Properties of plasmoid-dominated turbulent reconnection in a low-$\beta$ background plasma are investigated by resistive magnetohydrodynamic (MHD) simulations. In the $\beta_{\rm in} < 1$ regime, where $\beta_{\rm in}$ is plasma $\beta$ in the inflow region, the reconnection site is dominated by shocks and shock-related structures and plasma compression is significant. The effective reconnection rate increases from $0.01$ to $0.02$ as $\beta_{\rm in}$ decreases. We hypothesize that plasma compression allows faster reconnection rate, and then we estimate a speed-up factor, based on a compressible MHD theory. We validate our prediction by a series of MHD simulations. These results suggest that the plasmoid-dominated reconnection can be twice faster than expected in the $\beta \ll 1$ environment in a solar corona.

We develop a new numerical scheme for ideal magnetohydrodynamic (MHD) simulations, which is robust against one- and multi-dimensional shocks, and is accurate for low Mach number flows and discontinuities. The scheme belongs to a family of the advection upstream splitting method employed in computational aerodynamics, and it splits the inviscid flux in MHD equations into advection, pressure, and magnetic tension parts, and then individually evaluates mass, pressure, and magnetic tension fluxes at the interface of a computational cell. The mass flux is designed to avoid numerical shock instability in multidimension, while preserving contact discontinuity. The pressure flux possesses a proper scaling for low Mach number flows, allowing reliable simulations of nearly incompressible flows. The magnetic tension flux is built to be consistent with the HLLD approximate Riemann solver to preserve rotational discontinuity. We demonstrate various benchmark tests to verify the novel performance of the scheme. Our results indicate that the scheme must be a promising tool to tackle astrophysical systems that include both low and high Mach number flows, as well as magnetic field inhomogeneities.

We investigate the orbital motions of dwarf spheroidal galaxies (dSphs) in the halo of the Milky Way (MW) to understand their possible effects on the diversity of the star formation histories seen in these MW satellites. In this work, we explicitly consider a time-varying gravitational potential due to the growth of the MW's dark halo mass to calculate the long-term orbital evolutions of the dSphs, guided with {\it Gaia} DR2 proper motions, over the past 13.5 billion years. We find that the infall time of a satellite, defined at which the galaxy first crosses within the growing virial radius of the MW's halo, coincides well with the time when the star formation rate (SFR) is peaked for the sample of classical dSphs. On the other hand, ultra-faint dSphs already finished their SF activity prior to their infall times as already suggested in previous works, but there is a signature that their earlier SF histories are affected by interaction with the growing MW's halo to some extent. We also find, for classical dSphs, that the relative fraction of stars formed after the peak of the SFR to the current stellar mass is smaller for the smaller pericentric radius of the galaxy at its first infall. These results suggest that the infalling properties of the dSphs into the MW and the resultant environmental effects such as ram-pressure stripping and/or tidal disturbance in the MW's dark halo containing hot gas play important roles in their star formation histories.

A nonlinear force-free field (NLFFF) extrapolation is widely used to reconstruct the three-dimensional magnetic field in the solar corona from the observed photospheric magnetic field. However, the pressure gradient and gravitational forces are ignored in the NLFFF model, even though the photospheric and chromospheric magnetic fields are not in general force-free. Here we develop a magnetohydrodynamic (MHD) relaxation method that reconstructs the solar atmospheric (chromospheric and coronal) magnetic field as a non-force-free magnetic field (NFFF) in magnetohydrostatic equilibrium where the Lorentz, pressure gradient, and gravitational forces are balanced. The system of basic equations for the MHD relaxation method is derived, and mathematical properties of the system are investigated. A robust numerical solver for the system is constructed based on the modern high-order shock capturing scheme. Two-dimensional numerical experiments that include the pressure gradient and gravitational forces are also demonstrated.

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.

