Search

NinjaSat is an X-ray CubeSat designed for agile, long-term continuous observations of bright X-ray sources, with the size of 6U ($100\times200\times300$ mm$^3$) and a mass of 8 kg. NinjaSat is capable of pointing at X-ray sources with an accuracy of less than $0^{\circ}\hspace{-1.0mm}.1$ (2$\sigma$ confidence level) with 3-axis attitude control. The satellite bus is a commercially available NanoAvionics M6P, equipped with two non-imaging gas X-ray detectors covering an energy range of 2-50 keV. A total effective area of 32 cm$^2$ at 6 keV is capable of observing X-ray sources with a flux of approximately 10$^{-10}$ erg cm$^{-2}$ s$^{-1}$. The arrival time of each photon can be tagged with a time resolution of 61 $\mu$s. The two radiation belt monitors continuously measure the fluxes of protons above 5 MeV and electrons above 200 keV trapped in the geomagnetic field, alerting the X-ray detectors when the flux exceeds a threshold. The NinjaSat project started in 2020. Fabrication of the scientific payloads was completed in August 2022, and satellite integration and tests were completed in July 2023. NinjaSat was launched into a Sun-synchronous polar orbit at an altitude of about 530 km on 2023 November 11 by the SpaceX Transporter-9 mission. After about three months of satellite commissioning and payload verification, we observed the Crab Nebula on February 9, 2024, and successfully detected the 33.8262 ms pulsation from the neutron star. With this observation, NinjaSat met the minimum success criterion and stepped forward to scientific observations as initially planned. By the end of November 2024, we successfully observed 21 X-ray sources using NinjaSat. This achievement demonstrates that, with careful target selection, we can conduct scientific observations effectively using CubeSats, contributing to time-domain astronomy., Comment: 14 pages, 17 figures

Fuzzy dark matter (FDM) is an attractive dark matter candidate composed of ultralight particles. In this paper, toward a clear understanding of the core-halo relation in the FDM halos, we consider a simple model of the soliton-halo system, in which the self-gravitating soliton core is formed in the presence of Navarro-Frenk-White (NFW) halo potential as an external field. Solving numerically the Schr\"odinger-Poisson equation, the self-gravitating soliton is obtained as a ground-state solution, which is characterized by the two key parameters, i.e., size of soliton core and its strength of self-gravity relative to those of the NFW halo. Using our soliton-halo model, we investigate the properties of soliton cores found in cosmological simulation, and the key parameters characterizing these solitons are reconstructed in a self-consistent manner. Results suggest that (1) the soliton core properties depend critically on both the self-gravity of the soliton and the external potential of the host halo, and (2) the scatter observed in the core-halo relation cannot be explained solely by the one in the halo's concentration-mass relation, as previously suggested, but also significantly influenced by intrinsic features of the soliton core, potentially arising from local dynamics at the halo center. We also demonstrate that the FDM mass can be reconstructed from the simulation data characterizing the halo density profile, providing a basis for applying the model to observational studies., Comment: 30 pages, 17 figures

In February and March 2024, a series of many Type I X-ray bursts from the accreting neutron star SRGA J144459.2$-$604207, which has been identified by multiple X-ray satellites, with the first reports coming from INTEGRAL and NinjaSat. These observations reveal that after exhibiting very regular behavior as a ``clocked'' burster, the peak luminosity of the SRGA J144459.2$-$604207 X-ray bursts shows a gradual decline. The observed light curves exhibit a short plateau feature, potentially with a double peak, followed by a rapid decay in the tail-features unlike those seen in previously observed clocked bursters. In this study, we calculate a series of multizone X-ray burst models with various compositions of accreted matter, specifically varying the mass fractions of hydrogen ($X$), helium ($Y$), and heavier CNO elements or metallicity ($Z_{\rm CNO}$). We demonstrate that a model with higher $Z_{\rm CNO}$ and/or lower $X/Y$ compared to the solar values can reproduce the observed behavior of SRGA J144459.2$-$604207. Therefore, we propose that this new X-ray burster is likely the first clocked burster with non-solar elemental compositions. Moreover, based on the X-ray burst light curve morphology in the decline phase observed by NinjaSat, a He-enhanced model with $X/Y \approx 1.5$ seems preferred over high-metallicity cases. We also give a brief discussion on the implications for the neutron star mass, binary star evolution, inclination angle, and the potential for a high-metallicity scenario, the last of which is closely related to the properties of the hot CNO cycle., Comment: 7 pages, 5 figures, accepted for publication in PASJ Letter

