In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We present a systematic analysis of the phase lags associated with the type-C QPOs in GRS 1915+105 using RXTE data. Our sample comprises of 620 RXTE observations with type-C QPOs ranging from ~0.4 Hz to ~6.3 Hz. Based on our analysis, we confirm that the QPO phase lags decrease with QPO frequency, and change sign from positive to negative at a QPO frequency of ~2 Hz. In addition, we find that the slope of this relation is significantly different between QPOs below and above 2 Hz. The relation between the QPO lags and QPO rms can be well fitted with a broken line: as the QPO lags go from negative to positive, the QPO rms first increases, reaching its maximum at around zero lag, and then decreases. The phase-lag behaviour of the subharmonic of the QPO is similar to that of the QPO fundamental, where the subharmonic lags decrease with subharmonic frequency and change sign from positive to negative at a subharmonic frequency of ~1 Hz; on the contrary, the second harmonic of the QPO shows a quite different phase-lag behaviour, where all the second harmonics show hard lags that remain more or less constant. For both the QPO and its (sub)harmonics, the slope of the lag-energy spectra shows a similar evolution with frequency as the average phase lags. This suggests that the lag-energy spectra drives the average phase lags. We discuss the possibility for the change in lag sign, and the physical origin of the QPO lags.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows.

We present three simultaneous/quasi-simultaneous NuSTAR and Swift datasets of the black hole GRS 1716-249 in its hard intermediate state. The accretion disk in this state may have reached the innermost stable circular orbit, and the NuSTAR spectra show a broad relativistic iron line and a strong Compton hump. To measure the black hole spin, we construct a joint model consisting of a relativistic disk model kerrbb and a reflection model relxill, to fit the continuum and the reflection components, respectively. By applying this model to each dataset independently, a consistent result is obtained on the black hole spin and the disk inclination. The black hole spin is a* >~ 0.92, and the inclination angle (i) is around 40-50 degree, based on the measurements of all datasets. In the third dataset, a high black hole mass (M_BH) is strongly disfavored by the spectral fits. By unfreezing the black hole mass, we find a*>0.92, i=49.9^{+1.0}_{-1.3} degree and M_BH<8.0 Msun, at a 90% confidence level. Considering the lower limit derived from a previous optical constraint, M_BH is in a range of 4.9-8.0 Msun.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning and gamma ray burst (GRB) observations and , based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

We studied the long-term evolution of the spectral-temporal correlated properties of the black hole candidate Swift J1753.5-0127 from the onset of its outburst until 2011 with the Rossi X-ray Timing Explorer (RXTE). The source stayed most of its lifetime during hard state, with occasionally transitioned to the hard intermediate state. Similar to typical black hole transients, Swift J1753.5-0127 traces a clear hard line in absolute rms - intensity diagram during the low hard state, with expected highest absolute rms, while shows a clear turn during the hard intermediate state, accompanied by lower absolute rms. Different from Cyg X-1, we found that frequency-dependent time lag increased significantly in the 0.02 -- 3.2 Hz band during state transition in this source. The X-ray time lags in 0.02 -- 3.2 Hz can therefore be used as indicators of state transition in this source. Type-C QPO frequency is positively related with its fractional rms and X-ray photon index, suggesting a moving inwards disk/corona scenario. We discussed the physical interpretation of our results in our paper.

We report on the spectral behavior of the first Galactic ultraluminous X-ray pulsar Swift J0243.6+6124 with NuSTAR observations during its 2017-2018 outburst. At sub-Eddington levels, the source spectrum is characterized by three emission components, respectively from the accretion column, the hot spot, and a broad iron line emission region. When the source is above the Eddington limit, the hot spot temperature increases and the spectrum features two more blackbody components. One blackbody component has a radius of 10-20 km and is likely originated from the top of the accretion column. The other one saturates at a blackbody luminosity of (1 - 2)*10^38 erg/s, coincident with the Eddington limit of a neutron star. This is well consistent with the scenario that super-Eddington accretion onto compact objects will power optically-thick outflows and indicates an accretion rate 60-80 times the critical value. This suggests that super-Eddington accretion onto magnetized systems can also power massive winds. At super-Eddington levels, the iron line becomes more significant and blueshifted, and is argued to be associated with the ultrafast wind in the central funnel or jets. This source, if located in external galaxies, will appear like other ultraluminous pulsars.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

The black hole candidate MAXI J1535-571 is a recently discovered X-ray transient. We report on the monitoring observations of Swift Gamma-Ray Burst Mission during outburst. The source transits from the hard state to the intermediate state, and reaches the soft state; near the end of the outburst, it returns to the hard state passing through the low intermediate state, following a typical Q-shaped loop in the hardness-intensity diagram. During the high intermediate state, detailed spectral analyses using a multi-temperature disk model reveal that the disk luminosity is flatter than the inner disk temperature to the fourth power, and the disk temperature profile varies as the disk radius raised to the power of -0.5, consistent with the behavior of a slim disk, implying that the disk structure has been modified due to the high luminosity of ~ 10^{39} erg/s. Meanwhile, the column density increases with the Eddington ratio, suggesting that the outflow matters are driven by radiation pressure due to the high accretion rate.

