Showing 1–50 of 250 results for author: Blake, G

The high-contrast performance of the Keck Planet Imager and Characterizer

Authors: Jason J. Wang, Dimitri Mawet, Jerry W. Xuan, Chih-Chun Hsu, Jean-Baptiste Ruffio, Katelyn Horstman, Yinzi Xin, Jacques-Robert Delorme, Nemanja Jovanovic, Yapeng Zhang, Luke Finnerty, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Gregory W. Doppmann, Daniel Echeverri, Michael P. Fitzgerald, Joshua Liberman, Ronald Lopez, Evan Morris, Jacklyn Pezzato-Rovner, Ben Sappey, Tobias Schofield , et al. (3 additional authors not shown)

Abstract: The Keck Planet Imager and Characterizer (KPIC), a series of upgrades to the Keck II Adaptive Optics System and Instrument Suite, aims to demonstrate high-resolution spectroscopy of faint exoplanets that are spatially resolved from their host stars. In this paper, we measure KPIC's sensitivity to companions as a function of separation (i.e., the contrast curve) using on-sky data collected over fou… ▽ More The Keck Planet Imager and Characterizer (KPIC), a series of upgrades to the Keck II Adaptive Optics System and Instrument Suite, aims to demonstrate high-resolution spectroscopy of faint exoplanets that are spatially resolved from their host stars. In this paper, we measure KPIC's sensitivity to companions as a function of separation (i.e., the contrast curve) using on-sky data collected over four years of operation. We show that KPIC is able to reach contrasts of $1.3 \times 10^{-4}$ at 90 mas and $9.2 \times 10^{-6}$ at 420 mas separation from the star, and that KPIC can reach planet-level sensitivities at angular separations within the inner working angle of coronagraphic instruments such as GPI and SPHERE. KPIC is also able to achieve more extreme contrasts than other medium-/high-resolution spectrographs that are not as optimized for high-contrast performance. We decompose the KPIC performance budget into individual noise terms and discuss limiting factors. The fringing that results from combining a high-contrast imaging system with a high-resolution spectrograph is identified as an important source of systematic noise. After mitigation and correction, KPIC is able to reach within a factor of 2 of the photon noise limit at separations < 200 mas. At large separations, KPIC is limited by the background noise performance of NIRSPEC. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 16 pages, 6 figures, submitted to the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, 13096-69

Are these planets or brown dwarfs? Broadly solar compositions from high-resolution atmospheric retrievals of ~10-30 $M_\textrm{Jup}$ companions

Authors: Jerry W. Xuan, Chih-Chun Hsu, Luke Finnerty, Jason J. Wang, Jean-Baptiste Ruffio, Yapeng Zhang, Heather A. Knutson, Dimitri Mawet, Eric E. Mamajek, Julie Inglis, Nicole L. Wallack, Marta L. Bryan, Geoffrey A. Blake, Paul Mollière, Neda Hejazi, Ashley Baker, Randall Bartos, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Michael P. Fitzgerald, Nemanja Jovanovic, Joshua Liberman , et al. (10 additional authors not shown)

Abstract: Using Keck Planet Imager and Characterizer (KPIC) high-resolution ($R$~35000) spectroscopy from 2.29-2.49 $μ$m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ~10-30 $M_\textrm{Jup}$, orbital separations spanning ~50-360 au, and $T_\textrm{eff}$ between ~1500-2600 K. We find that all companions have solar C/O ratios, and metallicities, to within the… ▽ More Using Keck Planet Imager and Characterizer (KPIC) high-resolution ($R$~35000) spectroscopy from 2.29-2.49 $μ$m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ~10-30 $M_\textrm{Jup}$, orbital separations spanning ~50-360 au, and $T_\textrm{eff}$ between ~1500-2600 K. We find that all companions have solar C/O ratios, and metallicities, to within the 1-2$σ$ level, with the measurements clustered around solar composition. Stars in the same stellar associations as our systems have near-solar abundances, so these results indicate that this population of companions is consistent with formation via direct gravitational collapse. Alternatively, core accretion outside the CO snowline would be compatible with our measurements, though the high mass ratios of most systems would require rapid core assembly and gas accretion in massive disks. On a population level, our findings can be contrasted with abundance measurements for directly imaged planets with m<10 $M_\textrm{Jup}$, which show tentative atmospheric metal enrichment. In addition, the atmospheric compositions of our sample of companions are distinct from those of hot Jupiters, which most likely form via core accretion. For two companions with $T_\textrm{eff}$~1700-2000 K (kap And b and GSC 6214-210 b), our best-fit models prefer a non-gray cloud model with >3$σ$ significance. The cloudy models yield 2-3$σ$ lower $T_\textrm{eff}$ for these companions, though the C/O and [C/H] still agree between cloudy and clear models at the $1σ$ level. Finally, we constrain 12CO/13CO for three companions with the highest S/N data (GQ Lup b, HIP 79098 b, and DH Tau b), and report $v$sin($i$) and radial velocities for all companions. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted to ApJ, 52 pages, 20 figures

Rotation and Abundances of the Benchmark Brown Dwarf HD 33632 Ab from Keck/KPIC High-resolution Spectroscopy

Authors: Chih-Chun Hsu, Jason J. Wang, Jerry W. Xuan, Jean-Baptiste Ruffio, Daniel Echeverri, Yinzi Xin, Joshua Liberman, Luke Finnerty, Evan Morris, Katelyn Horstman, Ben Sappey, Gregory W. Doppmann, Dimitri Mawet, Nemanja Jovanovic, Michael P. Fitzgerald, Jacques-Robert Delorme, J. Kent Wallace, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Ronald A. López, Jacklyn Pezzato, Tobias Schofield , et al. (2 additional authors not shown)

Abstract: We present the projected rotational velocity and molecular abundances for HD 33632 Ab obtained via Keck Planet Imager and Characterizer high-resolution spectroscopy. HD 33632 Ab is a nearby benchmark brown dwarf companion at a separation of $\sim$20 au that straddles the L/T transition. Using a forward-modeling framework with on-axis host star spectra, self-consistent substellar atmospheric and re… ▽ More We present the projected rotational velocity and molecular abundances for HD 33632 Ab obtained via Keck Planet Imager and Characterizer high-resolution spectroscopy. HD 33632 Ab is a nearby benchmark brown dwarf companion at a separation of $\sim$20 au that straddles the L/T transition. Using a forward-modeling framework with on-axis host star spectra, self-consistent substellar atmospheric and retrieval models for HD 33632 Ab, we derive a projected rotational velocity of 53 $\pm$ 3 km/s and carbon/water mass fractions of log CO = $-$2.3 $\pm$ 0.3 and log H$_2$O = $-$2.7 $\pm$ 0.2. The inferred carbon-to-oxygen ratio (C/O = 0.58 $\pm$ 0.14), molecular abundances, and metallicity ([C/H] = 0.0 $\pm$ 0.2 dex) of HD 33632 Ab are consistent with its host star. Although detectable methane opacities are expected in L/T transition objects, we did not recover methane in our KPIC spectra, partly due to the high $v\sin{i}$ and to disequilibrium chemistry at the pressures we are sensitive to. We parameterize the spin as the ratio of rotation over break-up velocity, and compare HD 33632 Ab to a compilation of >200 very low-mass objects (M$\lesssim$0.1 M$_{\odot}$) that have spin measurements in the literature. There appears to be no clear trend for the isolated field low-mass objects versus mass, but a tentative trend is identified for low-mass companions and directly imaged exoplanets, similar to previous findings. A larger sample of close-in gas giant exoplanets and brown dwarfs will critically examine our understanding of their formation and evolution through rotation and chemical abundance measurements. △ Less

Submitted 18 June, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: Accepted for publication in the Astrophysical Journal. 36 pages, 15 figures, 5 tables

Fresh view of the hot brown dwarf HD 984 B through high-resolution spectroscopy

Authors: J. C. Costes, J. W. Xuan, A. Vigan, J. Wang, V. D'Orazi, P. Mollière, A. Baker, R. Bartos, G. A. Blake, B. Calvin, S. Cetre, J. Delorme, G. Doppmann, D. Echeveri, L. Finnerty, M. P. Fitzgerald, C. Hsu, N. Jovanovic, R. Lopez, D. Mawet, E. Morris, J. Pezzato, C. L. Phillips, J. Ruffio, B. Sappey , et al. (5 additional authors not shown)

Abstract: Context. High-resolution spectroscopy has the potential to drive a better understanding of the atmospheric composition, physics, and dynamics of young exoplanets and brown dwarfs, bringing clear insights into the formation channel of individual objects. Aims. Using the Keck Planet Imager and Characterizer (KPIC; R = 35,000), we aim to characterize a young brown dwarf HD 984 B. By measuring its C/O… ▽ More Context. High-resolution spectroscopy has the potential to drive a better understanding of the atmospheric composition, physics, and dynamics of young exoplanets and brown dwarfs, bringing clear insights into the formation channel of individual objects. Aims. Using the Keck Planet Imager and Characterizer (KPIC; R = 35,000), we aim to characterize a young brown dwarf HD 984 B. By measuring its C/O and 12CO/13CO ratios, we expect to gain new knowledge about its origin by confirming the difference in the formation pathways between brown dwarfs and super-Jupiters. Methods. We analysed the KPIC high-resolution spectrum (2.29-2.49 μm) of HD 984 B using an atmospheric retrieval framework based on nested sampling and petitRADTRANS, using both clear and cloudy models. Results. Using our best-fit model, we find C/O = 0.50+0.01-0.01 (0.01 is the statistical error) for HD 984 B which agrees with that of its host star within 1σ (0.40+0.20-0.20). We also retrieve an isotopolog 12CO/13CO ratio of 98+20-25 in its atmosphere, which is similar to that of the Sun. In addition, HD 984 B has a substellar metallicity with [Fe/H] = -0.62+0.02-0.02. Finally, we find that most of the retrieved parameters are independent of our choice of retrieval model. Conclusions. From our measured C/O and 12CO/13CO, the favored formation mechanism of HD 984 B seems to be via gravitational collapse or disk instability and not core accretion, which is a favored formation mechanism for giant exoplanets with m < 13 MJup and semimajor axis between 10 and 100 au. However, with only a few brown dwarfs with a measured 12CO/13CO ratio, similar analyses using high-resolution spectroscopy will become essential in order to determine planet formation processes more precisely. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 14 pages, 8 figures

Journal ref: A&A 686, A294 (2024)

Orbital and Atmospheric Characterization of the 1RXS J034231.8+121622 System Using High-Resolution Spectroscopy Confirms That The Companion is a Low-Mass Star

Authors: Clarissa R. Do Ó, Ben Sappey, Quinn M. Konopacky, Jean-Baptiste Ruffio, Kelly K. O'Neil, Tuan Do, Gregory Martinez, Travis S. Barman, Jayke S. Nguyen, Jerry W. Xuan, Christopher A. Theissen, Sarah Blunt, William Thompson, Chih-Chun Hsu, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Julie Inglis , et al. (11 additional authors not shown)

