Source author record

Yosuke Mizuno

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

39 published item(s)

preprint2022arXiv

Fiber-optic multimode interference sensing: comprehensive characterization and its potential for strain-insensitive temperature sensing

A strain-insensitive temperature sensor based on multimode interference using standard multimode fibers (MMFs) is proposed according to the comprehensive study of the characteristics of the MMFs. The temperature and strain dependences on the core diameter, numerical aperture (NA), and the length of the MMF section in the single-mode--multimode--single-mode (SMS) fiber structure are investigated experimentally. The results indicate that the larger core diameter of the MMF leads to higher temperature sensitivity but lower strain sensitivity (absolute values); the higher NA does not influence the temperature sensitivity but results in higher absolute value of strain sensitivity; the longer MMF section brings lower temperature sensitivity but does not have an impact on strain sensitivity. These findings also contribute to the theoretical analysis of the length dependence in the SMS fiber sensors. Besides, the results of the characterization study show that the strain sensitivity is relatively low, which brings a possibility to develop a strain-insensitive temperature sensor. The proposed sensor is used for temperature sensing while the strain is constantly applied from 0 to 1100 $με$ with steps of 100 $με$. The measured results are consistent with the comprehensive study. The mean temperature sensitivity is 6.14 pm/$^{\circ}$C with a standard deviation of 0.39 pm/$^{\circ}$C, which proves that the proposed temperature sensor exhibits good stability and is insensitive to strain. We expect that these results will provide a profound guideline to fiber sensors based on multimode interference.

preprint2022arXiv

GRMHD Simulations and Modeling for Jet Formation and Acceleration Region in AGNs

Relativistic jets are collimated plasma outflows with relativistic speeds. Astrophysical objects involving relativistic jets are a system comprising a compact object such as a black hole, surrounded by rotating accretion flows, with the relativistic jets produced near the central compact object. The most accepted models explaining the origin of relativistic jets involve magnetohydrodynamic (MHD) processes. Over the past few decades, many general relativistic MHD (GRMHD) codes have been developed and applied to model relativistic jet formation in various conditions. This short review provides an overview of the recent progress of GRMHD simulations in generating relativistic jets and their modeling for observations.

preprint2022arXiv

Particle acceleration in shearing flows: the self-generation of turbulent spine-sheath structures in relativistic MHD jet simulations

X-ray observations of several kiloparsec-scale extragalactic jets favour a synchrotron origin. The short cooling times of the emitting electrons requires distributed acceleration of electrons up to sub-PeV energies. In a previous paper, we found that this can be self-consistently explained by a shear acceleration model, where particles are accelerated to produce power-law spectra with a spectral index being determined mainly by the velocity profile and turbulence spectrum. In this paper, we perform 3D relativistic magneto-hydrodynamic simulations to investigate the formation of a spine-sheath structure and the development of turbulence for a relativistic jet propagating into a static cocoon. We explore different spine velocities and magnetic field profiles with values being chosen to match typical Fanaroff-Riley type I/II jets. We find that in all cases a sheath is generated on the interface of the spine and the cocoon mainly due to the Kelvin-Helmholtz instability. The large scale velocity profile in the sheath is close to linear. Turbulence develops in both the spine and the sheath, with a turbulent velocity spectrum consistent with Kolmogorov-scaling. The implications for shear particle acceleration are explored, with a focus on the particle spectral index.

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.

preprint2022arXiv

Stringent axion constraints with Event Horizon Telescope polarimetric measurements of M87$^\star$

