Showing 1–50 of 283 results for author: Blake, A

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 11 figures, accepted for publication in PRD

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

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

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

Submitted 21 June, 2024; originally announced June 2024.

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

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

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

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

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted to ApJ, 52 pages, 20 figures

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

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

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

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

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

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

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

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

Submitted 17 April, 2024; originally announced April 2024.

Comments: 14 pages, 8 figures

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

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

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

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

Submitted 11 April, 2024; originally announced April 2024.

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

Molecular Gas Tracers in Young and Old Protoplanetary Disks

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

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

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astrophysical Journal

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

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

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

Submitted 30 November, 2023; originally announced December 2023.

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

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

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

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

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

Comments: Accepted for publication in the Astrophysical Journal

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

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

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

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

Comments: 16 pages, 8 figures, accepted by ApJL

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

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

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

Submitted 19 September, 2023; originally announced September 2023.

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

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

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

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

Submitted 4 August, 2023; originally announced August 2023.

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

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

Keck/KPIC Emission Spectroscopy of WASP-33b

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

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

Submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted in AJ, 26 pages, 12 figures

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Exoplanet Volatile Carbon Content as a Natural Pathway for Haze Formation

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

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

Submitted 8 May, 2023; originally announced May 2023.

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

Simulated recovery of LEO objects using sCMOS blind stacking

Authors: Benjamin F. Cooke, Paul Chote, Don Pollacco, Richard West, James A. Blake, James McCormac, Robert Airey, Billy Shrive

Abstract: We present the methodology and results of a simulation to determine the recoverability of LEO objects using a blind stacking technique. The method utilises sCMOS and GPU technology to inject and recover LEO objects in real observed data. We explore the target recovery fraction and pipeline run-time as a function of three optimisation parameters; number of frames per data-set, exposure time, and bi… ▽ More We present the methodology and results of a simulation to determine the recoverability of LEO objects using a blind stacking technique. The method utilises sCMOS and GPU technology to inject and recover LEO objects in real observed data. We explore the target recovery fraction and pipeline run-time as a function of three optimisation parameters; number of frames per data-set, exposure time, and binning factor. Results are presented as a function of magnitude and velocity. We find that target recovery using blind stacking is significantly more complete, and can reach fainter magnitudes, than using individual frames alone. We present results showing that, depending on the combination of optimisation parameters, recovery fraction is up to 90% of detectable targets for magnitudes up to 13.5, and then falls off steadily up to a magnitude limit around 14.5. Run-time is shown to be a few multiples of the observing time for the best combinations of optimisation parameters, approaching real-time processing. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: 14 pages, 14 figures. Accepted for publication in Advances in Space Research (ASR)

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

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

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

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

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

Journal ref: ApJ, 943, 112 (2023)

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

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

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

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

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

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

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

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

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

Comments: Accepted for publication on AJ

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

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

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

Submitted 22 August, 2022; originally announced August 2022.

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

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

Authors: Olivia H. Wilkins, Geoffrey A. Blake

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

Submitted 17 August, 2022; originally announced August 2022.

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

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

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

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

Submitted 2 August, 2022; originally announced August 2022.

Comments: 33 pages, 16 figures, Accepted to ApJ

Airfall on Comet 67P/Churyumov-Gerasimenko

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

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

Submitted 1 August, 2022; originally announced August 2022.

Comments: 65 pages, 11 figures. Published manuscript

Journal ref: Icarus 354, 114004 (2021)

Gusts in the Headwind: Uncertainties in Direct Dark Matter Detection

Authors: Grace E. Lawrence, Alan R. Duffy, Chris A. Blake, Philip F. Hopkins

Abstract: We use high-resolution, hydrodynamic, galaxy simulations from the Latte suite of FIRE-2 simulations to investigate the inherent variation of dark matter in sub-sampled regions around the Solar Circle of a Milky Way-type analogue galaxy and its impact on direct dark matter detection. These simulations show that the baryonic backreaction, as well as the assembly history of substructures, has lasting… ▽ More We use high-resolution, hydrodynamic, galaxy simulations from the Latte suite of FIRE-2 simulations to investigate the inherent variation of dark matter in sub-sampled regions around the Solar Circle of a Milky Way-type analogue galaxy and its impact on direct dark matter detection. These simulations show that the baryonic backreaction, as well as the assembly history of substructures, has lasting impacts on the dark matter's spatial and velocity distributions. These are experienced as 'gusts' of dark matter wind around the Solar Circle, potentially complicating interpretations of direct detection experiments on Earth. We find that the velocity distribution function in the galactocentric frame shows strong deviations from the Maxwell Boltzmann form typically assumed in the fiducial Standard Halo Model, indicating the presence of high-velocity substructures. By introducing a new numerical integration technique which removes any dependencies on the Standard Halo Model, we generate event-rate predictions for both single-element Germanium and compound Sodium Iodide detectors, and explore how the variability of dark matter around the Solar Circle influences annual modulation signal predictions. We find that these velocity substructures contribute additional astrophysical uncertainty to the interpretation of event rates, although their impact on summary statistics such as the peak day of annual modulation is generally low. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 22 pages, 27 figures, 4 tables. Submitted to MNRAS. Comments welcome

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

