BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Joshua C. Dolence

Joshua C. Dolence appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRgr-qcastro-ph.GAastro-ph.IMmath.NANumerical AnalysisPerformancephysics.comp-ph

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

Not-Quite Transcendental Functions and their Applications

Transcendental functions, such as exponentials and logarithms, appear in a broad array of computational domains: from simulations in curvilinear coordinates, to interpolation, to machine learning. Unfortunately they are typically expensive to compute accurately. In this note, we argue that in many cases, the properties of the function matters more than the exact functional form. We present new functions, which are not transcendental, that can be used as drop-in replacements for the exponential and logarithm in many settings for a significant performance boost. We show that for certain applications using these functions result in no drop in the accuracy at all, as they are perfectly accurate representations of themselves, if not the original transcendental functions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Convection-Aided Explosions in One-Dimensional Core-Collapse Supernova Simulations I: Technique and Validation

Most one-dimensional core-collapse simulations fail to explode, yet multi-dimensional simulations often explode. A dominant multi-dimensional effect aiding explosion is neutrino-driven convection. We incorporate a convection model in approximate one-dimensional core-collapse supernova (CCSN) simulations. This is the 1D+ method. This convection model lowers the neutrino luminosity required for explosion by 30%, similar to the reduction observed in multi-dimensional simulations. The model is based upon the global turbulence model of Mabanta & Murphy (2018) and models the mean-field turbulent flow of neutrino-driven convection. In this preliminary investigation, we use simple neutrino heating and cooling algorithms to compare the critical condition in the 1D+ simulations with the critical condition observed in two-dimensional simulations. Qualitatively, the critical conditions in the 1D+ and the two-dimensional simulations are similar. The assumptions in the convection model affect the radial profiles of density, entropy, and temperature, and comparisons with the profiles of three dimensional simulations will help to calibrate these assumptions. These 1D+ simulations are consistent with the profiles and explosion conditions of equivalent two-dimensional CCSN simulations but are ~100 times faster, and the 1D+ prescription has the potential to be ~100,000 faster than three-dimensional CCSN simulations. The 1D+ technique will be ideally suited to test the explodability of thousands of progenitor models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

We present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress-energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Radiation four-fources are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way, with Compton scattering further approximated by the Kompaneets equation. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that MOCMC converges at least as $\sim N^{-1}$, rather than the canonical Monte Carlo $N^{-1/2}$, where $N$ is the number of samples per zone. On several problems we demonstrate substantial improvement over Eddington and M1 closures and gray opacities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting the Supernova Breakout Burst in Terrestrial Neutrino Detectors

We calculate the distance-dependent performance of a few representative terrestrial neutrino detectors in detecting and measuring the properties of the $ν_e$ breakout burst light curve in a Galactic core-collapse supernova. The breakout burst is a signature phenomenon of core collapse and offers a probe into the stellar core through collapse and bounce. We examine cases of no neutrino oscillations and oscillations due to normal and inverted neutrino-mass hierarchies. For the normal hierarchy, other neutrino flavors emitted by the supernova overwhelm the $ν_e$ signal, making a detection of the breakout burst difficult. For the inverted hierarchy (IH), some detectors at some distances should be able to see the $ν_e$ breakout burst peak and measure its properties. For the IH, the maximum luminosity of the breakout burst can be measured at 10 kpc to accuracies of $\sim$30% for Hyper-Kamiokande (Hyper-K) and $\sim$60% for the Deep Underground Neutrino Experiment (DUNE). Super-Kamiokande (Super-K) and Jiangmen Underground Neutrino Observatory (JUNO) lack the mass needed to make an accurate measurement. For the IH, the time of the maximum luminosity of the breakout burst can be measured in Hyper-K to an accuracy of $\sim$3 ms at 7 kpc, in DUNE $\sim$2 ms at 4 kpc, and JUNO and Super-K can measure the time of maximum luminosity to an accuracy of $\sim$2 ms at 1 kpc. Detector backgrounds in IceCube render a measurement of the $ν_e$ breakout burst unlikely. For the inverted hierarchy, a measurement of the maximum luminosity of the breakout burst could be used to differentiate between nuclear equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

