Showing 1–48 of 48 results for author: Desch, S J

The first white dwarf debris disk observed by JWST

Authors: Andrew Swan, Jay Farihi, Kate Y. L. Su, Steven J. Desch

Abstract: This letter reports the first JWST spectroscopy of a white dwarf debris disk, giving a preliminary assessment of the salient features, and recommendations for future observations. The polluted and dusty star WD 0145+234 experienced a major collisional event in its circumstellar disk in 2018, accompanied by an infrared outburst, and subsequently a gradual decrease in thermal emission. Time-series N… ▽ More This letter reports the first JWST spectroscopy of a white dwarf debris disk, giving a preliminary assessment of the salient features, and recommendations for future observations. The polluted and dusty star WD 0145+234 experienced a major collisional event in its circumstellar disk in 2018, accompanied by an infrared outburst, and subsequently a gradual decrease in thermal emission. Time-series NIRSPEC observations demonstrate that the circumstellar disk is returning to a quiescent state with a T~1000 K infrared excess similar to the bulk of known dusty white dwarfs. MIRI spectroscopy reveals a 9-12 micron solid-state emission feature consistent with silicate minerals as observed in debris disks observed with Spitzer IRS. The strength and morphology of the silicate feature appear unchanged relative to the continuum in spectra taken over a year apart, consistent with steady-state collisional evolution of the circumstellar debris. A tentative emission feature around 7 microns may be due to carbonates, and if confirmed would indicate aqueous alteration in the parent body. △ Less

Submitted 19 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS Letters

Planetesimals at DZ stars I: chondritic compositions and a massive accretion event

Authors: Andrew Swan, Jay Farihi, Carl Melis, Patrick Dufour, Steven J. Desch, Detlev Koester, Jincheng Guo

Abstract: There is a wealth of evidence to suggest that planetary systems can survive beyond the main sequence. Most commonly, white dwarfs are found to be accreting material from tidally disrupted asteroids, whose bulk compositions are reflected by the metals polluting the stellar photospheres. While many examples are known, most lack the deep, high-resolution data required to detect multiple elements, and… ▽ More There is a wealth of evidence to suggest that planetary systems can survive beyond the main sequence. Most commonly, white dwarfs are found to be accreting material from tidally disrupted asteroids, whose bulk compositions are reflected by the metals polluting the stellar photospheres. While many examples are known, most lack the deep, high-resolution data required to detect multiple elements, and thus characterise the planetesimals that orbit them. Here, spectra of seven DZ white dwarfs observed with Keck HIRES are analysed, where up to nine metals are measured per star. Their compositions are compared against those of solar system objects, working in a Bayesian framework to infer or marginalise over the accretion history. All of the stars have been accreting primitive material, similar to chondrites, with hints of a Mercury-like composition at one star. The most polluted star is observed several Myr after its last major accretion episode, in which a Moon-sized object met its demise. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: Accepted to MNRAS. 18 pages, 9 figures

Origin of Low-${}^{26}{\rm Al}/{}^{27}{\rm Al}$ Corundum/Hibonite Inclusions in Meteorites

Authors: Steven J. Desch, Emilie T. Dunham, Ashley K. Herbst, Cayman T. Unterborn, Thomas G. Sharp, Maitrayee Bose, Prajkta Mane, Curtis D. Williams

Abstract: Most meteoritic calcium-rich, aluminum-rich inclusions (CAIs) formed from a reservoir with ${}^{26}{\rm Al}/{}^{27}{\rm Al} \approx 5 \times 10^{-5}$, but some record lower $({}^{26}{\rm Al}/{}^{27}{\rm Al})_0$, demanding they sampled a reservoir without live ${}^{26}{\rm Al}$. This has been interpreted as evidence for "late injection" of supernova material into our protoplanetary disk. We instead… ▽ More Most meteoritic calcium-rich, aluminum-rich inclusions (CAIs) formed from a reservoir with ${}^{26}{\rm Al}/{}^{27}{\rm Al} \approx 5 \times 10^{-5}$, but some record lower $({}^{26}{\rm Al}/{}^{27}{\rm Al})_0$, demanding they sampled a reservoir without live ${}^{26}{\rm Al}$. This has been interpreted as evidence for "late injection" of supernova material into our protoplanetary disk. We instead interpret the heterogeneity as chemical, demonstrating that these inclusions are strongly associated with the refractory phases corundum or hibonite. We name them "Low-${}^{26}{\rm Al}/{}^{27}{\rm Al}$ Corundum/Hibonite Inclusions" (LAACHIs). We present a detailed astrophysical model for LAACHI formation in which they derive their Al from presolar corundum, spinel or hibonite grains $0.5 - 2 \, μ{\rm m}$ in size with no live ${}^{26}{\rm Al}$; live ${}^{26}{\rm Al}$ is carried on smaller ($<$50 nm) presolar chromium spinel grains from recent nearby Wolf-Rayet stars or supernovae. In hot ($\approx$ 1350-1425 K) regions of the disk these grains, and perovskite grains, would be the only survivors. These negatively charged grains would grow to sizes $1 - 10^3 \, μ{\rm m}$, even incorporating positively charged perovskite grains, but not the small, negatively charged ${}^{26}{\rm Al}$-bearing grains. Chemical and isotopic fractionations due to grain charging was a significant process in hot regions of the disk. Our model explains the sizes, compositions, oxygen isotopic signatures, and the large, correlated ${}^{48}{\rm Ca}$ and ${}^{50}{\rm Ti}$ anomalies (if carried by presolar perovskite) of LAACHIs, and especially how they incorporated no ${}^{26}{\rm Al}$ in a solar nebula with uniform, canonical ${}^{26}{\rm Al}/{}^{27}{\rm Al}$. A late injection of supernova material is obviated, although formation of the Sun in a high-mass star-forming region is demanded. △ Less

Submitted 15 July, 2023; originally announced July 2023.

Comments: Accepted to The Astrophysical Journal

The Nominal Range of Rocky Planet Masses, Radii, Surface Gravities and Bulk Densities

Authors: Cayman T. Unterborn, Steven J. Desch, Jonas Haldemann, Alejandro Lorenzo, Joseph G. Schulze, Natalie R. Hinkel, Wendy R. Panero

Abstract: The two primary observable quantities of an exoplanet--its mass and radius--alone are not sufficient to probe a rocky exoplanet's interior composition and mineralogy. To overcome this, host-star abundances of the primary planet-building elements (Mg, Si, Fe) are typically used as a proxy for the planet's bulk composition. The majority of small exoplanet hosts, however, do not have available abunda… ▽ More The two primary observable quantities of an exoplanet--its mass and radius--alone are not sufficient to probe a rocky exoplanet's interior composition and mineralogy. To overcome this, host-star abundances of the primary planet-building elements (Mg, Si, Fe) are typically used as a proxy for the planet's bulk composition. The majority of small exoplanet hosts, however, do not have available abundance data. Here we present the open-source ExoPlex mass-radius-composition solver. Unlike previous open-source mass-radius solvers, ExoPlex calculates the core chemistry and equilibrium mantle mineralogy for a bulk composition, including effects of mantle FeO content, core light elements and surface water/ice. We utilize ExoPlex to calculate the planetary radii, surface gravities and bulk densities for 10$^6$ model planets up to 2 R$_\oplus$ across these geochemistries, adopting the distribution of FGK stellar abundances to estimate of the range of bulk exoplanet compositions. We outline the $99.7\%$ distribution of radii, surface gravity and bulk densities that define planets as "nominally rocky." Planets outside this range require compositions outside those expected from stellar abundance data, likely making them either Fe-enriched super-Mercuries, or volatile-enriched mini-Neptunes. We apply our classification scheme to a sample of 85 well-resolved exoplanets without available host-star abundances. We estimate only 9 planets are within the "nominally rocky planet zone" at $>70\%$ confidence, while $\sim20\%$ and $\sim30\%$ of this sample can be reasonably classified as super-Mercuries or volatile-rich, respectively. Our results provide observers with a self-consistent way to broadly classify a planet as likely rocky, Mercury-like or volatile-enriched, using mass and radius measurements alone. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 21 figures, 2 tables. Accepted to ApJ

Statistical Chronometry of Meteorites. I. The Pb-Pb age of the Solar System's t=0

Authors: Steven J. Desch, Daniel R. Dunlap, Emilie T. Dunham, Curtis D. Williams, Prajkta Mane

Abstract: We use rapidly cooled achondrites to test the assumption of 26Al homogeneity in the solar nebula, by checking if there is a single value of tSS, the absolute "Pb-Pb" age of the Solar System's t=0, that makes concordant their ages from the Al-Mg and Pb-Pb systems. We find that values tSS = 4568.42 +/-0.24 Myr do make these ages concordant, and therefore the hypothesis of homogeneous 26Al is not fal… ▽ More We use rapidly cooled achondrites to test the assumption of 26Al homogeneity in the solar nebula, by checking if there is a single value of tSS, the absolute "Pb-Pb" age of the Solar System's t=0, that makes concordant their ages from the Al-Mg and Pb-Pb systems. We find that values tSS = 4568.42 +/-0.24 Myr do make these ages concordant, and therefore the hypothesis of homogeneous 26Al is not falsified. This age, defined to be when the solar nebula had (26Al/27Al) = 5.23 x 10^-5, is significantly older than the ~4567.3 Myr inferred from direct measurements of Pb-Pb ages in CAIs. Discrepancies between the Al-Mg and Pb-Pb chronometers in chondrules and CAIs have previously been interpreted as arising from heterogeneities in 26Al, under the presumption that the Al-Mg and Pb-Pb systems in CAIs closed simultaneously. We examine this assumption and show that resetting is to be expected in CAIs. In particular, we quantitatively demonstrate that it is plausible that Pb-Pb ages of CAIs were reset at late times, without resetting the earlier Al-Mg ages, if they were transiently heated in the same manner as chondrules. We critically examine Pb-Pb isochrons, refining data and suggesting best practices for their calculation and reporting. We advocate reporting chronometry as times of formation after t=0 rather than absolute ages, as only the former is useful for astrophysical models of the solar nebula. We strongly advocate averaging of multiple samples, rather than anchoring to individual meteorites, to improve precision. △ Less

Submitted 3 May, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: Accepted by Icarus

Statistical Chronometry of Meteorites. II. Initial Abundances and Homogeneity of Short-lived Radionuclides

Authors: Steven J. Desch, Daniel R. Dunlap, Curtis D. Williams, Prajkta Mane, Emilie T. Dunham