The hitherto unprecedented angular resolution of the Event Horizon Telescope (EHT) has created exciting opportunities in the search for new physics. Recently, the linear polarization of radiation emitted near the supermassive black hole M87$^\star$ was measured on four separate days, precisely enabling tests of the existence of a dense axion cloud produced by a spinning black hole. The presence of an axion cloud leads to a frequency-independent oscillation in the electric vector position angle (EVPA) of this linear polarization. For a nearly face-on M87$^\star$, this oscillation in the EVPA appears as a propagating wave along the photon ring. In this paper, we leverage the azimuthal distribution of EVPA measured by the EHT to study the axion-photon coupling. We propose a novel differential analysis procedure to reduce the astrophysical background, and derive stringent constraints on the existence of axions in the previously unexplored mass window $\sim (10^{-21}-10^{-20})$~eV.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2022arXiv

The science case and challenges of space-borne sub-millimeter interferometry

Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies.

preprint2022arXiv

Unravelling the Innermost Jet Structure of OJ 287 with the First GMVA+ALMA Observations

We present the first very-long-baseline interferometric (VLBI) observations of the blazar OJ287 carried out jointly with the Global Millimeter VLBI Array (GMVA) and the phased Atacama Large Millimeter/submillimeter Array (ALMA) at 3.5 mm on April 2, 2017. Participation of phased-ALMA not only has improved the GMVA north-south resolution by a factor of ~3, but also has enabled fringe detection with signal-to-noise ratios up to 300 at baselines longer than 2 Gλ. The high sensitivity has motivated us to image the data with the newly developed regularized maximum likelihood imaging methods, revealing the innermost jet structure with unprecedentedly high angular resolution. Our images reveal a compact and twisted jet extending along the northwest direction with two bends within the inner 200 μas that resembles a precessing jet in projection. The component at the southeastern end shows a compact morphology and high brightness temperature, and is identified as the VLBI core. An extended jet feature that lies at ~200 μas northwest of the core shows a conical shape in both total and linearly polarized intensity, and a bimodal distribution of the linear polarization electric vector position angle. We discuss the nature of this feature by comparing our observations with models and simulations of oblique and recollimation shocks with various magnetic field configurations. Our high-fidelity images also enabled us to search for possible jet features from the secondary supermassive black hole (SMBH) and test the SMBH binary hypothesis proposed for this source.

preprint2021arXiv

Impact of non-thermal particles on the spectral and structural properties of M87

The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations. The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87

preprint2020arXiv

Deep Horizon; a machine learning network that recovers accreting black hole parameters

The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from images of black hole shadows. We investigate the effects of a limited telescope resolution and observations at higher frequencies. We trained two convolutional deep neural networks on a large image library of simulated mock data. The first network is a Bayesian deep neural regression network and is used to recover the viewing angle $i$, and position angle, mass accretion rate $\dot{M}$, electron heating prescription $R_{\rm high}$ and the black hole mass $M_{\rm BH}$. The second network is a classification network that recovers the black hole spin $a$. We find that with the current resolution of the Event Horizon Telescope, it is only possible to accurately recover a limited number of parameters of a static image, namely the mass and mass accretion rate. Since potential future space-based observing missions will operate at frequencies above 230 GHz, we also investigated the applicability of our network at a frequency of 690 GHz. The expected resolution of space-based missions is higher than the current resolution of the Event Horizon Telescope, and we show that Deep Horizon can accurately recover the parameters of simulated observations with a comparable resolution to such missions.

preprint2020arXiv

How to tell an accreting boson star from a black hole

The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent with the properties of the accretion flow onto Sgr A*. To this end, we perform general relativistic magnetohydrodynamic simulations followed by general relativistic radiative transfer calculations in the boson star space-time. Synthetic reconstructed images considering realistic astronomical observing conditions show that, despite qualitative similarities, the differences in the appearance of a black hole -- either rotating or not -- and a boson star of the type considered here are large enough to be detectable. These differences arise from dynamical effects directly related to the absence of an event horizon, in particular, the accumulation of matter in the form of a small torus or a spheroidal cloud in the interior of the boson star, and the absence of an evacuated high-magnetization funnel in the polar regions. The mechanism behind these effects is general enough to apply to other horizonless and surfaceless black hole mimickers, strengthening confidence in the ability of the EHT to identify such objects via radio observations.

