Showing 1–50 of 291 results for author: Stone, J

The three-dimensional structure of black hole accretion flows within the plunging region

Authors: Andrew Mummery, James M. Stone

Abstract: We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e., at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry… ▽ More We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e., at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization "SANE" configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the "thin disc" limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 21 pages, 16 figures. Accepted for publication in MNRAS

Exploring the directly imaged HD 1160 system through spectroscopic characterization and high-cadence variability monitoring

Authors: Ben J. Sutlieff, Jayne L. Birkby, Jordan M. Stone, Annelotte Derkink, Frank Backs, David S. Doelman, Matthew A. Kenworthy, Alexander J. Bohn, Steve Ertel, Frans Snik, Charles E. Woodward, Ilya Ilyin, Andrew J. Skemer, Jarron M. Leisenring, Klaus G. Strassmeier, Ji Wang, David Charbonneau, Beth A. Biller

Abstract: The time variability and spectra of directly imaged companions provide insight into their physical properties and atmospheric dynamics. We present follow-up R~40 spectrophotometric monitoring of red companion HD 1160 B at 2.8-4.2 $μ$m using the double-grating 360° vector Apodizing Phase Plate (dgvAPP360) coronagraph and ALES integral field spectrograph on the Large Binocular Telescope Interferomet… ▽ More The time variability and spectra of directly imaged companions provide insight into their physical properties and atmospheric dynamics. We present follow-up R~40 spectrophotometric monitoring of red companion HD 1160 B at 2.8-4.2 $μ$m using the double-grating 360° vector Apodizing Phase Plate (dgvAPP360) coronagraph and ALES integral field spectrograph on the Large Binocular Telescope Interferometer. We use the recently developed technique of gvAPP-enabled differential spectrophotometry to produce differential light curves for HD 1160 B. We reproduce the previously reported ~3.2 h periodic variability in archival data, but detect no periodic variability in new observations taken the following night with a similar 3.5% level precision, suggesting rapid evolution in the variability of HD 1160 B. We also extract complementary spectra of HD 1160 B for each night. The two are mostly consistent, but the companion appears fainter on the second night between 3.0-3.2 $μ$m. Fitting models to these spectra produces different values for physical properties depending on the night considered. We find an effective temperature T$_{\text{eff}}$ = 2794$^{+115}_{-133}$ K on the first night, consistent with the literature, but a cooler T$_{\text{eff}}$ = 2279$^{+79}_{-157}$ K on the next. We estimate the mass of HD 1160 B to be 16-81 M$_{\text{Jup}}$, depending on its age. We also present R = 50,000 high-resolution optical spectroscopy of host star HD 1160 A obtained simultaneously with the PEPSI spectrograph. We reclassify its spectral type to A1 IV-V and measure its projected rotational velocity v sin i = 96$^{+6}_{-4}$ km s$^{-1}$. We thus highlight that gvAPP-enabled differential spectrophotometry can achieve repeatable few percent level precision and does not yet reach a systematic noise floor, suggesting greater precision is achievable with additional data or advanced detrending techniques. △ Less

Submitted 5 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 22 pages, 12 figures, published in MNRAS; typos corrected, references updated

Magnetized Accretion onto and Feedback from Supermassive Black Holes in Elliptical Galaxies

Authors: Minghao Guo, James M. Stone, Eliot Quataert, Chang-Goo Kim

Abstract: We present three-dimensional magnetohydrodynamic (MHD) simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent cooling medium on galactic scales, taking M87* as a typical case. We find that the mass accretion rate is increased by a factor of $\sim 10$ compared with analogous hydrodynamic simulations. The scaling of $\dot{M} \sim r^{1/2}$ roughly holds from… ▽ More We present three-dimensional magnetohydrodynamic (MHD) simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent cooling medium on galactic scales, taking M87* as a typical case. We find that the mass accretion rate is increased by a factor of $\sim 10$ compared with analogous hydrodynamic simulations. The scaling of $\dot{M} \sim r^{1/2}$ roughly holds from $\sim 10\,\mathrm{pc}$ to $\sim 10^{-3}\,\mathrm{pc}$ ($\sim 10\, r_\mathrm{g}$) with the accretion rate through the event horizon being $\sim 10^{-2}\, M_\odot\,\mathrm{yr^{-1}}$. The accretion flow on scales $\sim 0.03-3\,\mathrm{kpc}$ takes the form of magnetized filaments. Within $\sim 30\,\mathrm{pc}$, the cold gas circularizes, forming a highly magnetized ($β\sim 10^{-3}$) thick disk supported by a primarily toroidal magnetic field. The cold disk is truncated and transitions to a turbulent hot accretion flow at $\sim0.3\,\mathrm{pc}$ ($10^3\,r_\mathrm{g}$). There are strong outflows towards the poles driven by the magnetic field. The outflow energy flux increases with smaller accretor size, reaching $\sim 3\times10^{43}\,\mathrm{erg\,s^{-1}}$ for $r_\mathrm{in}=8\,r_\mathrm{g}$; this corresponds to a nearly constant energy feedback efficiency of $η\sim0.05-0.1$ independent of accretor size. The feedback energy is enough to balance the total cooling of the M87/Virgo hot halo out to $\sim 50$ kpc. The accreted magnetic flux at small radii is similar to that in magnetically arrested disk models, consistent with the formation of a powerful jet on horizon scales in M87. Our results motivate a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by $\sim (10r_\mathrm{g}/r_\mathrm{B})^{1/2}$. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 23 pages, 14 figures, submitted to ApJ

Nonthermal Signatures of Radiative Supernova Remnants

Authors: Rebecca Diesing, Minghao Guo, Chang-Goo Kim, James Stone, Damiano Caprioli

Abstract: The end of supernova remnant (SNR) evolution is characterized by a so-called "radiative" stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at t… ▽ More The end of supernova remnant (SNR) evolution is characterized by a so-called "radiative" stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to $γ$-rays, including radio and $γ$-ray brightening by nearly two orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array). △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 12 pages, 6 figures, submitted to ApJ

Dissipation of AGN jets in a clumpy interstellar medium

Authors: Riju Dutta, Prateek Sharma, Kartick C. Sarkar, James M. Stone

Abstract: Accreting supermassive black holes (SMBHs) frequently power jets that interact with the interstellar/circumgalactic medium (ISM/CGM), regulating star-formation in the galaxy. Highly supersonic jets launched by active galactic nuclei (AGN) power a cocoon that confines them and shocks the ambient medium. We build upon the models of narrow conical jets interacting with a smooth ambient medium, to inc… ▽ More Accreting supermassive black holes (SMBHs) frequently power jets that interact with the interstellar/circumgalactic medium (ISM/CGM), regulating star-formation in the galaxy. Highly supersonic jets launched by active galactic nuclei (AGN) power a cocoon that confines them and shocks the ambient medium. We build upon the models of narrow conical jets interacting with a smooth ambient medium, to include the effect of dense clouds that are an essential ingredient of a multiphase ISM. The key physical ingredient of this model is that the clouds along the supersonic jet-beam strongly decelerate the jet-head, but the subsonic cocoon easily moves around the clouds without much resistance. We propose scalings for important physical quantities -- cocoon pressure, head & cocoon speed, and jet radius. We obtain, for the first time, the analytic condition on clumpiness of the ambient medium for the jet to dissipate within the cocoon and verify it with numerical simulations of conical jets interacting with a uniform ISM with embedded spherical clouds. A jet is defined to be dissipated when the cocoon speed exceeds the speed of the jet-head. We compare our models to more sophisticated numerical simulations, direct observations of jet-ISM interaction (e.g., quasar J1316+1753), and discuss implications for the Fermi/eROSITA bubbles. Our work also motivates effective subgrid models for AGN jet feedback in a clumpy ISM unresolved by the present generation of cosmological galaxy formation simulations. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Comments: 23 pages, 12 figures, 3 tables; to be submitted; comments are welcome; accompanying video: http://youtu.be/DUpSwMMrGfk

Tearing-mediated reconnection in magnetohydrodynamic poorly ionized plasmas. I. Onset and linear evolution

Authors: Elizabeth A. Tolman, Matthew W. Kunz, James M. Stone, Lev Arzamasskiy