We report on a measurement of the size of charge clouds produced by X-ray photons in X-ray SOI (Silicon-On-Insulator) pixel sensor named XRPIX. We carry out a beam scanning experiment of XRPIX using a monochromatic X-ray beam at 5.0 keV collimated to $\sim 10$ $\mu$m with a 4-$\mu$m$\phi$ pinhole, and obtain the spatial distribution of single-pixel events at a sub-pixel scale. The standard deviation of charge clouds of 5.0 keV X-ray is estimated to be $\sigma_{\rm cloud} = 4.30 \pm 0.07$ $\mu$m. Compared to the detector response simulation, the estimated charge cloud size is well explained by a combination of photoelectron range, thermal diffusion, and Coulomb repulsion. Moreover, by analyzing the fraction of multi-pixel events in various energies, we find that the energy dependence of the charge cloud size is also consistent with the simulation.

We design a conservative finite difference scheme for ideal magnetohydrodynamic simulations that attains high-order accuracy, shock-capturing, and divergence-free condition of the magnetic field. The scheme interpolates pointwise physical variables from computational nodes to midpoints through a high-order nonlinear weighted average. The numerical flux is evaluated at the midpoint by a multi-state approximate Riemann solver for correct upwinding, and its spatial derivative is approximated by a high-order linear central difference to update the variables with designed order of accuracy and conservation. The magnetic and electric fields are defined at staggered grid points employed in the Constrained Transport (CT) method by Evans & Hawley (1988). We propose a new CT variant, in which the staggered electric field is evaluated so as to be consistent with the base one-dimensional Riemann solver and the staggered magnetic field is updated to be divergence-free as designed high-order finite difference representation. We demonstrate various benchmark tests to measure the performance of the present scheme. We discuss the effect of the choice of interpolation methods, Riemann solvers, and the treatment for the divergence-free condition on the quality of numerical solutions in detail.

A new Riemann solver is presented for the ideal magnetohydrodynamics (MHD) equations with the so-called Boris correction. The Boris correction is applied to reduce wave speeds, avoiding an extremely small timestep in MHD simulations. The proposed Riemann solver, Boris-HLLD, is based on the HLLD solver. As done by the original HLLD solver, (1) the Boris-HLLD solver has four intermediate states in the Riemann fan when left and right states are given, (2) it resolves the contact discontinuity, Alfv\'en waves, and fast waves, and (3) it satisfies all the jump conditions across shock waves and discontinuities except for slow shock waves. The results of a shock tube problem indicate that the scheme with the Boris-HLLD solver captures contact discontinuities sharply and it exhibits shock waves without any overshoot when using the minmod limiter. The stability tests show that the scheme is stable when $|u| \lesssim 0.5c$ for a low Alfv\'en speed ($V_A \lesssim c$), where $u$, $c$, and $V_A$ denote the gas velocity, speed of light, and Alfv\'en speed, respectively. For a high Alfv\'en speed ($V_A \gtrsim c$), where the plasma beta is relatively low in many cases, the stable region is large, $|u| \lesssim (0.6-1) c$. We discuss the effect of the Boris correction on physical quantities using several test problems. The Boris-HLLD scheme can be useful for problems with supersonic flows in which regions with a very low plasma beta appear in the computational domain.