preprint2020arXiv

Particle acceleration by relativistic magnetic reconnection driven by kink instability turbulence in Poynting flux dominated jets

Particle acceleration in magnetized relativistic jets still puzzles theorists, specially when one tries to explain the highly variable emission observed in blazar jets or gamma-ray bursts putting severe constraints on current models. In this work we investigate the acceleration of particles injected in a three-dimensional relativistic magnetohydrodynamical jet subject to current driven kink instability (CDKI), which drives turbulence and fast magnetic reconnection. Test protons injected in the nearly stationary snapshots of the jet, experience an exponential acceleration up to a maximum energy. For a background magnetic field of $B \sim 0.1$ G, this saturation energy is $\sim 10^{16}$ eV, while for $B \sim 10$ G it is $\sim 10^{18}$ eV. The simulations also reveal a clear association of the accelerated particles with the regions of fast reconnection. In the early stages of the development of the non-linear growth of CDKI in the jet, when there are still no sites of fast reconnection, injected particles are also efficiently accelerated, but by magnetic curvature drift in the wiggling jet spine. However, they have to be injected with an initial energy much larger than that required for particles to accelerate in reconnection sites. Finally, we have also obtained from the simulations an acceleration time due to reconnection with a weak dependence on the particles energy $E$, $t_A \propto E^{0.1}$. The energy spectrum of the accelerated particles develops a high energy tail with a power law index $p \sim$ -1.2 in the beginning of the acceleration, in agreement with earlier works. Our results provide an appropriate multi-dimensional framework for exploring this process in real systems and explain their complex emission patterns, specially in the very high energy bands and the associated neutrino emission recently detected in some blazars.

preprint2020arXiv

Plasmoid formation in global GRMHD simulations and AGN flares

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

preprint2020arXiv

Rapid Particle Acceleration due to Recollimation Shocks and Turbulent Magnetic Fields in Injected Jets with Helical Magnetic Fields

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

preprint2019arXiv

Modelling the polarised emission from black holes on event horizon-scales

Upcoming VLBI observations will resolve nearby supermassive black holes, most notably Sagittarius A* and M87, on event horizon-scales. Recent observations of Sagittarius A* with the Event Horizon Telescope have revealed horizon-scale structure. Accordingly, the detection and measurement of the back hole ``shadow" is expected to enable the existence of astrophysical black holes to be verified directly. Although the theoretical description of the shadow is straightforward, its observational appearance is largely determined by the properties of the surrounding accretion flow, which is highly turbulent. We introduce a new polarised general-relativistic radiative transfer code, \texttt{BHOSS}, which accurately solves the equations of polarised radiative transfer in arbitrary strong-gravity environments, providing physically-realistic images of astrophysical black holes on event horizon-scales, as well as also providing insight into the fundamental properties and nature of the surrounding accretion flow environment.

preprint2018arXiv

Particle acceleration and the origin of the very high energy emission around black holes and relativistic jets

Particle acceleration induced by fast magnetic reconnection may help to solve current puzzles related to the interpretation of the very high energy (VHE) and neutrino missions from AGNs and compact sources in general. Our general relativistic-MHD simulations of accretion disk-corona systems reveal the growth of turbulence driven by MHD instabilities that lead to the development of fast magnetic reconnection in the corona. In addition, our simulations of relativistic MHD jets reveal the formation of several sites of fast reconnection induced by current-driven kink turbulence. The injection of thousands of test particles in these regions cause acceleration up to energies of several PeVs, thus demonstrating the ability of this process to accelerate particles and produce VHE and neutrino emission, specially in blazars. Finally, we discuss how reconnection can also explain the observed VHE luminosity-black hole mass correlation, involving hundreds of non-blazar sources like Perseus A, and black hole binaries.

preprint2016arXiv

A new method for shadow calculations: application to parameterised axisymmetric black holes

