Source author record

A. K. H. Kong

A. K. H. Kong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

preprint2022arXiv

Multi-wavelength properties of 850-$μ$m selected sources from the North Ecliptic Pole SCUBA-2 survey

We present the multi-wavelength counterparts of 850-$μ$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$μ$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$μ$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

Searching for cataclysmic variable stars in unidentified X-ray sources

We carry out a photometric search for new cataclysmic variable stars (CVs), with the goal of identification for candidates of AR~Scorpii-type binary systems. We select GAIA sources that are likely associated with unidentified X-ray sources, and analyze the light curves taken by the Zwicky Transient Facility, Transiting Exoplanet Survey Satellite, and Lulin One-meter Telescope in Taiwan. We investigate eight sources as candidates for CVs, among which six sources are new identifications. Another two sources have been recognized as CVs in previous studies, but no detailed investigations have been done. We identify two eclipsing systems that are associated with an unidentified XMM-Newton or Swift source, and one promising candidate for polar associated with an unidentified ASKA source. Two polar candidates may locate in the so-called period gap of a CV, and the other six candidates have an orbital period shorter than that of the period gap. Although we do not identify a promising candidate for AR~Scorpii-type binary systems, our study suggests that CV systems that have X-ray emission and do not show frequent outbursts may have been missed in previous surveys.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Investigation of the timing and spectral properties of an Ultra-luminous X-ray pulsar NGC 7793 P13

We perform both timing and spectral analyses using the archival X-ray data taken with Swift, XMM-Newton, NICER, and NuSTAR from 2016 to 2020 to study an ultra-luminous pulsar, NGC 7793 P13, that showed a long period of super-Eddington accretion. We use the Rayleigh test to investigate the pulsation at different epochs, and confirm the variation of the pulse profile with the finite Gaussian mixture modelling and two-sample Kuiper test. Taking into account the periodic variation of the spin periods caused by the orbital Doppler effect, we further determine an orbital period of ~65 d and show that no significant correlation can be detected between the orbital phase and the pulsed fraction. The pulsed spectrum of NGC7793 P13 in 0.5-20keV can be simply described using a power-law with a high energy exponential cutoff, while the broad-band phase-averaged spectrum of the same energy range requires two additional components to account for the contribution of a thermal accretion disk and the Comptonization photons scattered into the hard X-rays. We find that NGC 7793 P13 stayed in the hard ultra-luminous state and the pulsed spectrum was relatively soft when the source was faint in the end of 2019. Moreover, an absorption feature close to 1.3 keV is marginally detected from the pulsed spectra and it is possibly associated with a cyclotron resonant scattering feature.

preprint2021arXiv

Overview of KAGRA: Calibration, detector characterization, physical environmental monitors, and the geophysics interferometer

KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu, Japan. In this series of articles, we present an overview of the baseline KAGRA, for which we finished installing the designed configuration in 2019. This article describes the method of calibration (CAL) used for reconstructing gravitational wave signals from the detector outputs, as well as the characterization of the detector (DET). We also review the physical environmental monitors (PEM) system and the geophysics interferometer (GIF). Both are used for characterizing and evaluating the data quality of the gravitational wave channel. They play important roles in utilizing the detector output for gravitational wave searches. These characterization investigations will be even more important in the near future, once gravitational wave detection has been achieved, and in using KAGRA in the gravitational wave astronomy era.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2020arXiv

A multi-wavelength study of the gamma-ray binary candidate HESS J1832-093

We investigate the nature of the unidentified very-high-energy (VHE) gamma-ray object, HESS J1832-093, in a multi-wavelength context. Based on X-ray variability and spectral index ($Γ_X\sim\,1.5$), and its broad-band spectrum (which was remarkably similar to HESS J0632+057, a confirmed gamma-ray binary), HESS J1832-093 has been considered to be a strong gamma-ray binary candidate in previous works. In this work, we provide further evidence for this scenario. We obtained a spectrum of its IR counterpart using Gemini/Flamingo, finding absorption lines that are usually seen in massive stars, in particular O stars. We also obtained a rather steep ATCA spectrum ($α=-1.18^{+1.04}_{-0.88}$) which prefers a gamma-ray binary over an AGN scenario. Based on spatial-spectral analysis and variability search, we found that 4FGL J1832.9-0913 is possible to be associated with SNR G22.7-0.2 rather than with HESS J1832-093 only.

preprint2020arXiv

A Variable X-ray Source Close to the Magnetar SGR 1935+2154

The recent discovery of a millisecond radio burst from the Galactic magnetar SGR 1935+2154 has initiated an intense discussion about the connection between magnetars and fast radio bursts (FRBs). Although some properties of the radio burst from SGR 1935+2154 are not the same as cosmological FRBs, there are theoretical models which propose a connection between magnetars and FRBs (see review by Margalit et al. 2020). In particular, the role of a magnetar wind nebula is included in some models, and therefore it is worthwhile to investigate the X-ray environment of SGR 1935+2154 in more detail. Here, we report on the discovery of an X-ray transient feature near SGR 1935+2154 using archival Chandra data and discuss its possible origin.

preprint2020arXiv

An X-ray study for white dwarf binary AR Scorpii

We report a study of the X-ray emission from the white dwarf/M-type star binary system AR Scorpii using archival data taken in 2016-2020. It has been known that the X-ray emission is dominated by the optically thin thermal plasma emission, and its flux level varies significantly over the orbital phase. The X-ray emission also contains a component that modulates with the beat frequency between the white dwarf's spin frequency and orbital frequency. In this new analysis, the 2020 data taken by NICER shows that the X-ray emission is modulating with the spin frequency as well as the beat frequency, indicating that part of the X-ray emission is coming from the white dwarf's magnetosphere. It is found that the signal of the spin frequency appears only at a specific orbital phase, while the beat signal appears over the orbital phase. We interpret the X-ray emission modulating with the spin frequency and the beat frequency as a result of the synchrotron emission from electrons with a smaller and larger pitch angle, respectively. In a long-term evolution, the beat pulse profile averaged over the orbital phase changed from a single-peak structure in 2016/2018 to a double-peak structure in 2020. The observed X-ray flux levels measured in 2016/2017 are higher than those measured in 2018/2020. The plasma temperature and amplitude of the orbital waveform might vary with time too. These results indicate that the X-ray emission from AR Scorpii evolves on a timescale of years. This long-term evolution would be explained by a super-orbital modulation related to, for example, a precession of the white dwarf, or a fluctuation of the system related to activity of the companion star.

preprint2020arXiv

Application of the independent component analysis to the iKAGRA data

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

preprint2020arXiv

Multi-Epoch X-ray Imaging of Globular Cluster M62 with Chandra

Using archival spectral-imaging data with a total exposure of $\sim144$ ks obtained by {\it Chandra}, 43 X-ray sources are detected within the half-light radius of globular cluster M62 (NGC6266). Based on the X-ray colour-luminosity diagram or the positional coincidences with known sources, we have classified these sources into different groups of compact binaries including cataclysmic variable (CV), quiescent low mass X-ray binary (qLMXB), millisecond pulsar (MSP) and black hole (BH). Candidates of the X-ray counterparts of 12 CVs, 4 qLMXBs, 2 MSPs and 1 BH are identified in our analysis. The data used in our analysis consist of two frames separated by 12 years, which enable us to search for the long-term variability as well as the short-term X-ray flux variability within each observation window. Evidence for the short-term variability and long-term variability have been found in 7 and 12 sources respectively. For a number of bright sources with X-ray luminosities $L_{x}\gtrsim 10^{32}$~erg/s, we have characterized their spectral properties in further details. By comparing the X-ray population in M62 with those in several other prototypical globular clusters, we found the proportion of bright sources is larger in M62 which can possibly be a result of their active dynamical formation processes.

preprint2020arXiv

