Showing 1–50 of 96 results for author: Berger, T

Revising Properties of Planet-Host Binary Systems. IV. The Radius Distribution of Small Planets in Binary Star Systems is Dependent on Stellar Separation

Authors: Kendall Sullivan, Adam L. Kraus, Travis A. Berger, Trent J. Dupuy, Elise Evans, Eric Gaidos, Daniel Huber, Michael J. Ireland, Andrew W. Mann, Erik A. Petigura, Pa Chia Thao, Mackenna L. Wood, Jingwen Zhang

Abstract: Small planets ($R_{p} \leq 4 R_{\oplus}$) are divided into rocky super-Earths and gaseous sub-Neptunes separated by a radius gap, but the mechanisms that produce these distinct planet populations remain unclear. Binary stars are the only main-sequence systems with an observable record of the protoplanetary disk lifetime and mass reservoir, and the demographics of planets in binaries may provide in… ▽ More Small planets ($R_{p} \leq 4 R_{\oplus}$) are divided into rocky super-Earths and gaseous sub-Neptunes separated by a radius gap, but the mechanisms that produce these distinct planet populations remain unclear. Binary stars are the only main-sequence systems with an observable record of the protoplanetary disk lifetime and mass reservoir, and the demographics of planets in binaries may provide insights into planet formation and evolution. To investigate the radius distribution of planets in binary star systems, we observed 207 binary systems hosting 283 confirmed and candidate transiting planets detected by the Kepler mission, then recharacterized the planets while accounting for the observational biases introduced by the secondary star. We found that the population of planets in close binaries ($ρ\leq 100$ au) is significantly different from the planet population in wider binaries ($ρ> 300$ au) or single stars. In contrast to planets around single stars, planets in close binaries appear to have a unimodal radius distribution with a peak near the expected super-Earth peak of $R_{p} \sim 1.3 R_{\oplus}$ and a suppressed population of sub-Neptunes. We conclude that we are observing the direct impact of a reduced disk lifetime, smaller mass reservoir, and possible altered distribution of solids reducing the sub-Neptune formation efficiency. Our results demonstrate the power of binary stars as a laboratory for exploring planet formation and as a controlled experiment of the impact of varied initial conditions on mature planet populations. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 19 pages, 9 figures, 6 tables. Accepted to AJ. Full tables available upon request to the first author prior to publication

PHANGS-JWST: Data Processing Pipeline and First Full Public Data Release

Authors: Thomas G. Williams, Janice C. Lee, Kirsten L. Larson, Adam K. Leroy, Karin Sandstrom, Eva Schinnerer, David A. Thilker, Francesco Belfiore, Oleg V. Egorov, Erik Rosolowsky, Jessica Sutter, Joseph DePasquale, Alyssa Pagan, Travis A. Berger, Gagandeep S. Anand, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Sinan Deger, Cosima Eibensteiner , et al. (33 additional authors not shown)

Abstract: The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as… ▽ More The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of $\sim$5-50pc. Here, we describe pjpipe, an image processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.9" (degrading the resolution by about 30 percent) reduces the noise by a factor of $\approx$3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2-21 $μ$m imaging for all nineteen galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program. △ Less

Submitted 9 May, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 49 pages (27 in Appendices), 54 Figures (39 in Appendices), 3 Tables. Accepted for publication in ApJS. Updated to match accepted version. Data available at https://archive.stsci.edu/hlsp/phangs/phangs-jwst

High-Cadence Thermospheric Density Estimation enabled by Machine Learning on Solar Imagery

Authors: Shreshth A. Malik, James Walsh, Giacomo Acciarini, Thomas E. Berger, Atılım Güneş Baydin

Abstract: Accurate estimation of thermospheric density is critical for precise modeling of satellite drag forces in low Earth orbit (LEO). Improving this estimation is crucial to tasks such as state estimation, collision avoidance, and re-entry calculations. The largest source of uncertainty in determining thermospheric density is modeling the effects of space weather driven by solar and geomagnetic activit… ▽ More Accurate estimation of thermospheric density is critical for precise modeling of satellite drag forces in low Earth orbit (LEO). Improving this estimation is crucial to tasks such as state estimation, collision avoidance, and re-entry calculations. The largest source of uncertainty in determining thermospheric density is modeling the effects of space weather driven by solar and geomagnetic activity. Current operational models rely on ground-based proxy indices which imperfectly correlate with the complexity of solar outputs and geomagnetic responses. In this work, we directly incorporate NASA's Solar Dynamics Observatory (SDO) extreme ultraviolet (EUV) spectral images into a neural thermospheric density model to determine whether the predictive performance of the model is increased by using space-based EUV imagery data instead of, or in addition to, the ground-based proxy indices. We demonstrate that EUV imagery can enable predictions with much higher temporal resolution and replace ground-based proxies while significantly increasing performance relative to current operational models. Our method paves the way for assimilating EUV image data into operational thermospheric density forecasting models for use in LEO satellite navigation processes. △ Less

Submitted 12 November, 2023; originally announced December 2023.

Comments: Accepted at the Machine Learning and the Physical Sciences workshop, NeurIPS 2023

Two Warm Super-Earths Transiting the Nearby M Dwarf TOI-2095

Authors: Elisa V. Quintana, Emily A. Gilbert, Thomas Barclay, Michele L. Silverstein, Joshua E. Schlieder, Ryan Cloutier, Samuel N. Quinn, Joseph E. Rodriguez, Andrew Vanderburg, Benjamin J. Hord, Dana R. Louie, Colby Ostberg, Stephen R. Kane, Kelsey Hoffman, Jason F. Rowe, Giada N. Arney, Prabal Saxena, Taran Richardson, Matthew S. Clement, Nicholas M. Kartvedt, Fred C. Adams, Marcus Alfred, Travis Berger, Allyson Bieryla, Paul Bonney , et al. (33 additional authors not shown)

Abstract: We report the detection and validation of two planets orbiting TOI-2095 (TIC 235678745). The host star is a 3700K M1V dwarf with a high proper motion. The star lies at a distance of 42 pc in a sparsely populated portion of the sky and is bright in the infrared (K=9). With data from 24 Sectors of observation during TESS's Cycles 2 and 4, TOI-2095 exhibits two sets of transits associated with super-… ▽ More We report the detection and validation of two planets orbiting TOI-2095 (TIC 235678745). The host star is a 3700K M1V dwarf with a high proper motion. The star lies at a distance of 42 pc in a sparsely populated portion of the sky and is bright in the infrared (K=9). With data from 24 Sectors of observation during TESS's Cycles 2 and 4, TOI-2095 exhibits two sets of transits associated with super-Earth-sized planets. The planets have orbital periods of 17.7 days and 28.2 days and radii of 1.30 and 1.39 Earth radii, respectively. Archival data, preliminary follow-up observations, and vetting analyses support the planetary interpretation of the detected transit signals. The pair of planets have estimated equilibrium temperatures of approximately 400 K, with stellar insolations of 3.23 and 1.73 times that of Earth, placing them in the Venus zone. The planets also lie in a radius regime signaling the transition between rock-dominated and volatile-rich compositions. They are thus prime targets for follow-up mass measurements to better understand the properties of warm, transition radius planets. The relatively long orbital periods of these two planets provide crucial data that can help shed light on the processes that shape the composition of small planets orbiting M dwarfs. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Submitted to AAS Journals

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Revising Properties of Planet-Host Binary Systems. III. There is No Observed Radius Gap For Kepler Planets in Binary Star Systems

Authors: Kendall Sullivan, Adam L. Kraus, Daniel Huber, Erik A. Petigura, Elise Evans, Trent Dupuy, Jingwen Zhang, Travis A. Berger, Eric Gaidos, Andrew W. Mann

Abstract: Binary stars are ubiquitous; the majority of solar-type stars exist in binaries. Exoplanet occurrence rate is suppressed in binaries, but some multiples do still host planets. Binaries cause observational biases in planet parameters, with undetected multiplicity causing transiting planets to appear smaller than they truly are. We have analyzed the properties of a sample of 119 planet-host binary s… ▽ More Binary stars are ubiquitous; the majority of solar-type stars exist in binaries. Exoplanet occurrence rate is suppressed in binaries, but some multiples do still host planets. Binaries cause observational biases in planet parameters, with undetected multiplicity causing transiting planets to appear smaller than they truly are. We have analyzed the properties of a sample of 119 planet-host binary stars from the Kepler mission to study the underlying population of planets in binaries that fall in and around the radius valley, which is a demographic feature in period-radius space that marks the transition from predominantly rocky to predominantly gaseous planets. We found no statistically significant evidence for a radius gap for our sample of 122 planets in binaries when assuming the primary stars are the planet hosts, with a low probability ($p < 0.05$) of the binary planet sample radius distribution being consistent with the single-star small planet population via an Anderson-Darling test. These results reveal demographic differences in the planet size distribution between planets in binary and single stars for the first time, showing that stellar multiplicity may fundamentally alter the planet formation process. A larger sample and further assessment of circumprimary versus circumsecondary transits is needed to either validate this non-detection or explore other scenarios, such as a radius gap with a location that is dependent on binary separation. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 18 pages, 6 figures, 4 tables, with full-length tables available in the supplemental files. Accepted to AJ

Evidence that Core-Powered Mass-Loss Dominates Over Photoevaporation in Shaping the Kepler Radius Valley