Abstract: The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a… ▽ More The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a temperature of 2510 $\pm$ 50 K and a metallicity of $0.13\substack{+0.12 \\ -0.11}$. Recent work suggests this system is associated with the Hyades, placing it an older age than previous estimates. Both metallicities agree with current $[Fe/H]$ Hyades measurements (0.11 -- 0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, of 0.30 $\pm$ 0.15 $M_\odot$ and 0.08 $\pm$ 0.01 $M_\odot$ (84 $\pm$ 11 $M_{Jup}$) respectively. Using the RVs and a new astrometry point from Keck/NIRC2, we find that the system is likely an edge-on, moderately eccentric ($0.41\substack{+0.27 \\ -0.08}$) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 $\pm$ 0.10 (primary) and 0.55 $\pm$ 0.10 (companion). From these results, we confirm that this system most likely went through a binary star formation process in the Hyades. The significant changes in this system's parameters since its discovery highlight the importance of high resolution spectroscopy for both orbital and atmospheric characterization of directly imaged companions. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 30 pages, 18 figures, accepted for publication in AJ

Molecular Gas Tracers in Young and Old Protoplanetary Disks

Authors: Dana E. Anderson, L. Ilsedore Cleeves, Geoffrey A. Blake, Chunhua Qi, Edwin A. Bergin, John M. Carpenter, Kamber R. Schwarz, Claire Thilenius, Ke Zhang

Abstract: Molecular emission is used to investigate both the physical and chemical properties of protoplanetary disks. Therefore, to accurately derive disk properties, we need a thorough understanding of the behavior of the molecular probes we rely on. Here we investigate how the molecular line emission of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O compare to other measured quantities in a set of 20 protoplane… ▽ More Molecular emission is used to investigate both the physical and chemical properties of protoplanetary disks. Therefore, to accurately derive disk properties, we need a thorough understanding of the behavior of the molecular probes we rely on. Here we investigate how the molecular line emission of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O compare to other measured quantities in a set of 20 protoplanetary disks. Overall, we find positive correlations between multiple line fluxes and the disk dust mass and radius. We also generally find strong positive correlations between the line fluxes of different molecular species. However, some disks do show noticeable differences in the relative fluxes of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O. These differences occur even within a single star-forming region. This results in a potentially large range of different disk masses and chemical compositions for systems of similar age and birth environment. While we make preliminary comparisons of molecular fluxes across different star-forming regions, more complete and uniform samples are needed in the future to search for trends with birth environment or age. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astrophysical Journal

Validation of elemental and isotopic abundances in late-M spectral types with the benchmark HIP 55507 AB system

Authors: Jerry W. Xuan, Jason J. Wang, Luke Finnerty, Katelyn Horstman, Simon Grimm, Anne Peck, Eric L. Nielsen, Heather A. Knutson, Dimitri Mawet, Howard Isaacson, Andrew W. Howard, Michael C. Liu, Sam Walker, Mark Phillips, Geoffrey Blake, Jean-Baptiste Ruffio, Yapeng Zhang, Julie Inglis, Nicole L. Wallack, Aniket Sanghi, Erica Gonzales, Fei Dai, Ashley Baker, Randall Bartos, Charlotte Bond , et al. (26 additional authors not shown)

Abstract: M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similar $T_{\rm eff}\sim2300-2800~K$. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution ($R\sim35,000$) $K$ band spectroscopy. First, by includi… ▽ More M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similar $T_{\rm eff}\sim2300-2800~K$. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution ($R\sim35,000$) $K$ band spectroscopy. First, by including KPIC relative radial velocities between the primary and secondary in the orbit fit, we improve the dynamical mass precision by 60% and find $M_B=88.0_{-3.2}^{+3.4}$ $M_{\rm Jup}$, putting HIP 55507 B above the stellar-substellar boundary. We also find that HIP 55507 B orbits its K6V primary star with $a=38^{+4}_{-3}$ AU and $e=0.40\pm0.04$. From atmospheric retrievals of HIP 55507 B, we measure $\rm [C/H]=0.24\pm0.13$, $\rm [O/H]=0.15\pm0.13$, and $\rm C/O=0.67\pm0.04$. Moreover, we strongly detect $\rm ^{13}CO$ ($7.8σ$ significance) and tentatively detect $\rm H_2^{18}O$ ($3.7σ$ significance) in companion's atmosphere, and measure $\rm ^{12}CO/^{13}CO=98_{-22}^{+28}$ and $\rm H_2^{16}O/H_2^{18}O=240_{-80}^{+145}$ after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure $\rm ^{12}CO/^{13}CO=79_{-16}^{+21}$ and $\rm C^{16}O/C^{18}O=288_{-70}^{+125}$ for the primary star. These results demonstrate that HIP 55507 A and B have consistent $\rm ^{12} C/^{13}C$ and $\rm ^{16}O/^{18}O$ to the $<1σ$ level, as expected for a chemically homogeneous binary system. Given the similar flux ratios and separations between HIP 55507 AB and systems with young, substellar companions, our results open the door to systematically measuring $\rm ^{13}CO$ and $\rm H_2^{18}O$ abundances in the atmospheres of substellar or even planetary-mass companions with similar spectral types. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: Accepted to ApJ, 28 pages, 14 figures

Atmospheric metallicity and C/O of HD 189733 b from high-resolution spectroscopy

Authors: Luke Finnerty, Jerry W. Xuan, Yinzi Xin, Joshua Liberman, Tobias Schofield, Michael P. Fitzgerald, Shubh Agrawal, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppman, Daniel Echeverri, Chih-Chun Hsu, Nemanja Jovanovic, Ronald A. López, Emily C. Martin, Dimitri Mawet, Evan Morris, Jacklyn Pezzato, Jean-Baptiste Ruffio, Ben Sappey, Andrew Skemer , et al. (5 additional authors not shown)

Abstract: We present high-resolution $K$-band emission spectra of the quintessential hot Jupiter HD 189733 b from the Keck Planet Imager and Characterizer (KPIC). Using a Bayesian retrieval framework, we fit the dayside pressure-temperature profile, orbital kinematics, mass-mixing ratios of H$_2$O, CO, CH$_4$, NH$_3$, HCN, and H$_2$S, and the $\rm ^{13}CO/^{12}CO$ ratio. We measure mass fractions of… ▽ More We present high-resolution $K$-band emission spectra of the quintessential hot Jupiter HD 189733 b from the Keck Planet Imager and Characterizer (KPIC). Using a Bayesian retrieval framework, we fit the dayside pressure-temperature profile, orbital kinematics, mass-mixing ratios of H$_2$O, CO, CH$_4$, NH$_3$, HCN, and H$_2$S, and the $\rm ^{13}CO/^{12}CO$ ratio. We measure mass fractions of $\rm \log H_2O = -2.0^{+0.4}_{-0.4}$ and $\rm \log CO = -2.2^{+0.5}_{-0.5}$, and place upper limits on the remaining species. Notably, we find $\rm \log CH_4 < -4.5$ at 99\% confidence, despite its anticipated presence at the equilibrium temperature of HD 189733 b assuming local thermal equilibrium. We make a tentative ($\sim3σ$) detection of $\rm ^{13}CO$, and the retrieved posteriors suggest a $\rm ^{12}C/^{13}C$ ratio similar to or substantially less than the local interstellar value. The possible $\rm ^{13}C$ enrichment would be consistent with accretion of fractionated material in ices or in the protoplanetary disk midplane. The retrieved abundances correspond to a substantially sub-stellar atmospheric $\rm C/O = 0.3\pm0.1$, while the carbon and oxygen abundances are stellar to slightly super-stellar, consistent with core-accretion models which predict an inverse correlation between C/O and metallicity. The specific combination of low C/O and high metallicity suggests significant accretion of solid material may have occurred late in the formation process of HD 189733 b. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: 17 pages, 7 figures, 2 tables, accepted in AJ

High-contrast JWST-MIRI spectroscopy of planet-forming disks for the JDISC Survey

Authors: Klaus M. Pontoppidan, Colette Salyk, Andrea Banzatti, Ke Zhang, Ilaria Pascucci, Karin I. Oberg, Feng Long, Carlos Munoz-Romero, John Carr, Joan Najita, Geoffrey A. Blake, Nicole Arulanantham, Sean Andrews, Nicholas P. Ballering, Edwin Bergin, Jenny Calahan, Douglas Cobb, Maria Jose Colmenares, Annie Dickson-Vandervelde, Anna Dignan, Joel Green, Phoebe Heretz, Greg Herczeg, Anusha Kalyaan, Sebastian Krijt , et al. (4 additional authors not shown)

Abstract: The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4… ▽ More The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4.9-28 micron spectral range. This is achieved using observations of Themis-family asteroids as precise empirical reference sources. High spectral contrast enables precise retrievals of physical parameters, searches for rare molecular species and isotopologues, and constraints on the inventories of carbon- and nitrogen-bearing species. JDISCS also offers significant improvements to the MRS wavelength and resolving power calibration. We describe the JDISCS calibrated data and demonstrate its quality using observations of the disk around the solar-mass young star FZ Tau. The FZ Tau MIRI spectrum is dominated by strong emission from warm water vapor. We show that the water and CO line emission originates from the disk surface and traces a range of gas temperatures of ~500-1500 K. We retrieve parameters for the observed CO and H2O lines, and show that they are consistent with a radial distribution represented by two temperature components. A high water abundance of n(H2O)~10^-4 fills the disk surface at least out to the 350 K isotherm at 1.5 au. We search the FZ Tau environs for extended emission detecting a large (radius of ~300 au) ring of emission from H2 gas surrounding FZ Tau, and discuss its origin. △ Less

Submitted 16 January, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal

Water-Rich Disks around Late M-stars Unveiled: Exploring the Remarkable Case of Sz114

Authors: Chengyan Xie, Ilaria Pascucci, Feng Long, Klaus M. Pontoppidan, Andrea Banzatti, Anusha Kalyaan, Colette Salyk, Yao Liu, Joan R. Najita, Paola Pinilla, Nicole Arulanantham, Gregory J. Herczeg, John Carr, Edwin A. Bergin, Nicholas P. Ballering, Sebastiaan Krijt, Geoffrey A. Blake, Ke Zhang, Karin I. Oberg, Joel D. Green, the JDISC collaboration

Abstract: We present an analysis of the JDISC JWST/MIRI-MRS spectrum of Sz~114, an accreting M5 star surrounded by a large dust disk with a shallow gap at $\sim 39$ au. The spectrum is molecular-rich: we report the detection of water, CO, CO$_2$, HCN, C$_2$H$_2$, and H$_2$. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water… ▽ More We present an analysis of the JDISC JWST/MIRI-MRS spectrum of Sz~114, an accreting M5 star surrounded by a large dust disk with a shallow gap at $\sim 39$ au. The spectrum is molecular-rich: we report the detection of water, CO, CO$_2$, HCN, C$_2$H$_2$, and H$_2$. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water lines with the 17.22 micron emission surpassing that of most mid-to-late M-star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, flux ratios of C$_2$H$_2$/H$_2$O and HCN/H$_2$O in Sz~114 also resemble those of earlier-type disks, with a slightly elevated CO$_2$/H$_2$O ratio. While accretional heating can boost all infrared lines, the unusual properties of Sz~114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and the presence of dust substructures. The latter delays the inward drift of icy pebbles and help preserve a lower C/O ratio over an extended period. In contrast, mid-to-late M-star disks--which are typically faint, small in size, and likely lack significant substructures--may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of mid-to-late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models. △ Less

Submitted 28 November, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 16 pages, 8 figures, accepted by ApJL

Recovering simulated planet and disk signals using SCALES aperture masking

Authors: Mackenzie Lach, Steph Sallum, Ravinder Banyal, Natalie Batalha, Geoff Blake, Tim Brandt, Zackery Briesemeister, Aditi Desai, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Isabel Kain, Charlie Kilpatrick, Katherine de Kleer, Michael Liu, Bruce Macintosh, Raquel Martinez, Dimitri Mawet, Brittany Miles, Caroline Morley, Imke de Pater, Diana Powell, Patrick Sheehan, Andrew Skemer , et al. (7 additional authors not shown)

Abstract: The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric t… ▽ More The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES' low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Journal ref: SPIE Proceedings Volume 12680, Techniques and Instrumentation for Detection of Exoplanets XI; 1268024 (2023)

Simulating medium-spectral-resolution exoplanet characterization with SCALES angular/reference differential imaging

Authors: Aditi Desai, Stephanie E. Sallum, Ravinder Banyal, Natalie Batalha, Natasha Batalha, Geoff Blake, Tim Brandt, Zack Briesemeister, Katherine de Kleer, Imke de Pater, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Isabel Kain, Charlie Kilpatrick, Mackenzie Lach, Mike Liu, Bruce Macintosh, Raquel A. Martinez, Dimitri Mawet, Brittany Miles, Caroline Morley, Diana Powell, Patrick Sheehan , et al. (8 additional authors not shown)

Abstract: SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution… ▽ More SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts. △ Less

Submitted 18 October, 2023; originally announced October 2023.

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES): driving science cases and expected outcomes

