Featured Paper of the Week
26 news

20200420
Featured Paper of the WeekBlack Hole as a Quantum Field Configuration
What's going on inside black holes? The observations to date have taught us the nature of black holes outside and around it, but it does not tell us anything about the inside yet. This is because black holes are "almost" black: gravity is so strong that physical signals from the inside "almost" never come out. So, how can we search for the interior of black holes? It is necessary to return to the basic principles of physics and reconsider what a black hole is. In this paper, we succeeded in describing the inside of black holes based on general relativity and quantum mechanics (quantum field theory). What is important here is a fact that black holes evaporate: According to quantum mechanics, black holes emit weak light, slowly lose energy, and eventually evaporate (Hawking radiation). By taking in this effect from the beginning and analyzing the collapsing process of stars, we obtained the following result: The black hole is a compact object with a surface (instead of horizon) that looks like a conventional black hole from the outside and eventually evaporates without a singularity. (See figure.) Indeed, we explained why this picture can be achieved in quantum field theory. In addition, we investigated how much information is contained inside and showed that the amount (entropy) is given by the surface area. We also described how the energy decreases. This work provides the first fieldtheoretical formulation of the black hole, and I think it will be the basis for investigating how the information that enters the black hole will go out.

20200124
Featured Paper of the WeekDistance between collapsing matter and trapping horizon in evaporating black holes
How deep can an apple fall into the horizon of a black hole? Of course, in the case of a classic black hole, the apple reaches the center. However, the real world is described by quantum theory, and black holes evaporate by Hawking radiation. This changes the fate of the apple. Negative energy due to vacuum quantum fluctuations is generated near black holes, which enters the black holes and reduces their energy. If there is a horizon, the negative energy will make the horizon smaller. Then, the horizon becomes timelike (while it is null in the classical case). In this paper, we solved approximately the Einstein equation in the vicinity of the horizon and examined the physical (proper) distance between the falling apple and the horizon. It turned out that the distance is a few of the Planck length (due to the exponentially delayed time inside the horizon). This means that essentially no apples fall into the black hole. This result suggest that we should reconsider whether the black hole really has a horizon.

20200109
Featured Paper of the WeekAkizukiNakano vanishing theorem and its application for globally Fsplit 3folds
In algebraic geometry, it is important to study spaces (=varieties) which may admit singular points. "Smoothing" is an effective way to deal with singularities. Roughly speaking, a smoothing of a (singular) variety helps us to find another variety which is nonsingular and is "similar" to the original variety. Namikawa proved that a complex Fano 3fold whose singular points are at worst ordinary double points admits a smoothing. In this paper, we give an analogy of this result in the case where the Fano variety is not complex, but is defined over a field of positive characteristic. The key ingredient of the proof is the vanishing of the second cohomology of the tangent space, which is a vector space containing an obstruction to construct a smoothing. We prove a positive characteristic analogue of the AkizukiNakano type vanishing for 3dimensional varieties and as a corollary, we conclude the result about a smoothing.

20190722
Featured Paper of the WeekSynthetic dimensions and topological chiral currents in mesoscopic rings
In condensed matter systems, the topological effect of chiral currents is intimately related to twodimensional systems. Recently emerging concept of "synthetic dimensions" provides a way to explore twodimensional physics using one dimensional setups by using internal degrees of freedom as an additional dimension. Synthetic dimensions have so far been predominantly developed in atomic and optical physics. In this work, we extend this concept to condensed matter systems of a mesoscopic ring coupled to a nanomagnet. Regarding the spin degrees of freedom of a nanomagnet as a dimension, we interpret a current in the ring, which is locked to the nanomagnet's spin, as chiral edge currents of a twodimensional quantum Hall system. We thus provide a broader conceptual setting for synthetic dimensions, extending it to a rich domain of potential practical applications in condensed matter systems.

