We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$\alpha$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$\alpha$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$\alpha$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$\alpha$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

This paper presents a data processing algorithm with machine learning for polarization extraction and event selection applied to photoelectron track images taken with X-ray polarimeters. The method uses a convolutional neural network (CNN) classification to predict the azimuthal angle and 2-D position of the initial photoelectron emission from a 2-D track image projected along the X-ray incident direction. Two CNN models are demonstrated with data sets generated by a Monte Carlo simulation: one has a commonly used loss function calculated by the cross entropy and the other has an additional loss term to penalize nonuniformity for an unpolarized modulation curve based on the $H$-test, which is used for periodic signal search in X-ray/$\gamma$-ray astronomy. The modulation curve calculated by the former model with unpolarized data has several irregular features, which can be canceled out by unfolding the angular response or simulating the detector rotation. On the other hand, the latter model can predict a flat modulation curve with a residual systematic modulation down to $\lesssim1$%. Both models show almost the same modulation factors and position accuracy of less than 2 pixel (or 240 $\mu$m) for all four test energies of 2.7, 4.5, 6.4, and 8.0 keV. In addition, event selection is performed based on probabilities from the CNN to maximize the polarization sensitivity considering a trade-off between the modulation factor and signal acceptance. The developed method with machine learning improves the polarization sensitivity by 10-20%, compared to that determined with the image moment method developed previously.

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.

The typical accreting neutron star, Aquila X-1, was observed with Suzaku seven times in the decay phase of an outburst in 2007 September-October. Among them, the second to the fourth observations were performed 10 to 22 days after the outburst peak, when the source was in the hard state with a luminosity of 2x10^36 erg/sec. A unified spectral model for this type of objects approximately reproduced the 0.8--100~keV spectra obtained in these 3 observations. However, the spectra all exhibited an enigmatic hump-like excess around 30 keV, above the hard X-ray continuum which is interpreted as arising via Comptonization. The excess feature was confirmed to be significant against statistical and systematic uncertainties. It was successfully represented by a Gaussian centered at ~32 keV, with a width (sigma) of ~6 keV and an equivalent width of ~8.6 keV. Alternatively, the feature can also be explained by a recombination edge model, that produces a quasi-continuum above an edge energy of ~27 keV with an electron temperature of ~11 keV and an equivalent width of ~6.3 keV. These results are discussed in the context of atomic features of heavy elements synthesized via rapid-proton capture process during thermonuclear flashes.

We present the third MAXI/GSC catalog in the high Galactic-latitude sky ($|b| > 10^\circ$) based on the 7-year data from 2009 August 13 to 2016 July 31, complementary to that in the low Galactic-latitude sky ($|b| < 10^\circ$; Hori et al. 2018). We compile 682 sources detected at significances of $s_{\rm D,4-10~keV} \geq 6.5$ in the 4--10 keV band. A two-dimensional image fit based on the Poisson likelihood algorithm ($C$-statistics) is adopted for the detections and constraints on their fluxes and positions. The 4--10 keV sensitivity reaches $\approx 0.48$ mCrab, or $\approx 5.9 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, over the half of the survey area. Compared with the 37-month catalog (Hiroi et al. 2013), which adopted a threshold of $s_{\rm D,4-10~keV} \geq 7$, the source number increases by a factor of $\sim$1.4. The fluxes in the 3--4 keV and 10--20 keV bands are further estimated, and hardness ratios (HRs) are calculated using the 3--4 keV, 4--10 keV, 3--10 keV, and 10--20 keV band fluxes. We also make the 4--10 keV lightcurves in one year bins for all the sources and characterize their variabilities with an index based on a likelihood function and the excess variance. Possible counterparts are found from five major X-ray survey catalogs by Swift, Uhuru, RXTE, XMM-Newton, and ROSAT, and an X-ray galaxy-cluster catalog (MCXC). Our catalog provides the fluxes, positions, detection significances, HRs, one-year bin lightcurves, variability indices, and counterpart candidates.

We report on the detection and follow-up multi-wavelength observations of the new X-ray transient MAXI J1807+132 with the MAXI/GSC, Swift, and ground-based optical telescopes. The source was first recognized with the MAXI/GSC on 2017 March 13. About a week later, it reached the maximum intensity ($\sim$10 mCrab in 2-10 keV), and then gradually faded in $\sim$10 days by more than one order of magnitude. Time-averaged Swift/XRT spectra in the decaying phase can be described by a blackbody with a relatively low temperature (0.1-0.5 keV), plus a hard power-law component with a photon index of $\sim$2. These spectral properties are similar to those of neutron star low-mass X-ray binaries (LMXBs) in their dim periods. The blackbody temperature and the radius of the emission region varied in a complex manner as the source became dimmer. The source was detected in the optical wavelength on March 27-31 as well. The optical flux decreased monotonically as the X-ray flux decayed. The correlation between the X-ray and optical fluxes is found to be consistent with those of known neutron star LMXBs, supporting the idea that the source is likely to be a transient neutron star LMXB.