Authors: Steph Sallum, Andrew Skemer, Deno Stelter, Ravinder Banyal, Natalie Batalha, Natasha Batalha, Geoff Blake, Tim Brandt, Zack Briesemeister, Katherine de Kleer, Imke de Pater, Aditi Desai, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Rebecca Jensen-Clem, Isabel Kain, Charlie Kilpatrick, Renate Kupke, Mackenzie Lach, Michael C. Liu, Bruce Macintosh, Raquel A. Martinez, Dimitri Mawet , et al. (12 additional authors not shown)

Abstract: The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~μ$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~μ$m… ▽ More The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~μ$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~μ$m) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES' unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements, and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in $\sim2025$. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 10 pages, 16 figures, submitted to Proceedings of the SPIE

The development of HISPEC for Keck and MODHIS for TMT: science cases and predicted sensitivities

Authors: Quinn M. Konopacky, Ashley D. Baker, Dimitri Mawet, Michael P. Fitzgerald, Nemanja Jovanovic, Charles Beichman, Garreth Ruane, Rob Bertz, Hiroshi Terada, Richard Dekany, Larry Lingvay, Marc Kassis, David Anderson, Motohide Tamura, Bjorn Benneke, Thomas Beatty, Tuan Do, Shogo Nishiyama, Peter Plavchan, Jason Wang, Ji Wang, Adam Burgasser, Jean-Baptiste Ruffio, Huihao Zhang, Aaron Brown , et al. (50 additional authors not shown)

Abstract: HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measureme… ▽ More HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 25 pages, 9 figures. To appear in the Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

CoCCoA: Complex Chemistry in hot Cores with ALMA. Selected oxygen-bearing species

Authors: Y. Chen, M. L. van Gelder, P. Nazari, C. L. Brogan, E. F. van Dishoeck, H. Linnartz, J. K. Jørgensen, T. R. Hunter, O. H. Wilkins, G. A. Blake, P. Caselli, K. -J. Chuang, C. Codella, I. Cooke, M. N. Drozdovskaya, R. T. Garrod, S. Ioppolo, M. Jin, B. M. Kulterer, N. F. W. Ligterink, A. Lipnicky, R. Loomis, M. G. Rachid, S. Spezzano, B. A. McGuire

Abstract: Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protos… ▽ More Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protostars using ALMA are still lacking. With the recent CoCCoA survey, we are able to determine the column density ratios of six O-COMs with respect to methanol (CH$_3$OH) in a sample of 14 high-mass protostellar sources to investigate their origin through ice and/or gas-phase chemistry. The selected species are: acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), dimethyl ether (DME, CH$_3$OCH$_3$), methyl formate (MF, CH$_3$OCHO), glycolaldehyde (GA, CH$_2$OHCHO), and ethylene glycol (EG, (CH$_2$OH)$_2$). DME and MF have the highest and most constant ratios within one order of magnitude, while the other four species have lower ratios and exhibit larger scatter by 1-2 orders of magnitude. We compare the O-COM ratios of high-mass CoCCoA sources with those of 5 low-mass protostars available from the literature, along with the results from experiments and simulations. We find that the O-COM ratios with respect to methanol are on the same level in both the high- and low-mass samples, which suggests that these species are mainly formed in similar environments during star formation, probably in ice mantles on dust grains during early pre-stellar stages. Current simulations and experiments can reproduce most observational trends with a few exceptions, and hypotheses exist to explain the differences between observations and simulations/experiments, such as the involvement of gas-phase chemistry and different emitting areas of molecules. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 17 pages and 7 figures in the main body, 15 pages and 8 figures in the appendix

Journal ref: A&A 678, A137 (2023)

Keck/KPIC Emission Spectroscopy of WASP-33b

Authors: Luke Finnerty, Tobias Schofield, Ben Sappey, Jerry W. Xuan, Jean-Baptiste Ruffio, Jason J. Wang, Jacques-Robert Delorme, Geoffrey A. Blake, Cam Buzard, Michael P. Fitzgerald, Ashley Baker, Randall Bartos, Charlotte Z. Bond, Benjamin Calvin, Sylvain Cetre, Greg Doppmann, Daniel Echeverri, Nemanja Jovanovic, Joshua Liberman, Ronald A. Lopez, Emily C. Martin, Dimitri Mawet, Evan Morris, Jacklyn Pezzato, Caprice L. Phillips , et al. (7 additional authors not shown)

Abstract: We present Keck/KPIC high-resolution ($R\sim35,000$) $K$-band thermal emission spectroscopy of the ultra-hot Jupiter WASP-33b. The use of KPIC's single-mode fibers greatly improves both blaze and line-spread stabilities relative to slit spectrographs, enhancing the cross-correlation detection strength. We retrieve the dayside emission spectrum with a nested sampling pipeline which fits for orbital… ▽ More We present Keck/KPIC high-resolution ($R\sim35,000$) $K$-band thermal emission spectroscopy of the ultra-hot Jupiter WASP-33b. The use of KPIC's single-mode fibers greatly improves both blaze and line-spread stabilities relative to slit spectrographs, enhancing the cross-correlation detection strength. We retrieve the dayside emission spectrum with a nested sampling pipeline which fits for orbital parameters, the atmospheric pressure-temperature profile, and molecular abundances.We strongly detect the thermally-inverted dayside and measure mass-mixing ratios for CO ($\log\rm CO_{MMR} = -1.1^{+0.4}_{-0.6}$), H$_2$O ($\log\rm H_2O_{MMR} = -4.1^{+0.7}_{-0.9}$) and OH ($\log\rm OH_{MMR} = -2.1^{+0.5}_{-1.1}$), suggesting near-complete dayside photodissociation of H$_2$O. The retrieved abundances suggest a carbon- and possibly metal-enriched atmosphere, with a gas-phase C/O ratio of $0.8^{+0.1}_{-0.2}$, consistent with the accretion of high-metallicity gas near the CO$_2$ snow line and post-disk migration or with accretion between the soot and H$_2$O snow lines. We also find tentative evidence for $\rm ^{12}CO/^{13}CO \sim 50$, consistent with values expected in protoplanetary disks, as well as tentative evidence for a metal-enriched atmosphere (2--15$\times$ solar). These observations demonstrate KPIC's ability to characterize close-in planets and the utility of KPIC's improved instrumental stability for cross-correlation techniques. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted in AJ, 26 pages, 12 figures

Exoplanet Volatile Carbon Content as a Natural Pathway for Haze Formation

Authors: E. A. Bergin, E. Kempton, M. Hirschmann, S. T. Bastelberger, D. J. Teal, G. A. Blake, F. Ciesla, J. Li

Abstract: We explore terrestrial planet formation with a focus on the supply of solid-state organics as the main source of volatile carbon. For the water-poor Earth, the water ice line, or ice sublimation front, within the planet-forming disk has long been a key focal point. We posit that the soot line, the location where solid-state organics are irreversibly destroyed, is also a key location within the dis… ▽ More We explore terrestrial planet formation with a focus on the supply of solid-state organics as the main source of volatile carbon. For the water-poor Earth, the water ice line, or ice sublimation front, within the planet-forming disk has long been a key focal point. We posit that the soot line, the location where solid-state organics are irreversibly destroyed, is also a key location within the disk. The soot line is closer to the host star than the water snowline and overlaps with the location of the majority of detected exoplanets. In this work, we explore the ultimate atmospheric composition of a body that receives a major portion of its materials from the zone between the soot line and water ice line. We model a silicate-rich world with 0.1% and 1% carbon by mass with variable water content. We show that as a result of geochemical equilibrium, the mantle of these planets would be rich in reduced carbon but have relatively low water (hydrogen) content. Outgassing would naturally yield the ingredients for haze production when exposed to stellar UV photons in the upper atmosphere. Obscuring atmospheric hazes appear common in the exoplanetary inventory based on the presence of often featureless transmission spectra (Kreidberg et al. 2014, Knutson et al. 2014, Libby-Roberts et al. 2020). Such hazes may be powered by the high volatile content of the underlying silicate-dominated mantle. Although this type of planet has no solar system counterpart, it should be common in the galaxy with potential impact on habitability. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 15 pages, 6 figures, accepted by Astrophysical Journal Letters. ApJL version has 5 figures. We include an extra figure (Figure 6) in this submission which is an artist rendering of a young disk including the soot and water ice lines. Image Credit for Fig. 6: Ari Gea/SayoStudio

Breaking Degeneracies in Formation Histories by Measuring Refractory Content in Gas Giants

Authors: Yayaati Chachan, Heather A. Knutson, Joshua Lothringer, Geoffrey A. Blake