Multi-wavelength observations of the BL Lac object Fermi J1544-0649: one year after its awakening

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($Γ_γ\sim$1.7) in the GeV band and at times also very hard (($Γ_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

preprint2020arXiv

Overview of KAGRA: Detector design and construction history

KAGRA is a newly built gravitational-wave telescope, a laser interferometer comprising arms with a length of 3\,km, located in Kamioka, Gifu, Japan. KAGRA was constructed under the ground and it is operated using cryogenic mirrors that help in reducing the seismic and thermal noise. Both technologies are expected to provide directions for the future of gravitational-wave telescopes. In 2019, KAGRA finished all installations with the designed configuration, which we call the baseline KAGRA. In this occasion, we present an overview of the baseline KAGRA from various viewpoints in a series of of articles. In this article, we introduce the design configurations of KAGRA with its historical background.

preprint2020arXiv

Periodicity search for Pulsar Binaries with TESS

Pulsar binaries, in particular redback systems, provide good sources to study the pulsar wind flow and its interaction with the companion stars. {\it Fermi}-LAT have proposed probable pulsar binary candidates in its catalogs. To identify pulsar binary sources from the catalog, orbital modulation search of binary candidates is an effective way. TESS observes in survey mode for a large part of the sky and thus provide an excellent data set to periodicity search of pulsar binary candidates by observing the flux variation, thought to mainly come from the stellar companion. Using TESS data we look for flux modulation of five pulsar binaries (or candidates) with reported orbital periods, including PSR J1023+0038, 3FGL J0523.3-2528, 3FGL J0212.1+5320, 3FGL J0744.1-2523 and PSR J1417-4402, demonstrating that TESS photometric data are very useful in identifying periodicities of redback-like systems. This method can be effective in searches for new pulsar binaries or similar binary systems in the future.

preprint2020arXiv

Repeated state change of variable gamma-ray pulsar, PSR J2021+4026

PSR J2021+4026 is a radio-quiet gamma-ray pulsar and the first pulsar that shows state change of the gamma-ray emission and spin-down rate. The state change of PSR J2021+4026 was first observed at 2011 October, at which the pulsar changes the state from high gamma-ray flux/low spin-down rate state to low gamma-ray flux/high spin-down rate st\ ate. In December 2014, PSR J2021+4026 recovered the state before the 2011 state change over a timescale of a few months. We report that the long term evolution of the gamma-ray flux and timing behavior suggests that PSR J2021+4026 changed the state near 2018 February 1st and entered a new low gamma-ray flux/high spin-down rate state. At the 2018 state change, the averaged flux dropped from $(1.29\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1}}$ to $(1.12\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1\ }}$, which has the similar behavior to the case of 2011 event. The spin-down rate has increased by $\sim 3\%$ in the new state since the 2018 state change. The shapes of pulse profile and spectrum in GeV bands also changed at the 2018 event, and they are consistent with behavior at the 2011 state change. Our results probably suggest that PSR J2021+4026 is switching between different states with a timescale of several years, like some radio pulsars (e.g. PSR~B1828-11). PSR J2021+4026 will provide a unique opportunity to study the mechanism of the state switching.

preprint2020arXiv

Searches for Pulsar-like Candidates from Unidentified Objects in the Third Catalog of Hard Fermi-LAT (3FHL) sources with Machine Learning Techniques

We report the results of searching pulsar-like candidates from the unidentified objects in the $3^{\rm rd}$ Catalog of Hard Fermi-LAT sources (3FHL). Using a machine-learning based classification scheme with a nominal accuracy of $\sim98\%$, we have selected 27 pulsar-like objects from 200 unidentified 3FHL sources for an identification campaign. Using archival data, X-ray sources are found within the $γ-$ray error ellipses of 10 3FHL pulsar-like candidates. Within the error circles of the much better constrained X-ray positions, we have also searched for the optical/infrared counterparts and examined their spectral energy distributions. Among our short-listed candidates, the most secure identification is the association of 3FHL J1823.3-1339 and its X-ray counterpart with the globular cluster Mercer 5. The $γ-$rays from the source can be contributed by a population of millisecond pulsars residing in the cluster. This makes Mercer 5 as one of the slowly growing hard $γ-$ray population of globular clusters with emission $>10$ GeV. Very recently, another candidate picked by our classification scheme, 3FHL J1405.1-6118, has been identified as a new $γ-$ray binary with an orbital period of $13.7$ days. Our X-ray analysis with a short Chandra observation has found a possible periodic signal candidate of $\sim1.4$ hrs and a putative extended X-ray tail of $\sim20$ arcsec long. Spectral energy distribution of its optical/infrared counterpart conforms with a blackbody of $T_{\rm bb}\sim40000$ K and $R_{\rm bb}\sim12R_{\odot}$ at a distance of 7.7 kpc. This is consistent with its identification as an early O star as found by infrared spectroscopy.

preprint2020arXiv

The X-ray Emissivity of Low-Density Stellar Populations

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.

preprint2019arXiv

A timing study of MAXI J1820+070 based on Swift/XRT and NICER monitoring in 2018/19

We present a detailed timing analysis of the bright black hole X-ray binary MAXI J1820+070 (ASASSN-18ey), during its first detected outburst lasting from March 2018 until October 2019 based on Swift/XRT window timing mode observations, corresponding UVOT data and NICER observations. The light curves clearly show four outbursts, with the source remaining in the hard state during its first outburst, while the rise of the second outburst corresponds with the transition to the soft state. A similar double outburst of GX339-4 has been observed in 2004. Here it is followed by two hard-state only outbursts. In many observations the power density spectra showed type-C quasi-periodic oscillations with characteristic frequencies below 1 Hz, which suggests that the source stayed in a state of low effective accretion for large parts of its outburst. The absence of other types of quasi-periodic oscillations hinders a precise determination of the state transitions, but from combining NICER and Swift/XRT data, we find that MAXI J1820+070 went from the hard-intermediate to the soft state in less than one day. The covariance ratios derived from NICER data show an increase towards lower energies, which indicate that the source should make a transition to the soft state. This transition finally took place, after MAXI J1820+070 stayed in the hard state at rather constant luminosity for about 116 days. The steepness of the increase of the covariance ratios is not correlated with the amount of rms variability and it does not show a monotonic evolution along the outburst.

preprint2019arXiv

An arm length stabilization system for KAGRA and future gravitational-wave detectors

Modern ground-based gravitational wave (GW) detectors require a complex interferometer configuration with multiple coupled optical cavities. Since achieving the resonances of the arm cavities is the most challenging among the lock acquisition processes, the scheme called arm length stabilization (ALS) had been employed for lock acquisition of the arm cavities. We designed a new type of the ALS, which is compatible with the interferometers having long arms like the next generation GW detectors. The features of the new ALS are that the control configuration is simpler than those of previous ones and that it is not necessary to lay optical fibers for the ALS along the kilometer-long arms of the detector. Along with simulations of its noise performance, an experimental test of the new ALS was performed utilizing a single arm cavity of KAGRA. This paper presents the first results of the test where we demonstrated that lock acquisition of the arm cavity was achieved using the new ALS and residual noise was measured to be $8.2\,\mathrm{Hz}$ in units of frequency, which is smaller than the linewidth of the arm cavity and thus low enough to lock the full interferometer of KAGRA in a repeatable and reliable manner.

preprint2019arXiv

Peculiar outbursts of an ultra luminous source likely signs of an aperiodic disc-wind

The metal rich globular cluster RZ 2109 in the massive Virgo elliptical galaxy NGC 4472 (M49) harbours the ultra luminous X-ray source XMMU 122939.7+075333. Previous studies showed that this source varies between bright and faint phases on timescales of just a few hours. Here, we report the discovery of two peculiar X-ray bursting events that last for about 8 and 3.5 hours separated by about 3 days. Such a recurring X-ray burst-like behaviour has never been observed before. We argue that type-I X-ray bursts or super bursts as well as outburst scenarios requiring a young stellar object are highly unlikely explanations for the observed light curve, leaving an aperiodic disc wind scenario driven by hyper-Eddington accretion as a viable explanation for this new type of X-ray flaring activities.

preprint2016arXiv

A likely inverse-Compton emission from the Type IIb SN 2013df

The inverse-Compton X-ray emission model for supernovae has been well established to explain the X-ray properties of many supernovae for over 30 years. However, no observational case has yet been found to connect the X-rays with the optical lights as they should be. Here, we report the discovery of a hard X-ray source that is associated with a Type II-b supernova. Simultaneous emission enhancements have been found in both the X-ray and optical light curves twenty days after the supernova explosion. While the enhanced X-rays are likely dominated by inverse-Compton scatterings of the supernova's lights from the Type II-b secondary peak, we propose a scenario of a high-speed supernova ejecta colliding with a low-density pre-supernova stellar wind that produces an optically thin and high-temperature electron gas for the Comptonization. The inferred stellar wind mass-loss rate is consistent with that of the supernova progenitor candidate as a yellow supergiant detected by the Hubble Space Telescope, providing an independent proof for the progenitor. This is also new evidence of the inverse-Compton emission during the early phase of a supernova.

