Showing 1–50 of 95 results for author: Jacobson, S

Commissioning results from the Robo-AO-2 facility for rapid visible and near-infrared AO imaging

Authors: Christoph Baranec, James Ou, Reed Riddle, Ruihan Zhang, Luke Mckay, Rachel Rampy, Morgan Bonnet, Iven Hamilton, Greg Ching, Jessica Young, Maıssa Salama, Paul Barnes, Shane Jacobson, Peter Onaka, Mark Chun, Zachary Werber, Keith Powell, Marcos A. van Dam, Benjamin Shappee

Abstract: We installed the next-generation automated laser adaptive optics system, Robo-AO-2, on the University of Hawaii 2.2-m telescope on Maunakea in 2023. We engineered Robo-AO-2 to deliver robotic, diffraction-limited observations at visible and near-infrared wavelengths in unprecedented numbers. This new instrument takes advantage of upgraded components, manufacturing techniques and control; and inclu… ▽ More We installed the next-generation automated laser adaptive optics system, Robo-AO-2, on the University of Hawaii 2.2-m telescope on Maunakea in 2023. We engineered Robo-AO-2 to deliver robotic, diffraction-limited observations at visible and near-infrared wavelengths in unprecedented numbers. This new instrument takes advantage of upgraded components, manufacturing techniques and control; and includes a parallel reconfigurable natural guide star wavefront sensor with which to explore hybrid wavefront sensing techniques. We present the results of commissioning in 2023 and 2024. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Comments: Proceedings of SPIE, 13097-15

BYORP and Dissipation in Binary Asteroids: Lessons from DART

Authors: Matija Ćuk, Harrison Agrusa, Rachel H. Cueva, Fabio Ferrari, Masatoshi Hirabayashi, Seth A. Jacobson, Jay McMahon, Patrick Michel, Paul Sánchez, Daniel J. Scheeres, Stephen Schwartz, Kevin J. Walsh, Yun Zhang

Abstract: The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary.… ▽ More The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary. DART revealed a largely oblate and ellipsoidal shape of Dimorphos before the impact, while the post-impact observations suggest that Dimorphos now has a prolate shape. Here we add those data points to the known properties of small binary asteroids and propose new paradigms of the radiative binary YORP (BYORP) effect as well as tidal dissipation in small binaries. We find that relatively spheroidal bodies like Dimorphos made of small debris may experience a weaker and more size-dependent BYORP effect than previously thought. This could explain the observed values of period drift in several well-characterized binaries. We also propose that energy dissipation in small binaries is dominated by relatively brief episodes of large-scale movement of (likely surface) materials, rather than long-term steady-state tidal dissipation. We propose that one such episode was triggered on Dimorphos by the DART impact. Depending on the longevity of this high-dissipation regime, it is possible that Dimorphos will be more dynamically relaxed in time for the Hera mission than it was in the weeks following the impact. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: Accepted for PSJ

Gliese 12 b: A temperate Earth-sized planet at 12 pc ideal for atmospheric transmission spectroscopy

Authors: M. Kuzuhara, A. Fukui, J. H. Livingston, J. A. Caballero, J. P. de Leon, T. Hirano, Y. Kasagi, F. Murgas, N. Narita, M. Omiya, Jaume Orell-Miquel, E. Palle, Q. Changeat, E. Esparza-Borges, H. Harakawa, C. Hellier, Yasunori Hori, Kai Ikuta, H. T. Ishikawa, T. Kodama, T. Kotani, T. Kudo, J. C. Morales, M. Mori, E. Nagel , et al. (81 additional authors not shown)

Abstract: Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We repor… ▽ More Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We report the discovery of an Earth-sized planet transiting the nearby (12 pc) inactive M3.0 dwarf Gliese 12 (TOI-6251) with an orbital period ($P_{\rm{orb}}$) of 12.76 days. The planet, Gliese 12b, was initially identified as a candidate with an ambiguous $P_{\rm{orb}}$ from TESS data. We confirmed the transit signal and $P_{\rm{orb}}$ using ground-based photometry with MuSCAT2 and MuSCAT3, and validated the planetary nature of the signal using high-resolution images from Gemini/NIRI and Keck/NIRC2 as well as radial velocity (RV) measurements from the InfraRed Doppler instrument on the Subaru 8.2 m telescope and from CARMENES on the CAHA 3.5 m telescope. X-ray observations with XMM-Newton showed the host star is inactive, with an X-ray-to-bolometric luminosity ratio of $\log L_{\rm X}/L_{\rm bol} \approx -5.7$. Joint analysis of the light curves and RV measurements revealed that Gliese 12b has a radius of 0.96 $\pm$ 0.05 $R_\oplus$, a 3$σ$ mass upper limit of 3.9 $M_\oplus$, and an equilibrium temperature of 315 $\pm$ 6 K assuming zero albedo. The transmission spectroscopy metric (TSM) value of Gliese 12b is close to the TSM values of the TRAPPIST-1 planets, adding Gliese 12b to the small list of potentially terrestrial, temperate planets amenable to atmospheric characterization with JWST. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 29 pages (20 pages in main body), 13 figures (10 figures in main body). Equal contributions from M. K. and A. F.. Accepted for Publication in ApJL at 2024 March 21

Journal ref: Published on 2024 May 23 by Astrophysical Journal Letters (ApJL) 967 L21

Forming Mercury from excited initial conditions

Authors: Jennifer Scora, Diana Valencia, Alessandro Morbidelli, Seth Jacobson

Abstract: Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found… ▽ More Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found more success forming Mercury analogues by adding additional embryos near Mercury's orbit. In this work, we aim to form Mercury by simulating the formation of the solar system in the presence of the giant planets Jupiter and Saturn. We test out the effect of an inner disk of embryos added on to the commonly-used narrow annulus of initial material. We form Mercury analogues with core-mass fractions (CMF) $> 0.4$ in $\sim 10\%$ of our simulations, and twice that number of Mercury analogues form during the formation process, but are unstable and do not last to the end of the simulations. Mercury analogues form at similar rates for both disks with and without an inner component, and most of our Mercury analogues have lower CMF than that of Mercury, $\sim 0.7$, due to significant accretion of debris material. We suggest that a more in-depth understanding of the fraction of debris mass that is lost to collisional grinding is necessary to understand Mercury's formation, or some additional mechanism is required to stop this debris from accreting. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 15 pages, accepted for publication in ApJ

Direct $N$-body simulations of satellite formation around small asteroids: insights from DART's encounter with the Didymos system

Authors: Harrison F. Agrusa, Yun Zhang, Derek C. Richardson, Petr Pravec, Matija Ćuk, Patrick Michel, Ronald-Louis Ballouz, Seth A. Jacobson, Daniel J. Scheeres, Kevin Walsh, Olivier Barnouin, R. Terik Daly, Eric Palmer, Maurizio Pajola, Alice Lucchetti, Filippo Tusberti, Joseph V. DeMartini, Fabio Ferrari, Alex J. Meyer, Sabina D. Raducan, Paul Sánchez

Abstract: We explore binary asteroid formation by spin-up and rotational disruption considering the NASA DART mission's encounter with the Didymos-Dimorphos binary, which was the first small binary visited by a spacecraft. Using a suite of $N$-body simulations, we follow the gravitational accumulation of a satellite from meter-sized particles following a mass-shedding event from a rapidly rotating primary.… ▽ More We explore binary asteroid formation by spin-up and rotational disruption considering the NASA DART mission's encounter with the Didymos-Dimorphos binary, which was the first small binary visited by a spacecraft. Using a suite of $N$-body simulations, we follow the gravitational accumulation of a satellite from meter-sized particles following a mass-shedding event from a rapidly rotating primary. The satellite's formation is chaotic, as it undergoes a series of collisions, mergers, and close gravitational encounters with other moonlets, leading to a wide range of outcomes in terms of the satellite's mass, shape, orbit, and rotation state. We find that a Dimorphos-like satellite can form rapidly, in a matter of days, following a realistic mass-shedding event in which only ${\sim}2-3\%$ of the primary's mass is shed. Satellites can form in synchronous rotation due to their formation near the Roche limit. There is a strong preference for forming prolate (elongated) satellites, although some simulations result in oblate spheroids like Dimorphos. The distribution of simulated secondary shapes is broadly consistent with other binary systems, measured through radar or lightcurves. Unless Dimorphos's shape is an outlier, and considering the observational bias against lightcurve-based determination of secondary elongations for oblate bodies, we suggest there could be a significant population of oblate secondaries. If these satellites initially form with elongated shapes, a yet-unidentified pathway is needed to explain how they become oblate. Finally, we show that this chaotic formation pathway occasionally forms asteroid pairs and stable triples, including co-orbital satellites and satellites in mean motion resonances. △ Less

Submitted 18 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 41 pages, 23 figures, accepted to PSJ, movies available at https://zenodo.org/doi/10.5281/zenodo.8387043

Ejecta Evolution Following a Planned Impact into an Asteroid: The First Five Weeks

Authors: Theodore Kareta, Cristina Thomas, Jian-Yang Li, Matthew M. Knight, Nicholas Moskovitz, Agata Rozek, Michele T. Bannister, Simone Ieva, Colin Snodgrass, Petr Pravec, Eileen V. Ryan, William H. Ryan, Eugene G. Fahnestock, Andrew S. Rivkin, Nancy Chabot, Alan Fitzsimmons, David Osip, Tim Lister, Gal Sarid, Masatoshi Hirabayashi, Tony Farnham, Gonzalo Tancredi, Patrick Michel, Richard Wainscoat, Rob Weryk , et al. (63 additional authors not shown)

Abstract: The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consis… ▽ More The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consistent dimming rates of 0.11-0.12 magnitudes/day in the first week, and 0.08-0.09 magnitudes/day over the entire study period. The system returned to its pre-impact brightness 24.3-25.3 days after impact through the primary ejecta tail remained. The dimming paused briefly eight days after impact, near in time to the appearance of the second tail. This was likely due to a secondary release of material after re-impact of a boulder released in the initial impact, through movement of the primary ejecta through the aperture likely played a role. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 16 pages, 5 Figures, accepted in the Astrophysical Journal Letters (ApJL) on October 16, 2023

An improved model of metal/silicate differentiation during Earth's accretion

Authors: K. I. Dale, D. C. Rubie, M. Nakajima, S. Jacobson, G. Nathan, G. J. Golabek, S. Cambioni, A. Morbidelli

Abstract: We improved the algorithm presented in Rubie et al. (2015) to model the chemical evolution of Earth driven by iron/silicate differentiation during the planet's accretion. The pressure at which the equilibration occurs during a giant impact is no longer a free parameter but is determined by the smooth particle hydrodynamic (SPH) simulations of Nakajima et al. (2021). Moreover, impacting planetesima… ▽ More We improved the algorithm presented in Rubie et al. (2015) to model the chemical evolution of Earth driven by iron/silicate differentiation during the planet's accretion. The pressure at which the equilibration occurs during a giant impact is no longer a free parameter but is determined by the smooth particle hydrodynamic (SPH) simulations of Nakajima et al. (2021). Moreover, impacting planetesimals are now assumed to be too small to cause melting and differentiation and thus their materials are stored in the crystalline upper mantle of the growing planet until a hydrostatically relaxed global magma ocean forms in the aftermath of a giant impact, whose depth is also estimated from Nakajima et al. (2021). With these changes, not all dynamical simulations lead to a satisfactory reproduction of the chemical composition of the bulk silicate Earth (BSE). Thus, the latter becomes diagnostic of the success of dynamical models. In the successful cases also the BSE abundances of W and Mo can be reproduced, that were previously hard to fit (Jennings et al., 2021). △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: Accepted in Icarus