Abstract: In high-Lundquist-number plasmas, reconnection proceeds via onset of tearing, followed by a nonlinear phase during which plasmoids continuously form, merge, and are ejected from the current sheet (CS). This process is understood in fully ionized, magnetohydrodynamic plasmas. However, many plasma environments, such as star-forming molecular clouds and the solar chromosphere, are poorly ionized. We… ▽ More In high-Lundquist-number plasmas, reconnection proceeds via onset of tearing, followed by a nonlinear phase during which plasmoids continuously form, merge, and are ejected from the current sheet (CS). This process is understood in fully ionized, magnetohydrodynamic plasmas. However, many plasma environments, such as star-forming molecular clouds and the solar chromosphere, are poorly ionized. We use theory and computation to study tearing-mediated reconnection in such poorly ionized systems. In this paper, we focus on the onset and linear evolution of this process. In poorly ionized plasmas, magnetic nulls on scales below $v_{\rm A,n0}/ν_{\rm ni0}$, with $v_{\rm A,n0}$ the neutral Alfvén speed and $ν_{\rm ni0}$ the neutral-ion collision frequency, will self-sharpen via ambipolar diffusion. This sharpening occurs at an increasing rate, inhibiting the onset of reconnection. Once the CS becomes thin enough, however, ions decouple from neutrals and thinning of the CS slows, allowing tearing to onset in a time of order $ν_{\rm ni0}^{-1}$. We find that the wavelength and growth rate of the mode that first disrupts the forming sheet can be predicted from a poorly ionized tearing dispersion relation; as the plasma recombination rate increases and ionization fraction decreases, the growth rate becomes an increasing multiple of $ν_{ni0}$ and the wavelength becomes a decreasing fraction of $v_{\rm A,n0}/ν_{\rm ni0}$. △ Less

Submitted 18 April, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 19 pages, 13 figures

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems V: Do Self-Consistent Atmospheric Models Represent JWST Spectra? A Showcase With VHS 1256 b

Authors: Simon Petrus, Niall Whiteford, Polychronis Patapis, Beth A. Biller, Andrew Skemer, Sasha Hinkley, Genaro Suárez, Anna Lueber, Paulina Palma-Bifani, Jordan M. Stone, Johanna M. Vos, Caroline V. Morley, Pascal Tremblin, Benjamin Charnay, Christiane Helling, Brittany E. Miles, Aarynn L. Carter, Jason J. Wang, Markus Janson, Eileen C. Gonzales, Ben Sutlieff, Kielan K. W. Hoch, Mickaël Bonnefoy, Gaël Chauvin, Olivier Absil , et al. (97 additional authors not shown)

Abstract: The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. W… ▽ More The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models. △ Less

Submitted 31 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 32 pages, 16 figures, 6 tables, 2 appendices

Hartree-Fock Formulation of the QMC Model at Finite Temperature

Authors: P A M Guichon, J R Stone, A W Thomas

Abstract: We present, for the first time, a detailed theory of high density matter including the entire baryon octet at finite temperature, based on a fully relativistic mean field model with a consistent treatment of exchange (Fock) terms, using the quark-meson-coupling model (QMC). It has been already demonstrated that the QMC equation of state is applicable in thermodynamic scenarios in stationary and ro… ▽ More We present, for the first time, a detailed theory of high density matter including the entire baryon octet at finite temperature, based on a fully relativistic mean field model with a consistent treatment of exchange (Fock) terms, using the quark-meson-coupling model (QMC). It has been already demonstrated that the QMC equation of state is applicable in thermodynamic scenarios in stationary and rotating isentropic proto-neutron stars, producing results in agreement with recent observation. It is also suitable for the simulation of the behavior following a binary neutron star merger [1]; https://compose.obspm.fr/eos/205. We develop a comprehensive demonstration of the impact of the Fock terms in the QMC energy density functional on properties of neutrinoless proto-neutron stars with cores containing the full hyperon octet with constant entropy, S/A=2kB. Given the interest in the properties of the proto-neutron star remaining after either a supernova explosion or the merger of two neutron stars, it is vital to develop modern equations of state at finite temperature. While much attention has been paid to relativistic mean-field calculations at finite temperature, it is crucial to explore the consequences of a consistent treatment of the Fock terms. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Report number: ADP-23-25/T1234

Recovering simulated planet and disk signals using SCALES aperture masking

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

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

Submitted 19 October, 2023; originally announced October 2023.

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

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

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

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

Submitted 18 October, 2023; originally announced October 2023.

The \textit{JWST} Early Release Science Program for Direct Observations of Exoplanetary Systems III: Aperture Masking Interferometric Observations of the star HIP\,65426 at $\boldsymbol{3.8\,\rm{μm}}$

Authors: Shrishmoy Ray, Steph Sallum, Sasha Hinkley, Anand Sivamarakrishnan, Rachel Cooper, Jens Kammerer, Alexandra Z. Greebaum, Deepashri Thatte, Cecilia Lazzoni, Andrei Tokovinin, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan , et al. (98 additional authors not shown)

Abstract: We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{JWST} Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of ${}0.5λ/D$ f… ▽ More We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{JWST} Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of ${}0.5λ/D$ for an interferometer), which are inaccessible with the classical inner working angles of the \textit{JWST} coronagraphs. When combined with \textit{JWST}'s unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a contrast of $Δm_{F380M}{\sim }7.8$\,mag relative to the host star at a separation of ${\sim}0.07\arcsec$ but detect no additional companions interior to the known companion HIP\,65426\,b. Our observations thus rule out companions more massive than $10{-}12\,\rm{M\textsubscript{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP\,65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on \textit{JWST} is sensitive to planetary mass companions orbiting at the water frost line, even for more distant stars at $\sim$100\,pc. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening essentially unexplored parameter space. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 15 pages, 9 figures, submitted to ApJ Letters

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems IV: NIRISS Aperture Masking Interferometry Performance and Lessons Learned

Authors: Steph Sallum, Shrishmoy Ray, Jens Kammerer, Anand Sivaramakrishnan, Rachel Cooper, Alexandra Z. Greebaum, Deepashri Thatte, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Sasha Hinkley, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan, Sarah K. Betti, Anthony Boccaletti , et al. (98 additional authors not shown)

Abstract: We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early… ▽ More We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim λ/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy. △ Less

Submitted 11 March, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

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

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

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

Submitted 10 October, 2023; originally announced October 2023.

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

A Systematic Study of Planetary Envelope Growth with 3D Radiation-Hydrodynamics Simulations

Authors: Avery Bailey, Jim Stone, Jeffrey Fung

Abstract: In the core accretion model of planet formation, envelope cooling regulates the accretion of material and ultimately sets the timescale to form a giant planet. Given the diversity of planet-forming environments, opacity uncertainties, and the advective transport of energy by 3-dimensional recycling flows, it is unclear whether 1D models can adequately describe envelope structure and accretion in a… ▽ More In the core accretion model of planet formation, envelope cooling regulates the accretion of material and ultimately sets the timescale to form a giant planet. Given the diversity of planet-forming environments, opacity uncertainties, and the advective transport of energy by 3-dimensional recycling flows, it is unclear whether 1D models can adequately describe envelope structure and accretion in all regimes. Even in 3D models, it is unclear whether approximate radiative transfer methods sufficiently model envelope cooling particularly at the planetary photosphere. To address these uncertainties, we present a suite of 3D radiation hydrodynamics simulations employing methods that directly solve the transfer equation. We perform a parameter space study, formulated in terms of dimensionless parameters, for a variety of envelope optical depths and cooling times. We find that the thermodynamic structure of the envelope ranges from adiabatic to isothermal based on the cooling time and by extension, the background disk temperature and density. Our models show general agreement with 1D static calculations, suggesting a limited role of recycling flows in determining envelope structure. By adopting a dimensionless framework, these models can be applied to a wide range of formation conditions and assumed opacities. In particular, we dimensionalize them to the case of a super-Earth and proto-Jupiter and place upper limits on the 3D mass accretion rates prior to runaway growth. Finally, we evaluate the fidelity of approximate radiative transfer methods and find that even in the most challenging cases, more approximate methods are sufficiently accurate and worth their savings in computational cost. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 15 pages, 16 figures

A Global 3-D Simulation of Magnetospheric Accretion: I. Magnetically Disrupted Disks and Surface Accretion

Authors: Zhaohuan Zhu, James M. Stone, Nuria Calvet

Abstract: We present a 3-D ideal MHD simulation of magnetospheric accretion onto a non-rotating star. The accretion process unfolds with intricate 3-D structures driven by various mechanisms. First, the disc develops filaments at the magnetospheric truncation radius ($R_T$) due to magnetic interchange instability. These filaments penetrate deep into the magnetosphere, form multiple accretion columns, and ev… ▽ More We present a 3-D ideal MHD simulation of magnetospheric accretion onto a non-rotating star. The accretion process unfolds with intricate 3-D structures driven by various mechanisms. First, the disc develops filaments at the magnetospheric truncation radius ($R_T$) due to magnetic interchange instability. These filaments penetrate deep into the magnetosphere, form multiple accretion columns, and eventually impact the star at $\sim$30$^o$ from the poles at nearly the free-fall speed. Over 50% (90%) of accretion occurs on just 5% (20%) of the stellar surface. Second, the disc region outside $R_T$ develops large-scale magnetically dominated bubbles, again due to magnetic interchange instability. These bubbles orbit at a sub-Keplerian speed, persisting for a few orbits while leading to asymmetric mass ejection. The disc outflow is overall weak because of mostly closed field lines. Third, magnetically-supported surface accretion regions appear above the disc, resembling a magnetized disc threaded by net vertical fields, a departure from traditional magnetospheric accretion models. Stellar fields are efficiently transported into the disc region due to above instabilities, contrasting with the ``X-wind'' model. The accretion rate onto the star remains relatively steady with a 23% standard deviation. The periodogram reveals variability occurring at around 0.2 times the Keplerian frequency at $R_T$, linked to the large-scale magnetic bubbles. The ratio of the spin-up torque to $\dot{M}(GM_*R_T)^{1/2}$ is around 0.8. Finally, after scaling the simulation, we investigate planet migration in the inner protoplanetary disc. The disc driven migration is slow in the MHD turbulent disc beyond $R_T$, while aerodynamic drag plays a significant role in migration within $R_T$. △ Less

