-
Circumnuclear Multi-phase Gas in the Circinus Galaxy. VI. Detectability of Molecular Inflow and Atomic Outflow
Authors:
Shunsuke Baba,
Keiichi Wada,
Takuma Izumi,
Yuki Kudoh,
Kosei Matsumoto
Abstract:
Recent submillimeter observations have revealed signs of pc-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multi-wavelength observations in t…
▽ More
Recent submillimeter observations have revealed signs of pc-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multi-wavelength observations in this paper series. The major axis position-velocity diagram (PVD) of CO(3-2) reproduces the observed faint emission at the systemic velocity, and our calculations confirm that this component originates from failed winds falling back to the disk plane. The minor-axis PVD of [CI]($^3P_1$-$^3P_0$), when created using only the gas with positive radial velocities, presents a sign of blue- and redshifted offset peaks similar to those in the observation, suggesting that the observed peaks indeed originate from the outflow, but that the model may lack outflows as strong as those in the Circinus galaxy. Similar to the observed HCN(3-2), the similar dense gas tracer HCO$^+$(3-2) can exhibit nuclear spectra with inverse P-Cygni profiles with $\sim$0.5 pc beams, but the line shape is azimuthally dependent. The corresponding continuum absorbers are inflowing clumps at 5-10 pc from the center. To detect significant absorption with a high probability, the inclination must be fairly edge-on ($\gtrsim$85$^\circ$), and the beam size must be small ($\lesssim$1 pc). These results suggest that HCN or HCO$^+$ and [CI] lines are effective for observing pc-scale inflows and outflows, respectively.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Formation of unipolar outflow and $\textit{protostellar rocket effect}$ in magnetized turbulent molecular cloud cores
Authors:
Daisuke Takaishi,
Yusuke Tsukamoto,
Miyu Kido,
Shigehisa Takakuwa,
Yoshiaki Misugi,
Yuki Kudoh,
Yasushi Suto
Abstract:
Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular m…
▽ More
Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular momentum transport to the protostars and protoplanetary disks. Here, we perform the three-dimensional non-ideal magnetohydrodynamics simulations to investigate the formation and early evolution of the asymmetrical protostellar outflows in magnetized turbulent isolated molecular cloud cores. We find, for the first time to our knowledge, that the unipolar outflow forms even in the single low-mass protostellar system. The results show that the unipolar outflow is driven in the weakly magnetized cloud cores with the dimensionless mass-to-flux ratios of $μ=8$ and $16$. Furthermore, we find the $\textit{protostellar rocket effect}$ of the unipolar outflow, which is similar to the launch and propulsion of a rocket. The unipolar outflow ejects the protostellar system from the central dense region to the outer region of the parent cloud core, and the ram pressure caused by its ejection suppresses the driving of additional new outflows. In contrast, the bending bipolar outflow is driven in the moderately magnetized cloud core with $μ=4$. The ratio of the magnetic to turbulent energies of a parent cloud core may play a key role in the formation of asymmetrical protostellar outflows.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
Circumnuclear Multi-phase Gas in the Circinus Galaxy. V. The Origin of the X-Ray Polarization in the Circinus Galaxy
Authors:
Atsushi Tanimoto,
Keiichi Wada,
Yuki Kudoh,
Hirokazu Odaka,
Ryosuke Uematsu,
Shoji Ogawa
Abstract:
The Imaging X-ray Polarimetry Explorer (IXPE) detected X-ray polarization in the nearest Seyfert 2 galaxy, the Circinus galaxy, for the first time. To reproduce the IXPE results, we computed the degree of polarization based on two types of radiative hydrodynamic simulations: a parsec-scale three-dimensional model and a sub-parsec-scale axisymmetric model with a higher spatial resolution. In a seri…
▽ More
The Imaging X-ray Polarimetry Explorer (IXPE) detected X-ray polarization in the nearest Seyfert 2 galaxy, the Circinus galaxy, for the first time. To reproduce the IXPE results, we computed the degree of polarization based on two types of radiative hydrodynamic simulations: a parsec-scale three-dimensional model and a sub-parsec-scale axisymmetric model with a higher spatial resolution. In a series of papers, we confirmed that these models naturally explain the multi-wavelength observations of the Circinus galaxy from radio to X-rays. We used a Monte Carlo Simulation for Astrophysics and Cosmology code to compute the linear polarization of continuum emission. We found that the degree of polarization based on the parsec-scale radiation-driven fountain model was smaller than that observed with the IXPE. The degree of polarization based on the sub-parsec-scale model depends on the hydrogen number density of the disk ($d$), and the degree of polarization obtained from our simulation is consistent with that observed with the IXPE in the case of $\log d/\mathrm{cm}^{-3} \geq 13$. We investigate where the photons are Compton scattered and imply that the origin of the X-ray polarization in the Circinus galaxy is the outflow inside $0.01 \ \mathrm{pc}$. In this case, the degree of polarization may change over a timescale of approximately ten years.
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
Supermassive black hole feeding and feedback observed on sub-parsec scales
Authors:
Takuma Izumi,
Keiichi Wada,
Masatoshi Imanishi,
Kouichiro Nakanishi,
Kotaro Kohno,
Yuki Kudoh,
Taiki Kawamuro,
Shunsuke Baba,
Naoki Matsumoto,
Yutaka Fujita,
Konrad R. W. Tristram
Abstract:
Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecul…
▽ More
Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecular gas is evident on sub-parsec scales. We calculate that less than 3% of this inflow is accreted by the black hole, with the rest being ejected by multiphase outflows, providing feedback to the host galaxy. The observations also reveal a dense gas disk surrounding the inflow; the disk is gravitationally unstable which drives the accretion into the central ~1 parsec.
△ Less
Submitted 13 November, 2023; v1 submitted 6 May, 2023;
originally announced May 2023.
-
Multiphase Gas Nature in the Sub-parsec Region of the Active Galactic Nuclei I: Dynamical Structures of Dusty and Dust-free Outflow
Authors:
Yuki Kudoh,
Keiichi Wada,
Nozomu Kawakatu,
Mariko Nomura
Abstract:
We investigated dusty and dust-free gas dynamics for a radiation-driven sub-pc scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole $10^7 M_\odot$ and bolometric luminosity $10^{44}$ erg s$^{-1}$ based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven ``lotus-like'' multi-shell outflow is launched from the inner part (…
▽ More
We investigated dusty and dust-free gas dynamics for a radiation-driven sub-pc scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole $10^7 M_\odot$ and bolometric luminosity $10^{44}$ erg s$^{-1}$ based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven ``lotus-like'' multi-shell outflow is launched from the inner part ($r \lesssim 0.04$ pc) of the geometrically thin disk, and it repeatedly and steadily produces shocks as mass accretion continues through the disk to the center. The shape of the dust sublimation radius is not spherical and depends on the angle ($θ$) from the disk plane, reflecting the non-spherical radiation field and nonuniform dust-free gas. Moreover, we found that the sublimation radius of $θ\sim 20$-$60$ deg varies on a timescale of several years. The ``inflow-induced outflow" contributes the obscuration of the nucleus in the sub-parsec region. The column density of the dust-free gas is $N_{\rm H} \gtrsim 10^{22}$ cm$^{-2}$ for $r \lesssim 0.04$ pc. Gases near the disk plane ($θ\lesssim 30$ degree) can be the origin of the Compton-thick component, which was suggested by the recent X-ray observations of AGNs. The dusty outflow from the sub-parsec region can be also a source of material for the radiation-driven fountain for a larger scale.
△ Less
Submitted 11 May, 2023; v1 submitted 12 April, 2023;
originally announced April 2023.
-
Multi-phase gas nature in the sub-pc region of the active galactic nuclei II: Possible origins of the changing-state AGNs
Authors:
Keiichi Wada,
Yuki Kudoh,
Tohru Nagao
Abstract:
Multi-wavelength observations of active galactic nuclei (AGNs) often reveal various time scales of variability. Among these phenomena, "changing-look AGNs" are extreme cases where broad emission lines become faint/bright or even disappear/emerge between multi-epoch observations, providing crucial information about AGN internal structures. We here focus on "changing-state" AGNs, specifically invest…
▽ More
Multi-wavelength observations of active galactic nuclei (AGNs) often reveal various time scales of variability. Among these phenomena, "changing-look AGNs" are extreme cases where broad emission lines become faint/bright or even disappear/emerge between multi-epoch observations, providing crucial information about AGN internal structures. We here focus on "changing-state" AGNs, specifically investigating the transition of optical spectra over years to tens of years. Based on the axisymmetric radiation-hydrodynamical simulations (Paper I) for the gas dynamics within the dust-sublimation radius, we investigate the spectral properties of ionized gas exposed to the radiation from an AGN with a 10^7 Msun supermassive black hole. We find significant time-dependent variations in the Balmer emission lines by utilizing post-process pseudo-three-dimensional calculations and the spectral synthesis code CLOUDY. The equivalent width of Halpha and Hbeta changes by a factor of 3, or the emission lines even disappear during 30 years for the same viewing angle. The time-dependent behaviour arises primarily from gas dynamics, particularly the formation of non-steady, radiation-driven outflows within the innermost region of the disc (r <10^-3 pc). The intricate interplay between non-spherical radiation sources at the core of AGNs and the dynamic behavior of gas within the dust sublimation radius gives rise to radiation-driven outflows. This non-steady outflow potentially contributes to the observed variability in Balmer line emissions over multi-year timescales in certain AGNs.
△ Less
Submitted 22 September, 2023; v1 submitted 28 March, 2023;
originally announced March 2023.
-
Probing dynamics and thermal properties inside molecular tori with CO rovibrational absorption lines
Authors:
Kosei Matsumoto,
Takao Nakagawa,
Keiichi Wada,
Shunsuke Baba,
Shusuke Onishi,
Taisei Uzuo,
Naoki Isobe,
Yuki Kudoh
Abstract:
A recent hydrodynamic model, "radiation-driven fountain model" (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers confirmed that this picture is consistent with multi-wavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovi…
▽ More
A recent hydrodynamic model, "radiation-driven fountain model" (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers confirmed that this picture is consistent with multi-wavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovibrational absorption lines ($ΔJ=\pm1$, $v=0-1$, $λ\sim 4.7$ $\mathrm{μm}$) could probe physical properties of the inside tori. However, the origin of the CO absorption lines has been under debate. In this paper, we investigate the origin of the absorption lines and conditions for detecting them by performing line radiative transfer calculations based on the radiation-driven fountain model. We find that CO rovibrational absorption lines are detected at inclination angles $θ_\mathrm{obs} = 50-80$ $^{\circ}$. At the inclination angle $θ_\mathrm{obs} = 77$ $^{\circ}$, we observe multi-velocity components: inflow ($v_\mathrm{LOS}=30$ $\mathrm{kms^{-1}}$), systemic ($v_\mathrm{LOS}=0 \, \mathrm{kms^{-1}}$), and outflows ($v_\mathrm{LOS}=-75,\, -95,$ and $-105$ $\mathrm{kms^{-1}}$). The inflow and outflow components ($v_\mathrm{LOS}= 30$ and $-95$ $\mathrm{kms^{-1}}$) are collisionally excited at the excitation temperature of $186$ and $380$ K up to $J=12$ and $4$, respectively. The inflow and outflow components originate from the accreting gas on the equatorial plane at $1.5$ pc from the AGN center and the outflowing gas driven by AGN radiation pressure at $1.0$ pc, respectively. These results suggest that CO rovibrational absorption lines can provide us with the velocities and kinetic temperatures of the inflow and outflow in the inner a-few-pc regions of AGN tori, and the observations can probe the gas circulation inside the tori.
△ Less
Submitted 1 June, 2022;
originally announced June 2022.
-
Circumnuclear Multi-phase Gas in Circinus Galaxy IV: absorption owing to high-$J$ CO rotational transitions
Authors:
Taisei Uzuo,
Keiichi Wada,
Takuma Izumi,
Shunsuke Baba,
Kosei Matsumoto,
Yuki Kudoh
Abstract:
We studied the absorption features of CO lines against the continuum originating from the heated dust in the obscuring tori around active galactic nuclei (AGNs). We investigated the formation of absorption lines corresponding to the CO rotational transitions using three-dimensional non-LTE line transfer simulations considering the dust thermal emission. As in Papers I--III of this series, we perfo…
▽ More
We studied the absorption features of CO lines against the continuum originating from the heated dust in the obscuring tori around active galactic nuclei (AGNs). We investigated the formation of absorption lines corresponding to the CO rotational transitions using three-dimensional non-LTE line transfer simulations considering the dust thermal emission. As in Papers I--III of this series, we performed post-processed radiative transfer calculations using the "radiation-driven fountain model" (wada2016}, which yields a geometrically thick obscuring structure around the nucleus. This model is consistent with the spectral energy distribution of the nearest type-2 Seyfert galaxy, the Circinus galaxy. We found that the continuum-subtracted channel maps of $J = 4-3$ and higher transitions show absorption regions along the disk mid-plane for an edge-on viewing angle. The spectra consist of multiple absorption and emission features, reflecting the internal inhomogeneous and turbulent structure of the torus. The deepest absorption feature is caused by the gas on the near-side of the torus between $r =10$ and 15 pc, which is located in front of the AGN-heated dust inside $r \simeq 5$ pc. We also found that a spatial resolution of 0.5--1.0 pc is necessary to resolve the absorption features. Moreover, the inclination angle must be close to the edge-on angle (i.e., $\sim 85^\circ$) to observe the absorption features. The findings of the present study imply that combining our radiation-hydrodynamic model with high-resolution observations of CO (7-6) by ALMA can provide new information about the internal structure of the molecular tori in nearby AGNs.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
Large-scale magnetic field structure of NGC 3627 based on magnetic vector map
Authors:
Kohei Kurahara,
Hiroyuki Nakanishi,
Yuki Kudoh
Abstract:
We analyzed the data of Stokes $I$, $Q$, and $U$ in C- and X-bands and investigated the large-scale magnetic field structure of NGC 3627. The polarization intensity and angle in each band were derived using Stokes $Q$ and $U$ maps. The rotation measure was calculated using the polarization-angle maps. Moreover, the magnetic field strength was calculated by assuming energy equipartition with the co…
▽ More
We analyzed the data of Stokes $I$, $Q$, and $U$ in C- and X-bands and investigated the large-scale magnetic field structure of NGC 3627. The polarization intensity and angle in each band were derived using Stokes $Q$ and $U$ maps. The rotation measure was calculated using the polarization-angle maps. Moreover, the magnetic field strength was calculated by assuming energy equipartition with the cosmic ray electrons. The structure of the magnetic field was well aligned with the spiral arms, which were consistent with those in the former studies. We applied the magnetic vector reconstruction method to NGC 3627 to derive a magnetic vector map, which showed that northern and southern disks were dominant with inward and outward magnetic vectors, respectively. Furthermore, we discussed the large-scale structure of the magnetic field in NGC 3627 and observed that the structure is bi-symmetric spiral in nature, and that the number of magnetic field mode is $ m_{\rm B} = 1 $ in outer region of galaxy. In addition, NGC 3627 has a mode of two spiral arms that were clearly visible in an optical image. The ratio of the mode of spiral arms to that of magnetic field is 2:1. In terms of NGC 3627, the large-scale magnetic field may be generated via the parametric resonance induced by the gravitational potential of the spiral arms.
△ Less
Submitted 17 December, 2020; v1 submitted 9 December, 2020;
originally announced December 2020.
-
Multiphase Circumnuclear Gas in a Low-$β$ Disk: Turbulence and Magnetic Field Reversals
Authors:
Yuki Kudoh,
Keiichi Wada,
Colin Norman
Abstract:
We studied the magnetic field structures and dynamics of magnetized multiphase gas on parsec scales around supermassive black holes by using global 3D magnetohydrodynamics (MHD) simulations. We considered the effect of radiative cooling and X-ray heating due to active galactic nuclei (AGNs). The gas disk consists of a multiphase gas with (1) cold ($\leq 10^3$ K) and thin, and (2) warm (…
▽ More
We studied the magnetic field structures and dynamics of magnetized multiphase gas on parsec scales around supermassive black holes by using global 3D magnetohydrodynamics (MHD) simulations. We considered the effect of radiative cooling and X-ray heating due to active galactic nuclei (AGNs). The gas disk consists of a multiphase gas with (1) cold ($\leq 10^3$ K) and thin, and (2) warm ($\sim 10^4$ K) and thick components with a wide range of number densities. The turbulent magnetic energy at maximum is comparable to the thermal and turbulent kinetic energies in the turbulent motion. We confirmed that the turbulent velocity of the warm gas in the ambient cold gas is caused by magnetoconvective instability. The turbulent magnetic field due to magnetorotational instability (MRI) is developed in the disk, but the mean toroidal magnetic field dominates and supports in a quasi-steady state, where the plasma-$β$, the ratio between gas pressure and magnetic pressure, is low ($β< 1$). As often seen in adiabatic MHD simulations of rotating disks, the direction of the mean toroidal field periodically reverses with time even in multiphase gas structures. The direction reversal is caused by magnetic flux vertically escaping from the disk and by the combination of the MRI and the Parker instability.
△ Less
Submitted 8 November, 2020; v1 submitted 16 August, 2020;
originally announced August 2020.
-
Two-temperature Magnetohydrodynamic simulations for sub-relativistic AGN jets:Dependence on the fraction of the electron heating
Authors:
T. Ohmura,
M. Machida,
K. E. Nakamura,
Y. Kudoh,
R. Matsumoto
Abstract:
We present the results of two-temperature magnetohydrodynamic simulations of the propagation of sub-relativistic jets of active galactic nuclei. The dependence of the electron and ion temperature distributions on the fraction of electron heating fe at the shock front is studied for fe=0, 0.05, and 0.2. Numerical results indicate that in sub-relativistic, rarefied jets, the jet plasma crossing the…
▽ More
We present the results of two-temperature magnetohydrodynamic simulations of the propagation of sub-relativistic jets of active galactic nuclei. The dependence of the electron and ion temperature distributions on the fraction of electron heating fe at the shock front is studied for fe=0, 0.05, and 0.2. Numerical results indicate that in sub-relativistic, rarefied jets, the jet plasma crossing the terminal shock forms a hot, two-temperature plasma in which the ion temperature is higher than the electron temperature. The two-temperature plasma expands and forms a backflow referred to as a cocoon, in which the ion temperature remains higher than the electron temperature for longer than 100 Myr. Electrons in the cocoon are continuously heated by ions through Coulomb collisions, and the electron temperature thus remains at Te > 10^9 K in the cocoon. X-ray emissions from the cocoon are weak because the electron number density is low. Meanwhile, soft X-rays are emitted from the shocked intracluster medium surrounding the cocoon. Mixing of the jet plasma and the shocked intracluster medium through the Kelvin--Helmholtz instability at the interface enhances X-ray emissions around the contact discontinuity between the cocoon and shocked intracluster medium.
△ Less
Submitted 2 March, 2020;
originally announced March 2020.
-
Magneto Rotational Instability in Magnetized AGN Tori
Authors:
Yuki Kudoh,
Keiichi Wada
Abstract:
It is widely believed that in active galactic nuclei (AGNs) a supermassive black hole with an accretion disk is surrounded by an optically and geometrically thick torus at sub-parsec scale. However, it is not clear how is the mass supply toward the central engine caused and how it is related with the internal structures of the tori. The magnetic field in the tori may contribute to the accretion pr…
▽ More
It is widely believed that in active galactic nuclei (AGNs) a supermassive black hole with an accretion disk is surrounded by an optically and geometrically thick torus at sub-parsec scale. However, it is not clear how is the mass supply toward the central engine caused and how it is related with the internal structures of the tori. The magnetic field in the tori may contribute to the accretion process via the magneto-rotational instability (MRI). Using global three dimensional magnetohydrodynamic (MHD) simulations taking the effects of X-ray heating and radiative cooling into account studied the numerical resolution for azimuthal direction for MRI driving. We found that a strongly magnetized disk consisted of a cold ($< 10^3$ K) and warm ($10^4$ K) gas is developed in about 30 rotational periods. We also found in high resolution model that the mean azimuthal magnetic fields reverse their direction quasi-periodically. We confirmed that the typical wave length of the MRI should be resolved with a least 20 azimuthal grid cells.
△ Less
Submitted 20 December, 2018; v1 submitted 20 December, 2018;
originally announced December 2018.
-
Cosmic Magnetism in Centimeter and Meter Wavelength Radio Astronomy
Authors:
Takuya Akahori,
Hiroyuki Nakanishi,
Yoshiaki Sofue,
Yutaka Fujita,
Kiyotomo Ichiki,
Shinsuke Ideguchi,
Osamu Kameya,
Takahiro Kudoh,
Yuki Kudoh,
Mami Machida,
Yoshimitsu Miyashita,
Hiroshi Ohno,
Takeaki Ozawa,
Keitaro Takahashi,
Motokazu Takizawa,
Dai G. Yamazaki
Abstract:
Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on…
▽ More
Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography.
△ Less
Submitted 7 September, 2017;
originally announced September 2017.
-
Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications
Authors:
Yosuke Matsumoto,
Yuta Asahina,
Yuki Kudoh,
Tomohisa Kawashima,
Jin Matsumoto,
Hiroyuki R. Takahashi,
Takashi Minoshima,
Seiji Zenitani,
Takahiro Miyoshi,
Ryoji Matsumoto
Abstract:
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 magne…
▽ More
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.
△ Less
Submitted 18 August, 2019; v1 submitted 6 November, 2016;
originally announced November 2016.
-
Approximate Riemann Solvers for the Cosmic Ray Magnetohydrodynamical Equations
Authors:
Yuki Kudoh,
Tomoyuki Hanawa
Abstract:
We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully conservation form, the CR energy equation is replaced with the CR "number" conservation, where the CR number density is defined as the three fourths power of the CR ener…
▽ More
We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully conservation form, the CR energy equation is replaced with the CR "number" conservation, where the CR number density is defined as the three fourths power of the CR energy density. The former contains an extra source term, while latter does not. An approximate Riemann solver is derived from the CR MHD equations in the fully conservation form. Based on the analysis, we propose a numerical scheme of which solutions satisfy the Rankine-Hugoniot relation at any shock. We demonstrate that it reproduces the Riemann solution derived by Pfrommer et al. (2006) for a 1D CR hydrodynamic shock tube problem. We compare the solution with those obtained by solving the CR energy equation. The latter solutions deviate from the Riemann solution seriously, when the CR pressure dominates over the gas pressure in the post-shocked gas. The former solutions converge to the Riemann solution and are of the second order accuracy in space and time. Our numerical examples include an expansion of high pressure sphere in an magnetized medium. Fast and slow shocks are sharply resolved in the example. We also discuss possible extension of the CR MHD equations to evaluate the average CR energy.
△ Less
Submitted 6 October, 2016; v1 submitted 10 August, 2016;
originally announced August 2016.
-
Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography: Japanese SKA Cosmic Magnetism Science
Authors:
Takuya Akahori,
Yutaka Fujita,
Kiyotomo Ichiki,
Shinsuke Ideguchi,
Takahiro Kudoh,
Yuki Kudoh,
Mami Machida,
Hiroyuki Nakanishi,
Hiroshi Ohno,
Takeaki Ozawa,
Keitaro Takahashi,
Motokazu Takizawa
Abstract:
Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides n…
▽ More
Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides numerous observation points in the cosmic space. Furthermore, the SKA should be designed to facilitate wideband polarimetry so as to allow us to examine sightline structures of magnetic fields by means of depolarization and Faraday Tomography. The SKA will be able to drive cosmic magnetism of the interstellar medium, the Milky Way, galaxies, AGN, galaxy clusters, and potentially the cosmic web which may preserve information of the primeval Universe. The Japan SKA Consortium (SKA-JP) Magnetism Science Working Group (SWG) proposes the project "Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography", which contains ten scientific use cases.
△ Less
Submitted 7 March, 2016; v1 submitted 7 March, 2016;
originally announced March 2016.