Chemical Equilibrium Calculations for Bulk Silicate Earth Material at High Temperatures

Authors: Bruce Fegley Jr., Katharina Lodders, Nathan S. Jacobson

Abstract: The chemical equilibrium distribution of 69 elements between gas and melt is modeled for bulk silicate Earth (BSE) material from 1000 - 4500 K and 1e-6 to 100 bar. The BSE melt is modeled as a non-ideal solution and the effects of different activity coefficients and ideal solution are studied. Results include 50% condensation temperatures, major gases of each element, and oxygen fugacity (fO2) of… ▽ More The chemical equilibrium distribution of 69 elements between gas and melt is modeled for bulk silicate Earth (BSE) material from 1000 - 4500 K and 1e-6 to 100 bar. The BSE melt is modeled as a non-ideal solution and the effects of different activity coefficients and ideal solution are studied. Results include 50% condensation temperatures, major gases of each element, and oxygen fugacity (fO2) of dry and wet BSE material. The dry BSE model excludes H, C, N, F, Cl, Br, I, S, Se, Te. The wet BSE model includes H and the other volatiles. Key conclusions are much higher condensation temperatures in silicate vapor than in solar composition gas at the same total P, a different condensation sequence in silicate vapor than in solar composition gas, good agreement between different activity coefficient models except for the alkalis, agreement, where overlap exists, with prior published work, condensation of Re, Mo, W, Ru, Os oxides instead of metals, a stability field for Ni-rich metal as reported by Lock et al. (2018), agreement between ideal solution (from this work and from Lock et al. 2018) and real solution condensation temperatures for elements with minor deviations from ideality in the oxide melt, similar 50% condensation temperatures, within a few degrees, in the dry and wet BSE models for the major elements Al, Ca, Fe, Mg, Si, and the minor elements Co, Cr, Li, Mn, Ti, V, and much lower 50 percent condensation temperatures for elements such as B, Cu, K, Na, Pb, Rb, which form halide, hydroxide, sulfide, selenide, telluride and oxyhalide gases. The latter results are preliminary because the poorly known solubilities and activities of volatile elements in silicate melts must be considered for the correct equilibrium distribution, condensation temperatures and mass balance of F, Cl, Br, I, H, S, Se and Te bearing species between melt and vapor (abridged). △ Less

Submitted 16 May, 2023; originally announced May 2023.

A spectroscopic thermometer: individual vibrational band spectroscopy with the example of OH in the atmosphere of WASP-33b

Authors: Sam O. M. Wright, Stevanus K. Nugroho, Matteo Brogi, Neale P. Gibson, Ernst J. W. de Mooij, Ingo Waldmann, Jonathan Tennyson, Hajime Kawahara, Masayuki Kuzuhara, Teruyuki Hirano, Takayuki Kotani, Yui Kawashima, Kento Masuda, Jayne L. Birkby, Chris A. Watson, Motohide Tamura, Konstanze Zwintz, Hiroki Harakawa, Tomoyuki Kudo, Klaus Hodapp, Shane Jacobson, Mihoko Konishi, Takashi Kurokawa, Jun Nishikawa, Masashi Omiya , et al. (4 additional authors not shown)

Abstract: Individual vibrational band spectroscopy presents an opportunity to examine exoplanet atmospheres in detail by distinguishing where the vibrational state populations of molecules differ from the current assumption of a Boltzmann distribution. Here, retrieving vibrational bands of OH in exoplanet atmospheres is explored using the hot Jupiter WASP-33b as an example. We simulate low-resolution spectr… ▽ More Individual vibrational band spectroscopy presents an opportunity to examine exoplanet atmospheres in detail by distinguishing where the vibrational state populations of molecules differ from the current assumption of a Boltzmann distribution. Here, retrieving vibrational bands of OH in exoplanet atmospheres is explored using the hot Jupiter WASP-33b as an example. We simulate low-resolution spectroscopic data for observations with the JWST's NIRSpec instrument and use high resolution observational data obtained from the Subaru InfraRed Doppler instrument (IRD). Vibrational band-specific OH cross section sets are constructed and used in retrievals on the (simulated) low and (real) high resolution data. Low resolution observations are simulated for two WASP-33b emission scenarios: under the assumption of local thermal equilibrium (LTE) and a toy non-LTE model for vibrational excitation of selected bands. We show that mixing ratios for individual bands can be retrieved with sufficient precision to allow the vibrational population distributions of the forward models to be reconstructed. A simple fit for the Boltzmann distribution in the LTE case shows that the vibrational temperature is recoverable in this manner. For high resolution, cross-correlation applications, we apply the individual vibrational band analysis to an IRD spectrum of WASP-33b, applying an 'un-peeling' technique. Individual detection significances for the two strongest bands are shown to be in line with Boltzmann distributed vibrational state populations consistent with the effective temperature of the WASP-33b atmosphere reported previously. We show the viability of this approach for analysing the individual vibrational state populations behind observed and simulated spectra including reconstructing state population distributions. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: Submitted for publication in AJ

Constraining the Origin of Mars via Simulations of Multi-Stage Core Formation

Authors: Gabriel Nathan, David C. Rubie, Seth A. Jacobson

Abstract: It remains an elusive goal to simultaneously model the astrophysics of Solar System accretion while reproducing the mantle chemistry of more than one inner terrestrial planet. Here, we used a multistage core-mantle differentiation model based on Rubie et al. (2011,2015) to track the formation and composition of Earth and Mars in various Grand Tack formation simulations. Prior studies showed that i… ▽ More It remains an elusive goal to simultaneously model the astrophysics of Solar System accretion while reproducing the mantle chemistry of more than one inner terrestrial planet. Here, we used a multistage core-mantle differentiation model based on Rubie et al. (2011,2015) to track the formation and composition of Earth and Mars in various Grand Tack formation simulations. Prior studies showed that in order to recreate Earth's mantle composition, it must grow first from reduced (Fe-metal rich and O-poor) building blocks and then from increasingly oxidized (FeO rich) material. This accretion chemistry occurs when an oxidation gradient exists across the disk so that the innermost solids are reduced and increasingly oxidized material is found at greater heliocentric distances. For a suite of Grand Tack simulations, we investigated whether Earth and Mars can be simultaneously produced by the same oxidation gradient. Our model did not find an oxidation gradient that simultaneously reproduces the mantle composition of Earth and Mars. Due to its small mass and rapid formation, the formation history of Mars-like planets is very stochastic which decreases the likelihood of compatibility with an Earth-producing oxidation gradient in any given realization. To reconcile the accretion history and ideal chemistry of the Mars-like planet with the oxidation gradient of an Earth-producing disk, we determined where in the Earth-producing disk Mars must have formed. We find that the FeO-rich composition of the Martian mantle requires that Mars' building blocks must originate exterior to 1.0 astronomical units (AU). △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 20 pages, 7 figures, Accepted for publication at ICARUS 04/2023

Smartphone screens as astrometric calibrators

Authors: Aidan Walk, Charles-Antoine Claveau, Michael Bottom, Mark Chun, Shane Jacobson, Maxwell Service, Jessica R. Lu

Abstract: Geometric optical distortion is a significant contributor to the astrometric error budget in large telescopes using adaptive optics. To increase astrometric precision, optical distortion calibration is necessary. We investigate using smartphone OLED screens as astrometric calibrators. Smartphones are low cost, have stable illumination, and can be quickly reconfigured to probe different spatial fre… ▽ More Geometric optical distortion is a significant contributor to the astrometric error budget in large telescopes using adaptive optics. To increase astrometric precision, optical distortion calibration is necessary. We investigate using smartphone OLED screens as astrometric calibrators. Smartphones are low cost, have stable illumination, and can be quickly reconfigured to probe different spatial frequencies of an optical system's geometric distortion. In this work, we characterize the astrometric accuracy of a Samsung S20 smartphone, with a view towards providing large format, flexible astrometric calibrators for the next generation of astronomical instruments. We find the placement error of the pixels to be 189 nm +/- 15 nm RMS. At this level of error, milliarcsecond astrometric accuracy can be obtained on modern astronomical instruments. △ Less

Submitted 6 July, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: 15 pages, 11 figures; accepted, Journal of Astronomical Instrumentation

Momentum Transfer from the DART Mission Kinetic Impact on Asteroid Dimorphos

Authors: Andrew F. Cheng, Harrison F. Agrusa, Brent W. Barbee, Alex J. Meyer, Tony L. Farnham, Sabina D. Raducan, Derek C. Richardson, Elisabetta Dotto, Angelo Zinzi, Vincenzo Della Corte, Thomas S. Statler, Steven Chesley, Shantanu P. Naidu, Masatoshi Hirabayashi, Jian-Yang Li, Siegfried Eggl, Olivier S. Barnouin, Nancy L. Chabot, Sidney Chocron, Gareth S. Collins, R. Terik Daly, Thomas M. Davison, Mallory E. DeCoster, Carolyn M. Ernst, Fabio Ferrari , et al. (44 additional authors not shown)

Abstract: The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of… ▽ More The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of asteroid deflection. Here we report the first determination of the momentum transferred to an asteroid by kinetic impact. Based on the change in the binary orbit period, we find an instantaneous reduction in Dimorphos's along-track orbital velocity component of 2.70 +/- 0.10 mm/s, indicating enhanced momentum transfer due to recoil from ejecta streams produced by the impact. For a Dimorphos bulk density range of 1,500 to 3,300 kg/m$^3$, we find that the expected value of the momentum enhancement factor, $β$, ranges between 2.2 and 4.9, depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have equal densities of 2,400 kg/m$^3$, $β$= 3.61 +0.19/-0.25 (1 $σ$). These $β$ values indicate that significantly more momentum was transferred to Dimorphos from the escaping impact ejecta than was incident with DART. Therefore, the DART kinetic impact was highly effective in deflecting the asteroid Dimorphos. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: accepted by Nature

Successful Kinetic Impact into an Asteroid for Planetary Defense

Authors: R. Terik Daly, Carolyn M. Ernst, Olivier S. Barnouin, Nancy L. Chabot, Andrew S. Rivkin, Andrew F. Cheng, Elena Y. Adams, Harrison F. Agrusa, Elisabeth D. Abel, Amy L. Alford, Erik I. Asphaug, Justin A. Atchison, Andrew R. Badger, Paul Baki, Ronald-L. Ballouz, Dmitriy L. Bekker, Julie Bellerose, Shyam Bhaskaran, Bonnie J. Buratti, Saverio Cambioni, Michelle H. Chen, Steven R. Chesley, George Chiu, Gareth S. Collins, Matthew W. Cox , et al. (76 additional authors not shown)

Abstract: While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority sp… ▽ More While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority space mission related to asteroid mitigation. NASA's Double Asteroid Redirection Test (DART) mission is the first full-scale test of kinetic impact technology. The mission's target asteroid was Dimorphos, the secondary member of the S-type binary near-Earth asteroid (65803) Didymos. This binary asteroid system was chosen to enable ground-based telescopes to quantify the asteroid deflection caused by DART's impact. While past missions have utilized impactors to investigate the properties of small bodies those earlier missions were not intended to deflect their targets and did not achieve measurable deflections. Here we report the DART spacecraft's autonomous kinetic impact into Dimorphos and reconstruct the impact event, including the timeline leading to impact, the location and nature of the DART impact site, and the size and shape of Dimorphos. The successful impact of the DART spacecraft with Dimorphos and the resulting change in Dimorphos's orbit demonstrates that kinetic impactor technology is a viable technique to potentially defend Earth if necessary. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: Accepted by Nature

Ejecta from the DART-produced active asteroid Dimorphos