Thermonuclear X-ray bursts from accreting neutron stars power brief but strong irradiation of their surroundings, providing a unique way to study accretion physics. We analyze MAXI/GSC and Swift/XRT spectra of a day-long flash observed from IGR J17062-6143 in 2015. It is a rare case of recurring bursts at a low accretion luminosity of 0.15% Eddington. Spectra from MAXI, Chandra, and NuSTAR observations taken between the 2015 burst and the previous one in 2012 are used to determine the accretion column. We find it to be consistent with the burst ignition column of 5x10^10 g cm^-2, which indicates that it is likely powered by burning in a deep helium layer. The burst flux is observed for over a day, and decays as a straight power law: F~t^-1.15. The burst and persistent spectra are well described by thermal emission from the neutron star, Comptonization of this emission in a hot optically thin medium surrounding the star, and reflection off the photoionized accretion disk. At the burst peak, the Comptonized component disappears, when the burst may dissipate the Comptonizing gas, and it returns in the burst tail. The reflection signal suggests that the inner disk is truncated at ~100 gravitational radii before the burst, but may move closer to the star during the burst. At the end of the burst, the flux drops below the burst cooling trend for 2 days, before returning to the pre-burst level.

The performance of the Time Projection Chamber (TPC) polarimeter for the Polarimeter for Relativistic Astrophysical X-ray Sources (PRAXyS) Small Explorer was evaluated using polarized and unpolarized X-ray sources. The PRAXyS mission will enable exploration of the universe through X-ray polarimetry in the 2-10 keV energy band. We carried out performance tests of the polarimeter at the Brookhaven National Laboratory, National Synchrotron Light Source (BNL-NSLS) and at NASA's Goddard Space Flight Center. The polarimeter was tested with linearly polarized, monochromatic X-rays at 11 different energies between 2.5 and 8.0 keV. At maximum sensitivity, the measured modulation factors at 2.7, 4.5 and 8.0 keV are 27%, 43% and 59%, respectively and the measured angle of polarization is consistent with the expected value at all energies. Measurements with a broadband, unpolarized X-ray source placed a limit of less than 1% on false polarization in the PRAXyS polarimeter.

We report nine long X-ray bursts from neutron stars, detected with Monitor of All-sky X-ray Image (MAXI). Some of these bursts lasted for hours, and hence are qualified as superbursts, which are prolonged thermonuclear flashes on neutron stars and are relatively rare events. MAXI observes roughly 85% of the whole sky every 92 minutes in the 2-20 keV energy band, and has detected nine bursts with a long e-folding decay time, ranging from 0.27 to 5.2 hours, since its launch in 2009 August until 2015 August. The majority of the nine events were found to originate from transient X-ray sources. The persistent luminosities of the sources, when these prolonged bursts were observed, were lower than 1% of the Eddington luminosity for five of them and lower than 20% for the rest. This trend is contrastive to the 18 superbursts observed before MAXI, all but two of which originated from bright persistent sources. The distribution of the total emitted energy, i.e., the product of e-folding time and luminosity, of these bursts clusters around $10^{41}$-$10^{42}$ erg, whereas either of the e-folding time and luminosity ranges for an order of magnitude. Among the nine events, two were from 4U 1850-086 during the phases of relatively low persistent-flux, whereas it usually exhibits standard short X-ray bursts during outbursts.

We report the results of Suzaku observations of the young supernova remnant, Vela Jr.\ (RX J0852.0$-$4622), which is known to emit synchrotron X-rays, as well as TeV gamma-rays. Utilizing 39 Suzaku mapping observation data from Vela Jr., a significant hard X-ray emission is detected with the hard X-ray detector (HXD) from the north-west TeV-emitting region. The X-ray spectrum is well reproduced by a single power-law model with the photon index of 3.15$^{+1.18}_{-1.14}$ in the 12--22 keV band. Compiling this with the soft X-ray spectrum simultaneously observed with the X-ray imaging spectrometer (XIS) onboard Suzaku, we find that the wide-band X-ray spectrum in the 2--22 keV band is reproduced with a single power-law or concave broken power-law model, which are statistically consistent with each other. Whichever the model of a single or broken power-law is appropriate, clearly the spectrum has no rolloff structure. Applying this result to the method introduced in \citet{yama2014}, we find that one-zone synchrotron model with electron spectrum having a power-law plus exponential cutoff may not be applicable to Vela Jr.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

The 2.1-s anomalous X-ray pulsar 1E 1547.0-5408 exhibited an X-ray outburst on 2009 January 22, emitting a large number of short bursts. The wide-band all- sky monitor (WAM) on-board Suzaku detected at least 254 bursts in the 160keV-6.2MeV band over the period of January 22 00:57-17:02 UT from the direction of 1E 1547.0-5408. One of these bursts, which occurred at 06:45:13, produced the brightest fluence in the 0.5-6.2MeV range, with an averaged 0.16-6.2MeV flux and extrapolated 25 keV-2 MeV fluence of about 3x10-6 erg cm-2 s-1 and about 3x10-4 erg cm-2, respectively. After pile-up corrections, the time-resolved WAM spectra of this burst were well-fitted in the 0.16-6.2MeV range by two-component models; specifically, a blackbody plus an optically thin thermal bremsstrahlung or a combination of a blackbody and a power-law component with an exponential cutoff. These results are compared with previous works reporting the persistent emission and weaker short bursts followed by the same outburst.