Submitted 14 December, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: Accepted to MNRAS

Direct images and spectroscopy of a giant protoplanet driving spiral arms in MWC 758

Authors: Kevin Wagner, Jordan Stone, Andrew Skemer, Steve Ertel, Ruobing Dong, Dániel Apai, Eckhart Spalding, Jarron Leisenring, Michael Sitko, Kaitlin Kratter, Travis Barman, Mark Marley, Brittany Miles, Anthony Boccaletti, Korash Assani, Ammar Bayyari, Taichi Uyama, Charles E. Woodward, Phil Hinz, Zackery Briesemeister, Kellen Lawson, François Ménard, Eric Pantin, Ray W. Russell, Michael Skrutskie , et al. (1 additional authors not shown)

Abstract: Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets t… ▽ More Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: Published in Nature Astronomy

Damping of MHD Turbulence in A Partially Ionized Medium

Authors: Yue Hu, Siyao Xu, Lev Arzamasskiy, James M. Stone, A. Lazarian

Abstract: The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions+neutrals) simulations generated with the AthenaK code. Specifically, we examine the velocity, densi… ▽ More The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions+neutrals) simulations generated with the AthenaK code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length scales. △ Less

Submitted 10 November, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 17 pages, 14 figures, accepted for publication in MNRAS

Implementation of chemistry in the Athena++ code

Authors: Munan Gong, Ka-Wai Ho, James M. Stone, Eve C. Ostriker, Paola Caselli, Tommaso Grassi, Chang-Goo Kim, Jeong-Gyu Kim, Goni Halevi

Abstract: Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating… ▽ More Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing the user to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ODEs with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photo-dissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++, aiming to provide a robust and flexible code for the astrochemical simulation community. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Mahakala: a Python-based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Space-times

Authors: Aniket Sharma, Lia Medeiros, Chi-kwan Chan, Goni Halevi, Patrick D. Mullen, James M. Stone, George N. Wong

Abstract: We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved space-times. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr metrics. JAX also enables Mahakala to run in parallel on both C… ▽ More We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved space-times. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr metrics. JAX also enables Mahakala to run in parallel on both CPUs and GPUs and achieve speeds comparable to C-based codes. Mahakala natively uses the Cartesian Kerr-Schild coordinate system, which avoids numerical issues caused by the "pole" of spherical coordinates. We demonstrate Mahakala's capabilities by simulating the 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature of Mahakala allows us to easily quantify the relative contribution of different regions of the flow to image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole. We also quantify the relative contribution of the disk, forward jet, and counter jet to 1.3 mm images. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 15 pages, 11 figures

The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity

Authors: Kengo Tomida, James M. Stone

Abstract: We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR hierarchy and TaskList framework of Athena++ for efficient parallelization. We adopt a conservative formulation for the Laplacian operator that avoids artificial accel… ▽ More We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR hierarchy and TaskList framework of Athena++ for efficient parallelization. We adopt a conservative formulation for the Laplacian operator that avoids artificial accelerations at level boundaries. Periodic, fixed, and zero-gradient boundary conditions are implemented, as well as open boundary conditions based on a multipole expansion. Hybrid parallelization using both MPI and OpenMP is adopted, and we present results of tests demonstrating the accuracy and scaling of the methods. On a uniform grid we show multigrid significantly outperforms methods based on FFTs, and requires only a small fraction of the compute time required by the (highly optimized) magnetohydrodynamic solver in Athena++. As a demonstration of the capabilities of the methods, we present the results of a test calculation of magnetized protostellar collapse on an adaptive mesh. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 28 pages, 13 figures, submitted to AAS Journals

An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization

Authors: Christopher J. White, Patrick D. Mullen, Yan-Fei Jiang, Shane W. Davis, James M. Stone, Viktoriya Morozova, Lizhong Zhang

Abstract: We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical… ▽ More We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems. △ Less

Submitted 28 March, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: To be published in ApJS

Report number: LA-UR-23-21167

Global Three-Dimensional Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar Mass Black Hole at Sub- and Near-critical Accretion Rates

Authors: Jiahui Huang, Yan-Fei Jiang, Hua Feng, Shane W. Davis, James M. Stone, Matthew J. Middleton

Abstract: We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundr… ▽ More We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 $r_{\rm g}$. The energy dissipation happens close to the mid-plane in the near-critical runs and near the disk surface in the low accretion rate run. The total radiative luminosity inside $\sim$20 $r_{\rm g}$ is about 1% to 30% the Eddington limit, with a radiative efficiency of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability (MRI) leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel with a terminal velocity of $\sim$0.1$c$ are seen only in the near-critical runs. We conclude that these magnetic pressure dominated disks are thermally stable and thicker than the $α$ disk, and the effective temperature profiles are much flatter than that in the $α$ disks. The magnetic pressure of these disks are comparable within an order of magnitude with the previous analytical magnetic pressure dominated disk model. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 17 pages, 13 figures, 3 tables, accepted for publication in ApJ

Measuring the variability of directly imaged exoplanets using vector Apodizing Phase Plates combined with ground-based differential spectrophotometry

Authors: Ben J. Sutlieff, Jayne L. Birkby, Jordan M. Stone, David S. Doelman, Matthew A. Kenworthy, Vatsal Panwar, Alexander J. Bohn, Steve Ertel, Frans Snik, Charles E. Woodward, Andrew J. Skemer, Jarron M. Leisenring, Klaus G. Strassmeier, David Charbonneau

Abstract: Clouds and other features in exoplanet and brown dwarf atmospheres cause variations in brightness as they rotate in and out of view. Ground-based instruments reach the high contrasts and small inner working angles needed to monitor these faint companions, but their small fields-of-view lack simultaneous photometric references to correct for non-astrophysical variations. We present a novel approach… ▽ More Clouds and other features in exoplanet and brown dwarf atmospheres cause variations in brightness as they rotate in and out of view. Ground-based instruments reach the high contrasts and small inner working angles needed to monitor these faint companions, but their small fields-of-view lack simultaneous photometric references to correct for non-astrophysical variations. We present a novel approach for making ground-based light curves of directly imaged companions using high-cadence differential spectrophotometric monitoring, where the simultaneous reference is provided by a double-grating 360° vector Apodizing Phase Plate (dgvAPP360) coronagraph. The dgvAPP360 enables high-contrast companion detections without blocking the host star, allowing it to be used as a simultaneous reference. To further reduce systematic noise, we emulate exoplanet transmission spectroscopy, where the light is spectrally-dispersed and then recombined into white-light flux. We do this by combining the dgvAPP360 with the infrared ALES integral field spectrograph on the Large Binocular Telescope Interferometer. To demonstrate, we observed the red companion HD 1160 B (separation ~780 mas) for one night, and detect $8.8\%$ semi-amplitude sinusoidal variability with a ~3.24 h period in its detrended white-light curve. We achieve the greatest precision in ground-based high-contrast imaging light curves of sub-arcsecond companions to date, reaching $3.7\%$ precision per 18-minute bin. Individual wavelength channels spanning 3.59-3.99 $μ$m further show tentative evidence of increasing variability with wavelength. We find no evidence yet of a systematic noise floor, hence additional observations can further improve the precision. This is therefore a promising avenue for future work aiming to map storms or find transiting exomoons around giant exoplanets. △ Less

Submitted 21 February, 2023; v1 submitted 20 January, 2023; originally announced January 2023.

Comments: 24 pages, 13 figures, in press in MNRAS; typos corrected, references updated

The Inner 2 pc of Sagittarius A*: Simulations of the Circumnuclear Disk and Multiphase Gas Accretion in the Galactic Center

Authors: Siddhant Solanki, Sean M. Ressler, Lena Murchikova, James M. Stone, Mark R. Morris

Abstract: We present hydrodynamic simulations of the inner few parsecs of the Milky Way's Galactic Center that, for the first time, combine a realistic treatment of stellar winds and the circumnuclear disk as they interact with the gravitational potential of the nuclear star cluster and Sagittarius~A*. We observe a complex interaction of the stellar winds with the inner edge of the circumnuclear disk, which… ▽ More We present hydrodynamic simulations of the inner few parsecs of the Milky Way's Galactic Center that, for the first time, combine a realistic treatment of stellar winds and the circumnuclear disk as they interact with the gravitational potential of the nuclear star cluster and Sagittarius~A*. We observe a complex interaction of the stellar winds with the inner edge of the circumnuclear disk, which leads to the growth of instabilities, induced accretion of cool gas from the inner edge of the disk, and the eventual formation of a small accretion disk of $\sim 10^4-10^5$ K within $r \sim 0.1$ pc. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 26 pages, 16 figures, (Submitted to ApJ)

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems: Best Practices for Data Collection in Cycle 2 and Beyond