Collaborative international efforts under the name of the Event Horizon Telescope project, using sub- mm very long baseline interferometry, are soon expected to provide the first images of the shadow cast by the candidate supermassive black hole in our Galactic center, Sagittarius A*. Observations of this shadow would provide direct evidence of the existence of astrophysical black holes. Although it is expected that astrophysical black holes are described by the axisymmetric Kerr solution, there also exist many other black hole solutions, both in general relativity and in other theories of gravity, which cannot presently be ruled out. To this end, we present calculations of black hole shadow images from various metric theories of gravity as described by our recent work on a general parameterisation of axisymmetric black holes [R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016)]. An algorithm to perform general ray-tracing calculations for any metric theory of gravity is first outlined and then employed to demonstrate that even for extremal metric deformation parameters of various black hole spacetimes, this parameterisation is both robust and rapidly convergent to the correct solution.

preprint2016arXiv

Magnetic Reconnection on Jet-Accretion disk Systems

Fast Magnetic Reconnection is currently regarded as an important process also beyond the solar system, specially in magnetically dominated regions of galactic and extragalactic sources like the surrounds of black holes and relativistic jets. In this lecture we discuss briefly the theory of fast magnetic reconnection, specially when driven by turbulence which is very frequent in Astrophysical flows, and its implications for relativistic particle acceleration. Then we discuss these processes in the context of the sources above, showing recent analytical and multidimensional numerical MHD studies that indicate that fast reconnection can be a powerful process to accelerate particles to relativistic velocities, produce the associated high energy non-thermal emission, and account for efficient conversion of magnetic into kinetic energy in these flows.

preprint2016arXiv

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

preprint2016arXiv

Spatial growth of current-driven instability in relativistic rotating jets and the search for magnetic reconnection

Using the three-dimensional relativistic magnetohydrodynamic code RAISHIN, we investigated the influence of radial density profile on the spatial development of the current-driven kink instability along magnetized rotating, relativistic jets. For the purpose of our study, we used a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers the growth of the kink instability. We studied light as well as heavy jets with respect to the environment depending on the density profile. Different angular velocity amplitudes have been also tested. The results show the propagation of a helically kinked structure along the jet and relatively stable configuration for the lighter jets. The jets appear to be collimated by the magnetic field and the flow is accelerated due to conversion of electromagnetic into kinetic energy. We also identify regions of high current density in filamentary current sheets, indicative of magnetic reconnection, which are associated to the kink unstable regions and correlated to the decrease of the sigma parameter of the flow. We discuss the implications of our findings for Poynting-flux dominated jets in connection with magnetic reconnection process. We find that fast magnetic reconnection may be driven by the kink-instability turbulence and govern the transformation of magnetic into kinetic energy thus providing an efficient way to power and accelerate particles in AGN and gamma-ray-burst relativistic jets.

preprint2015arXiv

High-speed optical correlation-domain reflectometry without using acousto-optic modulator

To achieve a distributed reflectivity measurement along an optical fiber, we develop a simplified cost-effective configuration of optical correlation- (or coherence-) domain reflectometry based on a synthesized optical coherence function by sinusoidal modulation. By excluding conventional optical heterodyne detection (practically, without using an acousto-optic modulator) and by exploiting the foot of the Fresnel reflection spectrum, the electrical bandwidth required for signal processing is lowered down to several megahertz. We evaluate the basic system performance and demonstrate its high-speed operation (10 ms for one scan) by tracking a moving reflection point in real time.

preprint2015arXiv

Jet Motion, Internal Working Surfaces, and Nested Shells in the Protostellar System HH 212

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.

preprint2015arXiv

Probing the innermost regions of AGN jets and their magnetic fields with RadioAstron. I. Imaging BL Lacertae at 21 microarcsecond resolution