We have been developing event-driven SOI Pixel Detectors, named `XRPIX' (X-Ray soiPIXel) based on the silicon-on-insulator (SOI) pixel technology, for the future X-ray astronomical satellite with wide band coverage from 0.5 keV to 40 keV. XRPIX has event trigger output function at each pixel to acquire a good time resolution of a few $\mu \rm s$ and has Correlated Double Sampling function to reduce electric noises. The good time resolution enables the XRPIX to reduce Non X-ray Background in the high energy band above 10\,keV drastically by using anti-coincidence technique with active shield counters surrounding XRPIX. In order to increase the soft X-ray sensitivity, it is necessary to make the dead layer on the X-ray incident surface as thin as possible. Since XRPIX1b, which is a device at the initial stage of development, is a front-illuminated (FI) type of XRPIX, low energy X-ray photons are absorbed in the 8 $\rm \mu$m thick circuit layer, lowering the sensitivity in the soft X-ray band. Therefore, we developed a back-illuminated (BI) device XRPIX2b, and confirmed high detection efficiency down to 2.6 keV, below which the efficiency is affected by the readout noise. In order to further improve the detection efficiency in the soft X-ray band, we developed a back-illuminated device XRPIX3b with lower readout noise. In this work, we irradiated 2--5 keV X-ray beam collimated to 4 $\rm \mu m \phi$ to the sensor layer side of the XRPIX3b at 6 $\rm \mu m$ pitch. In this paper, we reported the uniformity of the relative detection efficiency, gain and energy resolution in the subpixel level for the first time. We also confirmed that the variation in the relative detection efficiency at the subpixel level reported by Matsumura et al. has improved.

In low earth orbit, there are many cosmic rays composed primarily of high energy protons. These cosmic rays cause surface and bulk radiation effects, resulting in degradation of detector performance. Quantitative evaluation of radiation hardness is essential in development of X-ray detectors for astronomical satellites. We performed proton irradiation experiments on newly developed X-ray detectors called XRPIX based on silicon-on-insulator technology at HIMAC in National Institute of Radiological Sciences. We irradiated 6 MeV protons with a total dose of 0.5 krad, equivalent to 6 years irradiation in orbit. As a result, the gain increases by 0.2% and the energy resolution degrades by 0.5%. Finally we irradiated protons up to 20 krad and found that detector performance degraded significantly at 5 krad. With 5 krad irradiation corresponding to 60 years in orbit, the gain increases by 0.7% and the energy resolution worsens by 10%. By decomposing into noise components, we found that the increase of the circuit noise is dominant in the degradation of the energy resolution.

Nonlinear evolution of magnetic reconnection is investigated by means of magnetohydrodynamic simulations including uniform resistivity, uniform viscosity, and anisotropic thermal conduction. When viscosity exceeds resistivity (the magnetic Prandtl number Prm > 1), the viscous dissipation dominates outflow dynamics and leads to the decrease in the plasma density inside a current sheet. The low-density current sheet supports the excitation of the vortex. The thickness of the vortex is broader than that of the current for Prm > 1. The broader vortex flow more efficiently carries the upstream magnetic flux toward the reconnection region, and consequently boosts the reconnection. The reconnection rate increases with viscosity provided that thermal conduction is fast enough to take away the thermal energy increased by the viscous dissipation (the fluid Prandtl number Pr < 1). The result suggests the need to control the Prandtl numbers for the reconnection against the conventional resistive model.

Plasmoid structures in fast reconnection in low-beta plasmas are investigated by two-dimensional magnetohydrodynamic simulations. A high-resolution shock-capturing code enables us to explore a variety of shock structures: vertical slow shocks behind the plasmoid, another slow shocks in the outer-region, and the shock-reflection in the front side. The Kelvin-Helmholtz-like turbulence is also found inside the plasmoid. It is concluded that these shocks are rigorous features in reconnection in low-beta plasmas, where the reconnection jet speed or the upstream Alfven speed exceeds the sound speed.

We performed magnetohydrodynamic simulation of a formation process of coronal mass ejections (CMEs), focusing on interaction (reconnection) between an ejecting flux rope and its ambient field. We examined three cases with different ambient fields: no ambient field, and cases with dipole field of two opposite directions which are parallel and anti-parallel to that of the flux rope surface. As a result, while the flux rope disappears in the anti-parallel case, in other cases the flux ropes can evolve to CMEs and show different amounts of rotation of the flux rope. The results imply that the interaction between an ejecting flux rope and its ambient field is an important process for determining CME formation and CME orientation, and also show that the amount and direction of magnetic flux within the flux rope and the ambient field are key parameters for CME formation. Especially, the interaction (reconnection) plays a significant role to the rotation of the flux rope, with a process similar to "tilting instability" in a spheromak-type experiment of laboratory plasma.