The CubeSat X-ray observatory NinjaSat was launched on 2023 November 11 and has provided opportunities for agile and flexible monitoring of bright X-ray sources. On 2024 February 23, the NinjaSat team started long-term observation of the new X-ray source SRGA J144459.2$-$604207 as the first scientific target, which was discovered on 2024 February 21 and recognized as the sixth clocked X-ray burster. Our 25-day observation covered almost the entire decay of this outburst from two days after the peak at $\sim$100 mCrab on February 23 until March 18 at a few mCrab level. The Gas Multiplier Counter onboard NinjaSat successfully detected 12 Type-I X-ray bursts with a typical burst duration of $\sim$20 s, shorter than other clocked burster systems. As the persistent X-ray emission declined by a factor of five, X-ray bursts showed a notable change in its morphology: the rise time became shorter from 4.4(7) s to 0.3(3) s (1$\sigma$ errors), and the peak amplitude increased by 44%. The burst recurrence time $\Delta t_{\rm rec}$ also became longer from 2 hr to 10 hr, following the relation of $\Delta t_{\rm rec} \propto F_{\rm per}^{-0.84}$, where $F_{\rm per}$ is the persistent X-ray flux, by applying a Markov chain Monte Carlo method. The short duration of bursts is explained by the He-enhanced composition of accretion matter and the relation between $\Delta t_{\textrm{rec}}$ and $F_{\rm per}$ by a massive neutron star. This study demonstrated that CubeSat pointing observations can provide valuable astronomical X-ray data., Comment: 7 pages, 5 figures, submitted to PASJ Letter

Superradiant instability can form clouds around rotating black holes (BHs) composed of ultralight bosonic fields, such as axions. A BH with such a cloud in a binary system exhibits rich phenomena, and gravitational waves (GWs) from the BH merger provide a means to probe axions. For the first time, we study the evolution of axion clouds in a binary system during the inspiral phase, including axion self-interaction effects. When the self-interaction is significant, unlike in the negligible case, two types of clouds coexist through mode coupling. We examine the evolution of the system considering the effects of dissipation caused by both self-interaction and tidal interaction. For tidal interaction, in addition to the processes of emission to infinity and absorption by the BH, indirect emission via transitions (both resonant and off-resonant) is also considered as a second-order perturbation. Our results demonstrate that the signatures of axion self-interaction are imprinted in the modification of the GW phase. Furthermore, we find the possibility of a dynamical instability called bosenova during the binary inspiral phase., Comment: 19 pages, 13 figures, 1 table

Gravitational waves from condensates of ultra-light particles, such as axion, around rotating black holes are a promising probe to search for unknown physics. For this purpose, we need to characterize the signal to detect the gravitational waves, which requires tracking the evolution of the condensates, including various effects. The axion self-interaction causes the non-linear coupling between the superradiant modes, resulting in complicated branching of evolution. Most studies so far have considered evolution under the non-relativistic approximation or the two-mode approximation. In this paper, we numerically investigate the evolution of the axion condensate without these approximations, taking higher multipole modes into account. We also investigate the possible signature in gravitational waves from the condensate. We show that the higher multipole modes are excited, leading to the gravitational wave signal by the transition of the axion between different levels. The most prominent signal of gravitational waves arises from the transition between modes with their angular quantum numbers different by two. The gravitational wave signal is emitted in the deci-Hz band for stellar mass black holes, which might be observable with the future gravitational wave detectors., Comment: 19 pages, 17 figures

Superradiant instability of rotating black holes (BHs) leads to the formation of a cloud of ultralight bosons, such as axions. When the BH with the cloud belongs to a binary system and is in an inspiraling orbit, the resonant transition between the axion's bound states can occur. We study the history of the evolution of the binary system accompanying the cloud composed of the fastest growing mode, and its impact on the observational signatures, especially for small mass ratio cases. In this case, the hyperfine resonance, which has a very small resonance frequency, is relevant. Therefore, due to the long timescale, we should take into account the decaying process of axions in the transition destination mode, the backreaction to the orbital motion and the central BH, and gravitational emission from the cloud. We present a formulation to examine the evolution of the system around the resonance and useful expressions for the analysis. As a result, we found the mass of the cloud that can remain after the resonance is, at most, about $10^{-5}$ of the central BH. The maximum remaining cloud mass is achieved when the mass ratio of the binary is $q\sim10^{-3}$. In addition, we show that the resonant transition hardly changes the BH mass and spin distribution, while the associated modification of the gravitational wave frequency evolution when the binary pass through the resonance can be a signature of the presence of the cloud., Comment: 15 pages, 13 figures, 1 table

