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

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

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

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

Guidance of ultraviolet light down to 190 nm in a hollow-core optical fibre

Authors: Robbie Mears, Kerrianne Harrington, William J. Wadsworth, Jonathan C. Knight, James M. Stone, Tim A. Birks

Abstract: We report an anti-resonant hollow core fibre with ultraviolet transmission down to 190 nm, covering the entire UV-A, UV-B and much of the UV-C band. Guidance from 190 - 400 nm is achieved apart for a narrow high loss resonance band at 245 - 265 nm. The minimum attenuation is 0.13 dB/m at 235 nm and 0.16 dB/m at 325 nm. With an inscribed core diameter of ~ 12 $μ$m, the fibre's bend loss at 325 nm w… ▽ More We report an anti-resonant hollow core fibre with ultraviolet transmission down to 190 nm, covering the entire UV-A, UV-B and much of the UV-C band. Guidance from 190 - 400 nm is achieved apart for a narrow high loss resonance band at 245 - 265 nm. The minimum attenuation is 0.13 dB/m at 235 nm and 0.16 dB/m at 325 nm. With an inscribed core diameter of ~ 12 $μ$m, the fibre's bend loss at 325 nm was 0.22 dB per turn for a bend radius of 3 cm at 325 nm. △ Less

Submitted 24 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

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

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)

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

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

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

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

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

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.

Sub millimetre flexible fibre probe for background and fluorescence free Raman spectroscopy

Authors: Stephanos Yerolatsitis, András Kufcsák, Katjana Ehrlich, Harry A. C. Wood, Susan Fernandes, Tom Quinn, Vikki Young, Irene Young, Katie Hamilton, Ahsan R. Akram, Robert R. Thomson, Keith Finlayson, Kevin Dhaliwal, James M. Stone

Abstract: Using the shifted-excitation Raman difference spectroscopy technique and an optical fibre featuring a negative curvature excitation core and a coaxial ring of high numerical aperture collection cores, we have developed a portable, background and fluorescence free, endoscopic Raman probe. The probe consists of a single fibre with a diameter of less than 0.25 mm packaged in a sub-millimetre tubing,… ▽ More Using the shifted-excitation Raman difference spectroscopy technique and an optical fibre featuring a negative curvature excitation core and a coaxial ring of high numerical aperture collection cores, we have developed a portable, background and fluorescence free, endoscopic Raman probe. The probe consists of a single fibre with a diameter of less than 0.25 mm packaged in a sub-millimetre tubing, making it compatible with standard bronchoscopes. The Raman excitation light in the fibre is guided in air and therefore interacts little with silica, enabling an almost background free transmission of the excitation light. In addition, we used the shifted-excitation Raman difference spectroscopy technique and a tunable 785 nm laser to separate the fluorescence and the Raman spectrum from highly fluorescent samples, demonstrating the suitability of the probe for biomedical applications. Using this probe we also acquired fluorescence free human lung tissue data. △ Less

Submitted 16 December, 2020; originally announced December 2020.

High Contrast Thermal Infrared Spectroscopy with ALES: The 3-4$μ$m Spectrum of $κ$ Andromedae b

Authors: Jordan M. Stone, Travis Barman, Andrew J. Skemer, Zackery W. Briesemeister, Laci S. Brock, Philip M. Hinz, Jarron M. Leisenring, Charles E. Woodward, Michael F. Skrutskie, Eckhart Spalding

Abstract: We present the first $L-$band (2.8 to 4.1~$μ$m) spectroscopy of $κ$~Andromedae~b, a $\sim20~M_{\mathrm{Jup}}$ companion orbiting at $1^{\prime\prime}$ projected separation from its B9-type stellar host. We combine our Large Binocular Telescope ALES integral field spectrograph data with measurements from other instruments to analyze the atmosphere and physical characteristics of $κ$~And~b. We repor… ▽ More We present the first $L-$band (2.8 to 4.1~$μ$m) spectroscopy of $κ$~Andromedae~b, a $\sim20~M_{\mathrm{Jup}}$ companion orbiting at $1^{\prime\prime}$ projected separation from its B9-type stellar host. We combine our Large Binocular Telescope ALES integral field spectrograph data with measurements from other instruments to analyze the atmosphere and physical characteristics of $κ$~And~b. We report a discrepancy of $\sim20\%$ ($2σ$) in the $L^{\prime}$ flux of $κ$~And~b when comparing to previously published values. We add an additional $L^{\prime}$ constraint using an unpublished imaging dataset collected in 2013 using LBTI/LMIRCam, the instrument in which the ALES module has been built. The LMIRCam measurement is consistent with the ALES measurement, both suggesting a fainter $L$-band scaling than previous studies. The data, assuming the flux scaling measured by ALES and LMIRCam imaging, are well fit by an L3-type brown dwarf. Atmospheric model fits to measurements spanning 0.9-4.8~$μ$m reveal some tension with the predictions of evolutionary models, but the proper choice of cloud parameters can provide some relief. In particular, models with clouds extending to very-low pressures composed of grains $\leq1~μ$m appear to be necessary. If the brighter $L^{\prime}$ photometry is accurate, there is a hint that sub-solar metallicity may be required. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted for publication in AJ