Abstract: Astrophysical models of planet formation require accurate radiometric dating of meteoritic components by short-lived (Al-Mg, Mn-Cr, Hf-W) and long-lived (Pb-Pb) chronometers, to develop a timeline of such events in the solar nebula as formation of Ca-rich, Al-rich Inclusions (CAIs), chondrules, planetesimals, etc. CAIs formed mostly around a time ("t=0") when the short-lived radionuclide 26Al (t1/… ▽ More Astrophysical models of planet formation require accurate radiometric dating of meteoritic components by short-lived (Al-Mg, Mn-Cr, Hf-W) and long-lived (Pb-Pb) chronometers, to develop a timeline of such events in the solar nebula as formation of Ca-rich, Al-rich Inclusions (CAIs), chondrules, planetesimals, etc. CAIs formed mostly around a time ("t=0") when the short-lived radionuclide 26Al (t1/2 = 0.72 Myr) was present and presumably homogeneously distributed at a known level we define as (26Al/27Al)SS = 5.23 x10^-5. The time of formation after t=0 of another object can be found by determining its initial (26Al/27Al)0 ratio and comparing it to (26Al/27Al)SS. Dating of meteoritic objects using the Mn-Cr or Hf-W systems is hindered because the abundances (53Mn/55Mn)SS and (182Hf/180Hf)SS at t=0 are not known precisely. To constrain these quantities, we compile literature Al-Mg, Mn-Cr, Hf-W and Pb-Pb data for 14 achondrites and use novel statistical techniques to minimize the discrepancies between their times of formation across these systems. We find that for (53Mn/55Mn)SS = (8.09+/-0.65)x10^-6, (182Hf/180Hf)SS = (10.42+/-0.25)x10^-5, tSS = 4568.36+/-0.20 Myr, and a 53Mn half-life of 3.80+/-0.23 Myr, these four free parameters make concordant 37 formation times recorded by the different systems in 14 achondrites. These parameters also make concordant the ages derived for chondrules from CB/CH achondrites, formed simultaneously in an impact, and are apparently concordant with the I-Xe chronometer as well. Our findings provide very strong support for homogeneity of 26Al, 53Mn, and 182Hf in the solar nebula, and our approach offers a framework for more precise chronometry. △ Less

Submitted 8 May, 2023; v1 submitted 30 November, 2022; originally announced December 2022.

Comments: Accepted by Icarus

A primordial atmospheric origin of hydrospheric deuterium enrichment on Mars

Authors: Kaveh Pahlevan, Laura Schaefer, Linda T. Elkins-Tanton, Steven J. Desch, Peter R. Buseck

Abstract: The deuterium-to-hydrogen (D/H or 2H/1H) ratio of Martian atmospheric water (~6x standard mean ocean water, SMOW) is higher than that of known sources, requiring planetary enrichment. A recent measurement by NASA's Mars Science Laboratory rover Curiosity of >3 Gyr clays yields a D/H ratio ~3x SMOW, demonstrating that most enrichment occurs early in Mars's history. As on Venus, Mars's D/H enrichmen… ▽ More The deuterium-to-hydrogen (D/H or 2H/1H) ratio of Martian atmospheric water (~6x standard mean ocean water, SMOW) is higher than that of known sources, requiring planetary enrichment. A recent measurement by NASA's Mars Science Laboratory rover Curiosity of >3 Gyr clays yields a D/H ratio ~3x SMOW, demonstrating that most enrichment occurs early in Mars's history. As on Venus, Mars's D/H enrichment is thought to reflect preferential loss to space of 1H (protium) relative to 2H (deuterium), but the global environmental context of large and early hydrogen losses remain to be determined. Here, we apply a recent model of primordial atmosphere evolution to Mars, link the magma ocean of the accretion epoch with a subsequent water-ocean epoch, and calculate the behavior of deuterium for comparison with the observed record. We find that a ~2-3x hydrospheric deuterium-enrichment is produced if the Martian magma ocean is chemically reducing at last equilibration with the primordial atmosphere, making H2-CO the initially dominant species, with minor abundances of H2O-CO2. Reducing gases - in particular H2 - can cause greenhouse warming and prevent a water ocean from freezing immediately after the magma ocean epoch. Moreover, the pressure-temperature conditions are high enough to produce ocean-atmosphere H2O-H2 isotopic equilibrium such that surface H2O strongly concentrates deuterium relative to H2, which preferentially takes up protium and escapes from the primordial atmosphere. The proposed scenario of primordial H2-rich outgassing and escape suggests significant durations (>Myr) of chemical conditions on the Martian surface conducive to prebiotic chemistry immediately following Martian accretion. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 5 figures

Mantle Degassing Lifetimes through Galactic Time and the Maximum Age Stagnant-lid Rocky Exoplanets can Support Temperate Climates

Authors: Cayman T. Unterborn, Bradford J. Foley, Steven J. Desch, Patrick A. Young, Gregory Vance, Lee Chieffle, Stephen R. Kane

Abstract: The ideal exoplanets to search for life are those within a star's habitable zone. However, even within the habitable zone planets can still develop uninhabitable climate states. Sustaining a temperate climate over geologic ($\sim$Gyr) timescales requires a planet contain sufficient internal energy to power a planetary-scale carbon cycle. A major component of a rocky planet's energy budget is the h… ▽ More The ideal exoplanets to search for life are those within a star's habitable zone. However, even within the habitable zone planets can still develop uninhabitable climate states. Sustaining a temperate climate over geologic ($\sim$Gyr) timescales requires a planet contain sufficient internal energy to power a planetary-scale carbon cycle. A major component of a rocky planet's energy budget is the heat produced by the decay of radioactive elements, especially $^{40}$K, $^{232}$Th, $^{235}$U and $^{238}$U. As the planet ages and these elements decay, this radiogenic energy source dwindles. Here we estimate the probability distribution of the amount of these heat producing elements (HPEs) that enter into rocky exoplanets through Galactic history, by combining the system-to-system variation seen in stellar abundance data with the results from Galactic chemical evolution models. Using these distributions, we perform Monte-Carlo thermal evolution models that maximize the mantle cooling rate. This allows us to create a pessimistic estimate of lifetime a rocky, stagnant-lid exoplanet can support a global carbon cycle and temperate climate as a function of its mass and when it in Galactic history. We apply this framework to a sample of 17 likely rocky exoplanets with measured ages, 7 of which we predict are likely to be actively degassing today despite our pessimistic assumptions. For the remaining planets, including those orbiting TRAPPIST-1, we cannot confidently assume they currently contain sufficient internal heat to support mantle degassing at a rate sufficient to sustain a global carbon cycle or temperate climate without additional tidal heating or undergoing plate tectonics. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted to ApJ Letters

Short-Lived Radionuclides in Meteorites and the Sun's Birth Environment

Authors: Steven J. Desch, Edward D. Young, Emilie T. Dunham, Yusuke Fujimoto, Daniel R. Dunlap

Abstract: The solar nebula contained a number of short-lived radionuclides (SLRs) with half-lives of tens of Myr or less, comparable to the timescales for formation of protostars and protoplanetary disks. Therefore, determining the origins of SLRs would provide insights into star formation and the Sun's astrophysical birth environment. In this chapter, we review how isotopic studies of meteorites reveal the… ▽ More The solar nebula contained a number of short-lived radionuclides (SLRs) with half-lives of tens of Myr or less, comparable to the timescales for formation of protostars and protoplanetary disks. Therefore, determining the origins of SLRs would provide insights into star formation and the Sun's astrophysical birth environment. In this chapter, we review how isotopic studies of meteorites reveal the existence and abundances of these now-extinct radionuclides; and the evidence that the SLR ${}^{10}{\rm Be}$, which uniquely among the SLRs is not produced during typical stellar nucleosynthesis, was distributed homogeneously in the solar nebula. We review the evidence that the SLRs ${}^{26}{\rm Al}$, ${}^{53}{\rm Mn}$, and ${}^{182}{\rm Hf}$, and other radionuclides, were also homogeneously distributed and can be used to date events during the Solar System's planet-forming epoch. The homogeneity of the SLRs, especially ${}^{10}{\rm Be}$, strongly suggests they were all inherited from the Sun's molecular cloud, and that production by irradiation within the solar nebula was very limited, except for ${}^{36}{\rm Cl}$. We review astrophysical models for the origin of ${}^{10}{\rm Be}$, showing that it requires that the Sun formed in a spiral arm of the Galaxy with higher star formation rate than the Galaxy-wide average. Likewise, we review the astrophysical models for the origins of the other SLRs and show that they likely arose from contamination of the Sun's molecular cloud by massive stars over tens of Myr, most likely dominated by ejecta from Wolf-Rayet stars. The other SLRs also demand formation of the Sun in a spiral arm of the Galaxy with a star formation rate as high as demanded by the Solar System initial ${}^{10}{\rm Be}$ abundance. We discuss the astrophysical implications, and suggest further tests of these models and future directions for the field. △ Less

Submitted 15 January, 2023; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: 37 pages, 8 figures. Submitted as a chapter to Protostars and Planets VII

Matters Arising on "Breakup of a long-period comet as the origin of the dinosaur extinction" by Siraj & Loeb

Authors: Steven J. Desch, Alan P Jackson, Jessica L. Noviello, Ariel D. Anbar

Abstract: The recent publication by Siraj & Loeb (2021; Nature Scientific Reports 11, 3803) attempts to revive the debate over whether the Chicxulub impactor was a comet or an asteroid. They calculate that ~20% of long-period comets impacting Earth will have first been disrupted by passage inside the Sun's Roche limit, generating thousands of fragments, each the needed size of the Chicxulub impactor. This w… ▽ More The recent publication by Siraj & Loeb (2021; Nature Scientific Reports 11, 3803) attempts to revive the debate over whether the Chicxulub impactor was a comet or an asteroid. They calculate that ~20% of long-period comets impacting Earth will have first been disrupted by passage inside the Sun's Roche limit, generating thousands of fragments, each the needed size of the Chicxulub impactor. This would increase the impact rate of comets by a factor ~15, making them as likely to hit the Earth as an asteroid. They also argue that a comet would be a factor of 10 more likely to match the geochemical constraints, which indicate the Chicxulub impactor was carbonaceous chondrite-like. These conclusions are based on misinterpretations of the literature. Siraj & Loeb [1] overestimate the number of fragments produced during tidal disruption of a comet: tens of fragments are produced, not thousands. They also conflate 'carbonaceous chondrite' with specific types of carbonaceous chondrite, and ignore the evidence of iridium, making comets seem more likely than asteroids to match the Chicxulub impactor, when in fact they likely can be ruled out. Rather than a comet, an asteroidal impactor similar to CM or CR carbonaceous chondrites is strongly favored. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 4 pages. Submitted to Nature Scientific Reports

1I/'Oumuamua as an N2 ice fragment of an exo-pluto surface. II. Generation of N2 ice fragments and the origin of 'Oumuamua

Authors: Steven J. Desch, Alan P. Jackson

Abstract: The origin of the interstellar object 1I/'Oumuamua, has defied explanation. In a companion paper (Jackson & Desch, 2021), we show that a body of N2 ice with axes 45 m x 44 m x 7.5 m at the time of observation would be consistent with its albedo, non-gravitational acceleration, and lack of observed CO or CO2 or dust. Here we demonstrate that impacts on the surfaces of Pluto-like Kuiper belt objects… ▽ More The origin of the interstellar object 1I/'Oumuamua, has defied explanation. In a companion paper (Jackson & Desch, 2021), we show that a body of N2 ice with axes 45 m x 44 m x 7.5 m at the time of observation would be consistent with its albedo, non-gravitational acceleration, and lack of observed CO or CO2 or dust. Here we demonstrate that impacts on the surfaces of Pluto-like Kuiper belt objects (KBOs) would have generated and ejected ~10^14 collisional fragments--roughly half of them H2O ice fragments and half of them N2 ice fragments--due to the dynamical instability that depleted the primordial Kuiper belt. We show consistency between these numbers and the frequency with which we would observe interstellar objects like 1I/'Oumuamua, and more comet-like objects like 2I/Borisov, if other stellar systems eject such objects with efficiency like that of the Sun; we infer that differentiated KBOs and dynamical instabilities that eject impact-generated fragments may be near-universal among extrasolar systems. Galactic cosmic rays would erode such fragments over 4.5 Gyr, so that fragments are a small fraction (~0.1%) of long-period Oort comets, but C/2016 R2 may be an example. We estimate 'Oumuamua was ejected about 0.4-0.5 Gyr ago, from a young (~10^8 yr) stellar system, which we speculate was in the Perseus arm. Objects like 'Oumuamua may directly probe the surface compositions of a hitherto-unobserved type of exoplanet: "exo-plutos". 'Oumuamua may be the first sample of an exoplanet brought to us. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 29 pages, 1 figure. Companion to Jackson & Desch (2021), "1I/'Oumuamua as an N2 ice fragment of an exo-pluto surface. I. Size and compositional constraints"

1I/'Oumuamua as an N2 ice fragment of an exo-Pluto surface: I. size and compositional constraints

Authors: Alan P. Jackson, Steven J. Desch

Abstract: The origin of the interstellar object 1I/'Oumuamua has defied explanation. We perform calculations of the non-gravitational acceleration that would be experienced by bodies composed of a range of different ices and demonstrate that a body composed of N2 ice would satisfy the available constraints on the non-gravitational acceleration, size and albedo, and lack of detectable emission of CO or CO2 o… ▽ More The origin of the interstellar object 1I/'Oumuamua has defied explanation. We perform calculations of the non-gravitational acceleration that would be experienced by bodies composed of a range of different ices and demonstrate that a body composed of N2 ice would satisfy the available constraints on the non-gravitational acceleration, size and albedo, and lack of detectable emission of CO or CO2 or dust. We find that 'Oumuamua was small, with dimensions 45 m x 44 m x 7.5 m at the time of observation at 1.42 au from the Sun, with a high albedo of 0.64. This albedo is consistent with the N2 surfaces of bodies like Pluto and Triton. We estimate 'Oumuamua was ejected about 0.4-0.5 Gyr ago from a young stellar system, possibly in the Perseus arm. Objects like 'Oumuamua may directly probe the surface compositions of a hitherto-unobserved type of exoplanet: "exo-plutos". In a companion paper (Desch & Jackson, 2021) we demonstrate that dynamical instabilities like the one experienced by the Kuiper belt, in other stellar systems, plausibly could generate and eject large numbers of N2 ice fragments. 'Oumuamua may be the first sample of an exoplanet brought to us. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 24 pages, 3 figures, 2 tables. Accepted for publication in JGR: Planets. Code available at https://doi.org/10.5281/zenodo.45584

Icy Exomoons Evidenced by Spallogenic Nuclides in Polluted White Dwarfs

Authors: Alexandra E. Doyle, Steven J. Desch, Edward D. Young

Abstract: We present evidence that excesses in Be in polluted white dwarfs (WDs) are the result of accretion of icy exomoons that formed in the radiation belts of giant exoplanets. Here we use excess Be in the white dwarf GALEX J2339-0424 as an example. We constrain the parent body abundances of rock-forming elements in GALEX J2339-0424 and show that the overabundance of beryllium in this WD cannot be accou… ▽ More We present evidence that excesses in Be in polluted white dwarfs (WDs) are the result of accretion of icy exomoons that formed in the radiation belts of giant exoplanets. Here we use excess Be in the white dwarf GALEX J2339-0424 as an example. We constrain the parent body abundances of rock-forming elements in GALEX J2339-0424 and show that the overabundance of beryllium in this WD cannot be accounted for by differences in diffusive fluxes through the WD outer envelope nor by chemical fractionations during typical rock-forming processes. We argue instead that the Be was produced by energetic proton irradiation of ice mixed with rock. We demonstrate that the MeV proton fluence required to form the high Be/O ratio in the accreted parent body is consistent with irradiation of ice in the rings of a giant planet within its radiation belt, followed by accretion of the ices to form a moon that is later accreted by the WD. The icy moons of Saturn serve as useful analogs. Our results provide an estimate of spallogenic nuclide excesses in icy moons formed by rings around giant planets in general, including those in the solar system. While excesses in Be have been detected in two polluted WDs to date, including the WD described here, we predict that excesses in the other spallogenic elements Li and B, although more difficult to detect, should also be observed, and that such detections would also indicate pollution by icy exomoons formed in the ring systems of giant planets. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 13 pages, 4 figures, 1 table

Journal ref: ApJL 907 L35 (2021)

The Fundamental Connections Between the Solar System and Exoplanetary Science

Authors: Stephen R. Kane, Giada N. Arney, Paul K. Byrne, Paul A. Dalba, Steven J. Desch, Jonti Horner, Noam R. Izenberg, Kathleen E. Mandt, Victoria S. Meadows, Lynnae C. Quick

Abstract: Over the past several decades, thousands of planets have been discovered outside of our Solar System. These planets exhibit enormous diversity, and their large numbers provide a statistical opportunity to place our Solar System within the broader context of planetary structure, atmospheres, architectures, formation, and evolution. Meanwhile, the field of exoplanetary science is rapidly forging onw… ▽ More Over the past several decades, thousands of planets have been discovered outside of our Solar System. These planets exhibit enormous diversity, and their large numbers provide a statistical opportunity to place our Solar System within the broader context of planetary structure, atmospheres, architectures, formation, and evolution. Meanwhile, the field of exoplanetary science is rapidly forging onward towards a goal of atmospheric characterization, inferring surface conditions and interiors, and assessing the potential for habitability. However, the interpretation of exoplanet data requires the development and validation of exoplanet models that depend on in-situ data that, in the foreseeable future, are only obtainable from our Solar System. Thus, planetary and exoplanetary science would both greatly benefit from a symbiotic relationship with a two-way flow of information. Here, we describe the critical lessons and outstanding questions from planetary science, the study of which are essential for addressing fundamental aspects for a variety of exoplanetary topics. We outline these lessons and questions for the major categories of Solar System bodies, including the terrestrial planets, the giant planets, moons, and minor bodies. We provide a discussion of how many of these planetary science issues may be translated into exoplanet observables that will yield critical insight into current and future exoplanet discoveries. △ Less

Submitted 8 August, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 53 pages, 5 figures, invited review, accepted for publication in the Journal of Geophysical Research: Planets

Journal ref: Journal of Geophysical Research: Planets, 2021, Volume 126, Issue 2, article id. e06643

Exogeoscience and Its Role in Characterizing Exoplanet Habitability and the Detectability of Life

Authors: Cayman T. Unterborn, Paul K. Byrne, Ariel D. Anbar, Giada Arney, David Brain, Steve J. Desch, Bradford J. Foley, Martha S. Gilmore, Hilairy E. Hartnett, Wade G. Henning, Marc M. Hirschmann, Noam R. Izenberg, Stephen R. Kane, Edwin S. Kite, Laura Kreidberg, Kanani K. M. Lee, Timothy W. Lyons, Stephanie L. Olson, Wendy R. Panero, Noah J. Planavsky, Christopher T. Reinhard, Joseph P. Renaud, Laura K. Schaefer, Edward W. Schwieterman, Linda E. Sohl , et al. (2 additional authors not shown)

Abstract: The search for exoplanetary life must encompass the complex geological processes reflected in an exoplanet's atmosphere, or we risk reporting false positive and false negative detections. To do this, we must nurture the nascent discipline of "exogeoscience" to fully integrate astronomers, astrophysicists, geoscientists, oceanographers, atmospheric chemists and biologists. Increased funding for int… ▽ More The search for exoplanetary life must encompass the complex geological processes reflected in an exoplanet's atmosphere, or we risk reporting false positive and false negative detections. To do this, we must nurture the nascent discipline of "exogeoscience" to fully integrate astronomers, astrophysicists, geoscientists, oceanographers, atmospheric chemists and biologists. Increased funding for interdisciplinary research programs, supporting existing and future multidisciplinary research nodes, and developing research incubators is key to transforming true exogeoscience from an aspiration to a reality. △ Less

Submitted 23 July, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: Submitted as white paper to 2023-2033 Planetary Science and Astrobiology Decadal Survey; Updated to include all co-authors

The Dynamic Proto-atmospheres around Low-Mass Planets with Eccentric Orbits

Authors: Chuhong Mai, Steven J. Desch, Rolf Kuiper, Gabriel-Dominique Marleau, Cornelis Dullemond

Abstract: Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We inves… ▽ More Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We investigate the proto-atmospheres on low-mass planets with eccentric orbits with radiation-hydrodynamics simulations. A 2D radiative model of the proto-atmosphere is established with tabulated opacities for the gas and dust. The solutions reveal large-scale gas recycling inside a bow shock structure. The atmospheres on eccentric planets are typically three to four orders of magnitude less massive than those of planets with circular orbits. Overall, however, a supersonic environment is favorable for planets to keep an early stable atmosphere, rather than harmful, due to the steady gas supply through the recycling flow. We also quantitatively explore how such atmospheres are affected by the relative velocity of the planet to the gas, the planet mass, and the background gas density. Our time-dependent simulations track the orbital evolution of the proto-atmosphere with the planet--disk parameters changing throughout the orbit. Atmospheric properties show oscillatory patterns as the planet travels on an eccentric orbit, with a lag in phase. To sum up, low-mass eccentric planets can retain small proto-atmospheres despite the stripping effects of bow shocks. The atmospheres are always connected to and interacting with the disk gas. These findings provide important insights into the impacts of migration and scattering on planetary proto-atmospheres. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJ (23 pages, 15 figures, 2 tables)

A Geologically Robust Procedure For Observing Rocky Exoplanets to Ensure that Detection of Atmospheric Oxygen is an Earth-Like Biosignature

Authors: Carey M. Lisse, Steven J. Desch, Cayman T. Unterborn, Stephen R. Kane, Patrick R. Young, Hilairy E. Hartnett, Natalie R. Hinkel, Sang Heon Shim, Eric E. Mamajek, Noam R. Izenberg

Abstract: In the next decades, the astrobiological community will debate whether the first observations of oxygen in an exoplanet$'$s atmosphere signifies life, so it is critical to establish procedures now for collection and interpretation of such data. We present a step-by-step observational strategy for using oxygen as a robust biosignature, to prioritize exoplanet targets and design future observations.… ▽ More In the next decades, the astrobiological community will debate whether the first observations of oxygen in an exoplanet$'$s atmosphere signifies life, so it is critical to establish procedures now for collection and interpretation of such data. We present a step-by-step observational strategy for using oxygen as a robust biosignature, to prioritize exoplanet targets and design future observations. It is premised on avoiding planets lacking subaerial weathering of continents, which would imply geochemical cycles drastically different from Earth$'$s, precluding use of oxygen as a biosignature. The strategy starts with the most readily obtained data: semi-major axis and stellar luminosity to ensure residence in the habitable zone; stellar XUV flux, to ensure an exoplanet can retain a secondary (outgassed) atmosphere. Next, high-precision mass and radius information should be combined with high-precision stellar abundance data, to constrain the exoplanet$'$s water content; those incompatible with less than 0.1 wt % H$_{2}$O can be deprioritized. Then, reflectance photometry or low-resolution transmission spectroscopy should confirm an optically thin atmosphere. Subsequent long-duration, high-resolution transmission spectroscopy should search for oxygen and ensure that water vapor and CO$_{2}$ are present only at low (10$^{2}$-10$^{4}$ ppm levels). Assuming oxygen is found, attribution to life requires the difficult acquisition of a detailed, multispectral light curve of the exoplanet to ensure both surface land and water. Exoplanets failing some of these steps might be habitable, even have observable biogenic oxygen, but should be deprioritized because oxygen could not be attributed unambiguously to life. We show how this is the case for the Solar System, the 55 Cnc System, and the TRAPPIST-1 System, in which only the Earth and TRAPPIST-1e successfully pass through our procedure. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 27 Pages, 1 Figure, 0 Tables (accepted 09 June 2020, in press)

Journal ref: Astrophysical Journal Letters 2020

Detectability of Life Using Oxygen on Pelagic Planets and Water Worlds

Authors: Donald M Glaser, Hilairy Ellen Hartnett, Steven J. Desch, Cayman T. Unterborn, Ariel Anbar, Steffen Buessecker, Theresa Fisher, Steven Glaser, Stephen R. Kane, Carey M. Lisse, Camerian Millsaps, Susanne Neuer, Joseph G. ORourke, Nuno Santos, Sara Imari Walker, Mikhail Zolotov

Abstract: The search for life on exoplanets is one of the grand scientific challenges of our time. The strategy to date has been to find (e.g., through transit surveys like Kepler) Earth-like exoplanets in their stars habitable zone, then use transmission spectroscopy to measure biosignature gases, especially oxygen, in the planets atmospheres (e.g., using JWST, the James Webb Space Telescope). Already ther… ▽ More The search for life on exoplanets is one of the grand scientific challenges of our time. The strategy to date has been to find (e.g., through transit surveys like Kepler) Earth-like exoplanets in their stars habitable zone, then use transmission spectroscopy to measure biosignature gases, especially oxygen, in the planets atmospheres (e.g., using JWST, the James Webb Space Telescope). Already there are more such planets than can be observed by JWST, and missions like the Transiting Exoplanet Survey Satellite and others will find more. A better understanding of the geochemical cycles relevant to biosignature gases is needed, to prioritize targets for costly follow-up observations and to help design future missions. We define a Detectability Index to quantify the likelihood that a biosignature gas could be assigned a biological vs. non-biological origin. We apply this index to the case of oxygen gas, O2, on Earth-like planets with varying water contents. We demonstrate that on Earth-like exoplanets with 0.2 weight percent (wt%) water (i.e., no exposed continents) a reduced flux of bioessential phosphorus limits the export of photosynthetically produced atmospheric O2 to levels indistinguishable from geophysical production by photolysis of water plus hydrogen escape. Higher water contents >1wt% that lead to high-pressure ice mantles further slow phosphorus cycling. Paradoxically, the maximum water content allowing use of O2 as a biosignature, 0.2wt%, is consistent with no water based on mass and radius. Thus, the utility of an O2 biosignature likely requires the direct detection of both water and land on a planet. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Haumea's Shape and Composition

Authors: Emilie T. Dunham, Steven J. Desch, Luke Probst

Abstract: We have calculated the figure of equilibrium of a rapidly rotating, differentiated body to determine the shape, structure, and composition of the dwarf planet Haumea. Previous studies of Haumea's light curve have suggested Haumea is a uniform triaxial ellipsoid consistent with a Jacobi ellipsoid with axes $\approx 960 \times 774 \times 513$ km, and bulk density… ▽ More We have calculated the figure of equilibrium of a rapidly rotating, differentiated body to determine the shape, structure, and composition of the dwarf planet Haumea. Previous studies of Haumea's light curve have suggested Haumea is a uniform triaxial ellipsoid consistent with a Jacobi ellipsoid with axes $\approx 960 \times 774 \times 513$ km, and bulk density $\approx 2600 \, {\rm kg} \, {\rm m}^{-3}$. In contrast, observations of a recent stellar occultation by Haumea indicate its axes are $\approx 1161 \times 852 \times 523$ km and its bulk density $\approx 1885 \, {\rm kg} \, {\rm m}^{-3}$; these results suggest that Haumea cannot be a fluid in hydrostatic equilibrium and must be partially supported by interparticle forces. We have written a code to reconcile these contradictory results and to determine if Haumea is in fact a fluid in hydrostatic equilibrium. The code calculates the equilibrium shape, density, and ice crust thickness of a differentiated Haumea after imposing (semi-) axes lengths $a$ and $b$. We find Haumea is consistent with a differentiated triaxial ellipsoid fluid in hydrostatic equilibrium with axes of best fit $a$ = 1050 km, $b$ = 840 km, and $c$ = 537 km. This solution for Haumea has $ρ_{\rm avg} = 2018 \, {\rm kg} \, {\rm m}^{-3}$, $ρ_{\rm core} = 2680 \, {\rm kg} \, {\rm m}^{-3}$, and core axes $a_{\rm c} = 883$ km, $b_{\rm c} = 723$ km, and $c_{\rm c}=470$ km, which equates to an ice mantle comprising $\sim 17\%$ of Haumea's volume and ranging from 67 to 167 km in thickness. The thick ice crust we infer allows for Haumea's collisional family to represent only a small fraction of Haumea's pre-collisional ice crust. For a wide range of parameters, the core density we calculate for Haumea suggests that today the core is composed of hydrated silicates and likely underwent serpentinization in the past. △ Less

Submitted 31 March, 2019; originally announced April 2019.

Comments: 15 pages, 3 figures, accepted for publication in ApJ

Effect of Different Angular Momentum Transport mechanisms on the Distribution of Water in Protoplanetary Disks

Authors: Anusha Kalyaan, Steven J. Desch

Abstract: The snow line in a protoplanetary disk demarcates regions with H$_2$O ice from regions with H$_2$O vapor. Where a planet forms relative to this location determines how much water and other volatiles it forms with. Giant planet formation may be triggered at the water snow line if vapor diffuses outward and is cold-trapped beyond the snow line faster than icy particles can drift inward. In this stud… ▽ More The snow line in a protoplanetary disk demarcates regions with H$_2$O ice from regions with H$_2$O vapor. Where a planet forms relative to this location determines how much water and other volatiles it forms with. Giant planet formation may be triggered at the water snow line if vapor diffuses outward and is cold-trapped beyond the snow line faster than icy particles can drift inward. In this study we investigate the distribution of water across the snow line, considering three different radial profiles of the turbulence parameter $α(r)$, corresponding to three different angular momentum transport mechanisms. We consider the radial transport of water vapor and icy particles by diffusion, advection, and drift. We show that even for similar values of $α$, the gradient of $α$(r) across the snow line significantly changes the snow line location, the sharpness of the volatile gradient across the snow line, and the final water/rock ratio in planetary bodies. A profile of radially decreasing $α$, consistent with transport by hydrodynamic instabilities plus magnetic disk winds, appears consistent with the distribution of water in the solar nebula, with monotonically-increasing radial water content and a diverse population of asteroids with different water content. We argue that $Σ(r)$ and water abundance $N_{\rm H_2O}(r)/N_{\rm H_2}(r)$ are likely diagnostic of $α(r)$ and thus the mechanism for angular momentum transport in inner disks. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Accepted to ApJ (Mar 7, 2019)

Protoplanetary Disk Science Enabled by Extremely Large Telescopes

Authors: Hannah Jang-Condell, Sean Brittain, Alycia Weinberger, Michael Liu, Jacqueline Faherty, Jaehan Bae, Sean Andrews, Megan Ansdell, Til Birnstiel, Alan Boss, Laird Close, Thayne Currie, Steven J Desch, Sarah Dodson-Robinson, Chuanfei Dong, Gaspard Duchene, Catherine Espaillat, Kate Follette, Eric Gaidos, Peter Gao, Nader Haghighipour, Hilairy Hartnett, Yasuhiro Hasegawa, Mihkel Kama, Jinyoung Serena Kim , et al. (14 additional authors not shown)

Abstract: The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the ph… ▽ More The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Inferring the Composition of Disintegrating Planet Interiors from Dust Tails with Future James Webb Space Telescope Observations

Authors: Eva H. L. Bodman, Jason T. Wright, Steven J. Desch, Carey M. Lisse

Abstract: Disintegrating planets allow for the unique opportunity to study the composition of the interiors of small, hot, rocky exoplanets because the interior is evaporating and that material is condensing into dust, which is being blown away and then transiting the star. Their transit signal is dominated by dusty effluents forming a comet-like tail trailing the host planet (or leading it, in the case of… ▽ More Disintegrating planets allow for the unique opportunity to study the composition of the interiors of small, hot, rocky exoplanets because the interior is evaporating and that material is condensing into dust, which is being blown away and then transiting the star. Their transit signal is dominated by dusty effluents forming a comet-like tail trailing the host planet (or leading it, in the case of K2-22b), making these good candidates for transmission spectroscopy. To assess the ability of such observations to diagnose the dust composition, we simulate the transmission spectra from 5-14 $μ$m for the planet tail assuming an optically-thin dust cloud comprising a single dust species with a constant column density scaled to yield a chosen visible transit depth. We find that silicate resonant features near 10 $μ$m can produce transit depths that are at least as large as those in the visible. For the average transit depth of 0.55% in the Kepler band for K2-22b, the features in the transmission spectra can be as large as 1%, which is detectable with the JWST MIRI low-resolution spectrograph in a single transit. The detectability of compositional features is easier with an average grain size of 1 $μ$m despite features being more prominent with smaller grain sizes. We find most features are still detectable for transit depths of ~0.3% in the visible range. If more disintegrating planets are found with future missions such as the space telescope TESS, follow-up observations with JWST can explore the range of planetary compositions. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: 14 pages, 3 figures

Updated Compositional Models of the TRAPPIST-1 Planets

Authors: Cayman T. Unterborn, Natalie R. Hinkel, Steven J. Desch

Abstract: After publication of our initial mass-radius-composition models for the TRAPPIST-1 system in Unterborn et al. (2018), the planet masses were updated in Grimm et al. (2018). We had originally adopted the data set of Wang et al., 2017 who reported different densities than the updated values. The differences in observed density change the inferred volatile content of the planets. Grimm et al. (2018)… ▽ More After publication of our initial mass-radius-composition models for the TRAPPIST-1 system in Unterborn et al. (2018), the planet masses were updated in Grimm et al. (2018). We had originally adopted the data set of Wang et al., 2017 who reported different densities than the updated values. The differences in observed density change the inferred volatile content of the planets. Grimm et al. (2018) report TRAPPIST-1 b, d, f, g, and h as being consistent with <5 wt% water and TRAPPIST-1 c and e has having largely rocky interiors. Here, we present updated results recalculating water fractions and potential alternative compositions using the Grimm et al., 2018 masses. Overall, we can only reproduce the results of Grimm et al., 2018 of planets b, d and g having small water contents if the cores of these planets are small (<23 wt%). We show that, if the cores for these planets are roughly Earth-sized (33 wt%), significant water fractions up to 40 wt% are possible. We show planets c, e, f, and h can have volatile envelopes between 0-35 wt% that are also consistent with being totally oxidized and lacking an Fe-core entirely. We note here that a pure MgSiO$_3$ planet (Fe/Mg = 0) is not the true lowest density end-member mass-radius curve for determining the probability of a planet containing volatiles. All planets that are rocky likely contain some Fe, either within the core or oxidized in the mantle. We argue the true low density end-member for oxidizing systems is instead a planet with the lowest reasonable Fe/Mg and completely core-less. Using this logic, we assert that planets b, d and g likely must have significant volatile layers because the end-member planet models produce masses too high even when uncertainties in both mass and radius are taken into account. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: Accepted to RNAAS

Magnetic Fields Recorded by Chondrules Formed in Nebular Shocks

Authors: Chuhong Mai, Steven J. Desch, Aaron C. Boley, Benjamin P. Weiss

Abstract: Recent laboratory efforts (Fu et al., 2014) have constrained the remanent magnetizations of chondrules and the magnetic field strengths at which the chondrules were exposed to as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investi… ▽ More Recent laboratory efforts (Fu et al., 2014) have constrained the remanent magnetizations of chondrules and the magnetic field strengths at which the chondrules were exposed to as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investigate the amplification of the magnetic field above background values for two proposed chondrule formation mechanisms, large-scale nebular shocks and planetary bow shocks. Behind large-scale shocks, the magnetic field parallel to the shock front is amplified by factors $\sim 10-30$, regardless of the magnetic diffusivity. Therefore, chondrules melted in these shocks probably recorded an amplified magnetic field. Behind planetary bow shocks, the field amplification is sensitive to the magnetic diffusivity. We compute the gas properties behind a bow shock around a 3000 km-radius planetary embryo, with and without atmospheres, using hydrodynamics models. We calculate the ionization state of the hot, shocked gas, including thermionic emission from dust, and thermal ionization of gas-phase potassium atoms, and the magnetic diffusivity due to Ohmic dissipation and ambipolar diffusion. We find that the diffusivity is sufficiently large that magnetic fields have already relaxed to background values in the shock downstream where chondrules acquire magnetizations, and that these locations are sufficiently far from the planetary embryos that chondrules should not have recorded a significant putative dynamo field generated on these bodies. We conclude that, if melted in planetary bow shocks, chondrules probably recorded the background nebular field. △ Less

Submitted 13 March, 2018; originally announced March 2018.

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

A procedure for observing rocky exoplanets to maximize the likelihood that atmospheric oxygen will be a biosignature

Authors: Steven J. Desch, Stephen Kane, Carey M. Lisse, Cayman T. Unterborn, Hilairy E. Hartnett, Sang-Heon Shim

Abstract: Here we advocate an observational strategy to help prioritize exoplanet observations. It starts with more easily obtained observational data, and ranks exoplanets for more difficult follow-up observations based on the likelihood of avoiding planets for which oxygen is a false positives or even an inconclusive signature of life. We find that for oxygen to be a robust biosignature, both land and sur… ▽ More Here we advocate an observational strategy to help prioritize exoplanet observations. It starts with more easily obtained observational data, and ranks exoplanets for more difficult follow-up observations based on the likelihood of avoiding planets for which oxygen is a false positives or even an inconclusive signature of life. We find that for oxygen to be a robust biosignature, both land and surface water must be present. Landless exoplanets have much slower biogeochemical cycles, so while oxygenic photosynthesizing life could exist on such planets, it could not produce oxygen at a rate competitive with abiotic rates such as photolysis. These habitable planets, whose life would not be detectable, should be avoided. △ Less

Submitted 21 January, 2018; originally announced January 2018.

Comments: 6 pages, 1 figure, white paper submitted in response to the solicitation of feedback for the "Astrobiology Science Strategy for the Search for Life in the Universe" by the National Academy of Sciences

The Effect of Jupiter's Formation on the Distribution of Refractory Elements and Inclusions in Meteorites

Authors: Steven J. Desch, Anusha Kalyaan, Conel M. O'D. Alexander

Abstract: We present a comprehensive evolutionary model of the Sun's protoplanetary disk, constructed to resolve the "CAI Storage" problem of meteoritics. We predict the abundances of calcium-rich, aluminum-rich inclusions (CAIs) and refractory lithophile elements under the central assumption that Jupiter's $\sim 30 \, M_{\oplus}$ core forms at about 3 AU at around 0.6 Myr and opened a gap CAIs are trapped… ▽ More We present a comprehensive evolutionary model of the Sun's protoplanetary disk, constructed to resolve the "CAI Storage" problem of meteoritics. We predict the abundances of calcium-rich, aluminum-rich inclusions (CAIs) and refractory lithophile elements under the central assumption that Jupiter's $\sim 30 \, M_{\oplus}$ core forms at about 3 AU at around 0.6 Myr and opened a gap CAIs are trapped in the pressure maximum beyond Jupiter; carbonaceous chondrites formed there. Inside Jupiter's orbit, CAIs were depleted by aerodynamic drag; ordinary and enstatite chondrites formed there. For 16 chondrites and achondrites, we review meteoritic data on their CAI and refractory abundances and their times of formation, constrained by radiometric dating and thermal models. We predict their formation locations, finding excellent consistency with other location information (water content, asteroid spectra and parent bodies). We predict the size of particle concentrated by turbulence for each chondrite, finding excellent matches to each chondrites's mean chondrule diameter. These consistencies imply meteorite parent bodies assembled quickly from local materials concentrated by turbulence, and usually did not migrate far. We predict CI chondrites are depleted in refractory lithophile elements relative to the Sun, by about 12% (0.06 dex). We constrain the variation of turbulence parameter $α$ in the disk, and find a limited role for magnetorotational instability, favoring hydrodynamical instabilities in the outer disk, plus magnetic disk winds in the inner disk. Between 3 and 4 Myr at least, gas persisted outside Jupiter but was depleted inside it, and the solar nebula was a transition disk. △ Less

Submitted 8 August, 2018; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 44 pages, 10 figures, 5 tables, accepted for ApJ Supplement

A Catalog of Stellar Unified Properties (CATSUP) for 951 FGK-Stars Within 30 pc

Authors: Natalie R. Hinkel, Eric E. Mamajek, Margaret C. Turnbull, Ella Osby, Evgenya L. Shkolnik, Graeme H. Smith, Alexis Klimasewski, Garrett Somers, Steven J. Desch

Abstract: Almost every star in our Galaxy is likely to harbor a terrestrial planet, but accurate measurements of an exoplanet's mass and radius demands accurate knowledge of the properties of its host star. The imminent TESS and CHEOPS missions are slated to discover thousands of new exoplanets. Along with WFIRST, which will directly image nearby planets, these surveys make urgent the need to better charact… ▽ More Almost every star in our Galaxy is likely to harbor a terrestrial planet, but accurate measurements of an exoplanet's mass and radius demands accurate knowledge of the properties of its host star. The imminent TESS and CHEOPS missions are slated to discover thousands of new exoplanets. Along with WFIRST, which will directly image nearby planets, these surveys make urgent the need to better characterize stars in the nearby solar neighborhood (< 30 pc). We have compiled the CATalog of Stellar Unified Properties (CATSUP) for 951 stars, including such data as: Gaia astrometry; multiplicity within stellar systems; stellar elemental abundance measurements; standardized spectral types; Ca II H and K stellar activity indices; GALEX NUV and FUV photometry; and X-ray fluxes and luminosities from ROSAT, XMM, and Chandra. We use this data-rich catalog to find correlations, especially between stellar emission indices, colors, and galactic velocity. Additionally, we demonstrate that thick-disk stars in the sample are generally older, have lower activity, and have higher velocities normal to the galactic plane. We anticipate CATSUP will be useful in discerning other trends among stars within the nearby solar neighborhood, for comparing thin-disk vs. thick-disk stars, for comparing stars with and without planets, and for finding correlations between chemical and kinematic properties. △ Less

Submitted 13 September, 2017; originally announced September 2017.

Comments: 21 pages, 14 figures, 6 tables, full CATSUP database included in source data (datafile6.txt) per breakdown in Table 6 and will also be uploaded on Vizier, accepted for publication in ApJ

Inward Migration of the TRAPPIST-1 Planets as Inferred From Their Water-Rich Compositions

Authors: Cayman T. Unterborn, Steven J. Desch, Natalie R. Hinkel, Alejandro Lorenzo Jr

Abstract: Multiple planet systems provide an ideal laboratory for probing exoplanet composition, formation history and potential habitability. For the TRAPPIST-1 planets, the planetary radii are well established from transits (Gillon et al., 2016, Gillon et al., 2017), with reasonable mass estimates coming from transit timing variations (Gillon et al., 2017, Wang et al., 2017) and dynamical modeling (Quarle… ▽ More Multiple planet systems provide an ideal laboratory for probing exoplanet composition, formation history and potential habitability. For the TRAPPIST-1 planets, the planetary radii are well established from transits (Gillon et al., 2016, Gillon et al., 2017), with reasonable mass estimates coming from transit timing variations (Gillon et al., 2017, Wang et al., 2017) and dynamical modeling (Quarles et al., 2017). The low bulk densities of the TRAPPIST-1 planets demand significant volatile content. Here we show using mass-radius-composition models, that TRAPPIST-1f and g likely contain substantial ($\geq50$ wt\%) water/ice, with b and c being significantly drier ($\leq15$ wt\%). We propose this gradient of water mass fractions implies planets f and g formed outside the primordial snow line whereas b and c formed inside. We find that compared to planets in our solar system that also formed within the snow line, TRAPPIST-1b and c contain hundreds more oceans worth of water. We demonstrate the extent and timescale of migration in the TRAPPIST-1 system depends on how rapidly the planets formed and the relative location of the primordial snow line. This work provides a framework for understanding the differences between the protoplanetary disks of our solar system versus M dwarfs. Our results provide key insights into the volatile budgets, timescales of planet formation, and migration history of likely the most common planetary host in the Galaxy. △ Less

Submitted 6 February, 2018; v1 submitted 8 June, 2017; originally announced June 2017.

Comments: Accepted to Nature Astronomy, 30 pages, 6 figures

Formulas for Radial Transport in Protoplanetary Disks

Authors: Steven J. Desch, Paul R. Estrada, Anusha Kalyaan, Jeffrey N. Cuzzi

Abstract: Quantification of the radial transport of gaseous species and solid particles is important to many applications in protoplanetary disk evolution. An especially important example is determining the location of the water snow lines in a disk, which requires computing the rates of outward radial diffusion of water vapor and the inward radial drift of icy particles; however, the application is general… ▽ More Quantification of the radial transport of gaseous species and solid particles is important to many applications in protoplanetary disk evolution. An especially important example is determining the location of the water snow lines in a disk, which requires computing the rates of outward radial diffusion of water vapor and the inward radial drift of icy particles; however, the application is generalized to evaporation fronts of all volatiles. We review the relevant formulas using a uniform formalism. This uniform treatment is necessary because the literature currently contains at least six mutually exclusive treatments of radial diffusion of gas, only one of which is correct. We derive the radial diffusion equations from first principles, using Fick's law. For completeness, we also present the equations for radial transport of particles. These equations may be applied to studies of diffusion of gases and particles in protoplanetary and other accretion disks. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: 13 pages, 4 Figures, accepted for ApJ

A Comparison of Stellar Elemental Abundance Techniques and Measurements

Authors: Natalie R. Hinkel, Patrick A. Young, Michael D. Pagano, Steven J. Desch, Ariel D. Anbar, Vardan Adibekyan, Sergi Blanco-Cuaresma, Joleen K. Carlberg, Elisa Delgado Mena, Fan Liu, Thomas Nordlander, Sergio G. Sousa, Andreas Korn, Pieter Gruyters, Ulrike Heiter, Paula Jofre, Nuno C. Santos, Caroline Soubiran

Abstract: Stellar elemental abundances are important for understanding the fundamental properties of a star or stellar group, such as age and evolutionary history, as well as the composition of an orbiting planet. However, as abundance measurement techniques have progressed, there has been little standardization between individual methods and their comparisons. As a result, different stellar abundance proce… ▽ More Stellar elemental abundances are important for understanding the fundamental properties of a star or stellar group, such as age and evolutionary history, as well as the composition of an orbiting planet. However, as abundance measurement techniques have progressed, there has been little standardization between individual methods and their comparisons. As a result, different stellar abundance procedures determine measurements that vary beyond quoted error for the same elements within the same stars (Hinkel et al. 2014). The purpose of this paper is to better understand the systematic variations between methods and offer recommendations for producing more accurate results in the future. We have invited a number of participants from around the world (Australia, Portugal, Sweden, Switzerland, and USA) to calculate ten element abundances (C, O, Na, Mg, Al, Si, Fe, Ni, Ba, and Eu) using the same stellar spectra for four stars (HD361, HD10700, HD121504, HD202206). Each group produced measurements for each of the stars using: 1) their own autonomous techniques, 2) standardized stellar parameters, 3) standardized line list, and 4) both standardized parameters and line list. We present the resulting stellar parameters, absolute abundances, and a metric of data similarity that quantifies homogeneity of the data. We conclude that standardization of some kind, particularly stellar parameters, improves the consistency between methods. However, because results did not converge as more free parameters were standardized, it is clear there are inherent issues within the techniques that need to be reconciled. Therefore, we encourage more conversation and transparency within the community such that stellar abundance determinations can be reproducible as well as accurate and precise. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: 86 pages, 5 figures, 24 tables, accepted to Astrophysical Journal Supplement