bhlight: General Relativistic Radiation Magnetohydrodynamics with Monte Carlo Transport

We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics (GRRMHD) using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-Dimensional Core-Collapse Supernova Models with Multi-Dimensional Transport

We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino transport, including all relevant $\mathcal{O}(v/c)$ terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate $\mathcal{O}(v/c)$ terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 milliseconds after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ``ray-by-ray' approach employed by all other groups may be compromising their results. We show that ``ray-by-ray' calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Dominance of Neutrino-Driven Convection in Core-Collapse Supernovae

Multi-dimensional instabilities have become an important ingredient in core-collapse supernova (CCSN) theory. Therefore, it is necessary to understand the driving mechanism of the dominant instability. We compare our parameterized three-dimensional CCSN simulations with other buoyancy-driven simulations and propose scaling relations for neutrino-driven convection. Through these comparisons, we infer that buoyancy-driven convection dominates post-shock turbulence in our simulations. In support of this inference, we present four major results. First, the convective fluxes and kinetic energies in the neutrino-heated region are consistent with expectations of buoyancy-driven convection. Second, the convective flux is positive where buoyancy actively drives convection, and the radial and tangential components of the kinetic energy are in rough equipartition (i.e. K_r ~ K_θ + K_ϕ). Both results are natural consequences of buoyancy-driven convection, and are commonly observed in simulations of convection. Third, buoyant driving is balanced by turbulent dissipation. Fourth, the convective luminosity and turbulent dissipation scale with the driving neutrino power. In all, these four results suggest that in neutrino-driven explosions, the multi-dimensional motions are consistent with neutrino-driven convection.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Monte Carlo Method for Time-Dependent Neutrino Radiation Transport

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Investigation into the Character of Pre-Explosion Core-Collapse Supernova Shock Motion

We investigate the structure of the stalled supernova shock in both 2D and 3D and explore the differences in the effects of neutrino heating and the standing accretion shock instability (SASI). We find that early on the amplitude of the dipolar mode of the shock is factors of 2 to 3 smaller in 3D than in 2D. However, later in both 3D and 2D the monopole and dipole modes start to grow until explosion. Whereas in 2D the (l,m) = (1,0) mode changes sign quasi-periodically, producing the "up-and-down" motion always seen in modern 2D simulations, in 3D this almost never happens. Rather, in 3D when the dipolar mode starts to grow, it grows in magnitude and wanders stochastically in direction until settling before explosion to a particular patch of solid angle. In 2D we find that the amplitude of the dipolar shock deformation separates into two classes. For the first, identified with the SASI and for a wide range of "low" neutrino luminosities, this amplitude remains small and roughly constant. For the other, identified with higher luminosities and neutrino-driven convection, the dipolar amplitude grows sharply. Importantly, it is only for this higher luminosity class that we see neutrino-driven explosions within ~1 second of bounce. Moreover, for the "low" luminosity runs, the power spectra of these dipolar oscillations peak in the 30-50 Hz range associated with advection timescales, while for the high-luminosity runs the power spectra at lower frequencies are significantly more prominent. We associate this enhanced power at lower frequencies with slower convective effects and the secular growth of the dipolar shock amplitude. On the basis of our study, we hypothesize that neutrino-driven buoyant convection should almost always dominate the SASI when the supernova explosion is neutrino-driven.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dimensional Dependence of the Hydrodynamics of Core-Collapse Supernovae

The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two- and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Near-Infrared and X-ray Quasi-Periodic Oscillations in Numerical Models of Sgr A*

We report transient quasi-periodic oscillations (QPOs) on minute timescales in relativistic, radiative models of the galactic center source Sgr A*. The QPOs result from nonaxisymmetric $m=1$ structure in the accretion flow excited by MHD turbulence. Near-infrared (NIR) and X-ray power spectra show significant peaks at frequencies comparable to the orbital frequency at the innermost stable circular orbit (ISCO) $f_o$. The excess power is associated with inward propagating magnetic filaments inside the ISCO. The amplitudes of the QPOs are sensitive to the electron distribution function. We argue that transient QPOs appear at a range of frequencies in the neighborhood of $f_o$ and that the power spectra, averaged over long times, likely show a broad bump near $f_o$ rather than distinct, narrow QPO features.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Galactic Center Weather Forecast