Stratified Global MHD Models of Accretion Disks in Semi-Detached Binaries

Authors: Patryk Pjanka, James M. Stone

Abstract: We present results of the first global magnetohydrodynamic (MHD) simulations of accretion disks fed by Roche lobe overflow, including vertical stratification, in order to investigate the roles of spiral shocks, magnetorotational instability (MRI), and the accretion stream on disk structure and evolution. Our models include a simple treatment of gas thermodynamics, with orbital Mach numbers at the… ▽ More We present results of the first global magnetohydrodynamic (MHD) simulations of accretion disks fed by Roche lobe overflow, including vertical stratification, in order to investigate the roles of spiral shocks, magnetorotational instability (MRI), and the accretion stream on disk structure and evolution. Our models include a simple treatment of gas thermodynamics, with orbital Mach numbers at the inner edge of the disk $M_{\rm in}$ of 5 and 10. We find mass accretion rates to vary considerably on all time scales, with only the Mach 5 model reaching a clear quasi-stationary state. For Mach 10, the model undergoes an outside-in magnetically-driven accretion event occurring on a time scale of $\sim10$ orbital periods of the binary. Both models exhibit spiral shocks inclined with respect to the binary plane, with their position and inclination changing rapidly. However, the time-averaged location of these shocks in the equatorial plane is well-fit by simple linear models. MRI turbulence in the disk generates toroidal magnetic field patterns (butterfly diagrams) that are in some cases irregular, perhaps due to interaction with spiral structure. While many of our results are in good agreement with local studies, we find some features (most notably those related to spiral shocks) can only be captured in global models such as studied here. Thus, while global studies remain computationally expensive -- even as idealized models -- they are essential (along with more sophisticated treatment of radiation transport and disk thermodynamics) for furthering our understanding of accretion in binary systems. △ Less

Submitted 3 October, 2020; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 30 pages, 13 figures; Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 904:90 (2020)

Turbulent dissipation, CH$^+$ abundance, H$_2$ line luminosities, and polarization in the cold neutral medium

Authors: Eric R. Moseley, B. T. Draine, Kengo Tomida, James M. Stone

Abstract: In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We… ▽ More In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. We make two important simplifying assumptions: (i) the ${\rm H}_2/{\rm H}$ fraction is fixed everywhere, and (ii) we exclude from our analysis regions where the ion-neutral drift velocity is calculated to be greater than 5 km/s. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $μ$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. These findings contrast with those of Valdivia et al. (2017). Comparison of predicted dust polarization with observations by {\it Planck} suggests that the mean field $\gtrsim 5 μ$G, so that the turbulence is sub-Alfvénic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance. △ Less

Submitted 23 October, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 20 pages, 18 figures, accepted for publication to MNRAS

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

Authors: Sean M. Ressler, Christopher J. White, Eliot Quataert, James M. Stone

Abstract: We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained pro… ▽ More We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation. △ Less

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

Comments: Accepted by ApJL. 14 pages, 8 figures. Animations available at https://smressle.bitbucket.io/animations.html

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

Authors: James M. Stone, Kengo Tomida, Christopher J. White, Kyle G. Felker

Abstract: The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic exe… ▽ More The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: 50 pages, 41 figures, accepted for publication in American Astronomical Society journals

The HOSTS survey for exozodiacal dust: Observational results from the complete survey