The effect of multiple particle sizes on cooling rates of chondrules produced in large-scale shocks in the solar nebula

Authors: Melissa A. Morris, Stuart J. Weidenschilling, Steven J. Desch

Abstract: Chondrules represent one of the best probes of the physical conditions and processes acting in the early solar nebula. Proposed chondrule formation models are assessed based on their ability to match the meteoritic evidence, especially experimental constraints on their thermal histories. The model most consistent with chondrule thermal histories is passage through shock waves in the solar nebula.… ▽ More Chondrules represent one of the best probes of the physical conditions and processes acting in the early solar nebula. Proposed chondrule formation models are assessed based on their ability to match the meteoritic evidence, especially experimental constraints on their thermal histories. The model most consistent with chondrule thermal histories is passage through shock waves in the solar nebula. Existing models of heating by shocks generally yield a good first-order approximation to inferred chondrule cooling rates. However, they predict prolonged heating in the pre-shock region, which would cause volatile loss and isotopic fractionation, which are not observed. These models have typically included particles of a single (large) size, i.e., chondrule precursors, or at most, large particles accompanied by micron-sized grains. The size distribution of solids present during chondrule formation controls the opacity of the affected region, and significantly affects the thermal histories of chondrules. Micron-sized grains evaporate too quickly to prevent excessive heating of chondrule precursors. However, isolated grains in chondrule-forming regions would rapidly coagulate into fractal aggregates. Pre-shock heating by infrared radiation from the shock front would cause these aggregates to melt and collapse into intermediate-sized (tens of microns) particles. We show that inclusion of such particles yields chondrule cooling rates consistent with petrologic and isotopic constraints. △ Less