20190716
Featured Paper of the WeekThe KOvalued spectral flow for skewadjoint Fredholm operators
Spectral flow measures the number of eigenvalue crossings of continuous paths of Fredholm operators (which have finitedimensional kernels). Atiyah, Patodi and Singer first defined spectral flow as a way to study index theory for odd dimensional manifolds. Recent applications of index theory to topological phases, where antilinear symmetries may occur, have motivated the study of spectral flow on real Hilbert spaces. The paper develops a theory of spectral flow for skewadjoint Fredholm operators on real Hilbert spaces, which are a classifying space for KOtheory. Our construction applies to bounded and unbounded Fredholm operators, possibly with additional Clifford symmetries, and generalises all previous notions of analytic spectral flow. The results are also relevant to freefermionic topological phases, where the topological obstruction to two symmetric Hamiltonians having the same strong topological phase can be exactly measured by the KOvalued spectral flow.

20190709
Featured Paper of the WeekInstantons in Chiral Magnets
More than 50 years ago, Russian physicist Dzyaloshinsky has proposed the specific interaction in some kinds of magnets to explain weak ferromagnetic behavior of them, and soon after that, Professor Moriya has identified its microscopic origin as the spinorbit coupling. That interaction, now called the Dzyaloshinskii–Moriya (DM) interaction, is one of the most important topics in both experimental and theoretical study of magnetism. This is because it induces novel inhomogeneously ordered states in spin systems such as a chiral soliton lattice and Skyrmion crystal. These interesting states support nonvanishing "topological" charges and have been experimentally observed in some magnetic materials. In this paper, focussing on onedimensional antiferromagnetic spin chains with the DM interaction, we exhaustively construct instanton solutions with topological excitations and elucidate their properties based on the field theoretical (nonlinear sigma model) description of spin chains. As shown in Fig, the resulting spin configurations show helical ordering due to the DM interaction.

20190704
Featured Paper of the WeekRigidity of the mod 2 families SeibergWitten invariants and topology of families of spin 4manifolds
After the early 1980s, gauge theory has been used by mathematicians to study 4dimensional spaces (manifolds). The hottest mathematical subject studied by gauge theory is a sort of difference between the shape of topological (continuous) 4dimensional manifold and that of smooth 4dimensional manifold. On the other hand, in most of mathematical areas, it is quite significant to understand symmetry of a given mathematical object, such as a manifold. However, about the space of all symmetries of a 4dimensional manifold, only few things have been known. In this paper, we study the spaces of symmetries of 4dimensional manifolds. We revealed that gauge theory (more precisely SeibergWitten theory) can extract some difference between the shape of the space of symmetries of a topological 4dimensional manifold and that of a smooth 4dimensional manifold. In addition to gauge theory, we also used some classical but deep results in higherdimensional topology invented 50 years ago. Such a combination of gauge theory with classical results in topology is also a new and interesting direction of this paper.

20190613
Featured Paper of the WeekProbing localization and quantum geometry by spectroscopy
Localization has been one of the main topics of interest in condensed matter physics for more than half a century. Study of localizatoin has recently been entering into a new stage because of the developments in the concept of manybody localization, and also of the recent advances in the technology of quantumengineered systems such as ultracold atoms and trapped ions. In this paper, we propose a novel method to quantitatively study localization of a quantum state in the platform of quantumengineered systems. Our proposal makes use of the fluctuationdissipation theorem; we consider periodically modulating the quantum system and observing the exciation rate. We find that the localization of a quantum state is related to the integral of the excitation rate over the modulation frequency. Our proposal to probe localization does not require high resolutino microscopes to spatially resolve the quantum state. We apply our method to various examples which are of direct experimental relevance in ultracold atoms. Moreover, inspired by a relation between quantum fluctuations and the quantum metric, we describe how our scheme can be generalized to extract the full quantumgeometric tensor of manybody systems. The figure (left) shows how the wavefunction of two interacting particles in a harmonic trap spreads as one increases the interaction, and the figure (right) shows how such spreading can be probed through excitation rate measurements for various values of the interaction strengths.