Authors: Travis A. Berger, Joshua E. Schlieder, Daniel Huber, Thomas Barclay

Abstract: The dearth of planets with sizes around 1.8 $\mathrm{R_\oplus}$ is a key demographic feature discovered by the $Kepler$ mission. Two theories have emerged as potential explanations for this valley: photoevaporation and core-powered mass-loss. However, Rogers et al. (2021) shows that differentiating between the two theories is possible using the three-dimensional parameter space of planet radius, i… ▽ More The dearth of planets with sizes around 1.8 $\mathrm{R_\oplus}$ is a key demographic feature discovered by the $Kepler$ mission. Two theories have emerged as potential explanations for this valley: photoevaporation and core-powered mass-loss. However, Rogers et al. (2021) shows that differentiating between the two theories is possible using the three-dimensional parameter space of planet radius, incident flux, and stellar mass. We use homogeneously-derived stellar and planetary parameters to measure the $Kepler$ exoplanet radius gap in this three-dimensional space. We compute the slope of the gap as a function of incident flux at constant stellar mass ($α$ $\equiv$ $\left(\partial \log R_{\mathrm{gap}} / \partial \log S \right)_{M_\star}$) and the slope of the gap as a function of stellar mass at constant incident flux ($β$ $\equiv$ $\left(\partial \log R_{\mathrm{gap}} / \partial \log M_\star \right)_{S}$) and find $α$ = 0.069$^{+0.019}_{-0.023}$ and $β$ = $-$0.046$^{+0.125}_{-0.117}$. Given that Rogers et al. (2021) shows that core-powered mass-loss predicts $α$ $\approx$ 0.08 and $β$ $\approx$ 0.00 while photoevaporation predicts $α$ $\approx$ 0.12 and $β$ $\approx$ --0.17, our measurements are more consistent with core-powered mass-loss than photoevaporation. However, we caution that different gap-determination methods can produce systematic offsets in both $α$ and $β$; therefore, we motivate a comprehensive re-analysis of $Kepler$ light curves with modern, updated priors on eccentricity and mean stellar density to improve both the accuracy and precision of planet radii and subsequent measurements of the gap. △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: 10 pages, 4 figures. Submitted to ApJL

The Gaia-Kepler-TESS-Host Stellar Properties Catalog: Uniform Physical Parameters for 7993 Host Stars and 9324 Planets

Authors: Travis A. Berger, Joshua E. Schlieder, Daniel Huber

Abstract: We present the first homogeneous catalog of $Kepler$, $K2$, and $TESS$ host stars and the corresponding catalog of exoplanet properties, which contain 7993 stars and 9324 planets, respectively. We used isochrone fitting and $Gaia$ DR3 photometry, parallaxes, and spectrophotometric metallicities to compute precise, homogeneous $T_{\mathrm{eff}}$, $\log g$, masses, radii, mean stellar densities, lum… ▽ More We present the first homogeneous catalog of $Kepler$, $K2$, and $TESS$ host stars and the corresponding catalog of exoplanet properties, which contain 7993 stars and 9324 planets, respectively. We used isochrone fitting and $Gaia$ DR3 photometry, parallaxes, and spectrophotometric metallicities to compute precise, homogeneous $T_{\mathrm{eff}}$, $\log g$, masses, radii, mean stellar densities, luminosities, ages, distances, and V-band extinctions for 3248, 565, and 4180 $Kepler$, $K2$, and $TESS$ stars, respectively. We compared our stellar properties to studies using fundamental and precise constraints, such as interferometry and asteroseismology, and find residual scatters of 2.8%, 5.6%, 5.0%, and 31%, with offsets of 0.2%, 1.0%, 1.2%, and 0.7% between our $T_{\mathrm{eff}}$, radii, masses, and ages and those in the literature, respectively. In addition, we compute planet radii, semimajor axes, and incident fluxes for 4281, 676, and 4367 $Kepler$, $K2$, and $TESS$ planets, respectively, and find that the exoplanet radius gap is less prominent in the $K2$, $TESS$, and combined samples than it is in the $Kepler$ sample alone. We suspect this difference is largely due to heterogeneous planet-to-star radius ratios, shorter time baselines of $K2$ and $TESS$, and smaller sample sizes. Finally, we identify a clear radius inflation trend in our large sample of hot Jupiters and find 150 hot sub-Neptunian desert planets, in addition to a population of over 400 young host stars as potential opportunities for testing theories of planet formation and evolution. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 18 pages, 13 figures, and 5 tables. Tables 2, 4, and 5 are included as ancillary files in the sidebar. Submitted to AJ

Helio2024 Science White Paper: Solar and Heliospheric Magnetism in 5D

Authors: Alexei A. Pevtsov, T. Woods, V. Martinez-Pillet, D. Hassler, T. Berger, S. Gosain, T. Hoeksema, A. R. Jones, R. Kohnert, T. Y. Chen, L. Upton, A. Pulkkinen

Abstract: This White Paper argues for the urgent need for the multi-vantage/multi-point observations of the Sun and the heliosphere in the framework of six (6) key science objectives. We further emphasize the critical importance of 5D-``space'': three spatial, one temporal and the magnetic field components. The importance of such observations cannot be overstated both for scientific research and the operati… ▽ More This White Paper argues for the urgent need for the multi-vantage/multi-point observations of the Sun and the heliosphere in the framework of six (6) key science objectives. We further emphasize the critical importance of 5D-``space'': three spatial, one temporal and the magnetic field components. The importance of such observations cannot be overstated both for scientific research and the operational space weather forecast. △ Less

Submitted 12 November, 2022; originally announced November 2022.

Comments: White paper submitted to Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033, 6 pages

Planet Engulfment Detections are Rare According to Observations and Stellar Modeling

Authors: Aida Behmard, Fei Dai, John M. Brewer, Travis A. Berger, Andrew W. Howard

Abstract: Dynamical evolution within planetary systems can cause planets to be engulfed by their host stars. Following engulfment, the stellar photosphere abundance pattern will reflect accretion of rocky material from planets. Multi-star systems are excellent environments to search for such abundance trends because stellar companions form from the same natal gas cloud and are thus expected to share primord… ▽ More Dynamical evolution within planetary systems can cause planets to be engulfed by their host stars. Following engulfment, the stellar photosphere abundance pattern will reflect accretion of rocky material from planets. Multi-star systems are excellent environments to search for such abundance trends because stellar companions form from the same natal gas cloud and are thus expected to share primordial chemical compositions to within 0.03$-$0.05 dex. Abundance measurements have occasionally yielded rocky enhancements, but few observations targeted known planetary systems. To address this gap, we carried out a Keck-HIRES survey of 36 multi-star systems where at least one star is a known planet host. We found that only HAT-P-4 exhibits an abundance pattern suggestive of engulfment, but is more likely primordial based on its large projected separation (30,000 $\pm$ 140 AU) that exceeds typical turbulence scales in molecular clouds. To understand the lack of engulfment detections among our systems, we quantified the strength and duration of refractory enrichments in stellar photospheres using MESA stellar models. We found that observable signatures from 10 $M_{\oplus}$ engulfment events last for $\sim$90 Myr in 1 $M_{\odot}$ stars. Signatures are largest and longest lived for 1.1$-$1.2 $M_{\odot}$ stars, but are no longer observable $\sim$2 Gyr post-engulfment. This indicates that engulfment will rarely be detected in systems that are several Gyr old. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 15 pages, 12 figures; submitted to MNRAS

Seismic Monitoring of the Sun's Far Hemisphere: A Crucial Component in Future Space Weather Forecasting (A White Paper Submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) -- SSPH 2024-2033)

Authors: Kiran Jain, C. Lindsey, E. Adamson, C. N. Arge, T. E. Berger, D. C. Braun, R. Chen, Y. M. Collado-Vega, M. Dikpati, T. Felipe, C. J. Henney, J. T. Hoeksema, R. W. Komm, K. D. Leka, A. R. Marble, V. Martinez Pillet, M. Miesch, L. J. Nickisch, A. A. Pevtsov, V. J. Pizzo, W. K. Tobiska, S. C. Tripathy, J. Zhao

Abstract: The purpose of this white paper is to put together a coherent vision for the role of helioseismic monitoring of magnetic activity in the Sun's far hemisphere that will contribute to improving space weather forecasting as well as fundamental research in the coming decade. Our goal fits into the broader context of helioseismology in solar research for any number of endeavors when helioseismic monito… ▽ More The purpose of this white paper is to put together a coherent vision for the role of helioseismic monitoring of magnetic activity in the Sun's far hemisphere that will contribute to improving space weather forecasting as well as fundamental research in the coming decade. Our goal fits into the broader context of helioseismology in solar research for any number of endeavors when helioseismic monitors may be the sole synoptic view of the Sun's far hemisphere. It is intended to foster a growing understanding of solar activity, as realistically monitored in both hemispheres, and its relationship to all known aspects of the near-Earth and terrestrial environment. Some of the questions and goals that can be fruitfully pursued through seismic monitoring of farside solar activity in the coming decade include: What is the relationship between helioseismic signatures and their associated magnetic configurations, and how is this relationship connected to the solar EUV irradiance over the period of a solar rotation?; How can helioseismic monitoring contribute to data-driven global magnetic-field models for precise space weather forecasting?; What can helioseismic monitors tell us about prospects of a flare, CME or high-speed stream that impacts the terrestrial environment over the period of a solar rotation?; How does the inclusion of farside information contribute to forecasts of interplanetary space weather and the environments to be encountered by human crews in interplanetary space? Thus, it is crucial for the development of farside monitoring of the Sun be continued into the next decade either through ground-based or space-borne observations. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Is [Y/Mg] a Reliable Age Diagnostic for FGK Stars?