Authors: Steve Ertel, Denis Defrère, Philip M. Hinz, Bertrand Mennesson, Grant M. Kennedy, William C. Danchi, Christopher Gelino, John M. Hill, William F. Hoffmann, Johan Mazoyer, George Rieke, Andrew Shannon, Karl Stapelfeldt, Eckhart Spalding, Jordan M. Stone, Amali Vaz, Alycia J. Weinberger, Phil Willems, Olivier Absil, Paul Arbo, Vanessa P. Bailey, Charles Beichman, Geoffrey Bryden, Elwood C. Downey, Olivier Durney , et al. (21 additional authors not shown)

Abstract: The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model p… ▽ More The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates. △ Less

Submitted 6 March, 2020; originally announced March 2020.

Comments: accepted for publication in AJ

The Surprisingly Small Impact of Magnetic Fields On The Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds

Authors: Sean M. Ressler, Eliot Quataert, James M. Stone

Abstract: We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic… ▽ More We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $β$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $β$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: Accepted by MNRAS; 22 pages; 24 figures, comments welcome

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

Authors: Zhaohuan Zhu, Yan-Fei Jiang, James M. Stone

Abstract: FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion… ▽ More FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3. △ Less

Submitted 5 March, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 22 pages, 20 figures, accepted by MNRAS

Thermal Infrared Imaging of MWC 758 with the Large Binocular Telescope: Planetary Driven Spiral Arms?

Authors: Kevin Wagner, Jordan M. Stone, Eckhart Spalding, Daniel Apai, Ruobing Dong, Steve Ertel, Jarron Leisenring, Ryan Webster

Abstract: Theoretical studies suggest that a giant planet around the young star MWC 758 could be responsible for driving the spiral features in its circumstellar disk. Here, we present a deep imaging campaign with the Large Binocular Telescope with the primary goal of imaging the predicted planet. We present images of the disk in two epochs in the $L^{\prime}$ filter (3.8 $μm$) and a third epoch in the… ▽ More Theoretical studies suggest that a giant planet around the young star MWC 758 could be responsible for driving the spiral features in its circumstellar disk. Here, we present a deep imaging campaign with the Large Binocular Telescope with the primary goal of imaging the predicted planet. We present images of the disk in two epochs in the $L^{\prime}$ filter (3.8 $μm$) and a third epoch in the $M^{\prime}$ filter (4.8 $μm$). The two prominent spiral arms are detected in each observation, which constitute the first images of the disk at $M^\prime$, and the deepest yet in $L^\prime$ ($ΔL^\prime=$12.1 exterior to the disk at 5$σ$ significance). We report the detection of a S/N$\sim$3.9 source near the end of the Sourthern arm, and, from the source's detection at a consistent position and brightness during multiple epochs, we establish a $\sim$90% confidence-level that the source is of astrophysical origin. We discuss the possibilities that this feature may be a) an unresolved disk feature, and b) a giant planet responsible for the spiral arms, with several arguments pointing in favor of the latter scenario. We present additional detection limits on companions exterior to the spiral arms, which suggest that a $\lesssim$4 M$_{Jup}$ planet exterior to the spiral arms could have escaped detection. Finally, we do not detect the companion candidate interior to the spiral arms reported recently by Reggiani et al. (2018), although forward modelling suggests that such a source would have likely been detected. △ Less

Submitted 12 August, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 10 pages, 4 figures, accepted for publication in ApJ

Bondi-Hoyle-Lyttleton accretion in Supergiant X-ray binaries: stability and disk formation

Authors: Wenrui Xu, James M. Stone

Abstract: We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi-Hoyle-Lyttleton (BHL) accretion with and without transverse upstream gradients. We mainly focus on the regime of high (upstream) Mach number, weak upstream gradients and small accretor size, which is relevant to neutron star (NS) accretion in wind-fed Supergiant X-ray binaries (SgXBs). We present a systematic exploration of th… ▽ More We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi-Hoyle-Lyttleton (BHL) accretion with and without transverse upstream gradients. We mainly focus on the regime of high (upstream) Mach number, weak upstream gradients and small accretor size, which is relevant to neutron star (NS) accretion in wind-fed Supergiant X-ray binaries (SgXBs). We present a systematic exploration of the flow in this regime. When there are no upstream gradients, the flow is always stable regardless of accretor size or Mach number. For finite upstream gradients, there are three main types of behavior: stable flow (small upstream gradient), turbulent unstable flow without a disk (intermediate upstream gradient), and turbulent flow with a disk-like structure (relatively large upstream gradient). When the accretion flow is turbulent, the accretion rate decreases non-convergently as the accretor size decreases. The flow is more prone to instability and the disk is less likely to form than previously expected; the parameters of most observed SgXBs place them in the regime of a turbulent, disk-less accretion flow. Among the SgXBs with relatively well-determined parameters, we find OAO 1657-415 to be the only one that is likely to host a persistent disk (or disk-like structure); this finding is consistent with observations. △ Less