preprint2016arXiv

A Possible 55-day X-ray Period of the Ultraluminous Accreting Pulsar M82 X-2

We report a possible detection of a 55-day X-ray modulation for the ultraluminous accreting pulsar M82 X-2 from archival Chandra observations. Because M82 X-2 is known to have a 2.5-day orbital period, if the 55-day period is real, it will be the superorbital period of the system. We also investigated variabilities of other three nearby ultraluminous X-ray sources in the central region of M82 with the Chandra data and did not find any evidence of periodicities. Furthermore, we re-examined the previously reported 62-day periodicity near the central region of M82 by performing a systematic timing study with all the archival Rossi X-Ray Timing Explorer and Swift data. Using various dynamic timing analysis methods, we confirmed that the 62-day period is not stable, suggesting that it is not the orbital period of M82 X-1 in agreement with previous work.

preprint2016arXiv

Millihertz Quasi-periodic Oscillations in 4U 1636-536: Putting Possible Constraints on the Neutron Star Size

Based on previous studies of quasi-periodic oscillations in neutron star LMXBs, mHz quasi-periodic oscillations (QPO) are believed to be related to `marginally stable' burning on the neutron star (NS) surface. Our study of phase resolved energy spectra of these oscillations in 4U 1636-53 shows that the oscillations are not caused by variations in the blackbody temperature of the neutron star, but reveals a correlation between the change of the count rate during the mHz QPO pulse and the spatial extend of a region emitting blackbody emission. The maximum size of the emission area $R^2_{\mathrm{BB}}=216.7^{+93.2}_{-86.4}$km$^2$, provides the direct evidence that the oscillations originate from a variable surface area constrained on the NS and are therefore not related to instabilities in the accretion disk. The obtained lower limit on the size of the neutron star (11.0 km) rules out equations of state that prefer small NS radii. Observations of mHz QPOs in NS LMXBs with NICER and eXTP will reduce the statistical uncertainty in the lower limit on the NS radius, which together with better estimates of the hardening factor and distance, will allow improving discrimination between different equations of state and compact star models. Furthermore, future missions will allow us to measure the peak blackbody emission area for a single mHz QPO pulse, which will push the lower limit to larger radii.

preprint2016arXiv

Observing two dark accelerators around the Galactic Centre with Fermi Large Area Telescope

We report the results from a detailed $γ-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $γ-$ray spectrum can be modeled with a single power-law with a photon index of $Γ\sim2.5$ from few hundreds MeV to TeV. Moreover, an elongated feature, which extends from "Region A" toward northwest for $\sim1.3^{\circ}$, is discovered for the first time. The orientation of this feature is similar to that of a large scale atomic/molecular gas distribution. For HESS J1741-302, our analysis does not yield any MeV-GeV counterpart for this unidentified TeV source. On the other hand, we have detected a new point source, Fermi J1740.1-3013, serendipitously. Its spectrum is apparently curved which resembles that of a $γ-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023.

preprint2016arXiv

Recurring X-ray Outbursts in the Supernova Impostor SN~2010da in NGC~300

We present new observations of the "supernova impostor" SN 2010da using the Chandra X-ray Observatory and the Hubble Space Telescope. During the initial 2010 outburst, the 0.3-10 keV luminosity was observed by Swift to be $\sim5\times10^{38}$ erg s$^{-1}$ and faded by a factor of $\sim$25 in a four month period. Our two new Chandra observations show a factor of $\sim$10 increase in the 0.35-8 keV X-ray flux, from $\sim$4$\times10^{36}$ erg s$^{-1}$ to $4\times10^{37}$ erg s$^{-1}$ in $\sim$6 months, and the X-ray spectrum is consistent in both observations with a power law photon index of $Γ\sim0$. We find evidence of X-ray spectral state changes: when SN 2010da is in a high-luminosity state, the X-ray spectrum is harder ($Γ\sim0$) compared to the low-luminosity state ($Γ\sim1.2\pm0.8$). Using our Hubble observations, we fit the color magnitude diagram of the coeval stellar population to estimate a time since formation of the SN 2010da progenitor system of $\lesssim$5 Myr. Our observations are consistent with SN 2010da being a high-mass X-ray binary (HMXB) composed of a neutron star and a luminous blue variable-like companion, although we cannot rule out the possibility that SN 2010da is an unusually X-ray bright massive star. The $\lesssim$5 Myr age is consistent with the theoretically predicted delay time between the formation of a massive binary and the onset of the HMXB phase. It is possible that the initial 2010 outburst marked the beginning of X-ray production in the system, making SN 2010da possibly the first massive progenitor binary ever observed to evolve into an HMXB.

preprint2016arXiv

The 2015 hard-state only outburst of GS 1354-64

Since its outburst in 1997 GS 135-64 stayed in quiescence. In June 2015 renewed activity of GS 1354-64 was observed. Based on our analysis of energy spectra and timing properties obtained from Swift/XRT monitoring data we found that GS 1354-64 stayed in the hard state during the entire outburst. Such a hard state only (or "failed" outburst) has also been observed in 1997. In addition, we analysed an XMM-Newton observation taken on August 6th. We compared variability on long and short time scales using covariance ratio and found that the ratio showed a decrease towards lower energies instead of the increase that has been found in other black hole X-ray binaries. There are now two sources (H 1743-322 and GS 1354-64) that do not show an increase towards lower energies in their covariance ratio. Both sources have been observed during "failed" outbursts and showed photon indices much harder than what is usually observed in black hole X-ray binaries.

preprint2015arXiv

Exploring the intrabinary shock from the redback millisecond pulsar PSR J2129-0429

We have investigated the intrabinary shock emission from the redback millisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modulation in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak structure with a dip in between. The observed X-rays are non-thermal dominant which can be modeled by a power-law with a photon index of ~1.2. Intrabinary shock can be the origin of the observed X-rays. The UV light curve is resulted from the ellipsoidal modulation of the companion. Modeling the UV light curve prefers a large viewing angle. The heating effect of the UV light curve is found to be negligible which suggests the high energy radiation beam of PSR J2129-0429 does not direct toward its companion. On the other hand, no significant orbital modulation can be found in gamma-ray which suggests the majority of the gamma-rays come from the pulsar.

preprint2015arXiv

Long-term X-ray Variability of Ultraluminous X-ray Sources

Long-term X-ray modulations on timescales from tens to hundreds of days have been widely studied for X-ray binaries located in the Milky Way and the Magellanic Clouds. For other nearby galaxies, only the most luminous X-ray sources can be monitored with dedicated observations. We here present the first systematic study of long-term X-ray variability of four ultraluminous X-ray sources (ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, and NGC 5408 X-1) monitored with Swift. By using various dynamic techniques to analyse their light curves, we find several interesting low-frequency quasi-periodicities. Although the periodic signals may not represent any stable orbital modulations, these detections reveal that such long-term regular patterns may be related to superorbital periods and structure of the accretion discs. In particular, we show that the outburst recurrence time of ESO 243-49 HLX-1 varies over time and suggest that it may not be the orbital period. Instead, it may be due to some kinds of precession, and the true binary period is expected to be much shorter.

preprint2015arXiv

Searches for millisecond pulsar candidates among the unidentified Fermi objects

Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $γ-$ray sources in this catalog are identified as promising MSP candidates based on their $γ$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated their X-ray fluxes and compared with the corresponding $γ$-ray fluxes. The X-ray to $γ$-ray flux ratios for 2FGL J1653.6-0159 and 2FGL J1946.4-5402 are comparable with the typical value for pulsars. For 2FGL J1625.2-0020, 2FGL J1653.6-0159 and 2FGL J1946.4-5402, their candidate X-ray counterparts are bright enough for performing a detailed spectral and temporal analysis to discriminate their thermal/non thermal nature and search for the periodic signal. We have also searched for possible optical/IR counterparts at the X-ray positions. For the optical/IR source coincident with the brightest X-ray object that associated with 2FGL J1120.0-2204, its spectral energy distribution is comparable with a late-type star. Evidence for the variability has also been found by examining its optical light curve. All the aforementioned 2FGL sources resemble a pulsar in one or more aspects, which make them as the promising targets for follow-up investigations.

preprint2015arXiv

The Chandra Local Volume Survey I: The X-ray Point Source Populations of NGC 55, NGC 2403, and NGC 4214

We present comprehensive X-ray point source catalogs of NGC~55, NGC~2403, and NGC~4214 as part of the Chandra Local Volume Survey. The combined archival observations have effective exposure times of 56.5 ks, 190 ks, and 79 ks for NGC~55, NGC~2403, and NGC~4214, respectively. When combined with our published catalogs for NGC 300 and NGC 404, our survey contains 629 X-ray sources total down to a limiting unabsorbed luminosity of $\sim5\times10^{35}$ erg s$^{-1}$ in the 0.35-8 keV band in each of the five galaxies. We present X-ray hardness ratios, spectral analysis, radial source distributions, and an analysis of the temporal variability for the X-ray sources detected at high significance. To constrain the nature of each X-ray source, we carried out cross-correlations with multi-wavelength data sets. We searched overlapping Hubble Space Telescope observations for optical counterparts to our X-ray detections to provide preliminary classifications for each X-ray source as a likely X-ray binary, background AGN, supernova remnant, or foreground star.

preprint2014arXiv

A detailed X-ray investigation of PSR J2021+4026 and $γ$-Cygni supernova remnant

We have investigated the field around the radio-quiet $γ$-ray pulsar, PSR J2021+4026, with a ~140 ks XMM-Newton observation and a ~56 ks archival Chandra data. Through analyzing the pulsed spectrum, we show that the X-ray pulsation is purely thermal in nature which suggests the pulsation is originated from a hot polar cap with $T\sim3\times10^{6}$ K on the surface of a rotating neutron star. On the other hand, the power-law component that dominates the pulsar emission in the hard band is originated from off-pulse phases, which possibly comes from a pulsar wind nebula. In re-analyzing the Chandra data, we have confirmed the presence of bow-shock nebula which extends from the pulsar to west by ~10 arcsec. The orientation of this nebular feature suggests that the pulsar is probably moving eastward which is consistent with the speculated proper motion by extrapolating from the nominal geometrical center of the supernova remnant (SNR) G78.2+2.1 to the current pulsar position. For G78.2+2.1, our deep XMM-Newton observation also enables a study of the central region and part of the southeastern region with superior photon statistics. The column absorption derived for the SNR is comparable with that for PSR J2021+4026, which supports their association. The remnant emission in both examined regions are in an non-equilibrium ionization state. Also, the elapsed time of both regions after shock-heating is apparently shorter than the Sedov age of G78.2+2.1. This might suggest the reverse shock has reached the center not long ago. Apart from PSR J2021+4026 and G78.2+2.1, we have also serendipitously detected an X-ray flash-like event XMM J202154.7+402855 from this XMM-Newton observation.

preprint2014arXiv

Chandra detection of a new diffuse X-ray component from the globular cluster 47 Tucanae

In re-analyzing the archival Chandra data of the globular cluster 47 Tucanae, we have detected a new diffuse X-ray emission feature within the half-mass radius of the cluster. The spectrum of the diffuse emission can be described by a power-law model plus a plasma component with photon index $Γ\sim1.0$ and plasma temperature $kT\sim0.2$ keV. While the thermal component is apparently uniform, the non-thermal contribution falls off exponentially from the core. The observed properties could possibly be explained in the context of multiple shocks resulted from the collisions among the stellar wind in the cluster and the inverse Compton scattering between the pulsar wind and the relic photons.

preprint2014arXiv

Discovery of an ultracompact gamma-ray millisecond pulsar binary candidate

We report multi-wavelength observations of the unidentified Fermi object 2FGL J1653.6-0159. With the help of high-resolution X-ray observation, we have identified an X-ray and optical counterpart of 2FGL J1653.6-0159. The source exhibits a periodic modulation of 75 min in optical and possibly also in X-ray. We suggest that 2FGL J1653.6-0159 is a compact binary system with an orbital period of 75 min. Combining the gamma-ray and X-ray properties, 2FGL J1653.6-0159 is potentially a black widow/redback type gamma-ray millisecond pulsar (MSP). The optical and X-ray lightcurve profiles show that the companion is mildly heated by the high-energy emission and the X-rays are from intrabinary shock. Although no radio pulsation has been detected yet, we estimated that the spin period of the MSP is ~2ms based on a theoretical model. If pulsation can be confirmed in the future, 2FGL J1653.6-0159 will become the first ultracompact rotation-powered MSP.

preprint2014arXiv

Discovery of X-ray pulsations from "next Geminga" - PSR J1836+5925

We report the X-ray pulsation of ~173.3 ms for the "next Geminga", PSR J1836+5925, with recent XMM-Newton investigations. The X-ray periodicity is consistent wtih the gamma-ray ephemeris at the same epoch. The X-ray folded light curve has a sinusoidal structure which is different from the double-peaked gamma-ray pulse profile. We have also analysed the X-ray phase-averaged spectra which shows the X-ray emission from PSR J1836+5925 is thermal dominant. This suggests the X-ray pulsation mainly originates from the modulated hot spot on the stellar surface.

preprint2014arXiv

High-energy emissions from the gamma-ray binary LS 5039

We study mechanisms of multi-wavelength emissions (X-ray, GeV and TeV gamma-rays) from the gamma-ray binary LS~5039. This paper is composed of two parts. In the first part, we report on results of observational analysis using four year data of \fermi\ Large Area Telescope. Due to the improvement of instrumental response function and increase of the statistics, the observational uncertainties of the spectrum in $\sim$100-300 MeV bands and $>10$GeV bands are significantly improved. The present data analysis suggests that the 0.1-100GeV emissions from LS~5039 contain three different components; (i) the first component contributes to $<$1GeV emissions around superior conjunction, (ii) the second component dominates in 1-10GeV energy bands and (iii) the third component is compatible to lower energy tail of the TeV emissions. In the second part, we develop an emission model to explain the properties of the phase-resolved emissions in multi-wavelength observations. Assuming that LS~5039 includes a pulsar, we argue that both emissions from magnetospheric outer gap and inverse-Compton scattering process of cold-relativistic pulsar wind contribute to the observed GeV emissions. We assume that the pulsar is wrapped by two kinds of termination shock; Shock-I due to the interaction between the pulsar wind and the stellar wind and Shock-II due to the effect of the orbital motion. We propose that the X-rays are produced by the synchrotron radiation at Shock-I region and the TeV gamma-rays are produced by the inverse-Compton scattering process at Shock-II region.

preprint2014arXiv

High-energy observations of PSR B1259-63/LS 2883 through the 2014 periastron passage: connecting X-rays to the GeV flare

The binary system PSR B1259-63/LS 2883 is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital phase-dependent variability in these frequencies. The first detection of GeV gamma-rays from the system was made around the 2010 periastron passage. In this Letter, we present an analysis of X-ray and gamma-ray data obtained by the Swift/XRT, NuSTAR/FPM, and Fermi/LAT, through the recent periastron passage which occurred on 2014 May 4. While PSR B1259-63/LS 2883 was not detected by the LAT before and during this passage, we show that the GeV flares occurred at a similar orbital phase as in early 2011, thus establishing the repetitive nature of the post-periastron GeV flares. Multiple flares each lasting for a few days have been observed and short-term variability is seen as well. We also found X-ray flux variation contemporaneous with the GeV flare for the first time. A strong evidence of the keV-to-GeV connection came from the broadband high-energy spectra, which we interpret as synchrotron radiation from the shocked pulsar wind.

preprint2014arXiv

Identification campaign of supernova remnant candidates in the Milky Way. II. X-ray studies of G38.7-1.4

We report on XMM-Newton and Chandra observations of the Galactic supernova remnant candidate G38.7-1.4, together with complementary radio, infrared, and gamma-ray data. An approximately elliptical X-ray structure is found to be well correlated with radio shell as seen by the Very Large Array. The X-ray spectrum of G38.7-1.4 can be well-described by an absorbed collisional ionization equilibrium plasma model, which suggests the plasma is shock heated. Based on the morphology and the spectral behaviour, we suggest that G38.7-1.4 is indeed a supernova remnant belongs to a mix-morphology category.

preprint2014arXiv

Identification of the Optical Counterpart of Fermi Black Widow Millisecond Pulsar PSR J1544+4937

We report the optical identification of the companion to the {\it Fermi} black widow millisecond pulsar PSR J1544+4937. We find a highly variable source on Keck LRIS images at the nominal pulsar position, with 2 magnitude variations over orbital period in the B, g, R, and I bands. The nearly achromatic light curves are difficult to explain with a simply irradiated hemisphere model, and suggest that the optical emission is dominated by a nearly isothermal hot patch on the surface of the companion facing the pulsar. We roughly constrain the distance to PSR J1544+4937 to be between 2 and 5 kpc. A more reliable distance measurement is needed in order to constrain the composition of the companion.

preprint2014arXiv

Multi-wavelength emissions from the millisecond pulsar binary PSR J1023+0038 during an accretion active state

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

preprint2014arXiv

NuSTAR observations and broadband spectral energy distribution modeling of the millisecond pulsar binary PSR J1023+0038

We report the first hard X-ray (3-79 keV) observations of the millisecond pulsar (MSP) binary PSR J1023+0038 using NuSTAR. This system has been shown transiting between a low-mass X-ray binary (LMXB) state and a rotation-powered MSP state. The NuSTAR observations were taken in both LMXB state and rotation-powered state. The source is clearly seen in both states up to ~79 keV. During the LMXB state, the 3-79 keV flux is about a factor of 10 higher that in the rotation-powered state. The hard X-rays show clear orbital modulation during the X-ray faint rotation-powered state but the X-ray orbital period is not detected in the X-ray bright LMXB state. In addition, the X-ray spectrum changes from a flat power-law spectrum during the rotation-powered state to a steeper power-law spectrum in the LMXB state. We suggest that the hard X-rays are due to the intra-binary shock from the interaction between the pulsar wind and the injected material from the low-mass companion star. During the rotation-powered MSP state, the X-ray orbital modulation is due to Doppler boosting of the shocked pulsar wind. At the LMXB state, the evaporating matter of the accretion disk due to the gamma-ray irradiation from the pulsar stops almost all the pulsar wind, resulting the disappearance of the X-ray orbital modulation.

preprint2014arXiv

The birth of radio millisecond pulsars and their high-energy signature

Millisecond pulsars (MSPs) are thought to born in low-mass X-ray binaries when the neutron star has gained enough angular momentum from the accreting materials of its companion star. It is generally believed that a radio MSP is born when the neutron star stops accreting and enters a rotation-powered state. Exactly what happens during the transition time was poorly understood until a year ago. In the past year, observations have revealed a few objects that not only switched from one state to the other (as predicted in the above picture), but also have swung between the two states within weeks to years. In this work, we present observations of two of these transition objects (PSR J1023+0038 and XSS J12270-4859) and a theoretical framework that tries to explain their high-energy radiation.

preprint2013arXiv

Exploring the X-ray and gamma-ray properties of the redback millisecond pulsar PSR J1723-2837

We have investigated the X-ray and $γ$-ray properties of the redback millisecond pulsar PSR J1723-2837 with XMM-Newton, Chandra and Fermi. We have discovered the X-ray orbital modulation of this binary system with the minimum that coincides with the phases of radio eclipse. The X-ray emission is clearly non-thermal in nature which can be well described by a simple power-law with a photon index of $\sim1.2$. The phase-averaged luminosity is $\sim9\times10^{31}$ erg/s in 0.3-10 keV which consumes $\sim0.2\%$ of the spin-down power. We have detected the $γ-$ray emission in $0.1-300$ GeV from this system at a significance of $\sim6σ$ for the first time. The $γ-$rays in this energy range consumes $\sim2\%$ of the spin-down power and can be modeled by a power-law with a photon index of $\sim2.6$. We discuss the high energy properties of the new redback in the context of a intrabinary shock model.

preprint2013arXiv

Search for pulsed gamma-ray emission from globular cluster M28

Using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, we have searched for the gamma-ray pulsations from the direction of globular cluster M28 (NGC 6626). We report the discovery of a signal with the frequency consistent with that of the energetic millisecond pulsar (MSP) PSR B1821-24 in M28. A weighted H-test test statisic (TS) of 28.8 is attained which corresponds to a chance probability of ~1e-5 (4.3-sigma detection). With a phase-resolved analysis, the pulsed component is found to contribute ~25% of the total observed gamma-ray emission from the cluster. On the other hand, the unpulsed level provides a constraint for the underlying MSP population and the fundamental plane relations for the scenario of inverse Compton scattering. Follow-up timing observations in radio/X-ray are encouraged for further investigating this periodic signal candidate.

preprint2012arXiv

Discovery of an unidentified Fermi object as a black widow-like millisecond pulsar

The Fermi Gamma-ray Space Telescope has revolutionized our knowledge of the gamma-ray pulsar population, leading to the discovery of almost 100 gamma-ray pulsars and dozens of gamma-ray millisecond pulsars (MSPs). Although the outer-gap model predicts different sites of emission for the radio and gamma-ray pulsars, until now all of the known gamma-ray MSPs have been visible in the radio. Here we report the discovery of a "radio-quiet" gamma-ray emitting MSP candidate by using Fermi, Chandra, Swift, and optical observations. The X-ray and gamma-ray properties of the source are consistent with known gamma-ray pulsars. We also found a 4.63-hr orbital period in optical and X-ray data. We suggest that the source is a black widow-like MSP with a ~0.1 solar-mass late-type companion star. Based on the profile of the optical and X-ray light-curves, the companion star is believed to be heated by the pulsar while the X-ray emissions originate from pulsar magnetosphere and/or from intra-binary shock. No radio detection of the source has been reported yet and although no gamma-ray/radio pulsation has been found, we estimated that the spin period of the MSP is ~3-5 ms based on the inferred gamma-ray luminosity.

preprint2012arXiv

Discovery of γ-ray pulsation and X-ray emission from the black widow pulsar PSR J2051-0827

We report the discovery of pulsed γ-ray emission and X-ray emission from the black widow millisecond pulsar PSR J2051-0827 by using the data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope and the Advanced CCD Imaging Spectrometer array (ACIS-S) on the Chandra X-ray Observatory. Using 3 years of LAT data, PSR J2051-0827 is clearly detected in γ-ray with a signicance of \sim 8σ in the 0.2 - 20 GeV band. The 200 MeV - 20 GeV γ-ray spectrum of PSR J2051-0827 can be modeled by a simple power- law with a photon index of 2.46 \pm 0.15. Significant (\sim 5σ) γ-ray pulsations at the radio period were detected. PSR J2051-0827 was also detected in soft (0.3-7 keV) X-ray with Chandra. By comparing the observed γ-rays and X-rays with theoretical models, we suggest that the γ-ray emission is from the outer gap while the X-rays can be from intra-binary shock and pulsar magnetospheric synchrotron emissions.

preprint2012arXiv

Identification campaign of supernova remnant candidates in the Milky Way - I: Chandra observation of G308.3-1.4

ROSAT all-sky survey (RASS) data have provided another window to search for supernova remnants (SNRs). In reexamining this data archive, a list of unidentified extended X-ray objects have been suggested as promising SNR candidate. However, most of these targets have not yet been fully explored by the state-of-art X-ray observatories. For selecting a pilot target for a long-term identification campaign, we have observed the brightest candidate, G308.3-1.4, with Chandra X-ray observatory. An incomplete shell-like X-ray structure which well-correlated with the radio shell emission at 843 MHz has been revealed. The X-ray spectrum suggests the presence of a shock-heated plasma. All these evidences confirm G308.3-1.4 as a SNR. The brightest X-ray point source detected in this field-of-view is also the one locates closest to the geometrical center of G308.3-1.4, which has a soft spectrum. The intriguing temporal variability and the identification of optical/infrared counterpart rule out the possibility of an isolated neutron star. On the other hand, the spectral energy distribution from Ks band to R band suggests a late-type star. Together with a putative periodicity of \sim1.4 hrs, the interesting excesses in V, B bands and H-alpha suggest this source as a promising candidate of a compact binary survived in a supernova explosion (SN).

preprint2012arXiv

Orbital Phase Dependent Gamma-ray emissions from the Black Widow Pulsar

We report on evidence for orbital phase-dependence of the gamma-ray emission from PSR B1957+20 black widow system by using the data of the Fermi Large Area Telescope. We divide an orbital cycle into two regions: a region containing the inferior conjunction, and the other region containing rest of the orbital cycle. We show that the observed spectra for the different orbital regions are fitted by different functional forms. The spectrum of the orbital region containing inferior conjunction can be described by a power-law with an exponential cutoff (PLE) model, which gives the best-fit model for the orbital phase that does not contain the inferior conjunction, plus an extra component above ~2.7 GeV. The emission above 3 GeV in this region is detected with a ~7-sigma confidence level. The gamma-ray data above ~2.7 GeV are observed to be modulated at the orbital period at the ~2.3-sigma level. We anticipate that the PLE component dominating below ~2.7 GeV originates from the pulsar magnetosphere. We also show that the inverse-Compton scattering of the thermal radiation of the companion star off a "cold" ultra-relativistic pulsar wind can explain the extra component above ~2.7 GeV. The black widow pulsar PSR B1957+20 may be the member of a new class of object, in the sense that the system is showing gamma-ray emission with both magnetospheric and pulsar wind origins.

preprint2012arXiv

Pulsed γ-ray emission from magnetar 1E 2259+586

Anomalous x-ray pulsars (AXPs) are thought to be magnetars which are young isolated neutron stars with extremely strong magnetic fields of > 10^14Gauss. Their tremendous magnetic fields inferred from the spin parameters provide a huge energy reservoir to power the observed x-ray emission. High-energy emission above 0.3 MeV has never been detected despite intensive search. Here, we present the possible Fermi Large Area Telescope (LAT) detection of γ-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of γ-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy γ-rays originate from the outer magnetosphere of the magnetar.

preprint2012arXiv

The absence of radio emission from the globular cluster G1

The detections of both X-ray and radio emission from the cluster G1 in M31 have provided strong support for existing dynamical evidence for an intermediate mass black hole (IMBH) of mass 1.8 +/- 0.5 x 10^4 solar masses at the cluster center. However, given the relatively low significance and astrometric accuracy of the radio detection, and the non-simultaneity of the X-ray and radio measurements, this identification required further confirmation. Here we present deep, high angular resolution, strictly simultaneous X-ray and radio observations of G1. While the X-ray emission (L_X = 1.74^{+0.53}_{-0.44} x 10^{36} (d/750 kpc)^2 erg/s in the 0.5-10 keV band) remained fully consistent with previous observations, we detected no radio emission from the cluster center down to a 3-sigma upper limit of 4.7 microJy/beam. Our favored explanation for the previous radio detection is flaring activity from a black hole low mass X-ray binary (LMXB). We performed a new regression of the Fundamental Plane of black hole activity, valid for determining black hole mass from radio and X-ray observations of sub-Eddington black holes, finding log M_{BH} = (1.638 +/- 0.070)log L_R - (1.136 +/- 0.077)log L_X - (6.863 +/- 0.790), with an empirically-determined uncertainty of 0.44 dex. This constrains the mass of the X-ray source in G1, if a black hole, to be <9700 solar masses at 95% confidence, suggesting that it is a persistent LMXB. This annuls what was previously the most convincing evidence from radiation for an IMBH in the Local Group, though the evidence for an IMBH in G1 from velocity dispersion measurements remains unaffected by these results.

preprint2012arXiv

The Discovery of an X-ray/UV Stellar Flare from the Late-K/Early-M Dwarf LMC 335

We report the discovery of an X-ray/UV stellar flare from the source LMC 335, captured by XMM-Newton in the field of the Large Magellanic Cloud. The flare event was recorded continuously in X-ray for its first 10 hours from the precursor to the late decay phases. The observed fluxes increased by more than two orders of magnitude at its peak in X-ray and at least one in the UV as compared to quiescence. The peak 0.1-7.0 keV X-ray flux is derived from the two-temperature APEC model to be ~(8.4 +/- 0.6) x 10^-12 erg cm-2 s-1. Combining astrometric information from multiple X-ray observations in the quiescent and flare states, we identify the NIR counterpart of LMC 335 as the 2MASS source J05414534-6921512. The NIR color relations and spectroscopic parallax characterize the source as a Galactic K7-M4 dwarf at a foreground distance of (100 - 264) pc, implying a total energy output of the entire event of ~(0.4 - 2.9) x 10^35 erg. This report comprises detailed analyses of this late-K / early-M dwarf flare event that has the longest time coverage yet reported in the literature. The flare decay can be modeled with two exponential components with timescales of ~28 min and ~4 hours, with a single component decay firmly ruled out. The X-ray spectra during flare can be described by two components, a dominant high temperature component of ~40-60MK and a low temperature component of ~10MK, with a flare loop length of about 1.1-1.3 stellar radius.

preprint2012arXiv

X-ray studies of the Black Widow Pulsar PSR B1957+20

We report on Chandra observations of the black widow pulsar, PSR B1957+20. Evidence for a binary-phase dependence of the X-ray emission from the pulsar is found with a deep observation. The binary-phase resolved spectral analysis reveals non-thermal X-ray emission of PSR B1957+20, confirming the results of previous studies. This suggests that the X-rays are mostly due to intra-binary shock emission which is strongest when the pulsar wind interacts with the ablated material from the companion star. The geometry of the peak emission is determined in our study. The marginal softening of the spectrum of the non-thermal X-ray tail may indicate that particles injected at the termination shock is dominated by synchrotron cooling.

preprint2011arXiv

Chandra Detection of SN 2010da Four Months After Outburst: Evidence for a Supergiant X-ray Binary in NGC 300

We present the results of a 63 ks {\it Chandra} observation of the "supernova impostor" SN 2010da four months after it was first observed on 25 May 2010. We detect an X-ray source at $\sim7σ$ confidence coincident with the optical location of the SN 2010da outburst. Our imaging analysis has revealed a hard central point source, surrounded by soft diffuse emission extending as far as $\sim$8\asn north of the central source. The diffuse emission has a hardness ratio, 0.35-2 keV luminosity ($\sim6\times10^{35}$ erg s$^{-1}$), and size ($\sim20$ pc) consistent with that of a supernova remnant, although the low number of counts prohibits detailed spectral modeling. Our best-fit spectral model for the hard central source is a black body ($kT = 1.79^{+0.66}_{-0.43}$ keV) with no evidence for intrinsic absorption beyond the Galactic column. We estimate the 0.3-10 keV luminosity to be 1.7$^{+0.2}_{-0.5}\times10^{37}$ erg s$^{-1}$, a factor of $\sim$25 decrease since the initial outburst four months previously. The high X-ray luminosity and slow fading rate is not consistent with this object being a single massive star undergoing an outburst; instead, we favor the scenario where the massive star powering the SN 2010da optical transient is part of a wind-fed supergiant X-ray binary system with a compact companion powering the observed X-ray emission.

preprint2011arXiv

Discovery of gamma-ray emission from the supernova remnant Kes 17 with Fermi Large Area Telescope

We report the discovery of GeV emission at the position of supernova remnant Kes 17 by using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Kes 17 can be clearly detected with a significance of ~12 sigma in the 1 - 20 GeV range. Moreover, a number of gamma-ray sources were detected in its vicinity. The gamma-ray spectrum of Kes 17 can be well described by a simple power-law with a photon index of ~ 2.4. Together with the multi-wavelength evidence for its interactions with the nearby molecular cloud, the gamma-ray detection suggests that Kes 17 is a candidate acceleration site for cosmic-rays.

preprint2011arXiv

Exploring the dark accelerator HESS J1745-303 with Fermi Large Area Telescope

We present a detailed analysis of the gamma-ray emission from HESS J1745-303 with the data obtained by the Fermi Gamma-ray Space Telescope in the first ~29 months observation.The source can be clearly detected at the level of ~18-sigma and ~6-sigma in 1-20 GeV and 10-20 GeV respectively. Different from the results obtained by the Compton Gamma-ray Observatory, we do not find any evidence of variability. Most of emission in 10-20 GeV is found to coincide with the region C of HESS J1745-303. A simple power-law is sufficient to describe the GeV spectrum with a photon index of ~2.6. The power-law spectrum inferred in the GeV regime can be connected to that of a particular spatial component of HESS J1745-303 in 1-10 TeV without any spectral break. These properties impose independent constraints for understanding the nature of this "dark particle accelerator".

preprint2011arXiv

Fractal dimension and thermodynamic fluctuation properties of IDV light curves

Fractals are a basic tool to phenomenologically describe natural objects having a high degree of temporal or spatial variability. From a physical point of view the fractal properties of natural systems can also be interpreted by using the standard formalism of thermodynamical fluctuations. In the present paper we introduce and analyze the fractal dimension of the Intra - Day- Variability (IDV) light-curves of the BL Lac objects, in the optical, radio and X-ray bands, respectively. A general description of the fluctuation spectrum of these systems based on general thermodynamical principles is also proposed. Based on the general fractal properties of a given physical system, we also introduce the predictability index for the IDV light curves. We have explicitly determined the fractal dimension for the R-band observations of five blazars, as well as for the radio band observations of the compact extragalactic radio source J 1128+5925, and of several X-ray sources. Our results show that the fractal dimension of the optical and X-ray observations indicates an almost pure "Brownian noise" (random walk) spectrum, with a very low predictability index, while in the radio band the predictability index is much higher. We have also studied the spectral properties of the IDV light curves, and we have shown that their spectral index is very closely correlated with the corresponding fractal dimension.

preprint2011arXiv

Gamma-ray emission from the globular clusters Liller 1, M80, NGC 6139, NGC 6541, NGC 6624, and NGC 6752

Globular clusters (GCs) are emerging as a new class of gamma-ray emitters, thanks to the data obtained from the Fermi Gamma-ray Space Telescope. By now, eight GCs are known to emit gamma-rays at energies >100~MeV. Based on the stellar encounter rate of the GCs, we identify potential gamma-ray emitting GCs out of all known GCs that have not been studied in details before. In this paper, we report the discovery of a number of new gamma-ray GCs: Liller 1, NGC 6624, and NGC 6752, and evidence for gamma-ray emission from M80, NGC 6139, and NGC 6541, in which gamma-rays were found within the GC tidal radius. With one of the highest metallicity among all GCs in the Milky Way, the gamma-ray luminosity of Liller 1 is found to be the highest of all known gamma-ray GCs. In addition, we confirm a previous report of significant gamma-ray emitting region next to NGC 6441. We briefly discuss the observed offset of gamma-rays from some GC cores. The increasing number of known gamma-ray GCs at distances out to ~10 kpc is important for us to understand the gamma-ray emitting mechanism and provides an alternative probe to the underlying millisecond pulsar populations of the GCs.

preprint2011arXiv

The deep XMM-Newton Survey of M 31

Aims. M 31 has been completely imaged for the first time obtaining a luminosity lower limit ~10^35 erg s-1 in the 0.2-4.5 keV band. Our XMM-Newton EPIC survey combines archival observations with observations taken between June 2006 and February 2008. The main goal of the paper is a study of the X-ray source population of M 31. Methods. An X-ray catalogue of 1897 sources was created; 914 were detected for the first time. Source classification and identification were based on X-ray hardness ratios, spatial extent of the sources, and by cross correlating with catalogues in the X-ray, optical, infrared and radio wavelengths. We also analysed the long-term variability of the X-ray sources. Results. Fourteen of the 30 supersoft source candidates represent supersoft emission of optical novae. Many of the 25 supernova remnants (SNRs) and 31 SNR candidates lie within the 10 kpc dust ring and other star forming regions in M 31; this connection between SNRs and star forming regions implies that most of the remnants originate in type II supernovae. The brightest sources in X-rays in M 31 belong to the class of X-ray binaries (XRBs). Ten low mass XRBs (LMXBs) and 26 LMXB candidates were identified based on their temporal variability. In addition 36 LMXBs and 17 LMXB candidates were identified due to correlations with globular clusters and globular cluster candidates. From optical and X-ray colour-colour diagrams, possible high mass XRB (HMXB) candidates were selected. Two of these candidates have an X-ray spectrum as is expected for an HMXB containing a neutron star primary. Conclusions. While our survey has greatly improved our understanding of the X-ray source populations in M 31, at this point 65% of the sources can still only be classified as "hard" sources. Deeper observations in X-ray and at other wavelengths would help classify these sources. ABRIDGED

preprint2011arXiv

The Fundamental Plane of Gamma-ray Globular Clusters

We have investigated the properties of a group of $γ$-ray emitting globular clusters (GCs) which have recently been uncovered in our Galaxy. By correlating the observed $γ$-ray luminosities $L_γ$ with various cluster properties, we probe the origin of the high energy photons from these GCs. We report $L_γ$ is positively correlated with the encounter rate $Γ_{c}$ and the metalicity $\left[{\rm Fe/H}\right]$ which place an intimate link between the gamma-ray emission and the millisecond pulsar population. We also find a tendency that $L_γ$ increase with the energy densities of the soft photon at the cluster location. Furthermore, the two-dimensional regression analysis suggests that $L_γ$, soft photon densities, and $Γ_{c}$/$\left[{\rm Fe/H}\right]$ possibly span fundamental planes which potentially provide better predictions for the $γ$-ray properties of GCs.

preprint2010arXiv

Cataclysmic Variables and Other Compact Binaries in the Globular Cluster NGC 362: Candidates from Chandra and HST

Highly sensitive and precise X-ray imaging from Chandra, combined with the superb spatial resolution of HST optical images, dramatically enhances our empirical understanding of compact binaries such as cataclysmic variables and low mass X-ray binaries, their progeny, and other stellar X-ray source populations deep into the cores of globular clusters. Our Chandra X-ray images of the globular cluster NGC 362 reveal 100 X-ray sources, the bulk of which are likely cluster members. Using HST color-magnitude and color-color diagrams, we quantitatively consider the optical content of the NGC 362 Chandra X-ray error circles, especially to assess and identify the compact binary population in this condensed-core globular cluster. Despite residual significant crowding in both X-rays and optical, we identify an excess population of Hα-emitting objects that is statistically associated with the Chandra X-ray sources. The X-ray and optical characteristics suggest that these are mainly cataclysmic variables, but we also identify a candidate quiescent low mass X-ray binary. A potentially interesting and largely unanticipated use of observations such as these may be to help constrain the macroscopic dynamic state of globular clusters.

preprint2010arXiv

Fermi Discovery of Gamma-ray Emission from the Globular Cluster Terzan 5

We report the discovery of gamma-ray emission from the Galactic globular cluster Terzan 5 using data taken with the Fermi Gamma-ray Space Telescope, from 2008 August 8 to 2010 January 1. Terzan 5 is clearly detected in the 0.5-20 GeV band by Fermi at ~27σlevel. This makes Terzan 5 as the second gamma-ray emitting globular cluster seen by Fermi after 47 Tuc. The energy spectrum of Terzan 5 is best represented by an exponential cutoff power-law model, with a photon index of ~1.9 and a cutoff energy at ~3.8 GeV. By comparing to 47 Tuc, we suggest that the observed gamma-ray emission is associated with millisecond pulsars, and is either from the magnetospheres or inverse Compton scattering between the relativistic electrons/positrons in the pulsar winds and the background soft photons from the Galactic plane. Furthermore, it is suggestive that the distance to Terzan 5 is less than 10 kpc and > 10 GeV photons can be seen in the future.

preprint2010arXiv

Localization of the X-ray Source in the Globular Cluster G1 with Chandra

We report the most accurate X-ray position of the X-ray source in the giant globular cluster G1 in M31 by using the Chandra X-ray Observatory, Hubble Space Telescope (HST), and Canada-France-Hawaii Telescope (CFHT). G1 is clearly detected with Chandra and by cross-registering with HST and CFHT images, we derive a 1-sigma error radius of 0.15 arcsec, significantly smaller than the previous measurement by XMM-Newton. We conclude that the X-ray emission of G1 is likely to come from within the core radius of the cluster. We have considered a number of possibilities for the origin of the X-ray emission but can rule all but two scenarios out: it could be due to either accretion onto a central intermediate-mass black hole (IMBH), or an ordinary low-mass X-ray binary (LMXB). Based on the X-ray luminosity and the Bondi accretion rate, an IMBH accreting from the cluster gas seems unlikely and we suggest that the X-rays are due to accretion from a companion. Alternatively, the probability that a 1.5 solar masses cluster LMXB lies within the 95 per cent X-ray error circle is about 0.7. Therefore we cannot rule out a single LMXB as the origin of the X-ray emission. While we cannot distinguish between different models with current observations, future high-resolution and high-sensitivity radio imaging observations will reveal whether there is an IMBH at the centre of G1.

preprint2010arXiv

Swift observations of the ultraluminous X-ray source Holmberg IX X-1

Holmberg IX X-1 is a well-known ultraluminous X-ray source with an X-ray luminosity of ~1e40 erg/s. The source has been monitored by the X-ray Telescope of Swift regularly. Since 2009 April, the source has been in an extended low luminosity state. We utilize the co-added spectra taken at different luminosity states to study the spectral behavior of the source. Simple power-law and multi-color disk blackbody models can be ruled out. The best overall fits, however, are provided by a dual thermal model with a cool blackbody and a warm disk blackbody. This suggests that Holmberg IX X-1 may be a 10 solar-mass black hole accreting at 7 times above the Eddington limit or a 100 solar-mass maximally rotating black hole accreting at the Eddington limit, and we are observing both the inner regions of the accretion disk and outflows from the compact object.

preprint2010arXiv

The Origin of Gamma-Rays from Globular Clusters

Fermi has detected gamma-ray emission from eight globular clusters. We suggest that the gamma-ray emission from globular clusters may result from the inverse Compton scattering between relativistic electrons/positrons in the pulsar wind of MSPs in the globular clusters and background soft photons including cosmic microwave/relic photons, background star lights in the clusters, the galactic infrared photons and the galactic star lights. We show that the gamma-ray spectrum from 47 Tuc can be explained equally well by upward scattering of either the relic photons, the galactic infrared photons or the galactic star lights whereas the gamma-ray spectra from other seven globular clusters are best fitted by the upward scattering of either the galactic infrared photons or the galactic star lights. We also find that the observed gamma-ray luminosity is correlated better with the combined factor of the encounter rate and the background soft photon energy density. Therefore the inverse Compton scattering may also contribute to the observed gamma-ray emission from globular clusters detected by Fermi in addition to the standard curvature radiation process. Furthermore, we find that the emission region of high energy photons from globular cluster produced by inverse Compton scattering is substantially larger than the core of globular cluster with a radius >10pc. The diffuse radio and X-rays emitted from globular clusters can also be produced by synchrotron radiation and inverse Compton scattering respectively. We suggest that future observations including radio, X-rays, and gamma-rays with energy higher than 10 GeV and better angular resolution can provide better constraints for the models.

preprint2009arXiv

The first two transient supersoft X-ray sources in M 31 globular clusters and the connection to classical novae

Classical novae (CNe) have been found to represent the major class of supersoft X-ray sources (SSS) in our neighbour galaxy M 31. We determine properties and evolution of the two first SSSs ever discovered in the M 31 globular cluster (GC) system. We have used XMM-Newton, Chandra and Swift observations of the centre region of M 31 to discover both SSS and to determine their X-ray light curves and spectra. We performed detailed analysis of XMM-Newton EPIC PN spectra of the source in Bol 111 (SS1) using blackbody and NLTE white dwarf (WD) atmosphere models. For the SSS in Bol 194 (SS2) we used optical monitoring data to search for an optical counterpart. Both GC X-ray sources were classified as SSS. We identify SS1 with the CN M31N 2007-06b recently discovered in the M 31 GC Bol 111. For SS2 we did not find evidence for a recent nova outburst and can only provide useful constraints on the time of the outburst of a hypothetical nova. The only known CN in a M 31 GC can be identified with the first SSS found in a M31 GC. We discuss the impact of our observations on the nova rate for the M 31 GC system.

preprint2008arXiv

An extremely luminous X-ray outburst at the birth of a supernova

Massive stars end their short lives in spectacular explosions, supernovae, that synthesize new elements and drive galaxy evolution. Throughout history supernovae were discovered chiefly through their delayed optical light, preventing observations in the first moments (hours to days) following the explosion. As a result, the progenitors of some supernovae and the events leading up to their violent demise remain intensely debated. Here we report the serendipitous discovery of a supernova at the time of explosion, marked by an extremely luminous X-ray outburst. We attribute the outburst to the break-out of the supernova shock-wave from the progenitor, and show that the inferred rate of such events agrees with that of all core-collapse supernovae. We forecast that future wide-field X-ray surveys will catch hundreds of supernovae each year in the act of explosion, and thereby enable crucial neutrino and gravitational wave detections that may ultimately unravel the explosion mechanism.

preprint2007arXiv

GRB 070201: A possible Soft Gamma Ray Repeater in M31

The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures.

preprint2007arXiv

The Spectrum of the Black Hole X-ray Nova V404 Cygni in Quiescence as Measured by XMM-Newton

We present XMM observations of the black hole X-ray nova V404 Cygni in quiescence. Its quiescent spectrum can be best fitted by a simple power-law with slope 2. The spectra are consistent with that expected for the advection-dominated accretion flow (ADAF). V404 Cyg was roughly equal in luminosity compared to the previous observation of Chandra. We see variability of a factor of 4 during the observation. We find no evidence for the presence of fluorescent or H-like/He-like iron emission, with upper limits of 52 eV and 110 eV respectively. The limit on the fluorescent emission is improved by a factor of 15 over the previous estimate, and the restriction on H-like/He-like emission is lower than predicted from models by a factor of roughly 2.

preprint2006arXiv

Chandra and Hubble Space Telescope Study of the Globular Cluster NGC 288

We report on the Chandra X-ray Observatory observations of the globular cluster NGC 288. We detect four X-ray sources within the core radius and seven additional sources within the half-mass radius down to a limiting luminosity of Lx=7e30 erg/s (assuming cluster membership) in the 0.3-7 keV band. We also observed the cluster with the Hubble Space Telescope Advanced Camera for Surveys and identify optical counterparts to seven X-ray sources out of the nine sources within the HST field-of-view. Based on the X-ray and optical properties, we find 2-5 candidates of cataclysmic variables (CVs) or chromospherically active binaries, and 2-5 background galaxies inside the half-mass radius. Since the core density of NGC 288 is very low, the faint X-ray sources of NGC 288 found in the Chandra and HST observations is higher than the prediction on the basis of the collision frequency. We suggest that the CVs and chromospherically active binaries are primordial in origin, in agreement with theoretical expectation.

preprint2005arXiv

A "Combination Nova" Outburst in Z Andromedae: Nuclear Shell Burning Triggered by a Disk Instability

We describe observational evidence for a new kind of interacting-binary-star outburst that involves both an accretion instability and an increase in thermonuclear shell burning on the surface of an accreting white dwarf. We refer to this new type of eruption as a combination nova. In late 2000, the prototypical symbiotic star Z Andromedae brightened by roughly two magnitudes in the optical. We observed the outburst in the radio with the VLA and MERLIN, in the optical both photometrically and spectroscopically, in the far ultraviolet with FUSE, and in the X-rays with both Chandra and XMM. The two-year-long event had three distinct stages. During the first stage, the optical rise closely resembled an earlier, small outburst that was caused by an accretion-disk instability. In the second stage, the hot component ejected an optically thick shell of material. In the third stage, the shell cleared to reveal a white dwarf whose luminosity remained on the order of 10^4 Lsun for approximately one year. The eruption was thus too energetic to have been powered by accretion alone. We propose that the initial burst of accretion was large enough to trigger enhanced nuclear burning on the surface of the white dwarf and the ejection of an optically thick shell of material. This outburst therefore combined elements of both a dwarf nova and a classical nova. Our results have implications for the long-standing problem of producing shell flashes with short recurrence times on low-mass white dwarfs in symbiotic stars.

preprint2002arXiv

Bright X-Ray Sources in M31 Globular Clusters

We have conducted Chandra observations of ~2560 square arcmin (~131 kpc^2) of of M31, and find that the most luminous X-ray sources in our fields are in globular clusters (GCs). Of the 28 GC X-ray sources in our fields, 17 are newly discovered. Approximately 1/3 of all the sources have L_x([0.5-7] keV) > 10^{37} ergs/s; approximately 1/10 of all the sources have L_x([0.5-7] keV) close to or above 10^{38} erg/s. The most luminous source, in the globular cluster Bo 375, is consistently observed to have L_x greater than 2 x 10^{38} ergs/s. (1) We present data on the spectra and/or light curves of the 5 most luminous M31 GC sources. (2)We explore possible explanations for the high X-ray luminosities of the brightest sources. These include that the X-ray sources may be composites, the radiation we receive may be beamed, metallicity effects could be at work, or the sources may be accreting black holes. We weigh each of these possibilities against the data. In addition, we introduce a neutron star model in which mass transfer proceeds on the thermal time scale of the donor star. Our model can produce luminosities of several times 10^{38} ergs/s, and leads to a set of well-defined predictions. (3) We compute the X-ray luminosity function and the distribution of counts in wavebands that span the range of energies to which Chandra is sensitive. We find the peak X-ray luminosity is higher and that systems with L_x > 10^{37} ergs/s constitute a larger fraction of all GC sources than in the MW. (4) We study the possible reasons for this difference between M31 and Galactic GC X-ray sources and identify three promising explanations.