Source author record

Takayoshi Sano

Takayoshi Sano 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

13works
7topics
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

13 published item(s)

preprint2026arXiv

Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

Recent observations by the IceCube Neutrino Observatory have revealed a significant excess of high-energy neutrinos from nearby Seyfert galaxies, such as NGC~1068, without a corresponding flux of high-energy gamma-rays. This suggests that neutrinos are produced via hadronic interactions in a region opaque to gamma-rays, likely a hot corona surrounding the central supermassive black hole. However, the mechanism responsible for accelerating the parent protons to the required energies ($\sim 100$ TeV) remains an open question. In this study, we investigate diffusive shock acceleration (DSA) in active galactic nucleus (AGN) coronae using a suite of one-dimensional Particle-in-cell (PIC) simulations spanning a broad range of plasma parameters. We find that DSA is a robust and efficient mechanism for proton acceleration, consistently channeling approximately 10\% of the shock's kinetic energy into non-thermal ions, even for shocks with sonic Mach number as low as $ M_s \approx 2$. In contrast, the efficiency of electron acceleration is highly variable and less efficient ($<1\%$) in our parameter survey. These findings provide strong, first-principles support for the hadronic models of neutrino production in AGN and offer quantitative constraints that can explain the observed gamma-ray deficit.

preprint2021arXiv

Optab: Public code for generating gas opacity tables for radiation hydrodynamics simulations

We have developed a public code, Optab, that outputs Rosseland, Planck, and two-temperature Planck mean gas opacity tables for radiation hydrodynamics simulations in astrophysics. The code is developed for modern high-performance computing, being written in Fortran 90 and using Message Passing Interface and Hierarchical Data Format, Version 5. The purpose of this work is to provide a platform on which users can generate opacity tables for their own research purposes. Therefore, the code has been designed so that a user can easily modify, change, or add opacity sources in addition to those already implemented, which include bremsstrahlung, photoionization, Rayleigh scattering, line absorption, and collision-induced absorption. In this paper, we provide details of the opacity calculations in our code and present validation tests to evaluate the performance of our code.

preprint2020arXiv

Suppression of the Richtmyer-Meshkov Instability due to a Density Transition Layer at the Interface

We have investigated the effects of a smooth transition layer at the contact discontinuity on the growth of the Richtmyer-Meshkov instability (RMI) by hydrodynamic numerical simulations and derived an empirical condition for the suppression of the instability. The transition layer has little influence on the RMI when the thickness $L$ is narrower than the wavelength of an interface modulation $λ$. However, if the transition layer becomes broader than $λ$, the perturbed velocity associated with the RMI is reduced considerably. The suppression condition is interpreted as the cases that the shock transit time through the transition layer is longer than the sound crossing time of the modulation wavelength. The fluctuation kinetic energy decreases as $L^{-p}$ with $p = 2.5$, which indicates that the growth velocity of the RMI decreases in proportion to $L^{-p/2}$ by the presence of the transition layer. This feature is found to be quite universal and appeared in a wide range of shock-interface interactions.

preprint2020arXiv

Thermonuclear Fusion Triggered by Collapsing Standing Whistler Waves in Magnetized Overdense Plasmas

Thermal fusion plasmas initiated by standing whistler waves are investigated numerically by two- and one-dimensional Particle-in-Cell simulations. When a standing whistler wave collapses due to the wave breaking of ion plasma waves, the energy of the electromagnetic waves transfers directly to the ion kinetic energy. Here we find that the ion heating by the standing whistler wave is operational even in multi-dimensional simulations of multi-ion species targets, such as deuterium-tritium (DT) ices and solid ammonia borane (H$_6$BN). The energy conversion efficiency to ions becomes as high as 15% of the injected laser energy, which depends significantly on the target thickness and laser pulse duration. The ion temperature could reach a few tens of keV or much higher if appropriate laser-plasma conditions are selected. DT fusion plasmas generated by this method must be useful as efficient neutron sources. Our numerical simulations suggest that the neutron generation efficiency exceeds 10$^9$ n/J per steradian, which is beyond the current achievements of the state-of-the-art laser experiments. The standing whistler wave heating would expand the experimental possibility for an alternative ignition design of magnetically confined laser fusion, and also for more difficult fusion reactions including the aneutronic proton-boron reaction.