20190528
Featured Paper of the WeekUncovering critical properties of the human respiratory syncytial virus by combining in vitro assays and in silico analyses
When a virion is left in a liquid medium over some time, it loses its ability to infect cells. This is often because the virion surface proteins become deformed, e.g. because of heat or the ionic content of the medium, such that they can no longer effectively bind cells and infect them. For many viruses (e.g. influenza, HIV and hepatitis C) this loss of infectivity over time follows an exponential. In this paper, we evaluated the rate at which the respiratory syncytial virus (RSV) loses infectivity over time. We were surprised when we found that the loss of infectivity of RSV in our data did not follow an exponential. We used a Weibull distribution to describe the decay and estimated its shape parameter to be k=0.30 [95% credible region: 0.270.36], which means the rate of decay gets smaller over time. This is typical of a bad production system wherein a lot of defective products are made which fail soon after production, leaving only wellmade products which become defective at a lower rate. In a 2013 paper, Liljeroos and others (doi: 10.1073/pnas.1309070110) showed that as RSV buds from the surface of the cell, bud pinching (closing the bud to free the virion from the cell) is physically difficult, and often leads to misformed buds. The bad structure leads the filamentous (wormshaped) RSV bud to “plop” outwards (become spherical), flip its surface F proteins (responsible for RSV attachment to host cells) from their pretrigger (very infectious) to their posttrigger (poorly infectious) state. This is consistent with our finding of k<1, wherein RSV that don’t pinch their bud properly upon cell exit fail (lose infectivity) very early on, and the remaining RSV population which pinched properly will decay following a much less steep exponential.

20190515
Featured Paper of the WeekComments on holographic entanglements in cutoff AdS
AdS/CFT correspondence tells us that a gravity in AdS space is equivalent to CFT which lives on the boundary of the AdS. Recently it has been proposed that an integrable TT deformation of CFTs produces a deformed AdS/CFT. According to the conjecture, the AdS boundary is located at a finite radial cutoff. To investigate the deformed correspondence, we studied holographic entanglement entropy in the cutoff AdS, which is well known as a probe of information theoretic quantity. In particular, we explored its phase transitions. For twointerval entanglement entropy, the transition point monotonically decreases with a deformation parameter, which means that by the TT deformation the degrees of freedom in subsystems are decreasing. Our result implies that the effect of the TT deformation can be regarded as the rescaling of the energy scale.

20190312
Featured Paper of the WeekDiffusive NambuGoldstone modes in quantum timecrystals
Symmetry and its breaking is one of the most important notions in modern physics. It is known that when a continuum symmetry is spontaneously broken, there appears a gapless excitation called the Nambu–Goldstone (NG) mode. Such a mode appears in several length scales. In the femtometer scale, pions that are responsible for nuclear force, are the NG modes associated with chiral symmetry breaking. In the meter scale, the sound mode in solid is the NG mode associated with spontaneous breaking of translational symmetry. In this paper, we study the NG modes associated with spontaneous breaking of the continuous timetranslation symmetry in an open quantum system. For this purpose, we study a system with a van der Pol type nonlinearfriction. By considering small fluctuations around a timeperiodic meanfield solution, we find that a gapless collective mode, which is a diffusive one, necessarily appears; this is nothing but the NG mode of a time crystal. We also show that noncommutative breaking of the timetranslation and U(1) symmetries results in mixing of the NG modes, and the (typically) propagating NG mode appears.