Authors: Jian-Yang Li, Masatoshi Hirabayashi, Tony L. Farnham, Jessica M. Sunshine, Matthew M. Knight, Gonzalo Tancredi, Fernando Moreno, Brian Murphy, Cyrielle Opitom, Steve Chesley, Daniel J. Scheeres, Cristina A. Thomas, Eugene G. Fahnestock, Andrew F. Cheng, Linda Dressel, Carolyn M. Ernst, Fabio Ferrari, Alan Fitzsimmons, Simone Ieva, Stavro L. Ivanovski, Teddy Kareta, Ludmilla Kolokolova, Tim Lister, Sabina D. Raducan, Andrew S. Rivkin , et al. (39 additional authors not shown)

Abstract: Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstra… ▽ More Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstrated the activation process of an asteroid from an impact under precisely known impact conditions. Here we report the observations of the DART impact ejecta with the Hubble Space Telescope (HST) from impact time T+15 minutes to T+18.5 days at spatial resolutions of ~2.1 km per pixel. Our observations reveal a complex evolution of ejecta, which is first dominated by the gravitational interaction between the Didymos binary system and the ejected dust and later by solar radiation pressure. The lowest-speed ejecta dispersed via a sustained tail that displayed a consistent morphology with previously observed asteroid tails thought to be produced by impact. The ejecta evolution following DART's controlled impact experiment thus provides a framework for understanding the fundamental mechanisms acting on asteroids disrupted by natural impact. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: accepted by Nature

Direct Imaging Explorations for Companions around Mid-Late M Stars from the Subaru/IRD Strategic Program

Authors: Taichi Uyama, Charles Beichman, Masayuki Kuzuhara, Markus Janson, Takayuki Kotani, Dimitri Mawet, Bun'ei Sato, Motohide Tamura, Hiroyuki Tako Ishikawa, Bryson Cale, Thayne Currie, Hiroki Harakawa, Thomas Henning, Teruyuki Hirano, Klaus Hodapp, Yasunori Hori, Masato Ishizuka, Shane Jacobson, Yui Kasagi, Eiichiro Kokubo, Mihoko Konishi, Tomoyuki Kudo, Takashi Kurokawa, Nobuhiko Kusakabe, Jungmi Kwon , et al. (10 additional authors not shown)

Abstract: The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this… ▽ More The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this RV information remains useful for studying the stellar multiplicity of nearby M~dwarfs. To search for companions around such ``deprioritized" M~dwarfs, we observed 14 IRD-SSP targets using Keck/NIRC2 observations with pyramid wavefront sensing at NIR wavelengths, leading to high sensitivity to substellar-mass companions within a few arcseconds. We detected two new companions (LSPM~J1002+1459~B and LSPM~J2204+1505~B) and two new candidates that are likely companions (LSPM~J0825+6902~B and LSPM~J1645+0444~B) as well as one known companion. Including two known companions resolved by the IRD fiber injection module camera, we detected seven (four new) companions at projected separations between $\sim2-20$~au in total. A comparison of the colors with the spectral library suggests that LSPM~J2204+1505~B and LSPM~J0825+6902~B are located at the boundary between late-M and early-L spectral types. Our deep high-contrast imaging for targets where no bright companions were resolved did not reveal any additional companion candidates. The NIRC2 detection limits could constrain potential substellar-mass companions ($\sim10-75\ M_{\rm Jup}$) at 10~au or further. The failure with Keck/NIRC2 around the IRD-SSP stars having significant RV trends makes these objects promising targets for further RV monitoring or deeper imaging with JWST to search for smaller-mass companions below the NIRC2 detection limits. △ Less

Submitted 26 February, 2023; originally announced February 2023.

Comments: 16 pages, 8 figures, accepted for publication in AJ

Comparisons of the core and mantle compositions of earth analogs from different terrestrial planet formation scenarios

Authors: Jesse T. Gu, Rebecca A. Fischer, Matthew C. Brennan, Matthew S. Clement, Seth A. Jacobson, Nathan A. Kaib, David P. O'Brien, Sean N. Raymond

Abstract: The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, G… ▽ More The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores and the initial embryo masses in N-body simulations. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Absence of extended atmospheres in low-mass star radius-gap planets GJ 9827 b, GJ 9827 d and TOI-1235 b

Authors: Vigneshwaran Krishnamurthy, Teruyuki Hirano, Eric Gaidos, Bunei Sato, Ravi Kopparapu, Thomas Barclay, Katherine Garcia-Sage, Hiroki Harakawa, Klaus Hodapp, Shane Jacobson, Mihoko Konishi, Takayuki Kotani, Tomoyuki Kudo, Takashi Kurokawa, Masayuki Kuzuhara, Eric Lopez, Jun Nishikawa, Masashi Omiya, Joshua E. Schlieder, Takuma Serizawa, Motohide Tamura, Akitoshi Ueda, Sebastien Vievard

Abstract: \textit{Kepler} showed a paucity of planets with radii of 1.5 - 2 $\mathrm R_{\oplus}$ around solar mass stars but this radius-gap has not been well studied for low-mass star planets. Energy-driven escape models like photoevaporation and core-powered mass-loss predict opposing transition regimes between rocky and non-rocky planets when compared to models depicting planets forming in gas-poor envir… ▽ More \textit{Kepler} showed a paucity of planets with radii of 1.5 - 2 $\mathrm R_{\oplus}$ around solar mass stars but this radius-gap has not been well studied for low-mass star planets. Energy-driven escape models like photoevaporation and core-powered mass-loss predict opposing transition regimes between rocky and non-rocky planets when compared to models depicting planets forming in gas-poor environments. Here we present transit observations of three super-Earth sized planets in the radius-gap around low-mass stars using high-dispersion InfraRed Doppler (IRD) spectrograph on the Subaru 8.2m telescope. The planets GJ 9827 b and d orbit around a K6V star and TOI-1235 b orbits a M0.5 star. We limit any planet-related absorption in the 1083.3 nm lines of triplet He I by placing an upper-limit on the equivalent width of 14.71 mÅ, 18.39 mÅ, and 1.44 mÅ, for GJ 9827 b (99% confidence), GJ 9827 d (99% confidence) and TOI-1235 b (95% confidence) respectively. Using a Parker wind model, we cap the mass-loss at $>$0.25 $\mathrm M_{\oplus}$ Gyr$^{-1}$ and $>$0.2 $\mathrm M_{\oplus}$ Gyr$^{-1}$ for GJ 9827 b and d, respectively (99% confidence), and $>$0.05 $\mathrm M_{\oplus}$ Gyr$^{-1}$ for TOI-1235 b (95\% confidence) for a representative wind temperature of 5000 K. Our observed results for the three planets are more consistent with the predictions from photoevaporation and/or core-powered mass-loss models than the gas-poor formation models. However, more planets in the radius-gap regime around the low-mass stars are needed to robustly predict the atmospheric evolution in planets around low-mass stars. △ Less

Submitted 17 March, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: Published for MNRAS. 12 pages, 15 figures

Rocky histories: The effect of high excitations on the formation of rocky planets

Authors: Jennifer Scora, Diana Valencia, Alessandro Morbidelli, Seth Jacobson

Abstract: Rocky planets both in and outside of our solar system are observed to have a range of core-mass fractions (CMFs). Imperfect collisions can preferentially strip mantle material from a planet, changing its CMF, and are therefore thought to be the most likely cause of this observed CMF variation. However, previous work that implements these collisions into N-body simulations of planet formation has s… ▽ More Rocky planets both in and outside of our solar system are observed to have a range of core-mass fractions (CMFs). Imperfect collisions can preferentially strip mantle material from a planet, changing its CMF, and are therefore thought to be the most likely cause of this observed CMF variation. However, previous work that implements these collisions into N-body simulations of planet formation has struggled to reliably form high CMF super-Earths. In this work, we specify our initial conditions and simulation parameters to maximize the prevalence of high-energy, CMF-changing collisions in order to form planets with highly diverse CMFs. High-energy collisions have a large $v_{imp}/v_{esc}$ ratio, so we maximize this ratio by starting simulations with high-eccentricity and inclination disks to increase the difference in their orbital velocities, maximizing $v_{imp}$. Additionally, we minimize $v_{esc}$ by starting with small embryos. The final planets undergo more high-energy, debris-producing collisions, and experience significant CMF change over their formation. However, we find that a number of processes work together to average out the CMF of a planet over time, therefore we do not consistently form high-CMF, high mass planets. We do form high-CMF planets below 0.5 $M_{\oplus}$. Additionally, we find in these highly eccentric environments, loss of debris mass due to collisional grinding has a significant effect on final planet masses and CMFs, resulting in smaller planets and a higher average planet CMF. This work highlights the importance of improving measurements of high-density planets to better constrain their CMFs. △ Less

Submitted 2 December, 2022; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 19 pages. This paper has been accepted to ApJ

Planet(esimal)s Around Stars with TESS (PAST) III: A Search for Triplet He I in the Atmospheres of Two 200 Myr-old Planets

Authors: Eric Gaidos, Teruyuki Hirano, Rena A. Lee, Hiroki Harakawa, Klaus Hodapp, Shane Jacobson, Takayuki Kotani, Tomoyuki Kudo, Takashi Kurokawa, Masayuki Kuzuhara, Jun Nishikawa, Masashi Omiya, Takuma Serizawa, Motohide Tamura, Akitoshi Ueda, Sebastien Vievard

Abstract: We report a search for excess absorption in the 1083.2 nm line of ortho (triplet) helium during transits of TOI-1807b and TOI-2076b, 1.25 and 2.5R$_{\rm Earth}$ planets on 0.55- and 10.4-day orbits around nearby $\sim$200~Myr-old K dwarf stars. We limit the equivalent width of any transit-associated absorption to $<$4 and $<$8 mA, respectively. We limit the escape of solar-composition atmospheres… ▽ More We report a search for excess absorption in the 1083.2 nm line of ortho (triplet) helium during transits of TOI-1807b and TOI-2076b, 1.25 and 2.5R$_{\rm Earth}$ planets on 0.55- and 10.4-day orbits around nearby $\sim$200~Myr-old K dwarf stars. We limit the equivalent width of any transit-associated absorption to $<$4 and $<$8 mA, respectively. We limit the escape of solar-composition atmospheres from TOI-1807b and TOI-2076b to $\lesssim$1 and $\lesssim$0.1M$_{\rm Earth}$ Gyr$^{-1}$, respectively, depending on wind temperature. The absence of a H/He signature for TOI-1807b is consistent with a measurement of mass indicating a rocky body and the prediction by a hydrodynamic model that any H-dominated atmosphere would be unstable and already have been lost. Differential spectra obtained during the transit of TOI-2076b contain a He I-like feature, but this closely resembles the stellar line and extends beyond the transit interval. Until additional transits are observed, we suspect this to be the result of variation in the stellar He I line produced by rotation of active regions and/or flaring on the young, active host star. Non-detection of escape could mean that TOI-2076b is more massive than expected, the star is less EUV-luminous, the models overestimate escape, or the planet has a H/He-poor atmosphere that is primarily molecules such as H$_2$O. Photochemical models of planetary winds predict a semi-major axis at which triplet He I observations are most sensitive to mass loss: TOI-2076b orbits near this optimum. Future surveys could use a distance criterion to increase the yield of detections. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Comments: Accepted to MNRAS

Metal-silicate partitioning of W and Mo and the role of carbon in controlling their abundances in the Bulk Silicate Earth

Authors: E. S. Jennings, S. A. Jacobson, D. C. Rubie, Y. Nakajima, A. K. Vogel, L. A. Rose-Weston, D. J. Frost