Submitted 13 July, 2019; originally announced July 2019.

Comments: 25 pages, 28 figures, 2 tables. Accepted by MNRAS

The X-ray Halo Scaling Relations of Supermassive Black Holes

Authors: M. Gaspari, D. Eckert, S. Ettori, P. Tozzi, L. Bassini, E. Rasia, F. Brighenti, M. Sun, S. Borgani, S. D. Johnson, G. R. Tremblay, J. M. Stone, P. Temi, H. -Y. K. Yang, F. Tombesi, M. Cappi

Abstract: We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by… ▽ More We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-σ_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions). △ Less

Submitted 22 October, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: 40 pages, 27 figures, 3 tables; ApJ accepted version (minor revision) - Reviewer: "The authors are to be congratulated on a very interesting paper that deserves to be read by anyone interested in the connection between the SMBH and host galaxy."

Journal ref: ApJ, 884, 169 (2019)

Non-ideal MHD simulation of HL Tau disk: formation of rings

Authors: Xiao Hu, Zhaohuan Zhu, Satoshi Okuzumi, Xue-Ning Bai, Lile Wang, Kengo Tomida, James M. Stone

Abstract: Recent high resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D non-ideal MHD simulations including effects from both Ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The non-ideal MHD diffusion profiles are calculated based on the… ▽ More Recent high resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D non-ideal MHD simulations including effects from both Ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The non-ideal MHD diffusion profiles are calculated based on the global dust evolution calculation including sintering effects. Disk ionization structure changes dramatically across the snow line due to the change of dust size distribution close to snow line of major volatiles. We find that accretion is mainly driven by disk wind. Gaps and rings can be quickly produced from different accretion rates across snow line. Furthermore, ambipolar diffusion (AD) leads to highly preferential accretion at midplane, followed by magnetic reconnection. This results a local zone of decretion that drains of mass in the field reconnection area, which leaves a gap and an adjacent ring just outside it. Overall, under the favorable condition, both snow lines and non-ideal MHD effects can lead to gaseous gaps and rings in protoplanetary disks. △ Less

Submitted 30 October, 2019; v1 submitted 18 April, 2019; originally announced April 2019.

Comments: 13 pages, 10 figures, published on ApJ

High Spatial Resolution Thermal-Infrared Spectroscopy with ALES: Resolved Spectra of the Benchmark Brown Dwarf Binary HD 130948BC

Authors: Zackery W. Briesemeister, Andrew J. Skemer, Jordan M. Stone, Travis S. Barman, Philip Hinz, Jarron Leisenring, Michael F. Skrutskie, Charles E. Woodward, Eckhart Spalding

Abstract: We present 2.9-4.1 micron integral field spectroscopy of the L4+L4 brown dwarf binary HD 130948BC, obtained with the Arizona Lenslets for Exoplanet Spectroscopy (ALES) mode of the Large Binocular Telescope Interferometer (LBTI). The HD 130948 system is a hierarchical triple system, in which the G2V primary is joined by two co-orbiting brown dwarfs. By combining the age of the system with the dynam… ▽ More We present 2.9-4.1 micron integral field spectroscopy of the L4+L4 brown dwarf binary HD 130948BC, obtained with the Arizona Lenslets for Exoplanet Spectroscopy (ALES) mode of the Large Binocular Telescope Interferometer (LBTI). The HD 130948 system is a hierarchical triple system, in which the G2V primary is joined by two co-orbiting brown dwarfs. By combining the age of the system with the dynamical masses and luminosities of the substellar companions, we can test evolutionary models of cool brown dwarfs and extra-solar giant planets. Previous near-infrared studies suggest a disagreement between HD 130948BC luminosities and those derived from evolutionary models. We obtained spatially-resolved, low-resolution (R~20) L-band spectra of HD 130948B and C to extend the wavelength coverage into the thermal infrared. Jointly using JHK photometry and ALES L-band spectra for HD 130948BC, we derive atmospheric parameters that are consistent with parameters derived from evolutionary models. We leverage the consistency of these atmospheric quantities to favor a younger age (0.50 \pm 0.07 Gyr) of the system compared to the older age (0.79 \pm 0.22 Gyr) determined with gyrochronology in order to address the luminosity discrepancy. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 17 pages, 9 figures, Accepted to ApJ

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