Ultra-light particles, such as axions, form a macroscopic condensate around a highly spinning black hole by the superradiant instability. Due to its macroscopic nature, the condensate opens the possibility of detecting the axion through gravitational wave observations. However, the precise evolution of the condensate must be known for the actual detection. For future observation, we numerically study the influence of the self-interaction, especially interaction between different modes, on the evolution of the condensate in detail. First, we focus on the case when condensate starts with the smallest possible angular quantum number. For this case, we perform the non-linear calculation and show that the dissipation induced by the mode interaction is strong enough to saturate the superradiant instability, even if the secondary cloud starts with quantum fluctuations. Our result indicates that explosive phenomena such as bosenova do not occur in this case. We also show that the condensate settles to a quasi-stationary state mainly composed of two modes, one with the smallest angular quantum number for which the superradiant instability occurs and the other with the adjacent higher angular quantum number. We also study the case when the condensate starts with the dominance of the higher angular quantum number. We show that the dissipation process induced by the mode coupling does not occur for small gravitational coupling. Therefore, bosenova might occur in this case., Comment: 29 pages, 25 figures, 1 table

Ultra light axion fields, motivated by the string theory, form a large condensate (axion cloud) around rotating black holes through superradiant instability. Several effects due to the axion cloud, such as the spin-down of black holes and the emission of monochromatic gravitational waves, open a new window to search for axions by astrophysical observations. When the axion self-interaction is considered, the evolution of cloud is altered significantly, and an explosive phenomenon called bosenova can happen. Thus, it is necessary to understand the precise evolution of self-interacting clouds for the detection of axions by astrophysical observations. In this paper, we propose a new method to track the whole process of the growth of self-interacting axion clouds employing the adiabatic approximation. We emphasize that our method relies neither on the non-relativistic approximation nor on perturbative treatment of the self-interaction, which is often used in literature. Our main finding is that the evolution of cloud in the strongly self-interacting regime depends on the strength of the gravitational coupling between the axion and the black hole. For a large coupling, the cloud evolves into a quasi-stationary state where the superradiant energy gain is balanced with the energy dissipation to infinity by the self-interaction. On the other hand, when one decreases the size of coupling, clouds become unstable at some energy, which would be interpreted as the onset of bosenova., Comment: 22pages, 15 figures, 1 table

Ultralight scalar fields such as axions can form clouds around rotating black holes (BHs) by the superradiant instability. It is important to consider the evolution of clouds associated with BH binaries for the detectability of the presence of clouds through gravitational wave signals and observations of the mass and spin parameters of BHs. We re-examine the impact on the axion cloud due to the tidal perturbation from the companion in a binary system taking into account the following points. First, we study the influence of higher-multipole moments. Second, we consider the backreaction due to the angular momentum transfer between the cloud and the orbital motion. This angular momentum transfer further causes the backreaction to the hyperfine split through the change in geometry. Finally, we calculate the particle number flux to the infinity induced by the tidal interaction. As a result, we found that the scalar field is not reabsorbed by the BH. Instead, the scalar particles are radiated away to evaporate during the inspiral, irrespectively of the direction of the orbital motion, for almost equal mass binaries., Comment: 21 pages, 7 figures

Gravitational wave observation has the potential of probing ultralight bosonic fields such as axion. Axion forms a cloud around a rotating black hole (BH) by superradiant instability and should affect the gravitational waveform from binary BHs. On the other hand, considering the cloud associated with a BH in a binary system, tidal interaction depletes the cloud in some cases during the inspiral phase. We made the exhaustive study of cloud depletion numerically in a wide parameter range for equal mass binaries, assuming only the quadrupolar tidal perturbation is at work. We found that clouds can avoid disappearing due to the tidal effect only when $l=1$ mode is the fastest growing mode and when the binary orbit is counter-rotating in the non-relativistic parameter region., Comment: 15 pages, 4 figures

