Researcher profile

Christopher J. White

Christopher J. White contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Blacklight: A General-Relativistic Ray-Tracing and Analysis Tool

We describe the Blacklight code, intended for post-processing general-relativistic magnetohydrodynamic simulation data. Beyond polarized ray tracing of synchrotron radiation, it can produce a number of outputs that aid in analyzing data sets, such as maps of auxiliary quantities and false-color renderings. Additional features include support for adaptive mesh refinement input, slow-light calculations, and adaptive ray tracing. The code is written with ease of use, readability, and transparency as primary objectives, while it still achieves high performance. Blacklight is publicly available and released into the public domain.

preprint2022arXiv

Observational Signatures of Black Hole Accretion: Rotating vs. Spherical Flows with Tilted Magnetic Fields

We study the observational signatures of magnetically arrested black hole accretion with non-rotating inflow onto a rotating black hole; we consider a range of angles between the black hole spin and the initial magnetic field orientation. We compare the results of our General Relativistic Magneto-Hydrodynamic simulations to more commonly used rotating initial conditions and to the Event Horizon Telescope (EHT) observations of M87. We find that the mm intensity images, polarization images, and synchrotron emission spectra are very similar among the different simulations when post-processed with the same electron temperature model; observational differences due to different electron temperature models are significantly larger than those due to the different realizations of magnetically arrested accretion. The orientation of the mm synchrotron polarization is particularly insensitive to the initial magnetic field orientation, the electron temperature model, and the rotation of the inflowing plasma. The largest difference among the simulations with different initial rotation and magnetic tilt is in the strength and stability of the jet; spherical inflow leads to kink-unstable jets. We discuss the implications of our results for current and future EHT observations and for theoretical models of event-horizon-scale black hole accretion.

preprint2022arXiv

Remarkable correspondence of Sagittarius A* submillimeter variability with a stellar-wind-fed accretion flow model

We compare the 230 GHz near-horizon emission from Sagittarius A* to simulations representing three classes of accretion flows. Using the structure function to capture the variability statistics of the light curve, we find a noticeable discrepancy between the observations and models based on torus-fed accretion disks, whether those disks bring in a small or large amount of net magnetic flux. On the other hand, the simulations that are fed more realistically by stellar winds match the observed structure function very well. We describe the differences between models, arguing that feeding by stellar winds may be a critical component in constructing theoretical models for accretion in the Galactic Center.

preprint2022arXiv

The Physical Effects of Progenitor Rotation: Comparing Two Long-Duration 3D Core-Collapse Supernova Simulations

We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despite this difference, the asymptotic explosion energies ($\sim$10$^{50}$ ergs) and residual neutron star baryon masses ($\sim$1.3 M$_{\odot}$) are similar. We find that the proto-neutron star (PNS) core can deleptonize and cool significantly more quickly. Soon into the evolution of the rotating model, we witness more vigorous and extended PNS core convection that early in its evolution envelopes the entire inner sphere, not just a shell. Moreover, we see a corresponding excursion in both the $ν_e$ luminosity and gravitational-wave strain that may be diagnostic of this observed dramatic phenomenon. In addition, after bounce the innermost region of the rotating model seems to execute meridional circulation. The rotationally-induced growth of the convective PNS region may facilitate the growth of core B-fields by the dynamo mechanism by facilitating the achievement of the critical Rossby number condition for substantial growth of a dipole field, obviating the need for rapid rotation rates to create dipole fields of significance. The next step is to explore the progenitor-mass and spin dependencies across the progenitor continuum of the supernova explosion, dynamics, and evolution of PNS convection and its potential role in the generation of magnetar and pulsar magnetic fields.

preprint2021arXiv

Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87.

preprint2021arXiv

On the Origin of Pulsar and Magnetar Magnetic Fields

In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective proto-neutron stars at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 solar masses to 25 solar masses, including one with initial rotation, have sufficiently vigorous proto-neutron-star convection to generate dipole fields of order ~10^15 gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.

preprint2021arXiv

The Effects of Tilt on the Time Variability of Millimeter and Infrared Emission from Sagittarius A*

Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24 degrees) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A*), though our results apply more generally to optically thin flows. For most quantities, tilted accretion is more variable, primarily due to a significantly hotter and denser "coronal" region well off the disk midplane. We find (1) a 150% increase in millimeter light curve variability when adding tilt to the flow; (2) the tilted image centroid in the millimeter shifts on a scale of 3.7 microarcseconds over 28 hours (5000 gravitational times) for some electron temperature models; (3) tilted disk image diameters in the millimeter can be 10% larger (52 versus 47 microarcseconds) than those of aligned disks at certain viewing angles; (4) the tilted models produce significant IR flux, similar to that seen in Sgr A*, with comparable or even greater variability than observed; (5) for some electron models, the tilted IR centroid moves by more than 50 microarcseconds over several hours, in a similar fashion to the centroid motion detected by the GRAVITY interferometer.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The Effects of Resolution on Black Hole Accretion Simulations of Jets

We perform three general-relativistic magnetohydrodynamic simulations of black hole accretion designed to test how sensitive results are to grid resolution in the jet region. The cases differ only in numerics, modelling the same physical scenario of a radiatively inefficient, geometrically thick, magnetically arrested flow onto a rapidly spinning black hole. Properties inferred with the coarsest grid generally agree with those found with higher resolutions, including total jet power and its decomposition into different forms, velocity structure, nonaxisymmetric structure, and the appearance of resolved millimetre images. Some measures of variability and magnetization are sensitive to resolution. We conclude that most results obtained by limiting resolution near the jet for computational expediency should still be reliable, at least insofar as they would not be improved with a finer grid.

preprint2020arXiv

The Effects of Tilt on the Images of Black Hole Accretion Flows

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

preprint2020arXiv

The Structure of Radiatively Inefficient Black Hole Accretion Flows

We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting.

preprint2019arXiv

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motivated by the recent discovery of fast-moving objects that occupy unusual locations on the HR diagram, we explore the impact of the thermonuclear supernovae on the donors in this specific double detonation scenario. We use MESA to model the binary up to the moment of detonation, then 3D Athena++ to model the hydrodynamic interaction of the supernova ejecta with the donor star, calculating the amount of mass that is stripped and the entropy deposited in the deep stellar interior by the strong shock that traverses it. We show that these donor remnants are ejected with velocities primarily set by their orbital speeds: $700-900\ \rm km\ s^{-1}$. We model the long-term thermal evolution of remnants by introducing the shock entropy into MESA models. In response to this entropy change, donor remnants expand and brighten for timescales ranging from $10^6-10^8$ years, giving ample time for these runaway stars to be observed in their inflated state before they leave the galaxy. Even after surface layers are stripped, some donors retain enough mass to resume core helium burning and further delay fading for more than $10^8$ years.