Abstract: The liquid metal-liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf-W isotopic system is used to date early planetary evolution. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous exper… ▽ More The liquid metal-liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf-W isotopic system is used to date early planetary evolution. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous experiments to re-examine the systematics of Mo and W partitioning. W partitioning is particularly sensitive to silicate and metallic melt compositions and becomes more siderophile with increasing temperature. We show that W has a 6+ oxidation state in silicate melts over the full experimental fO2 range of $Δ$IW -1.5 to -3.5. Mo has a 4+ oxidation state and its partitioning is less sensitive to silicate melt composition, but also depends on metallic melt composition. DMo stays approximately constant with increasing depth in Earth. Both W and Mo become more siderophile with increasing C content of the metal, so we fit epsilon interaction parameters. W and Mo along with C are incorporated into a combined N-body accretion and core-mantle differentiation model. We show that W and Mo require the early accreting Earth to be sulfur-depleted and carbon-enriched so that W and Mo are efficiently partitioned into Earth's core and do not accumulate in the mantle. If this is the case, the produced Earth-like planets possess mantle compositions matching the BSE for all simulated elements. However, there are two distinct groups of estimates of the bulk mantle's C abundance in the literature: low (100 ppm), and high (800 ppm), and all models are consistent with the higher estimated carbon abundance. The low BSE C abundance would be achievable when the effects of the segregation of dispersed metal droplets produced in deep magma oceans by the disproportionation of Fe2+ to Fe3+ plus metallic Fe is considered. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Comments: 61 pages including supplemental material

Journal ref: Geochimica et Cosmochimica Acta. 293, 40-69 (2021)

Noble gases and stable isotopes track the origin and early evolution of the Venus atmosphere

Authors: Guillaume Avice, Rita Parai, Seth Jacobson, Jabrane Labidi, Melissa G. Trainer, Mihail P. Petkov

Abstract: The composition the atmosphere of Venus results from the integration of many processes entering into play over the entire geological history of the planet. Determining the elemental abundances and isotopic ratios of noble gases (He, Ne, Ar, Kr, Xe) and stable isotopes (H, C, N, O, S) in the Venus atmosphere is a high priority scientific target since it could open a window on the origin and early e… ▽ More The composition the atmosphere of Venus results from the integration of many processes entering into play over the entire geological history of the planet. Determining the elemental abundances and isotopic ratios of noble gases (He, Ne, Ar, Kr, Xe) and stable isotopes (H, C, N, O, S) in the Venus atmosphere is a high priority scientific target since it could open a window on the origin and early evolution of the entire planet. This chapter provides an overview of the existing dataset on noble gases and stable isotopes in the Venus atmosphere. The current state of knowledge on the origin and early and long-term evolution of the Venus atmosphere deduced from this dataset is summarized. A list of persistent and new unsolved scientific questions stemming from recent studies of planetary atmospheres (Venus, Earth and Mars) are described. Important mission requirements pertaining to the measurement of volatile elements in the atmosphere of Venus as well as potential technical difficulties are outlined. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 40 pages, 10 figures, 1 table

First tests of a 1 megapixel near-infrared avalanche photodiode array for ultra-low background space astronomy

Authors: Charles-Antoine Claveau, Michael Bottom, Shane Jacobson, Klaus Hodapp, Aidan Walk, Markus Loose, Ian Baker, Egle Zemaityte, Matthew Hicks, Keith Barnes, Richard Powell, Ryan Bradley, Eric Moore

Abstract: Spectroscopy of Earth-like exoplanets and ultra-faint galaxies are priority science cases for the coming decades. Here, broadband source flux rates are measured in photons per square meter per hour, imposing extreme demands on detector performance, including dark currents lower than 1 e-/pixel/kilosecond, read noise less than 1 e-/pixel/frame, and large formats. There are currently no infrared det… ▽ More Spectroscopy of Earth-like exoplanets and ultra-faint galaxies are priority science cases for the coming decades. Here, broadband source flux rates are measured in photons per square meter per hour, imposing extreme demands on detector performance, including dark currents lower than 1 e-/pixel/kilosecond, read noise less than 1 e-/pixel/frame, and large formats. There are currently no infrared detectors that meet these requirements. The University of Hawaii and industrial partners are developing one promising technology, linear mode avalanche photodiodes (LmAPDs), using fine control over the HgCdTe bandgap structure to enable noise-free charge amplification and minimal glow. Here we report first results of a prototype megapixel format LmAPD operated in our cryogenic testbed. At 50 Kelvin, we measure a dark current of about 3 e-/pixel/kilosecond, which is due to an intrinsic dark current consistent with zero (best estimate of 0.1 e-/pixel/kilosecond) and a ROIC glow of 0.08 e-/pixel/frame. The read noise of these devices is about 10 e-/pixel/frame at 3 volts, and decreases by 30% with each additional volt of bias, reaching 2 e- at 8 volts. Upcoming science-grade devices are expected to substantially improve upon these figures, and address other issues uncovered during testing. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: To be published in Proc. SPIE 12191, X-Ray, Optical, and Infrared Detectors for Astronomy X, 12191-34 (July 18, 2022)

Impact Induced Atmosphere-Mantle Exchange Sets the Volatile Elemental Ratios on Primitive Earths

Authors: Howard Chen, Seth A. Jacobson

Abstract: Conventional planet formation theory suggests that chondritic materials have delivered crucial atmospheric and hydrospheric elements such as carbon (C), nitrogen (N), and hydrogen (H) onto primitive Earth. However, recent measurements highlight the significant elemental ratio discrepancies between terrestrial parent bodies and the supposed planet building blocks. Here we present a volatile evoluti… ▽ More Conventional planet formation theory suggests that chondritic materials have delivered crucial atmospheric and hydrospheric elements such as carbon (C), nitrogen (N), and hydrogen (H) onto primitive Earth. However, recent measurements highlight the significant elemental ratio discrepancies between terrestrial parent bodies and the supposed planet building blocks. Here we present a volatile evolution model during the assembly of Earth and Earth-like planets. Our model includes impact losses, atmosphere-mantle exchange, and time dependent effects of accretion and outgassing calculated from dynamical modeling outcomes. Exploring a wide range of planetesimal properties (i.e., size and composition) as well as impact history informed by N-body accretion simulations, we find that while the degree of CNH fractionation has inherent stochasticity, the evolution of C/N and C/H ratios can be traced back to certain properties of the protoplanet and projectiles. Interestingly, the majority of our Earth-like planets acquire superchondritic final C/N ratios, implying that the volatile elemental ratios on terrestrial planets are driven by the complex interplay between delivery, atmospheric ablation, and mantle degassing. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 24 pages, 10 figures, 3 tables; accepted to Earth and Planetary Science Letters

Stochastic accretion of the Earth

Authors: Paolo A. Sossi, Ingo L. Stotz, Seth A. Jacobson, Alessandro Morbidelli, Hugh St. C. O'Neill

Abstract: Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to… ▽ More Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to that expected from loss by partial vaporisation/condensation. Here we reconcile these observations by showing, using N-body simulations, that Earth accreted stochastically from many precursor bodies whose variable compositions reflect the temperatures at which they formed. Impact-induced atmospheric loss was efficient only when the proto-Earth was small, and elements that accreted thereafter retain near-chondritic isotope ratios. Earth's composition is reproduced when initial temperatures of planetesimal- to embryo-sized bodies are set by disk accretion rates of (1.08 $\pm$ 0.17) $\times$ 10$^{-7}$ solar masses/yr, although they may be perturbed by $^{26}$Al heating on bodies formed at different times. The model implies a heliocentric gradient in composition and rapid planetesimal formation within $\sim$ 1 Myr, in accord with radiometric volatile depletion ages of Earth. △ Less

Submitted 17 July, 2022; originally announced July 2022.

Comments: 13 pages, 4 figures. Nat Astron (2022)

Predictions for the Dynamical States of the Didymos System before and after the Planned DART Impact

Authors: Derek C. Richardson, Harrison F. Agrusa, Brent Barbee, William F. Bottke, Andrew F. Cheng, Siegfried Eggl, Fabio Ferrari, Masatoshi Hirabayashi, Özgür Karatekin, Jay McMahon, Stephen R. Schwartz, Ronald-Louis Ballouz, Adriano Campo Bagatin, Elisabetta Dotto, Eugene G. Fahnestock, Oscar Fuentes-Muñoz, Ioannis Gkolias, Douglas P. Hamilton, Seth A. Jacobson, Martin Jutzi, Josh Lyzhoft, Rahil Makadia, Alex J. Meyer, Patrick Michel, Ryota Nakano , et al. (11 additional authors not shown)

Abstract: NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a reali… ▽ More NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $β$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, 13 figures, published in PSJ

Journal ref: Planet. Sci. J. 3 157 (2022)

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

Authors: Hiroki Harakawa, Takuya Takarada, Yui Kasagi, Teruyuki Hirano, Takayuki Kotani, Masayuki Kuzuhara, Masashi Omiya, Hajime Kawahara, Akihiko Fukui, Yasunori Hori, Hiroyuki Tako Ishikawa, Masahiro Ogihara, John Livingston, Timothy D. Brandt, Thayne Currie, Wako Aoki, Charles A. Beichman, Thomas Henning, Klaus Hodapp, Masato Ishizuka, Hideyuki Izumiura, Shane Jacobson, Markus Janson, Eiji Kambe, Takanori Kodama , et al. (24 additional authors not shown)

Abstract: We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass… ▽ More We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: Accepted for publication in PASJ (May 23, 2022)

Early Solar System instability triggered by dispersal of the gaseous disk

Authors: Beibei Liu, Sean N. Raymond, Seth A. Jacobson

Abstract: The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets. However, the instability trigger and timing have not been clearly established. Hydrodynamical modeling has shown that while the Sun's gaseous protoplanetary disk was present the giant planets migrated into a compact orbital configuration in a chain of resonances. He… ▽ More The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets. However, the instability trigger and timing have not been clearly established. Hydrodynamical modeling has shown that while the Sun's gaseous protoplanetary disk was present the giant planets migrated into a compact orbital configuration in a chain of resonances. Here we use dynamical simulations to show that the giant planets' instability was likely triggered by the dispersal of the gaseous disk. As the disk evaporated from the inside-out, its inner edge swept successively across and dynamically perturbed each planet's orbit in turn. The associated orbital shift caused a dynamical compression of the exterior part of the system, ultimately triggering instability. The final orbits of our simulated systems match those of the Solar System for a viable range of astrophysical parameters. The giant planet instability therefore took place as the gaseous disk dissipated, constrained by astronomical observations to be a few to ten million years after the birth of the Solar System. Terrestrial planet formation would not complete until after such an early giant planet instability; the growing terrestrial planets may even have been sculpted by its perturbations, explaining the small mass of Mars relative to Earth. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: Authors' version, 12 pages, 4 figures

Validation and atmospheric exploration of the sub-Neptune TOI-2136b around a nearby M3 dwarf

Authors: K. Kawauchi, F. Murgas, E. Palle, N. Narita, A. Fukui, T. Hirano, H. Parviainen, H. T. Ishikawa, N. Watanabe, E. Esparaza-Borges, M. Kuzuhara, J. Orell-Miquel, V. Krishnamurthy, M. Mori, T. Kagetani, Y. Zou, K. Isogai, J. H. Livingston, S. B. Howell, N. Crouzet, J. P. de Leon, T. Kimura, T. Kodama, J. Korth, S. Kurita , et al. (29 additional authors not shown)