Authors: Sasha Hinkley, Beth Biller, Andrew Skemer, Aarynn L. Carter, Julien Girard, Dean Hines, Jens Kammerer, Jarron Leisenring, William Balmer, Elodie Choquet, Maxwell A. Millar-Blanchaer, Marshall Perrin, Laurent Pueyo, Jason Wang, Kimberly Ward-Duong, Anthony Boccaletti, Brittany Miles, Polychronis Patapis, Isabel Rebollido, Emily Rickman, B. Sargent, Kadin Worthen, Kielan Hoch, Christine Chen, Stephanie Sallum , et al. (13 additional authors not shown)

Abstract: We present a set of recommended best practices for JWST data collection for members of the community focussed on the direct imaging and spectroscopy of exoplanetary systems. These findings and recommendations are based on the early analysis of the JWST Early Release Science Program 1386, "High-Contrast Imaging of Exoplanets and Exoplanetary Systems with JWST." Our goal is for this information to b… ▽ More We present a set of recommended best practices for JWST data collection for members of the community focussed on the direct imaging and spectroscopy of exoplanetary systems. These findings and recommendations are based on the early analysis of the JWST Early Release Science Program 1386, "High-Contrast Imaging of Exoplanets and Exoplanetary Systems with JWST." Our goal is for this information to be useful for observers in preparation of JWST proposals for Cycle 2 and beyond. In addition to compiling a set of best practices from our ERS program, in a few cases we also draw on the expertise gained within the instrument commissioning programs, as well as include a handful of data processing best practices. We anticipate that this document will be regularly updated and resubmitted to arXiv.org to ensure that we have distributed our knowledge of best-practices for data collection as widely and efficiently as possible. △ Less

Submitted 25 January, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Not yet submitted for publication. Intended only to be a community resource for JWST Cycle 2 proposals

Toward Horizon-scale Accretion Onto Supermassive Black Holes in Elliptical Galaxies

Authors: Minghao Guo, James M. Stone, Chang-Goo Kim, Eliot Quataert

Abstract: We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and span more than 6 orders of magnitude in radius, reaching scales similar to the black hole horizon. The key p… ▽ More We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and span more than 6 orders of magnitude in radius, reaching scales similar to the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk and relatively rare chaotic stages ($\lesssim 10\%$ of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as $\dot{M}\propto r^{1/2}$ when hot gas dominates and we obtain $\dot{M}\simeq10^\mathrm{-4}-10^\mathrm{-3}\,M_\odot\,\mathrm{yr^{-1}}$ near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by $\sim (r_\mathrm{g}/r_{\rm Bondi})^{1/2}$. Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale. △ Less

Submitted 29 March, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Comments: 28 pages, 25 figures, 2 tables, accepted for publication in ApJ

Prograde and Retrograde Gas Flow around Disk-embedded Companions: Dependence on Eccentricity, Mass and Disk Properties

Authors: Yi-Xian Chen, Avery Bailey, James M. Stone, Zhaohuan Zhu

Abstract: We apply 3D hydrodynamical simulations to study the rotational aspect of gas flow patterns around eccentric companions embedded in an accretion disk around its primary host. We sample a wide range of companion mass ratio q and disk aspect ratio h, and confirm a generic transition from prograde (steady tidal interaction dominated) to retrograde (background Keplerian shear dominated) circum-companio… ▽ More We apply 3D hydrodynamical simulations to study the rotational aspect of gas flow patterns around eccentric companions embedded in an accretion disk around its primary host. We sample a wide range of companion mass ratio q and disk aspect ratio h, and confirm a generic transition from prograde (steady tidal interaction dominated) to retrograde (background Keplerian shear dominated) circum-companion flow when orbital eccentricity exceeds a critical value et. We find et \sim h for sub-thermal companions while et \sim (q/h)^1/3 for super-thermal companions, and propose an empirical formula to unify the two scenarios. Our results also suggest that et is insensitive to modest levels of turbulence, modeled in the form of a kinematic viscosity term. In the context of stellar-mass Black Holes (sBHs) embedded in AGN accretion disks, the bifurcation of their circum-stellar disk (CSD) rotation suggest the formation of a population of nearly anti-aligned sBHs, whose relevance to low spin gravitational waves (GW) events can be probed in more details with future population models of sBH evolution in AGN disks, making use of our quantitative scaling for et; In the context of circum-planetary disks (CPDs), our results suggest the possibility of forming retrograde satellites in-situ in retrograde CPDs around eccentric planets. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: Accepted to ApJL, 10 pages, 5 figures

A Tale of a Tail: A Tidally-Disrupting Ultra-Diffuse Galaxy in the M81 Group

Authors: Rokas Žemaitis, Annette M. N. Ferguson, Sakurako Okamoto, Jean-Charles Cuillandre, Connor J. Stone, Nobuo Arimoto, Mike J. Irwin

Abstract: We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper… ▽ More We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper Suprime-Cam survey of the M81 Group. It has an average surface brightness of $μ_g \sim 32$ mag arcsec$^{-2}$ and can be traced for over a degree on the sky (60 kpc at the distance of F8D1) with our current imagery. We revisit the main body properties of F8D1 using deep multiband imagery acquired with MegaCam on CFHT and measure effective radii of $1.7-1.9$ kpc, central surface brightnesses of $24.7-25.7$ mag and a stellar mass of $\sim7 \times 10^7 M_{\odot}$. Assuming a symmetric feature on the other side of the galaxy, we calculate that $30-36$% of F8D1's present-day luminosity is contained in the tail. We argue that the most likely origin of F8D1's disruption is a recent close passage to M81, which would have stripped its gas and quenched its star formation. As the only UDG that has so far been studied to such faint surface brightness depths, the unveiling of F8D1's tidal disruption is important. It leaves open the possibility that many other UDGs could be the result of similar processes, with the most telling signatures of this lurking below current detection limits. △ Less

Submitted 27 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. 16 pages, 12 figures, 3 tables

Neutron Tagging following Atmospheric Neutrino Events in a Water Cherenkov Detector

Authors: K. Abe, Y. Haga, Y. Hayato, K. Hiraide, K. Ieki, M. Ikeda, S. Imaizumi, K. Iyogi, J. Kameda, Y. Kanemura, Y. Kataoka, Y. Kato, Y. Kishimoto, S. Miki, S. Mine, M. Miura, T. Mochizuki, S. Moriyama, Y. Nagao, M. Nakahata, T. Nakajima, Y. Nakano, S. Nakayama, T. Okada, K. Okamoto , et al. (281 additional authors not shown)

Abstract: We present the development of neutron-tagging techniques in Super-Kamiokande IV using a neural network analysis. The detection efficiency of neutron capture on hydrogen is estimated to be 26%, with a mis-tag rate of 0.016 per neutrino event. The uncertainty of the tagging efficiency is estimated to be 9.0%. Measurement of the tagging efficiency with data from an Americium-Beryllium calibration agr… ▽ More We present the development of neutron-tagging techniques in Super-Kamiokande IV using a neural network analysis. The detection efficiency of neutron capture on hydrogen is estimated to be 26%, with a mis-tag rate of 0.016 per neutrino event. The uncertainty of the tagging efficiency is estimated to be 9.0%. Measurement of the tagging efficiency with data from an Americium-Beryllium calibration agrees with this value within 10%. The tagging procedure was performed on 3,244.4 days of SK-IV atmospheric neutrino data, identifying 18,091 neutrons in 26,473 neutrino events. The fitted neutron capture lifetime was measured as 218 \pm 9 μs. △ Less

Submitted 20 September, 2022; v1 submitted 18 September, 2022; originally announced September 2022.

Journal ref: JINST 17 P10029 (2022)

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems II: A 1 to 20 Micron Spectrum of the Planetary-Mass Companion VHS 1256-1257 b

Authors: Brittany E. Miles, Beth A. Biller, Polychronis Patapis, Kadin Worthen, Emily Rickman, Kielan K. W. Hoch, Andrew Skemer, Marshall D. Perrin, Niall Whiteford, Christine H. Chen, B. Sargent, Sagnick Mukherjee, Caroline V. Morley, Sarah E. Moran, Mickael Bonnefoy, Simon Petrus, Aarynn L. Carter, Elodie Choquet, Sasha Hinkley, Kimberly Ward-Duong, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Laurent Pueyo, Shrishmoy Ray, Karl R. Stapelfeldt , et al. (79 additional authors not shown)

Abstract: We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude… ▽ More We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $μ$m to 20 $μ$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion. △ Less

Submitted 4 July, 2024; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems I: High Contrast Imaging of the Exoplanet HIP 65426 b from 2-16 $μ$m

Authors: Aarynn L. Carter, Sasha Hinkley, Jens Kammerer, Andrew Skemer, Beth A. Biller, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Simon Petrus, Jordan M. Stone, Kimberly Ward-Duong, Jason J. Wang, Julien H. Girard, Dean C. Hines, Marshall D. Perrin, Laurent Pueyo, William O. Balmer, Mariangela Bonavita, Mickael Bonnefoy, Gael Chauvin, Elodie Choquet, Valentin Christiaens, Camilla Danielski, Grant M. Kennedy, Elisabeth C. Matthews, Brittany E. Miles , et al. (86 additional authors not shown)

Abstract: We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exo… ▽ More We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $μ$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$σ$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $μ$m and MIRI at 11.3 $μ$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $μ$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail. △ Less

Submitted 3 May, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Design of SCALES: A 2-5 Micron Coronagraphic Integral Field Spectrograph for Keck Observatory

Authors: Andrew Skemer, R. Deno Stelter, Stephanie Sallum, Nicholas MacDonald, Renate Kupke, Christopher Ratliffe, Ravinder Banyal, Amirul Hasan, Hari Mohan Varshney, Arun Surya, Ajin Prakash, Sivarani Thirupathi, Ramya Sethuraman, Govinda K. V., Michael P. Fitzgerald, Eric Wang, Marc Kassis, Olivier Absil, Carlos Alvarez, Natasha Batalha, Marc-Andre Boucher, Cyril Bourgenot, Timothy Brandt, Zackery Briesemeister, Katherine de Kleer , et al. (27 additional authors not shown)

Abstract: We present the design of SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) a new 2-5 micron coronagraphic integral field spectrograph under construction for Keck Observatory. SCALES enables low-resolution (R~50) spectroscopy, as well as medium-resolution (R~4,000) spectroscopy with the goal of discovering and characterizing cold exoplanets that are brightest in the thermal… ▽ More We present the design of SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) a new 2-5 micron coronagraphic integral field spectrograph under construction for Keck Observatory. SCALES enables low-resolution (R~50) spectroscopy, as well as medium-resolution (R~4,000) spectroscopy with the goal of discovering and characterizing cold exoplanets that are brightest in the thermal infrared. Additionally, SCALES has a 12x12" field-of-view imager that will be used for general adaptive optics science at Keck. We present SCALES's specifications, its science case, its overall design, and simulations of its expected performance. Additionally, we present progress on procuring, fabricating and testing long lead-time components. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 19 pages, 12 Figures, SPIE Astronomical Telescopes + Instrumentation, (#12184-18), 2022, Montreal, Canada