Looking out for a sustainable space

Authors: James A. Blake

Abstract: October 1957, and the successful launch of Sputnik 1 into Earth orbit, marked the dawn of the Space Age. The first of the 'fellow travellers' - humanity's first artificial satellite - orbited for a mere three months before re-entering the Earth's atmosphere, though its mission paved the way for an era of exploration that has endured to the present day. For many, a world without satellites would be… ▽ More October 1957, and the successful launch of Sputnik 1 into Earth orbit, marked the dawn of the Space Age. The first of the 'fellow travellers' - humanity's first artificial satellite - orbited for a mere three months before re-entering the Earth's atmosphere, though its mission paved the way for an era of exploration that has endured to the present day. For many, a world without satellites would be a difficult one to imagine. As a society, we have become reliant on them for a vast array of services and applications. With a divine view of large swathes of the Earth's surface, and the ability to relay signals around its curvature, satellites have enabled the fast transfer of data on a global scale, bypassing the challenges associated with ground-based broadcasting, long-distance wiring, and so on. Positioning, Navigation and Timing (PNT) satellites have revolutionised transportation by land, air, and sea, while weather satellites enable scientists to monitor and warn of large-scale phenomena as they develop in near real-time. Satellites have extended the frontiers of observation: looking outwards, astronomers are able to circumvent the Earth's atmosphere to look deeper into the cosmos than ever before; looking inwards, patterns and processes that feed into general circulation models can be monitored on a range of timescales, improving our understanding of climate change. Satellites, and the services they provide, are not going to disappear any time soon. That said, threats to satellite safety do exist, and it is important that they be addressed as soon as possible to avoid long-lasting damage to operations in the space domain. Nearly sixty-five years on from the advent of human activity in space, I chart the evolution of the orbital debris environment and review latest efforts to make space operations more sustainable. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 9 pages, 6 figures - Invited review for Astronomy & Geophysics (RAS Journals) - Accepted for publication in Volume 63, Issue 2

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

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

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

Submitted 4 February, 2022; originally announced February 2022.

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

New Constraints on Protoplanetary Disk Gas Masses in Lupus

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

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

Submitted 1 February, 2022; originally announced February 2022.

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

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

Abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre… ▽ More For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs. △ Less

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

Authors: R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, M. Hayashi , et al. (140 additional authors not shown)

Abstract: We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of t… ▽ More We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Indications of a Cosmic Ray Source in the Perseus-Pisces Supercluster

Authors: Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon , et al. (135 additional authors not shown)

Abstract: The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A… ▽ More The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with $20^\circ$-radius circles indicates an excess in the arrival direction of events with a local significance of about 4 standard deviations. The probability of having such excess close to the PPSC by chance is estimated to be 3.5 standard deviations. This result indicates that a cosmic ray source likely exists in that supercluster. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 8 pages, 4 figures, 1 table

Mapping physical parameters in Orion KL at high spatial resolution

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

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

Submitted 26 October, 2021; originally announced October 2021.

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

Searching for solar KDAR with DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti, M. P. Andrews , et al. (1157 additional authors not shown)

Abstract: The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search.… ▽ More The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 19 pages, 13 figures

Report number: FERMILAB-PUB-21-322-LBNF-ND

Journal ref: JCAP10(2021)065

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

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

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

Submitted 14 July, 2021; originally announced July 2021.

Comments: 31 pages, 12 figures, Accepted to AJ

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

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

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

Submitted 6 April, 2021; originally announced April 2021.

Comments: 50 pages including supplementary materials

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

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

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

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

Submitted 6 April, 2021; originally announced April 2021.

Comments: 21 pages including main article and supplementary materials

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

Surface detectors of the TAx4 experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (124 additional authors not shown)

Abstract: Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with… ▽ More Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with energies greater than 57 EeV. In order to confirm this evidence with more data, it is necessary to increase the data collection rate.We have begun building an expansion of TA that we call TAx4. In this paper, we explain the motivation, design, technical features, and expected performance of the TAx4 SD. We also present TAx4's current status and examples of the data that have already been collected. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 26 pages, 17 figures, submitted to Nuclear Inst. and Methods in Physics Research, A

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

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

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

Submitted 18 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, Accepted in ApJ

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

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

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

Submitted 8 January, 2021; originally announced January 2021.

Comments: Accepted for publication in The Astrophysical Journal

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

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

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

Submitted 27 December, 2020; originally announced December 2020.

Comments: 22 pages, 8 figures, accepted to AJ

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

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

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

Submitted 10 December, 2020; originally announced December 2020.

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

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

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

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

Submitted 13 November, 2020; originally announced November 2020.

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

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

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

Submitted 13 November, 2020; originally announced November 2020.