Authors: Travis A. Berger, Jennifer L. van Saders, Daniel Huber, Eric Gaidos, Joshua E. Schlieder, Zachary R. Claytor

Abstract: Current spectroscopic surveys are producing large catalogs of chemical abundances for stars of all types. The yttrium to magnesium ratio, [Y/Mg], has emerged as a candidate age indicator for solar twins in the local stellar neighborhood. However, it is unclear whether it is a viable age diagnostic for more diverse stellar types, so we investigate [Y/Mg] as an age indicator for the FGK-type planet… ▽ More Current spectroscopic surveys are producing large catalogs of chemical abundances for stars of all types. The yttrium to magnesium ratio, [Y/Mg], has emerged as a candidate age indicator for solar twins in the local stellar neighborhood. However, it is unclear whether it is a viable age diagnostic for more diverse stellar types, so we investigate [Y/Mg] as an age indicator for the FGK-type planet host stars observed by $Kepler$. We find that the [Y/Mg] "Clock" is most precise for solar twins, with a [Y/Mg]/Age slope of $m$ = $-$0.0370 $\pm$ 0.0071 dex/Gyr and $σ_{\mathrm{Age}}$ = 2.6 Gyr. We attribute the lower precision compared to literature results to non-solar twins contaminating our solar twin sample and recommend a 1.5 Gyr systematic uncertainty for stellar ages derived with any [Y/Mg]-Age relation. We also analyzed the [Y/Mg] Clock as a function of $T_{\mathrm{eff}}$, $\log g$, and metallicity individually and find no strong trends, but compute statistically significant [Y/Mg]-Age relations for subsamples defined by ranges in $T_{\mathrm{eff}}$, $\log g$, and metallicity. Finally, we compare [Y/Mg] and rotation ages and find statistically similar trends as for isochrone ages, although we find that rotation ages perform better for GK dwarfs while isochrones perform better for FG subgiants. We conclude that the [Y/Mg] Clock is most precise for solar twins and analogs but is also a useful age diagnostic for FGK stars. △ Less

Submitted 4 August, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: 18 pages, 13 figures, 2 tables. Accepted by ApJ

The Discovery of a Planetary Companion Interior to Hot Jupiter WASP-132 b

Authors: Benjamin J. Hord, Knicole D. Colón, Travis A. Berger, Veselin Kostov, Michele L. Silverstein, Keivan G. Stassun, Jack J. Lissauer, Karen A. Collins, Richard P. Schwarz, Ramotholo Sefako, Carl Ziegler, César Briceño, Nicholas Law, Andrew W. Mann, George R. Ricker, David W. Latham, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Luke G. Bouma, Ben Falk, Guillermo Torres, Joseph D. Twicken, Andrew Vanderburg

Abstract: Hot Jupiters are generally observed to lack close planetary companions, a trend that has been interpreted as evidence for high-eccentricity migration. We present the discovery and validation of WASP-132 c (TOI-822.02), a 1.85 $\pm$ 0.10 $R_{\oplus}$ planet on a 1.01 day orbit interior to the hot Jupiter WASP-132 b. Transiting Exoplanet Survey Satellite (TESS) and ground-based follow-up observation… ▽ More Hot Jupiters are generally observed to lack close planetary companions, a trend that has been interpreted as evidence for high-eccentricity migration. We present the discovery and validation of WASP-132 c (TOI-822.02), a 1.85 $\pm$ 0.10 $R_{\oplus}$ planet on a 1.01 day orbit interior to the hot Jupiter WASP-132 b. Transiting Exoplanet Survey Satellite (TESS) and ground-based follow-up observations, in conjunction with vetting and validation analysis, enable us to rule out common astrophysical false positives and validate the observed transit signal produced by WASP-132 c as a planet. Running the validation tools \texttt{vespa} and \texttt{triceratops} on this signal yield false positive probabilities of $9.02 \times 10^{-5}$ and 0.0107, respectively. Analysis of archival CORALIE radial velocity data leads to a 3$σ$ upper limit of 28.23 ms$^{-1}$ on the amplitude of any 1.01-day signal, corresponding to a 3$σ$ upper mass limit of 37.35 $M_{\oplus}$. Dynamical simulations reveal that the system is stable within the 3$σ$ uncertainties on planetary and orbital parameters for timescales of $\sim$100 Myr. The existence of a planetary companion near the hot Jupiter WASP-132 b makes the giant planet's formation and evolution via high-eccentricity migration highly unlikely. Being one of just a handful of nearby planetary companions to hot Jupiters, WASP-132 c carries with it significant implications for the formation of the system and hot Jupiters as a population. △ Less

Submitted 11 May, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

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

Rotation Distributions around the Kraft Break with TESS and Kepler: The Influences of Age, Metallicity, and Binarity

Authors: Ellis A. Avallone, Jamie N. Tayar, Jennifer L. van Saders, Travis A. Berger, Zachary R. Claytor, Rachael L. Beaton, Johanna Teske, Diego Godoy-Rivera, Kaike Pan

Abstract: Stellar rotation is a complex function of mass, metallicity, and age and can be altered by binarity. To understand the importance of these parameters in main sequence stars, we have assembled a sample of observations that spans a range of these parameters using a combination of observations from The Transiting Exoplanet Survey Satellite (TESS) and the Kepler Space Telescope. We find that while we… ▽ More Stellar rotation is a complex function of mass, metallicity, and age and can be altered by binarity. To understand the importance of these parameters in main sequence stars, we have assembled a sample of observations that spans a range of these parameters using a combination of observations from The Transiting Exoplanet Survey Satellite (TESS) and the Kepler Space Telescope. We find that while we can measure rotation periods and identify other classes of stellar variability (e.g., pulsations) from TESS lightcurves, instrument systematics prevent the detection of rotation signals longer than the TESS orbital period of 13.7 days. Due to this detection limit, we also utilize rotation periods constrained using rotational velocities measured by the APOGEE spectroscopic survey and radii estimated using the Gaia mission for both TESS and Kepler stars. From these rotation periods, we 1) find we can track rotational evolution along discrete mass tracks as a function of stellar age, 2) find we are unable to recover trends between rotation and metallicity that were observed by previous studies, and 3) note that our sample reveals that wide binary companions do not affect rotation, while close binary companions cause stars to exhibit more rapid rotation than single stars. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 19 pages, 13 figures, Accepted for publication in the Astrophysical Journal

Machine Learning Approaches to Solar-Flare Forecasting: Is Complex Better?

Authors: Varad Deshmukh, Srinivas Baskar, Elizabeth Bradley, Thomas Berger, James D. Meiss

Abstract: Recently, there has been growing interest in the use of machine-learning methods for predicting solar flares. Initial efforts along these lines employed comparatively simple models, correlating features extracted from observations of sunspot active regions with known instances of flaring. Typically, these models have used physics-inspired features that have been carefully chosen by experts in orde… ▽ More Recently, there has been growing interest in the use of machine-learning methods for predicting solar flares. Initial efforts along these lines employed comparatively simple models, correlating features extracted from observations of sunspot active regions with known instances of flaring. Typically, these models have used physics-inspired features that have been carefully chosen by experts in order to capture the salient features of such magnetic field structures. Over time, the sophistication and complexity of the models involved has grown. However, there has been little evolution in the choice of feature sets, nor any systematic study of whether the additional model complexity is truly useful. Our goal is to address these issues. To that end, we compare the relative prediction performance of machine-learning-based, flare-forecasting models with varying degrees of complexity. We also revisit the feature set design, using topological data analysis to extract shape-based features from magnetic field images of the active regions. Using hyperparameter training for fair comparison of different machine-learning models across different feature sets, we show that simpler models with fewer free parameters \textit{generally perform better than more-complicated models}, ie., powerful machinery does not necessarily guarantee better prediction performance. Secondly, we find that \textit{abstract, shape-based features contain just as much useful information}, for the purposes of flare prediction, as the set of hand-crafted features developed by the solar-physics community over the years. Finally, we study the effects of dimensionality reduction, using principal component analysis, to show that streamlined feature sets, overall, perform just as well as the corresponding full-dimensional versions. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Journal ref: A&A 674, A159 (2023)

Decreasing False Alarm Rates in ML-based Solar Flare Prediction using SDO/HMI Data

Authors: Varad Deshmukh, Natasha Flyer, Kiera Van Der Sande, Thomas Berger

