Engineering topics

Krolik, Julian H.

Publications and source records attributed to Krolik, Julian H..

At least 19 records

DOE OSTI2024

Accretion onto a Supermassive Black Hole Binary before Merger

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20 M (G = c = 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate $\dot{M}$ onto the black holes first decreases and then plateaus, dropping by only a factor of ~3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, "sloshing" streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless level $\phi$ ~ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PATCHWORK MHD, decreasing the computational expense to ~3% of conventional single-grid methods' cost.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

Disk Emission from Magnetohydrodynamic Simulations of Spinning Black Holes

We present the results of a new series of global, three-dimensional, relativistic magnetohydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of H/R approx. 0.05 and spin parameters of a/M = 0, 0.5, 0.9, and 0.99. Using the ray-tracing code Pandurata, we generate broadband thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Finally, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well known to the continuum-fitting method of measuring black hole spin.

accretion↗

NASA NTRS2013

A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes

We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption mechanisms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code and also present results from a number of numerical tests to assess its accuracy and convergence.

accretion disks↗

NASA NTRS2012

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, nonrotating black hole. For the first time, we are able to explain from first principles in a self-consistent way the X-ray spectra observed from stellar-mass black holes, including a thermal peak, Compton reflection hump, power-law tail, and broad iron line. Varying only the mass accretion rate, we are able to reproduce the low/hard, steep power-law, and thermal-dominant states seen in most galactic black hole sources. The temperature in the corona is T(sub e) ~ 10 keV in a boundary layer near the disk and rises smoothly to T(sub e) greater than or approximately 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to approximately equal to 6M as the accretion rate decreases, we find that the shape of the Fe Ka line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

Schnittman, Jeremy D.↗

NASA NTRS2011

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present new global calculations of X-ray spectra from fully relativistic magneto-hydrodynamic (MHO) simulations of black hole (BH) accretion disks. With a self consistent radiative transfer code including Compton scattering and returning radiation, we can reproduce the predominant spectral features seen in decades of X-ray observations of stellar-mass BHs: a broad thermal peak around 1 keV, power-law continuum up to >100 keV, and a relativistically broadened iron fluorescent line. By varying the mass accretion rate, different spectral states naturally emerge: thermal-dominant, steep power-law, and low/hard. In addition to the spectral features, we briefly discuss applications to X-ray timing and polarization.

Schnittman, Jeremy D.↗

NASA NTRS2005

Physics-Based Spectra of Accretion Disks around Black Holes

The purpose of this grant was to begin the process of deriving the light output of accretion disks around black holes directly from the actual processes that inject heat into the accreting matter, rather than from guessed dependences of heating rate on physical parameters. At JHU, the effort has focussed so far on models of accretion onto "intermediate mass black holes", a possible class of black holes, examples of which may have recently been discovered in nearby galaxies. There, Krolik and his student (Yawei Hui) have computed stellar atmospheres for uniformly-heated disks around this class of black holes. Their models serve two purposes: they are the very first serious attempts to compute the spectrum from accreting black holes in this mass range; and a library of such models can be used later in this program as contrasts for those computed on the basis of real disk dynamics. The output from these local disk calculations has also been successfully coupled to a program that applies the appropriate relativistic transformations and computes photon trajectories in order to predict the spectrum received by observers located at different polar angles. The principal new result of these calculations is the discovery of potentially observable ionization edges of H-like C and O at frequencies near the peak in flux from these objects. Most of the grant money at UCSB was spent on supporting graduate student Shane Davis. In addition. some money was spent on supporting two other students: Ari Socrates (now a Hubble Fellow at Princeton), and Laura Melling. Davis spent the year constructing stellar atmosphere models of accretion disks appropriate for the high/soft (thermal) state of black hole X-ray binaries. As with AGN models published previously by our collaboration with NASA support. our models include a complete general relativistic treatment of both the disk structure and the propagation of photons from the disk to a distant observer. They also include all important continuum opacity sources, including Compton scattering and bound-free opacity from abundant metal species. The principal new result is that bound-free opacity is very significant in altering the continuum spectral shape, resulting for example in quite different "color correction factors" compared to those predicted previously. In addition, the models predict a relationship between luminosity and inner disk temperature that is, for the first time, in accord with that observed. The primary purpose of the grant was to incorporate more realistic accretion disk physics, learned largely from simulations, into such spectral models. The Davis et al. paper includes consideration of a vertical dissipation profile computed from radiation magneto-hydrodynamic simulations of MRI turbulence by N. J. Turner (2004). So long as the disk is effectively thick, such dissipation profiles do not affect the predicted spectrum significantly. (More work needs to be done on these simulations, however.) A potentially more serious issue is that MRI turbulence produces substantial inhomogeneities, as do photon bubble instabilities. These inhomogeneities can affect the spectra by enhancing the effects of absorption opacity over scattering opacity. We have done some preliminary Monte Carlo calculations to explore these effects.

Krolik, Julian H.↗

NASA NTRS2004

[Predicting Spectra of Accretion Disks Around Galactic Black Holes]

The purpose of this grant was to construct detailed atmosphere solutions in order to predict the spectra of accretion disks around Galactic black holes. Our plan of action was to take an existing disk atmosphere code (TLUSTY, created by Ivan Hubeny) and introduce those additional physical processes necessary to make it applicable to disks of this variety. These modifications include: treating Comptonization; introducing continuous opacity due to heavy elements; incorporating line opacity due to heavy elements; adopting a disk structure that reflects readjustments due to radiation pressure effects; and injecting heat via a physically-plausible vertical distribution.

Krolik, Julian H.↗

NASA NTRS2002

The X-ray Reflectors in the Nucleus of the Seyfert Galaxy NGC 1068

Based on observations of the Seyfert nucleus in NGC 1068 with ASCA, RXTE and BeppoSAX, we report the discovery of a flare (increase in flux by a factor of approximately 1.6) in the 6.7 keV Fe K line component between observations obtained four months apart, with no significant change in the other (6.21, 6.4, and 6.97 keV) Fe Kalpha line components. During this time, the continuum flux decreased by approximately 20%. The RXTE spectrum requires an Fe K absorption edge near 8.6 keV (Fe XXIII- XXV). The spectral data indicate that the 2-10 keV continuum emission is dominated (approximately 2/3 of the luminosity) by reflection from a previously unidentified region of warm, ionized gas located approximately or less than 0.2 pc from the AGN. The remaining approximately 1/3 of the observed X-ray emission is reflected from optically thick, neutral gas. The coronal gas in the inner Narrow-Line Region (NLR) and/or the cold gas at the inner surface of the obscuring 'torus' are possible cold reflectors. The inferred properties of the warm reflector are: size (diameter) approximately or less than 0.2 pc, gas density n approximately or greater than 10(exp 5.5)/cu cm, ionization parameter xi is approximately 10(exp 3.5) erg cm s(exp -1), and covering fraction 0.003 (L(sub 0)/ 10(exp 43.5) erg s(exp -1)(exp -1) less than (omega/4pi) less than 0.024 (L(sub 0)/ 10(exp 43.5) erg s(exp -1) (exp -1) where L(sub 0) is the intrinsic 2-10 keV X-ray luminosity of the AGN. We suggest that the warm reflector gas is the source of the (variable) 6.7 keV Fe line emission, and the 6.97 keV Fe line emission. The 6.7 keV line flare is assumed to be due to an increase in the emissivity of the warm reflector gas from a decrease (by 20-30%) in L(sub 0). The properties of the warm reflector are most consistent with an intrinsically X-ray weak AGN with L(sub 0) approximately equals 10(exp 43.0) erg s(exp -1). The optical and UV emission that scatters from the warm reflector into our line of sight is required to suffer strong extinction, which can be reconciled if the line-of-sight skims the outer surface of the torus. Thermal bremsstrahlung radio emission from the warm reflector may be detectable in VLBA radio maps of the NGC 1068 nucleus.

Colbert, Edward J. M.↗

NASA NTRS1998

Variability Analysis and the Structure of Active Galactic Nuclei

This five-year Long-Term Space Astrophysics grant provided the support for several major steps in advancing our knowledge of the internal structure of active galactic nuclei. The single largest portion of this program had to do with the development and application of techniques for "reverberation mapping", the use of spectral monitoring of several different bands related by radiation reprocessing to infer the internal geometry of sources. Major steps were taken in this regard, particularly in establishing the distribution in radius of emission line material, and in relating the apparent reprocessing of continuum bands to the underlying structure of the accretion disk. Another major effort built directly upon these results. Once the case for continuum reprocessing was made by the monitoring, it next behooved us to understand the spectral output of AGN as a result of this reprocessing. As a result, our view of continuum production in AGN is now much better focussed on the key problems. A third focus of effort had to do with the nature of X-ray variability in AGN, and what it can tell us about the dynamics of extremely hot material in the immediate outskirts of the supermassive black holes that form the central engines of active galactic nuclei. In addition to these primary efforts, this grant also supported many other, smaller projects. Several of these were demonstrations of how the material spewed out of AGN in relativistic.ets generate the radiation by which we observe them. J Finally, the portion of this study that had to do with continuum production by accretion disks in AGN led naturally to several papers in which new developments were presented having to do with "advection-dominated accretion disks", those disks in which accretion appears to proceed at a substantial rate, but in which radiation processes are weak.

Krolik, Julian H.↗

NASA NTRS1998

What is the Nature of Accretion in Active Galactic Nuclei?

The purpose of this grant was to support theoretical research on the nature of accretion in active galactic nuclei. In the brief time of the award, four papers that appeared in refereed journals were written, as well as two invited reviews in conference proceedings. These papers significantly advanced our understanding of the structure of the most important parts of bright accretion disks around accreting black holes, such as active galactic nuclei.

Krolik, Julian H.↗

NASA NTRS1998

Variability Analysis and the Structure of Active Galactic Nuclei

This five-year Long-Term Space Astrophysics grant provided the support for several major steps in advancing our knowledge of the internal structure of active galactic nuclei. The single largest portion of this program had to do with the development and application of techniques for 'reverberation mapping, the use of spectral monitoring of several different bands related by radiation reprocessing to infer the internal geometry of sources. Major steps were taken in d-ds regard, particularly in establishing the distribution in radius of emission line material, and in relating the apparent reprocessing of continuum bands to the underlying structure of the accretion disk. Another major effort built directly upon these results. Once the case for continuum reprocessing was made by the monitoring, it next behooved us to understand the spectral output of AGN as a result of this reprocessing. As a result, our view of continuum production in AGN is now much better focussed on the key problems. A third focus of effort had to do with the nature of X-ray variability in AGN, and what it can tell us about the dynamics of extremely hot material in the immediate outskirts of the supermassive black holes that form the central engines of active galactic nuclei. In addition to these primary efforts, this grant also supported many other, smaller projects. Several of these were demonstrations of how the material spewed out of AGN in relativistic jets generate the radiation by which we observe them. Finally, the portion of this study that had to do with continuum production by accretion disks in AGN led naturally to several papers in which new developments were presented having to do with 'advection-dominated accretion disks', those disks in which accretion appears to proceed at a substantial rate, but in which radiation processes are weak.

Krolik, Julian H.↗

NASA NTRS1995

Reverberation mapping by regularized linear inversion

Reverberation mapping of active galactic nucleus (AGN) emission-line regions requires the numerical deconvolution of two time series. We suggest the application of a new method, regularized linear inversion, to the solution of this problem. This method possesses many good features; it imposes no restrictions on the sign of the response function; it can provide clearly defined uncertainty estimates; it involves no guesswork about unmeasured data; it can give a clear indication of when the underlying convolution model is inadequate; and it is computationally very efficient. Using simulated data, we find the minimum S/N and length of the time series in order for this method to work satisfactorily. We also define guidelines for choosing the principal tunable parameter of the method and for interpreting the results. Finally, we reanalyze published data from the 1989 NGC 5548 campaign using this new method and compare the results to those previously obtained by maximum entropy analysis. For some lines we find good agreement, but for others, especially C III lambda(1909) and Si IV lambda(1400), we find significant differences. These can be attributed to the inability of the maximum entropy method to find negative values of the response function, but also illustrate the nonuniqueness of any deconvolution technique. We also find evidence that certain line light curves (e.g., C IV lambda(1549)) cannot be fully described by the simple linear convolution model.

Krolik, Julian H.↗

NASA NTRS1994

X-ray refelection from photoionized media in active galactic nuclei

We calculate the spectrum of X-ray radiation and reprocessed by a partly ionized optically thick medium in an active galactic nucleus. We self-consistently calculate the ionization balance and thermal balance in the medium along with the distribution of X-ray intensity with optical depth. In addition to absorption or scattering of the incident X-rays, we also compute the spectrum of X-rays emitted by the material, including lines, edges, and bremsstrahlung. The albedo of the medium depends primarily on the X-ray ionization parameter (ratio of incident flux to gas density, zeta(sub Chi), and secondarily on the UV flux generated by dissipation inside the disk; we locate the critical range of zeta(sub Chi) over which the albedo increases from small to nearly unity. While the continuum reflection is very weak below 10 keV when zeta(sub Chi) is small, significnat fluxes are emitted in atomic lines and edges in this energy range. In the limit of large zeta(sub Chi), the albedo below 10 keV increases, but reflection in this band is never gray: some photoelectric absorption remains up to rather large values of zeta(sub Chi), while at still higher values, inverse Compton scattering amplifies the soft X-ray flux. These features are sufficiently sharp that current and near-future X-ray experiments should permit diagnostic measures of zeta(sub Chi).

Zycki, Piotr T.↗

NASA NTRS1994

Hubble Space Telescope imaging of the Seyfert 2 galaxy NGC 2110

Hubble Space Telescope (HST) images of the Seyfert 2 galaxy NGC 2110 in the emission lines of (O III) lambda lambda-4959, 5007 and H-alpha + (N II) lambda lambda-6548, 6584 and the nearby continuum are presented. After deconvolution, the resolution of these images is approximately 0.1 sec, corresponding to a linear scale of approximately 23 pc at the assumed distance (46.8 Mpc) to this galaxy. The dominant line emission is associated with a strongly curved, narrow, 1 sec (230 pc) long feature emanating from the nucleus to the north. The morphology of this feature suggests that it may represent gas flowing outward from the nucleus (perhaps as a 'jet'). While this optical emission-line feature and the jetlike radio emission extend along a similar P.A., their detailed brightness distributions are loosely anticorrelated. In particular, there is no line emission associated with the off-nuclear 6 cm peaks and only weak radio emission along the optical 'jet'. On a larger scale, there is an S-shaped distribution of ionized gas (seen best in H-alpha) extending approximately 4 sec north and south of the nucleus. Ground-based kinematical data suggest this larger scale component represents normally rotating, ambient, interstellar gas photoionized by the central source. A continuum color map derived from HST observations shows evidence for an unresolved region (scale less than 23 pc) of very strong reddening at the location of the optical and radio nucleus. The location and size of the reddened region are consistent with obscuration by a dusty torus, but the observed reddening provides only a lower limit (A(sub V) greater than or equal to 4.6(sub -1.7 sup +1.2) mag) on the amount of obscuration to the nucleus. Nevertheless, this obscuration is sufficient to account for the 'energy deficit' previously reported in this object.

Mulchaey, John S.↗

NASA NTRS1994

Relativistic induced Compton scattering in synchrotron self-absorbed sources

We present a model of the radio emission from synchrotron self-absorbed sources, including the effects of induced Compton scattering by the relativistic electrons in the source. Order of magnitude estimates show that stimulated scattering becomes the dominant absorption process when (kT(sub B)/m(sub e) c(exp 2))Tau(sub T) greater than or approximately 1.0. Numerical simulations demonstrate that relativistic induced Compton scattering limits the brightness temperature of a self-absorbed synchroton sources to T(sub B) less than or approximately 2 x 10(exp 11) nu(exp -1/p + 3) (sub zero, 9) gamma(exp p + 2/p + 5) (sub min) K, where gamma(sub min) is the low-energy cutoff to the relativistic electron distribution with a power-law index of p. It can also significantly flatten the radio spectrum. The radio spectrum of the core of 3C 279 is well matched by a model in which stimulated scattering is important, and the additional constraint T(sub B) less than 2 x 10(exp 11)K may be important to the interpretation of the broadband spectra in variable extragalactic compact radio sources. Stimulated scattering reduces the amplitude of the radio frequency variability relative to the X-ray variability, an effect which can be detected by multiwavelength variability studies. Data for a sample of resolved compact radio cores indicate that it is inconsistent to neglect induced Compton scattering when inferring the physical parameters of the sources. The necessary generalizations to the standard synchrotron self-Compton theory are p resented.

Sincell, Mark W.↗

NASA NTRS1994

Hubble Space Telescope images of the Seyfert galaxies NGC 5929 and MCG 8-11-11

We present the initial results of a program to obtain high resolution images of Seyfert galaxies with the Planetary Camera aboard the Hubble Space Telescope (HST). In this paper we discuss the images of the type 2 Seyfert NGC 5929 and the type 1.5 Seyfert MCG 8-11-11 (= UGC 3374). The images were obtained in the emission lines of (O III) lambda lambda 4959 and 5007 A and H alpha + (N II) lambda lambda 6548 and 6583 A and their adjacent continua. The high-excitation gas in the narrow line region (NLR) of NGC 5929 is resolved into individual clouds in the central 1 sec .5. Although the (O III) emission is clearly not spherically symmetric with respect to the nucleus, it does not define a distinct 'bicone' morphology, as observed by the HST in a few other Seyfert galaxies. We find no direct evidence for the reddening and/or obscuration effects characteristic of a dusty torus, which, in the context of 'unified models', is expected to obscure the active galactic nucleus (AGN) in type 2 Seyfert galaxies. The correspondence between the emission line gas and the radio morphology suggests that the structure of the NLR in NGC 5929 is governed by matter ejected from the AGN. A comparison of the recombination rate of hydrogen in the brightest emission line cloud with an upper limit on the ionizing luminosity emitted toward Earth provides no evidence that the central ionizing source radiates anisotropically. The images of MCG 8-11-11 show only an unresolved nuclear source. No emission line gas associated with the extended radio source is detected. We estimate upper limits on the intensity of extended line emission in this galaxy and examine their significance.

Bower, Gary A.↗

NASA NTRS1994

X-ray bumps, iron K-alpha lines, and X-ray suppression by obscuring tori in Seyfert galaxies

We investigate the X-ray spectral properties of unobscured type 1 and obscured type 2 Seyferts as predicted by the unified Seyfert scheme. We consider the reprocessing of X-ray photons by photoelectric absorption, iron fluorescence, and Compton downscattering in the obscuring tori surrounding these active nuclei, and compute by Monte Carlo methods the reprocessed spectra as a function of the viewing angle. Depending on the optical depth and shape of the torus, and on the viewing angle, the X-ray flux can be suppressed by substantial factors when our line of sight is obscured. We show that an immediate consequence of the existence of an obscuring thick torus is the production in the spectra of type 1 Seyfert galaxies of a bump in the continuum above 10-20 keV and an Fe K-alpha line with significant equivalent width. In those type 2 Seyferts for which the hard X-ray spectrum has been substantially suppressed, the equivalent width of the Fe K-alpha line in the transmitted spectrum can be very large.

Krolik, Julian H.↗

NASA NTRS1993

Active galactic nuclei make the cosmic X-ray background

While it has long been known that active galactic nuclei (AGNs) make up a significant part of the X-ray background near 1 keV, it has been generally thought that the typical AGN X-ray spectrum did not fit the higher energy background spectrum. More recent data (from Ginga, SIGMA, and OSSE) show that the X-ray spectra of radio-quiet AGNs do in fact resemble the X-ray background spectrum. These data are well described (from 2 to 100 keV) by a model in which the X-ray photons are made by thermal Comptonization, and a fraction (around a half) of the X-rays are reflected by a cool surface. Using the observed evolution of X-ray-bright AGNs, we show that the known AGN population, radiating the spectrum observed at low redshift and integrated in redshift back to z about 4, is capable of producing the entire X-ray background from 2 to 100 keV.

Zdziarski, Andrzej A.↗