preprint2019arXiv

Peta-Pascal Pressure Driven by Fast Isochoric Heating with Multi-Picosecond Intense Laser Pulse

Fast isochoric laser heating is a scheme to heat a matter with relativistic-intensity ($>$ 10$^{18}$ W/cm$^2$) laser pulse or X-ray free electron laser pulse. The fast isochoric laser heating has been studied for creating efficiently ultra-high-energy-density (UHED) state. We demonstrate an fast isochoric heating of an imploded dense plasma using a multi-picosecond kJ-class petawatt laser with an assistance of externally applied kilo-tesla magnetic fields for guiding fast electrons to the dense plasma.The UHED state with 2.2 Peta-Pascal is achieved experimentally with 4.6 kJ of total laser energy that is one order of magnitude lower than the energy used in the conventional implosion scheme. A two-dimensional particle-in-cell simulation reveals that diffusive heating from a laser-plasma interaction zone to the dense plasma plays an essential role to the efficient creation of the UHED state.

preprint2016arXiv

Spontaneous Formation of Surface Magnetic Structure from Large-scale Dynamo in Strongly-stratified Convection

We report the first successful simulation of spontaneous formation of surface magnetic structures from a large-scale dynamo by strongly-stratified thermal convection in Cartesian geometry. The large-scale dynamo observed in our strongly-stratified model has physical properties similar to those in earlier weakly-stratified convective dynamo simulations, indicating that the $α^2$-type mechanism is responsible for it. Additionally to the large-scale dynamo, we find that large-scale structures of the vertical magnetic field are spontaneously formed in the convection zone surface only for the case of strongly-stratified atmosphere. The organization of the vertical magnetic field proceeds in the upper convection zone within tens of convective turn-over time and band-like bipolar structures are recurrently-appeared in the dynamo-saturated stage. We examine possibilities of several candidates as the origin of the surface magnetic structure formation, and then suggest the existence of an as-yet-unknown mechanism for the self-organization of the large-scale magnetic structure, which should be inherent in the strongly-stratified convective atmosphere.

preprint2015arXiv

Dependence of the saturation level of magnetorotational instability on gas pressure and magnetic Prandtl number

A large set of numerical simulations of magnetohydrodynamic (MHD) turbulence induced by the magnetorotational instability (MRI) is presented. Revisiting the previous survey conducted by Sano et al. (2004), we investigate the gas pressure dependence of the saturation level. In ideal MHD simulations, the gas pressure dependence is found to be very sensitive to the choice of a numerical scheme. This is because the numerical magnetic Prandtl number varies according to the scheme as well as the pressure, which considerably affects the results. The saturation level is more sensitive to the numerical magnetic Prandtl number than the pressure. In MHD simulations with explicit viscosity and resistivity, the saturation level increases with the physical magnetic Prandtl number, and it is almost independent of the gas pressure when the magnetic Prandtl number is constant. This is indicative of the incompressible turbulence saturated by the secondary tearing instability.

preprint2014arXiv

Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $α$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $α$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

preprint2014arXiv

Long-term Evolution of Large-scale Magnetic Fields in Rotating Stratified Convection

Convective dynamo simulations are performed in local Cartesian geometry. We report the first successful simulation of a large-scale oscillatory dynamo in rigidly rotating convection without stably stratified layers. A key requirement for exciting the large-scale dynamo is a sufficiently long integration time comparable to the ohmic diffusion time. By comparing two models with and without stably stratified layers, their effect on the large-scale dynamo is also studied. The spatiotemporal evolution of the large-scale magnetic field is similar in both models. However, it is intriguing that the magnetic cycle is much shorter in the model without the stable layer than with the stable layer. This suggests that the stable layer impedes the cyclic variations of the large-scale magnetic field.