Abstract: A hybrid two-stage machine learning architecture that addresses the problem of excessive false positives (false alarms) in solar flare prediction systems is investigated. The first stage is a convolutional neural network (CNN) model based on the VGG-16 architecture that extracts features from a temporal stack of consecutive Solar Dynamics Observatory (SDO) Helioseismic and Magnetic Imager (HMI) ma… ▽ More A hybrid two-stage machine learning architecture that addresses the problem of excessive false positives (false alarms) in solar flare prediction systems is investigated. The first stage is a convolutional neural network (CNN) model based on the VGG-16 architecture that extracts features from a temporal stack of consecutive Solar Dynamics Observatory (SDO) Helioseismic and Magnetic Imager (HMI) magnetogram images to produce a flaring probability. The probability of flaring is added to a feature vector derived from the magnetograms to train an extremely randomized trees (ERT) model in the second stage to produce a binary deterministic prediction (flare/no flare) in a 12-hour forecast window. To tune the hyperparameters of the architecture a new evaluation metric is introduced, the "scaled True Skill Statistic". It specifically addresses the large discrepancy between the true positive rate and the false positive rate in the highly unbalanced solar flare event training datasets. Through hyperparameter tuning to maximize this new metric, our two-stage architecture drastically reduces false positives by $\approx$ $48\%$ without significantly affecting the true positives (reduction by $\approx$ $12\%$), when compared with predictions from the first stage CNN alone. This, in turn, improves various traditional binary classification metrics sensitive to false positives such as the precision, F1 and the Heidke Skill Score. The end result is a more robust 12-hour flare prediction system that could be combined with current operational flare forecasting methods. Additionally, using the ERT-based feature ranking mechanism, we show that the CNN output probability is highly ranked in terms of flare prediction relevance. △ Less

Submitted 20 November, 2021; originally announced November 2021.

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

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

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

Submitted 16 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS 14 October 2021

Measuring the magnetic origins of solar flares, CMEs and Space Weather

Authors: Philip Judge, Matthias Rempel, Rana Ezzeddine, Lucia Kleint, Ricky Egeland, Svetlana Berdyugina, Thomas Berger, Joan Burkepile, Rebecca Centeno, Giuliana de Toma, Mausumi Dikpati, Yuhong Fan, Holly Gilbert, Daniela Lacatus

Abstract: We take a broad look at the problem of identifying the magnetic solar causes of space weather. With the lackluster performance of extrapolations based upon magnetic field measurements in the photosphere, we identify a region in the near UV part of the spectrum as optimal for studying the development of magnetic free energy over active regions. Using data from SORCE, Hubble Space Telescope, and SKY… ▽ More We take a broad look at the problem of identifying the magnetic solar causes of space weather. With the lackluster performance of extrapolations based upon magnetic field measurements in the photosphere, we identify a region in the near UV part of the spectrum as optimal for studying the development of magnetic free energy over active regions. Using data from SORCE, Hubble Space Telescope, and SKYLAB, along with 1D computations of the near-UV (NUV) spectrum and numerical experiments based on the MURaM radiation-MHD and HanleRT radiative transfer codes, we address multiple challenges. These challenges are best met through a combination of near UV lines of bright \ion{Mg}{2}, and lines of \ion{Fe}{2} and \ion{Fe}{1} (mostly within the $4s-4p$ transition array) which form in the chromosphere up to $2\times10^4$ K. Both Hanle and Zeeman effects can in principle be used to derive vector magnetic fields. However, for any given spectral line the $τ=1$ surfaces are generally geometrically corrugated owing to fine structure such as fibrils and spicules. By using multiple spectral lines spanning different optical depths, magnetic fields across nearly-horizontal surfaces can be inferred in regions of low plasma $β$, from which free energies, magnetic topology and other quantities can be derived. Based upon the recently-reported successful suborbital space measurements of magnetic fields with the CLASP2 instrument, we argue that a modest space-borne telescope will be able to make significant advances in the attempts to predict solar eruptions. Difficulties associated with blended lines are shown to be minor in an Appendix. △ Less

Submitted 16 June, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: Accepted in ApJ, V2 fixed author spelling error

Slow Cooling and Fast Reinflation for Hot Jupiters

Authors: Daniel P. Thorngren, Jonathan J. Fortney, Eric D. Lopez, Travis A. Berger, Daniel Huber

Abstract: The unexpectedly large radii of hot Jupiters are a longstanding mystery whose solution will provide important insights into their interior physics. Many potential solutions have been suggested, which make diverse predictions about the details of inflation. In particular, although any valid model must allow for maintaining large planetary radii, only some allow for radii to increase with time. This… ▽ More The unexpectedly large radii of hot Jupiters are a longstanding mystery whose solution will provide important insights into their interior physics. Many potential solutions have been suggested, which make diverse predictions about the details of inflation. In particular, although any valid model must allow for maintaining large planetary radii, only some allow for radii to increase with time. This reinflation process would potentially occur when the incident flux on the planet is increased. In this work, we examine the observed population of hot Jupiters to see if they grow as their parent stars brighten along the main sequence. We consider the relation between radius and other observables, including mass, incident flux, age, and fractional age (age over main sequence lifetime), and show that main sequence brightening is often sufficient to produce detectable reinflation. We further argue that these provide strong evidence for the relatively rapid reinflation of giant planets, and discuss the implications for proposed heating mechanisms. In our population analysis we also find evidence for a "delayed-cooling effect", wherein planets cool and contract far more slowly than expected. While not capable of explaining the observed radii alone, it may represent an important component of the effect. Finally, we identify a weak negative relationship between stellar metallicity and planet radius which is presumably the result of enhanced planetary bulk metallicity around metal-rich stars and has important implications for planet formation theory. △ Less

Submitted 23 February, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 10 pages, 3 figures, 1 table. Revision improves clarity, but scientific results are unchanged. Accepted to The Astrophysical Journal Letters

The Stars Kepler Missed: Investigating the Kepler Target Selection Function Using Gaia DR2

Authors: Linnea M. Wolniewicz, Travis A. Berger, Daniel Huber

Abstract: The Kepler Mission revolutionized exoplanet science and stellar astrophysics by obtaining highly precise photometry of over 200,000 stars over 4 years. A critical piece of information to exploit Kepler data is its selection function, since all targets had to be selected from a sample of half a million stars on the Kepler CCDs using limited information. Here we use Gaia DR2 to reconstruct the Keple… ▽ More The Kepler Mission revolutionized exoplanet science and stellar astrophysics by obtaining highly precise photometry of over 200,000 stars over 4 years. A critical piece of information to exploit Kepler data is its selection function, since all targets had to be selected from a sample of half a million stars on the Kepler CCDs using limited information. Here we use Gaia DR2 to reconstruct the Kepler selection function and explore possible biases with respect to evolutionary state, stellar multiplicity, and kinematics. We find that the Kepler target selection is nearly complete for stars brighter than $Kp < 14$ mag and was effective at selecting main-sequence stars, with the fraction of observed stars decreasing from 95% to 60% between $14 < Kp < 16$ mag. We find that the observed fraction for subgiant stars is only 10% lower, confirming that a significant number of subgiants selected for observation were believed to be main-sequence stars. Conversely we find a strong selection bias against low-luminosity red giant stars ($R \approx 3-5 R_\odot$, $T_{eff} \approx 5500$K), dropping from 90% at $Kp = 14$ mag to below 30% at $Kp = 16$ mag, confirming that the target selection was efficient at distinguishing dwarfs from giants. We compare the Gaia Re-normalized Unit Weight Error (RUWE) values of the observed and non-observed main sequence stars and find a difference in elevated ($>$ 1.2) RUWE values at $\sim\,5\,σ$ significance, suggesting that the Kepler target selection shows some bias against either close or wide binaries. We furthermore use the Gaia proper motions to show that the Kepler selection function was unbiased with respect to kinematics. △ Less

Submitted 8 January, 2021; originally announced January 2021.

Comments: 15 pages, 11 figures. Submitted to the AAS Astrophysical Journal. Presented at AAS 237

Shape-based Feature Engineering for Solar Flare Prediction

Authors: Varad Deshmukh, Thomas Berger, James Meiss, Elizabeth Bradley

Abstract: Solar flares are caused by magnetic eruptions in active regions (ARs) on the surface of the sun. These events can have significant impacts on human activity, many of which can be mitigated with enough advance warning from good forecasts. To date, machine learning-based flare-prediction methods have employed physics-based attributes of the AR images as features; more recently, there has been some w… ▽ More Solar flares are caused by magnetic eruptions in active regions (ARs) on the surface of the sun. These events can have significant impacts on human activity, many of which can be mitigated with enough advance warning from good forecasts. To date, machine learning-based flare-prediction methods have employed physics-based attributes of the AR images as features; more recently, there has been some work that uses features deduced automatically by deep learning methods (such as convolutional neural networks). We describe a suite of novel shape-based features extracted from magnetogram images of the Sun using the tools of computational topology and computational geometry. We evaluate these features in the context of a multi-layer perceptron (MLP) neural network and compare their performance against the traditional physics-based attributes. We show that these abstract shape-based features outperform the features chosen by the human experts, and that a combination of the two feature sets improves the forecasting capability even further. △ Less

Submitted 28 December, 2020; originally announced December 2020.

Comments: To be published in Proceedings for Innovative Applications of Artificial Intelligence Conference 2021

Journal ref: AAAI Conference on Artificial Intelligence, 35(17), 2021, 15293-15300

TESS Asteroseismology of $α$ Mensae: Benchmark Ages for a G7 Dwarf and its M-dwarf Companion