Submitted 11 March, 2016; originally announced March 2016.

Comments: Accepted for publication in Meteoritics and Planetary Science

Magneto-rotational instability in the protolunar disk

Authors: Augusto Carballido, Steven J. Desch, G. Jeffrey Taylor

Abstract: (Abridged) We perform the first study of magnetohydrodynamic processes in the protolunar disk (PLD). With the use of published data on the chemical composition of the PLD, along with existing analytical models of the disk structure, we show that the high temperatures that were prevalent in the disk would have led to ionization of Na, K, SiO, Zn and, to a lesser extent, O$_2$. We assume that the di… ▽ More (Abridged) We perform the first study of magnetohydrodynamic processes in the protolunar disk (PLD). With the use of published data on the chemical composition of the PLD, along with existing analytical models of the disk structure, we show that the high temperatures that were prevalent in the disk would have led to ionization of Na, K, SiO, Zn and, to a lesser extent, O$_2$. We assume that the disk has a vapor structure. The resulting ionization fractions, together with a relatively weak magnetic field, would have been sufficient to trigger the magneto-rotational instability, or MRI. We calculate the intensity of the resulting magnetohydrodynamic turbulence, as parameterized by the dimensionless ratio $α$ of turbulent stresses to gas pressure, and obtain maximum values $α\sim 10^{-2}$ along most of the vertical extent of the disk, and at different orbital radii. Under these conditions, turbulent mixing within the PLD due to the MRI was likely capable of transporting isotopic and chemical species efficiently. To test these results in a conservative manner, we carry out a numerical magnetohydrodynamic simulation of a small, rectangular patch of the PLD, located at 4 Earth radii ($r_{\rm E}$) from the center of the Earth, and assuming once again that the disk is completely gaseous. We obtain turbulence with an average intensity $α\sim 7\times 10^{-6}$ over the course of 280 orbital periods (133 days at 4$r_{\rm{E}}$). Despite this relatively low value of $α$, the effective turbulent diffusivity of a passive tracer introduced in the flow is large enough to allow the tracer to spread across a radial distance of $10r_{\rm{E}}$ in $\sim13-129$ yr, less than the estimated cooling time of the PLD of $\sim250$ yr. Further improvements to our model will need to incorporate the energy balance in the disk, a complete two-phase structure, and a more realistic equation of state. △ Less