We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope (HST) observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope (VLT) and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4-28) x 10^10 M_sun, and the mass of the companion galaxy of is estimated to be M_* = (6.8 +/- 2.4) x 10^9 M_sun, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2-5) x 10^8 M_sun and an Eddington ratio of about 0.01-0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142.

We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions.

We report on XMM-Newton/Chandra/Swift/HST observations of the ultraluminous X-ray source (ULX) in NGC 247, which is found to make transitions between the supersoft ultraluminous (SSUL) regime with a spectrum dominated by a cool ($\sim 0.1$ keV) blackbody component and the soft ultraluminous (SUL) regime with comparable luminosities shared by the blackbody and power-law components. Multi-epoch observations revealed an anti-correlation between the blackbody radius and temperature, $R_{\rm bb} \propto T_{\rm bb}^{-2.8 \pm 0.3}$, ruling out a standard accretion disk as the origin of the soft X-ray emission. The soft X-ray emission is much more variable on both short and long timescales in the SSUL regime than in the SUL regime. We suggest that the SSUL regime may be an extension of the ultraluminous state toward the high accretion end, being an extreme case of the SUL regime, with the blackbody emission arising from the photosphere of thick outflows and the hard X-rays being emission leaked from the embedded accretion disk via the central low-density funnel or advected through the wind. However, the scenario that the supersoft ULXs are standard ULXs viewed nearly edge-on cannot be ruled out. Flux dips on a timescale of 200 s were observed. The dips cannot be explained by an increase of absorption, but could be due to the change of accretion rate or related to thermal fluctuations in the wind or disk. The optical emission of NGC 247 ULX exhibits a blackbody spectrum at a temperature of 19,000 K with a radius of 20 $R_\odot$, likely arising from an OB supergiant companion star.

We report the detection of a $78.1\pm0.5$ day period in the X-ray lightcurve of the extreme ultraluminous X-ray source NGC 5907 ULX1 ($L_{\rm{X,peak}}\sim5\times10^{40}$ erg s$^{-1}$), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of $\sim$3-4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3$\sigma$ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

We present a multi-epoch hard X-ray analysis of Cygnus X-1 in its soft state based on four observations with NuSTAR. Despite the basic similarity of the observed spectra, there is clear spectral variability between epochs. To investigate this variability, we construct a model incorporating both the standard disk-corona continuum and relativistic reflection from the accretion disk, based on prior NuSTAR work on Cygnus X-1, and apply this model to each epoch independently. We find excellent consistency for the black hole spin, and the iron abundance of the accretion disk, which are expected to remain constant on observational timescales. In particular, we confirm that Cygnus X-1 hosts a rapidly rotating black hole, 0.93<a*<0.96, in broad agreement with the majority of prior studies of the relativistic disk reflection and constraints on the spin obtained through studies of the thermal accretion disk continuum. Our work also confirms the apparent misalignment between the inner disk and the orbital plane of the binary system reported previously, finding the magnitude of this warp to be ~10-15deg. This level of misalignment does not significantly change (and may even improve) the agreement between our reflection results and the thermal continuum results regarding the black hole spin. The spectral variability observed by NuSTAR is dominated by the primary continuum, implying variability in the temperature of the scattering electron plasma. Finally, we consistently observe absorption from ionized iron at ~6.7 keV, which varies in strength as a function of orbital phase in a manner consistent with the absorbing material being an ionized phase of the focused stellar wind from the supergiant companion star.

The black hole candidate XTE J1908+094 went into outburst for the first time since 2003 in October 2013. We report on an observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) and monitoring observations with Swift during the outburst. NuSTAR caught the source in the soft state: the spectra show a broad relativistic iron line, and the light curves reveal a ~40 ks flare with the count rate peaking about 40% above the non-flare level and with significant spectral variation. A model combining a multi-temperature thermal component, a power-law, and a reflection component with an iron line provides a good description of the NuSTAR spectrum. Although relativistic broadening of the iron line is observed, it is not possible to constrain the black hole spin with these data. The variability of the power-law component, which can also be modeled as a Comptonization component, is responsible for the flux and spectral change during the flare, suggesting that changes in the corona (or possibly continued jet activity) are the likely cause of the flare.