Authors: Ashley Chontos, Daniel Huber, Travis A. Berger, Hans Kjeldsen, Aldo M. Serenelli, Victor Silva Aguirre, Warrick H. Ball, Sarbani Basu, Timothy R. Bedding, William J. Chaplin, Zachary R. Claytor, Enrico Corsaro, Rafael A. García, Steve B. Howell, Mia S. Lundkvist, Savita Mathur, Travis S. Metcalfe, Martin B. Nielsen, Jia Mian Joel Ong, Zeynep Çelik Orhan, Sibel Örtel, Maïssa Salama, Keivan G. Stassun, R. H. D. Townsend, Jennifer L. van Saders , et al. (5 additional authors not shown)

Abstract: Asteroseismology of bright stars has become increasingly important as a method to determine fundamental properties (in particular ages) of stars. The Kepler Space Telescope initiated a revolution by detecting oscillations in more than 500 main-sequence and subgiant stars. However, most Kepler stars are faint, and therefore have limited constraints from independent methods such as long-baseline int… ▽ More Asteroseismology of bright stars has become increasingly important as a method to determine fundamental properties (in particular ages) of stars. The Kepler Space Telescope initiated a revolution by detecting oscillations in more than 500 main-sequence and subgiant stars. However, most Kepler stars are faint, and therefore have limited constraints from independent methods such as long-baseline interferometry. Here, we present the discovery of solar-like oscillations in $α$ Men A, a naked-eye (V=5.1) G7 dwarf in TESS's Southern Continuous Viewing Zone. Using a combination of astrometry, spectroscopy, and asteroseismology, we precisely characterize the solar analog alpha Men A (Teff = 5569 +/- 62 K, R = 0.960 +/- 0.016 Rsun, M = 0.964 +/- 0.045 Msun). To characterize the fully convective M dwarf companion, we derive empirical relations to estimate mass, radius, and temperature given the absolute Gaia magnitude and metallicity, yielding M = 0.169 +/- 0.006, R = 0.19 +/- 0.01 and Teff = 3054 +/- 44 K. Our asteroseismic age of 6.2 +/- 1.4 (stat) +/- 0.6 (sys) Gyr for the primary places $α$ Men B within a small population of M dwarfs with precisely measured ages. We combined multiple ground-based spectroscopy surveys to reveal an activity cycle of 13.1 +/- 1.1 years, a period similar to that observed in the Sun. We used different gyrochronology models with the asteroseismic age to estimate a rotation period of ~30 days for the primary. Alpha Men A is now the closest (d=10pc) solar analog with a precise asteroseismic age from space-based photometry, making it a prime target for next-generation direct imaging missions searching for true Earth analogs. △ Less

Submitted 4 December, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: Accepted to The Astrophysical Journal; 15 pages, 10 figures

The Occurrence of Rocky Habitable Zone Planets Around Solar-Like Stars from Kepler Data

Authors: Steve Bryson, Michelle Kunimoto, Ravi K. Kopparapu, Jeffrey L. Coughlin, William J. Borucki, David Koch, Victor Silva Aguirre, Christopher Allen, Geert Barentsen, Natalie. M. Batalha, Travis Berger, Alan Boss, Lars A. Buchhave, Christopher J. Burke, Douglas A. Caldwell, Jennifer R. Campbell, Joseph Catanzarite, Hema Chandrasekharan, William J. Chaplin, Jessie L. Christiansen, Jorgen Christensen-Dalsgaard, David R. Ciardi, Bruce D. Clarke, William D. Cochran, Jessie L. Dotson , et al. (57 additional authors not shown)

Abstract: We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orb… ▽ More We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orbiting stars with effective temperatures between 4800 K and 6300 K. We find that $η_\oplus$ for the conservative HZ is between $0.37^{+0.48}_{-0.21}$ (errors reflect 68\% credible intervals) and $0.60^{+0.90}_{-0.36}$ planets per star, while the optimistic HZ occurrence is between $0.58^{+0.73}_{-0.33}$ and $0.88^{+1.28}_{-0.51}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates using both a Poisson likelihood Bayesian analysis and Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with $95\%$ confidence that, on average, the nearest HZ planet around G and K dwarfs is about 6 pc away, and there are about 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun. △ Less

Submitted 3 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: To appear in The Astronomical Journal

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Projected WIMP Sensitivity of the XENONnT Dark Matter Experiment

Authors: The XENON collaboration, E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, V. C. Antochi, E. Angelino, J. R. Angevaare, F. Arneodo, D. Barge, L. Baudis, B. Bauermeister, L. Bellagamba, M. L. Benabderrahmane, T. Berger, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon , et al. (115 additional authors not shown)

Abstract: XENONnT is a dark matter direct detection experiment, utilizing 5.9 t of instrumented liquid xenon, located at the INFN Laboratori Nazionali del Gran Sasso. In this work, we predict the experimental background and project the sensitivity of XENONnT to the detection of weakly interacting massive particles (WIMPs). The expected average differential background rate in the energy region of interest, c… ▽ More XENONnT is a dark matter direct detection experiment, utilizing 5.9 t of instrumented liquid xenon, located at the INFN Laboratori Nazionali del Gran Sasso. In this work, we predict the experimental background and project the sensitivity of XENONnT to the detection of weakly interacting massive particles (WIMPs). The expected average differential background rate in the energy region of interest, corresponding to (1, 13) keV and (4, 50) keV for electronic and nuclear recoils, amounts to $12.3 \pm 0.6$ (keV t y)$^{-1}$ and $(2.2\pm 0.5)\times 10^{-3}$ (keV t y)$^{-1}$, respectively, in a 4 t fiducial mass. We compute unified confidence intervals using the profile construction method, in order to ensure proper coverage. With the exposure goal of 20 t$\,$y, the expected sensitivity to spin-independent WIMP-nucleon interactions reaches a cross-section of $1.4\times10^{-48}$ cm$^2$ for a 50 GeV/c$^2$ mass WIMP at 90% confidence level, more than one order of magnitude beyond the current best limit, set by XENON1T. In addition, we show that for a 50 GeV/c$^2$ WIMP with cross-sections above $2.6\times10^{-48}$ cm$^2$ ($5.0\times10^{-48}$ cm$^2$) the median XENONnT discovery significance exceeds 3$σ$ (5$σ$). The expected sensitivity to the spin-dependent WIMP coupling to neutrons (protons) reaches $2.2\times10^{-43}$ cm$^2$ ($6.0\times10^{-42}$ cm$^2$). △ Less

Submitted 17 November, 2020; v1 submitted 17 July, 2020; originally announced July 2020.

Journal ref: JCAP11(2020)031

Excess Electronic Recoil Events in XENON1T

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, V. C. Antochi, E. Angelino, J. R. Angevaare, F. Arneodo, D. Barge, L. Baudis, B. Bauermeister, L. Bellagamba, M. L. Benabderrahmane, T. Berger, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, B. Cimmino , et al. (114 additional authors not shown)

Abstract: We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector. With an exposure of 0.65 t-y and an unprecedentedly low background rate of $76\pm2$ events/(t y keV) between 1 and 30 keV, the data enables sensitive searches for solar axions, an enhanced neutrino magnetic moment, and bosonic dark matter. An excess over known backgrounds is o… ▽ More We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector. With an exposure of 0.65 t-y and an unprecedentedly low background rate of $76\pm2$ events/(t y keV) between 1 and 30 keV, the data enables sensitive searches for solar axions, an enhanced neutrino magnetic moment, and bosonic dark matter. An excess over known backgrounds is observed at low energies and most prominent between 2 and 3 keV. The solar axion model has a 3.4$σ$ significance, and a 3D 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae}<3.8 \times 10^{-12}$, $g_{ae}g_{an}^{eff}<4.8\times 10^{-18}$, and $g_{ae}g_{aγ}<7.7\times10^{-22} GeV^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{aγ}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$σ$ and a confidence interval of $μ_ν \in (1.4,2.9)\times10^{-11}μ_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $β$ decays of tritium at 3.2$σ$ with a trace amount that can neither be confirmed nor excluded with current knowledge of its production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are reduced to 2.0$σ$ and 0.9$σ$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3\pm0.2$) keV (68% C.L.) with a 3.0$σ$ global (4.0$σ$ local) significance. We also consider the possibility that $^{37}$Ar may be present in the detector and yield a 2.82 keV peak. Contrary to tritium, the $^{37}$Ar concentration can be tightly constrained and is found to be negligible. △ Less

Submitted 16 October, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 26 pages, 15 figures. v2 added Ar37 background discussion and best-fit mass of bosonic dark matter, v3 updated Ar37 discussion, tritium estimation, and solar axion energy spectrum. Data in Fig. 2, 4, and 15, including unbinned energy points in Fig. 4, are available in 10.5281/zenodo.4088778

Journal ref: Phys. Rev. D 102, 072004 (2020)

The Gaia-Kepler Stellar Properties Catalog. II. Planet Radius Demographics as a Function of Stellar Mass and Age

Authors: Travis A. Berger, Daniel Huber, Eric Gaidos, Jennifer L. van Saders, Lauren M. Weiss