Submitted 19 December, 2015; originally announced December 2015.

Comments: Accepted for publication in Icarus

Report number: CASPER-15-5

External Photoevaporation of the Solar Nebula II: Effects on Disk Structure and Evolution with Non-Uniform Turbulent Viscosity due to the Magnetorotational Instability

Authors: Anusha Kalyaan, Steven J Desch, Nikhil Monga

Abstract: The structure and evolution of protoplanetary disks, especially the radial flows of gas through them, are sensitive to a number of factors. One that has been considered only occasionally in the literature is external photoevaporation by far-ultraviolet (FUV) radiation from nearby, massive stars, despite the fact that nearly half of all disks will experience photoevaporation. Another effect apparen… ▽ More The structure and evolution of protoplanetary disks, especially the radial flows of gas through them, are sensitive to a number of factors. One that has been considered only occasionally in the literature is external photoevaporation by far-ultraviolet (FUV) radiation from nearby, massive stars, despite the fact that nearly half of all disks will experience photoevaporation. Another effect apparently not considered in the literature is a spatially and temporally varying value of $α$ in the disk [where the turbulent viscosity $ν$ is $α$ times the sound speed C times the disk scale height H]. Here we use the formulation of Bai \& Stone (2011) to relate $α$ to the ionization fraction in the disk, assuming turbulent transport of angular momentum is due to the magnetorotational instability. We find that disk evolution is most sensitive to the surface area of dust. Typically $α\lesssim 10^{-5}$ in the inner disk ($< 2$ AU), rising to $\sim 10^{-1}$ beyond 20 AU. This drastically alters the structure of the disk and the flow of mass through it: while the outer disk rapidly viscously spreads, the inner disk hardly evolves; this leads to a steep surface density profile with a slope < p > $\approx$ 2 - 5 in the 5-30 AU region) that is made steeper by external photoevaporation. We also find that the combination of variable $α$ and external photoevaporation eventually causes gas as close as 3 AU, previously accreting inward, to be drawn outward to the photoevaporated outer edge of the disk. These effects have drastic consequences for planet formation and volatile transport in protoplanetary disks. △ Less

