Normalized to: Kawashima, T.
[1]
oai:arXiv.org:2004.01161 [pdf] - 2074562
SYMBA: An end-to-end VLBI synthetic data generation pipeline
Roelofs, F.;
Janssen, M.;
Natarajan, I.;
Deane, R.;
Davelaar, J.;
Olivares, H.;
Porth, O.;
Paine, S. N.;
Bouman, K. L.;
Tilanus, R. P. J.;
van Bemmel, I. M.;
Falcke, H.;
Akiyama, K.;
Alberdi, A.;
Alef, W.;
Asada, K.;
Azulay, R.;
Baczko, A.;
Ball, D.;
Baloković, M.;
Barrett, J.;
Bintley, D.;
Blackburn, L.;
Boland, W.;
Bower, G. C.;
Bremer, M.;
Brinkerink, C. D.;
Brissenden, R.;
Britzen, S.;
Broderick, A. E.;
Broguiere, D.;
Bronzwaer, T.;
Byun, D.;
Carlstrom, J. E.;
Chael, A.;
Chan, C.;
Chatterjee, S.;
Chatterjee, K.;
Chen, M.;
Chen, Y.;
Cho, I.;
Christian, P.;
Conway, J. E.;
Cordes, J. M.;
Crew, G. B.;
Cui, Y.;
De Laurentis, M.;
Dempsey, J.;
Desvignes, G.;
Dexter, J.;
Doeleman, S. S.;
Eatough, R. P.;
Fish, V. L.;
Fomalont, E.;
Fraga-Encinas, R.;
Friberg, P.;
Fromm, C. M.;
Gómez, J. L.;
Galison, P.;
Gammie, C. F.;
García, R.;
Gentaz, O.;
Georgiev, B.;
Goddi, C.;
Gold, R.;
Gu, M.;
Gurwell, M.;
Hada, K.;
Hecht, M. H.;
Hesper, R.;
Ho, L. C.;
Ho, P.;
Honma, M.;
Huang, C. L.;
Huang, L.;
Hughes, D. H.;
Ikeda, S.;
Inoue, M.;
Issaoun, S.;
James, D. J.;
Jannuzi, B. T.;
Jeter, B.;
Jiang, W.;
Johnson, M. D.;
Jorstad, S.;
Jung, T.;
Karami, M.;
Karuppusamy, R.;
Kawashima, T.;
Keating, G. K.;
Kettenis, M.;
Kim, J.;
Kim, J.;
Kim, J.;
Kino, M.;
Koay, J. Yi;
Koch, P. M.;
Koyama, S.;
Kramer, M.;
Kramer, C.;
Krichbaum, T. P.;
Kuo, C.;
Lauer, T. R.;
Lee, S.;
Li, Y.;
Li, Z.;
Lindqvist, M.;
Lico, R.;
Liu, K.;
Liuzzo, E.;
Lo, W.;
Lobanov, A. P.;
Loinard, L.;
Lonsdale, C.;
Lu, R.;
MacDonald, N. R.;
Mao, J.;
Markoff, S.;
Marrone, D. P.;
Marscher, A. P.;
Martí-Vidal, I.;
Matsushita, S.;
Matthews, L. D.;
Medeiros, L.;
Menten, K. M.;
Mizuno, Y.;
Mizuno, I.;
Moran, J. M.;
Moriyama, K.;
Moscibrodzka, M.;
Müller, C.;
Nagai, H.;
Nagar, N. M.;
Nakamura, M.;
Narayan, R.;
Narayanan, G.;
Neri, R.;
Ni, C.;
Noutsos, A.;
Okino, H.;
Ortiz-León, G. N.;
Oyama, T.;
Özel, F.;
Palumbo, D. C. M.;
Patel, N.;
Pen, U.;
Pesce, D. W.;
Piétu, V.;
Plambeck, R.;
PopStefanija, A.;
Prather, B.;
Preciado-López, J. A.;
Psaltis, D.;
Pu, H.;
Ramakrishnan, V.;
Rao, R.;
Rawlings, M. G.;
Raymond, A. W.;
Rezzolla, L.;
Ripperda, B.;
Rogers, A.;
Ros, E.;
Rose, M.;
Roshanineshat, A.;
Rottmann, H.;
Roy, A. L.;
Ruszczyk, C.;
Ryan, B. R.;
Rygl, K. L. J.;
Sánchez, S.;
Sánchez-Arguelles, D.;
Sasada, M.;
Savolainen, T.;
Schloerb, F. Peter;
Schuster, K.;
Shao, L.;
Shen, Z.;
Small, D.;
Sohn, B. Won;
SooHoo, J.;
Tazaki, F.;
Tiede, P.;
Titus, M.;
Toma, K.;
Torne, P.;
Trent, T.;
Trippe, S.;
Tsuda, S.;
van Langevelde, H. J.;
van Rossum, D. R.;
Wagner, J.;
Wardle, J.;
Ward-Thompson, D.;
Weintroub, J.;
Wex, N.;
Wharton, R.;
Wielgus, M.;
Wong, G. N.;
Wu, Q.;
Young, A.;
Young, K.;
Younsi, Z.;
Yuan, F.;
Yuan, Y.;
Zensus, J. A.;
Zhao, G.;
Zhao, S.;
Zhu, Z.
Submitted: 2020-04-02
Realistic synthetic observations of theoretical source models are essential
for our understanding of real observational data. In using synthetic data, one
can verify the extent to which source parameters can be recovered and evaluate
how various data corruption effects can be calibrated. These studies are
important when proposing observations of new sources, in the characterization
of the capabilities of new or upgraded instruments, and when verifying
model-based theoretical predictions in a comparison with observational data. We
present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a
novel synthetic data generation pipeline for Very Long Baseline Interferometry
(VLBI) observations. SYMBA takes into account several realistic atmospheric,
instrumental, and calibration effects. We used SYMBA to create synthetic
observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has
recently captured the first image of a black hole shadow. After testing SYMBA
with simple source and corruption models, we study the importance of including
all corruption and calibration effects. Based on two example general
relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed
case studies to assess the attainable image quality with the current and future
EHT array for different weather conditions. The results show that the effects
of atmospheric and instrumental corruptions on the measured visibilities are
significant. Despite these effects, we demonstrate how the overall structure of
the input models can be recovered robustly after performing calibration steps.
With the planned addition of new stations to the EHT array, images could be
reconstructed with higher angular resolution and dynamic range. In our case
study, these improvements allowed for a distinction between a thermal and a
non-thermal GRMHD model based on salient features in reconstructed images.
[2]
oai:arXiv.org:2003.07569 [pdf] - 2065521
Pulsed fraction of super-critical column accretion flows onto neutron
stars: modeling of ultraluminous X-ray pulsars
Submitted: 2020-03-17
We calculate the pulsed fraction (PF) of the super-critical column accretion
flows onto magnetized neutron stars (NSs), of which the magnetic axis is
misaligned with the rotation axis, based on the simulation results by Kawashima
et al.(2016, PASJ, 68, 83). Here, we solve the geodesic equation for light in
the Schwarzschild spacetime in order to take into account the light bending
effect. The gravitational redshift and the relativistic doppler effect from gas
motions of the accretion columns are also incorporated. The pulsed emission
appears since the observed luminosity, which exceeds the Eddington luminosity
for the stellar-mass black holes, periodically changes via precession of the
column caused by the rotation of the NS. The PF tends to increase as
$\theta_{\rm obs}$ approaching to $\theta_{\rm B}$, where $\theta_{\rm obs}$
and $\theta_{\rm B}$ are the observer's viewing angle and the polar angle of
the magnetic axis measured from the rotation axis. The maximum PF is around 50
%. Also, we find that the PF becomes less than 5 % for $\theta_{\rm obs}
\lesssim 5^\circ$ or for $\theta_{\rm B} \lesssim 5^\circ$. Our results are
consistent with observations of ultraluminous X-ray pulsars (ULXPs) with few
exceptions, since the ULXPs mostly exhibit the PF of $\lesssim$ 50 %. Our
present study supports the hypothesis that the ULXPs are powered by the
super-critical column accretion onto NSs.
[3]
oai:arXiv.org:2002.00954 [pdf] - 2065403
Polarization Imaging of M87 Jets by General Relativistic Radiative
Transfer Calculation based on GRMHD Simulations
Submitted: 2020-02-03, last modified: 2020-02-05
The spectacular images of the M87 black hole taken by the Event Horizon
Telescope (EHT) have opened a new era of black hole research. One of the next
issues is to take polarization images around the central black hole (BH). Since
radio emission is produced by synchrotron process, polarization properties
should vividly reflect the magnetic field structure at the jet base and thus
provide good information regarding the magnetic mechanism of jet formation.
With this kept in mind we perform general relativistic (GR) radiative transfer
calculations of polarized light based on the GR magnetohydrodynamic (MHD)
simulation data of accretion flow and outflow in M87, to obtain their linear
and circular polarization images in the horizon-scale. We found that the linear
polarization components from the jet base and inner accretion flow should
experience Faraday rotation when passing through magnetized plasmas around the
BH, thus sensitively depending on the BH spin. Through the comparison with
intensity image at 1.3mm by EHT and the polarized degree and the rotation
measure (RM) at 1.3mm with the Submillimeter Array, the model with the spin
parameter of a=0.9M_BH (with M_BH being the BH mass) is favored over other
models with a=0.5M_BH or 0.99M_BH, though we need further systematic studies
for confirmation. We also find in low-temperature models clear ring-like image
in the circular polarization map, which arises because of Faraday conversion of
the linearly polarized synchrotron emission and is thus indicative of magnetic
field. This occurs only when the emission region is threaded with well-ordered
magnetic fields and hence no clear images are expected in high-temperature disk
models, in which disk emission is appreciable. We will be able to elucidate the
field configuration through the comparison between the simulated polarization
images and future polarimetry with EHT and other VLBI observations.
[4]
oai:arXiv.org:1909.01411 [pdf] - 1996114
Studying Black Holes on Horizon Scales with VLBI Ground Arrays
Blackburn, Lindy;
Doeleman, Sheperd;
Dexter, Jason;
Gómez, José L.;
Johnson, Michael D.;
Palumbo, Daniel C.;
Weintroub, Jonathan;
Bouman, Katherine L.;
Chael, Andrew A.;
Farah, Joseph R.;
Fish, Vincent;
Loinard, Laurent;
Lonsdale, Colin;
Narayanan, Gopal;
Patel, Nimesh A.;
Pesce, Dominic W.;
Raymond, Alexander;
Tilanus, Remo;
Wielgus, Maciek;
Akiyama, Kazunori;
Bower, Geoffrey;
Broderick, Avery;
Deane, Roger;
Fromm, Christian M.;
Gammie, Charles;
Gold, Roman;
Janssen, Michael;
Kawashima, Tomohisa;
Krichbaum, Thomas;
Marrone, Daniel P.;
Matthews, Lynn D.;
Mizuno, Yosuke;
Rezzolla, Luciano;
Roelofs, Freek;
Ros, Eduardo;
Savolainen, Tuomas K.;
Yuan, Feng;
Zhao, Guangyao
Submitted: 2019-09-03, last modified: 2019-11-08
High-resolution imaging of supermassive black holes is now possible, with new
applications to testing general relativity and horizon-scale accretion and
relativistic jet formation processes. Over the coming decade, the EHT will
propose to add new strategically placed VLBI elements operating at 1.3mm and
0.87mm wavelength. In parallel, development of next-generation backend
instrumentation, coupled with high throughput correlation architectures, will
boost sensitivity, allowing the new stations to be of modest collecting area
while still improving imaging fidelity and angular resolution. The goal of
these efforts is to move from imaging static horizon scale structure to dynamic
reconstructions that capture the processes of accretion and jet launching in
near real time.
[5]
oai:arXiv.org:1909.01405 [pdf] - 1955406
Studying black holes on horizon scales with space-VLBI
Haworth, Kari;
Johnson, Michael D.;
Pesce, Dominic W.;
Palumbo, Daniel C. M.;
Blackburn, Lindy;
Akiyama, Kazunori;
Boroson, Don;
Bouman, Katherine L.;
Farah, Joseph R.;
Fish, Vincent L.;
Honma, Mareki;
Kawashima, Tomohisa;
Kino, Motoki;
Raymond, Alexander;
Silver, Mark;
Weintroub, Jonathan;
Wielgus, Maciek;
Doeleman, Sheperd S.;
Gomez, Jose L.;
Kauffmann, Jens;
Keating, Garrett K.;
Krichbaum, Thomas P.;
Loinard, Laurent;
Narayanan, Gopal;
James, Akihiro Doi David J.;
Marrone, Daniel P.;
Mizuno, Yosuke;
Nagai, Hiroshi
Submitted: 2019-09-03
The Event Horizon Telescope (EHT) recently produced the first horizon-scale
image of a supermassive black hole. Expanding the array to include a 3-meter
space telescope operating at >200 GHz enables mass measurements of many black
holes, movies of black hole accretion flows, and new tests of general
relativity that are impossible from the ground.
[6]
oai:arXiv.org:1611.01775 [pdf] - 1944769
Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications
Matsumoto, Yosuke;
Asahina, Yuta;
Kudoh, Yuki;
Kawashima, Tomohisa;
Matsumoto, Jin;
Takahashi, Hiroyuki R.;
Minoshima, Takashi;
Zenitani, Seiji;
Miyoshi, Takahiro;
Matsumoto, Ryoji
Submitted: 2016-11-06, last modified: 2019-08-18
We present a new magnetohydrodynamic (MHD) simulation code with the aim of
providing accurate numerical solutions to astrophysical phenomena where
discontinuities, shock waves, and turbulence are inherently important. The code
implements the HLLD approximate Riemann solver, the
fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the
hyperbolic divergence cleaning method for a magnetic field. This choice of
schemes significantly improved numerical accuracy and stability, and saved
computational costs in multidimensional problems. Numerical tests of one- and
two-dimensional problems showed the advantages of using the high-order scheme
by comparing with results from a standard second-order TVD MUSCL scheme. The
present code enabled us to explore long-term evolution of a three-dimensional
accretion disk around a black hole, in which compressible MHD turbulence caused
continuous mass accretion via nonlinear growth of the magneto-rotational
instability (MRI). Numerical tests with various computational cell sizes
exhibited a convergent picture of the early nonlinear growth of the MRI in a
global model, and indicated that the MP5 scheme has more than twice the
resolution of the MUSCL scheme in practical applications.
[7]
oai:arXiv.org:1905.10717 [pdf] - 1901970
Black Hole Spin Signature in the Black Hole Shadow of M87 in the Flaring
State
Submitted: 2019-05-25
Imaging the immediate vicinity of supermassive black holes (SMBH) and
extracting a BH-spin signature is one of the grand challenges in astrophysics.
M87 is known as one of the best targets for imaging the BH shadow and it can be
partially thick against synchrotron self-absorption (SSA), particularly in a
flaring state with high mass accretion rate. However, little is known about
influences of the SSA-thick region on BH shadow images. Here we investigate BH
shadow images of M87 at 230 GHz properly taking into account the SSA-thick
region. When the BH has a high spin value, the corresponding BH shadow image
shows the positional offset between the center of the photon ring and that of
the SSA-thick ring at the innermost stable circular orbit (ISCO) due to the
frame-dragging effect in the Kerr spacetime. As a result, we find that a
dark-crescent structure is generally produced between the photon ring and the
SSA-thick ISCO ring in the BH shadow image. The scale size of the dark-crescent
increase with BH spin: its width reaches up to $\sim 2$ gravitational radius
when the BH spin is 99.8% of its maximum value. The dark crescent is regarded
as a new signature of a highly spinning BH. This feature is expected to appear
in flaring states with relatively high mass accretion rate rather than the
quiescent states. We have simulated the image reconstruction of our theoretical
image by assuming the current and future Event Horizon Telescope (EHT) array,
and have found that the future EHT including space-very long baseline
interferometry in 2020s can detect the dark crescent.
[8]
oai:arXiv.org:1809.01151 [pdf] - 1779631
Systematic two-dimensional radiation hydrodynamic simulations of
super-Eddington accretion flow and outflow: Comparison with the slim disk
model
Submitted: 2018-09-04
To what extent can the one-dimensional slim disk model reproduce the
multi-dimensional results of global radiation-hydrodynamic simulations of
super-Eddington accretion? With this question in mind, we perform systematic
simulation study of accretion flow onto a non-spinning black hole for a variety
of black hole masses of $(10 - 10^7) M_{\odot}$ and mass accretion rates of
$(1.4 \times 10^2 - 5.6 \times 10^3)L_{\rm Edd}/c^2$ (with $L_{\rm Edd}$ and
$c$ being the Eddington luminosity and speed of light). In order to adequately
resolve large-scale outflow structure, we extensively expand a simulation box
to cover the space of $3000 r_{\rm S}$ (with $r_{\rm S}$ being the
Schwarzschild radius), larger than those in most previous studies, so that we
can put relatively large angular momentum to the gas injected from the outer
simulation boundary. The adopted Keplerian radius, at which the centrifugal
force balances with the gravitational force, is $r_{\rm K}=300 r_{\rm S}$. The
injected mass first falls and is accumulated at around this radius and then
slowly accretes towards the central black hole via viscosity. We simulate such
accretion processes, taking inverse and bulk Compton scattering into account.
The simulated accretion flow is in a quasi-steady state inside $r_{\rm qss}\sim
200 r_{\rm S}$. Within this radius the flow properties are, as a whole, in good
agreement with those described by the slim disk model except that the radial
density profile of the underlying disk is much flatter, $\rho\propto r^{-0.73}$
(cf. $\rho\propto r^{-3/2}$ in the slim disk model), due probably to efficient
convection. We find very weak outflow from inside $r\sim 200 r_{\rm S}$ unlike
the past studies.
[9]
oai:arXiv.org:1802.00567 [pdf] - 1641480
Three-dimensional structure of clumpy outflow from supercritical
accretion flow onto black holes
Submitted: 2018-02-02
We perform global three-dimensional (3D) radiation-hydrodynamic (RHD)
simulations of out- flow from supercritical accretion flow around a 10 Msun
black hole. We only solve the outflow part, starting from the axisymmetric 2D
simulation data in a nearly steady state but with small perturbations in a
sinusoidal form being added in the azimuthal direction. The mass accretion rate
onto the black hole is ~10^2 L_E/c^2 in the underlying 2D simulation data and
the outflow rate is ~10 L_E/c^2 (with LE and c being the Eddington luminosity
and speed of light, respectively). We first confirm the emergence of clumpy
outflow, which was discovered by the 2D RHD simulations, above the photosphere
located at a few hundreds of Schwarzschild radii (r_S) from the central black
hole. As prominent 3D features we find that the clumps have the shape of a torn
sheet, rather than a cut string, and that they are rotating around the central
black hole with a sub-Keplerian velocity at a distance of ~10^3 r_S from the
center. The typical clump size is ~30 r_S or less in the radial direction, and
is more elongated in the angular directions, ~hundreds of r_S at most. The
sheet separation ranges from 50 to 150 r_S. We expect stochastic time
variations when clumps pass across the line of the sight of a distant observer.
Variation timescales are estimated to be several seconds for a black hole with
mass of ten to several tens of Msun, in rough agreement with the observations
of some ultra-luminous X-ray sources.
[10]
oai:arXiv.org:1709.01531 [pdf] - 1604894
Theoretical modeling of Comptonized X-ray spectra of super-Eddington
accretion flow: origin of hard excess in Ultraluminous X-Ray Sources
Submitted: 2017-09-05
X-ray continuum spectra of super-Eddington accretion flow are studied by
means of Monte Carlo radiative transfer simulations based on the radiation
hydrodynamic simulation data, in which both of thermal and bulk Compton
scatterings are taken into account. We compare the calculated spectra of
accretion flow around black holes with masses of $M_{\rm BH} = 10, 10^2, 10^3$,
and $10^4 M_\odot$ for a fixed mass injection rate (from the computational
boundary at $10^3 r_{\rm s}$) of $10^3 L_{\rm Edd}/c^2$ (with $r_{\rm s}$,
$L_{\rm Edd}$, and $c$ being the Schwarzschild radius, the Eddington
luminosity, and the speed of light, respectively). The soft X-ray spectra
exhibit mass dependence in accordance with the standard-disk relation; the
maximum surface temperature is scaled as $T \propto M_{\rm BH}^{-1/4}$. The
spectra in the hard X-ray bands, by contrast, look quite similar among
different models, if we normalize the radiation luminosity by $M_{\rm BH}$.
This reflects that the hard component is created by thermal and bulk Compton
scattering of soft photons originating from an accretion flow in the
over-heated and/or funnel regions, the temperatures of which have no mass
dependence. The hard X-ray spectra can be reproduced by a Wien spectrum with
temperature of $T\sim 3$ keV accompanied by a hard excess at photon energy
above several keV. The excess spectrum can be well fitted with a power law with
a photon index of $\Gamma \sim 3$. This feature is in good agreement with that
of the recent NuSTAR observations of ULX (Ultra-Luminous X-ray sources).
[11]
oai:arXiv.org:1705.04776 [pdf] - 1577577
White Paper on East Asian Vision for mm/submm VLBI: Toward Black Hole
Astrophysics down to Angular Resolution of 1~R$_{S}$
Asada, K.;
Kino, M.;
Honma, M.;
Hirota, T.;
Lu, R. -S.;
Inoue, M.;
Sohn, B. -W.;
Shen, Z. -Q.;
Ho, P. T. P.;
Akiyama, K.;
Algaba, J-C.;
An, T.;
Bower, G.;
Byun, D-Y.;
Dodson, R.;
Doi, A.;
Edwards, P. G.;
Fujisawa, K.;
Gu, M-F.;
Hada, K.;
Hagiwara, Y.;
Jaroenjittichai, P.;
Jung, T.;
Kawashima, T.;
Koyama, S.;
Lee, S-S.;
Matsushita, S.;
Nagai, H.;
Nakamura, M.;
Niinuma, K.;
Phillips, C.;
Park, J-H.;
Pu, H-Y.;
Ro, H-W.;
Stevens, J.;
Trippe, S.;
Wajima, K.;
Zhao, G-Y.
Submitted: 2017-05-12
This White Paper details the intentions and plans of the East Asian Very Long
Baseline Interferometry (VLBI) community for pushing the frontiers of
millimeter/submillimeter VLBI. To this end, we shall endeavor to actively
promote coordinated efforts in the East Asia region. Our goal is to establish
firm collaborations among the East Asia VLBI community in partnership with
related institutes in North America and Europe and to expand existing global
mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an
ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b)
investigating the dynamics of circumstellar gas in star-forming regions and
late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In
the first half of this White Paper, we highlight scientific accomplishments of
the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr
A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are
described, and future visions of our VLBI science are briefly presented. Maser
science of star formation, stellar evolution, and physics of accretion disks
around AGNs are also discussed. A new vision for conducting multi-transition
maser studies using mm/submm VLBI together with the Atacama Large
Millimeter/Submillimeter Array (ALMA) is described. In the second half of this
White Paper, we describe the EA community's vision for using mm/submm VLBI
arrays in the framework or extended version of the Event Horizon Telescope
(EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the
aforementioned goal will maximize the overall scientific outcomes of mm/submm
VLBI in the world.
[12]
oai:arXiv.org:1702.07903 [pdf] - 1560350
A Possible Time-Delayed Brightening of the Sgr A* Accretion Flow after
the Pericenter Passage of G2 Cloud
Submitted: 2017-02-25, last modified: 2017-04-10
A possibility of time-delayed radio brightenings of Sgr A* triggered by the
pericenter passage of the G2 cloud is studied by carrying out global
three-dimensional magnetohydrodynamic simulations taking into account the
radiative cooling of the tidal debris of the G2 cloud. Magnetic fields in the
accretion flow are strongly perturbed and re-organized after the passage of G2.
We have found that the magnetic energy in the accretion flow increases by a
factor 3-4 in 5-10 years after the pericenter passage of G2 by a dynamo
mechanism driven by the magneto-rotational instability. Since this B-field
amplification enhances the synchrotron emission from the disk and the outflow,
the radio and the infrared luminosity of Sgr A* is expected to increase around
A.D. 2020. The time-delay of the radio brightening enables us to determine the
rotation axis of the preexisting disk.
[13]
oai:arXiv.org:1702.07519 [pdf] - 1537934
Magnetohydrodynamic Simulations of the Formation of Molecular Clouds
toward the Stellar Cluster Westerlund 2: Interaction of a Jet with Clumpy
Interstellar Medium
Submitted: 2017-02-24
The formation mechanism of CO clouds observed with NANTEN2 and Mopra
telescope toward the stellar cluster Westerlund 2 is studied by
three-dimensional magnetohydrodynamic (MHD) simulations taking into account the
interstellar cooling. These molecular clouds show a peculiar shape composing of
an arc-shaped cloud in one side of a TeV-ray source HESS J1023-575 and a linear
distribution of clouds (jet clouds) in another side. We propose that these
clouds are formed by the interaction of a jet with interstellar neutral
hydrogen (HI) clumps. By studying the dependence of the shape of dense cold
clouds formed by shock compression and cooling on the filling factor of HI
clumps, we found that the density distribution of HI clumps determines the
shape of molecular clouds formed by the jet-cloud interaction; arc-clouds are
formed when the filling factor is large. On the other hand, when the filling
factor is small, molecular clouds align with the jet. The jet propagates faster
in models with small filling factors.
[14]
oai:arXiv.org:1608.04211 [pdf] - 1477665
A radiation-hydrodynamic model of accretion columns for ultra-luminous
X-ray pulsars
Submitted: 2016-08-15, last modified: 2016-08-29
Prompted by the recent discovery of pulsed emission from an ultra-luminous
X-ray source, M82 X-2 ("ULX-pulsar"), we perform a two-dimensional
radiation-hydrodynamic simulation of a super-critical accretion flow onto a
neutron star through a narrow accretion column. We set an accretion column with
a cone shape filled with tenuous gas with density of $10^{-4} {\rm g}~ {\rm
cm}^{-3}$ above a neutron star and solve the two dimensional gas motion and
radiative transfer within the column. The side boundaries are set such that
radiation can freely escape, while gas cannot. Since the initial gas layer is
not in a hydrostatic balance, the column gas falls onto the neutron-star
surface, thereby a shock being generated. As a result, the accretion column is
composed of two regions: an upper, nearly free-fall region and a lower settling
region, as was noted by Basko \& Sunyaev (1976). The average accretion rate is
very high; ${\dot M}\sim 10^{2-3} L_{\rm E}/c^2$ (with $L_{\rm E}$ being the
Eddington luminosity), and so radiation energy dominates over gas internal
energy entirely within the column. Despite the high accretion rate, the
radiation flux in the laboratory frame is kept barely below $L_{\rm E}/(4\pi
r^2)$ at a distance $r$ in the settling region so that matter can slowly
accrete. This adjustment is made possible, since large amount of photons
produced via dissipation of kinetic energy of matter can escape through the
side boundaries. The total luminosity can greatly exceed $L_{\rm E}$ by several
orders of magnitude, whereas the apparent luminosity observed from the top of
the column is much less. Due to such highly anisotropic radiation fields,
observed flux should exhibit periodic variations with the rotation period,
provided that the rotation and magnetic axes are misaligned.
[15]
oai:arXiv.org:1607.02443 [pdf] - 1580049
TARGET 5: a new multi-channel digitizer with triggering capabilities for
gamma-ray atmospheric Cherenkov telescopes
Albert, A.;
Funk, S.;
Kawashima, T.;
Murphy, M.;
Okumura, A.;
Quagliani, R.;
Sapozhnikov, L.;
Tajima, H.;
Tibaldo, L.;
Vandenbroucke, J.;
Varner, G.;
Wu, T.
Submitted: 2016-07-08, last modified: 2016-07-11
TARGET~5 is a new application-specific integrated circuit (ASIC) of the
TARGET family, designed for the readout of signals from photosensors in the
cameras of imaging atmospheric Cherenkov telescopes (IACTs) for ground-based
gamma-ray astronomy. TARGET~5 combines sampling and digitization on 16 signal
channels with the formation of trigger signals based on the analog sum of
groups of four channels. We describe the ASIC architecture and performance.
TARGET~5 improves over the performance of the first-generation TARGET ASIC,
achieving: tunable sampling frequency from $<0.4$~GSa/s to $>1$~GSa/s; a
dynamic range on the data path of 1.2 V with effective dynamic range of 11~bits
and DC noise of ${\sim}0.6$~mV; 3-dB bandwidth of 500 MHz; crosstalk between
adjacent channels $<1.3\%$; charge resolution improving from 40\% to $<4\%$
between 3 photoelectrons (p.e.) and $>100$~p.e. (assuming 4 mV per p.e.); and
minimum stable trigger threshold of 20 mV (5 p.e.) with trigger noise of 5 mV
(1.2 p.e.), which is mostly limited by interference between trigger and
sampling operations. TARGET~5 is the first ASIC of the TARGET family used in an
IACT prototype, providing one development path for readout electronics in the
forthcoming Cherenkov Telescope Array (CTA).
[16]
oai:arXiv.org:1605.04992 [pdf] - 1443976
Formation of Overheated Regions and Truncated Disks around Black Holes;
Three-dimensional General Relativistic Radiation-magnetohydrodynamics
Simulations
Submitted: 2016-05-16
Using three-dimensional general relativistic radiation magnetohydrodynamics
simulations of accretion flows around stellar mass black holes, we report that
the relatively cold disk ($\gtrsim 10^{7}$K) is truncated near the black hole.
Hot and less-dense regions, of which the gas temperature is $ \gtrsim 10^9$K
and more than ten times higher than the radiation temperature (overheated
regions), appear within the truncation radius. The overheated regions also
appear above as well as below the disk, and sandwich the cold disk, leading to
the effective Compton upscattering. The truncation radius is $\sim 30 r_{\rm
g}$ for $\dot{M} \sim L_{\rm Edd}/c^2$, where $r_{\rm g}, \dot M,
L_\mathrm{Edd}, c$ are the gravitational radius, mass accretion rate, Eddington
luminosity, and light speed. Our results are consistent with observations of
very high state, whereby the truncated disk is thought to be embedded in the
hot rarefied regions. The truncation radius shifts inward to $\sim 10 r_{\rm
g}$ with increasing mass accretion rate $\dot{M} \sim 100 L_{\rm Edd}/c^2$,
which is very close to an innermost stable circular orbit. This model
corresponds to the slim disk state observed in ultra luminous X-ray sources.
Although the overheated regions shrink if the Compton cooling effectively
reduces the gas temperature, the sandwich-structure does not disappear at the
range of $\dot{M} \lesssim 100L_{\rm Edd}/c^2$. Our simulations also reveal
that the gas temperature in the overheated regions depends on black hole spin,
which would be due to efficient energy transport from black hole to disks
through the Poynting flux, resulting gas heating.
[17]
oai:arXiv.org:1508.06296 [pdf] - 1938070
TARGET: toward a solution for the readout electronics of the Cherenkov
Telescope Array
Tibaldo, L.;
Vandenbroucke, J. A.;
Albert, A. M.;
Funk, S.;
Kawashima, T.;
Kraus, M.;
Okumura, A.;
Sapozhnikov, L.;
Tajima, H.;
Varner, G. S.;
Wu, T.;
Zink, A.
Submitted: 2015-08-25, last modified: 2015-12-18
TARGET is an application specific integrated circuit (ASIC) designed to read
out signals recorded by the photosensors in cameras of very-high-energy
gamma-ray telescopes exploiting the imaging of Cherenkov radiation from
atmospheric showers. TARGET capabilities include sampling at a high rate
(typically 1 GSample/s), digitization, and triggering on the sum of four
adjacent pixels. The small size, large number of channels read out per ASIC
(16), low cost per channel, and deep buffer for trigger latency (~16 $\mu$s at
1 GSample/s) make TARGET ideally suited for the readout in systems with a large
number of telescopes instrumented with compact photosensors like multi-anode or
silicon photomultipliers combined with dual-mirror optics. The possible
advantages of such systems are better sensitivity, a larger field of view, and
improved angular resolution. The two latest generations of TARGET ASICs, TARGET
5 and TARGET 7, are soon to be used for the first time in two prototypes of
small-sized and medium-sized dual-mirror telescopes proposed in the framework
of the Cherenkov Telescope Array (CTA) project. In this contribution we report
on the performance of the TARGET ASICs and discuss future developments.
[19]
oai:arXiv.org:1412.4521 [pdf] - 1223640
Radiation hydrodynamics simulations of wide-angle outflows from
super-critical accretion disks around black holes
Submitted: 2014-12-15
By performing two-dimensional radiation hydrodynamics simulations with large
computational domain of 5000 Schwarzschild radius, we revealed that wide-angle
outflow is launched via the radiation force from the super-critical accretion
flows around black holes. The angular size of the outflow, of which the radial
velocity (v_r) is over the escape velocity (v_esc), increases with an increase
of the distance from the black hole. As a result, the mass is blown away with
speed of v_r > v_esc in all direction except for the very vicinity of the
equatorial plane, theta=0-85^circ, where theta is the polar angle. The mass
ejected from the outer boundary per unit time by the outflow is larger than the
mass accretion rate onto the black hole, ~150L_Edd/c^2, where L_Edd and c are
the Eddington luminosity and the speed of light. Kinetic power of such
wide-angle high-velocity outflow is comparable to the photon luminosity and is
a few times larger than the Eddington luminosity. This corresponds to
~10^39-10^40 erg/s for the stellar mass black holes. Our model consistent with
the observations of shock excited bubbles observed in some ultra-luminous X-ray
sources (ULXs), supporting a hypothesis that ULXs are powered by the
super-critical accretion onto stellar mass black holes.
[20]
oai:arXiv.org:1407.2381 [pdf] - 1215543
Magneto-hydrodynamic Simulations of a Jet Drilling an HI Cloud: Shock
Induced Formation of Molecular Clouds and Jet Breakup
Submitted: 2014-07-09
The formation mechanism of the jet-aligned CO clouds found by NANTEN CO
observations is studied by magnetohydrodynamical (MHD) simulations taking into
account the cooling of the interstellar medium. Motivated by the association of
the CO clouds with the enhancement of HI gas density, we carried out MHD
simulations of the propagation of a supersonic jet injected into the dense HI
gas. We found that the HI gas compressed by the bow shock ahead of the jet is
cooled down by growth of the cooling instability triggered by the density
enhancement. As a result, cold dense sheath is formed around the interface
between the jet and the HI gas. The radial speed of the cold, dense gas in the
sheath is a few km/s almost independent of the jet speed. Molecular clouds can
be formed in this region. Since the dense sheath wrapping the jet reflects
waves generated in the cocoon, the jet is strongly perturbed by the vortices of
the warm gas in the cocoon, which breaks up the jet and forms a secondary shock
in the HI-cavity drilled by the jet. The particle acceleration at the shock can
be the origin of radio and X-ray filaments observed near the eastern edge of
W50 nebula surrounding the galactic jet source SS433.
[21]
oai:arXiv.org:1305.4943 [pdf] - 1171495
Recurrent Outbursts and Jet Ejections Expected in Swift J1644+57:
Limit-Cycle Activities in a Supermassive Black Hole
Submitted: 2013-05-21
The tidal disruption event by a supermassive black hole in Swift J1644+57 can
trigger limit-cycle oscillations between a supercritically accreting X-ray
bright state and a subcritically accreting X-ray dim state. Time evolution of
the debris gas around a black hole with mass $M=10^{6} {\MO}$ is studied by
performing axisymmetric, two-dimensional radiation hydrodynamic simulations. We
assumed the $\alpha$-prescription of viscosity, in which the viscous stress is
proportional to the total pressure. The mass supply rate from the outer
boundary is assumed to be ${\dot M}_{\rm supply}=100L_{\rm Edd}/c^2$, where
$L_{\rm Edd}$ is the Eddington luminosity, and $c$ is the light speed. Since
the mass accretion rate decreases inward by outflows driven by radiation
pressure, the state transition from a supercritically accreting slim disk state
to a subcritically accreting Shakura-Sunyaev disk starts from the inner disk
and propagates outward in a timescale of a day. The sudden drop of the X-ray
flux observed in Swift J1644+57 in August 2012 can be explained by this
transition. As long as ${\dot M}_{\rm supply}$ exceeds the threshold for the
existence of a radiation pressure dominant disk, accumulation of the accreting
gas in the subcritically accreting region triggers the transition from a gas
pressure dominant Shakura-Sunyaev disk to a slim disk. This transition takes
place at $t {\sim}~50/({\alpha}/0.1)$ days after the X-ray darkening. We expect
that if $\alpha > 0.01$, X-ray emission with luminosity $\gtrsim 10^{44}$ ${\rm
erg}{\cdot}{\rm s}^{-1}$ and jet ejection will revive in Swift J1644+57 in
2013--2014.
[22]
oai:arXiv.org:0904.4123 [pdf] - 1002032
New Spectral State of Supercritical Accretion Flow with Comptonizing
Outflow
Submitted: 2009-04-27
Supercritical accretion flows inevitably produce radiation-pressure driven
outflows, which will Compton up-scatter soft photons from the underlying
accretion flow, thereby making hard emission. We perform two dimensional
radiation hydrodynamic simulations of supercritical accretion flows and
outflows, incorporating such Compton scattering effects, and demonstrate that
there appears a new hard spectral state at higher photon luminosities than that
of the slim-disk state. In this state, as the photon luminosity increases, the
photon index decreases and the fraction of the hard emission increases. The
Compton $y$-parameter is of the order of unity (and thus the photon index will
be $\sim 2$) when the apparent photon luminosity is ${\sim}30L_{\rm E}$ (with
$L_{\rm E}$ being the Eddington luminosity) for nearly face-on sources. This
explains the observed spectral hardening of the ULX NGC1313 X-2 in its
brightening phase and thus supports the model of supercritical accretion onto
stellar mass black holes in this ULX.