[Plasma 2020 Decadal] The Material Properties of Weakly Collisional, High-Beta Plasmas

Authors: M. W. Kunz, J. Squire, S. A. Balbus, S. D. Bale, C. H. K. Chen, E. Churazov, S. C. Cowley, C. B. Forest, C. F. Gammie, E. Quataert, C. S. Reynolds, A. A. Schekochihin, L. Sironi, A. Spitkovsky, J. M. Stone, I. Zhuravleva, E. G. Zweibel

Abstract: This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood:… ▽ More This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

A Resolution Study of Magnetically Arrested Disks

Authors: Christopher J. White, James M. Stone, Eliot Quataert

Abstract: We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree ac… ▽ More We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford-Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: 19 pages, 18 figures, to be published in ApJ

Journal ref: ApJ 874 (2019) 168

Hydrodynamic Torques in Circumbinary Accretion Disks

Authors: Mackenzie S. L. Moody, Ji-Ming Shi, James M. Stone

Abstract: Gaseous disks have been proposed as a mechanism for facilitating mergers of binary black holes. We explore circumbinary disk systems to determine the evolution of the central binary. To do so, we perform 3D, hydrodynamic, locally isothermal simulations of circumbinary disks on a Cartesian grid. We focus on binaries of equal mass ratios on fixed circular orbits. To investigate the orbital evolution… ▽ More Gaseous disks have been proposed as a mechanism for facilitating mergers of binary black holes. We explore circumbinary disk systems to determine the evolution of the central binary. To do so, we perform 3D, hydrodynamic, locally isothermal simulations of circumbinary disks on a Cartesian grid. We focus on binaries of equal mass ratios on fixed circular orbits. To investigate the orbital evolution of the binary, we examine the various torques exerted on the system. For the case where the disk plane and binary orbital plane are aligned, we find that the total torque is positive so that the semi-major axis of the binary increases. For the misaligned case, we run simulations with the binary orbital plane and disk midplane misaligned by 45\degree and find the same results - the binary grows. The timescale for the circumbinary disk to realign to the plane of the binary is consistent with the global viscous timescale of the disk. △ Less

Submitted 28 February, 2019; originally announced March 2019.

Comments: Accepted to ApJ. 13 pages, 19 figures

Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution

Authors: Xue-Ning Bai, Eve C. Ostriker, Illya Plotnikov, James M. Stone