There are strong interests in considering ultra-light scalar fields (especially axion) around a rapidly rotating black hole because of the possibility of observing gravitational waves from axion condensate (axion cloud) around black holes. Motivated by this consideration, we propose a new method to study the dynamics of an ultra-light scalar field with self-interaction around a rapidly rotating black hole, which uses the dynamical renormalization group method. We find that for relativistic clouds, saturation of the superradiant instability by the scattering of the axion due to the self-interaction does not occur in the weakly non-linear regime when we consider the adiabatic growth of the cloud from a single superradiant mode. This may suggest that for relativistic axion clouds, an explosive phenomenon called the Bosenova may inevitably happen, at least once in its evolutionary history., Comment: 31 pages, 7 figures, 2 tables

A. Sannami constructed an example of the differentiable Cantor set embedded in the real line whose difference set has a positive measure. In this paper, we generalize the definition of the difference sets for sets of the two dimensional Euclidean space as the sets of vectors between two sets, and estimate their measures. For the quadratic map Q_c(z)=z^2+c, we obtain that the measure of the difference set of its Julia set vanishes if |c|>3+sqrt{3}., Comment: 7 pages, 2 figures

Spin-state crossover beyond a conventional ligand-field theory has been a fundamental issue in condensed matter physics. Here, we report microscopic observations of spin states and low-energy dynamics through orbital-resolved NMR spectroscopy in the prototype compound LaCoO$_3$. The $^{59}$Co NMR spectrum shows the preserved crystal symmetry across the crossover, inconsistent with $d$ orbital ordering due to the Jahn-Teller distortion. The orbital degeneracy results in a pseudospin ($\tilde{J} = 1$) excited state with an orbital moment observed as $^{59}$Co hyperfine coupling tensors. We found that the population of the excited state evolves above the heart crossover temperature. The crossover involves critical spin-state fluctuations emerging under the magnetic field. These results suggest that the spin-state crossover can be mapped into a statistical problem, analogous to the supercritical liquid in liquid-gas transition., Comment: 5 pages, 5 figures

We propose a mechanism for quasi-periodic oscillations of both coronal mass ejections (CMEs) and flare loops as related to magnetic reconnection in eruptive solar flares. We perform two-dimensional numerical MHD simulations of magnetic flux rope eruption, with three different values of the global Lundquist number. In the low Lundquist number run, no oscillatory behavior is found. In the moderate Lunquist number run, on the other hand, quasi-periodic oscillations are excited both at the bottom of the flux rope and at the flare loop-top. In the high Lundquist number run, quasi-periodic oscillations are also excited; in the meanwhile, the dynamics become turbulent due to the formation of multiple plasmoids in the reconnection current sheet. In high and moderate Lundquist number runs, thin reconnection jet collide with the flux rope bottom or flare loop-top and dig them deeply. Steep oblique shocks are formed as termination shocks where reconnection jet is bent (rather than decelerated) in horizontal direction, resulting in supersonic back-flows. The structure becomes unstable, and quasi-periodic oscillation of supersonic back-flows appear at locally confined high-beta region at both the flux rope bottom and flare loop-top. We compare the observational characteristics of quasi-periodic oscillations in erupting flux ropes, post-CME current sheets, flare ribbons and light curves, with corresponding dynamical structures found in our simulation., Comment: 30 pages, 20 figures, accepted for publication in ApJ