Turbulent Magnetic Field Amplification by the Interaction of Shock Wave and Inhomogeneous Medium

Authors: Yue Hu, Siyao Xu, James M. Stone, Alex Lazarian

Abstract: Magnetic fields on the order of 100 $μ$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven b… ▽ More Magnetic fields on the order of 100 $μ$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven by the strongest preshock density contrast and follows the Kolmogorov scaling. The resulting turbulence amplifies the postshock magnetic field. The time evolution of the magnetic fields agrees with the prediction of the nonlinear turbulent dynamo theory in Xu & Lazarian (2016). When the initial weak magnetic field is perpendicular to the shock normal, the maximum amplification of the field's strength achieves a factor of $\approx200$, which is twice larger than that for a parallel shock. We find that the perpendicular shock exhibits a smaller turbulent Alfvén Mach number in the vicinity of the shock front than the parallel shock. However, the strongest magnetic field has a low volume filling factor and is limited by the turbulent energy due to the reconnection diffusion taking place in a turbulent and magnetized fluid. The magnetic field strength averaged along the $z$-axis is reduced by a factor $\gtrsim10$. We decompose the turbulent velocity and magnetic field into solenoidal and compressive modes. The solenoidal mode is dominant and evolves to follow the Kolmogorov scaling, even though the preshock density distribution has a shallow spectrum. When the preshock density distribution has a Kolmogorov spectrum, the turbulent velocity's compressive component increases. △ Less

Submitted 9 November, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 19 pages, 15 figures, accepted for publication in ApJ

Search for supernova bursts in Super-Kamiokande IV

Authors: The Super-Kamiokande collaboration, :, M. Mori, K. Abe, Y. Hayato, K. Hiraide, K. Ieki, M. Ikeda, S. Imaizumi, J. Kameda, Y. Kanemura, R. Kaneshima, Y. Kashiwagi, Y. Kataoka, S. Miki, S. Mine, M. Miura, S. Moriyama, Y. Nagao, M. Nakahata, Y. Nakano, S. Nakayama, Y. Noguchi, T. Okada, K. Okamoto , et al. (223 additional authors not shown)

Abstract: Super-Kamiokande has been searching for neutrino bursts characteristic of core-collapse supernovae continuously, in real time, since the start of operations in 1996. The present work focuses on detecting more distant supernovae whose event rate may be too small to trigger in real time, but may be identified using an offline approach. The analysis of data collected from 2008 to 2018 found no eviden… ▽ More Super-Kamiokande has been searching for neutrino bursts characteristic of core-collapse supernovae continuously, in real time, since the start of operations in 1996. The present work focuses on detecting more distant supernovae whose event rate may be too small to trigger in real time, but may be identified using an offline approach. The analysis of data collected from 2008 to 2018 found no evidence of distant supernovae bursts. This establishes an upper limit of 0.29 year$^{-1}$ on the rate of core-collapse supernovae out to 100 kpc at 90% C.L.. For supernovae that fail to explode and collapse directly to black holes the limit reaches to 300 kpc. △ Less

Submitted 2 June, 2022; originally announced June 2022.

The JWST Early Release Science Program for the Direct Imaging & Spectroscopy of Exoplanetary Systems

Authors: Sasha Hinkley, Aarynn L. Carter, Shrishmoy Ray, Andrew Skemer, Beth Biller, Elodie Choquet, Maxwell A. Millar-Blanchaer, Stephanie Sallum, Brittany Miles, Niall Whiteford, Polychronis Patapis, Marshall D. Perrin, Laurent Pueyo, Glenn Schneider, Karl Stapelfeldt, Jason Wang, Kimberly Ward-Duong, Brendan P. Bowler, Anthony Boccaletti, Julien H. Girard, Dean Hines, Paul Kalas, Jens Kammerer, Pierre Kervella, Jarron Leisenring , et al. (61 additional authors not shown)

Abstract: The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe e… ▽ More The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$μ$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$μ$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles. △ Less

Submitted 12 September, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 21 pages, 7 figures. Accepted for Publication in PASP

Testing Non-Standard Interactions Between Solar Neutrinos and Quarks with Super-Kamiokande

Authors: Super-Kamiokande Collaboration, :, P. Weatherly, K. Abe, C. Bronner, Y. Hayato, K. Hiraide, M. Ikeda, K. Iyogi, J. Kameda, Y. Kanemura, Y. Kataoka, Y. Kato, Y. Kishimoto, S. Miki, M. Miura, S. Moriyama, T. Mochizuki, M. Nakahata, Y. Nakano, S. Nakayama, T. Okada, K. Okamoto, A. Orii, G. Pronost , et al. (248 additional authors not shown)

Abstract: Non-Standard Interactions (NSI) between neutrinos and matter affect the neutrino flavor oscillations. Due to the high matter density in the core of the Sun, solar neutrinos are suited to probe these interactions. Using the $277$ kton-yr exposure of Super-Kamiokande to $^{8}$B solar neutrinos, we search for the presence of NSI. Our data favors the presence of NSI with down quarks at 1.8$σ$, and wit… ▽ More Non-Standard Interactions (NSI) between neutrinos and matter affect the neutrino flavor oscillations. Due to the high matter density in the core of the Sun, solar neutrinos are suited to probe these interactions. Using the $277$ kton-yr exposure of Super-Kamiokande to $^{8}$B solar neutrinos, we search for the presence of NSI. Our data favors the presence of NSI with down quarks at 1.8$σ$, and with up quarks at 1.6$σ$, with the best fit NSI parameters being ($ε_{11}^{d},ε_{12}^{d}$) = (-3.3, -3.1) for $d$-quarks and ($ε_{11}^{u},ε_{12}^{u}$) = (-2.5, -3.1) for $u$-quarks. After combining with data from the Sudbury Neutrino Observatory and Borexino, the significance increases by 0.1$σ$. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Author: Pierce Weatherly 25 pages. To be submitted to Physical Review D

L-band Integral Field Spectroscopy of the HR 8799 Planetary System

Authors: David S. Doelman, Jordan M. Stone, Zackery W. Briesemeister, Andrew J. I. Skemer, Travis Barman, Laci S. Brock, Philip M. Hinz, Alexander Bohn, Matthew Kenworthy, Sebastiaan Y. Haffert, Frans Snik, Steve Ertel, Jarron M. Leisenring, Charles E. Woodward, Michael F. Skrutskie