Abstract: Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using Gaia Data Release 2 by Berger et al. (2020), we re-compute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to b… ▽ More Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using Gaia Data Release 2 by Berger et al. (2020), we re-compute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to be $d \log R_{\mathrm{p}}$/$d \log M_\star = 0.26^{+0.21}_{-0.16}$, consistent with the slope predicted by a planet mass dependence on stellar mass ($0.24-0.35$) and core-powered mass-loss (0.33). We also find first evidence of a stellar age dependence of the planet populations straddling the radius valley. Specifically, we determine that the fraction of super-Earths ($1-1.8 \mathrm{R_\oplus}$) to sub-Neptunes ($1.8-3.5 \mathrm{R_\oplus}$) increases from $0.61 \pm 0.09$ at young ages (< 1 Gyr) to $1.00 \pm 0.10$ at old ages (> 1 Gyr), consistent with the prediction by core-powered mass-loss that the mechanism shaping the radius valley operates over Gyr timescales. Additionally, we find a tentative decrease in the radii of relatively cool ($F_{\mathrm{p}} < 150 \mathrm{F_\oplus}$) sub-Neptunes over Gyr timescales, which suggests that these planets may possess H/He envelopes instead of higher mean molecular weight atmospheres. We confirm the existence of planets within the hot sub-Neptunian "desert" ($2.2 < R_{\mathrm{p}} < 3.8 \mathrm{R_\oplus}$, $F_{\mathrm{p}} > 650 \mathrm{F_\oplus}$) and show that these planets are preferentially orbiting more evolved stars compared to other planets at similar incident fluxes. In addition, we identify candidates for cool ($F_{\mathrm{p}} < 20 \mathrm{F_\oplus}$) inflated Jupiters, present a revised list of habitable zone candidates, and find that the ages of single- and multiple-transiting planet systems are statistically indistinguishable. △ Less

Submitted 29 June, 2020; v1 submitted 29 May, 2020; originally announced May 2020.

Comments: 22 pages, 13 figures, 1 table. Accepted by AJ. The electronic version of Table 1 is available as an ancillary file (sidebar on the right). For a brief video explaining this paper, see https://youtu.be/d42ffDq1D94

Exploring the evolution of stellar rotation using Galactic kinematics

Authors: Ruth Angus, Angus Beane, Adrian M. Price-Whelan, Elisabeth Newton, Jason L. Curtis, Travis Berger, Jennifer van Saders, Rocio Kiman, Daniel Foreman-Mackey, Yuxi Lu, Lauren Anderson, Jacqueline K. Faherty

Abstract: The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the… ▽ More The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: Accepted for publication in the Astronomical Journal

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

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

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

Submitted 28 March, 2020; originally announced March 2020.

Comments: Accepted to MNRAS

Leveraging the Mathematics of Shape for Solar Magnetic Eruption Prediction

Authors: V. Deshmukh, T. E. Berger, E. Bradley, J. D. Meiss

Abstract: Current operational forecasts of solar eruptions are made by human experts using a combination of qualitative shape-based classification systems and historical data about flaring frequencies. In the past decade, there has been a great deal of interest in crafting machine-learning (ML) flare-prediction methods to extract underlying patterns from a training set---e.g., a set of solar magnetogram ima… ▽ More Current operational forecasts of solar eruptions are made by human experts using a combination of qualitative shape-based classification systems and historical data about flaring frequencies. In the past decade, there has been a great deal of interest in crafting machine-learning (ML) flare-prediction methods to extract underlying patterns from a training set---e.g., a set of solar magnetogram images, each characterized by features derived from the magnetic field and labeled as to whether it was an eruption precursor. These patterns, captured by various methods (neural nets, support vector machines, etc.), can then be used to classify new images. A major challenge with any ML method is the \textit{featurization} of the data: pre-processing the raw images to extract higher-level properties, such as characteristics of the magnetic field, that can streamline the training and use of these methods. It is key to choose features that are informative, from the standpoint of the task at hand. To date, the majority of ML-based solar eruption methods have used physics-based magnetic and electric field features such as the total unsigned magnetic flux, the gradients of the fields, the vertical current density, etc. In this paper, we extend the relevant feature set to include characteristics of the magnetic field that are based purely on the geometry and topology of 2D magnetogram images and show that this improves the prediction accuracy of a neural-net based flare-prediction method. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Journal ref: Space Weather and Space Climate 10: 13 (2020)

Energy resolution and linearity of XENON1T in the MeV energy range

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, V. C. Antochi, E. Angelino, J. Angevaare, F. Arneodo, D. Barge, L. Baudis, B. Bauermeister, L. Bellagamba, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon , et al. (113 additional authors not shown)