Fast interplanetary coronal mass ejections (interplanetary CMEs, or ICMEs) are the drivers of strongest space weather storms such as solar energetic particle events and geomagnetic storms. The connection between space weather impacting solar wind disturbances associated with fast ICMEs at Earth and the characteristics of causative energetic CMEs observed near the Sun is a key question in the study of space weather storms as well as in the development of practical space weather prediction. Such shock-driving fast ICMEs usually expand at supersonic speed during the propagation, resulting in the continuous accumulation of shocked sheath plasma ahead. In this paper, we propose the "sheath-accumulating propagation" (SAP) model that describe the coevolution of the interplanetary sheath and decelerating ICME ejecta by taking into account the process of upstream solar wind plasma accumulation within the sheath region. Based on the SAP model, we discussed (1) ICME deceleration characteristics, (2) the fundamental condition for fast ICME at Earth, (3) thickness of interplanetary sheath, (4) arrival time prediction and (5) the super-intense geomagnetic storms associated with huge solar flares. We quantitatively show that not only speed but also mass of the CME are crucial in discussing the above five points. The similarities and differences among the SAP model, the drag-based model and the`snow-plough' model proposed by \citet{tappin2006} are also discussed., Comment: 20 pages, 5 figures, accepted for publication in ApJL

Flare associated coronal shock waves sometimes interact with solar prominences leading to large amplitude prominence oscillations. Such prominence activation gives us unique opportunity to track time evolution of shock-cloud interaction in cosmic plasmas. Although the dynamics of interstellar shock-cloud interaction is extensively studied, coronal shock-solar prominence interaction is rarely studied in the context of shock-cloud interaction. Associated with X5.4 class solar flare occurred on 7 March, 2012, a globally propagated coronal shock wave interacted with a polar prominence leading to large amplitude prominence oscillation. In this paper, we studied bulk acceleration and excitation of internal flow of the shocked prominence using three-dimensional MHD simulations. We studied eight magnetohydrodynamic (MHD) simulation runs with different mass density structure of the prominence, and one hydrodynamic simulation run, and compared the result. In order to compare observed motion of activated prominence with corresponding simulation, we also studied prominence activation by injection of triangular shaped coronal shock. We found that magnetic tension force mainly accelerate (and then decelerate) the prominence. The internal flow, on the other hand, is excited during the shock front sweeps through the the prominence and damps almost exponentially. We construct phenomenological model of bulk momentum transfer from shock to the prominence, which agreed quantitatively with all the simulation results. Based on the phenomenological prominence-activation model, we diagnosed physical parameters of coronal shock wave. The estimated energy of the coronal shock is several percent of total energy released during the X5.4 flare., Comment: 48 pages, 20 figures, accepted for publication in ApJ

In order to discuss the potential impact of solar 'superflares' on space weather, we investigated statistical relations among energetic proton peak flux with energy higher than $ 10 \rm MeV$ ($F_p$), CME speed near the Sun ($V_{CME}$) obtained by {\it SOHO}/LASCO coronagraph and flare soft X-ray peak flux in 1-8\AA band ($F_{SXR}$) during 110 major solar proton events (SPEs) recorded from 1996 to 2014. The linear regression fit results in the scaling relations $V_{CME} \propto F_{SXR}^\alpha$, $F_p\propto F_{SXR}^\beta$ and $F_p\propto V_{CME}^\gamma$ with $\alpha = 0.30\pm 0.04$, $\beta = 1.19 \pm 0.08$ and $\gamma = 4.35 \pm 0.50$, respectively. On the basis of simple physical assumptions, on the other hand, we derive scaling relations expressing CME mass ($M_{CME}$), CME speed and energetic proton flux in terms of total flare energy ($E_{flare}$) as, $M_{CME}\propto E_{flare}^{2/3}$, $V_{CME}\propto E_{flare}^{1/6}$ and $F_{p}\propto E_{flare}^{5/6}\propto V_{CME}^5$, respectively. We then combine the derived scaling relations with observation, and estimated the upper limit of $V_{CME}$ and $F_p$ to be associated with possible solar superflares., Comment: 17 pages, 5 figures, accepted for publication in Astrophysical Journal Letters

An X5.4 class flare occurred in active region (AR) NOAA11429 on 2012 March 7. The flare was associated with very fast coronal mass ejection (CME) with its velocity of over 2500 km/s. In the images taken with STEREO-B/COR1, a dome-like disturbance was seen to detach from expanding CME bubble and propagated further. A Type-II radio burst was also observed at the same time. On the other hand, in EUV images obtained by SDO/AIA, expanding dome-like structure and its foot print propagating to the north were observed. The foot print propagated with its average speed of about 670 km/s and hit a prominence located at the north pole and activated it. While the activation, the prominence was strongly brightened. On the basis of some observational evidence, we concluded that the foot print in AIA images and the ones in COR1 images are the same, that is MHD fast mode shock front. With the help of a linear theory, the fast mode mach number of the coronal shock is estimated to be between 1.11 and 1.29 using the initial velocity of the activated prominence. Also, the plasma compression ratio of the shock is enhanced to be between 1.18 and 2.11 in the prominence material, which we consider to be the reason of the strong brightening of the activated prominence. The applicability of linear theory to the shock problem is tested with nonlinear MHD simulation.