Abstract: Relating planet formation to atmospheric composition has been a long-standing goal of the planetary science community. So far, most modeling studies have focused on predicting the enrichment of heavy elements and the C/O ratio in giant planet atmospheres. Although this framework provides useful constraints on the potential formation locations of gas giant exoplanets, carbon and oxygen measurements… ▽ More Relating planet formation to atmospheric composition has been a long-standing goal of the planetary science community. So far, most modeling studies have focused on predicting the enrichment of heavy elements and the C/O ratio in giant planet atmospheres. Although this framework provides useful constraints on the potential formation locations of gas giant exoplanets, carbon and oxygen measurements alone are not enough to determine where a given gas giant planet originated. Here, we show that characterizing the abundances of refractory elements (e.g., silicon, iron) can break these degeneracies. Refractory elements are present in the solid phase throughout most of the disk and their atmospheric abundances therefore reflect the solid-to-gas accretion ratio during formation. We introduce a new framework that parameterizes the atmospheric abundances of gas giant exoplanets in the form of three ratios: Si/H, O/Si, and C/Si. Si/H traces the solid-to-gas accretion ratio of a planet and is loosely equivalent to earlier notions of 'metallicity'. For O/Si and C/Si, we present a global picture of their variation with distance and time based on what we know from the solar system meteorites and an updated understanding of the variations of thermal processing within protoplanetary disks. We show that ultra-hot Jupiters are ideal targets for atmospheric characterization studies using this framework, as we can measure the abundances of refractories, oxygen and carbon in the gas phase. Finally, we propose that hot Jupiters with silicate clouds and low water abundances might have accreted their envelopes between the soot line and the water snowline. △ Less

Submitted 1 February, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: accepted to ApJ on Nov 22, 2022. Added DOI and journal ref

Journal ref: ApJ, 943, 112 (2023)

Retrieving C and O Abundance of HR 8799 c by Combining High- and Low-Resolution Data

Authors: Ji Wang, Jason J. Wang, Jean-Baptiste Ruffio, Geoffrey A. Blake, Dimitri Mawet, Ashley Baker, Randall Bartos, Charlotte Z. Bond, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Nemanja Jovanovic, Ronald Lopez, Emily C. Martin, Evan Morris, Jacklyn Pezzato, Sam Ragland, Garreth Ruane, Ben Sappey, Tobias Schofield, Andrew Skemer , et al. (7 additional authors not shown)

Abstract: The formation and evolution pathway for the directly-imaged multi-planetary system HR 8799 remains mysterious. Accurate constraints on the chemical composition of the planetary atmosphere(s) are key to solving the mystery. We perform a detailed atmospheric retrieval on HR 8799~c to infer the chemical abundances and abundance ratios using a combination of photometric data along with low- and high-r… ▽ More The formation and evolution pathway for the directly-imaged multi-planetary system HR 8799 remains mysterious. Accurate constraints on the chemical composition of the planetary atmosphere(s) are key to solving the mystery. We perform a detailed atmospheric retrieval on HR 8799~c to infer the chemical abundances and abundance ratios using a combination of photometric data along with low- and high-resolution spectroscopic data (R$\sim$20-35,000). We specifically retrieve [C/H], [O/H], and C/O and find them to be 0.55$^{+0.36}_{-0.39}$, 0.47$^{+0.31}_{-0.32}$, and 0.67$^{+0.12}_{-0.15}$ at 68\% confidence. The super-stellar C and O abundances, yet a stellar C/O ratio, reveal a potential formation pathway for HR 8799~c. Planet c, and likely the other gas giant planets in the system, formed early on (likely within $\sim$1 Myr), followed by further atmospheric enrichment in C and O through the accretion of solids beyond the CO iceline. The enrichment either preceded or took place during the early phase of the inward migration to the planet current locations. △ Less

Submitted 26 October, 2022; v1 submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 6 figures, 3 tables, accepted to AAS journals

The kinematics and excitation of infrared water vapor emission from planet-forming disks: results from spectrally-resolved surveys and guidelines for JWST spectra

Authors: Andrea Banzatti, Klaus M. Pontoppidan, José Pérez Chávez, Colette Salyk, Lindsey Diehl, Simon Bruderer, Greg J. Herczeg, Andres Carmona, Ilaria Pascucci, Sean Brittain, Stanley Jensen, Sierra Grant, Ewine F. van Dishoeck, Inga Kamp, Arthur D. Bosman, Karin I. Öberg, Geoff A. Blake, Michael R. Meyer, Eric Gaidos, Adwin Boogert, John T. Rayner, Caleb Wheeler

Abstract: This work presents ground-based spectrally-resolved water emission at R = 30000-100000 over infrared wavelengths covered by JWST (2.9-12.8 $μ$m). Two new surveys with iSHELL and VISIR are combined with previous spectra from CRIRES and TEXES to cover parts of multiple ro-vibrational and rotational bands observable within telluric transmission bands, for a total of $\approx160$ spectra and 85 disks… ▽ More This work presents ground-based spectrally-resolved water emission at R = 30000-100000 over infrared wavelengths covered by JWST (2.9-12.8 $μ$m). Two new surveys with iSHELL and VISIR are combined with previous spectra from CRIRES and TEXES to cover parts of multiple ro-vibrational and rotational bands observable within telluric transmission bands, for a total of $\approx160$ spectra and 85 disks (30 of which are JWST targets in Cycle 1). The general expectation of a range of regions and excitation conditions traced by infrared water spectra is for the first time supported by the combined kinematics and excitation as spectrally resolved at multiple wavelengths. The main findings from this analysis are: 1) water lines are progressively narrower from the ro-vibrational bands at 2-9 $μ$m to the rotational lines at 12 $μ$m, and partly match a broad (BC) and narrow (NC) emission components, respectively, as extracted from ro-vibrational CO spectra; 2) rotation diagrams of resolved water lines from upper level energies of 4000-9500 K show vertical spread and curvatures indicative of optically thick emission ($\approx 10^{18}$ cm$^{-2}$) from a range of excitation temperatures ($\approx 800$-1100 K); 3) the new 5 $μ$m spectra demonstrate that slab model fits to the rotational lines at $> 10$ $μ$m strongly over-predict the ro-vibrational emission bands at $< 9$ $μ$m, implying non-LTE vibrational excitation. We discuss these findings in the context of emission from a disk surface and a molecular inner disk wind, and provide a list of guidelines to support the analysis of spectrally-unresolved JWST spectra. △ Less

Submitted 16 November, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Accepted for publication on AJ

An unusual reservoir of water emission in the VV CrA A protoplanetary disk

Authors: Colette Salyk, Klaus M. Pontoppidan, Andrea Banzatti, Ulrich Käufl, Cassandra Hall, Ilaria Pascucci, Andrés Carmona, Geoffrey A. Blake, Richard Alexander, Inga Kamp

Abstract: We present an analysis of an unusual pattern of water vapor emission from the $\sim$2 Myr-old low-mass binary system VV CrA, as observed in infrared spectra obtained with VLT-CRIRES, VLT-VISIR, and Spitzer-IRS. Each component of the binary shows emission from water vapor in both the L ($\sim3\,μ$m) and N ($\sim 12\,μ$m) bands. The N-band and Spitzer spectra are similar to those previously observed… ▽ More We present an analysis of an unusual pattern of water vapor emission from the $\sim$2 Myr-old low-mass binary system VV CrA, as observed in infrared spectra obtained with VLT-CRIRES, VLT-VISIR, and Spitzer-IRS. Each component of the binary shows emission from water vapor in both the L ($\sim3\,μ$m) and N ($\sim 12\,μ$m) bands. The N-band and Spitzer spectra are similar to those previously observed from young stars with disks, and are consistent with emission from an extended protoplanetary disk. Conversely, the CRIRES L-band data of VV CrA A show an unusual spectrum, which requires the presence of a water reservoir with high temperature ($T\gtrsim1500$ K), column density ($N_\mathrm{H2O}\sim 3\times10^{20}\ \mathrm{cm}^{-2}$), and turbulent broadening ($v\sim 10$ km s$^{-1}$), but very small emitting area ($A\lesssim0.005$ AU$^2$). Similarity with previously observed water emission from V1331 Cyg (Doppmann et al. 2011) and SVS 13 (Carr et al. 2004) suggests that the presence of such a reservoir may be linked to evolutionary state, perhaps related to the presence of high accretion rates or winds. While the inner disk may harbor such a reservoir, simple Keplerian models do not match well with emitting line shapes, and alternative velocity fields must be considered. We also present a new idea, that the unusual emission could arise in a circumplanetary disk, embedded within the larger VV CrA A protoplanetary disk. Additional data are likely required to determine the true physical origin of this unusual spectral pattern. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 32 pages, 17 figures, 5 appendix figures, 2 tables. Accepted for publication in ApJ

On the Relationship between CH$_3$OD Abundance and Temperature in the Orion KL Nebula

Authors: Olivia H. Wilkins, Geoffrey A. Blake

Abstract: The relative abundances of singly-deuterated methanol isotopologues, [CH$_{2}$DOH]/[CH$_{3}$OD], in star-forming regions deviate from the statistically expected ratio of 3. In Orion KL, the nearest high-mass star-forming region to Earth, the singly-deuterated methanol ratio is about 1, and the cause for this observation has been explored through theory for nearly three decades. We present high-ang… ▽ More The relative abundances of singly-deuterated methanol isotopologues, [CH$_{2}$DOH]/[CH$_{3}$OD], in star-forming regions deviate from the statistically expected ratio of 3. In Orion KL, the nearest high-mass star-forming region to Earth, the singly-deuterated methanol ratio is about 1, and the cause for this observation has been explored through theory for nearly three decades. We present high-angular resolution observations of Orion KL using the Atacama Large Millimeter/submillimeter Array to map small-scale changes in CH$_{3}$OD column density across the nebula, which provide a new avenue to examine the deuterium chemistry during star and planet formation. By considering how CH$_{3}$OD column densities vary with temperature, we find evidence of chemical processes that can significantly alter the observed column densities. The astronomical data are compared with existing theoretical work and support D-H exchange between CH$_{3}$OH and heavy water (i.e., HDO and D$_{2}$O) at methanol's hydroxyl site in the icy mantles of dust grains. The enhanced CH$_{3}$OD column densities are localized to the Hot Core-SW region, a pattern that may be linked to the coupled evolution of ice mantel chemistry and star formation in giant molecular clouds. This work provides new perspectives on deuterated methanol chemistry in Orion KL and informs considerations that may guide future theoretical, experimental, and observational work. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted for publication in the Journal of Physical Chemistry A

A Clear View of a Cloudy Brown Dwarf Companion from High-Resolution Spectroscopy

Authors: Jerry W. Xuan, Jason Wang, Jean-Baptiste Ruffio, Heather Knutson, Dimitri Mawet, Paul Mollière, Jared Kolecki, Arthur Vigan, Sagnick Mukherjee, Nicole Wallack, Ji Wang, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Charlotte Z. Bond, Marta Bryan, Benjamin Calvin, Sylvain Cetre, Mark Chun, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Katelyn Horstman , et al. (15 additional authors not shown)

