BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Julian H. Krolik

Julian H. Krolik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcphysics.comp-phphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
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

10 published item(s)

preprint2022arXiv

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by an $m=1$ overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases, we find most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are $8$ times stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

What powers the radio emission in TDE AT2019dsg: a long-lived jet or the disruption itself?

The tidal disruption event AT2019dsg was observed from radio to X-rays and was possibly accompanied by a high-energy neutrino. Previous interpretations have focused on continued injection by a central engine as the source of energy for radio emission. We show that continuous energy injection is unnecessary; the radio data can be explained by a single ejection of plasma that supplies all the energy needed. To support this assertion, we analyze the synchrotron self-absorbed spectra in terms of the equipartition model. Similar to previous analyses, we find that the energy in the radio-emitting region increases approximately $\propto t^{0.7}$ and the lengthscale of this region grows $\propto t$ at a rate $\simeq0.06c$. This event resembles the earliest stage of a supernova remnant: because the ejected mass is much greater than the shocked external mass, its velocity remains unchanged, while the energy in shocked gas grows with time. The radio-emitting material gains energy from the outflow, not continuing energy injection by the central object. Although energy injection from an accreting BH cannot be completely excluded, the energy injection rate is very different from the fallback luminosity, and maintaining constant outflow velocity requires fine-tuning demanding further physical explanation. If the neutrino association is real, the energy injection needed is much greater than for the radio emission, suggesting that the detected neutrino did not arise from the radio-emitting region.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-energy emission from tidal disruption events in active galactic nuclei

Tidal disruption events (TDEs) taking place in active galactic nuclei (AGNs) are different from ordinary TDEs. In these events, the returning tidal debris stream drills through the pre-existing AGN accretion disk near the stream pericenter, destroying the inner disk in the process, and then intersects with the disk a second time at radii ranging from a few to hundreds of times the pericenter distance. The debris dynamics of such TDEs, and hence their appearance, are distinct from ordinary TDEs. Here we explore the observational signatures of this "second impact" of the stream with the disk. Strong shocks form as the dilute stream is stopped by the denser disk. Compton cooling of the shocked material produces hard X-rays, even soft gamma-rays, with most of the energy emitted between ~10 keV and 1 MeV. The luminosity follows the mass-return rate, peaking between ~$10^{42}$ and $10^{44}$ erg/s. The X-ray hardness and the smoothness of the light curve provide possible means for distinguishing the second impact from ordinary AGN flares, which exhibit softer spectra and more irregular light curves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Future Simulations of Tidal Disruption Events

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprised of an arbitrary number of moving, local meshes, or patches, which are free to employ their own resolution, coordinate system/topology, physics equations, reference frame, and in our new approach, numerical method. Each local patch exchanges boundary data with a single global patch on which all other patches reside through a client-router-server parallelization model. In generalizing Patchwork to be compatible with arbitrary order time integration, PatchworkMHD and PatchworkWave have significantly improved the interpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in the original Patchwork code. Furthermore, we extend Patchwork to be multimethod by allowing multiple state vectors to be updated simultaneously, with each state vector providing its own interpatch interpolation and transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partial differential equations. We demonstrate our changes through the implementation of a scalar wave toy-model that is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Radiative Emission Mechanisms of Tidal Disruption Events

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tidal disruption events in active galactic nuclei

A fraction of tidal disruption events (TDEs) occur in active galactic nuclei (AGNs) whose black holes possess accretion disks; these TDEs can be confused with common AGN flares. The disruption itself is unaffected by the disk, but the evolution of the bound debris stream is modified by its collision with the disk when it returns to pericenter. The outcome of the collision is largely determined by the ratio of the stream mass current to the azimuthal mass current of the disk rotating underneath the stream footprint, which in turn depends on the mass and luminosity of the AGN. To characterize TDEs in AGNs, we simulated a suite of stream--disk collisions with various mass current ratios. The collision excites shocks in the disk, leading to inflow and energy dissipation orders of magnitude above Eddington; however, much of the radiation is trapped in the inflow and advected into the black hole, so the actual bolometric luminosity may be closer to Eddington. The emergent spectrum may not be thermal, TDE-like, or AGN-like. The rapid inflow causes the disk interior to the impact point to be depleted within a fraction of the mass return time. If the stream is heavy enough to penetrate the disk, part of the outgoing material eventually hits the disk again, dissipating its kinetic energy in the second collision; another part becomes unbound, emitting synchrotron radiation as it shocks with surrounding gas.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Magnetorotational instability in eccentric disks

Eccentric disks arise in such astrophysical contexts as tidal disruption events, but it is unknown whether the magnetorotational instability (MRI), which powers accretion in circular disks, operates in eccentric disks as well. We examine the linear evolution of unstratified, incompressible MRI in an eccentric disk orbiting a point mass. We consider vertical modes of wavenumber $k$ on a background flow with uniform eccentricity $e$ and vertical Alfvén speed $v_\mathrm A$ along an orbit with mean motion $n$. We find two mode families, one with dominant magnetic components, the other with dominant velocity components; the former is unstable at $(1-e)^3f^2\lesssim3$, where $f\equiv kv_\mathrm A/n$, the latter at $e\gtrsim0.8$. For $f^2\lesssim3$, MRI behaves much like in circular disks, but the growth per orbit declines slowly with increasing $e$; for $f^2\gtrsim3$, modes grow by parametric amplification, which is resonant for $0<e\ll1$. MRI growth and the attendant angular momentum and energy transport happen chiefly near pericenter, where orbital shear dominates magnetic tension.

0 citations0 reviewsNo code linkedOpen access