20190306
Featured Paper of the WeekNonsinusoidal Waveform in TemperatureCompensated Circadian Oscillations
We have autonomous daily rhythm in our body, socalled “circadian clock” for which we can wake up without alarm clock. The rhythm governs the timing of physiological events such as hormone secretion. Curiously, the period of our daily rhythms is stable to temperature, a phenomenon known as “temperature compensation.” Temperature compensation has been a mystery for many years because underlying reaction processes tend to accelerate with temperature. If our internal clock was sensitive to temperature, the clock would not function properly any more. To uncover the mystery, we developed a simpler model of circadian clock based on experimental evidences whilst many realistic simulators have been proposed. Using the model, we derived the period formula. The formula states that “the more nonsinusoidal, the longer period.” From the formula, we predicted that more nonsinusoidal waveform at higher temperature should be observed in reality which is necessary to cancel out the period shortening effect at higher temperature. Unexpectedly, we obtained a similar period formula for a limit cycle oscillation of electrical circuit (van der Pol model), indicating that nonsinusoidal waveform lengthens period. How can circadian clock and electrical circuit behave similarly? We wish to understand the reason in a future.

20190207
Featured Paper of the WeekQuantum annealing for polynomial systems of equations
The advent of Noisy IntermediateScale Quantum (NISQ) quantum computers has galvanized efforts towards discovering nearterm applications. An algorithm for solving polynomial systems of equations was proposed by Dr. Chang (iTHEMS) in collaboration with Dr. Gambhir (Lawrence Livermore National Lab), Dr. Humble (Oak Ridge National Lab), and Dr. Sota (RCCS) and implemented a linear solver on a DWave quantum annealer. While the problems are currently limited to sizes that are easily solved by classical computers, the team showed that the quantum algorithm exhibits constant scaling with increasing condition number, in direct contrast with classical methods. Additionally, the quantum algorithm may also be applied iteratively to exponentially decrease the relative residual, allowing for the classical solution to be reproduced by the quantum computer to single precision. However, the scaling with problem size is unfortunately exponentially bad, reflecting limitations of current quantum computers. Fortunately, there is a great amount of interest and effort put fourth by the greater quantum annealing community geared towards tackling this problem, including using inhomogeneous driving fields, reverse annealing, and even hardware developments towards universal quantum annealers.

20190131
Featured Paper of the WeekBack Reaction of 4D Conformal Fields on Static Geometry
In classical mechanics, a black hole is described by a vacuum solution with the horizon of the Einstein equation. For spherical case, it is the Schwarzschild metric, and the location of the horizon is given by the Schwarzschild radius. (Note that the Schwarzschild radius can also be defined even for a star without horizon.) In quantum mechanics, a black hole evaporates and information inside it seems to be lost, which is contradict to the principle of quantum mechanics. An effective way to address this problem is to consider again “What is the black hole in quantum mechanics?” In this paper, we examined how robust the use of the Schwarzschild metric to represent a black hole is in quantum mechanics. We consider conformal matters (e.g. electromagnetic field) and introduce the quantum effect (4D conformal anomaly) into the Einstein equation, which necessarily makes the equation nonvacuum. We start from the Schwarzschild metric, add the quantum effect perturbatively, and solve the Einstein equation in a selfconsistent manner. Then, we showed that the quantum effect can play a crucial role in shaping the static geometry near the Schwarzschild radius. The geometry depends on a parameter corresponding to a boundary condition, and the existence of the horizon requires the finetuning. Therefore, in quantum mechanics, a typical static spherical solution does not have a horizon.

20190122
Featured Paper of the WeekIndex theory and topological phases of aperiodic lattices
The paper studies the interactions between dynamical systems, represented via groupoid operator algebras, and index theory, described using Kasparov's bivariant Kfunctor. The key motivating example is the dynamical systems associated to Delone sets. Delone sets give a model of aperiodic atomic configurations, where the material may be amorphous or quasicrystalline. The atomic configuration does not have a group structure, but key properties can be encoded in a discrete groupoid. We construct bivariant Kcycles, a generalisation of Ktheory classes, that encode the dynamics of the aperiodic lattice. For the case of quasicrystalline configurations, we obtain finer information coming from the longrange aperiodic order of the material. We then use our constructed Kcycles to define bulk topological phases of aperiodic materials. A bulkboundary correspondence is also shown, which arises from factorisation properties of the constructed Kcycles.