NGC 247 is a nearby late-type bulgeless spiral galaxy that contains an inactive nucleus. We report a serendipitous discovery of an X-ray flare from the galaxy center with a luminosity up to 2*10^39 erg/s in the 0.3-10 keV band with XMM-Newton. A Chandra observation confirms that the new X-ray source is spatially coincident with the galaxy nucleus. The XMM-Newton data revealed a hard power-law spectrum with a spectral break near 3-4 keV, no pulsations on timescales longer than 150 ms, and a flat power spectrum consistent with Poisson noise from 1 mHz to nearly 10 Hz. Follow-up observations with Swift detected a second flux peak followed by a luminosity drop by factor of almost 20. The spectral and temporal behaviors of the nuclear source are well consistent with the scenario that the flare was due to an outburst of a low-mass X-ray binary that contains a stellar-mass black hole emitting near its Eddington limit at the peak. However, it cannot be ruled out that the sudden brightening in the nucleus was due to accretion onto a possible low-mass nuclear black hole, fed by a tidally disrupted star or a gas cloud; the MAXI observations limit the peak luminosity of the flare to less than ~10^43 erg/s, suggesting that it is either a low mass black hole or an inefficient tidal disruption event (TDE).

The positron fraction in cosmic rays was found to be a steadily increasing in function of energy, above $\sim$ 10 GeV. This behaviour contradicts standard astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic rays during the propagation in the interstellar medium. The observed anomaly in the positron fraction triggered a lot of excitement, as it could be interpreted as an indirect signature of the presence of dark matter species in the Galaxy. Alternatively, it could be produced by nearby astrophysical sources, such as pulsars. Both hypotheses are probed in this work in light of the latest AMS-02 positron fraction measurements. The transport of the primary and secondary positrons in the Galaxy is described using a semi-analytic two-zone model. MicrOMEGAs is used to model the positron flux generated by dark matter species. The description of the positron fraction from astrophysical sources is based on the pulsar observations included in the ATNF catalogue. We find that the mass of the favoured dark matter candidates is always larger than 500 GeV. The only dark matter species that fulfils the numerous gamma ray and cosmic microwave background bounds is a particle annihilating into four leptons through a light scalar or vector mediator, with a mixture of tau (75%) and electron (25%) channels, and a mass between 0.5 and 1 TeV. The positron anomaly can also be explained by a single astrophysical source and a list of five pulsars from the ATNF catalogue is given. Those results are obtained with the cosmic ray transport parameters that best fit the B/C ratio. Uncertainties in the propagation parameters turn out to be very significant. In the WIMP annihilation cross section to mass plane for instance, they overshadow the error contours derived from the positron data.

We report on multiwavelength observations of the supersoft ultraluminous X-ray source (ULX) in NGC 247 made with the Chandra X-ray Observatory and Hubble Space Telescope (HST). We aligned the X-ray and optical images using three objects present on both and identified a unique, point-like optical counterpart to the ULX. The X-ray to optical spectrum is well fitted with an irradiated disk model if the extinction measured for Cepheids in NGC 247 is used. Assuming only Galactic extinction, then the spectrum can be modeled as a standard thin accretion disk. Either result leads to the conclusion that a disk interpretation of the X-ray spectrum is valid, thus the source may be in the X-ray thermal state and contain an intermediate mass black hole of at least 600 solar masses. In contrast to other supersoft ULXs which are transient and exhibit a luminosity temperature relation inconsistent with a disk interpretation of the X-ray emission, the NGC 247 ULX has a relatively steady flux and all available X-ray data are consistent with emission from a disk in the thermal state.

We report on UV and X-ray spectroscopy and broad-band optical observations of the ultraluminous X-ray source in Holmberg II. Fitting various stellar spectral models to the combined, non-simultaneous data set, we find that normal metallicity stellar spectra are ruled out by the data, while low metallicity, Z = 0.1 Z_{\odot}, late O-star spectra provide marginally acceptable fits, if we allow for the fact that X-ray ionization from the compact object may reduce or eliminate UV absorption/emission lines from the stellar wind. By contrast, an irradiated disk model fits both UV and optical data with chi^2/dof=175.9/178, and matches the nebular extinction with a reddening of E(B-V)=0.05^{+0.05}_{-0.04}. These results suggest that the UV/optical flux of Holmberg II X-1 may be dominated by X-ray irradiated disk emission.

We obtained 3 epochs of simultaneous Hubble Space Telescope (HST)/Wide Field Camera 3 and Chandra observations of the ultraluminous X-ray source NGC 5408 X-1. The counterpart of the X-ray source is seen in all HST filters, from the UV through the NIR, and for the first time, we resolve the optical nebula around the ULX. We identified a small OB association near the ULX that may be the birthplace of the system. The stellar association is young, ~ 5 Myr, contains massive stars up to 40 solar masses, and is thus similar to associations seen near other ULXs, albeit younger. The UV/Optical/NIR spectral energy distribution (SED) of the ULX counterpart is consistent with that of a B0I supergiant star. We are also able to fit the whole SED from the X-rays to the NIR with an irradiated disk model. The three epochs of data show only marginal variability and thus, we cannot firmly conclude on the nature of the optical emission.

