-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
[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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Methane in Analogs of Young Directly Imaged Exoplanets
Authors:
Brittany E. Miles,
Andrew J. Skemer,
Travis S. Barman,
Katelyn N. Allers,
Jordan M. Stone
Abstract:
We present detections of methane in R of $\sim$1300, L band spectra of VHS 1256 b and PSO 318.5, two low gravity, red, late L dwarfs that share the same colors as the HR 8799 planets. These spectra reveal shallow methane features, which indicate VHS 1256 b and PSO 318.5 have photospheres that are out of chemical equilibrium. Directly imaged exoplanets usually have redder near infrared colors than…
▽ More
We present detections of methane in R of $\sim$1300, L band spectra of VHS 1256 b and PSO 318.5, two low gravity, red, late L dwarfs that share the same colors as the HR 8799 planets. These spectra reveal shallow methane features, which indicate VHS 1256 b and PSO 318.5 have photospheres that are out of chemical equilibrium. Directly imaged exoplanets usually have redder near infrared colors than the field-age population of brown dwarfs on a color magnitude diagram. These objects along the L to T transition show reduced methane absorption and evidence of photospheric clouds. Compared to the H and K bands, L band (3 micron - 4 micron) spectroscopy provides stronger constraints on the methane abundances of brown dwarfs and directly imaged exoplanets that have similar effective temperatures as L to T transition objects. When combined with near infrared spectra, the L band extends our conventional wavelength coverage, increasing our understanding of atmospheric cloud structure. Our model comparisons show relatively strong vertical mixing and photospheric clouds can explain the molecular absorption features and continua of VHS 1256 b and PSO 318.5. We also discuss the implications of this work for future exoplanet focused instruments and observations with the James Webb Space Telescope.
△ Less
Submitted 21 November, 2018; v1 submitted 10 October, 2018;
originally announced October 2018.
-
Bound Outflows, Unbound Ejecta, and the Shaping of Bipolar Remnants during Stellar Coalescence
Authors:
Morgan MacLeod,
Eve C. Ostriker,
James M. Stone
Abstract:
Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe ove…
▽ More
Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae.
△ Less
Submitted 16 November, 2018; v1 submitted 17 August, 2018;
originally announced August 2018.