20181130
Featured Paper of the WeekEngaging the Public with Supernova and Supernova Remnant Research Using Virtual Reality
On April 21, 2018, residents of Wako came to RIKEN for the Open Day. Over one hundred people visited the exhibit of the Astrophysical Big Bang Laboratory (ABBL) and interacted with real scientific data using virtual reality (VR). With the help of the entire ABBL and several students from Rikkyo U, people as young as 2 and as old as 80 got to explore a supernova remnant in 3D, as if it was in the room with them. (A supernova remnant results from the explosion of a star; ABBL does numerical simulations to understand the explosion mechanism.) In this report, Gilles Ferrand (ABBL/iTHEMS) and Don Warren (iTHEMS) explain what makes VR exciting to explore your data and communicate your research, and how to make your own VR booth at next year's Open Day — or any other event! We are happy to discuss this exhibit, and answer any questions you have, at a coffee meeting, wine meeting, or another occasion.

20181112
Featured Paper of the WeekGamowTeller transitions from highspin isomers in $N=Z$ nuclei
Atomic nuclei are composed of protons and neutrons. We may imagine like a certain number of (quantum) balls putting in a (selfbound) container. On the one hand, the nuclear magic numbers show the singleparticle features of nuclei; On the other hand, the giant resonances show the emergent collective features. For example, an intuitive image of giant dipole resonance is: All protons form a group while all neutrons form another group, and these two groups oscillate against to each other collectively. GamowTeller resonance is an important kind of giant resonances involving not only the spin but also the isospin degrees of freedom in nuclei. So far, most, if not all, of experimental and theoretical studies on this direction focus on the resonances from the nuclei being in their ground states, largely due to the experimental limitations. However, it is expected that, in the near future, our worldleading nuclear facility in RIKEN can generate highquality radioactive beams with nuclei being in their highspin isomeric states. Motived by this perspective, we predict that the GamowTeller transitions from the highspin isomers could be much more collective than the corresponding transitions from their ground states, by taking nuclei 52Fe and 94Ag as examples. This would be one of the promising collaborations between theorists and experimentalists in RIKEN.

20181102
Featured Paper of the WeekEffective field theory of timetranslational symmetry breaking in nonequilibrium open system
Although we know "all things are in a state of flux. (Patna rhei.)", it has been still an interesting question whether we have "universal" features in such timeevolution of all things. Considering the nonequilibrium open systems like the Brownian particles, we have tackled this problem from the point of view of timetranslational symmetry. We regard the transient diffusive (or timeoscillating) behaviors as "symmetrybroken phases of matter", and develop a way to describe their spacetime evolution. Our formalism is based on the socalled effective field theory (EFT) originally developed in the context of particle physics and now used in almost every branches of physics. The advantage of the EFT is that it enables us to capture the possible universal properties in the largescale (lowenergy) behaviors of systems. As a result, we find the existence of the NambuGoldstonelike mode during transient timeevolution, which gives the universal diffusive behaviors of all things in a state of flux.

20181018
Featured Paper of the WeekNonequilibrium Chiral Magnetic/Vortical Effects in Viscous Fluids
Chiral magnetic and vortical effects are novel effects in quantum fluids with chiral anomaly. We utilize the chiral kinetic theory in a relaxationtime approximation to investigate the nonlinear anomalous responses of chiral fluids with viscous effects. Unlike the cases in equilibrium, it is found that the chiral magnetic effect and chiral vortical effect are modified by the shear and bulk strengths. Particularly, the shear strength could result in charged Hall currents for chiral magnetic and chiral vortical effects, which propagate perpendicular to applied magnetic fields and vorticity. These quantum corrections stemming from side jumps and anomalies are dissipative and pertinent to interactions.