Abstract: Direct imaging studies have mainly used low-resolution spectroscopy ($R\sim20-100$) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g. C/O, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer… ▽ More Direct imaging studies have mainly used low-resolution spectroscopy ($R\sim20-100$) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g. C/O, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer (KPIC) uses adaptive optics and single-mode fibers to transport light into NIRSPEC ($R\sim35,000$ in $K$ band), and aims to address these challenges with high-resolution spectroscopy. Using an atmospheric retrieval framework based on petitRADTRANS, we analyze KPIC high-resolution spectrum ($2.29-2.49~μ$m) and archival low-resolution spectrum ($1-2.2~μ$m) of the benchmark brown dwarf HD 4747 B ($m=67.2\pm1.8~M_{\rm{Jup}}$, $a=10.0\pm0.2$ au, $T_{\rm eff}\approx1400$ K). We find that our measured C/O and metallicity for the companion from the KPIC high-resolution spectrum agree with that of its host star within $1-2σ$. The retrieved parameters from the $K$ band high-resolution spectrum are also independent of our choice of cloud model. In contrast, the retrieved parameters from the low-resolution spectrum are highly sensitive to our chosen cloud model. Finally, we detect CO, H$_2$O, and CH$_4$ (volume mixing ratio of log(CH$_4$)=$-4.82\pm0.23$) in this L/T transition companion with the KPIC data. The relative molecular abundances allow us to constrain the degree of chemical disequilibrium in the atmosphere of HD 4747 B, and infer a vertical diffusion coefficient that is at the upper limit predicted from mixing length theory. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 33 pages, 16 figures, Accepted to ApJ

Airfall on Comet 67P/Churyumov-Gerasimenko

Authors: B. J. R. Davidsson, S. Birch, G. A. Blake, D. Bodewits, J. P. Dworkin, D. P. Glavin, Y. Furukawa, J. I. Lunine, J. L. Mitchell, A. N. Nguyen, S. Squyres, A. Takigawa, J. -B. Vincent, K. Zacny

Abstract: We here study the transfer process of material from one hemisphere to the other (deposition of airfall material) on an active comet nucleus, specifically 67P/Churyumov-Gerasimenko. Our goals are to: 1) quantify the thickness of the airfall debris layers and how it depends on the location of the target area, 2) determine the amount of $\mathrm{H_2O}$ and $\mathrm{CO_2}$ ice that are lost from icy d… ▽ More We here study the transfer process of material from one hemisphere to the other (deposition of airfall material) on an active comet nucleus, specifically 67P/Churyumov-Gerasimenko. Our goals are to: 1) quantify the thickness of the airfall debris layers and how it depends on the location of the target area, 2) determine the amount of $\mathrm{H_2O}$ and $\mathrm{CO_2}$ ice that are lost from icy dust assemblages of different sizes during transfer through the coma, and 3) estimate the relative amount of vapor loss in airfall material after deposition in order to understand what locations are expected to be more active than others on the following perihelion approach. We use various numerical simulations, that include orbit dynamics, thermophysics of the nucleus and of individual coma aggregates, coma gas kinetics and hydrodynamics, as well as dust dynamics due to gas drag, to address these questions. We find that the thickness of accumulated airfall material varies substantially with location, and typically is of the order $0.1$-$1\,\mathrm{m}$. The airfall material preserves substantial amounts of water ice even in relatively small (cm-sized) coma aggregates after a rather long ($12\,\mathrm{h}$) residence in the coma. However, $\mathrm{CO_2}$ is lost within a couple of hours even in relatively large (dm-sized) aggregates, and is not expected to be an important component in airfall deposits. We introduce reachability and survivability indices to measure the relative capacity of different regions to simultaneously collect airfall and to preserve its water ice until the next perihelion passage, thereby grading their potential of contributing to comet activity during the next perihelion passage. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 65 pages, 11 figures. Published manuscript

Journal ref: Icarus 354, 114004 (2021)

Near-infrared spectra of the inflated post-common envelope brown dwarf NLTT 5306B

Authors: Cam Buzard, Sarah Casewell, Joshua Lothringer, Geoffrey Blake

Abstract: NLTT 5306 is a post-common envelope binary made up of a white dwarf host and brown dwarf companion that has shown evidence of inflation and active mass donation despite not filling its Roche lobe. Two proposed mechanisms for the brown dwarf's inflation are magnetic interactions and a high metallicity, cloudy atmosphere. We present moderate resolution ($R\lesssim2000$) $J$ band Keck/NIRSPEC observa… ▽ More NLTT 5306 is a post-common envelope binary made up of a white dwarf host and brown dwarf companion that has shown evidence of inflation and active mass donation despite not filling its Roche lobe. Two proposed mechanisms for the brown dwarf's inflation are magnetic interactions and a high metallicity, cloudy atmosphere. We present moderate resolution ($R\lesssim2000$) $J$ band Keck/NIRSPEC observations of this system. These phase-resolved data allow us to constrain differences between atmospheric parameters of the day- and night-side of the brown dwarf. Our day- and night-side effective temperature measurements are consistent, in agreement with the brightness temperatures measurements from Casewell et al. 2021a. The day-side favors a slightly lower surface gravity, perhaps stemming from the material streaming between the two objects. Finally, our data show a preference for low metallicity models. This would be expected from the system's old age, but provides direct evidence that a high metallicity, cloudy brown dwarf atmosphere is not responsible for the witnessed inflation. These results strengthen the case for magnetic interactions leading to inflation of NLTT 5306 B. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 12 pages, 7 figures, Accepted for publication in AJ

Retrieving the C and O Abundances of HR 7672~AB: a Solar-Type Primary Star with a Benchmark Brown Dwarf

Authors: Ji Wang, Jared R. Kolecki, Jean-Baptiste Ruffio, Jason J. Wang, Dimitri Mawet, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Charlotte Z. Bond, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Nemanja Jovanovic, Michael C. Liu, Ronald Lopez, Evan Morris, Anusha Pai Asnodkar, Jacklyn Pezzato, Sam Ragland, Arpita Roy, Garreth Ruane , et al. (8 additional authors not shown)

Abstract: A benchmark brown dwarf (BD) is a BD whose properties (e.g., mass and chemical composition) are precisely and independently measured. Benchmark BDs are valuable in testing theoretical evolutionary tracks, spectral synthesis, and atmospheric retrievals for sub-stellar objects. Here, we report results of atmospheric retrieval on a synthetic spectrum and a benchmark BD -- HR 7672~B -- with \petit. Fi… ▽ More A benchmark brown dwarf (BD) is a BD whose properties (e.g., mass and chemical composition) are precisely and independently measured. Benchmark BDs are valuable in testing theoretical evolutionary tracks, spectral synthesis, and atmospheric retrievals for sub-stellar objects. Here, we report results of atmospheric retrieval on a synthetic spectrum and a benchmark BD -- HR 7672~B -- with \petit. First, we test the retrieval framework on a synthetic PHOENIX BT-Settl spectrum with a solar composition. We show that the retrieved C and O abundances are consistent with solar values, but the retrieved C/O is overestimated by 0.13-0.18, which is $\sim$4 times higher than the formal error bar. Second, we perform retrieval on HR 7672~B using high spectral resolution data (R=35,000) from the Keck Planet Imager and Characterizer (KPIC) and near infrared photometry. We retrieve [C/H], [O/H], and C/O to be $-0.24\pm0.05$, $-0.19\pm0.04$, and $0.52\pm0.02$. These values are consistent with those of HR 7672~A within 1.5-$σ$. As such, HR 7672~B is among only a few benchmark BDs (along with Gl 570~D and HD 3651~B) that have been demonstrated to have consistent elemental abundances with their primary stars. Our work provides a practical procedure of testing and performing atmospheric retrieval, and sheds light on potential systematics of future retrievals using high- and low-resolution data. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 29 pages, 17 figures, 5 tables, resubmitted to AAS journals after first revision

New Constraints on Protoplanetary Disk Gas Masses in Lupus

Authors: Dana E. Anderson, L. Ilsedore Cleeves, Geoffrey A. Blake, Edwin A. Bergin, Ke Zhang, John M. Carpenter, Kamber R. Schwarz

Abstract: Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H$_2$ content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologu… ▽ More Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H$_2$ content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologues of CO. However, without additional constraints, the degeneracy between H$_2$ mass and the elemental composition of the gas leads to large uncertainties in such estimates. Here we explore the gas compositions of seven disks from the Lupus sample representing a range of CO-to-dust ratios. With Band 6 and 7 ALMA observations, we measure line emission for HCO$^+$, HCN, and N$_2$H$^+$. We find a tentative correlation among the line fluxes for these three molecular species across the sample, but no correlation with $^{13}$CO or sub-mm continuum fluxes. For the three disks where N$_2$H$^+$ is detected, we find that a combination of high disk gas masses and sub-interstellar C/H and O/H are needed to reproduce the observed values. We find increases of $\sim$10-100$\times$ previous mass estimates are required to match the observed line fluxes. This study highlights how multi-molecular studies are essential for constraining the physical and chemical properties of the gas in populations of protoplanetary disks and that CO isotopologues alone are not sufficient for determining the mass of many observed disks. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: Accepted for publication in The Astrophysical Journal, 19 pages, 15 figures, 2 tables

Mapping physical parameters in Orion KL at high spatial resolution

Authors: Olivia H. Wilkins, P. Brandon Carroll, Geoffrey A. Blake

Abstract: The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high angular resolution (sub-arcsecond) observations of $^13$CH$_{3}$OH $ν=0$ ($\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$) and CH$_3$CN $ν_8=1$ ($\sim$0.2$^{\prime\prime}$ and $\sim$0.9$^{\prime\prime}$) line emission with the Atacama Large Millimet… ▽ More The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high angular resolution (sub-arcsecond) observations of $^13$CH$_{3}$OH $ν=0$ ($\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$) and CH$_3$CN $ν_8=1$ ($\sim$0.2$^{\prime\prime}$ and $\sim$0.9$^{\prime\prime}$) line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) to investigate Orion KL's structure on small spatial scales (${\le}350$ au). Gas kinematics, excitation temperatures, and column densities were derived from the molecular emission via a pixel-by-pixel spectral line fitting of the image cubes, enabling us to examine the small-scale variation of these parameters. Sub-regions of the Hot Core have a higher excitation temperature in a 0.2$^{\prime\prime}$ beam than a 0.9$^{\prime\prime}$ beam, indicative of possible internal sources of heating. Furthermore, the velocity field includes a bipolar ${\sim}7{-}8$ km s$^{-1}$ feature with a southeast-northwest orientation against the surrounding ${\sim}4{-}5$ km s$^{-1}$ velocity field, which may be due to an outflow. We also find evidence of a possible source of internal heating toward the Northwest Clump, since the excitation temperature there is higher in a smaller beam versus a larger beam. Finally, the region southwest of the Hot Core (Hot Core-SW) presents itself as a particularly heterogeneous region bridging the Hot Core and Compact Ridge. Additional studies to identify the (hidden) sources of luminosity and heating within Orion KL are necessary to better understand the nebula and its chemistry. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 28 pages, 15 figures, to be published in ApJ