The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August 2009. It is the fourth in a series of normal outbursts associated with the periastron, but is unusual by presenting a double-peaked light curve. The two peaks reached a flux of ~450 mCrab in the 15-50 keV range. We present results of the timing and spectral analysis of INTEGRAL, RXTE, and Suzaku observations of the outburst. The energy dependent pulse profiles and their evolution during the outburst are studied. No significant differences with respect to other normal outbursts are observed. The centroid energy of the fundamental cyclotron line shows no significant variation during the outburst. A spectral hardening with increasing luminosity is observed. We conclude that the source is accreting in the sub-critical regime. We discuss possible explanations for the double-peaked outburst.

We present analysis of 4U 1626-67, a 7.7 s pulsar in a low-mass X-ray binary system, observed with the hard X-ray detector of the Japanese X-ray satellite Suzaku in March 2006 for a net exposure of \sim88 ks. The source was detected at an average 10-60 keV flux of \sim4 x10^-10 erg cm^-2 s^-1. The phase-averaged spectrum is reproduced well by combining a negative and positive power-law times exponential cutoff (NPEX) model modified at \sim 37 keV by a cyclotron resonance scattering feature (CRSF). The phase-resolved analysis shows that the spectra at the bright phases are well fit by the NPEX with CRSF model. On the other hand, the spectrum in the dim phase lacks the NPEX high-energy cutoff component, and the CRSF can be reproduced by either an emission or an absorption profile. When fitting the dim phase spectrum with the NPEX plus Gaussian model, we find that the feature is better described in terms of an emission rather than an absorption profile. The statistical significance of this result, evaluated by means of an F-test, is between 2.91 x 10^-3 and 1.53 x 10^-5, taking into account the systematic errors in the background evaluation of HXD-PIN. We find that, the emission profile is more feasible than the absorption one for comparing the physical parameters in other phases. Therefore, we have possibly detected an emission line at the cyclotron resonance energy in the dim phase.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

The Be/X-ray binary A0535+26 showed a giant outburst in December 2009 that reached ~5.14 Crab in the 15-50 keV range. Unfortunately, due to Sun constraints it could not be observed by most X-ray satellites. The outburst was preceded by four weaker outbursts associated with the periastron passage of the neutron star. The fourth of them, in August 2009, presented a peculiar double-peaked light curve, with a first peak lasting about 9 days that reached a (15-50 keV) flux of 440 mCrab. The flux then decreased to less than 220 mCrab, and increased again reaching 440 Crab around the periastron. The outburst was monitored with INTEGRAL, RXTE, and Suzaku TOO observations. One orbital period (~111 days) after the 2009 giant outburst, a new and unexpectedly bright outburst took place (~1.4Crab in the 15-50 keV range). It was monitored with TOO obs ervations with INTEGRAL, RXTE, Suzaku, and Swift. First results of the spectral and timing analysis of these observations are presented, with a specific focus on the cyclotron lines present in the system and its variation with the mass accretion rate.

The WIDeField telescope for Gamma-ray burst Early Timing (WIDGET) is used for a fully automated, ultra-wide-field survey aimed at detecting the prompt optical emission associated with Gamma-ray Bursts (GRBs). WIDGET surveys the HETE-2 and Swift/BAT pointing directions covering a total field of view of 62 degree x 62 degree every 10 secounds using an unfiltered system. This monitoring survey allows exploration of the optical emission before the gamma-ray trigger. The unfiltered magnitude is well converted to the SDSS r' system at a 0.1 mag level. Since 2004, WIDGET has made a total of ten simultaneous and one pre-trigger GRB observations. The efficiency of synchronized observation with HETE-2 is four times better than that of Swift. There has been no bright optical emission similar to that from GRB 080319B. The statistical analysis implies that GRB080319B is a rare event. This paper summarizes the design and operation of the WIDGET system and the simultaneous GRB observations obtained with this instrument.

A normal outburst of the Be/X-ray binary system A0535+26 has taken place in August 2009. It is the fourth in a series of normal outbursts that have occured around the periastron passage of the source, but is unusual by starting at an earlier orbital phase and by presenting a peculiar double-peaked light curve. A first "flare" (lasting about 9 days from MJD 55043 on) reached a flux of 440 mCrab. The flux then decreased to less than 220 mCrab, and increased again reaching 440 mCrab around the periastron at MJD 55057. Target of Opportunity observations have been performed with INTEGRAL, RXTE and Suzaku. First results of these observations are presented, with special emphasis on the cyclotron lines present in the X-ray spectrum of the source, as well as in the pulse period and energy dependent pulse profiles of the source.