20181009
Featured Paper of the WeekOnline Model Parameter Estimation With Ensemble Data Assimilation in the Real Global Atmosphere
When and where will it rain, and how heavy? – this is a central question that meteorology tries to address. Numerical weather prediction (NWP) is a major approach using data assimilation with mathematical models to predict the weather. The NWP models discretize the earth atmosphere into 3dimensional grids, and compute the evolution of the atmospheric states based on physical processes (e.g., fluid dynamics, radiation, and water phase changes). The NWP models cannot include subgridscale and complex processes explicitly, and these complex phenomena are represented by simplified equations socalled “parameterization”. Parameterization usually contains tunable model parameters, and manual tuning of these parameters is a tedious but important task. This study explores an objective and autonomous approach to optimizing these parameters using data assimilation. We chose a parameter of the subgridscale parameterization of raindrop initiation processes of a global NWP model. We successfully mitigated the overproduced precipitation of the NWP model by estimating the model parameter with satelliteobserved precipitation data. Figure: Global precipitation forecasts (mm 6h1) at 0000 UTC on 16 June 2014 by an NWP model. (left) control experiment with the default model setting, (middle) test experiment with the model parameter estimation, and (right) satellite observation, respectively. Overproduced precipitation over ocean in the control experiment is successfully mitigated by the model parameter estimation.

20181001
Featured Paper of the WeekSelflearning Monte Carlo method with BehlerParrinello neural networks
Self Learning MonteCarlo (SLMC) method is one of the recent promising applications of machine learning techniques to computational physics especially MarcovChain MonteCarlo (MCMC) simulations in statistical physics systems. In SLMC, we prepare an effective Hamiltonian with tunable coupling constants and try to reproduce the value of the real Hamiltonian by adjusting the coupling constants. This procedure corresponds to the supervised machine learning. After the learning, we employ it for global updates in MCMC, and it drastically reduces the autocorrelation, i.e. similarity of configurations in the MarcovChain, thanks to the local structure of the effective Hamiltonian. In our paper, we present two novel techniques for SLMC in a quantum MonteCarlo simulation. First technique is the use of neural networks. Based on the neural net architecture by Behler and Parrinello in the context of molecular dynamics machine learning, we develop how to construct general purposed effective Hamiltonian for SLMC. Second contribution is a new proposal for regularization in SLMC that we call batchatom normalization. It is a generalization of the wellknown technique in deep learning, and we observe it drastically improve the learning procedure for effective Hamiltonian represented by neural network.

20180919
Featured Paper of the WeekThe relative MishchenkoFomenko higher index and almost flat bundles
Existence of a positive scalar curvature (PSC) metric has been one of the central topics in differential topology of higher dimensional manifolds. Index theory provides a topological invariant whose vanishing is an effective necessary condition of the existence of a PSC metric. Recently, ChangWeinbergerYu introduce a new index theoretic invariant called the relative higher index detecting the nonexistence of a PSC metric on spin manifolds with boundary. The aim of this paper is to understand this invariant from Riemannian gemetric point of view. Firstly, we introduce a new definition of the relative higher index and prove that it is equivalent to the existing definitions. Secondly, by using our new definition, we relate the relative higher index with a Riemannian geometry of vector bundles; the index pairing with almost flat vector bundles. This is a relative analogue of the theory of GromovLawson for closed spin manifolds.

20180828
Featured Paper of the WeekTwo Ck1δ transcripts regulated by m6A methylation code for two antagonistic kinases in the control of the circadian clock
Thanks to the accumulation of molecular knowledges about biological clock, we can predict a missing link “X” in the system by using a realistic simulator. Mealtime is known to influence our sleepwake timing. But the mechanism has been unknown. This time, the group of Dr. Okamura at Kyoto University discovered a new enzyme (named, Ck1δ2) that possibly links between our metabolic system and biological clock. Based on the experimental results by Dr. Okamura, two iTHEMS members (Gibo and Kurosawa) predicted the specific biochemical process, activated by the new enzyme. The prediction was conducted by using a 190 variables model with all genes and proteins. Fortunately, the prediction was confirmed experimentally. The point for the prediction was the activation of the new enzyme “slowed” the biological clock in the experiment while the activation of most enzymes “accelerate” the clock in the simulator. Now, the mechanism is studied by using a simpler model. Then, a possible scenario is that aminoacids in our food influence our sleepwake timing via the new enzyme.