Submitted 17 November, 2015; originally announced November 2015.

Comments: 65 pages, 23 figures, Accepted for publication in the Astrophysical Journal

High-Temperature Ionization in Protoplanetary Disks

Authors: Steven J. Desch, Neal J. Turner

Abstract: We calculate the abundances of electrons and ions in the hot (> 500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show t… ▽ More We calculate the abundances of electrons and ions in the hot (> 500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show that the Saha equation often does not hold, because free charges are produced by thermionic and ion emission and destroyed when they stick to grain surfaces. This means the ionization state depends not on the first ionization potential of the alkali atoms, but rather on the grains' work functions. The charged species' abundances typically rise abruptly above about 800 K, with little qualitative dependence on the work function, gas density, or dust-to-gas mass ratio. Applying our results, we find that protoplanetary disks' dead zone, where high diffusivities stifle magnetorotational turbulence, has its inner edge located where the temperature exceeds a threshold value ~1000 K. The threshold is set by ambipolar diffusion except at the highest densities, where it is set by Ohmic resistivity. We find that the disk gas can be diffusively loaded onto the stellar magnetosphere at temperatures below a similar threshold. We investigate whether the "short-circuit" instability of current sheets can operate in disks and find that it cannot, or works only in a narrow range of conditions; it appears not to be the chondrule formation mechanism. We also suggest that thermionic emission is important for determining the rate of Ohmic heating in hot Jupiters. △ Less

Submitted 31 August, 2015; originally announced August 2015.

Comments: 68 pages, 20 figures

Short GRB 130603B: Discovery of a jet break in the optical and radio afterglows, and a mysterious late-time X-ray excess

Authors: Wen-fai Fong, Edo Berger, Brian D. Metzger, Raffaella Margutti, Ryan Chornock, Giulia Migliori, Ryan J. Foley, B. Ashley Zauderer, Ragnhild Lunnan, Tanmoy Laskar, Steven J. Desch, Karen J. Meech, Sarah Sonnett, Claire M. Dickey, Anne M. Hedlund, Paul Harding

Abstract: We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio ba… ▽ More We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio band. We model the temporal evolution of the spectral energy distribution to determine the burst explosion properties and find an isotropic-equivalent kinetic energy of (0.6-1.7) x 10^51 erg and a circumburst density of 5 x 10^-3-30 cm^-3. From the inferred opening angle of GRB 130603B, we calculate beaming-corrected energies of Egamma (0.5-2) x 10^49 erg and EK (0.1-1.6) x 10^49 erg. Along with previous measurements and lower limits we find a median short GRB opening angle of 10 deg. Using the all-sky observed rate of 10 Gpc^-3 yr^-1, this implies a true short GRB rate of 20 yr^-1 within 200 Mpc, the Advanced LIGO/VIRGO sensitivity range for neutron star binary mergers. Finally, we uncover evidence for significant excess emission in the X-ray afterglow of GRB 130603B at >1 d and conclude that the additional energy component could be due to fall-back accretion or spin-down energy from a magnetar formed following the merger. △ Less

Submitted 28 September, 2013; originally announced September 2013.

Comments: Submitted to ApJ; emulateapj style; 10 pages, 1 table, 3 figures

High-Temperature Processing of Solids Through Solar Nebular Bow Shocks: 3D Radiation Hydrodynamics Simulations with Particles

Authors: A. C. Boley, M. A. Morris, S. J. Desch

Abstract: A fundamental, unsolved problem in Solar System formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks has been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals an… ▽ More A fundamental, unsolved problem in Solar System formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks has been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through 3D radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H$_2$ is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities $\sim$few$\times10^{-8}$ L$_{\odot}$. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids. △ Less

Submitted 27 August, 2013; v1 submitted 13 August, 2013; originally announced August 2013.

Comments: Accepted for publication in ApJ. Images have been resized to conform to arXiv limits, but are all readable upon adjusting the zoom. Changes from v1: Corrected typos discovered in proofs. Most changes are in the appendix

Phyllosilicate Emission from Protoplanetary Disks: Is the Indirect Detection of Extrasolar Water Possible?

Authors: Melissa A. Morris, Steven J. Desch

Abstract: Phyllosilicates are hydrous minerals formed by interaction between rock and liquid water and are commonly found in meteorites originating in the asteroid belt. Collisions between asteroids contribute to zodiacal dust, which therefore reasonably could include phyllosilicates. Collisions between planetesimals in protoplanetary disks may also produce dust containing phyllosilicates. These minerals po… ▽ More Phyllosilicates are hydrous minerals formed by interaction between rock and liquid water and are commonly found in meteorites originating in the asteroid belt. Collisions between asteroids contribute to zodiacal dust, which therefore reasonably could include phyllosilicates. Collisions between planetesimals in protoplanetary disks may also produce dust containing phyllosilicates. These minerals possess characteristic emission features in the mid-infrared and could be detectable in extrasolar protoplanetary disks. Here we determine whether phyllosilicates in protoplanetary disks are detectable in the infrared using instruments such as those on board the Spitzer Space Telescope and SOFIA (Stratospheric Observatory for Infrared Astronomy). We calculate opacities for the phyllosilicates most common in meteorites and compute the emission of radiation from a protoplanetary disk using a 2-layer radiative transfer model. We find that phyllosilicates present at the 3% level lead to observationally significant differences in disk spectra, and should therefore be detectable using infrared observations and spectral modeling. Detection of phyllosilicates in a protoplanetary disk would be diagnostic of liquid water in planetesimals in that disk and would demonstrate similarity to our own Solar System. We also discuss use of phyllosilicate emission to test the "waterworlds" hypothesis, which proposes that liquid water in planetesimals should correlate with the inventory of short-lived radionuclides in planetary systems, especially 26Al. △ Less

Submitted 9 April, 2013; originally announced April 2013.

Journal ref: ASTROBIOLOGY, Volume 9, Number 10, 2009

Mixing of Clumpy Supernova Ejecta into Molecular Clouds

Authors: Liubin Pan, Steven J. Desch, Evan Scannapieco, F. X. Timmes

Abstract: Several lines of evidence, from isotopic analyses of meteorites to studies of the Sun's elemental and isotopic composition, indicate that the solar system was contaminated early in its evolution by ejecta from a nearby supernova (SN). Previous models have invoked SN material being injected into an extant protoplanetary disk, or isotropically expanding ejecta sweeping over a distant (>10 pc) cloud… ▽ More Several lines of evidence, from isotopic analyses of meteorites to studies of the Sun's elemental and isotopic composition, indicate that the solar system was contaminated early in its evolution by ejecta from a nearby supernova (SN). Previous models have invoked SN material being injected into an extant protoplanetary disk, or isotropically expanding ejecta sweeping over a distant (>10 pc) cloud core, simultaneously enriching it and triggering its collapse. Here we consider a new astrophysical setting: the injection of clumpy SN ejecta, as observed in the Cas A SN remnant, into the molecular gas at the periphery of an HII region created by the SN's progenitor star. To track these interactions we have conducted a suite of high-resolution (1500^3 effective) 3D simulations that follow the evolution of individual clumps as they move into molecular gas. Even at these high resolutions, our simulations do not quite achieve numerical convergence, due to the challenge of properly resolving the small-scale mixing of ejecta and molecular gas, although they do allow some robust conclusions to be drawn. Isotropically exploding ejecta do not penetrate into the molecular cloud, but, if cooling is properly accounted for, clumpy ejecta penetrate to distances ~10^18 cm and mix effectively with star-forming molecular gas. The ~2 M_\odot high-metallicity ejecta from a core-collapse SN is likely to mix with ~2 \times 10^4 M_\odot of molecular gas. Thus all stars forming late (~5 Myr) in the evolution of an HII region may be contaminated by SN ejecta at a level ~10^-4. This level of contamination is consistent with the abundances of short-lived radionuclides and possibly some stable isotopic shifts in the early solar system, and is potentially consistent with the observed variability in stellar elemental abundances. SN contamination of forming planetary systems may be a common, universal process. △ Less

Submitted 27 June, 2012; originally announced June 2012.

Comments: 21 pages, 10 figures, accepted by ApJ

Chondrule Formation in Bow Shocks around Eccentric Planetary Embryos

Authors: Melissa A. Morris, Aaron C. Boley, Steven J. Desch, Themis Athanassiadou

Abstract: Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules - millimeter-sized igneous inclusions found in meteorites - were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We n… ▽ More Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules - millimeter-sized igneous inclusions found in meteorites - were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We numerically model the hydrodynamics of the flow, and find that such large bodies retain an atmosphere, with Kelvin-Helmholtz instabilities allowing mixing of this atmosphere with the gas and particles flowing past the embryo. We calculate the trajectories of chondrules flowing past the body, and find that they are not accreted by the protoplanet, but may instead flow through volatiles outgassed from the planet's magma ocean. In contrast, chondrules are accreted onto smaller planetesimals. We calculate the thermal histories of chondrules passing through the bow shock. We find that peak temperatures and cooling rates are consistent with the formation of the dominant, porphyritic texture of most chondrules, assuming a modest enhancement above the likely solar nebula average value of chondrule densities (by a factor of 10), attributable to settling of chondrule precursors to the midplane of the disk or turbulent concentration. We calculate the rate at which a planetary embryo's eccentricity is damped and conclude that a single planetary embryo scattered into an eccentric orbit can, over ~ 10e5 years, produce ~ 10e24 g of chondrules. In principle, a small number (1-10) of eccentric planetary embryos can melt the observed mass of chondrules in a manner consistent with all known constraints. △ Less

Submitted 3 April, 2012; originally announced April 2012.

Comments: Accepted for publication in The Astrophysical Journal

A Critical Examination of the X-Wind Model for Chondrule and Calcium-rich, Aluminum-rich Inclusion Formation and Radionuclide Production

Authors: S. J. Desch, M. A. Morris, H. C. Connolly, Jr., Alan P. Boss

Abstract: Meteoritic data, especially regarding chondrules and calcium-rich, aluminum-rich inclusions (CAIs), and isotopic evidence for short-lived radionuclides (SLRs) in the solar nebula, potentially can constrain how planetary systems form. Intepretation of these data demands an astrophysical model, and the "X-wind" model of Shu et al. (1996) and collaborators has been advanced to explain the origin of c… ▽ More Meteoritic data, especially regarding chondrules and calcium-rich, aluminum-rich inclusions (CAIs), and isotopic evidence for short-lived radionuclides (SLRs) in the solar nebula, potentially can constrain how planetary systems form. Intepretation of these data demands an astrophysical model, and the "X-wind" model of Shu et al. (1996) and collaborators has been advanced to explain the origin of chondrules, CAIs and SLRs. It posits that chondrules and CAIs were thermally processed < 0.1 AU from the protostar, then flung by a magnetocentrifugal outflow to the 2-3 AU region to be incorporated into chondrites. Here we critically examine key assumptions and predictions of the X-wind model. We find a number of internal inconsistencies: theory and observation show no solid material exists at 0.1 AU; particles at 0.1 AU cannot escape being accreted into the star; particles at 0.1 AU will collide at speeds high enough to destroy them; thermal sputtering will prevent growth of particles; and launching of particles in magnetocentrifugal outflows is not modeled, and may not be possible. We also identify a number of incorrect predictions of the X-wind model: the oxygen fugacity where CAIs form is orders of magnitude too oxidizing; chondrule cooling rates are orders of magnitude lower than those experienced by barred olivine chondrules; chondrule-matrix complementarity is not predicted; and the SLRs are not produced in their observed proportions. We conclude that the X-wind model is not relevant to chondrule and CAI formation and SLR production. We discuss more plausible models for chondrule and CAI formation and SLR production. △ Less

Submitted 15 November, 2010; originally announced November 2010.

Comments: Accepted for publication in The Astrophysical Journal

Collateral Effects on Solar Nebula Oxygen Isotopes due to Injection of 26Al by a Nearby Supernova

Authors: C. I. Ellinger, P. A. Young, S. J. Desch

Abstract: Injection of material from a core-collapse supernova into the solar system's already-formed disk is one proposed mechanism for producing the short-lived radionuclides, such as 26Al and 41Ca, inferred from isotopic studies of meteorites to have existed in the solar nebula. This hypothesis has recently been challenged on the basis that the injection of enough supernova material to match the meteorit… ▽ More Injection of material from a core-collapse supernova into the solar system's already-formed disk is one proposed mechanism for producing the short-lived radionuclides, such as 26Al and 41Ca, inferred from isotopic studies of meteorites to have existed in the solar nebula. This hypothesis has recently been challenged on the basis that the injection of enough supernova material to match the meteoritic abundances of 26Al and 41Ca would produce large, measureable, and unobserved collateral effects on oxygen isotopes. Here we calculate again the shifts in oxygen isotopes due to injection of supernova material in the solar nebula, using a variety of nucleosynthetic conditions of our own progenitor explosions. Unlike previous studies of this type, we also consider the effect of non-homogeneity in abundance distribution of the nucleosynthesis products after the explosion. We calculate the shifts in oxygen isotopes due to injection of sufficient supernova material to produce the meteoritic abundances of 26Al and 41Ca, and analyze the predicted shifts in detail for compatibility with meteoritic data. We find that the range in possible isotopic shifts is considerable and sensitive to parameters such as progenitor mass and anisotropy of the explosion; however, a small number of compatible scenarios do exist. Because of the wide range of outcomes and the sensitivity of isotopic yields to assumed conditions, it is difficult to constrain the supernova that may have led to injection of 26Al in the solar nebula. Conversely, we argue that the existence of viable counterexamples demonstrates that it is premature to use oxygen isotopes to rule out injection of 26Al and 41Ca into the solar nebula protoplanetary disk by a nearby supernova. △ Less

Submitted 22 October, 2010; originally announced October 2010.

Comments: 35 pages, 4 figures, 5 tables

Journal ref: Ap.J. 725 (2010) 1495

Thermal Histories of Chondrules in Solar Nebula Shocks

Authors: M. A. Morris, S. J. Desch

Abstract: Chondrules are important early Solar System materials that can provide a wealth of information on conditions in the solar nebula, if their formation mechanism can be understood. The theory most consistent with observational constraints, especially thermal histories, is the so-called "shock model", in which chondrules were melted in solar nebula shocks. However, several problems have been identifie… ▽ More Chondrules are important early Solar System materials that can provide a wealth of information on conditions in the solar nebula, if their formation mechanism can be understood. The theory most consistent with observational constraints, especially thermal histories, is the so-called "shock model", in which chondrules were melted in solar nebula shocks. However, several problems have been identified with previous shock models. These problems all pertained to the treatment of the radiation field, namely the input boundary condition to the radiation field, the proper treatment of the opacity of solids, and the proper treatment of molecular line cooling. In this paper, we present the results of our updated shock model, which corrects for the problems listed above. Our new hydrodynamic shock code includes a complete treatment of molecular line cooling due to H2O. Previously, shock models including line cooling predicted chondrule cooling rates exceeding 100,000 K/hr. Contrary to these expectations, we have found that the effect of line cooling is minimal; after the inclusion of line cooling, the cooling rates of chondrules are 10-1000 K/hr. The reduction in the otherwise rapid cooling rates attributable to line cooling is due to a combination of factors, including buffering due to hydrogen recombination/dissociation, high column densities of water, and backwarming. Our model demonstrates that the shock model for chondrule formation remains consistent with observational constraints. △ Less

Submitted 16 August, 2010; originally announced August 2010.

Comments: Accepted for publication in The Astrophysical Journal, a journal of The American Astronomical Society published by IOP Publishing Ltd which trades as IOP Publishing Ltd of Dirac House, Temple Back, Bristol BS1 6BE, UK

Cooling of Dense Gas by H2O Line Emission and an Assessment of its Effects in Chondrule-Forming Shocks

Authors: M. A. Morris, S. J. Desch, F. J. Ciesla

Abstract: We consider gas at densities appropriate to protoplanetary disks and calculate its ability to cool due to line radiation emitted by H2O molecules within the gas. Our work follows that of Neufeld & Kaufman (1993; ApJ, 418, 263), expanding on their work in several key aspects, including use of a much expanded line database, an improved escape probability formulism, and the inclusion of dust grains… ▽ More We consider gas at densities appropriate to protoplanetary disks and calculate its ability to cool due to line radiation emitted by H2O molecules within the gas. Our work follows that of Neufeld & Kaufman (1993; ApJ, 418, 263), expanding on their work in several key aspects, including use of a much expanded line database, an improved escape probability formulism, and the inclusion of dust grains, which can absorb line photons. Although the escape probabilities formally depend on a complicated combination of optical depth in the lines and in the dust grains, we show that the cooling rate including dust is well approximated by the dust-free cooling rate multiplied by a simple function of the dust optical depth. We apply the resultant cooling rate of a dust-gas mixture to the case of a solar nebula shock pertinent to the formation of chondrules, millimeter-sized melt droplets found in meteorites. Our aim is to assess whether line cooling can be neglected in chondrule-forming shocks or if it must be included. We find that for typical parameters, H2O line cooling shuts off a few minutes past the shock front; line photons that might otherwise escape the shocked region and cool the gas will be absorbed by dust grains. During the first minute or so past the shock, however, line photons will cool the gas at rates ~ 10,000 K/hr, dropping the temperature of the gas (and most likely the chondrules within the gas) by several hundred K. Inclusion of H2O line cooling therefore must be included in models of chondrule formation by nebular shocks. △ Less

Submitted 29 September, 2008; originally announced September 2008.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: Astrophysical Journal 691 (2009) 320-331

Spitzer Observations of the HII Region NGC 2467: An Analysis of Triggered Star Formation

Authors: Keely D. Snider, J. Jeff Hester, Steven J. Desch, Kevin R. Healy, Jon Bally

Abstract: We present new Spitzer Space Telescope observations of the region NGC 2467, and use these observations to determine how the environment of an HII region affects the process of star formation. Our observations comprise IRAC (3.6, 4.5, 5.8, and 8.0 um) and MIPS (24 um) maps of the region, covering approximately 400 square arcminutes. The images show a region of ionized gas pushing out into the sur… ▽ More We present new Spitzer Space Telescope observations of the region NGC 2467, and use these observations to determine how the environment of an HII region affects the process of star formation. Our observations comprise IRAC (3.6, 4.5, 5.8, and 8.0 um) and MIPS (24 um) maps of the region, covering approximately 400 square arcminutes. The images show a region of ionized gas pushing out into the surrounding molecular cloud, powered by an O6V star and two clusters of massive stars in the region. We have identified as candidate Young Stellar Objects (YSOs) 45 sources in NGC 2467 with infrared excesses in at least two mid-infrared colors. We have constructed color-color diagrams of these sources and have quantified their spatial distribution within the region. We find that the YSOs are not randomly distributed in NGC 2467; rather, over 75% of the sources are distributed at the edge of the HII region, along ionization fronts driven by the nearby massive stars. The high fraction of YSOs in NGC 2467 that are found in proximity to gas that has been compressed by ionization fronts supports the hypothesis that a significant fraction of the star formation in NGC 2467 is triggered by the massive stars and the expansion of the HII region. At the current rate of star formation, we estimate at least 25-50% of the total population of YSOs formed by this process. △ Less

Submitted 23 June, 2009; v1 submitted 9 November, 2007; originally announced November 2007.

Comments: Accepted to the Astrophysical Journal, set to appear in Volume 701; 18 pages, 13 figures, 4 tables. This version reflects a few major changes made in the accepted version, including new figures

Journal ref: Astrophys.J.700:506-522,2009

Interaction of Supernova Ejecta with Nearby Protoplanetary Disks

Authors: N. Ouellette, S. J. Desch, J. J. Hester

Abstract: The early Solar System contained short-lived radionuclides such as 60Fe (t1/2 = 1.5 Myr) whose most likely source was a nearby supernova. Previous models of Solar System formation considered a supernova shock that triggered the collapse of the Sun's nascent molecular cloud. We advocate an alternative hypothesis, that the Solar System's protoplanetary disk had already formed when a very close (<… ▽ More The early Solar System contained short-lived radionuclides such as 60Fe (t1/2 = 1.5 Myr) whose most likely source was a nearby supernova. Previous models of Solar System formation considered a supernova shock that triggered the collapse of the Sun's nascent molecular cloud. We advocate an alternative hypothesis, that the Solar System's protoplanetary disk had already formed when a very close (< 1 pc) supernova injected radioactive material directly into the disk. We conduct the first numerical simulations designed to answer two questions related to this hypothesis: will the disk be destroyed by such a close supernova; and will any of the ejecta be mixed into the disk? Our simulations demonstrate that the disk does not absorb enough momentum from the shock to escape the protostar to which it is bound. Only low amounts (< 1%) of mass loss occur, due to stripping by Kelvin-Helmholtz instabilities across the top of the disk, which also mix into the disk about 1% of the intercepted ejecta. These low efficiencies of destruction and injectation are due to the fact that the high disk pressures prevent the ejecta from penetrating far into the disk before stalling. Injection of gas-phase ejecta is too inefficient to be consistent with the abundances of radionuclides inferred from meteorites. On the other hand, the radionuclides found in meteorites would have condensed into dust grains in the supernova ejecta, and we argue that such grains will be injected directly into the disk with nearly 100% efficiency. The meteoritic abundances of the short-lived radionuclides such as 60Fe therefore are consistent with injection of grains condensed from the ejecta of a nearby (< 1 pc) supernova, into an already-formed protoplanetary disk. △ Less

Submitted 12 April, 2007; originally announced April 2007.

Comments: 57 pages, 16 figures

Journal ref: Astrophys.J.662:1268-1281,2007

Understanding our Origins: Star Formation in H II Region Environments

Authors: J. Jeff Hester, Steven J. Desch

Abstract: Recent analysis of the decay products of short-lived radiounclides (SLRs) in meteorites, in particular the confirmation of the presence of live 60Fe in the early Solar System, provides unambiguous evidence that the Sun and Solar System formed near a massive star. We consider the question of the formation of low-mass stars in environments near massive stars, presenting a scenario for the evolutio… ▽ More Recent analysis of the decay products of short-lived radiounclides (SLRs) in meteorites, in particular the confirmation of the presence of live 60Fe in the early Solar System, provides unambiguous evidence that the Sun and Solar System formed near a massive star. We consider the question of the formation of low-mass stars in environments near massive stars, presenting a scenario for the evolution of a star and its disk around the periphery of an expanding H II region. The stages in this scenario begin with compression of molecular gas around the edge of an H II region, continue as forming stars are overrun by the advancing ionization front, and culminate when ejecta from one or more nearby supernova explosions sweeps over YSO disks located in the low density interior of the H II region, injecting SLRs including 26Al and 60Fe. We review the evidence that this mode of star formation is more characteristic of formation of low-mass stars than is the mode of star formation seen in regions such as the Taurus-Auriga molecular cloud. We discuss the implications of this scenario for our understanding of star formation, as well as the effects of the young Sun's astrophysical environment on the formation and evolution of the Solar System. We conclude that low-mass stars and their accompanying disks form and evolve very differently near massive stars than they do in regions like Taurus-Auriga, and that these differences have profound implications for our understanding of our origins. △ Less

Submitted 8 June, 2005; originally announced June 2005.

Comments: 25 pages, 10 figures. Refereed and accepted submission of an invited presentation to the Proceedings from the Workshop on Chondrites and the Protoplanetary Disk, Kaua'i Hawai'i, November 8-11, 2004. ASP conference series, editors A. Krot, E. Scott, & B. Reipurth. Higher quality copy available at http://eagle.la.asu.edu/hester/kauai/kauai_hester_final.pdf

Annealing of Silicate Dust by Nebular Shocks at 10 AU

Authors: David E. Harker, Steven J. Desch

Abstract: Silicate dust grains in the interstellar medium are known to be mostly amorphous, yet crystalline silicate grains have been observed in many long-period comets and in protoplanetary disks. Annealing of amorphous silicate grains into crystalline grains requires temperatures > 1000 K, but exposure of dust grains in comets to such high temperatures is incompatible with the generally low temperature… ▽ More Silicate dust grains in the interstellar medium are known to be mostly amorphous, yet crystalline silicate grains have been observed in many long-period comets and in protoplanetary disks. Annealing of amorphous silicate grains into crystalline grains requires temperatures > 1000 K, but exposure of dust grains in comets to such high temperatures is incompatible with the generally low temperatures experienced by comets. This has led to the proposal of models in which dust grains were thermally processed near the protoSun, then underwent considerable radial transport until they reached the gas giant planet region where the long-period comets originated. We hypothesize instead that silicate dust grains were annealed in situ, by shock waves triggered by gravitational instabilities. We assume a shock speed of 5 km/s, a plausible value for shocks driven by gravitational instabilities. We calculate the peak temperatures of micron and submicron amorphous pyroxene grains of chondritic composition under conditions typical in protoplanetary disks at 5 - 10 AU. Our results also apply to chondritic amorphous olivine grains. We show that {\it in situ} thermal annealing of submicron and micron-sized silicate dust grains can occur, obviating the need for large-scale radial transport. △ Less

Submitted 20 December, 2001; originally announced December 2001.

Comments: 12 pages; includes 1 figure, 1 table; accepted by ApJ Letters

Ambipolar Diffusion and Far-Infrared Polarization from the Galactic Circumnuclear Disk

Authors: S. J. Desch, W. G. Roberge

Abstract: We describe an implicit prediction of the accretion disk models constructed by Wardle and Konigl (1990) for the circumnuclear disk (CND) of gas and dust near the Galactic center: supersonic ambipolar diffusion, an essential dynamical ingredient of the Wardle-Konigl disks, will cause the alignment of dust grains due to a process described by Roberge, Hanany, & Messinger (1995). We calculate synth… ▽ More We describe an implicit prediction of the accretion disk models constructed by Wardle and Konigl (1990) for the circumnuclear disk (CND) of gas and dust near the Galactic center: supersonic ambipolar diffusion, an essential dynamical ingredient of the Wardle-Konigl disks, will cause the alignment of dust grains due to a process described by Roberge, Hanany, & Messinger (1995). We calculate synthetic maps of the polarized thermal emission which would be caused by ambipolar alignment in the preferred Wardle-Konigl model. Our maps are in reasonable agreement with 100 micron polarimetry of the CND if we assume that the grains have shapes similar to those of grains in nearby molecular clouds and that the CND contains a disordered magnetic field in energy equipartition with its ordered field. △ Less

Submitted 20 November, 1996; originally announced November 1996.

Comments: 14 pages (LaTeX+2 eps figs); accepted by ApJ Letters 11/19/96