Using archival Hubble Space Telescope (HST) imaging data, we report the multiband photometric properties of 13 ultraluminous X-ray sources (ULXs) that have a unique compact optical counterpart. Both magnitude and color variation are detected at time scales of days to years. The optical color, variability, and X-ray to optical flux ratio indicate that the optical emission of most ULXs is dominated by X-ray reprocessing on the disk, similar to that of low mass X-ray binaries. For most sources, the optical spectrum is a power-law, $F_{\nu} \propto \nu^{\alpha}$ with $\alpha$ in the range 1.0 to 2.0 and the optically emitting region has a size on the order of 1e12 cm. Exceptions are NGC 2403 X-1 and M83 IXO 82, which show optical spectra consistent with direct emission from a standard thin disk, M101 ULX-1 and M81 ULS1, which have X-ray to optical flux ratios more similar to high-mass X-ray binaries, and IC 342 X-1, in which the optical light may be dominated by the companion star. Inconsistent extinction between the optical counterpart of NGC 5204 X-1 and the nearby optical nebulae suggests that they may be unrelated.

We present Hubble Space Telescope observations in the far UV of the ultraluminous X-ray source in NGC 6946 associated with the optical nebula MF 16. Both a point-like source coincident with the X-ray source and the surrounding nebula are detected in the FUV. The point source has a flux of 5E-16 erg s^-1 cm^-2 Ang^-1 and the nebula has a flux of 1.6E-15 erg s^-1 cm^-2 Ang^-1, quoted at 1533 Ang and assuming an extinction of A_V = 1.54. Thus, MF 16 appears to host the first directly detected ultraluminous UV source (ULUV). The flux of the point-like source is consistent with a blackbody with T ~ 30,000 K, possibly from a massive companion star, but this spectrum does not create sufficient ionizing radiation to produce the nebular HeII flux and a second, hotter emission component would be required. A multicolor disk blackbody spectrum truncated with an outer disk temperature of ~16,000 K provides an adequate fit to the FUV, B, V, I, and HeII fluxes and can produce the needed ionizing radiation. Additional observations are required to determine the physical nature of the source.

One of a family previously proposed ``central engines'' for cosmic gamma-ray burst sources (Klu\'zniak & Ruderman 1998) is considered in some detail. A steadily accreting $10^6$ Gauss magnetic white dwarf should ultimately collapse to a strongly differentially rotating, millisecond-rotation-period neutron star for a wide range of steady accretion rates and initial masses if the accreting white dwarf has an evolved O-Ne-Mg composition. A similar neutron star could also result from an initial C-O white dwarf but only for more constrained accretion rates. Because the collapsing white dwarf begins as a $\gamma=4/3$ polytrope, the final neutron star's spin-rate increases strongly with cylindrical radius. A stable wind-up of the neutron star's poloidal magnetic field then produces buoyant magnetic toroids which grow, break loose, rise, and partly penetrate the neutron star surface to form a transient, $B\approx 10^{17}$ G millisecond spin-period pulsar with a powerful pulsar wind (Usov 1992). This pulsar wind emission is then rapidly suppressed by the surface shear motion from the strong stellar differential rotation. This wind-up and transient pulsar formation can occur at other times on different cylinders and/or repeat on the same one, with (re-)wind up and surface penetration time scales hugely longer han the neutron star's millisecond spin period. In this way, differential rotation both opens and closes the doors which allow neutron star spin-energy to be emitted in powerful bursts of pulsar wind. Predictions of this model compare favorably to needed central engine properties of gamma-ray burst sources (total energy, duration, sub-burst fluctuations and time scales, variability among burst events, and baryon loading).

Beginning from a relatively simple set of dynamical equations for a fluid permeated by a radiative field strong enough to produce significant forces, we find the structure of plane-parallel equilibria and study their stability to small acoustic disturbances. In doing this, we neglect viscous effects and complications of nongreyness. We find that acoutic instabilities occur over a wide range of conditions below the Eddington limit. This result is in line with findings reported twenty years ago but it contradicts some more recent reports of the absence of instabilities. We briefly attempt to identify the causes of the discrepancies and then close with a discussion of the possible astrophysical interest of such instabilities.

We study the weakly nonlinear evolution of acoustic instability of a plane- parallel polytrope with thermal dissipation in the form of Newton's law of cooling. The most unstable horizontal wavenumbers form a band around zero and this permits the development of a nonlinear pattern theory leading to a complex Ginzburg-Landau equation (CGLE). Numerical solutions for a subcritical, quintic CGLE produce vertically oscillating, localized structures that resemble the oscillons observed in recent experiments of vibrated granular material.