Abstract: Understanding the physical processes sculpting the appearance of young gas-giant planets is complicated by degeneracies confounding effective temperature, surface gravity, cloudiness, and chemistry. To enable more detailed studies, spectroscopic observations covering a wide range of wavelengths is required. Here we present the first L-band spectroscopic observations of HR 8799 d and e and the firs… ▽ More Understanding the physical processes sculpting the appearance of young gas-giant planets is complicated by degeneracies confounding effective temperature, surface gravity, cloudiness, and chemistry. To enable more detailed studies, spectroscopic observations covering a wide range of wavelengths is required. Here we present the first L-band spectroscopic observations of HR 8799 d and e and the first low-resolution wide bandwidth L-band spectroscopic measurements of HR 8799 c. These measurements were facilitated by an upgraded LMIRCam/ALES instrument at the LBT, together with a new apodizing phase plate coronagraph. Our data are generally consistent with previous photometric observations covering similar wavelengths, yet there exists some tension with narrowband photometry for HR 8799 c. With the addition of our spectra, each of the three innermost observed planets in the HR 8799 system have had their spectral energy distributions measured with integral field spectroscopy covering $\sim0.9$ to $4.1~μ\mathrm{m}$. We combine these spectra with measurements from the literature and fit synthetic model atmospheres. We demonstrate that the bolometric luminosity of the planets is not sensitive to the choice of model atmosphere used to interpolate between measurements and extrapolate beyond them. Combining luminosity with age and mass constraints, we show that the predictions of evolutionary models are narrowly peaked for effective temperature, surface gravity, and planetary radius. By holding these parameters at their predicted values, we show that more flexible cloud models can provide good fits to the data while being consistent with the expectations of evolutionary models. △ Less

Submitted 8 April, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 19 pages, 11 figures, accepted for publication in The Astronomical Journal; added reference, updated figure 6 and table 4

Parthenon -- a performance portable block-structured adaptive mesh refinement framework

Authors: Philipp Grete, Joshua C. Dolence, Jonah M. Miller, Joshua Brown, Ben Ryan, Andrew Gaspar, Forrest Glines, Sriram Swaminarayan, Jonas Lippuner, Clell J. Solomon, Galen Shipman, Christoph Junghans, Daniel Holladay, James M. Stone, Luke F. Roberts

Abstract: On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived… ▽ More On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multi-dimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of $1.7\times10^{13}$ zone-cycles/s on 9,216 nodes (73,728 logical GPUs) at ~92% weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR. △ Less

Submitted 21 November, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Comments: 17 pages, 11 figures, accepted for publication in IJHPCA, Codes available at https://github.com/parthenon-hpc-lab

Report number: LA-UR-22-21270

Planet populations inferred from debris discs: insights from 178 debris systems in the ISPY, LEECH and LIStEN planet-hunting surveys

Authors: Tim D. Pearce, Ralf Launhardt, Robert Ostermann, Grant M. Kennedy, Mario Gennaro, Mark Booth, Alexander V. Krivov, Gabriele Cugno, Thomas K. Henning, Andreas Quirrenbach, Arianna Musso Barcucci, Elisabeth C. Matthews, Henrik L. Ruh, Jordan M. Stone