Abstract: The NASA space telescope $TESS$ is currently in the extended mission of its all-sky search for new transiting planets. Of the thousands of candidates that TESS is expected to deliver, transiting planets orbiting nearby M dwarfs are particularly interesting targets since they provide a great opportunity to characterize their atmospheres by transmission spectroscopy. We aim to validate and character… ▽ More The NASA space telescope $TESS$ is currently in the extended mission of its all-sky search for new transiting planets. Of the thousands of candidates that TESS is expected to deliver, transiting planets orbiting nearby M dwarfs are particularly interesting targets since they provide a great opportunity to characterize their atmospheres by transmission spectroscopy. We aim to validate and characterize the new sub-Neptune-sized planet candidate TOI-2136.01 orbiting a nearby M dwarf ($d = 33.36 \pm 0.02$ pc, $T_{eff} = 3373 \pm 108$ K) with an orbital period of 7.852 days. We use TESS data, ground-based multicolor photometry, and radial velocity measurements with the InfraRed Doppler (IRD) instrument on the Subaru Telescope to validate the planetary nature of TOI-2136.01 and estimate the stellar and planetary parameters. We also conduct high-resolution transmission spectroscopy to search for helium in its atmosphere. We confirmed that TOI-2136.01 (now named as TOI-2136b) is a bona fide planet with a planetary radius of $R_p = 2.2 \pm 0.07$ $R_{Earth}$ and a mass of $M_p = 4.7^{+3.1}_{-2.6}$ $M_{Earth}$. We also search for helium 10830 Å absorption lines and place an upper limit on the equivalent width of $<$ 7.8 mÅ and on the absorption signal of $<$ 1.44 % with 95 % confidence. TOI-2136b is a sub-Neptune transiting a nearby and bright star (J=10.8) and is a potentially hycean planet, which is a new class of habitable planets with large oceans under a H$_2$-rich atmosphere, making it an excellent target for atmospheric studies to understand the formation, evolution, and habitability of the small planets. △ Less

Submitted 11 July, 2022; v1 submitted 21 February, 2022; originally announced February 2022.

Comments: 19 pages, 16 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 666, A4 (2022)

The Metal-Silicate Partitioning of Carbon During Earth's Accretion and its Distribution in the Early Solar System

Authors: I. Blanchard, D. C. Rubie, E. S. Jennings, I. A. Franchi, X. Zhao, S. Petitgirard, N. Miyajima, S. A. Jacobson, A. Morbidelli

Abstract: Carbon is an essential element for the existence and evolution of life on Earth. Its abundance in Earth's crust and mantle (the Bulk Silicate Earth, BSE) is surprisingly high given that carbon is strongly siderophile (metal-loving) at low pressures and temperatures, and hence should have segregated almost completely into Earth's core during accretion. Estimates of the concentration of carbon in th… ▽ More Carbon is an essential element for the existence and evolution of life on Earth. Its abundance in Earth's crust and mantle (the Bulk Silicate Earth, BSE) is surprisingly high given that carbon is strongly siderophile (metal-loving) at low pressures and temperatures, and hence should have segregated almost completely into Earth's core during accretion. Estimates of the concentration of carbon in the BSE lie in the range 100-260 ppm and are much higher than expected based on simple models of core-mantle differentiation. Here we show through experiments at the putative conditions of Earth's core formation (49-71 GPa and 3600-4000 K) that carbon is significantly less siderophile at these conditions than at the low pressures ($\leq$13 GPa) and temperatures ($\leq$2500 K) of previous large volume press studies, but at least an order of magnitude more siderophile than proposed recently based on an experimental approach that is similar to ours. Using our new data along with previously published results, we derive a new parameterization of the pressure-temperature dependence of the metal-silicate partitioning of carbon. We apply this parameterization in a model that combines planet formation and core-mantle differentiation that is based on astrophysical N-body accretion simulations. Because differentiated planetesimals were almost completely depleted in carbon due to sublimation at high temperatures, almost all carbon in the BSE was added by the accretion of fully-oxidized carbonaceous chondrite material from the outer solar system. Carbon is added to the mantle continuously throughout accretion and its concentration reaches values within the BSE range (e.g. 140$\pm$40 ppm) at the end of accretion. The corresponding final core and bulk Earth carbon concentrations are 1270$\pm$300 ppm and 495$\pm$125 ppm respectively. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Supplementary material included in pdf

Journal ref: Earth and Planetary Science Letters. 580, 117374 (2022)

Silicon and Strontium abundances of very metal-poor stars determined from near-infrared spectra

Authors: Wako Aoki, Timothy C. Beers, Satoshi Honda, Hiroyuki T. Ishikawa, Tadafumi Matsuno, Vinicius M. Placco, Jinmi Yoon, Hiroki Harakawa, Teruyuki Hirano, Klaus Hodapp, Masato Ishizuka, Shane Jacobson, Takayuki Kotani, Tomoyuki Kudo, Takashi Kurokawa, Masayuki Kuzuhara, Jun Nishikawa, Masashi Omiya, Takuma Serizawa, Motohide Tamura, Akitoshi Ueda, Sebastien Vievard

Abstract: Silicon and Strontium are key elements to explore the nucleosynthesis and chemical evolution of the Galaxy by measurements of very metal-poor stars. There are, however, only a few useful spectral lines of these elements in the optical range that are measurable for such low-metallicity stars. Here we report on abundances of these two elements determined from near-infrared high-resolution spectra ob… ▽ More Silicon and Strontium are key elements to explore the nucleosynthesis and chemical evolution of the Galaxy by measurements of very metal-poor stars. There are, however, only a few useful spectral lines of these elements in the optical range that are measurable for such low-metallicity stars. Here we report on abundances of these two elements determined from near-infrared high-resolution spectra obtained with the Subaru Telescope Infrared Doppler instrument (IRD). Si abundances are determined for as many as 26 Si lines for six very and extremely metal-poor stars (-4.0<[Fe/H]<-1.5), which significantly improves the reliability of the abundance measurements. All six stars, including three carbon-enhanced objects, show over-abundances of Si ([Si/Fe]~+0.5). Two stars with [Fe/H]~-1.5 have relatively small over-abundances. The [Mg/Si] ratios agree with the solar value, except for one metal-poor star with carbon excess. Strontium abundances are determined from the triplet lines for four stars, including two for the first time. The consistency of the Sr abundances determined from near-infrared and optical spectra require further examination from additional observations. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 3 figures, 6 tables, PASJ, in press

The "Breaking The Chains" migration model for super-Earths formation: the effect of collisional fragmentation

Authors: Leandro Esteves, André Izidoro, Bertram Bitsch, Seth A. Jacobson, Sean N. Raymond, Rogerio Deienno, Othon C. Winter

Abstract: Planets between 1-4 Earth radii with orbital periods <100 days are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disk inner edge in chains of mean motion resonances. After gas disk dispersal, simulations show that super-Earth's gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts.… ▽ More Planets between 1-4 Earth radii with orbital periods <100 days are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disk inner edge in chains of mean motion resonances. After gas disk dispersal, simulations show that super-Earth's gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts. The instability phase is key to matching the orbital spacing of observed systems. Yet, most previous simulations have modelled collisions as perfect accretion events, ignoring fragmentation. In this work, we investigate the impact of imperfect accretion on the breaking the chains scenario. We performed N-body simulations starting from distributions of planetary embryos and modelling the effects of pebble accretion and migration in the gas disk. Our simulations also follow the long-term dynamical evolution of super-Earths after the gas disk dissipation. We compared the results of simulations where collisions are treated as perfect merging events with those where imperfect accretion and fragmentation are allowed. We concluded that the perfect accretion is a suitable approximation in this regime, from a dynamical point of view. Although fragmentation events are common, only ~10% of the system mass is fragmented during a typical "late instability phase", with fragments being mostly reacreted by surviving planets. This limited total mass in fragments proved to be insufficient to alter qualitatively the final system dynamical configuration -- e.g. promote strong dynamical friction or residual migration -- compared to simulations where fragmentation is neglected. △ Less

Submitted 29 October, 2021; originally announced November 2021.

Comments: 13 pages, 8 figures and 1 table. Accepted for publication in MNRAS

MSC Class: 85-10

Zodiacal Exoplanets in Time (ZEIT) XIII: Planet Orbits and Atmospheres in the V1298 Tau System, a Keystone in Studies of Early Planetary Evolution

Authors: E. Gaidos, T. Hirano, C. Beichman, J. Livingston, H. Harakawa, K. W. Hodapp, M. Ishizuka, S. Jacobson, M. Konishi, T. Kotani, T. Kudo, T. Kurokawa, M. Kuzuhara, J. Nishikawa, M. Omiya, T. Serizawa, M. Tamura, A. Ueda, S. Vievard

Abstract: Studies of planetary systems of stars in star-forming regions and young clusters open a window on the formative stages of planetary evolution. We obtained high-cadence high-resolution infrared spectroscopy of the solar-mass Taurus association-member V1298 Tau during a transit of its 10R$_{\oplus}$-size "b" planet. We measured the systemic radial velocity and find that the kinematics of V1298 Tau s… ▽ More Studies of planetary systems of stars in star-forming regions and young clusters open a window on the formative stages of planetary evolution. We obtained high-cadence high-resolution infrared spectroscopy of the solar-mass Taurus association-member V1298 Tau during a transit of its 10R$_{\oplus}$-size "b" planet. We measured the systemic radial velocity and find that the kinematics of V1298 Tau suggest an affiliation with a $\gtrsim$6 Myr-old subgroup. A comparison of V1298 Tau and the nearby, co-moving star 2M0405 with stellar evolution models suggests an age of $\sim$10-25 Myr. We measured the projected spin-orbit angle of "b" as $λ=15_{-16}^{+15}$ and $λ= 2_{-4}^{+12}$ degrees using the apparent RV shift and change in line profile, respectively, induced by the transient occultation of the rotating star by the planet. These values indicate a prograde orbit like that of the interior "c" planet of V1298 Tau and point to a co-planar multi-planet system that formed within a disk. We also measured variation in the strength of the 1083 nm triplet of neutral orthohelium as a probe of any extended/escaping atmosphere around "b". We detect a steady decrease in absorption over the transit which appears to arise from the star or its planetary system. While this variation could be ascribed to "b" or possibly to the immediately preceding transit of "d", we cannot rule out that this is due to rapid variation in the stellar disk-integrated flux in the triplet. The amplitude of variation ($\sim$0.04 nm) is consistent with moderate estimates of atmospheric escape driven by XUV radiation from the central star. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS 20 October 2021

Barrel Instability in Binary Asteroids

Authors: Matija Ćuk, Seth A. Jacobson, Kevin J. Walsh

Abstract: Most close-in planetary satellites are in synchronous rotation, which is usually the stable end-point of tidal despinning. Saturn's moon Hyperion is a notable exception by having a chaotic rotation. Hyperion's dynamical state is a consequence of its high eccentricity and its highly prolate shape (Wisdom et al. 1984). As many binary asteroids also have elongated secondaries, chaotic rotation is exp… ▽ More Most close-in planetary satellites are in synchronous rotation, which is usually the stable end-point of tidal despinning. Saturn's moon Hyperion is a notable exception by having a chaotic rotation. Hyperion's dynamical state is a consequence of its high eccentricity and its highly prolate shape (Wisdom et al. 1984). As many binary asteroids also have elongated secondaries, chaotic rotation is expected for moons in eccentric binaries (Ćuk and Nesvorný, 2010), and a minority of asteroidal secondaries may be in that state (Pravec et al. 2016). The question of secondary rotation is also important for the action of the BYORP effect, which can quickly evolve orbits of synchronous (but not non-synchronous) secondaries (Ćuk and Burns, 2005). Here we report results of a large set of short numerical simulations which indicate that, apart from synchronous and classic chaotic rotation, close-in irregularly-shaped asteroidal secondaries can occupy an additional, intermediate rotational state. In this ``barrel instability'' the secondary slowly rolls along its long axis, while the longest axis is staying largely aligned with the primary-secondary line. This behavior may be more difficult to detect through lightcurves than a fully chaotic rotation, but would likewise shut down BYORP. We show that the binary's eccentricity, separation measured in secondary's radii and the secondary's shape are all important for determining whether the system settles in synchronous rotation, chaotic tumbling, or barrel instability. We compare our results for synthetic asteroids with known binary pairs to determine which of these behaviors may be present in the Near-Earth Asteroid binary population. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: Accepted for PSJ