We present the first polarimetric space VLBI imaging observations at 22 GHz. BL Lacertae was observed in 2013 November 10 with the RadioAstron space VLBI mission, including a ground array of 15 radio telescopes. The instrumental polarization of the space radio telescope is found to be within 9%, demonstrating the polarimetric imaging capabilities of RadioAstron at 22 GHz. Ground-space fringes were obtained up to a projected baseline distance of 7.9 Earth's diameters in length, allowing us to image the jet in BL Lacertae with a maximum angular resolution of 21 $μ$as, the highest achieved to date. We find evidence for emission upstream of the radio core, which may correspond to a recollimation shock at about 40 $μ$as from the jet apex, in a pattern that includes other recollimation shocks at approximately 100 $μ$as and 250 $μ$as from the jet apex. Polarized emission is detected in two components within the innermost 0.5 mas from the core, as well as in some knots 3 mas downstream. Faraday rotation analysis, obtained from combining RadioAstron 22 GHz and ground-based 15 GHz and 43 GHz images, shows a gradient in rotation measure and Faraday corrected polarization vector as a function of position angle with respect to the core, suggesting that the jet in BL Lacertae is threaded by a helical magnetic field. The intrinsic de-boosted brightness temperature in the unresolved core exceeds $3\!\times\!10^{12}$ K, suggesting at the very least departure from equipartition of energy between the magnetic field and radiating particles.

preprint2015arXiv

Recollimation Shocks in Magnetized Relativistic Jets

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

preprint2015arXiv

Steady General Relativistic Magnetohydrodynamic Inflow/Outflow Solution along Large-Scale Magnetic Fields that Thread a Rotating Black Hole

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

preprint2015arXiv

Studies of Relativistic Jets in Active Galactic Nuclei with SKA

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

preprint2014arXiv

Alternative implementation of simplified Brillouin optical correlation-domain reflectometry

We developed an alternative configuration of simplified Brillouin optical correlation-domain reflectometry, which can overcome the drawbacks of the original configuration. This system uses, as reference light, the light that is Fresnel reflected at a partial reflection point artificially produced near an optical circulator. We show that the influence of the 0th correlation peak fixed at the partial reflection point can be suppressed by replacing the nearby fibers with other fibers having different Brillouin frequency shift values (here, multi-mode fibers are used). Finally, we demonstrate a distributed measurement for detecting a 1.46-m-long strained section with a high signal-to-noise ratio.

preprint2014arXiv

Brillouin scattering in multi-core optical fibers

We measure the Brillouin gain spectra in two cores (the central core and one of the outer cores) of a ~3-m-long, silica-based, 7-core multi-core fiber (MCF) with incident light of 1550 nm wavelength, and investigate the Brillouin frequency shift (BFS) and its dependence on strain and temperature. The BFSs of both the cores are ~10.92 GHz, and the strain- and temperature-dependence coefficients of the BFS in the central core are 484.8 MHz/% and 1.08 MHz/K, respectively, whereas those in the outer core are 516.9 MHz/% and 1.03 MHz/K. All of these values are not largely different from those in a silica single-mode fiber, which is expected because the cores are composed of the same material (silica). The difference in the BFS dependence of strain between the two cores may originate from the difference in structural deformation when strain is applied to the fiber, which is a unique characteristic to MCFs. The future prospect on distributed strain and temperature sensing based on Brillouin scattering in MCFs is finally presented.

preprint2014arXiv

Fiber-Optic Interferometry Using Narrowband Light Source and Electrical Spectrum Analyzer: Influence on Brillouin Measurement

We observe an interference pattern using a simple fiber-optic interferometer consisting of an electrical spectrum analyzer and a narrowband light source, which is commonly employed for observing the Brillouin gain spectrum. This interference pattern expands well beyond the frequency range corresponding to the Brillouin frequency shift in silica fibers (approximately 11 GHz at 1550 nm). Using both silica single-mode and polymer optical sensing fibers, we then experimentally prove that the distinctive noise in a self-heterodyne-based Brillouin measurement with an unoptimized polarization state originates from the interference between the reference light and the Fresnel-reflected light. This noise can be almost completely suppressed by employing a delay line that is longer than the coherence length of the light source and by artificially applying a high loss near the open end of the sensing fiber.