Abstract: Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XEN… ▽ More Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{ββ}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $σ/μ$ is as low as (0.80$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{ββ}$. We also present a new signal correction method to rectify the saturation effects of the signal readout system, resulting in more accurate position reconstruction and indirectly improving the energy resolution. The very good result achieved in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events. △ Less

Submitted 9 September, 2020; v1 submitted 8 March, 2020; originally announced March 2020.

Comments: 9 pages, 7 figures

Journal ref: Eur. Phys. J. C 80, 785 (2020)

Calibrating GONG Magnetograms with End-to-end Instrument Simulation II: Theory of Calibration

Authors: Joseph Plowman, Thomas Berger

Abstract: This is the second of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. In the first paper, we described the GONG instrument and our `end-to-end' simulation of its measurement process. In this paper, we consider the theory of calibration, and magnetograph comparison in general, identifying some of the significant i… ▽ More This is the second of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. In the first paper, we described the GONG instrument and our `end-to-end' simulation of its measurement process. In this paper, we consider the theory of calibration, and magnetograph comparison in general, identifying some of the significant issues and pitfalls. The calibration of a magnetograph is a function of whether or not it preserves flux, independent of its spatial resolution. However, we find that the one-dimensional comparison methods most often used for magnetograph calibration and comparison will show dramatic differences between two magnetograms with differing spatial resolution, even if they both preserve flux. Some of the apparent disagreement between magnetograms found in the literature are likely a result of these instrumental resolution differences rather than any intrinsic calibration differences. To avoid them, spatial resolution must be carefully matched prior to comparing magnetograms or making calibration curves. In the third paper, we apply the lessons learned here to absolute calibration of GONG using our `end-to-end' measurement simulation. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: Submitted to Solar Physics, second of three papers (the other two are arXiv submissions 3034511 and 3034660)

Calibrating GONG Magnetograms with End-to-end Instrument Simulation I: Background, the GONG Instrument, and End-to-end Simulation

Authors: Joseph Plowman, Thomas Berger

Abstract: This is the first of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. The input to this simulation is a MURaM 3D MHD photospheric simulation and the output is the corresponding set of simulated data numbers which would be recorded by the GONG detectors. These simulated data numbers are then used to produce `synthe… ▽ More This is the first of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. The input to this simulation is a MURaM 3D MHD photospheric simulation and the output is the corresponding set of simulated data numbers which would be recorded by the GONG detectors. These simulated data numbers are then used to produce `synthetic magnetograms' which can be compared with the simulation inputs. This paper describes the GONG instrument, the MURaM datacube, our instrument simulator, and calculation of synthetic magnetograms, setting the stage for the subsequent two papers. These will first lay groundwork for calibration (and magnetogram comparison in general), then apply them to calibration of GONG using the simulation results. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: Submitted to Solar Physics, first of three papers (the other two are arXiv submissions 3034641 and 3034660)

Calibrating GONG Magnetograms with End-to-end Instrument Simulation III: Comparison, Calibration, and Results

Authors: Joseph Plowman, Thomas Berger

Abstract: This is the last of three papers describing an `absolute' calibration of the GONG magnetograph using and end-to-end simulation of its measurement process. The simulation begins with a MURaM 3D MHD datacube and ends with a `synthetic magnetogram' of the corresponding magnetic field values as they would be observed by GONG. We determine a calibration by comparing the synthetic magnetic field measure… ▽ More This is the last of three papers describing an `absolute' calibration of the GONG magnetograph using and end-to-end simulation of its measurement process. The simulation begins with a MURaM 3D MHD datacube and ends with a `synthetic magnetogram' of the corresponding magnetic field values as they would be observed by GONG. We determine a calibration by comparing the synthetic magnetic field measurements with the MURaM magnetic field values that produced them. The previous two papers have described the GONG measurement process (both instrument and data processing), our simulation of it, and the theory of magnetogram comparison and calibration. In this paper, we address some final points on calibration, combine all of this work into a set of calibration curves, and consider the results. We also review the results of the previous two papers for locality of reference. Our calibration indicates that GONG magnetograms underestimate weak flux by a factor of over 2 near disk center, but that factor decreases to about 1 as the line-of-sight approaches the limb. A preliminary investigation of the generalizability of these results suggests other instruments will be affected in a similar way. We also find that some differences in previous magnetograph comparisons are artifacts of instrumental resolution which do not reflect an intrinsic calibration difference, and the measurements are more similar than sometimes thought. These results are directly applicable to question of solar wind prediction model accuracies, particularly in the search for the cause of the common discrepancy between predicted solar wind magnetic flux at 1 AU and values measured in situ by current satellite missions. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: Submitted to Solar Physics, third of three papers (the other two are arXiv submissions 3034511 and 3034641)

The Lunar Lander Neutron and Dosimetry (LND) Experiment on Chang'E 4

Authors: Robert F. Wimmer-Schweingruber, Jia Yu, Stephan I. Böttcher, Shenyi Zhang, Sönke Burmeister, Henning Lohf, Jingnan Guo, Zigong Xu, Björn Schuster, Lars Seimetz, Johan L. Freiherr von Forstner, Ali Ravanbakhsh, Violetta Knierim, Stefan Kolbe, Hauke Woyciechowsky, Shrinivasrao R. Kulkarni, Bin Yuan, Guohong Shen, Chunqing Wang, Zheng Chang, Thomas Berger, Christine E. Hellweg, Daniel Matthiä, Donghui Hou, Alke Knappmann , et al. (4 additional authors not shown)

Abstract: Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von Kármán crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientifi… ▽ More Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von Kármán crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 38 pages, submitted to Space Science Reviews

The Gaia-Kepler Stellar Properties Catalog. I. Homogeneous Fundamental Properties for 186,301 Kepler Stars

Authors: Travis A. Berger, Daniel Huber, Jennifer L. van Saders, Eric Gaidos, Jamie Tayar, Adam L. Kraus

Abstract: An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia-Kepler Stellar Properties Catalog, a set of stellar p… ▽ More An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia-Kepler Stellar Properties Catalog, a set of stellar properties of 186,301 Kepler stars, homogeneously derived from isochrones and broadband photometry, Gaia Data Release 2 parallaxes, and spectroscopic metallicities, where available. Our photometric effective temperatures, derived from $g-K_s$ colors, are calibrated on stars with interferometric angular diameters. Median catalog uncertainties are 112 K for $T_{\mathrm{eff}}$, 0.05 dex for $\log g$, 4% for $R_\star$, 7% for $M_\star$, 13% for $ρ_\star$, 10% for $L_\star$, and 56% for stellar age. These precise constraints on stellar properties for this sample of stars will allow unprecedented investigations into trends in stellar and exoplanet properties as a function of stellar mass and age. In addition, our homogeneous parameter determinations will permit more accurate calculations of planet occurrence and trends with stellar properties. △ Less

Submitted 2 June, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: 22 pages, 13 figures, 2 tables. Published in AJ. Electronic versions of Tables 1 & 2 are now, in this submission, included as ancillary files (sidebar on the right)

The full $\mathit{Kepler}$ phase curve of the eclipsing hot white dwarf binary system KOI-964

Authors: Ian Wong, Avi Shporer, Juliette C. Becker, Benjamin J. Fulton, Travis A. Berger, Nevin N. Weinberg, Phil Arras, Andrew W. Howard, Björn Benneke

Abstract: We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011… ▽ More We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011. A periodogram of the residuals from the phase curve fit reveals an additional stellar variability signal from the host star with a characteristic period of $0.620276\pm0.000011$ days and a full amplitude of $24\pm2$ ppm. We also present new Keck/HIRES radial velocity observations which we use to measure the orbit and obtain a mass ratio of $q=0.106\pm0.012$. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are $2.23\pm0.12\,M_{\odot}$ and $0.236^{+0.028}_{-0.027}\,M_{\odot}$, respectively. The effective temperatures of the two components are $9940^{+260}_{-230}$ K and $15080\pm400$ K, respectively, and we determine the age of the system to be $0.21^{+0.11}_{-0.08}$ Gyr. We use the measured system properties to compute predicted phase curve amplitudes and find that while the measured Doppler boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible non-synchronous rotation of the host star. △ Less

Submitted 22 November, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: Accepted for publication in AJ, 18 pages, 5 figures

Sensitivity Analyses of Exoplanet Occurrence Rates from Kepler and Gaia

Authors: Megan Shabram, Natalie Batalha, Susan E. Thompson, Danley C. Hsu, Eric B. Ford, Jessie L. Christiansen, Daniel Huber, Travis Berger, Joseph Catanzarite, Benjamin E. Nelson, Steve Bryson, Ruslan Belikov, Chris Burke, Doug Caldwell

Abstract: We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the genera… ▽ More We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation and evolution around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog). We attribute the smaller mean to many of the instrumental false positives at longer orbital periods being removed from the $DR25$ catalog. We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space. △ Less

Submitted 26 September, 2020; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: 16 pages, 3 figures, published in the Astrophysical Journal

Journal ref: Astrophys. J., 160, 1 (2020)

Search for Light Dark Matter Interactions Enhanced by the Migdal effect or Bremsstrahlung in XENON1T

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, V. C. Antochi, E. Angelino, F. Arneodo, D. Barge, L. Baudis, B. Bauermeister, L. Bellagamba, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre , et al. (109 additional authors not shown)

Abstract: Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, whic… ▽ More Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach. △ Less

Submitted 18 August, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. Lett. 123, 241803 (2019)

Light Dark Matter Search with Ionization Signals in XENON1T

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, V. C. Antochi, E. Angelino, F. Arneodo, D. Barge, L. Baudis, B. Bauermeister, L. Bellagamba, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre , et al. (108 additional authors not shown)

Abstract: We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe… ▽ More We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_χ$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_χ> 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_χ$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$. △ Less

Submitted 17 December, 2019; v1 submitted 26 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. Lett. 123, 251801 (2019)

Planetesimals Around Stars with TESS (PAST): I. Transient Dimming of a Binary Solar Analog at the End of the Planet Accretion Era

Authors: E. Gaidos, T. Jacobs, D. LaCourse, A. Vanderburg, S. Rappaport, T. Berger, L. Pearce, A. W. Mann, L. Weiss, B. Fulton, A. Behmard, A. W. Howard, M. Ansdell, G. R. Ricker, R. K. Vanderspek, D. W. Latham, S. Seager, J. N. Winn, J. M. Jenkins

Abstract: We report detection of quasi-periodic (1.5 day) dimming of HD 240779, the solar-mass primary in a 5" visual binary (also TIC 284730577), by the Transiting Exoplanet Survey Satellite. This dimming, as has been shown for other "dipper" stars, is likely due to occultation by circumstellar dust. The barycentric space motion, lithium abundance, rotation, and chromospheric emission of the stars in this… ▽ More We report detection of quasi-periodic (1.5 day) dimming of HD 240779, the solar-mass primary in a 5" visual binary (also TIC 284730577), by the Transiting Exoplanet Survey Satellite. This dimming, as has been shown for other "dipper" stars, is likely due to occultation by circumstellar dust. The barycentric space motion, lithium abundance, rotation, and chromospheric emission of the stars in this system point to an age of ~125 Myr, and possible membership in the AB Doradus moving group. As such it occupies an important but poorly explored intermediate regime of stars with transient dimming between young stellar objects in star forming regions and main sequence stars, and between UX Orionis-type Ae/Be stars and M-type "dippers". HD 240779, but not its companion BD+10714B, has WISE-detected excess infrared emission at 12 and 22 microns indicative of circumstellar dust. We propose that infrared emission is produced by collisions of planetesimals during clearing of a residual disk at the end of rocky planet formation, and that quasi-periodic dimming is produced by the rapid disintegration of a 100 km planetesimal near the silicate evaporation radius. Further studies of this and similar systems will illuminate a poorly understood final phase of rocky planet formation like that which produced the inner Solar System. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: to appear in MNRAS

Growth Model Interpretation of Planet Size Distribution

Authors: Li Zeng, Stein B. Jacobsen, Dimitar D. Sasselov, Michail I. Petaev, Andrew Vanderburg, Mercedes Lopez-Morales, Juan Perez-Mercader, Thomas R. Mattsson, Gongjie Li, Matthew Z. Heising, Aldo S. Bonomo, Mario Damasso, Travis A. Berger, Hao Cao, Amit Levi, Robin D. Wordsworth

Abstract: The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observatio… ▽ More The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: PNAS link: https://www.pnas.org/content/116/20/9723 Complete data and mass-radius tables are available at: https://www.cfa.harvard.edu/~lzeng/planetmodels.html

Journal ref: PNAS May 2019, 116 (20) 9723-9728

A Probabilistic Approach to Kepler Completeness and Reliability for Exoplanet Occurrence Rates

Authors: Steve Bryson, Jeffrey Coughlin, Natalie M. Batalha, Travis Berger, Dan Huber, Christopher Burke, Jessie Dotson, Susan E. Mullally

Abstract: Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uni… ▽ More Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uniformly vetted planet candidate catalog and data products that facilitate corrections. We present a new probabilistic approach to the characterization of Kepler completeness and reliability, making full use of the Kepler DR25 products. We illustrate the impact of completeness and reliability corrections with a Poisson-likelihood occurrence rate method, using a recent stellar properties catalog that incorporates Gaia stellar radii and essentially uniform treatment of the stellar population. Correcting for reliability has a significant impact: the exoplanet occurrence rate for orbital period and radius within 20% of Earth's around GK dwarf stars, corrected for reliability, is 0.015+0.011-0.007, whereas not correcting results in 0.034+0.018-0.012 - correcting for reliability reduces this occurrence rate by more than a factor of two. We further show that using Gaia-based vs. DR25 stellar properties impacts the same occurrence rate by a factor of two. We critically examine the the DR25 catalog and the assumptions behind our occurrence rate method. We propose several ways in which confidence in both the Kepler catalog and occurrence rate calculations can be improved. This work provides an example of how the community can use the DR25 completeness and reliability products. △ Less

Submitted 6 February, 2020; v1 submitted 9 June, 2019; originally announced June 2019.

Comments: Accepted for publication in The Astronomical Journal

The Curious Case of KOI 4: Confirming Kepler's First Exoplanet

Authors: Ashley Chontos, Daniel Huber, David W. Latham, Allyson Bieryla, Vincent Van Eylen, Timothy R. Bedding, Travis Berger, Lars A. Buchhave, Tiago L. Campante, William J. Chaplin, Isabel L. Colman, Jeff L. Coughlin, Guy Davies, Teruyuki Hirano, Andrew W. Howard, Howard Isaacson

Abstract: The discovery of thousands of planetary systems by Kepler has demonstrated that planets are ubiquitous. However, a major challenge has been the confirmation of Kepler planet candidates, many of which still await confirmation. One of the most enigmatic examples is KOI 4.01, Kepler's first discovered planet candidate detection (as KOI 1.01, 2.01, and 3.01 were known prior to launch). Here we present… ▽ More The discovery of thousands of planetary systems by Kepler has demonstrated that planets are ubiquitous. However, a major challenge has been the confirmation of Kepler planet candidates, many of which still await confirmation. One of the most enigmatic examples is KOI 4.01, Kepler's first discovered planet candidate detection (as KOI 1.01, 2.01, and 3.01 were known prior to launch). Here we present the confirmation and characterization of KOI 4.01 (now Kepler-1658), using a combination of asteroseismology and radial velocities. Kepler-1658 is a massive, evolved subgiant (Mstar = 1.45 +/- 0.06 Msun, Rstar = 2.89 +/- 0.12 Rsun) hosting a massive (Mp = 5.88 +/- 0.47 MJ, Rp = 1.07 +/- 0.05 RJ) hot Jupiter that orbits every 3.85 days. Kepler-1658 joins a small population of evolved hosts with short-period (<=100 days) planets and is now the closest known planet in terms of orbital period to an evolved star. Because of its uniqueness and short orbital period, Kepler-1658 is a new benchmark system for testing tidal dissipation and hot Jupiter formation theories. Using all 4 years of Kepler data, we constrain the orbital decay rate to be Pdot <= -0.42 s/yr, corresponding to a strong observational limit of Qstar >= 4.826 x 10^3 for the tidal quality factor in evolved stars. With an effective temperature Teff ~6200 K, Kepler-1658 sits close to the spin-orbit misalignment boundary at ~6250 K, making it a prime target for follow-up observations to better constrain its obliquity and to provide insight into theories for hot Jupiter formation and migration. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: 18 pages, 11 figures, 5 tables; accepted for publication in AJ

Constraining the Spin-Dependent WIMP-Nucleon Cross Sections with XENON1T

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, M. Anthony, V. C. Antochi, F. Arneodo, L. Baudis, B. Bauermeister, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre, A. P. Colijn , et al. (105 additional authors not shown)

Abstract: We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set ex… ▽ More We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator. △ Less

Submitted 30 April, 2019; v1 submitted 8 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. Lett. 122, 141301 (2019)

First results on the scalar WIMP-pion coupling, using the XENON1T experiment

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, M. Anthony, V. C. Antochi, F. Arneodo, L. Baudis, B. Bauermeister, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre, A. P. Colijn , et al. (107 additional authors not shown)

Abstract: We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in… ▽ More We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass. △ Less

Submitted 22 February, 2019; v1 submitted 29 November, 2018; originally announced November 2018.

Report number: INT-PUB-18-058

Journal ref: Phys. Rev. Lett. 122, 071301 (2019)

Roadmap for Reliable Ensemble Forecasting of the Sun-Earth System

Authors: Gelu Nita, Rafal Angryk, Berkay Aydin, Juan Banda, Tim Bastian, Tom Berger, Veronica Bindi, Laura Boucheron, Wenda Cao, Eric Christian, Georgia de Nolfo, Edward DeLuca, Marc DeRosa, Cooper Downs, Gregory Fleishman, Olac Fuentes, Dale Gary, Frank Hill, Todd Hoeksema, Qiang Hu, Raluca Ilie, Jack Ireland, Farzad Kamalabadi, Kelly Korreck, Alexander Kosovichev , et al. (22 additional authors not shown)

Abstract: The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the und… ▽ More The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed. △ Less

Submitted 29 October, 2018; v1 submitted 19 October, 2018; originally announced October 2018.

Comments: Workshop Report

Skynet Algorithm for Single-Dish Radio Mapping I: Contaminant-Cleaning, Mapping, and Photometering Small-Scale Structures

Authors: J. R. Martin, D. E. Reichart, D. A. Dutton, M. P. Maples, T. A. Berger, F. D. Ghigo, J. B. Haislip, O. H. Shaban, A. S. Trotter, L. M. Barnes, M. L. Paggen, R. L. Gao, C. P. Salemi, G. I. Langston, S. Bussa, J. A. Duncan, S. White, S. A. Heatherly, J. B. Karlik, E. M. Johnson, J. E. Reichart, A. C. Foster, V. V. Kouprianov, S. Mazlin, J. Harvey

Abstract: We present a single-dish mapping algorithm with a number of advantages over traditional techniques. (1) Our algorithm makes use of weighted modeling, instead of weighted averaging, to interpolate between signal measurements. This smooths the data, but without blurring the data beyond instrumental resolution. Techniques that rely on weighted averaging blur point sources sometimes as much as 40%. (2… ▽ More We present a single-dish mapping algorithm with a number of advantages over traditional techniques. (1) Our algorithm makes use of weighted modeling, instead of weighted averaging, to interpolate between signal measurements. This smooths the data, but without blurring the data beyond instrumental resolution. Techniques that rely on weighted averaging blur point sources sometimes as much as 40%. (2) Our algorithm makes use of local, instead of global, modeling to separate astronomical signal from instrumental and/or environmental signal drift along the telescope's scans. Other techniques, such as basket weaving, model this drift with simple functional forms (linear, quadratic, etc.) across the entirety of scans, limiting their ability to remove such contaminants. (3) Our algorithm makes use of a similar, local modeling technique to separate astronomical signal from radio-frequency interference (RFI), even if only continuum data are available. (4) Unlike other techniques, our algorithm does not require data to be collected on a rectangular grid or regridded before processing. (5) Data from any number of observations, overlapping or not, may be appended and processed together. (6) Any pixel density may be selected for the final image. We present our algorithm, and evaluate it using both simulated and real data. We are integrating it into the image-processing library of the Skynet Robotic Telescope Network, which includes optical telescopes spanning four continents, and now also Green Bank Observatory's 20-meter diameter radio telescope in West Virginia. Skynet serves hundreds of professional users, and additionally tens of thousands of students, of all ages. Default data products are generated on the fly, but will soon be customizable after the fact. △ Less

Submitted 18 August, 2018; originally announced August 2018.

Comments: 50 pages, 54 figures, 4 tables, accepted for publication in ApJS

Robust Chauvenet Outlier Rejection

Authors: M. P. Maples, D. E. Reichart, N. C. Konz, T. A. Berger, A. S. Trotter, J. R. Martin, D. A. Dutton, M. L. Paggen, R. E. Joyner, C. P. Salemi

Abstract: Sigma clipping is commonly used in astronomy for outlier rejection, but the number of standard deviations beyond which one should clip data from a sample ultimately depends on the size of the sample. Chauvenet rejection is one of the oldest, and simplest, ways to account for this, but, like sigma clipping, depends on the sample's mean and standard deviation, neither of which are robust quantities:… ▽ More Sigma clipping is commonly used in astronomy for outlier rejection, but the number of standard deviations beyond which one should clip data from a sample ultimately depends on the size of the sample. Chauvenet rejection is one of the oldest, and simplest, ways to account for this, but, like sigma clipping, depends on the sample's mean and standard deviation, neither of which are robust quantities: Both are easily contaminated by the very outliers they are being used to reject. Many, more robust measures of central tendency, and of sample deviation, exist, but each has a tradeoff with precision. Here, we demonstrate that outlier rejection can be both very robust and very precise if decreasingly robust but increasingly precise techniques are applied in sequence. To this end, we present a variation on Chauvenet rejection that we call "robust" Chauvenet rejection (RCR), which uses three decreasingly robust/increasingly precise measures of central tendency, and four decreasingly robust/increasingly precise measures of sample deviation. We show this sequential approach to be very effective for a wide variety of contaminant types, even when a significant -- even dominant -- fraction of the sample is contaminated, and especially when the contaminants are strong. Furthermore, we have developed a bulk-rejection variant, to significantly decrease computing times, and RCR can be applied both to weighted data, and when fitting parameterized models to data. We present aperture photometry in a contaminated, crowded field as an example. RCR may be used by anyone at https://skynet.unc.edu/rcr, and source code is available there as well. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 62 pages, 48 figures, 7 tables, accepted for publication in ApJS

Dark Matter Search Results from a One Tonne$\times$Year Exposure of XENON1T

Authors: E. Aprile, J. Aalbers, F. Agostini, M. Alfonsi, L. Althueser, F. D. Amaro, M. Anthony, F. Arneodo, L. Baudis, B. Bauermeister, M. L. Benabderrahmane, T. Berger, P. A. Breur, A. Brown, A. Brown, E. Brown, S. Bruenner, G. Bruno, R. Budnik, C. Capelli, J. M. R. Cardoso, D. Cichon, D. Coderre, A. P. Colijn, J. Conrad , et al. (95 additional authors not shown)

Abstract: We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra… ▽ More We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. △ Less

Submitted 13 September, 2018; v1 submitted 31 May, 2018; originally announced May 2018.

Comments: 8 pages, 5 figures, v2 limit points as csv

Journal ref: Phys. Rev. Lett. 121, 111302 (2018)