Comments: Accepted for publication in The Astronomical Journal

Observations of the Origin of Downward Terrestrial Gamma-Ray Flashes

Authors: J. W. Belz, P. R. Krehbiel, J. Remington, M. A. Stanley, R. U. Abbasi, R. LeVon, W. Rison, D. Rodeheffer, the Telescope Array Scientific Collaboration, :, T. Abu-Zayyad, M. Allen, E. Barcikowski, D. R. Bergman, S. A. Blake, M. Byrne, R. Cady, B. G. Cheon, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich , et al. (116 additional authors not shown)

Abstract: In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (I… ▽ More In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud-to-ground and low-altitude intracloud flashes, and that the IBPs are produced by a newly-identified streamer-based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark-like transient conducting events (TCEs) within the fast streamer system, and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub-pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub-pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP. △ Less

Submitted 12 October, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Typo fixed and reference added. Manuscript is 36 pages. Supplemental Information is 42 pages. This paper is to be published in the Journal of Geophysical Research: Atmospheres. Online data repository: Open Science Framework DOI: 10.17605/OSF.IO/Z3XDA

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

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

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

Submitted 18 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

DebrisWatch I: A survey of faint geosynchronous debris

Authors: James A. Blake, Paul Chote, Don Pollacco, William Feline, Grant Privett, Andrew Ash, Stuart Eves, Arthur Greenwood, Nick Harwood, Thomas R. Marsh, Dimitri Veras, Christopher Watson

Abstract: Recent anomalies exhibited by satellites and rocket bodies have highlighted that a population of faint debris exists at geosynchronous (GEO) altitudes, where there are no natural removal mechanisms. Despite previous optical surveys probing to around 10-20 cm in size, regular monitoring of faint sources at GEO is challenging, thus our knowledge remains sparse. It is essential that we continue to ex… ▽ More Recent anomalies exhibited by satellites and rocket bodies have highlighted that a population of faint debris exists at geosynchronous (GEO) altitudes, where there are no natural removal mechanisms. Despite previous optical surveys probing to around 10-20 cm in size, regular monitoring of faint sources at GEO is challenging, thus our knowledge remains sparse. It is essential that we continue to explore the faint debris population using large telescopes to better understand the risk posed to active GEO satellites. To this end, we present photometric results from a survey of the GEO region carried out with the 2.54 m Isaac Newton Telescope in La Palma, Canary Islands. We probe to 21st visual magnitude (around 10 cm, assuming Lambertian spheres with an albedo of 0.1), uncovering 129 orbital tracks with GEO-like motion across the eight nights of dark-grey time comprising the survey. The faint end of our brightness distribution continues to rise until the sensitivity limit of the sensor is reached, suggesting that the modal brightness could be even fainter. We uncover a number of faint, uncatalogued objects that show photometric signatures of rapid tumbling, many of which straddle the limiting magnitude of our survey over the course of a single exposure, posing a complex issue when estimating object size. This work presents the first instalment of DebrisWatch, an ongoing collaboration between the University of Warwick and the Defence Science and Technology Laboratory (UK) investigating the faint population of GEO debris. △ Less

Submitted 28 August, 2020; originally announced August 2020.

Comments: 24 pages, 8 figures. Accepted for publication in Adv. Space Res

Supernova Neutrino Burst Detection with the Deep Underground Neutrino Experiment

Authors: DUNE collaboration, B. Abi, R. Acciarri, M. A. Acero, G. Adamov, D. Adams, M. Adinolfi, Z. Ahmad, J. Ahmed, T. Alion, S. Alonso Monsalve, C. Alt, J. Anderson, C. Andreopoulos, M. P. Andrews, F. Andrianala, S. Andringa, A. Ankowski, M. Antonova, S. Antusch, A. Aranda-Fernandez, A. Ariga, L. O. Arnold, M. A. Arroyave, J. Asaadi , et al. (949 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The gen… ▽ More The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The general capabilities of DUNE for neutrino detection in the relevant few- to few-tens-of-MeV neutrino energy range will be described. As an example, DUNE's ability to constrain the $ν_e$ spectral parameters of the neutrino burst will be considered. △ Less

Submitted 29 May, 2021; v1 submitted 15 August, 2020; originally announced August 2020.

Comments: 29 pages, 17 figures; paper based on DUNE Technical Design Report. arXiv admin note: substantial text overlap with arXiv:2002.03005

Report number: FERMILAB-PUB-20-380-LBNF

Search for Large-scale Anisotropy on Arrival Directions of Ultra-high-energy Cosmic Rays Observed with the Telescope Array Experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (121 additional authors not shown)

Abstract: Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is… ▽ More Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger. △ Less

Submitted 27 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 6 pages, 3 figures, 1 table, Proofed title. Added journal reference and DOI

Journal ref: The Astrophysical Journal Letters 898, L28 (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.

Comments: 20 pages, 11 figures

Journal ref: ApJ 893, 2 (2020)