Reinvestigation of the Multi-Epoch Direct Detections of HD 88133 b and Upsilon Andromedae b

Authors: Cam Buzard, Danielle Piskorz, Alexandra C. Lockwood, Geoffrey Blake, Travis S. Barman, Björn Benneke, Chad F. Bender, John S. Carr

Abstract: We reanalyze the multi-epoch direct detections of HD 88133 b and ups And b that were published in Piskorz et al. 2016 and Piskorz et al. 2017, respectively. Using simulations to attempt to reproduce the detections, we find that with the 6 and 7 $L$ band Keck/NIRSPEC epochs analyzed in the original works, the planets would not have been detectable unless they had unreasonably large radii. HD88133 a… ▽ More We reanalyze the multi-epoch direct detections of HD 88133 b and ups And b that were published in Piskorz et al. 2016 and Piskorz et al. 2017, respectively. Using simulations to attempt to reproduce the detections, we find that with the 6 and 7 $L$ band Keck/NIRSPEC epochs analyzed in the original works, the planets would not have been detectable unless they had unreasonably large radii. HD88133 and ups And both have fairly large stellar radii, which contributed to the difficulty in detecting the planets. We take this opportunity to consider how these planets may have been detectable with the small number of epochs originally presented by running simulations both with the upgraded NIRSPEC instrument and with near-zero primary velocities, as recommended by Buzard et al. 2021. While 7 $L$ band NIRSPEC2.0 epochs with near-zero primary velocities would have allowed a strong ($10.8σ$) detection of ups And b, many more than 6 $L$ band epochs would have been required for a strong detection of HD88133b, which could be due in part to both this system's large stellar radius and low stellar temperature. This work stresses the importance of careful analytic procedures and the usefulness of simulations in understanding the expected sensitivity of high-resolution spectroscopic data. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 11 pages, 8 figures, Accepted for publication in AJ

Detection and Bulk Properties of the HR 8799 Planets with High Resolution Spectroscopy

Authors: Jason J. Wang, Jean-Baptiste Ruffio, Evan Morris, Jacques-Robert Delorme, Nemanja Jovanovic, Jacklyn Pezzato, Daniel Echeverri, Luke Finnerty, Callie Hood, J. J. Zanazzi, Marta L. Bryan, Charlotte Z. Bond, Sylvain Cetre, Emily C. Martin, Dimitri Mawet, Andy Skemer, Ashley Baker, Jerry W. Xuan, J. Kent Wallace, Ji Wang, Randall Bartos, Geoffrey A. Blake, Andy Boden, Cam Buzard, Benjamin Calvin , et al. (27 additional authors not shown)

Abstract: Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral reso… ▽ More Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured $v\sin(i)$ values of $10.1^{+2.8}_{-2.7}$~km/s for HR 8799 d and $15.0^{+2.3}_{-2.6}$~km/s for HR 8799 e, and placed an upper limit of $< 14$~km/s of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anti-correlated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower mass planets. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 31 pages, 12 figures, Accepted to AJ

Primary Velocity and Orbital Phase Effects on Planetary Detectability from Small Epoch Number Data Sets

Authors: Cam Buzard, Stefan Pelletier, Danielle Piskorz, Björn Benneke, Geoffrey Blake

Abstract: Cross correlation analyses of high resolution spectroscopic data have recently shown great success in directly detecting planetary signals and enabling the characterization of their atmospheres. One such technique aims to observe a system at multiple epochs and combine the measured planetary radial velocities from each epoch into a measurement of the planetary Keplerian orbital velocity $K_p$, con… ▽ More Cross correlation analyses of high resolution spectroscopic data have recently shown great success in directly detecting planetary signals and enabling the characterization of their atmospheres. One such technique aims to observe a system at multiple epochs and combine the measured planetary radial velocities from each epoch into a measurement of the planetary Keplerian orbital velocity $K_p$, constituting a direct detection of the planetary signal. Recent work has shown that in few-epoch ($\sim$5) data sets, unintended structure can arise at a high level, obscuring the planetary detection. In this work, we look to simulations to examine whether there are ways to reduce this structured noise in few-epoch data sets by careful planning of observations. The choice of observation date allows observers to select the primary (stellar) velocity - through a set systemic velocity and chosen barycentric velocity - and the planetary orbital phase, and so we focus on the effects of these two parameters. We find that epochs taken when the primary velocity is near zero, and the stellar lines remain relatively fixed to the telluric rest-frame, greatly reduce the level of structured noise and allow for much stronger planetary detections, on average more than twice the significance of detections made with epochs using randomly selected primary velocities. Following these results, we recommend that observers looking to build up high-resolution multi-epoch data sets target nights when their system has a near-zero primary velocity. △ Less

Submitted 15 April, 2021; originally announced April 2021.

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

Early volatile depletion on planetesimals inferred from C-S systematics of iron meteorite parent bodies

Authors: Marc M. Hirschmann, Edwin A. Bergin, Geoff A. Blake, Fred J. Ciesla, Jie Li

Abstract: During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partition… ▽ More During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration and inferred parent body C concentrations range from 0.0004 to 0.11 wt.%. Parent bodies fall into 2 compositional clusters characterized by cores with medium, and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed system extraction to degassing of a wholly molten body, show that significant open system silicate melting and volatile loss is required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss, and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 50 pages including supplementary materials

Journal ref: Publications of the National Academy of Sciences, March 30, 2021 118 (13)

Earth's carbon deficit caused by early loss through irreversible sublimation

Authors: Jie Li, Edwin A. Bergin, Geoffrey A. Blake, Fred J. Ciesla, Marc M. Hirschmann

Abstract: Carbon is an essential element for life but its behavior during Earth's accretion is not well understood. Carbonaceous grains in meteoritic and cometary materials suggest that irreversible sublimation, and not condensation, governs carbon acquisition by terrestrial worlds. Through astronomical observations and modeling we show that the sublimation front of carbon carriers in the solar nebula, or t… ▽ More Carbon is an essential element for life but its behavior during Earth's accretion is not well understood. Carbonaceous grains in meteoritic and cometary materials suggest that irreversible sublimation, and not condensation, governs carbon acquisition by terrestrial worlds. Through astronomical observations and modeling we show that the sublimation front of carbon carriers in the solar nebula, or the soot line, moved inward quickly so that carbon-rich ingredients would be available for accretion at 1 au after the first million years. On the other hand, geological constraints firmly establish a severe carbon deficit in Earth, requiring the destruction of inherited carbonaceous organics in the majority of its building blocks. The carbon-poor nature of the Earth thus implies carbon loss in its precursor material through sublimation within the first million years. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 21 pages including main article and supplementary materials

Journal ref: Science Advances published 02 Apr 2021: Vol. 7, no. 14, DOI: 10.1126/sciadv.abd3632

The TW Hya Rosetta Stone Project IV: A hydrocarbon rich disk atmosphere

Authors: L. Ilsedore Cleeves, Ryan A. Loomis, Richard Teague, Edwin A. Bergin, David J. Wilner, Jennifer B. Bergner, Geoffrey A. Blake, Jenny K. Calahan, Paolo Cazzoletti, Ewine F. van Dishoeck, Viviana V. Guzman, Michiel R. Hogerheijde, Jane Huang, Mihkel Kama, Karin I. Oberg, Chunhua Qi, Jeroen Terwisscha van Scheltinga, Catherine Walsh

Abstract: Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project… ▽ More Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at $15-20$ au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ comes from a layer near the surface that extends no deeper than $z/r=0.25$. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, Accepted in ApJ

Observing Carbon & Oxygen Carriers in Protoplanetary Disks at Mid-infrared Wavelengths

Authors: D. E. Anderson, G. A. Blake, L. I. Cleeves, E. A. Bergin, K. Zhang, K. R. Schwarz, C. Salyk, A. D. Bosman

Abstract: Infrared observations probe the warm gas in the inner regions of planet-forming disks around young sun-like, T Tauri stars. In these systems, H$_2$O, OH, CO, CO$_2$, C$_2$H$_2$, and HCN have been widely observed. However, the potentially abundant carbon carrier CH$_4$ remains largely unconstrained. The James Webb Space Telescope (JWST) will be able to characterize mid-infrared fluxes of CH$_4$ alo… ▽ More Infrared observations probe the warm gas in the inner regions of planet-forming disks around young sun-like, T Tauri stars. In these systems, H$_2$O, OH, CO, CO$_2$, C$_2$H$_2$, and HCN have been widely observed. However, the potentially abundant carbon carrier CH$_4$ remains largely unconstrained. The James Webb Space Telescope (JWST) will be able to characterize mid-infrared fluxes of CH$_4$ along with several other carriers of carbon and oxygen. In anticipation of the JWST mission, we model the physical and chemical structure of a T Tauri disk to predict the abundances and mid-infrared fluxes of observable molecules. A range of compositional scenarios are explored involving the destruction of refractory carbon materials and alterations to the total elemental (volatile and refractory) C/O ratio. Photon-driven chemistry in the inner disk surface layers largely destroys the initial carbon and oxygen carriers. This causes models with the same physical structure and C/O ratio to have similar steady state surface compositions, regardless of the initial chemical abundances. Initial disk compositions are better preserved in the shielded inner disk midplane. The degree of similarity between the surface and midplane compositions in the inner disk will depend on the characteristics of vertical mixing at these radii. Our modeled fluxes of observable molecules respond sensitively to changes in the disk gas temperature, inner radius, and the total elemental C/O ratio. As a result, mid-infrared observations of disks will be useful probes of these fundamental disk parameters, including the C/O ratio, which can be compared to values determined for planetary atmospheres. △ Less

Submitted 8 January, 2021; originally announced January 2021.

Comments: Accepted for publication in The Astrophysical Journal

Contrast and Temperature Dependence of Multi-Epoch High-Resolution Cross-Correlation Exoplanet Spectroscopy

Authors: Luke Finnerty, Cam Buzard, Stefan Pelletier, Danielle Piskorz, Alexandra C. Lockwood, Chad F. Bender, Björn Benneke, Geoffrey A. Blake

