-
Retrieval of Thermally-Resolved Water Vapor Distributions in Disks Observed with JWST-MIRI
Authors:
Carlos E. Muñoz-Romero,
Andrea Banzatti,
Karin I. Öberg,
Klaus M. Pontoppidan,
Colette Salyk,
Joan Najita,
Geoffrey A. Blake,
Sebastiaan Krijt,
Nicole Arulanantham,
Paola Pinilla,
Feng Long,
Giovanni Rosotti,
Sean M. Andrews,
David J. Wilner,
Jenny Calahan,
The JDISCS Collaboration
Abstract:
The mid-infrared water vapor emission spectrum provides a novel way to characterize the delivery of icy pebbles towards the innermost ($<5$ au) regions of planet-forming disks. Recently, JWST MIRI-MRS showed that compact disks exhibit an excess of low-energy water vapor emission relative to extended multi-gapped disks, suggesting that icy pebble drift is more efficient in the former. We carry out…
▽ More
The mid-infrared water vapor emission spectrum provides a novel way to characterize the delivery of icy pebbles towards the innermost ($<5$ au) regions of planet-forming disks. Recently, JWST MIRI-MRS showed that compact disks exhibit an excess of low-energy water vapor emission relative to extended multi-gapped disks, suggesting that icy pebble drift is more efficient in the former. We carry out detailed emission line modeling to retrieve the excitation conditions of rotational water vapor emission in a sample of four compact and three extended disks within the JDISC Survey. We present two-temperature H$_2$O slab model retrievals and, for the first time, constrain the spatial distribution of water vapor by fitting parametric radial temperature and column density profiles. Such models statistically outperform the two-temperature slab fits. We find a correlation between the observable hot water vapor mass and stellar mass accretion rate, as well as an anti-correlation between cold water vapor mass and sub-mm dust disk radius, confirming previously reported water line flux trends. We find that the mid-IR spectrum traces H$_2$O with temperatures down to 180-300 K, but the coldest 150-170 K gas remains undetected. Furthermore the H$_2$O temperature profiles are generally steeper and cooler than the expected `super-heated' dust temperature in passive irradiated disks. The column density profiles are used to estimate icy pebble mass fluxes, which suggest that compact and extended disks may produce markedly distinct inner-disk exoplanet populations if local feeding mechanisms dominate their assembly.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
RV measurements of directly imaged brown dwarf GQ Lup B to search for exo-satellites
Authors:
Katelyn Horstman,
Jean-Baptiste Ruffio,
Konstantin Batygin,
Dimitri Mawet,
Ashley Baker,
Chih-Chun Hsu,
Jason J. Wang,
Ji Wang,
Sarah Blunt,
Jerry W. Xuan,
Yinzi Xin,
Joshua Liberman,
Shubh Agrawal,
Quinn M. Konopacky,
Geoffrey A. Blake,
Clarissa R. Do O,
Randall Bartos,
Charlotte Z. Bond,
Benjamin Calvin,
Sylvain Cetre,
Jacques-Robert Delorme,
Greg Doppmann,
Daniel Echeverri,
Luke Finnerty,
Michael P. Fitzgerald
, et al. (13 additional authors not shown)
Abstract:
GQ Lup B is one of the few substellar companions with a detected cicumplanetary disk, or CPD. Observations of the CPD suggest the presence of a cavity, possibly formed by an exo-satellite. Using the Keck Planet Imager and Characterizer (KPIC), a high contrast imaging suite that feeds a high resolution spectrograph (1.9-2.5 microns, R$\sim$35,000), we present the first dedicated radial velocity (RV…
▽ More
GQ Lup B is one of the few substellar companions with a detected cicumplanetary disk, or CPD. Observations of the CPD suggest the presence of a cavity, possibly formed by an exo-satellite. Using the Keck Planet Imager and Characterizer (KPIC), a high contrast imaging suite that feeds a high resolution spectrograph (1.9-2.5 microns, R$\sim$35,000), we present the first dedicated radial velocity (RV) observations around a high-contrast, directly imaged substellar companion, GQ Lup B, to search for exo-satellites. Over 11 epochs, we find a best and median RV error of 400-1000 m/s, most likely limited by systematic fringing in the spectra due to transmissive optics within KPIC. With this RV precision, KPIC is sensitive to exomoons 0.6-2.8% the mass of GQ Lup B ($\sim 30 M_{\text{Jup}}$) at separations between the Roche limit and $65 R_{\text{Jup}}$, or the extent of the cavity inferred within the CPD detected around GQ Lup B. Using simulations of HISPEC, a high resolution infrared spectrograph planned to debut at W.M. Keck Observatory in 2026, we estimate future exomoon sensitivity to increase by over an order of magnitude, providing sensitivity to less massive satellites potentially formed within the CPD itself. Additionally, we run simulations to estimate the amount of material that different masses of satellites could clear in a CPD to create the observed cavity. We find satellite-to-planet mass ratios of $q > 2 \times 10^{-4}$ can create observable cavities and report a maximum cavity size of $\sim 51 \, R_{\text{Jup}}$ carved from a satellite.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
Fringing analysis and forward modeling of Keck Planet Imager and Characterizer (KPIC) spectra
Authors:
Katelyn A. Horstman,
Jean-Baptiste Ruffio,
Jason J. Wang,
Chih-Chun Hsu,
Ashley Baker,
Luke Finnerty,
Jerry Xuan,
Daniel Echeverri,
Dimitri Mawet,
Geoffrey A. Blake,
Randall Bartos,
Charlotte Z. Bond,
Benjamin Calvin,
Sylvain Cetre,
Jacques-Robert Delorme,
Greg Doppmann,
Michael P. Fitzgerald,
Nemanja Jovanovic,
Ronald Lopez,
Emily C. Martin,
Evan Morris,
Jacklyn Pezzato,
Garreth Ruane,
Ben Sappey,
Tobias Schofield
, et al. (5 additional authors not shown)
Abstract:
The Keck Planet Imager and Characterizer (KPIC) combines high contrast imaging with high resolution spectroscopy (R$\sim$35,000 in K band) to study directly imaged exoplanets and brown dwarfs in unprecedented detail. KPIC aims to spectrally characterize substellar companions through measurements of planetary radial velocities, spins, and atmospheric composition. Currently, the dominant source of s…
▽ More
The Keck Planet Imager and Characterizer (KPIC) combines high contrast imaging with high resolution spectroscopy (R$\sim$35,000 in K band) to study directly imaged exoplanets and brown dwarfs in unprecedented detail. KPIC aims to spectrally characterize substellar companions through measurements of planetary radial velocities, spins, and atmospheric composition. Currently, the dominant source of systematic noise for KPIC is fringing, or oscillations in the spectrum as a function of wavelength. The fringing signal can dominate residuals by up to 10% of the continuum for high S/N exposures, preventing accurate wavelength calibration, retrieval of atmospheric parameters, and detection of planets with flux ratios less than 1% of the host star. To combat contamination from fringing, we first identify its three unique sources and adopt a physically informed model of Fabry-Pérot cavities to apply to post-processed data. We find this strategy can effectively model the fringing in observations of A0V/F0V stars, reducing the residual systematics caused by fringing by a factor of 2. Next, we wedge two of the transmissive optics internal to KPIC to eliminate two sources of fringing and confirm the third source as the entrance window to the spectrograph. Finally, we apply our previous model of the Fabry-Pérot cavity to new data taken with the wedged optics to reduce the amplitude of the residuals by a factor of 10.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
Atmospheric characterization of the super-Jupiter HIP 99770 b with KPIC
Authors:
Yapeng Zhang,
Jerry W. Xuan,
Dimitri Mawet,
Jason J. Wang,
Chih-Chun Hsu,
Jean-Bapiste Ruffio,
Heather A. Knutson,
Julie Inglis,
Geoffrey A. Blake,
Yayaati Chachan,
Katelyn Horstman,
Ashley Baker,
Randall Bartos,
Benjamin Calvin,
Sylvain Cetre,
Jacques-Robert Delorme,
Greg Doppmann,
Daniel Echeverri,
Luke Finnerty,
Michael P. Fitzgerald,
Nemanja Jovanovic,
Joshua Liberman,
Ronald A. López,
Evan Morris,
Jacklyn Pezzato
, et al. (6 additional authors not shown)
Abstract:
Young, self-luminous super-Jovian companions discovered by direct imaging provide a challenging test of planet formation and evolution theories. By spectroscopically characterizing the atmospheric compositions of these super-Jupiters, we can constrain their formation histories. Here we present studies of the recently discovered HIP 99770 b, a 16 MJup high-contrast companion on a 17 au orbit, using…
▽ More
Young, self-luminous super-Jovian companions discovered by direct imaging provide a challenging test of planet formation and evolution theories. By spectroscopically characterizing the atmospheric compositions of these super-Jupiters, we can constrain their formation histories. Here we present studies of the recently discovered HIP 99770 b, a 16 MJup high-contrast companion on a 17 au orbit, using the fiber-fed high-resolution spectrograph KPIC (R~35,000) on the Keck II telescope. Our K-band observations led to detections of H2O and CO in the atmosphere of HIP 99770 b. We carried out free retrieval analyses using petitRADTRANS to measure its chemical abundances, including the metallicity and C/O ratio, projected rotation velocity (vsini), and radial velocity (RV). We found that the companion's atmosphere has C/O=0.55(-0.04/+0.06) and [M/H]=0.26(-0.23/+0.24) (1σ confidence intervals), values consistent with those of the Sun and with a companion formation via gravitational instability or core accretion. The projected rotation velocity < 7.8 km/s is small relative to other directly imaged companions with similar masses and ages. This may imply a near pole-on orientation or effective magnetic braking by a circumplanetary disk. In addition, we added the companion-to-primary relative RV measurement to the orbital fitting and obtained updated constraints on orbital parameters. Detailed characterization of super-Jovian companions within 20 au like HIP 99770 b is critical for understanding the formation histories of this population.
△ Less
Submitted 30 July, 2024;
originally announced July 2024.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Anisotropic Lattice Expansion of Monolayer WSe2 Revealed by Ultrafast Electron Diffraction
Authors:
Tony E. Karam,
Geoffrey A. Blake
Abstract:
Bulk layered MX2 transition metal chalcogenides (M = Mo, W and X = S, Se) are known to exhibit an indirect to direct band gap transition as the number of layers is reduced. Previous time-resolved work has principally focused on the investigation of the transient evolution of the band structure after photo-excitation, but additional information on the dynamics of the concomitant lattice rearrangeme…
▽ More
Bulk layered MX2 transition metal chalcogenides (M = Mo, W and X = S, Se) are known to exhibit an indirect to direct band gap transition as the number of layers is reduced. Previous time-resolved work has principally focused on the investigation of the transient evolution of the band structure after photo-excitation, but additional information on the dynamics of the concomitant lattice rearrangement is needed to fully understand this phenomenon. Here, ultrafast electron diffraction is used to probe the atomic motion and bond dilation in monolayer WSe2 with femtosecond temporal resolution. The change in the intensity of the Bragg diffraction spots is characterized by single-exponential dynamics, consistent with a collective response of the lattice during electron-phonon and phonon-phonon equilibration that is repeatable over many hours of illumination with ultrafast pulses. Moreover, a transient and highly anisotropic lattice expansion is observed. A possible explanation for this behavior is axial strain induced by the laser excitation that breaks the degeneracy of the in-plane phonon modes, and that strongly influences the electronic band structure. Such degeneracy-breaking induced by distortions in the lattice are well characterized in static strain measurements, and the results presented here provide valuable insights into the nature and time scales of the structural rearrangement that occurs following optical photo-excitation in monolayer tungsten diselenide.
△ Less
Submitted 5 February, 2023;
originally announced February 2023.
-
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.
-
Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality
Authors:
Tomislav Begušić,
Geoffrey A. Blake
Abstract:
Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals…
▽ More
Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.
△ Less
Submitted 20 April, 2023; v1 submitted 31 October, 2022;
originally announced October 2022.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy
Authors:
Tomislav Begušić,
Xuecheng Tao,
Geoffrey A. Blake,
Thomas F. Miller III
Abstract:
Two-dimensional Raman and hybrid terahertz/Raman spectroscopic techniques provide invaluable insight into molecular structure and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer mole…
▽ More
Two-dimensional Raman and hybrid terahertz/Raman spectroscopic techniques provide invaluable insight into molecular structure and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.
△ Less
Submitted 15 March, 2022; v1 submitted 2 February, 2022;
originally announced February 2022.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Sum-Frequency Signals in 2D-Terahertz-Terahertz-Raman Spectroscopy
Authors:
Griffin Mead,
Haw-Wei Lin,
Ioan-Bogdan Magdău,
Thomas F. Miller III,
Geoffrey A. Blake
Abstract:
We demonstrate that halogenated methane 2D-Terahertz Terahertz Raman (2D-TTR) spectra are determined by the complicated structure of the instrument response function (IRF) along $f_1$ and by the molecular coherences along $f_2$. Experimental improvements have helped increase the resolution and dynamic range of the measurements, including accurate THz pulse shape characterization. Sum-frequency exc…
▽ More
We demonstrate that halogenated methane 2D-Terahertz Terahertz Raman (2D-TTR) spectra are determined by the complicated structure of the instrument response function (IRF) along $f_1$ and by the molecular coherences along $f_2$. Experimental improvements have helped increase the resolution and dynamic range of the measurements, including accurate THz pulse shape characterization. Sum-frequency excitations convolved with the IRF are found to quantitatively reproduce the 2D-TTR signal. A new Reduced Density Matrix model which incorporates sum-frequency pathways, with linear and harmonic operators fully supports this (re)interpretation of the 2D-TTR spectra.
△ Less
Submitted 16 July, 2020;
originally announced July 2020.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.