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

Chi-Kwan Chan

Chi-Kwan Chan 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.IMastro-phastro-ph.EPastro-ph.GAastro-ph.SRgr-qcphysics.flu-dyn

Catalog footprint

What is connected

18works
8topics
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

18 published item(s)

preprint2022arXiv

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

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
preprint2022arXiv

Resolving the inner parsec of the blazar J1924-2914 with the Event Horizon Telescope

The blazar J1924-2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic Center's black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 $μ$as resolution of the EHT. J1924-2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5-11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100-parsec scales. We combine the multi-frequency images of J1924-2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90 degrees between 2.3 and 230 GHz. Linearly polarized intensity images of J1924-2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Brightness Asymmetry of Black Hole Images as a Probe of Observer Inclination

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discretization and Filtering Effects on Black Hole Images Obtained with the Event Horizon Telescope

Interferometers, such as the Event Horizon Telescope (EHT), do not directly observe the images of sources but rather measure their Fourier components at discrete spatial frequencies up to a maximum value set by the longest baseline in the array. Construction of images from the Fourier components or analysis of them with high-resolution models requires careful treatment of fine source structure nominally beyond the array resolution. The primary EHT targets, Sgr A* and M87, are expected to have black-hole shadows with sharp edges and strongly filamentary emission from the surrounding plasma on scales much smaller than those probed by the currently largest baselines. We show that for aliasing not to affect images reconstructed with regularized maximum likelihood methods and model images that are directly compared to the data, the sampling of these images (i.e., their pixel spacing) needs to be significantly finer than the scale probed by the largest baseline in the array. Using GRMHD simulations of black-hole images, we estimate the maximum allowable pixel spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT targets, this corresponds to an angular pixel size of <0.5 microarcseconds. With aliasing under control, we then advocate use of the second-order Butterworth filter with a cut-off scale equal to the maximum array baseline as optimal for visualizing the reconstructed images. In contrast to the traditional Gaussian filters, this Butterworth filter retains most of the power at the scales probed by the array while suppressing the fine image details for which no data exist.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bayesian techniques for comparing time-dependent GRMHD simulations to variable Event Horizon Telescope observations

The Event Horizon Telescope (EHT) is a millimeter-wavelength, very-long-baseline interferometry (VLBI) experiment that is capable of observing black holes with horizon-scale resolution. Early observations have revealed variable horizon-scale emission in the Galactic Center black hole, Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools that explicitly consider the variability in both the data and the models. We develop here a Bayesian method to compare time-resolved simulation images to variable VLBI data, in order to infer model parameters and perform model comparisons. We use mock EHT data based on GRMHD simulations to explore the robustness of this Bayesian method and contrast it to approaches that do not consider the effects of variability. We find that time-independent models lead to offset values of the inferred parameters with artificially reduced uncertainties. Moreover, neglecting the variability in the data and the models often leads to erroneous model selections. We finally apply our method to the early EHT data on Sgr A*.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Particle Acceleration and the Origin of X-ray Flares in GRMHD simulations of Sgr A*

Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform GRMHD simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A General Relativistic Null Hypothesis Test with Event Horizon Telescope Observations of the black-hole shadow in Sgr A*

The half opening angle of a Kerr black-hole shadow is always equal to (5+-0.2)GM/Dc^2, where M is the mass of the black hole and D is its distance from the Earth. Therefore, measuring the size of a shadow and verifying whether it is within this 4% range constitutes a null hypothesis test of General Relativity. We show that the black hole in the center of the Milky Way, Sgr A*, is the optimal target for performing this test with upcoming observations using the Event Horizon Telescope. We use the results of optical/IR monitoring of stellar orbits to show that the mass-to-distance ratio for Sgr A* is already known to an accuracy of +-4%. We investigate our prior knowledge of the properties of the scattering screen between Sgr A and the Earth, the effects of which will need to be corrected for in order for the black-hole shadow to appear sharp against the background emission. Finally, we explore an edge detection scheme for interferometric data and a pattern matching algorithm based on the Hough/Radon transform and demonstrate that the shadow of the black hole at 1.3 mm can be localized, in principle, to within ~9%. All these results suggest that our prior knowledge of the properties of the black hole, of scattering broadening, and of the accretion flow can only limit this General Relativistic null hypothesis test with Event Horizon Telescope observations of Sgr A* to 10%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Variability and mm/IR flares in GRMHD Models of Sgr A* from Strong-Field Gravitational Lensing

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Power of Imaging: Constraining the Plasma Properties of GRMHD Simulations using EHT Observations of Sgr A*

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black-holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast GPU-based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all twelve models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black-holes and their thermodynamic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRay: a Massively Parallel GPU-Based Code for Ray Tracing in Relativistic Spacetimes