Abstract: While high-resolution cross-correlation spectroscopy (HRCCS) techniques have proven effective at characterizing the atmospheres of transiting and non-transiting hot Jupiters, the limitations of these techniques are not well understood. We present a series of simulations of one HRCCS technique, which combines the cross-correlation functions from multiple epochs, to place temperature and contrast li… ▽ More While high-resolution cross-correlation spectroscopy (HRCCS) techniques have proven effective at characterizing the atmospheres of transiting and non-transiting hot Jupiters, the limitations of these techniques are not well understood. We present a series of simulations of one HRCCS technique, which combines the cross-correlation functions from multiple epochs, to place temperature and contrast limits on the accessible exoplanet population for the first time. We find that planets approximately Saturn-size and larger within $\sim$0.2 AU of a Sun-like star are likely to be detectable with current instrumentation in the $L$-band, a significant expansion compared with the previously-studied population. Cooler ($ \rm T_{eq} \leq 1000$ K) exoplanets are more detectable than suggested by their photometric contrast alone as a result of chemical changes which increase spectroscopic contrast. The $L$-band CH$_4$ spectrum of cooler exoplanets enables robust constraints on the atmospheric C/O ratio at $\rm T_{eq} \sim 900K$, which have proven difficult to obtain for hot Jupiters. These results suggest that the multi-epoch approach to HRCCS can detect and characterize exoplanet atmospheres throughout the inner regions of Sun-like systems with existing high-resolution spectrographs. We find that many epochs of modest signal-to-noise ($\rm S/N_{epoch} \sim 1500$) yield the clearest detections and constraints on C/O, emphasizing the need for high-precision near-infrared telluric correction with short integration times. △ Less

Submitted 27 December, 2020; originally announced December 2020.

Comments: 22 pages, 8 figures, accepted to AJ

The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

Authors: Jenny Calahan, Edwin Bergin, Ke Zhang, Richard Teague, Ilsedore Cleeves, Jennifer Bergner, Geoffrey A. Blake, Paolo Cazzoletti, Viviana Guzman, Michiel R. Hogerheijde, Jane Huang, Mihkel Kama, Ryan Loomis, Karin Oberg, Ewine F. van Dishoeck, Jeroen Terwisscha van Scheltinga, Catherine Walsh, David Wilner, Charlie Qi

Abstract: The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations o… ▽ More The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 28 pages, 12 main figures, 14 appendix figures. Accepted to ApJ

The TW Hya Rosetta Stone Project II: Spatially resolved emission of formaldehyde hints at low-temperature gas-phase formation

Authors: Jeroen Terwisscha van Scheltinga, Michiel R. Hogerheijde, L. Ilsedore Cleeves, Ryan A. Loomis, Catherine Walsh, Karin I. Öberg, Edwin A. Bergin, Jennifer B. Bergner, Geoffrey A. Blake, Jenny K. Calahan, Paolo Cazzoletti, Ewine F. van Dishoeck, Viviana V. Guzmán, Jane Huang, Mihkel Kama, Chunhua Qi, Richard Teague, David J. Wilner

Abstract: Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pa… ▽ More Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pathway dominates, gaseous H$_2$CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and $<3$ linked to cold formation in the solid-state. We present spatially resolved ALMA observations of multiple ortho- and para-H$_2$CO transitions in the TW Hya protoplanetary disk to test H$_2$CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of $33\pm2$ K, ($1.1\pm0.1)\times10^{12}$ cm$^{-2}$ and $25\pm2$ K, $(4.4\pm0.3)\times10^{11}$ cm$^{-2}$ for ortho- and para-H$_2$CO, respectively, and an OPR of $2.49\pm0.23$. A radially resolved analysis shows that the observed H$_2$CO emits mostly at rotational temperatures of 30-40 K, corresponding to a layer with $z/R\ge0.25$. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to $2.0\pm0.5$. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, lead us to speculate that gas-phase formation is responsible for the observed H$_2$CO across the TW Hya disk. △ Less

Submitted 13 November, 2020; originally announced November 2020.

The TW Hya Rosetta Stone Project I: Radial and vertical distributions of DCN and DCO+

Authors: Karin I. Oberg, L. Ilsedore Cleeves, Jennifer B. Bergner, Joseph Cavanaro, Richard Teague, Jane Huang, Ryan A. Loomis, Edwin A. Bergin, Geoffrey A. Blake, Jenny Calahan, Paolo Cazzoletti, Viviana Veloso Guzman, Michiel R. Hogerheijde, Mihkel Kama, Jeroen Terwisscha van Scheltinga, Chunhua Qi, Ewine van Dishoeck, Catherine Walsh, David J. Wilner

Abstract: Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project,… ▽ More Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30~au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. There appears to be an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature, or something peculiar to the old TW Hya disk. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: Accepted for publication in The Astronomical Journal

The Evolution of Disk Winds from a Combined Study of Optical and Infrared Forbidden Lines

Authors: I. Pascucci, A. Banzatti, U. Gorti, M. Fang, K. Pontoppidan, R. Alexander, G. Ballabio, S. Edwards, C. Salyk, G. Sacco, E. Flaccomio, G. A. Blake, A. Carmona, C. Hall, I. Kamp, H. U. Kaufl, G. Meeus, M. Meyer, T. Pauly, S. Steendam, M. Sterzik

Abstract: We analyze high-resolution (dv=<10km/s) optical and infrared spectra covering the [OI] 6300 angstrom and [NeII] 12.81 micron lines from a sample of 31 disks in different evolutionary stages. Following work at optical wavelengths, we use Gaussian profiles to fit the [NeII] lines and classify them into HVC (LVC) if the line centroid is more (less) blueshifted than 30 km/s with respect to the stellar… ▽ More We analyze high-resolution (dv=<10km/s) optical and infrared spectra covering the [OI] 6300 angstrom and [NeII] 12.81 micron lines from a sample of 31 disks in different evolutionary stages. Following work at optical wavelengths, we use Gaussian profiles to fit the [NeII] lines and classify them into HVC (LVC) if the line centroid is more (less) blueshifted than 30 km/s with respect to the stellar radial velocity. Unlike for the [OI] where a HVC is often accompanied by a LVC, all 17 sources with a [NeII] detection have either a HVC or a LVC. [NeII] HVCs are preferentially detected toward high accretors (Macc > 10$^{-8}$ Msun/yr) while LVCs are found in sources with low Macc, low [OI] luminosity, and large infrared spectral index (n13-31). Interestingly, the [NeII] and [OI] LVC luminosities display an opposite behaviour with n13-31: as the inner dust disk depletes (higher n13-31) the [NeII] luminosity increases while the [OI] weakens. The [NeII] and [OI] HVC profiles are generally similar with centroids and FWHMs showing the expected behaviour from shocked gas in micro-jets. In contrast, the [NeII] LVC profiles are typically more blueshifted and narrower than the [OI] profiles. The FWHM and centroid vs. disk inclination suggest that the [NeII] LVC predominantly traces unbound gas from a slow, wide-angle wind that has not lost completely the Keplerian signature from its launching region. We sketch an evolutionary scenario that could explain the combined [OI] and [NeII] results and includes screening of hard (~1keV) X-rays in inner, mostly molecular, MHD winds. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks

Authors: Ryan A. Loomis, Karin I. Öberg, Sean M. Andrews, Edwin Bergin, Jennifer Bergner, Geoffrey A. Blake, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, Romane Le Gal, Francois Menard, Jamila Pegues, Chunhua Qi, Catherine Walsh, Jonathan P. Williams, David J. Wilner

Abstract: The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We p… ▽ More The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 20 pages, 11 figures

Journal ref: ApJ 893, 2 (2020)

An evolutionary study of volatile chemistry in protoplanetary disks

Authors: Jennifer B. Bergner, Karin I. Oberg, Edwin A. Bergin, Sean M. Andrews, Geoffrey A. Blake, John M. Carpenter, L. Ilsedore Cleeves, Viviana V. Guzman, Jane Huang, Jes K. Jorgensen, Chunhua Qi, Kamber R. Schwarz, Jonathan P. Williams, David J. Wilner

Abstract: The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class… ▽ More The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: Accepted to ApJ

Carbon-grain sublimation: a new top-down component of protostellar chemistry

Authors: Merel L. R. van 't Hoff, Edwin A. Bergin, Jes K. Jørgensen, Geoffrey A. Blake

Abstract: Earth's carbon deficit has been an outstanding problem in our understanding of the formation of our Solar System. A possible solution would be the sublimation of carbon grains at the so-called soot line (~300 K) early in the planet-formation process. Here, we argue that the most likely signatures of this process are an excess of hydrocarbons and nitriles inside the soot line, and a higher excitati… ▽ More Earth's carbon deficit has been an outstanding problem in our understanding of the formation of our Solar System. A possible solution would be the sublimation of carbon grains at the so-called soot line (~300 K) early in the planet-formation process. Here, we argue that the most likely signatures of this process are an excess of hydrocarbons and nitriles inside the soot line, and a higher excitation temperature for these molecules compared to oxygen-bearing complex organics that desorb around the water snowline (~100 K). Such characteristics have been reported in the literature, for example, in Orion KL, although not uniformly, potentially due to differences in observational settings and analysis methods of different studies or related to the episodic nature of protostellar accretion. If this process is active, this would mean that there is a heretofore unknown component to the carbon chemistry during the protostellar phase that is acting from the top down - starting from the destruction of larger species - instead of from the bottom up from atoms. In the presence of such a top-down component, the origin of organic molecules needs to be re-explored. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJL. 14 pages, 2 figures, 1 table (4 pages)

Simulating the Multi-Epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-Transiting Hot Jupiter HD187123B

Authors: Cam Buzard, Luke Finnerty, Danielle Piskorz, Stefan Pelletier, Björn Benneke, Chad F. Bender, Alexandra C. Lockwood, Nicole L. Wallack, Olivia H. Wilkins, Geoffrey A. Blake

Abstract: We report the 6.5$σ$ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity $K_p$ of 53 $\pm$ 13 km/s. This high confidence detection is made using a multi-epoch, high resolution, cross correlation technique, and corresponds to a planetary mass of 1.4$^{+0.5}_{-0.3}$ $M_J$ and an orbital inclination of 21 $\pm$ 5$^{\circ}$. The technique works by treating the planet/st… ▽ More We report the 6.5$σ$ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity $K_p$ of 53 $\pm$ 13 km/s. This high confidence detection is made using a multi-epoch, high resolution, cross correlation technique, and corresponds to a planetary mass of 1.4$^{+0.5}_{-0.3}$ $M_J$ and an orbital inclination of 21 $\pm$ 5$^{\circ}$. The technique works by treating the planet/star system as a spectroscopic binary and obtaining high signal-to-noise, high resolution observations at multiple points across the planet's orbit to constrain the system's binary dynamical motion. All together, seven epochs of Keck/NIRSPEC $L$-band observations were obtained, with five before the instrument upgrade and two after. Using high resolution SCARLET planetary and PHOENIX stellar spectral models, along with a line-by-line telluric absorption model, we were able to drastically increase the confidence of the detection by running simulations that could reproduce, and thus remove, the non-random structured noise in the final likelihood space well. The ability to predict multi-epoch results will be extremely useful for furthering the technique. Here, we use these simulations to compare three different approaches to combining the cross correlations of high resolution spectra and find that the Zucker 2003 log(L) approach is least affected by unwanted planet/star correlation for our HD187123 data set. Furthermore, we find that the same total S/N spread across an orbit in many, lower S/N epochs rather than fewer, higher S/N epochs could provide a more efficient detection. This work provides a necessary validation of multi-epoch simulations which can be used to guide future observations and will be key to studying the atmospheres of further separated, non-transiting exoplanets. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: Accepted to AJ, 14 pages, 10 figures