In accretion-based models for Sgr A* the X-ray, infrared, and millimeter emission arise in a hot, geometrically thick accretion flow close to the black hole. The spectrum and size of the source depend on the black hole mass accretion rate $\dot{M}$. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* that will arrive in summer 2013, $\dot{M}$ may increase from its present value $\dot{M}_0$. We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., which is based on a general relativistic MHD flow model and fully relativistic radiative transfer, for a range of $\dot{M}$. We find that for modest increases in $\dot{M}$ the characteristic ring of emission due to the photon orbit becomes brighter, more extended, and easier to detect by the planned Event Horizon Telescope submm VLBI experiment. If $\dot{M} \gtrsim 8 \dot{M}_0$ this "silhouette of the black hole will be hidden beneath the synchrotron photosphere at 230 GHz, and for $\dot{M} \gtrsim 16 \dot{M}_0$ the silhouette is hidden at 345 GHz. We also find that for $\dot{M} > 2 \dot{M}_0$ the near-horizon accretion flow becomes a persistent X-ray and mid-infrared source, and in the near-infrared Sgr A* will acquire a persistent component that is brighter than currently observed flares.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global GRMHD Simulations of Black Hole Accretion Flows: a Convergence Study

Global, general relativistic magnetohydrodynamic (GRMHD) simulations of nonradiative, magnetized disks are widely used to model accreting black holes. We have performed a convergence study of GRMHD models computed with HARM3D. The models span a factor of 4 in linear resolution, from 96x96x64 to 384x384x256. We consider three diagnostics of convergence: (1) dimensionless shell-averaged quantities such as plasma β; (2) the azimuthal correlation length of fluid variables; and (3) synthetic spectra of the source including synchrotron emission, absorption, and Compton scattering. Shell-averaged temperature is, except for the lowest resolution run, nearly independent of resolution; shell-averaged plasma β decreases steadily with resolution but shows signs of convergence. The azimuthal correlation lengths of density, internal energy, and temperature decrease steadily with resolution but show signs of convergence. In contrast, the azimuthal correlation length of magnetic field decreases nearly linearly with grid size. We argue by analogy with local models, however, that convergence should be achieved with another factor of 2 in resolution. Synthetic spectra are, except for the lowest resolution run, nearly independent of resolution. The convergence behavior is consistent with that of higher physical resolution local model (shearing box) calculations and with the recent nonrelativistic global convergence studies of Hawley et al. (2011).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pair Production in Low Luminosity Galactic Nuclei

Electron-positron pairs may be produced near accreting black holes by a variety of physical processes, and the resulting pair plasma may be accelerated and collimated into a relativistic jet. Here we use a self-consistent dynamical and radiative model to investigate pair production by γγcollisions in weakly radiative accretion flows around a black hole of mass M and accretion rate \dot{M}. Our flow model is drawn from general relativistic magnetohydrodynamic simulations, and our radiation field is computed by a Monte Carlo transport scheme assuming the electron distribution function is thermal. We argue that the pair production rate scales as r^{-6} M^{-1} \dot{M}^{6}. We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in M, \dot{M} parameter space. If \dot{M} is too low the implied pair density over the poles of the black hole is below the Goldreich-Julian density and γγpair production is relatively unimportant; if \dot{M} is too high the models are radiatively efficient. We also argue that for a power-law spectrum the pair production rate should scale with the observables L_X \equiv X-ray luminosity and M as L_X^2 M^{-4}. We confirm this numerically and argue that this relation likely holds even for radiatively efficient flows. The pair production rates are sensitive to black hole spin and to the ion-electron temperature ratio which are fixed in this exploratory calculation. We finish with a brief discussion of the implications for Sgr A* and M87.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

grmonty: a Monte Carlo Code for Relativistic Radiative Transport

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radiative Models of Sgr A* from GRMHD Simulations

Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron temperature ratio Ti/Te. We find that: (1) models with Ti/Te=1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a*; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=10 and i=85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

0 citations0 reviewsNo code linkedOpen access