We introduce GRay, a massively parallel integrator designed to trace the trajectories of billions of photons in a curved spacetime. This GPU-based integrator employs the stream processing paradigm, is implemented in CUDA C/C++, and runs on nVidia graphics cards. The peak performance of GRay using single precision floating-point arithmetic on a single GPU exceeds 300 GFLOP (or 1 nanosecond per photon per time step). For a realistic problem, where the peak performance cannot be reached, GRay is two orders of magnitude faster than existing CPU-based ray tracing codes. This performance enhancement allows more effective searches of large parameter spaces when comparing theoretical predictions of images, spectra, and lightcurves from the vicinities of compact objects to observations. GRay can also perform on-the-fly ray tracing within general relativistic magnetohydrodynamic algorithms that simulate accretion flows around compact objects. Making use of this algorithm, we calculate the properties of the shadows of Kerr black holes and the photon rings that surround them. We also provide accurate fitting formulae of their dependencies on black hole spin and observer inclination, which can be used to interpret upcoming observations of the black holes at the center of the Milky Way, as well as M87, with the Event Horizon Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of Saturated Energy Condensation in Two-Dimensional Turbulence

In two-dimensional forced Navier-Stokes turbulence, energy cascades to the largest scales in the system to form a pair of coherent vortices known as the Bose condensate. We show, both numerically and analytically, that the energy condensation saturates and the system reaches a statistically stationary state. The time scale of saturation is inversely proportional to the viscosity and the saturation energy level is determined by both the viscosity and the force. We further show that, without sufficient resolution to resolve the small-scale enstrophy spectrum, numerical simulations can give a spurious result for the saturation energy level. We also find that the movement of the condensate is similar to the motion of an inertial particle with an effective drag force. Furthermore, we show that the profile of the saturated coherent vortices can be described by a Gaussian core with exponential wings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MHD Turbulence in Accretion Disk Boundary Layers

The physical modeling of the accretion disk boundary layer, the region where the disk meets the surface of the accreting star, usually relies on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity, widely adopted in astrophysics, satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability is inefficient in this inner disk region. I will discuss the results of a recent study on the generation of hydromagnetic stresses and energy density in the boundary layer around a weakly magnetized star. Our findings suggest that although magnetic energy density can be significantly amplified in this region, angular momentum transport is rather inefficient. This seems consistent with the results obtained in numerical simulations and suggests that the detailed structure of turbulent MHD boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Hydromagnetic Stresses in Accretion Disk Boundary Layers

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light into physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Class of Physically Motivated Closures for Radiation Hydrodynamics

Radiative transfer and radiation hydrodynamics use the relativistic Boltzmann equation to describe the kinetics of photons. It is difficult to solve the six-dimensional time-dependent transfer equation unless the problem is highly symmetric or in equilibrium. When the radiation field is smooth, it is natural to take angular moments of the transfer equation to reduce the degrees of freedom. However, low order moment equations contain terms that depend on higher order moments. To close the system of moment equations, approximations are made to truncate this hierarchy. Popular closures used in astrophysics include flux limited diffusion and the M1 closure, which are rather ad hoc and do not necessarily capture the correct physics. In this paper, we propose a new class of closures for radiative transfer and radiation hydrodynamics. We start from a different perspective and highlight the consistency of a fully relativistic formalism. We present a generic framework to approximate radiative transfer based on relativistic Grad's moment method. We then derive a 14-field method that minimizes unphysical photon self-interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Oscillations of the Inner Regions of Viscous Accretion Disks

Although quasi-periodic oscillations (QPOs) have been discovered in different X-ray sources, their origin is still a matter of debate. Analytical studies of hydrodynamic accretion disks have shown three types of trapped global modes with properties that appear to agree with the observations. However, these studies take only linear effects into account and do not address the issues of mode excitation and decay. Moreover, observations suggest that resonances between modes play a crucial role. A systematic, numerical study of this problem is therefore needed. In this paper, we use a pseudo-spectral algorithm to perform a parameter study of the inner regions of hydrodynamic disks. By assuming alpha-viscosity, we show that steady state solutions rarely exist. The inner edges of the disks oscillate and excite axisymmetric waves. In addition, the flows inside the inner edges are sometimes unstable to non-axisymmetric perturbations. One-armed, or even two-armed, spirals are developed, which provides a plausible explanation for the high-frequency QPOs observed from accreting black holes. When the Reynolds numbers are above certain critical values, the inner disks go through some transient turbulent states characterized by strong trailing spirals; while large-scale leading spirals developed in the outer disks. We compared our numerical results with standard thin disk oscillation models. Although the non-axisymmetric features have their analytical counterparts, more careful study is needed to explain the axisymmetric oscillations.

0 citations0 reviewsNo code linkedOpen access