An ALMA Survey of H$_2$CO in Protoplanetary Disks

Authors: Jamila Pegues, Karin I. Öberg, Jennifer B. Bergner, Ryan A. Loomis, Chunhua Qi, Romane Le Gal, L. Ilsedore Cleeves, Viviana V. Guzmán, Jane Huang, Jes K. Jørgensen, Sean M. Andrews, Geoffrey A. Blake, John M. Carpenter, Kamber R. Schwarz, Jonathan P. Williams, David J. Wilner

Abstract: H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and c… ▽ More H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and centrally-depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H$_2$CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H$_2$CO excitation temperatures and column densities for four disks with multiple H$_2$CO line detections. The temperatures are between 20-50K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H$_2$CO emission in disks is generally emerging from the warm molecular layer, with some contributions from the colder midplane. Applying the same H$_2$CO excitation temperatures to all disks in the survey, we find that H$_2$CO column densities span almost three orders of magnitude ($\sim 5 \times 10^{11} - 5 \times 10^{14} \mathrm{cm}^{-2}$). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H$_2$CO compared to T Tauri disks, possibly because of less CO freeze-out. More H$_2$CO observations towards Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H$_2$CO and other oxygen-bearing organics efficiently form during planet formation. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: 20 pages, 14 figures, 3 figure sets, 9 tables. Published in ApJ (February 2020)

Journal ref: The Astrophysical Journal, 890, 142 (2020)

Probing planet formation and disk substructures in the inner disk of Herbig Ae stars with CO rovibrational emission

Authors: Arthur D. Bosman, Andrea Banzatti, Simon Bruderer, Alexander G. G. M. Tielens, Geoffrey A. Blake, Ewine F. van Dishoeck

Abstract: [abridged]CO rovibrational lines are efficient probes of warm molecular gas and can give unique insights into the inner 10 AU of proto-planetary disks. Recent studies have found a relation between the ratio of lines originating from the second and first vibrationally excited state, denoted as $v2/v1$, and the emitting radius of CO. In disks around Herbig Ae stars the vibrational excitation is low… ▽ More [abridged]CO rovibrational lines are efficient probes of warm molecular gas and can give unique insights into the inner 10 AU of proto-planetary disks. Recent studies have found a relation between the ratio of lines originating from the second and first vibrationally excited state, denoted as $v2/v1$, and the emitting radius of CO. In disks around Herbig Ae stars the vibrational excitation is low when CO lines come from close to the star, and high when lines only probe gas at large radii (more than 5 AU). We aim to find explanations for the observed trends between CO vibrational ratio, emitting radii, and NIR excess, and identify their implications in terms of the physical and chemical structure of inner disks around Herbig stars. Slab models and full disk thermo chemical models are calculated. Simulated observations from the models are directly compared to the data. Broad CO lines with low vibrational ratios are best explained by a warm (400-1300 K) inner disk surface with gas-to-dust ratios below 1000; no CO is detected within/at the inner dust rim, due to dissociation at high temperatures. In contrast, explaining the narrow lines with high vibrational ratios requires an inner cavity of a least 5 AU in both dust and gas, followed by a cool (100-300 K) molecular gas reservoir with gas-to-dust ratios greater than 10000 at the cavity wall. In all cases the CO gas must be close to thermalization with the dust. The high gas-to-dust ratios needed to explain high $v2/v1$ in narrow CO lines for a subset of group I disks can naturally be interpreted as due to the dust traps that have been proposed to explain millimeter dust cavities. The broad lines seen in most group II objects indicate a very flat disk in addition to inner disk substructures within 10 AU that can be related to the substructures recently observed with ALMA. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 17 pages, 16 figures, accepted by A&A

Journal ref: A&A 631, A133 (2019)

Searches for Interstellar HCCSH and H$_2$CCS

Authors: Brett A. McGuire, Christopher N. Shingledecker, Eric R. Willis, Kin Long Kelvin Lee, Marie-Aline Martin-Drumel, Geoffrey A. Blake, Crystal L. Brogan, Andrew M. Burkhardt, Paola Caselli, Ko-Ju Chuang, Samer El-Abd, Todd R. Hunter, Sergio Ioppolo, Harold Linnartz, Anthony J. Remijan, Ci Xue, Michael C. McCarthy

Abstract: A long standing problem in astrochemistry is the inability of many current models to account for missing sulfur content. Many relatively simple species that may be good candidates to sequester sulfur have not been measured experimentally at the high spectral resolution necessary to enable radioastronomical identification. On the basis of new laboratory data, we report searches for the rotational l… ▽ More A long standing problem in astrochemistry is the inability of many current models to account for missing sulfur content. Many relatively simple species that may be good candidates to sequester sulfur have not been measured experimentally at the high spectral resolution necessary to enable radioastronomical identification. On the basis of new laboratory data, we report searches for the rotational lines in the microwave, millimeter, and sub-millimeter regions of the sulfur-containing hydrocarbon HCCSH. This simple species would appear to be a promising candidate for detection in space owing to the large dipole moment along its $b$-inertial axis, and because the bimolecular reaction between two highly abundant astronomical fragments (CCH and SH radicals) may be rapid. An inspection of multiple line surveys from the centimeter to the far-infrared toward a range of sources from dark clouds to high-mass star-forming regions, however, resulted in non-detections. An analogous search for the lowest-energy isomer, H$_2$CCS, is presented for comparison, and also resulted in non-detections. Typical upper limits on the abundance of both species relative to hydrogen are $10^{-9}$-$10^{-10}$. We thus conclude that neither isomer is a major reservoir of interstellar sulfur in the range of environments studied. Both species may still be viable candidates for detection in other environments or at higher frequencies, providing laboratory frequencies are available. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: Accepted in the Astrophysical Journal

High-resolution Infrared Spectrograph for Exoplanet Characterization with the Keck and Thirty Meter Telescopes

Authors: Dimitri Mawet, Michael Fitzgerald, Quinn Konopacky, Charles Beichman, Nemanja Jovanovic, Richard Dekany, David Hover, Eric Chisholm, David Ciardi, Etienne Artigau, Ravinder Banyal, Thomas Beatty, Bjorn Benneke, Geoffrey A. Blake, Adam Burgasser, Gabriela Canalizo, Guo Chen, Tuan Do, Greg Doppmann, Rene Doyon, Courtney Dressing, Min Fang, Thomas Greene, Lynne Hillenbrand, Andrew Howard , et al. (24 additional authors not shown)

Abstract: HISPEC (High-resolution Infrared Spectrograph for Exoplanet Characterization) is a proposed diffraction-limited spectrograph for the W.M. Keck Observatory, and a pathfinder for the MODHIS facility project (Multi-Object Diffraction-limited High-resolution Infrared Spectrograph) on the Thirty Meter Telescope. HISPEC/MODHIS builds on diffraction-limited spectrograph designs which rely on adaptively c… ▽ More HISPEC (High-resolution Infrared Spectrograph for Exoplanet Characterization) is a proposed diffraction-limited spectrograph for the W.M. Keck Observatory, and a pathfinder for the MODHIS facility project (Multi-Object Diffraction-limited High-resolution Infrared Spectrograph) on the Thirty Meter Telescope. HISPEC/MODHIS builds on diffraction-limited spectrograph designs which rely on adaptively corrected single-mode fiber feeds. Seeing-limited high-resolution spectrographs, by virtue of the conservation of beam etendue, grow in volume following a D^3 power law (D is the telescope diameter), and are subject to daunting challenges associated with their large size. Diffraction-limited spectrographs fed by single mode fibers are decoupled from the telescope input, and are orders of magnitude more compact and have intrinsically stable line spread functions. Their efficiency is directly proportional to the performance of the adaptive optics (AO) system. AO technologies have matured rapidly over the past two decades and are baselined for future extremely large telescopes. HISPEC/MODHIS will take R>100,000 spectra of a few objects in a 10" field-of-view sampled at the diffraction limit (~10-50 mas), simultaneously from 0.95 to 2.4 microns (y-K). The scientific scope ranges from exoplanet infrared precision radial velocities, spectroscopy of transiting, close-in, and directly imaged exoplanets (atmospheric composition and dynamics, RM effect, spin measurements, Doppler imaging), brown dwarf characterization, stellar physics/chemistry, proto-planetary disk kinematics/composition, Solar system, extragalactic science, and cosmology. HISPEC/MODHIS features a compact, cost-effective design optimized to fully exploit the existing Keck-AO and future TMT-NFIRAOS infrastructures and boost the scientific reach of Keck Observatory and TMT soon after first light. △ Less

Submitted 9 August, 2019; originally announced August 2019.

Journal ref: Astro2020 APC White Paper: Optical and Infrared Observations from the Ground

Probing CO and N$_2$ Snow Surfaces in Protoplanetary Disks with N$_2$H$^+$ Emission

Authors: Chunhua Qi, Karin I. Öberg, Catherine C. Espaillat, Connor E. Robinson, Sean M. Andrews, David J. Wilner, Geoffrey A. Blake, Edwin A. Bergin, L. Ilsedore Cleeves

Abstract: Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of b… ▽ More Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$''$-0.4$''$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: Accepted for publication in ApJ

Unlocking CO Depletion in Protoplanetary Disks II. Primordial C/H Predictions Inside the CO Snowline

Authors: Kamber R. Schwarz, Edwin A. Bergin, L. Ilsedore Cleeves, Ke Zhang, Karin I. Öberg, Geoffrey A. Blake, Dana E. Anderson

Abstract: CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemicall… ▽ More CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemically reprocessed into less volatile species, which are frozen on dust grain surfaces as ice. This mechanism has the potential to change the primordial C/H ratio in the gas. However, current observations primarily probe the upper layers of the disk. It is not clear if the low volatile carbon abundances extend to the midplane, where planets form. We have run a grid of 198 chemical models, exploring how the chemical reprocessing of CO depends on disk mass, dust grain size distribution, temperature, cosmic ray and X-ray ionization rate, and initial water abundance. Building on our previous work focusing on the warm molecular layer, here we analyze the results for our grid of models in the disk midplane at 12 au. We find that either an ISM level cosmic ray ionization rate or the presence of UV photons due to a low dust surface density are needed to chemically reduce the midplane CO gas abundance by at least an order of magnitude within 1 Myr. In the majority of our models CO does not undergo substantial reprocessing by in situ chemistry and there is little change in the gas phase C/H and C/O ratios over the lifetime of the typical disk. However, in the small sub-set of disks where the disk midplane is subject to a source of ionization or photolysis, the gas phase C/O ratio increases by up to nearly 9 orders of magnitude due to conversion of CO into volatile hydrocarbons. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: Accepted for publication in ApJ, 15 pages, 10 figures, 3 tables