20180730
Featured Paper of the WeekStructural Bifurcation Analysis in Chemical Reaction Networks
In living cells, large numbers of chemical reactions form complicated networks (e.g. metabolic networks). It is considered that dynamics arising from such networks are the origin of biological functions. In this study, we developed a theoretical method to analyze bifurcation properties, i.e. properties of structural changes of steadystate solutions, from topology of reaction networks. Such a topological approach is advantageous because reaction kinetics is not well known in detail. Specifically, we showed that a reaction network can be decomposed into subnetworks based on a certain topological criteria, and showed that any bifurcation point is associated with a particular subnetwork ((3) in the figure). Further, we identified the parameters inducing bifurcations ((4) in the figure) and chemical concentrations exhibiting bifurcating behaviors ((5) in the figure). Our method is particularly useful to study large networks, since it allows us network decomposition. Biologically, bifurcation properties are related with plasticity of living systems, i.e. qualitative changes of physiological states induced by environmental conditions or geneexpressions. Our result suggests that biological plasticity arises from network topology.

20180723
Featured Paper of the WeekHilbert schemes of two points on K3 surfaces and certain rational cubic fourfolds
The goal of birational geometry is the classification of algebraic varieties up to birational equivalence. An algebraic variety is called rational if it is birationally equivalent to the projective space. In this paper, Genki Ouchi studied the rationality problem of (complex) cubic fourfolds, that is four dimensional complex hypersurface defined by a polynomial of degree 3. Conjecturelly, very general cubic fourfolds are irrational. However, no cubic fourfold has been proven to be irrational so far. On the other hand, there are five known examples of rational cubic fourfolds. They expect that the mysterious relation between rational cubic fourfolds and K3 surfaces is a key to solve the rationality problem. There are two inconsistent conjectures about it. So they should modify one of them at least. Together with previous works, he proved that known rational cubic fourfolds satisfy both conjectures. To modify the conjectures, we have to find a new rational cubic fourfold.

20180717
Featured Paper of the WeekProduction of intense episodic Alfvén pulses: GRMHD simulation of black hole accretion disks
The episodic dynamics of the magnetic eruption of spinning black hole (BH) accretion discs and the associated intense shapeup of their jets are studied via threedimensional generalrelativistic magnetohydrodynamics (GRMHD). The embedded magnetic fields in the disc are amplified by a magnetorotational instability (MRI) so large as to cause an eruption of the magnetic field (reconnection) and large chunks of matter accrete episodically toward the roots of the jets upon such an event. We also find that the eruption events produce intensive Alfvén pulses, which propagate through the jets. After the eruption, the disc returns to the weakly magnetic state. Such disc activities cause shorttime variabilities in mass accretion rate at the event horizon, as well as electromagnetic luminosity inside the jet. Since the dimensionless strength parameter a0 = eE/m_eωc of these Alfvén wave pulses is extremely high for a substantial fraction of Eddington accretion rate accretion flows on to supermassive black holes, the Alfvén shocks turn into ultrarelativistic (a0 >> 1) bow wake acceleration, manifesting as ultrahighenergy cosmic rays and electrons, which finally emit gammarays. Since our GRMHD model has universality in its spatial and temporal scales, it is applicable to a wide range of astrophysical objects, ranging from active galactic nuclei (AGNs, the primary target of this research) to microquasars. Properties such as the time variabilities of blazar gammaray flares and the spectrum observed by the Fermi Gammaray Observatory are explained well by linear acceleration of electrons by a bow wake. Figure caption: A snap shot of accreting gas onto a spinning black hole and Poynting flux dominated jet formation by 3dimensional general relativistic magnetohydrodynamic simulation. Logscaled inverse of plasma beta (magnetic pressure / thermal pressure) (xy plane), mass density (yz plane), and magnetic pressure (xz plane) are shown.
26 news