Abstract: We know little about the outermost exoplanets in planetary systems, because our detection methods are insensitive to moderate-mass planets on wide orbits. However, debris discs can probe the outer-planet population, because dynamical modelling of observed discs can reveal properties of perturbing planets. We use four sculpting and stirring arguments to infer planet properties in 178 debris-disc sy… ▽ More We know little about the outermost exoplanets in planetary systems, because our detection methods are insensitive to moderate-mass planets on wide orbits. However, debris discs can probe the outer-planet population, because dynamical modelling of observed discs can reveal properties of perturbing planets. We use four sculpting and stirring arguments to infer planet properties in 178 debris-disc systems from the ISPY, LEECH and LIStEN planet-hunting surveys. Similar analyses are often conducted for individual discs, but we consider a large sample in a consistent manner. We aim to predict the population of wide-separation planets, gain insight into the formation and evolution histories of planetary systems, and determine the feasibility of detecting these planets in the near future. We show that a `typical' cold debris disc likely requires a Neptune- to Saturn-mass planet at 10-100 au, with some needing Jupiter-mass perturbers. Our predicted planets are currently undetectable, but modest detection-limit improvements (e.g. from JWST) should reveal many such perturbers. We find that planets thought to be perturbing debris discs at late times are similar to those inferred to be forming in protoplanetary discs, so these could be the same population if newly formed planets do not migrate as far as currently thought. Alternatively, young planets could rapidly sculpt debris before migrating inwards, meaning that the responsible planets are more massive (and located further inwards) than debris-disc studies assume. We combine self-stirring and size-distribution modelling to show that many debris discs cannot be self-stirred without having unreasonably high masses; planet- or companion-stirring may therefore be the dominant mechanism in many (perhaps all) debris discs. Finally, we provide catalogues of planet predictions, and identify promising targets for future planet searches. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 41 pages, 16 figures, accepted for publication in A&A

New Methods and Simulations for Cosmogenic Induced Spallation Removal in Super-Kamiokande-IV

Authors: Super-Kamiokande Collaboration, :, S. Locke, A. Coffani, K. Abe, C. Bronner, Y. Hayato, M. Ikeda, S. Imaizumi, H. Ito, J. Kameda, Y. Kataoka, M. Miura, S. Moriyama, Y. Nagao, M. Nakahata, Y. Nakajima, S. Nakayama, T. Okada, K. Okamoto, A. Orii, G. Pronost, H. Sekiya, M. Shiozawa, Y. Sonoda , et al. (196 additional authors not shown)

Abstract: Radioactivity induced by cosmic muon spallation is a dominant source of backgrounds for $\mathcal{O}(10)~$MeV neutrino interactions in water Cherenkov detectors. In particular, it is crucial to reduce backgrounds to measure the solar neutrino spectrum and find neutrino interactions from distant supernovae. In this paper we introduce new techniques to locate muon-induced hadronic showers and effici… ▽ More Radioactivity induced by cosmic muon spallation is a dominant source of backgrounds for $\mathcal{O}(10)~$MeV neutrino interactions in water Cherenkov detectors. In particular, it is crucial to reduce backgrounds to measure the solar neutrino spectrum and find neutrino interactions from distant supernovae. In this paper we introduce new techniques to locate muon-induced hadronic showers and efficiently reject spallation backgrounds. Applying these techniques to the solar neutrino analysis with an exposure of $2790\times22.5$~kton.day increases the signal efficiency by $12.6\%$, approximately corresponding to an additional year of detector running. Furthermore, we present the first spallation simulation at SK, where we model hadronic interactions using FLUKA. The agreement between the isotope yields and shower pattern in this simulation and in the data gives confidence in the accuracy of this simulation, and thus opens the door to use it to optimize muon spallation removal in new data with gadolinium-enhanced neutron capture detection. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Diffuse Supernova Neutrino Background Search at Super-Kamiokande

Authors: Super-Kamiokande Collaboration, :, K. Abe, C. Bronner, Y. Hayato, K. Hiraide, M. Ikeda, S. Imaizumi, J. Kameda, Y. Kanemura, Y. Kataoka, S. Miki, M. Miura, S. Moriyama, Y. Nagao, M. Nakahata, S. Nakayama, T. Okada, K. Okamoto, A. Orii, G. Pronost, H. Sekiya, M. Shiozawa, Y. Sonoda, Y. Suzuki , et al. (197 additional authors not shown)

Abstract: A new search for the diffuse supernova neutrino background (DSNB) flux has been conducted at Super-Kamiokande (SK), with a $22.5\times2970$-kton$\cdot$day exposure from its fourth operational phase IV. The new analysis improves on the existing background reduction techniques and systematic uncertainties and takes advantage of an improved neutron tagging algorithm to lower the energy threshold comp… ▽ More A new search for the diffuse supernova neutrino background (DSNB) flux has been conducted at Super-Kamiokande (SK), with a $22.5\times2970$-kton$\cdot$day exposure from its fourth operational phase IV. The new analysis improves on the existing background reduction techniques and systematic uncertainties and takes advantage of an improved neutron tagging algorithm to lower the energy threshold compared to the previous phases of SK. This allows for setting the world's most stringent upper limit on the extraterrestrial $\barν_e$ flux, for neutrino energies below 31.3 MeV. The SK-IV results are combined with the ones from the first three phases of SK to perform a joint analysis using $22.5\times5823$ kton$\cdot$days of data. This analysis has the world's best sensitivity to the DSNB $\barν_e$ flux, comparable to the predictions from various models. For neutrino energies larger than 17.3 MeV, the new combined $90\%$ C.L. upper limits on the DSNB $\barν_e$ flux lie around $2.7$ cm$^{-2}$$\cdot$$\text{sec}^{-1}$, strongly disfavoring the most optimistic predictions. Finally, potentialities of the gadolinium phase of SK and the future Hyper-Kamiokande experiment are discussed. △ Less

Submitted 2 November, 2021; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: 42 pages, 37 figures, 14 tables

First Gadolinium Loading to Super-Kamiokande

Authors: K. Abe, C. Bronner, Y. Hayato, K. Hiraide, M. Ikeda, S. Imaizumi, J. Kameda, Y. Kanemura, Y. Kataoka, S. Miki, M. Miura, S. Moriyama, Y. Nagao, M. Nakahata, S. Nakayama, T. Okada, K. Okamoto, A. Orii, G. Pronost, H. Sekiya, M. Shiozawa, Y. Sonoda, Y. Suzuki, A. Takeda, Y. Takemoto , et al. (192 additional authors not shown)

Abstract: In order to improve Super-Kamiokande's neutron detection efficiency and to thereby increase its sensitivity to the diffuse supernova neutrino background flux, 13 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ (gadolinium sulfate octahydrate) was dissolved into the detector's otherwise ultrapure water from July 14 to August 17, 2020, marking the start of the SK-Gd phase of operations. During the loa… ▽ More In order to improve Super-Kamiokande's neutron detection efficiency and to thereby increase its sensitivity to the diffuse supernova neutrino background flux, 13 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ (gadolinium sulfate octahydrate) was dissolved into the detector's otherwise ultrapure water from July 14 to August 17, 2020, marking the start of the SK-Gd phase of operations. During the loading, water was continuously recirculated at a rate of 60 m$^3$/h, extracting water from the top of the detector and mixing it with concentrated $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ solution to create a 0.02% solution of the Gd compound before injecting it into the bottom of the detector. A clear boundary between the Gd-loaded and pure water was maintained through the loading, enabling monitoring of the loading itself and the spatial uniformity of the Gd concentration over the 35 days it took to reach the top of the detector. During the subsequent commissioning the recirculation rate was increased to 120 m$^3$/h, resulting in a constant and uniform distribution of Gd throughout the detector and water transparency equivalent to that of previous pure-water operation periods. Using an Am-Be neutron calibration source the mean neutron capture time was measured to be $115\pm1$ $μ$s, which corresponds to a Gd concentration of $111\pm2$ ppm, as expected for this level of Gd loading. This paper describes changes made to the water circulation system for this detector upgrade, the Gd loading procedure, detector commissioning, and the first neutron calibration measurements in SK-Gd. △ Less

Submitted 15 December, 2021; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: 37 pages, 19 Figures, Accepted for publication in Nucl. Instrum. Meth. A

Journal ref: Nuclear Inst. and Methods in Physics Research, A 1027 (2022) 166248

Imaging low-mass planets within the habitable zones of nearby stars with ground-based mid-infrared imaging

Authors: Kevin Wagner, Steve Ertel, Jordan Stone, Jarron Leisenring, Dániel Apai, Markus Kasper, Olivier Absil, Laird Close, Denis Defrère, Olivier Guyon, Jared Males

Abstract: Giant exoplanets on 10-100 au orbits have been directly imaged around young stars. The peak of the thermal emission from these warm young planets is in the near-infrared (~1-5 microns), whereas mature, temperate exoplanets (i.e., those within their stars' habitable zones) radiate primarily in the mid-infrared (mid-IR: ~10 microns). If the background noise in the mid-IR can be mitigated, then exopl… ▽ More Giant exoplanets on 10-100 au orbits have been directly imaged around young stars. The peak of the thermal emission from these warm young planets is in the near-infrared (~1-5 microns), whereas mature, temperate exoplanets (i.e., those within their stars' habitable zones) radiate primarily in the mid-infrared (mid-IR: ~10 microns). If the background noise in the mid-IR can be mitigated, then exoplanets with low masses--including rocky exoplanets--can potentially be imaged in very deep exposures. Here, we review the recent results of the Breakthrough Watch/New Earths in the Alpha Centauri Region (NEAR) program on the Very Large Telescope (VLT) in Chile. NEAR pioneered a ground-based mid-IR observing approach designed to push the capabilities for exoplanet imaging with a specific focus on the closest stellar system, Alpha Centauri. NEAR combined several new optical technologies--including a mid-IR optimized coronagraph, adaptive optics system, and rapid chopping strategy to mitigate noise from the central star and thermal background within the habitable zone. We focus on the lessons of the VLT/NEAR campaign to improve future instrumentation--specifically on strategies to improve noise mitigation through chopping. We also present the design and commissioning of the Large Binocular Telescope's Exploratory Survey for Super-Earths Orbiting Nearby Stars (LESSONS), an experiment in the Northern hemisphere that is building on what was learned from NEAR to further push the sensitivity of mid-IR imaging. Finally, we briefly discuss some of the possibilities that mid-IR imaging will enable for exoplanet science. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Submitted proceedings to SPIE Optical Engineering + Applications 2021, Techniques and Instrumentation for Detection of Exoplanets X

Cloud Properties of Brown Dwarf Binaries Across the L/T Transition

Authors: Laci Shea Brock, Travis Barman, Quinn M. Konopacky, Jordan M. Stone

Abstract: We present a new suite of atmosphere models with flexible cloud parameters to investigate the effects of clouds on brown dwarfs across the L/T transition. We fit these models to a sample of 13 objects with well-known masses, distances, and spectral types spanning L3-T5. Our modelling is guided by spatially-resolved photometry from the Hubble Space Telescope and the W. M. Keck Telescopes covering v… ▽ More We present a new suite of atmosphere models with flexible cloud parameters to investigate the effects of clouds on brown dwarfs across the L/T transition. We fit these models to a sample of 13 objects with well-known masses, distances, and spectral types spanning L3-T5. Our modelling is guided by spatially-resolved photometry from the Hubble Space Telescope and the W. M. Keck Telescopes covering visible to near-infrared wavelengths. We find that, with appropriate cloud parameters, the data can be fit well by atmospheric models with temperature and surface gravity in agreement with the predictions of evolutionary models. We see a clear trend in the cloud parameters with spectral type, with earlier-type objects exhibiting higher-altitude clouds with smaller grains (0.25-0.50 micron) and later-type objects being better fit with deeper clouds and larger grains ($\geq$1 micron). Our results confirm previous work that suggests L dwarfs are dominated by submicron particles, whereas T dwarfs have larger particle sizes. △ Less

Submitted 17 June, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Glassy quantum nuclear pasta in neutron star crusts

Authors: William G. Newton, Mark Alexander Kaltenborn, Sarah Cantu, Shuxi Wang, Amber Stinson, Jirina Rikovska Stone

Abstract: We investigate the nuclear pasta phases in neutron star crusts by conducting a large number of three-dimensional Hartree-Fock+BCS calculations at densities leading to the crust-core transition. We survey the shape parameter space of pasta at constant pressure. Spaghetti, waffles, lasagna, bi-continuous phases and cylindrical holes occupy local minima in the resulting Gibbs energy surfaces. The bi-… ▽ More We investigate the nuclear pasta phases in neutron star crusts by conducting a large number of three-dimensional Hartree-Fock+BCS calculations at densities leading to the crust-core transition. We survey the shape parameter space of pasta at constant pressure. Spaghetti, waffles, lasagna, bi-continuous phases and cylindrical holes occupy local minima in the resulting Gibbs energy surfaces. The bi-continuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions and therefore protons are delocalized, appears over a large range of depths. Our results support the idea that nuclear pasta is a glassy system. Multiple pasta configurations coexist in a given layer of the crust. At a characteristic temperature, of order $10^8$-$10^9$K, different phases become frozen into domains whose sizes we estimate to be 1-50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, there is very little long-range order and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature dependent transport properties in the deep layers of the crust. We identify 4 distinct regions: (1) nuclear pasta first appears as a local minima, but spherical nuclei are the ground state; (2) nuclear pasta become the absolute minimum, but spherical nuclei are still a local minimum (3) only nuclear pasta appears in local minima, and protons are still localized in at least one dimension (4) only pasta appears, and protons are delocalized. The whole pasta region can occupy up to 70% of the crust by mass and 40% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 25% of its thickness. △ Less

Submitted 23 April, 2021; originally announced April 2021.

Comments: 32 pages, 23 figures, submitted to Physical Review C

The vector-apodizing phase plate coronagraph: design, current performance, and future development

Authors: D. S. Doelman, F. Snik, E. H. Por, S. P. Bos, G. P. P. L. Otten, M. Kenworthy, S. Y. Haffert, M. Wilby, A. J. Bohn, B. J. Sutlieff, K. Miller, M. Ouellet, J. de Boer, C. U. Keller, M. J. Escuti, S. Shi, N. Z. Warriner, K. J. Hornburg, J. L. Birkby, J. Males, K. M. Morzinski, L. M. Close, J. Codona, J. Long, L. Schatz , et al. (28 additional authors not shown)

Abstract: Over the last decade, the vector-apodizing phase plate (vAPP) coronagraph has been developed from concept to on-sky application in many high-contrast imaging systems on 8-m class telescopes. The vAPP is an geometric-phase patterned coronagraph that is inherently broadband, and its manufacturing is enabled only by direct-write technology for liquid-crystal patterns. The vAPP generates two coronagra… ▽ More Over the last decade, the vector-apodizing phase plate (vAPP) coronagraph has been developed from concept to on-sky application in many high-contrast imaging systems on 8-m class telescopes. The vAPP is an geometric-phase patterned coronagraph that is inherently broadband, and its manufacturing is enabled only by direct-write technology for liquid-crystal patterns. The vAPP generates two coronagraphic PSFs that cancel starlight on opposite sides of the point spread function (PSF) and have opposite circular polarization states. The efficiency, that is the amount of light in these PSFs, depends on the retardance offset from half-wave of the liquid-crystal retarder. Using different liquid-crystal recipes to tune the retardance, different vAPPs operate with high efficiencies ($>96\%$) in the visible and thermal infrared (0.55 $μ$m to 5 $μ$m). Since 2015, seven vAPPs have been installed in a total of six different instruments, including Magellan/MagAO, Magellan/MagAO-X, Subaru/SCExAO, and LBT/LMIRcam. Using two integral field spectrographs installed on the latter two instruments, these vAPPs can provide low-resolution spectra (R$\sim$30) between 1 $μ$m and 5 $μ$m. We review the design process, development, commissioning, on-sky performance, and first scientific results of all commissioned vAPPs. We report on the lessons learned and conclude with perspectives for future developments and applications. △ Less

Submitted 4 November, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 38 pages, 17 figures, accepted for publication in Applied Optics, added NSF grant acknowledgement

Search for neutrinos in coincidence with gravitational wave events from the LIGO-Virgo O3a Observing Run with the Super-Kamiokande detector

Authors: The Super-Kamiokande collaboration, :, K. Abe, C. Bronner, Y. Hayato, M. Ikeda, S. Imaizumi, J. Kameda, Y. Kanemura, Y. Kataoka, S. Miki, M. Miura, S. Moriyama, Y. Nagao, M. Nakahata, S. Nakayama, T. Okada, K. Okamoto, A. Orii, G. Pronost, H. Sekiya, M. Shiozawa, Y. Sonoda, Y. Suzuki, A. Takeda , et al. (189 additional authors not shown)

Abstract: The Super-Kamiokande detector can be used to search for neutrinos in time coincidence with gravitational waves detected by the LIGO-Virgo Collaboration (LVC). Both low-energy ($7-100$ MeV) and high-energy ($0.1-10^5$ GeV) samples were analyzed in order to cover a very wide neutrino spectrum. Follow-ups of 36 (out of 39) gravitational waves reported in the GWTC-2 catalog were examined; no significa… ▽ More The Super-Kamiokande detector can be used to search for neutrinos in time coincidence with gravitational waves detected by the LIGO-Virgo Collaboration (LVC). Both low-energy ($7-100$ MeV) and high-energy ($0.1-10^5$ GeV) samples were analyzed in order to cover a very wide neutrino spectrum. Follow-ups of 36 (out of 39) gravitational waves reported in the GWTC-2 catalog were examined; no significant excess above the background was observed, with 10 (24) observed neutrinos compared with 4.8 (25.0) expected events in the high-energy (low-energy) samples. A statistical approach was used to compute the significance of potential coincidences. For each observation, p-values were estimated using neutrino direction and LVC sky map ; the most significant event (GW190602_175927) is associated with a post-trial p-value of $7.8\%$ ($1.4σ$). Additionally, flux limits were computed independently for each sample and by combining the samples. The energy emitted as neutrinos by the identified gravitational wave sources was constrained, both for given flavors and for all-flavors assuming equipartition between the different flavors, independently for each trigger and by combining sources of the same nature. △ Less

Submitted 13 September, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 16 pages, 5 figures. v2: adding corrections from The Astrophysical Journal review

Journal ref: The Astrophysical Journal, Volume 918, Number 2 (2021)

The HOSTS survey: evidence for an extended dust disk and constraints on the presence of giant planets in the Habitable Zone of $β$ Leo

Authors: D. Defrère, P. M. Hinz, G. M. Kennedy, J. Stone, J. Rigley, S. Ertel, A. Gaspar, V. P. Bailey, W. F. Hoffmann, B. Mennesson, R. Millan-Gabet, W. C. Danchi, O. Absil, P. Arbo, C. Beichman, M. Bonavita, G. Brusa, G. Bryden, E. C. Downey, S. Esposito, P. Grenz, C. Haniff, J. M. Hill, J. M. Leisenring, J. R. Males , et al. (16 additional authors not shown)

Abstract: The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Tel… ▽ More The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47\%$\pm$0.050\% within the central 1.5 au, rising to 0.81\%$\pm$0.026\% within 4.5 au, outside the habitable zone of $β$~Leo. This dust level is 50 $\pm$ 10 times greater than in the solar system's zodiacal cloud. Poynting-Robertson drag on the cold dust detected by Spitzer and Herschel under-predicts the dust present in the habitable zone of $β$~Leo, suggesting an additional delivery mechanism (e.g.,~comets) or an additional belt at $\sim$5.5 au. A model of these dust components is provided which implies the absence of planets more than a few Saturn masses between $\sim$5 au and the outer belt at $\sim$40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8~$μ$m wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5 au to 50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: 11 pages, 9 figures, accepted for publication in Astronomical Journal

Exoplanet Imaging Data Challenge: benchmarking the various image processing methods for exoplanet detection

Authors: F. Cantalloube, C. Gomez-Gonzalez, O. Absil, C. Cantero, R. Bacher, M. J. Bonse, M. Bottom, C. -H. Dahlqvist, C. Desgrange, O. Flasseur, T. Fuhrmann, Th. Henning, R. Jensen-Clem, M. Kenworthy, D. Mawet, D. Mesa, T. Meshkat, D. Mouillet, A. Mueller, E. Nasedkin, B. Pairet, S. Pierard, J. -B. Ruffio, M. Samland, J. Stone , et al. (1 additional authors not shown)

Abstract: The Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is host… ▽ More The Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is hosted on the CodaLab competition platform, where participants can upload their results. The specifications of the data challenge are published on our website. The first phase, launched on the 1st of September 2019 and closed on the 1st of October 2020, consisted in detecting point sources in two types of common data-set in the field of high-contrast imaging: data taken in pupil-tracking mode at one wavelength (subchallenge 1, also referred to as ADI) and multispectral data taken in pupil-tracking mode (subchallenge 2, also referred to as ADI mSDI). In this paper, we describe the approach, organisational lessons-learnt and current limitations of the data challenge, as well as preliminary results of the participants submissions for this first phase. In the future, we plan to provide permanent access to the standard library of data sets and metrics, in order to guide the validation and support the publications of innovative image processing algorithms dedicated to high-contrast imaging of planetary systems. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 36 pages, 27 figures. Proceedings of SPIE "Astronomical Telescopes and Instrumentation" 2020

Journal ref: SPIE Proceedings Volume 11448, Adaptive Optics Systems VII; 114485A (2020)

Search for Tens of MeV Neutrinos associated with Gamma-Ray Bursts in Super-Kamiokande

Authors: The Super-Kamiokande Collaboration, A. Orii, K. Abe, C. Bronner, Y. Hayato, M. Ikeda, S. Imaizumi, H. Ito, J. Kameda, Y. Kataoka, Y. Kato, Y. Kishimoto, M. Miura, S. Moriyama, T. Mochizuki, Y. Nagao, M. Nakahata, Y. Nakajima, S. Nakayama, T. Okada, K. Okamoto, G. Pronost, H. Sekiya, M. Shiozawa, Y. Sonoda , et al. (195 additional authors not shown)

Abstract: A search for neutrinos produced in coincidence with Gamma-Ray Bursts(GRB) was conducted with the Super-Kamiokande (SK) detector. Between December 2008 and March 2017, the Gamma-ray Coordinates Network recorded 2208 GRBs that occurred during normal SK operation. Several time windows around each GRB were used to search for coincident neutrino events. No statistically significant signal in excess of… ▽ More A search for neutrinos produced in coincidence with Gamma-Ray Bursts(GRB) was conducted with the Super-Kamiokande (SK) detector. Between December 2008 and March 2017, the Gamma-ray Coordinates Network recorded 2208 GRBs that occurred during normal SK operation. Several time windows around each GRB were used to search for coincident neutrino events. No statistically significant signal in excess of the estimated backgrounds was detected. The $\barν_e$ fluence in the range from 8 MeV to 100 MeV in positron total energy for $\barν_e+p\rightarrow e^{+}+n$ was found to be less than $\rm 5.07\times10^5$ cm$^{-2}$ per GRB in 90\% C.L. Upper bounds on the fluence as a function of neutrino energy were also obtained. △ Less

Submitted 26 June, 2021; v1 submitted 10 January, 2021; originally announced January 2021.

Comments: 29 pages, 16 figures, 5 tables, accepted by PTEP