preprint2014arXiv

Mean-Field Modeling of $α^2$-Dynamo Coupled with Direct Numerical Simulations of Rigidly Rotating Convection

The mechanism of large-scale dynamos in rigidly rotating stratified convection is explored by direct numerical simulations (DNS) in Cartesian geometry. A mean-field dynamo model is also constructed using turbulent velocity profiles consistently extracted from the corresponding DNS results. By quantitative comparison between the DNS and our mean-field model, it is demonstrated that the oscillatory $α^2$ dynamo wave, excited and sustained in the convection zone, is responsible for large-scale magnetic activities such as cyclic polarity reversal and spatiotemporal migration. The results provide strong evidence that a nonuniformity of the $α$-effect, which is a natural outcome of rotating stratified convection, can be an important prerequisite for large-scale stellar dynamos, even without the $Ω$-effect.

preprint2013arXiv

Critical Magnetic Field Strength for Suppression of the Richtmyer-Meshkov Instability in Plasmas

The critical strength of a magnetic field required for the suppression of the Richtmyer-Meshkov instability (RMI) is investigated numerically by using a two-dimensional single-mode analysis. For the cases of MHD parallel shocks, the RMI can be stabilized as a result of the extraction of vorticity from the interface. A useful formula describing a critical condition for MHD RMI has been introduced, and which is successfully confirmed by the direct numerical simulations. The critical field strength is found to be largely depending on the Mach number of the incident shock. If the shock is strong enough, even low-$β$ plasmas can be subject to the growth of the RMI.

preprint2012arXiv

Magnetic Field Amplification Associated with the Richtmyer-Meshkov Instability

The amplification of a magnetic field due to the Richtmyer-Meshkov instability (RMI) is investigated by two-dimensional MHD simulations. Single-mode analysis is adopted to reveal definite relation between the nonlinear evolution of RMI and the field enhancement. It is found that an ambient magnetic field is stretched by fluid motions associated with the RMI, and the strength is amplified significantly by more than two orders of magnitude. The saturation level of the field is determined by a balance between the amplified magnetic pressure and the thermal pressure after shock passage. This effective amplification can be achieved in a wide range of the conditions for the RMI such as the Mach number of an incident shock and the density ratio at a contact discontinuity. The results suggest that the RMI could be a robust mechanism of the amplification of interstellar magnetic fields and cause the origin of localized strong fields observed at the shock of supernova remnants.

preprint2007arXiv

Dead Zone Formation and Nonsteady Hyperaccretion in Collapsar Disks : A Possible Origin of Short-Term Variability in the Prompt Emission of Gamma-Ray Bursts

The central engine of gamma-ray bursts (GRBs) is believed to be a hot and dense disk with hyperaccretion onto a few solar-mass black hole. We investigate where the magnetorotational instability (MRI) actively operates in the hyperaccretion disk, which can cause angular momentum transport in the disk. The inner region of hyperaccretion disks can be neutrino opaque, and the energy- and momentum-transport by neutrinos could affect the growth of the MRI significantly. Assuming reasonable disk models and a weak magnetic field $B \lesssim 10^{14} \rm{G}$, it is found that the MRI is strongly suppressed by the neutrino viscosity in the inner region of hyperaccretion disks. On the other hand, the MRI can drive active MHD turbulence in the outer neutrino-transparent region regardless of the field strength. This suggests that the baryonic matter is accumulated into the inner dead zone where the MRI grows inactively and the angular momentum transport is inefficient. When the dead zone gains a large amount of mass and becomes gravitationally unstable, intense mass accretion onto the central black hole would occur episodically through the gravitational torque. This process can be a physical mechanism of the short-term variability in the prompt emission of GRBs. Finally, the origin of flaring activities in the X-ray afterglow is predicted in the context of our episodic accretion scenario.