Abstract: The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfvén waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($λ\ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$)… ▽ More The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfvén waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($λ\ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $δf$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfvén modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 34 pages, 28 figures, submitted to ApJ

Polygram Stars: Resonant Tidal Excitation of Fundamental Oscillation Modes in Asynchronous Stellar Coalescence

Authors: Morgan MacLeod, Michelle Vick, Dong Lai, James M. Stone

Abstract: The prevalence of binary stars at close separations implies that many of these systems will interact or merge during the binary's lifetime. This paper presents hydrodynamic simulations of the scenario of binary coalescence through unstable mass transfer, which drives the pair to closer separations. When the donor star does not rotate synchronously with respect to the orbit, dynamical tidal waves a… ▽ More The prevalence of binary stars at close separations implies that many of these systems will interact or merge during the binary's lifetime. This paper presents hydrodynamic simulations of the scenario of binary coalescence through unstable mass transfer, which drives the pair to closer separations. When the donor star does not rotate synchronously with respect to the orbit, dynamical tidal waves are excited in its envelope. We show that resonance crossings with high azimuthal-order $(m\sim3$ to $6$) fundamental modes induce a visible "polygram" distortion to the star. As the binary orbit tightens, the system sweeps through resonance with modes of decreasing azimuthal order, which are selectively excited. We compare our hydrodynamic simulations to predictions from linear theory of resonant mode excitation. The linear theory provides an estimate of mode amplitudes to within a factor of two, even as the oscillations become quite non-linear as the stars coalesce. We estimate that a wave with 10% radial amplitude generates approximately 1% photometric variability; this may be detectible if such a binary coalescence is caught in action by future photometric all-sky surveys. △ Less

Submitted 18 December, 2018; originally announced December 2018.

Comments: Submitted to AAS Journals. Comments for revision are welcome

The LEECH Exoplanet Imaging Survey: Limits on Planet Occurrence Rates Under Conservative Assumptions

Authors: Jordan M. Stone, Andrew J. Skemer, Philip M. Hinz, Mariangela Bonavita, Kaitlin M. Kratter, Anne-Lise Maire, Denis Defrere, Vanessa P. Bailey, Eckhart Spalding, Jarron M. Leisenring, S. Desidera, M. Bonnefoy, Beth Biller, Charles E. Woodward, Th. Henning, Michael F. Skrutskie, J. A. Eisner, Justin R. Crepp, Jennifer Patience, Gerd Weigelt, Robert J. De Rosa, Joshua Schlieder, Wolfgang Brandner, Dániel Apai, Kate Su , et al. (11 additional authors not shown)

Abstract: We present the results of the largest $L^{\prime}$ ($3.8~μ$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass,… ▽ More We present the results of the largest $L^{\prime}$ ($3.8~μ$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass, old, and cold-start giant planets. We emphasize proximity over youth in our target selection, probing physical separations smaller than other direct imaging surveys. For FGK stars, LEECH outperforms many previous studies, placing tighter constraints on the hot-start planet occurrence frequency interior to $\sim20$ au. For less luminous, cold-start planets, LEECH provides the best constraints on giant-planet frequency interior to $\sim20$ au around FGK stars. Direct imaging survey results depend sensitively on both the choice of evolutionary model (e.g., hot- or cold-start) and assumptions (explicit or implicit) about the shape of the underlying planet distribution, in particular its radial extent. Artificially low limits on the planet occurrence frequency can be derived when the shape of the planet distribution is assumed to extend to very large separations, well beyond typical protoplanetary dust-disk radii ($\lesssim50$ au), and when hot-start models are used exclusively. We place a conservative upper limit on the planet occurrence frequency using cold-start models and planetary population distributions that do not extend beyond typical protoplanetary dust-disk radii. We find that $\lesssim90\%$ of FGK systems can host a 7 to 10 $M_{\mathrm{Jup}}$ planet from 5 to 50 au. This limit leaves open the possibility that planets in this range are common. △ Less

Submitted 6 December, 2018; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 31 pages, 13 figures, accepted to AJ

Accretion of Magnetized Stellar Winds in the Galactic Centre: Implications for Sgr A* and PSR J1745-2900

Authors: Sean M. Ressler, Eliot Quataert, James M. Stone

Abstract: The observed rotation measures (RMs) towards the galactic centre magnetar and towards Sagittarius A* provide a strong constraint on MHD models of the galactic centre accretion flow, probing distances from the black hole separated by many orders of magnitude. We show, using 3D simulations of accretion via magnetized stellar winds of the Wolf-Rayet stars orbiting the black hole, that the large, time… ▽ More The observed rotation measures (RMs) towards the galactic centre magnetar and towards Sagittarius A* provide a strong constraint on MHD models of the galactic centre accretion flow, probing distances from the black hole separated by many orders of magnitude. We show, using 3D simulations of accretion via magnetized stellar winds of the Wolf-Rayet stars orbiting the black hole, that the large, time-variable RM observed for the pulsar PSR J1745-2900 can be explained by magnetized wind-wind shocks of nearby stars in the clockwise stellar disc. In the same simulation, both the total X-ray luminosity integrated over 2-10$''$, the time variability of the magnetar's dispersion measure, and the RM towards Sagittarius A* are consistent with observations. We argue that (in order for the large RM of the pulsar to not be a priori unlikely) the pulsar should be on an orbit that keeps it near the clockwise disc of stars. We present a 2D RM map of the central 1/2 parsec of the galactic centre that can be used to test our models. Our simulations predict that Sgr A* is typically accreting a significantly ordered magnetic field that ultimately could result in a strongly magnetized flow with flux threading the horizon at $\sim$ 10$\%$ of the magnetically arrested limit. △ Less

Submitted 22 November, 2018; v1 submitted 19 October, 2018; originally announced October 2018.

Comments: Accepted into MNRAS Letters. Comments welcome