Zodiacal Exoplanets in Time (ZEIT) XII: A Directly-Imaged Planetary-Mass Companion to a Young Taurus M Dwarf Star

Authors: E. Gaidos, T. Hirano, A. L. Kraus, M. Kuzuhara, Z. Zhang, R. A. Lee, M. Salama, T. A. Berger, S. K. Grunblatt, M. Ansdell, M. C. Liu, H. Harakawa, K. W. Hodapp, S. Jacobson, M. Konishi, T. Kotani, T. Kudo, T. Kurokawa, J. Nishikawa, M. Omiya, T. Serizawa, M. Tamura, A. Ueda, S. Vievard

Abstract: We report the discovery of a resolved (0".9) substellar companion to a member of the 1-5 Myr Taurus star-forming region. The host star (2M0437) is a single mid-M type ($T_{eff}\approx$3100K) dwarf with a position, space motion, and color-magnitude that support Taurus membership, and possible affiliation with a $\sim$2.5 Myr-old sub-group. A comparison with stellar models suggests a 2-5 Myr age and… ▽ More We report the discovery of a resolved (0".9) substellar companion to a member of the 1-5 Myr Taurus star-forming region. The host star (2M0437) is a single mid-M type ($T_{eff}\approx$3100K) dwarf with a position, space motion, and color-magnitude that support Taurus membership, and possible affiliation with a $\sim$2.5 Myr-old sub-group. A comparison with stellar models suggests a 2-5 Myr age and a mass of 0.15-0.18M$_{\odot}$. Although K2 detected quasi-periodic dimming from close-in circumstellar dust, the star lacks detectable excess infrared emission from a circumstellar disk and its H$α$ emission is not commensurate with accretion. Astrometry based on three years of AO imaging shows that the companion (2M0437b) is co-moving, while photometry of two other sources at larger separation indicates they are likely heavily-reddened background stars. A comparison of the luminosity of 2M0437b with models suggests a mass of 3-5$M_{\rm jup}$, well below the deuterium burning limit, and an effective temperature of 1400-1500K, characteristic of a late L spectral type. The $H$-$K$ color is redder than the typical L dwarf, but comparable to other directly detected young planets, e.g. those around HR 8799. The discovery of a super-Jupiter around a very young, very low mass star challenges models of planet formation by either core accretion (which requires time) or disk instability (which requires mass). We also detected a second, co-moving, widely-separated (75") object which appears to be a heavily-extincted star. This is certainly a fellow member of this Taurus sub-group and statistically likely to be a bound companion. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS 14 October 2021

The Excited Spin State of Dimorphos Resulting from the DART Impact

Authors: Harrison F. Agrusa, Ioannis Gkolias, Kleomenis Tsiganis, Derek C. Richardson, Alex J. Meyer, Daniel J. Scheeres, Matija Ćuk, Seth A. Jacobson, Patrick Michel, Özgür Karatekin, Andrew F. Cheng, Masatoshi Hirabayashi, Yun Zhang, Eugene G. Fahnestock, Alex B. Davis

Abstract: The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and t… ▽ More The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and two numerical) to investigate the post-impact attitude stability of Dimorphos as a function of its axial ratios, $a/b$ and $b/c$ ($a \ge b \ge c$), and the momentum transfer efficiency $β$. The first method uses a novel analytic approach in which we assume a circular orbit and a point-mass primary that identifies four fundamental frequencies of motion corresponding to the secondary's mean motion, libration, precession, and nutation frequencies. At resonance locations among these four frequencies, we find that attitude instabilities are possible. Using two independent numerical codes, we recover many of the resonances predicted by the analytic model and indeed show attitude instability. With one code, we use fast Lyapunov indicators to show that the secondary's attitude can evolve chaotically near the resonance locations. Then, using a high-fidelity numerical model, we find that Dimorphos enters a chaotic tumbling state near the resonance locations and is especially prone to unstable rotation about its long axis, which can be confirmed by ESA's Hera mission arriving at Didymos in late 2026. We also show that a fully coupled treatment of the spin and orbital evolution of both bodies is crucial to accurately model the long-term evolution of the secondary's spin state and libration amplitude. Finally, we discuss the implications of a post-impact tumbling or rolling state, including the possibility of terminating BYORP evolution if Dimorphos is no longer in synchronous rotation. △ Less

Submitted 29 July, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 38 pages, 13 figures, Accepted for publication in Icarus

Non-detection of Helium in the upper atmospheres of TRAPPIST-1b, e and f

Authors: Vigneshwaran Krishnamurthy, Teruyuki Hirano, Guðmundur Stefánsson, Joe P. Ninan, Suvrath Mahadevan, Eric Gaidos, Ravi Kopparapu, Bunei Sato, Yasunori Hori, Chad F. Bender, Caleb I. Cañas, Scott A. Diddams, Samuel Halverson, Hiroki Harakawa, Suzanne Hawley, Fred Hearty, Leslie Hebb, Klaus Hodapp, Shane Jacobson, Shubham Kanodia, Mihoko Konishi, Takayuki Kotani, Adam Kowalski, Tomoyuki Kudo, Takashi Kurokawa , et al. (15 additional authors not shown)

Abstract: We obtained high-resolution spectra of the ultra-cool M-dwarf TRAPPIST-1 during the transit of its planet `b' using two high dispersion near-infrared spectrographs, IRD instrument on the Subaru 8.2m telescope and HPF instrument on the 10m Hobby-Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet `b' using the APO/ARCTIC, which assiste… ▽ More We obtained high-resolution spectra of the ultra-cool M-dwarf TRAPPIST-1 during the transit of its planet `b' using two high dispersion near-infrared spectrographs, IRD instrument on the Subaru 8.2m telescope and HPF instrument on the 10m Hobby-Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet `b' using the APO/ARCTIC, which assisted us to capture the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets `b', `e' and `f' from IRD, we attempt to constrain the atmospheric escape of the planet using the He I triplet 10830 Å absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet related absorption, we place an upper limit on the equivalent widths of <7.754 mÅ for planet `b', <10.458 mÅ for planet `e', and <4.143 mÅ for planet `f' at 95% confidence from the IRD data, and <3.467 mÅ for planet `b' at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This non-detection of extended atmospheres having low mean-molecular weight in all three planets aids in further constraining their atmospheric composition by steering the focus towards the search of high molecular weight species in their atmospheres. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: 11 pages; 4 figures; Accepted for publication in AJ

The Effect of Inefficient Accretion on Planetary Differentiation

Authors: Saverio Cambioni, Seth A. Jacobson, Alexandre Emsenhuber, Erik Asphaug, David C. Rubie, Travis S. J. Gabriel, Stephen R. Schwartz, Roberto Furfaro

Abstract: Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core… ▽ More Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 21 pages, 2 tables, 7 figures. Published open access on PSJ: https://iopscience.iop.org/article/10.3847/PSJ/abf0ad

Journal ref: Planet. Sci. J. 2 93 (2021)

Large Adaptive Optics Survey for Substellar Objects (LASSO) Around Young, Nearby, Low-mass Stars with Robo-AO

Authors: Maissa Salama, James Ou, Christoph Baranec, Michael C. Liu, Brendan P. Bowler, Paul Barnes, Morgan Bonnet, Mark Chun, Dmitry A. Duev, Sean Goebel, Don Hall, Shane Jacobson, Rebecca Jensen-Clem, Nicholas M. Law, Charles Lockhart, Reed Riddle, Heather Situ, Eric Warmbier, Zhoujian Zhang

Abstract: We present results from the Large Adaptive optics Survey for Substellar Objects (LASSO), where the goal is to directly image new substellar companions (<70 M$_{Jup}$) at wide orbital separations ($\gtrsim$50 AU) around young ($\lesssim$300 Myrs), nearby (<100 pc), low-mass ($\approx$0.1-0.8 M$_{\odot}$) stars. We report on 427 young stars imaged in the visible (i') and near-infrared (J or H) simul… ▽ More We present results from the Large Adaptive optics Survey for Substellar Objects (LASSO), where the goal is to directly image new substellar companions (<70 M$_{Jup}$) at wide orbital separations ($\gtrsim$50 AU) around young ($\lesssim$300 Myrs), nearby (<100 pc), low-mass ($\approx$0.1-0.8 M$_{\odot}$) stars. We report on 427 young stars imaged in the visible (i') and near-infrared (J or H) simultaneously with Robo-AO on the Kitt Peak 2.1-m telescope and later the Maunakea University of Hawaii 2.2-m telescope. To undertake the observations, we commissioned a new infrared camera for Robo-AO that uses a low-noise high-speed SAPHIRA avalanche photodiode detector. We detected 121 companion candidates around 111 stars, of which 62 companions are physically associated based on Gaia DR2 parallaxes and proper motions, another 45 require follow-up observations to confirm physical association, and 14 are background objects. The companion separations range from 2-1101 AU and reach contrast ratios of 7.7 magnitudes in the near infrared compared to the primary. The majority of confirmed and pending candidates are stellar companions, with ~5 being potentially substellar and requiring follow-up observations for confirmation. We also detected a 43$\pm$9 M$_{Jup}$ and an 81$\pm$5 M$_{Jup}$ companion that were previously reported. We found 34 of our targets have acceleration measurements detected using Hipparcos-Gaia proper motions. Of those, 58$^{+12}_{-14}$% of the 12 stars with imaged companion candidates have significant accelerations ($χ^2 >11.8$), while only 23$^{+11}_{-6}$% of the remaining 22 stars with no detected companion have significant accelerations. The significance of the acceleration decreases with increasing companion separation. These young accelerating low-mass stars with companions will eventually yield dynamical masses with future orbit monitoring. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: 44 pages: 20 pages of text with 14 figures, 24 pages of tables in the appendix. Accepted for publication in the Astronomical Journal

How Sublimation Delays the Onset of Dusty Debris Disk Formation Around White Dwarf Stars

Authors: Jordan K. Steckloff, John Debes, Amy Steele, Brandon Johnson, Elisabeth R. Adams, Seth A. Jacobson, Alessondra Springmann

Abstract: Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host m… ▽ More Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host more dynamically unstable planetary systems. Here we apply a sophisticated material sublimation model to white dwarf systems to show that these statistics are actually a natural result of the interplay of thermal and tidal forces, and show how they define the circumstellar regions where dusty debris disks can form. We demonstrate that these processes tend to prevent stability against both sublimative destruction and reaccretion into planetesimals for rocky materials until white dwarfs cool to below ~25,000-32,000 K, in agreement with the observed limit of ~27,000 K. For pure water ice, this critical temperature is less than 2,700 K (requiring a cooling age older the universe); this precludes pure water ice-rich debris disks forming through the accepted two-step mechanism. The critical temperature is size-dependent; more massive white dwarfs could potentially host dusty debris disks at warmer temperatures. Our model suggests that the location of the disks within the PG 0010+280, GD 56, GD 362, and PG 1541+651 systems are consistent with a forsterite-dominated olivine composition. We also find that very cool white dwarfs may simultaneously host multiple, independently formed dusty debris disks, consistent with observations of the LSPM J0207+3331 system. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: Astrophysical Journal Letters. Accepted for Publication