preprint2014arXiv

Magnetic Field Amplification and Saturation in Turbulence Behind a Relativistic Shock

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

preprint2014arXiv

Simplified configuration of Brillouin optical correlation-domain reflectometry

We develop a simple and cost-efficient configuration of Brillouin optical correlation-domain reflectometry (BOCDR), the setup of which does not include an additional reference path used in standard BOCDR systems. The Fresnel-reflected light from an open end of a sensing fiber is used as reference light. The limitations of spatial resolution, measurement range, and their ratio are theoretically clarified, and then a distributed strain measurement with a < 100-mm spatial resolution and a 4.1-m measurement range is demonstrated with a high signal-to-noise ratio.

preprint2014arXiv

Spatial Growth of the Current-Driven Instability in Relativistic Jets

We have investigated the influence of velocity shear and a radial density profile on the spatial development of the current driven kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static non-propagating current driven kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of maximum toroidal magnetic field have a relatively stable configuration.

preprint2013arXiv

Polymer optical fiber fuse

Although high-transmission-capacity optical fibers are in demand, the problem of the fiber fuse phenomenon needs to be resolved to prevent the destruction of fibers. As polymer optical fibers become more prevalent, clarifying their fuse properties has become important. Here, we experimentally demonstrate a fuse propagation velocity of 21.9 mm/s, which is 1 to 2 orders of magnitude slower than that in standard silica fibers. The achieved threshold power density and proportionality constant between the propagation velocity and the power density are respectively 1/186 of and 16.8 times the values for silica fibers. An oscillatory continuous curve instead of periodic voids is formed after the passage of the fuse. An easy fuse termination method is presented herein, along with its potential plasma applications.

preprint2013arXiv

The Role of the Equation of State in Resistive Relativistic Magnetohydrodynamics

We have investigated the role of the equation of state in resistive relativistic magnetohydrodynamics using a newly developed resistive relativistic magnetohydrodynamic code. A number of numerical tests in one-dimension and multi-dimensions are carried out in order to check the robustness and accuracy of the new code. The code passes all the tests in situations involving both small and large uniform conductivities. Equations of state which closely approximate the single-component perfect relativistic gas are introduced. Results from selected numerical tests using different equations of state are compared. The main conclusion is that the choice of the equation of state as well as the value of the electric conductivity can result in considerable dynamical differences in simulations involving shocks, instabilities, and magnetic reconnection.

preprint2012arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. III. Rotating Relativistic Jets

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

preprint2012arXiv

Two Types of Ergospheric Jets from Accreting Black Holes: The Dichotomy of Fanaroff-Riley Galaxies

We investigate the extraction of the rotational energy of a black hole under different accreting environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can sustain a magnetic field with moderate strength at the event horizon, leading to the formation of relativistic jets with moderate luminosity and speed via the magnetohrodynamic Penrose process. When the accretion rate increases enough, on the other hand, the disk becomes geometrically thin near the horizon. In this slim disk, the denser plasmas can sustain a stronger magnetic field than that in a combined disk, leading to the formation of jets with greater luminosity and speed via the Blandford-Znajek processs. It is discussed that the former jets are associated with the Fanaroff-Riley (FR) I galaxies and the latter with FR II galaxies.

preprint2011arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability with A Sub-Alfvenic Jet: Temporal Properties

We have investigated the influence of a velocity shear surface on the linear and non-linear development of the CD kink instability of force-free helical magnetic equilibria in 3D. In this study we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than occur for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and non-linear development can be similar to the development of a static plasma column.

preprint2010arXiv

Magnetic Field Amplification by Turbulence in A Relativistic Shock Propagating through An Inhomogeneous Medium

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

preprint2010arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. II. Relaxation of Pulsar Wind Nebula

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter σ, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.