Two Bright M Dwarfs Hosting Ultra-Short-Period Super-Earths with Earth-like Compositions

Authors: Teruyuki Hirano, John H. Livingston, Akihiko Fukui, Norio Narita, Hiroki Harakawa, Hiroyuki Tako Ishikawa, Kohei Miyakawa, Tadahiro Kimura, Akifumi Nakayama, Naho Fujita, Yasunori Hori, Keivan G. Stassun, Allyson Bieryla, Charles Cadieux, David R. Ciardi, Karen A. Collins, Masahiro Ikoma, Andrew Vanderburg, Thomas Barclay, C. E. Brasseur, Jerome P. de Leon, John P. Doty, René Doyon, Emma Esparza-Borges, Gilbert A. Esquerdo , et al. (36 additional authors not shown)

Abstract: We present observations of two bright M dwarfs (TOI-1634 and TOI-1685: $J=9.5-9.6$) hosting ultra-short period (USP) planets, identified by the TESS mission. The two stars are similar in temperature, mass, and radius ($T_\mathrm{eff}\,\approx\,3500$ K, $M_\star\,\approx\,0.45-0.46\,M_\odot$, and $R_\star\approx 0.45-0.46\,R_\odot$), and the planets are both super-Earth-sized (… ▽ More We present observations of two bright M dwarfs (TOI-1634 and TOI-1685: $J=9.5-9.6$) hosting ultra-short period (USP) planets, identified by the TESS mission. The two stars are similar in temperature, mass, and radius ($T_\mathrm{eff}\,\approx\,3500$ K, $M_\star\,\approx\,0.45-0.46\,M_\odot$, and $R_\star\approx 0.45-0.46\,R_\odot$), and the planets are both super-Earth-sized ($1.25\,R_\oplus<R_p<2.0\,R_\oplus$). For both systems, light curves from the ground-based photometry exhibit planetary transits, whose depths are consistent with those by the TESS photometry. We also refine the transit ephemerides based on the ground-based photometry, finding the orbital periods of $P=0.9893436\pm0.0000020$ day and $P=0.6691416\pm0.0000019$ day for TOI-1634b and TOI-1685b, respectively. Through intensive radial velocity (RV) observations using IRD on the Subaru 8.2m telescope, we confirm the planetary nature of the TOIs, and measure their masses: $10.14\pm0.95\,M_\oplus$ and $3.43\pm0.93\,M_\oplus$ for TOI-1634b and TOI-1685b, respectively, when the observed RVs are fitted with a single-planet circular-orbit model. Combining those with the planet radii of $R_p=1.749\pm 0.079\,R_\oplus$ (TOI-1634b) and $1.459\pm0.065\,R_\oplus$ (TOI-1685b), we find that both USP planets have mean densities consistent with an Earth-like internal composition, which is typical for small USP planets. TOI-1634b is currently the most massive USP planet in this category, and it resides near the radius valley, which makes it a benchmark planet in the context of discussing the size limit of rocky planet cores as well as testing the formation scenarios for USP planets. Excess scatter in the RV residuals for TOI-1685 suggests the presence of a possible secondary planet or unknown activity/instrumental noise in the RV data, but further observations are required to check those possibilities. △ Less

Submitted 5 July, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 29 pages, 21 figures, 4 tables; AJ in press

First Detection of Hydroxyl Radical Emission from an Exoplanet Atmosphere: High-dispersion Characterization of WASP-33b using Subaru/IRD

Authors: Stevanus K. Nugroho, Hajime Kawahara, Neale P. Gibson, Ernst J. W. de Mooij, Teruyuki Hirano, Takayuki Kotani, Yui Kawashima, Kento Masuda, Matteo Brogi, Jayne L. Birkby, Chris A. Watson, Motohide Tamura, Konstanze Zwintz, Hiroki Harakawa, Tomoyuki Kudo, Masayuki Kuzuhara, Klaus Hodapp, Masato Ishizuka, Shane Jacobson, Mihoko Konishi, Takashi Kurokawa, Jun Nishikawa, Masashi Omiya, Takuma Serizawa, Akitoshi Ueda , et al. (1 additional authors not shown)

Abstract: We report the first detection of a hydroxyl radical (OH) emission signature in the planetary atmosphere outside the solar system, in this case, in the day-side of WASP-33b. We analyze high-resolution near-infrared emission spectra of WASP-33b taken using the InfraRed Doppler spectrograph on the 8.2-m Subaru telescope. The telluric and stellar lines are removed using a de-trending algorithm, SysRem… ▽ More We report the first detection of a hydroxyl radical (OH) emission signature in the planetary atmosphere outside the solar system, in this case, in the day-side of WASP-33b. We analyze high-resolution near-infrared emission spectra of WASP-33b taken using the InfraRed Doppler spectrograph on the 8.2-m Subaru telescope. The telluric and stellar lines are removed using a de-trending algorithm, SysRem. The residuals are then cross-correlated with OH and H$_{2}$O planetary spectrum templates produced using several different line-lists. We check and confirm the accuracy of OH line-lists by cross-correlating with the spectrum of GJ 436. As a result, we detect the emission signature of OH at $K_\mathrm{p}$ of 230.9$^{+6.9}_{-7.4}$ km s$^{-1}$ and $v_{\mathrm{sys}}$ of $-$0.3$^{+5.3}_{-5.6}$ km s$^{-1}$ with S/N of 5.4 and significance of 5.5$σ$. Additionally, we marginally detect H$_{2}$O emission in the H-band with S/N of 4.0 and significance of 5.2$σ$ using the POKAZATEL line-list. However, no significant signal is detected using the HITEMP 2010, which might be due to differences in line positions and strengths, as well as the incompleteness of the line-lists. Nonetheless, this marginal detection is consistent with the prediction that H$_{2}$O is mostly thermally dissociated in the upper atmosphere of the ultra-hot Jupiters. Therefore, along with CO, OH is expected to be one of the most abundant O-bearing molecules in the day-side atmosphere of ultra-hot Jupiters and should be considered when studying their atmosphere. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: Accepted for publication on ApJL; 12 pages, 6 figures

Quantitative estimates of impact induced crustal erosion during accretion and its influence on the Sm/Nd ratio of the Earth

Authors: Laetitia Allibert, Sébastien Charnoz, Julien Siebert, Seth A. Jacobson, Sean N. Raymond

Abstract: Dynamical scenarios of terrestrial planets formation involve strong perturbations of the inner part of the solar system by the giant-planets, leading to enhanced impact velocities and subsequent collisional erosion. We quantitatively estimate the effect of collisional erosion on the resulting composition of Earth, and estimate how it may provide information on the dynamical context of its formatio… ▽ More Dynamical scenarios of terrestrial planets formation involve strong perturbations of the inner part of the solar system by the giant-planets, leading to enhanced impact velocities and subsequent collisional erosion. We quantitatively estimate the effect of collisional erosion on the resulting composition of Earth, and estimate how it may provide information on the dynamical context of its formation. We simulate and quantify the erosion of Earth's crust in the context of Solar System formation scenarios, including the classical model and Grand Tack scenario that invokes orbital migration of Jupiter during the gaseous disk phase (Walsh et al., 2011; Raymond et al., 2018). We find that collisional erosion of the early crust is unlikely to produce an excess of about 6% of the Sm/Nd ratio in terrestrial rock samples compared to chondrites for most simulations. Only Grand Tack simulations in which the last giant impact on Earth occurred later than 50 million years after the start of Solar System formation can account for such an offset. However, this time frame is consistent with current cosmochemical and dynamical estimates of the Moon forming impact (Chyba, 1991; Walker, 2009; Touboul et al.,2007, 2009, 2015; Pepin and Porcelli, 2006; Norman et al., 2003; Nyquist et al., 2006; Boyet et al.,2015). Such a late fractionation in the Sm/Nd ratio is unlikely to be responsible for a 20-ppm $^{142}$Nd excess in terrestrial rocks due to the half life of the radiogenic system. Additionally, such a large and late fractionation in the Sm/Nd ratio would accordingly induce non-observed anomalies in the $^{143}$Nd/$^{144}$Nd ratio. Considering our results, the Grand Tack model with a late Moon-forming impact cannot be easily reconciled with the Nd isotopic Earth contents. △ Less

Submitted 9 March, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 22 pages, 3 figures, supplementary, to be published

Zodiacal Exoplanets in Time. XI. The Orbit and Radiation Environment of the Young M Dwarf-Hosted Planet K2-25b

Authors: E. Gaidos, T. Hirano, D. J. Wilson, K. France, K. Rockcliffe, E. Newton, G. Feiden, V. Krishnamurthy, H. Harakawa, K. W. Hodapp, M. Ishizuka, S. Jacobson, M. Konishi, T. Kotani, T. Kudo, T. Kurokawa, M. Kuzuhara, J. Nishikawa, M. Omiya, T. Serizawa, M. Tamura, A. Ueda, S. Vievard

Abstract: M dwarf stars are high-priority targets for searches for Earth-size and potentially Earth-like planets, but their planetary systems may form and evolve in very different circumstellar environments than those of solar-type stars. To explore the evolution of these systems, we obtained transit spectroscopy and photometry of the Neptune-size planet orbiting the ~650 Myr-old Hyades M dwarf K2-25. An an… ▽ More M dwarf stars are high-priority targets for searches for Earth-size and potentially Earth-like planets, but their planetary systems may form and evolve in very different circumstellar environments than those of solar-type stars. To explore the evolution of these systems, we obtained transit spectroscopy and photometry of the Neptune-size planet orbiting the ~650 Myr-old Hyades M dwarf K2-25. An analysis of the variation in spectral line shape induced by the Doppler "shadow" of the planet indicate that the planet's orbit is closely aligned with the stellar equator (lambda = -1.7+5.8/-3.7 deg), and that an eccentric orbit found by previous work could arise from perturbations by another planet on a co-planar orbit. We detect no significant variation in the depth of the He I line at 1083 nm during transit. A model of atmospheric escape as a isothermal Parker wind with a solar composition show that this non-detection is not constraining compared to escape rate predictions of ~0.1 Mearth/Gyr; at such rates, at least several Gyr are required for a Neptune-like planet to evolve into a rocky super-Earth. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Comments: to appear in MNRAS Letters

Precision radial velocity measurements by the forward-modeling technique in the near-infrared

Authors: Teruyuki Hirano, Masayuki Kuzuhara, Takayuki Kotani, Masashi Omiya, Tomoyuki Kudo, Hiroki Harakawa, Sébastien Vievard, Takashi Kurokawa, Jun Nishikawa, Motohide Tamura, Klaus Hodapp, Masato Ishizuka, Shane Jacobson, Mihoko Konishi, Takuma Serizawa, Akitoshi Ueda, Eric Gaidos, Bun'ei Sato

Abstract: Precision radial velocity (RV) measurements in the near-infrared are a powerful tool to detect and characterize exoplanets around low-mass stars or young stars with higher magnetic activity. However, the presence of strong telluric absorption lines and emission lines in the near infrared that significantly vary in time can prevent extraction of RV information from these spectra by classical techni… ▽ More Precision radial velocity (RV) measurements in the near-infrared are a powerful tool to detect and characterize exoplanets around low-mass stars or young stars with higher magnetic activity. However, the presence of strong telluric absorption lines and emission lines in the near infrared that significantly vary in time can prevent extraction of RV information from these spectra by classical techniques, which ignore or mask the telluric lines. We present a methodology and pipeline to derive precision RVs from near-infrared spectra using a forward-modeling technique. We applied this to spectra with a wide wavelength coverage (Y, J, and H bands, simultaneously), taken by the InfraRed Doppler (IRD) spectrograph on the Subaru 8.2-m telescope. Our pipeline extracts the instantaneous instrumental profile of the spectrograph for each spectral segment, based on a reference spectrum of the laser-frequency comb that is injected into the spectrograph simultaneously with the stellar light. These profiles are used to derive the intrinsic stellar template spectrum, which is free from instrumental broadening and telluric features, as well as model and fit individual observed spectra in the RV analysis. Implementing a series of numerical simulations using theoretical spectra that mimic IRD data, we test the pipeline and show that IRD can achieve $<2$ m s$^{-1}$ precision for slowly rotating mid-to-late M dwarfs with a signal-to-noise ratio $> 100$ per pixel at 1000 nm. Dependences of RV precision on various stellar parameters (e.g., $T_{\rm eff}$, $v\sin i$, [Fe/H]) and the impact of telluric-line blendings on the RV accuracy are discussed through the mock spectra analyses. We also apply the RV-analysis pipeline to the observed spectra of GJ 699 and TRAPPIST-1, demonstrating that the spectrograph and the pipeline are capable of an RV accuracy of $<3$ m s$^{-1}$ at least on a time scale of a few months. △ Less

Submitted 28 August, 2020; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: 21 pages, 20 figures; accepted for publication in PASJ; typos corrected

Small Bodies Tell the Story of the Solar System: A Scientific Rationale for a Multi-Target Small Body Sample Return Program including the Earth-based Laboratory Analysis of Returned Samples

Authors: Seth A. Jacobson, Maitrayee Bose, Dennis Bodewits, Marc Fries, Devanshu Jha, Prajkta Mane, Larry Nittler, Scott Sandford, Michelle Thompson

Abstract: Small bodies are time-capsules of different eras of solar system history from the most primitive materials within the solar system to evolved pieces of larger bodies. A small body sample return program is an essential component of small body exploration, and such a program should include opportunities for both missions and laboratory analysis. Small bodies are time-capsules of different eras of solar system history from the most primitive materials within the solar system to evolved pieces of larger bodies. A small body sample return program is an essential component of small body exploration, and such a program should include opportunities for both missions and laboratory analysis. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Comments: A White Paper submitted to the Planetary Science and Astrobiology Decadal Survey 2023-2032

Limits on the Spin-Orbit Angle and Atmospheric Escape for the 22 Myr-old Planet AU Mic b

Authors: Teruyuki Hirano, Vigneshwaran Krishnamurthy, Eric Gaidos, Heather Flewelling, Andrew W. Mann, Norio Narita, Peter Plavchan, Takayuki Kotani, Motohide Tamura, Hiroki Harakawa, Klaus Hodapp, Masato Ishizuka, Shane Jacobson, Mihoko Konishi, Tomoyuki Kudo, Takashi Kurokawa, Masayuki Kuzuhara, Jun Nishikawa, Masashi Omiya, Takuma Serizawa, Akitoshi Ueda, Sébastien Vievard

Abstract: We obtained spectra of the pre-main sequence star AU Microscopii during a transit of its Neptune-sized planet to investigate its orbit and atmosphere. We used the high-dispersion near-infrared spectrograph IRD on the Subaru telescope to detect the Doppler "shadow" from the planet and constrain the projected stellar obliquity. Modeling of the observed planetary Doppler shadow suggests a spin-orbit… ▽ More We obtained spectra of the pre-main sequence star AU Microscopii during a transit of its Neptune-sized planet to investigate its orbit and atmosphere. We used the high-dispersion near-infrared spectrograph IRD on the Subaru telescope to detect the Doppler "shadow" from the planet and constrain the projected stellar obliquity. Modeling of the observed planetary Doppler shadow suggests a spin-orbit alignment of the system ($λ=-4.7_{-6.4}^{+6.8}$ degrees), but additional observations are needed to confirm this finding. We use both the IRD data and spectra obtained with NIRSPEC on Keck-II to search for absorption in the 1083 nm line of metastable triplet He I by the planet's atmosphere and place an upper limit for the equivalent width of 3.7 mÅat 99 $\%$ confidence. With this limit and a Parker wind model we constrain the escape rate from the atmosphere to $<0.15-0.45\, M_{\oplus}$ Gyr$^{-1}$, comparable to the rates predicted by an XUV energy-limited escape calculation and hydrodynamic models, but refinement of the planet mass is needed for rigorous tests. △ Less

Submitted 20 July, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 10 pages, 5 figures, accepted for publication in ApJ Letters

Scaling laws for the geometry of an impact-induced magma ocean

Authors: Miki Nakajima, Gregor J. Golabek, Kai Wünnemann, David C. Rubie, Christoph Burger, Henry J. Melosh, Seth A. Jacobson, Lukas Manske, Scott D. Hull

Abstract: Here, we develop scaling laws for (1) the distribution of impact-induced heat within the mantle and (2) shape of the impact-induced melt based on more than 100 smoothed particle hydrodynamic (SPH) simulations. We use Legendre polynomials to describe these scaling laws and determine their coefficients by linear regression, minimizing the error between our model and SPH simulations. The input parame… ▽ More Here, we develop scaling laws for (1) the distribution of impact-induced heat within the mantle and (2) shape of the impact-induced melt based on more than 100 smoothed particle hydrodynamic (SPH) simulations. We use Legendre polynomials to describe these scaling laws and determine their coefficients by linear regression, minimizing the error between our model and SPH simulations. The input parameters are the impact angle $θ$ ($0^{\circ}, 30^{\circ}, 60^{\circ}$, and $90^{\circ}$), total mass $M_T$ ($1M_{\rm Mars}-53M_{\rm Mars}$, where $M_{\rm Mars}$ is the mass of Mars), impact velocity $v_{\rm imp}$ ($v_{\rm esc} - 2v_{\rm esc}$, where $v_{\rm esc}$ is the mutual escape velocity), and impactor-to-total mass ratio $γ$ ($0.03-0.5$). We find that the equilibrium pressure at the base of a melt pool can be higher (up to $\approx 80 \%$) than those of radially-uniform global magma ocean models. This could have a significant impact on element partitioning. These melt scaling laws are publicly available on GitHub ($\href{https://github.com/mikinakajima/MeltScalingLaw}{https://github.com/mikinakajima/MeltScalingLaw}$). △ Less

Submitted 2 June, 2021; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: Published in Earth and Planetary Science Letters

Zodiacal Exoplanets in Time. X. The Orbit and Atmosphere of the Young "Neptune Desert"-Dwelling Planet K2-100b

Authors: E. Gaidos, T. Hirano, A. W. Mann, D. A. Owens, T. A. Berger, K. France, A. Vanderburg, H. Harakawa, K. W. Hodapp, M. Ishizuka, S. Jacobson, M. Konishi, T. Kotani, T. Kudo, T. Kurokawa, J. Nishikawa, M. Omiya, T. Serizawa, M. Tamura, A. Oeda

Abstract: We obtained high-resolution infrared spectroscopy and short-cadence photometry of the 600-800 Myr Praesepe star K2-100 during transits of its 1.67-day planet. This Neptune-size object, discovered by the NASA K2 mission, is an interloper in the "desert" of planets with similar radii on short period orbits. Our observations can be used to understand its origin and evolution by constraining the orbit… ▽ More We obtained high-resolution infrared spectroscopy and short-cadence photometry of the 600-800 Myr Praesepe star K2-100 during transits of its 1.67-day planet. This Neptune-size object, discovered by the NASA K2 mission, is an interloper in the "desert" of planets with similar radii on short period orbits. Our observations can be used to understand its origin and evolution by constraining the orbital eccentricity by transit fitting, measuring the spin-orbit obliquity by the Rossiter-McLaughlin effect, and detecting any extended, escaping hydrogen-helium envelope with the 10830A line of neutral helium in the 2s3S triplet state. Transit photometry with 1-min cadence was obtained by the K2 satellite during Campaign 18 and transit spectra were obtained with the IRD spectrograph on the Subaru telescope. While the elevated activity of K2-100 prevented us from detecting the Rossiter-McLaughlin effect, the new photometry combined with revised stellar parameters allowed us to constrain the eccentricity to e < 0.15/0.28 with 90%/99% confidence. We modeled atmospheric escape as an isothermal, spherically symmetric Parker wind, with photochemistry driven by UV radiation that we estimate by combining the observed spectrum of the active Sun with calibrations from observations of K2-100 and similar young stars in the nearby Hyades cluster. Our non-detection (<5.7mA) of a transit-associated He I line limits mass loss of a solar-composition atmosphere through a T<10000K wind to <0.3Me/Gyr. Either K2-100b is an exceptional desert-dwelling planet, or its mass loss is occurring at a lower rate over a longer interval, consistent with a core accretion-powered scenario for escape. △ Less

Submitted 28 March, 2020; originally announced March 2020.

Comments: Accepted to MNRAS

Subsolar Al/Si and Mg/Si ratios of non-carbonaceous chondrites reveal planetesimal formation during early condensation in the protoplanetary disk

Authors: A. Morbidelli, G. Libourel, H. Palme, S. A. Jacobson, D. C. Rubie

Abstract: The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the s… ▽ More The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the subsequent formation of residual condensates with strongly sub-solar Al/Si and Mg/Si ratios. The mixing of residual condensates with different amounts of material with solar refractory/Si element ratios explains the Al/Si and Mg/Si values of non-carbonaceous chondrites. To match quantitatively the observed ratios, we find that the first-planetesimals should have accreted when the disk temperature was ~1,330-1,400 K depending on pressure and assuming a solar C/O ratio of the disk. We discuss how this model relates to our current understanding of disk evolution, grain dynamics, and planetesimal formation. We also extend the discussion to moderately volatile elements (e.g., Na), explaining how it may be possible that the depletion of these elements in non-carbonaceous chondrites is correlated with the depletion of refractory elements (e.g., Al). Extending the analysis to Cr, we find evidence for a higher than solar C/O ratio in the protosolar disk's gas when/where condensation from a fractionated gas occurred. Finally, we discuss the possibility that the supra-solar Al/Si and Mg/Si ratios of the Earth are due to the accretion of ~40% of the mass of our planet from the first-generation of refractory-rich planetesimals. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: Accepted for publication in EPSL

Chemical Diversity of Super-Earths As a Consequence of Formation

Authors: Jennifer Scora, Diana Valencia, Alessandro Morbidelli, Seth A. Jacobson

Abstract: Recent observations of rocky super-Earths have revealed an apparent wider distribution of Fe/Mg ratios, or core to mantle ratios, than the planets in our Solar System. This study aims to understand how much of the chemical diversity in the super-Earth population can arise from giant impacts during planetary formation. Planet formation simulations have only recently begun to treat collisions more r… ▽ More Recent observations of rocky super-Earths have revealed an apparent wider distribution of Fe/Mg ratios, or core to mantle ratios, than the planets in our Solar System. This study aims to understand how much of the chemical diversity in the super-Earth population can arise from giant impacts during planetary formation. Planet formation simulations have only recently begun to treat collisions more realistically in an attempt to replicate the planets in our Solar System. We investigate planet formation more generally by simulating the formation of rocky super-Earths with varying initial conditions using a version of SyMBA, a gravitational N-body code, that incorporates realistic collisions. We track the maximum plausible change in composition after each impact. The final planets span a range of Fe/Mg ratios similar to the Solar System planets, but do not completely match the distribution in super-Earth data. We only form a few planets with minor iron-depletion, suggesting other mechanisms are at work. The most iron-rich planets have a lower Fe/Mg ratio than Mercury, and are less enriched than planets such as Kepler-100b. This indicates that further work on our understanding of planet formation and further improvement of precision of mass and radius measurements are required to explain planets at the extremes of this Fe/Mg distribution. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 15 pages, 14 figures. Accepted for publication in MNRAS