BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shuang-Nan Zhang

Shuang-Nan Zhang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.IMastro-ph.SRastro-ph.COastro-ph.GAgr-qchep-phastro-phphysics.ins-detastro-ph.EPhep-exhep-thphysics.ed-ph

Catalog footprint

What is connected

111works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

111 published item(s)

preprint2026arXiv

Brightest GRB flare observed in GRB 221009A: bridge the last gap between flare and prompt emission in GRB

Flares are usually observed during the afterglow phase of Gamma-Ray Bursts (GRBs) in soft X-ray, optical and radio bands, but rarely in gamma-ray band. Despite the extraordinary brightness, GECAM-C has accurately measured both the bright prompt emission and flare emission of GRB 221009A without instrumental effects, offering a good opportunity to study the relation between them. In this work, we present a comprehensive analysis of flare emission of GRB 221009A, which is composed of a series of flares. Among them, we identify an exceptionally bright flare with a record-breaking isotropic energy $E_{\rm iso} = 1.82 \times 10^{53}$ erg of GRB flares. It exhibits the highest peak energy ever detected in GRB flares, $E_{\rm peak} \sim 300$ keV, making it a genuine gamma-ray flare. It also shows rapid rise and decay timescales, significantly shorter than those of typical X-ray flares observed in soft X-ray or optical band, but comparable to those observed in prompt emissions. Despite these exceptional properties, the flare shares several common properties with typical GRB flares. We note that this is the first observation of a GRB flare in the keV-MeV band with sufficiently high temporal resolution and high statistics, which bridges the last gap between prompt emission and flare.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Detection of Oscillations in a Type I X-Ray Burst of 4U 0614+091 with SVOM/ECLAIRs

On 2025 January 10, a thermonuclear (Type I) X-ray burst from the neutron star low-mass X-ray binary \textit{4U~0614+091} was detected with the ECLAIRs instrument on board the \textit{SVOM} mission. We present here a time-resolved spectroscopic analysis of the burst, along with the detection of burst oscillations within a 51-second interval during the decay phase. The oscillation frequency is measured to be $ν= 413.674 \pm 0.002\,\mathrm{Hz}$, consistent with previous reports. However, we detect a significant downward frequency drift over the burst duration, characterized by $\dotν = (-4.7 \pm 0.3) \times 10^{-3}\,\mathrm{Hz\,s^{-1}}$. This frequency evolution is atypical compared to those observed in similar burst oscillation sources. We tentatively attribute the observed drift to a Doppler shift induced by orbital motion. Under this interpretation, the inferred orbital period must be shorter than 20 minutes, placing \textit{4U~0614+091} among the most compact known low-mass X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

On the Ultra-Long Gamma-Ray Transient GRB 250702B/EP250702

GRB 250702B/EP250702a is an interesting long-duration gamma-ray transient whose nature is in debate. To obtain a full picture in gamma-ray band, we implement a comprehensive targeted search of burst emission in a wide window of 30 days jointly with Insight-HXMT, GECAM and Fermi/GBM data within the ETJASMIN framework. In gamma-ray band, we find there is a 50-second precursor about 25 hours before the 4-hour main burst, which generally consists of 4 emission episodes. Remarkably, we find that the soft X-ray emission (after the main burst) decays as a power-law with start time aligning with the last episode of main emission and index of -5/3 perfectly consistent with the canonical prediction of fallback accretion. We conclude that the properties of precursor, main burst and the following soft X-ray emission strongly support the atypical collapsar Ultra-Long Gamma-Ray Burst (ULGRB) scenario rather than the Tidal Disruption Event (TDE), and all these gamma-ray and soft X-ray emission probably originate from relativistic jet whose luminosity is dominated by the fallback accretion rate during the death collapse of a supergiant star.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

X-ray and radio observations of the AMXP MAXI J1957+032 covering the 2022-2025 outbursts

We presented a comprehensive multi-epoch timing and multiwavelength analysis of the accreting millisecond X-ray pulsar MAXI J1957+032, covering two major outbursts in 2022 and 2025. By reanalyzing the 2022 outburst data from the Neutron Star Interior Composition Explorer (NICER), we found the spin frequency and orbital parameters from the observations in 0.3-5 keV. For the 2025 outburst, we reported the detection of pulsations with the Einstein Probe (EP). Based on the $\sim$3-year baseline between these two outbursts, we measured a significant long-term spin-down rate of $\dotν= (-5.73 \pm 0.28) \times 10^{-14}~{\rm Hz~s^{-1}}$. Assuming that the quiescent spin-down is driven by magnetic dipole radiation, we inferred a spin-down luminosity of $L \approx 1.1 \times 10^{36}~{\rm erg~s^{-1}}$ and a surface dipolar magnetic field of $B \approx (7.3 - 10.4) \times 10^8$ G. Furthermore, we conducted a deep radio pulsation search with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during the X-ray quiescent state in 2024, resulting in a non-detection with a 7$σ$ flux density upper limit of 12.3 $μ$Jy. This corresponds to a radio efficiency upper limit of $ξ< 2.8 \times 10^{-10}$, which is significantly lower than that of typical millisecond pulsars with a similar spin-down power. This profound radio pulsation faintness can be explained by two primary scenarios: either a geometric effect, wherein the pulsar's radio beam is directed away from our line of sight, or a physical suppression of the emission mechanism, potentially caused by a persistent low-level accretion flow during the X-ray quiescent state.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Detection of disk-jet co-precession in a tidal disruption event

Theories and simulations predict that intense spacetime curvature near black holes bends the trajectories of light and matter, driving disk and jet precession under relativistic torques. However, direct observational evidence of disk-jet co-precession remains elusive. Here, we report the most compelling case to date: a tidal disruption event (TDE) exhibiting unprecedented 19.6-day quasi-periodic variations in both X-rays and radio, with X-ray amplitudes exceeding an order of magnitude. The nearly synchronized X-ray and radio variations suggest a shared mechanism regulating the emission regions. We demonstrate that a disk-jet Lense-Thirring precession model successfully reproduces these variations while requiring a low-spin black hole. This study uncovers previously uncharted short-term radio variability in TDEs, highlights the transformative potential of high-cadence radio monitoring, and offers profound insights into disk-jet physics.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

RMS-flux slope in MAXI J1820+070: a measure of the disk-corona coupling

Linear RMS-flux relation has been well established in different spectral states of all accreting systems. In this work, we study the evolution of the frequency-dependent RMS-flux relation of MAXI J1820+070 during the initial decaying phase of the 2018 outburst with Insight-HXMT over a broad energy range 1-150 keV. As the flux decreases, we first observe a linear RMS-flux relation at frequencies from 2 mHz to 10 Hz, while such a relation breaks at varying times for different energies, leading to a substantial reduction in the slope. Moreover, we find that the low-frequency variability exhibits the highest sensitivity to the break, which occurs prior to the hard-to-hard state transition time determined through time-averaged spectroscopy, and the time deviation increases with energy. The overall evolution of the RMS-flux slope and intercept suggests the presence of a two-component Comptonization system. One component is radially extended, explaining the strong disk-corona coupling before the break, while the other component extends vertically, contributing to the reduction of the disk-corona coupling after the break. A further vertical expansion of the latter component is required to accommodate the dynamic evolution observed in the RMS-flux slope. In conclusion, we suggest that the RMS-flux slope in 1-150 keV band can be employed as an indicator of the disk-corona coupling and the hard-to-hard state transition in MAXI J1820+070 could be partially driven by the changes in the corona geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A 2-hr binary period for the black hole transient MAXI J0637-430

We revisit various sets of published results from X-ray and optical studies of the Galactic black hole (BH) candidate MAXI J0637-430, which went into outburst in 2019. Combining the previously reported values of peak outburst luminosity, best-fitting radii of inner and outer accretion disk, viewing angle, exponential decay timescale and peak-to-peak separation of the He II 4686 disk emission line, we improve the constraints on the system parameters. We estimate a heliocentric distance d = (8.7 +/- 2.3) kpc, a projected Galactocentric distance R = (13.2 +/- 1.8) kpc and a height |z| = (3.1 +/- 0.8) kpc from the Galactic plane. It is the currently known Milky Way BH candidate located farthest from the Galactic Centre. We infer a BH mass M_1 = (5.1 +/- 1.6) M_{sun}, a spin parameter a* <~ 0.25, a donor star mass M_2 = (0.25 +/- 0.07) M_{sun}, a peak Eddington ratio lambda = 0.17 +/- 0.11 and a binary period P_{orb} = 2.2^{+0.8}_{-0.6} hr. This is the shortest period measured or estimated so far for any Galactic BH X-ray binary. If the donor star is a main-sequence dwarf, such a period corresponds to the evolutionary stage where orbital shrinking is driven by gravitational radiation and the star has regained contact with its Roche lobe (low end of the period gap). The three Galactic BHs with the shortest period (<~3 hr) are also those with the highest vertical distance from the Galactic plane (>~2 kpc). This is probably because binaries with higher binding energies can survive faster natal kicks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CATCH: Chasing All Transients Constellation Hunters Space Mission

In time-domain astronomy, a substantial number of transients will be discovered by multi-wavelength and multi-messenger observatories, posing a great challenge for follow-up capabilities. We have thus proposed an intelligent X-ray constellation, the Chasing All Transients Constellation Hunters (CATCH) space mission. Consisting of 126 micro-satellites in three types, CATCH will have the capability to perform follow-up observations for a large number of different types of transients simultaneously. Each satellite in the constellation will carry lightweight X-ray optics and use a deployable mast to increase the focal length. The combination of different optics and detector systems enables different types of satellites to have multiform observation capabilities, including timing, spectroscopy, imaging, and polarization. Controlled by the intelligent system, different satellites can cooperate to perform uninterrupted monitoring, all-sky follow-up observations, and scanning observations with a flexible field of view (FOV) and multi-dimensional observations. Therefore, CATCH will be a powerful mission to study the dynamic universe. Here, we present the current design of the spacecraft, optics, detector system, constellation configuration and observing modes, as well as the development plan.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Determination of QPO properties in the presence of strong broad-band noise: a case study on the data of MAXI J1820+070

Accurate calculation of the phase lags of quasi-periodic oscillations (QPOs) will provide insight into their origin. In this paper we investigate the phase lag correction method which has been applied to calculate the intrinsic phase lags of the QPOs in MAXI J1820+070. We find that the traditional additive model between BBN and QPOs in the time domain is rejected, but the convolution model is accepted. By introducing a convolution mechanism in the time domain, the Fourier cross-spectrum analysis shows that the phase lags between QPOs components in different energy bands will have a simple linear relationship with the phase lags between the total signals, so that the intrinsic phase lags of the QPOs can be obtained by linear correction. The power density spectrum (PDS) thus requires a multiplicative model to interpret the data. We briefly discuss a physical scenario for interpreting the convolution. In this scenario, the corona acts as a low-pass filter, the Green's function containing the noise is convolved with the QPOs to form the low-frequency part of the PDS, while the high-frequency part requires an additive component. We use a multiplicative PDS model to fit the data observed by Insight-HXMT. The overall fitting results are similar compared to the traditional additive PDS model. Neither the width nor the centroid frequency of the QPOs obtained from each of the two PDS models were significantly different, except for the r.m.s. of the QPOs. Our work thus provides a new perspective on the coupling of noise and QPOs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gamma-Ray Polarimetry of the Crab Pulsar Observed by POLAR

The X/$γ$ ray polarimetry of the Crab pulsar/nebula is believed to hold crucial information on their emission models. In the past, several missions have shown evidence of polarized emission from the Crab. The significance of these measurements remains however limited. New measurements are therefore required. POLAR is a wide Field of View Compton-scattering polarimeter (sensitive in 50-500 keV) onboard the Chinese spacelab Tiangong-2 which took data from September 2016 to April 2017. Although not designed to perform polarization measurements of pulsars, we present here a novel method which can be applied to POLAR as well as that of other wide Field of View polarimeters. The novel polarimetric joint-fitting method for the Crab pulsar observations with POLAR, allows us to obtain constraining measurements of the pulsar component. The best fitted values and corresponding 1$σ$ deviations for the averaged phase interval: (PD=$14\substack{+15 \\ -10}$\%, PA=$108\substack{+33 \\ -54} ^{\circ}$), for Peak 1: (PD=$17\substack{+18 \\ -12}$\%, PA=$174\substack{+39 \\ -36} ^{\circ}$) and for Peak 2: (PD=$16\substack{+16 \\ -11}$\%, PA=$78\substack{+39 \\ -30} ^{\circ}$). Further more, the 3$σ$ upper limits on the polarization degree are for the averaged phase interval (55\%), Peak 1 (66\%) and Peak 2 (57\%). Finally, to illustrate the capabilities of this method in the future, we simulated two years observation to the Crab pulsar with POLAR-2. The results show that POLAR-2 is able to confirm the emission to be polarized with $5σ$ and $4σ$ confidence level if the Crab pulsar is polarized at $20\%$ and $10\%$ respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

GRB 210121A: Observation of photospheric emissions from different regimes and the evolution of the outflow

GRB 210121A was observed by Insight-HXMT, Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM), Fermi Gamma-ray Burst Monitor (Fermi/GBM) on Jan 21st 2021. In this work, photospheric emission from a structured jet is preferred to interpret the prompt emission phase of GRB 210121A and emissions from different regimes are observed on-axis. Particularly, the emission from the intermediate photosphere is first observed in the first 1.3 s of the prompt emission, while those from the other part are dominant by the saturated regime, and offers an alternative explanation compared with the previous work. Moreover, the emissions with considering the intermediate photosphere can well interpret the changes on low-energy photon index $α$ during the pulses. Besides, the evolution of the outflow is extracted from time-resolved analysis, and a correlation of $Γ_0 \propto L^{0.25\pm0.05}_0$ is obtained, which implies that the jet may be mainly launched by neutrino annihilation in a hyper-accretion disk.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Insight-HXMT discovery of the highest energy CRSF from the first Galactic ultra-luminous X-ray pulsar Swift J0243.6+6124

The detection of cyclotron resonance scattering features (CRSFs) is the only way to directly and reliably measure the magnetic field near the surface of a neutron star (NS). The broad energy coverage and large collection area of \emph{Insight}-HXMT in the hard X-ray band allowed us to detect the CRSF with the highest energy known to date, reaching about 146 keV during the 2017 outburst of the first galactic pulsing ultraluminous X-ray source (pULX) Swift J0243.6+6124. During this outburst, the CRSF was only prominent close to the peak luminosity $\sim 2\times10^{39}$ erg s$^{-1}$, the highest to date in any of the Galactic pulsars. The CRSF is most significant in the spin phase region corresponding to the main pulse of the pulse profile, and its centroid energy evolves with phase from 120 to 146 keV. We identify this feature as the fundamental CRSF, since no spectral feature exists at $60-70$ keV. This is the first unambiguous detection of an electron CRSF from an ULX. We also estimate a surface magnetic field $\sim1.6\times10^{13}$ G for Swift J0243.6+6124. Considering that the dipole magnetic field strengths, inferred from several independent estimates of magnetosphere radius, are at least an order of magnitude lower than our measurement, we argue that the detection of the highest energy CRSF reported here unambiguously proves the presence of multipole field components close to the surface of the neutron star. Such a scenario has previously been suggested for several pulsating ULXs, including Swift J0243.6+6124, and our result represents the first direct confirmation of this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Insight-HXMT observation on 4U~1608--52: evidence of interplay between thermonuclear burst and accretion environments

A type-I burst could influence the accretion process through radiation pressure and Comptonization both for the accretion disk and the corona/boundary layer of an X-ray binary, and vice versa. We investigate the temporal evolution of a bright photospheric radius expansion (PRE) burst of 4U 1608-52 detected by Insight-HXMT in 1-50 keV, with the aim of studying the interplay between the burst and persistent emission. Apart from the emission from the neutron star (NS) surface, we find the residuals both in the soft (<3 keV) and hard (>10 keV) X-ray band. Time-resolved spectroscopy reveals that the excess can be attributed to an enhanced pre-burst/persistent emission or the Comptonization of the burst emission by the corona/boundary layer. The Comptonization model is a convolution thermal-Comptonization model (thcomp in XSPEC) and the Comptonization parameters are fixed at the values derived from the persistent emission. We find, during the PRE phase, after the enhanced pre-burst/persistent emission or the Comptonization of the burst emission is removed, the NS surface emission shows a plateau, and then a rise until the photosphere touches down to the NS surface, resulting in a flux peak at that moment. We speculate that the findings above correspond to that the obscured lower part of the NS surface by the disk is exposed to the line of sight due to the inner disk evaporation by the burst emission. The consistency between the fa model and convolution thermal-Comptonization model indicates the interplay between thermonuclear bursts and accretion environments. These phenomena did not usually show up by the conventional blackbody model fitting, which may be due to low count rate and narrow energy coverage in previous observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Phase dependent evolution within large luminosity range of 1A 0535+262 observed by Insight-HXMT during 2020 giant outburst

We have performed phase-resolved spectral analysis of the accreting pulsar 1A~0535+262 based on observations of Insight-HXMT during the 2020 Type-II outburst of the source. We focus on the two-dimensional dependence of the cyclotron resonance scattering features (CRSFs) along the outburst time and at different phases. The fundamental CRSF line (f-CRSF) shows different time- and phase-dependent behaviors. At higher luminosity, the phase profile of the f-CRSF energy changes from a single peak to double peaks, with the transition occurring at MJD 59185. On the contrary, the first harmonic CRSF (1-st CRSF) at $\sim$ 100 keV is only detected within a narrow phase range (0.8$-$1.0) accompanied by a shallow f-CRSF line. Based on these results, we speculate that when the source enters the supercritical regime, the higher accretion column can significantly enhance the harmonic line at a narrow phase through an "anti-pencil" beam at a higher energy band. At the same time, it will also affect the behavior of the fundamental line.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quasi-periodic oscillations of the X-ray burst from the magnetar SGR J1935+2154 and associated with the fast radio burst FRB 200428

The origin(s) and mechanism(s) of fast radio bursts (FRBs), which are short radio pulses from cosmological distances, have remained a major puzzle since their discovery. We report a strong Quasi-Periodic Oscillation(QPO) of 40 Hz in the X-ray burst from the magnetar SGR J1935+2154 and associated with FRB 200428, significantly detected with the Hard X-ray Modulation Telescope (Insight-HXMT) and also hinted by the Konus-Wind data. QPOs from magnetar bursts have only been rarely detected; our 3.4 sigma (p-value is 2.9e-4) detection of the QPO reported here reveals the strongest QPO signal observed from magnetars (except in some very rare giant flares), making this X-ray burst unique among magnetar bursts. The two X-ray spikes coinciding with the two FRB pulses are also among the peaks of the QPO. Our results suggest that at least some FRBs are related to strong oscillation processes of neutron stars. We also show that we may overestimate the significance of the QPO signal and underestimate the errors of QPO parameters if QPO exists only in a fraction of the time series of a X-ray burst which we use to calculate the Leahy-normalized periodogram.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of the corona in MAXI J1535-571 through type-C quasi-periodic oscillations with Insight-HXMT

Type-C quasi-periodic oscillations (QPOs) in black hole X-ray transients can appear when the source is in the low-hard and hard-intermediate states. The spectral-timing evolution of the type-C QPO in MAXI J1535-571 has been recently studied with Insight-HXMT. Here we fit simultaneously the time-averaged energy spectrum, using a relativistic reflection model, and the fractional rms and phase-lag spectra of the type-C QPOs, using a recently developed time-dependent Comptonization model when the source was in the intermediate state. We show, for the first time, that the time-dependent Comptonization model can successfully explain the X-ray data up to 100 keV. We find that in the hard-intermediate state the frequency of the type-C QPO decreases from 2.6 Hz to 2.1 Hz, then increases to 3.3 Hz, and finally increases to ~ 9 Hz. Simultaneously with this, the evolution of corona size and the feedback fraction (the fraction of photons up-scattered in the corona that return to the disc) indicates the change of the morphology of the corona. Comparing with contemporaneous radio observations, this evolution suggests a possible connection between the corona and the jet when the system is in the hard-intermediate state and about to transit into the soft-intermediate state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The First Insight-HXMT Gamma-Ray Burst Catalog: The First Four Years

The Hard X-ray Modulation Telescope (Insight-HXMT), is China's first X-ray astronomy satellite launched on June 15, 2017. The anti-coincidence CsI detectors of the High Energy X-ray telescope (HE) onboard Insight-HXMT could serve as an all-sky gamma-ray monitor in about 0.2-3 MeV. In its first four years of operation, Insight-HXMT has detected 322 Gamma-Ray Bursts (GRBs) by offline search pipeline including blind search and targeted search. For the GOLDEN sample of Insight-HXMT GRBs, joint analyses were performed with other GRB missions, including Fermi Gamma-ray Burst Monitor (Fermi/GBM), Swift Burst Alert Telescope (Swift/BAT) and Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM). It shows that Insight-HXMT can provide better constraint on GRB spectrum at higher energy band. The properties of Insight-HXMT GRBs are reported in detail, including their trigger time, duration, spectral parameters, peak fluxes of different time scales and fluence. This catalog is an official product of the Insight-HXMT GRB team.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

In-orbit timing calibration of the Insight-Hard X-ray Modulation Telescope

We describe the timing system and the timing calibration results of the three payloads on-board the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). These three payloads are the High Energy X-ray telescope (HE, 20-250 keV), the Medium Energy X-ray telescope (ME, 5-30 keV) and the low Energy X-ray telescope (LE, 1-10 keV). We present a method to correct the temperature-dependent period response and the long-term variation of the on-board crystal oscillator, especially for ME that does not carry a temperature-compensated crystal oscillator. The time of arrivals (ToAs) of the Crab pulsar are measured to evaluate the accuracy of the timing system. As the ephemeris of the Crab pulsar given by Jodrell Bank observatory has systematic errors around 40 μs (Rots et al. 2014), we use the quasi-simultaneous observations of the X-ray Timing Instrument (XTI) on-board the Neutron star Interior Composition Explorer (NICER) to produce the Crab ephemerides and to verify the timing system of Insight-HXMT. The energy-dependent ToAs' offsets relative to the NICER measurements including physical and instrumental origins are about 24.7μs, 10.1μs and 864.7μs, and the systematic errors of the timing system are determined as 12.1μs, 8.6μs, and 15.8μs, for HE, ME and LE respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Synchronous X-ray/Optical QPOs from the Black Hole LMXB MAXI J1820+070

We present high-speed optical photometry from SAAO and SALT on the black hole LMXB \MAXI (ASSASN-18ey), some of it simultaneous with NICER, Swift and Insight-HXMT X-ray coverage. We detect optical Quasi-Periodic Oscillations (QPOs) that move to higher frequencies as the outburst progresses, tracking both the frequency and evolution of similar X-ray QPOs previously reported. Correlated X-ray/optical data reveal a complex pattern of lags, including an anti-correlation and a sub-second lag that evolve over the first few weeks of outburst. They also show correlated components separated by a lag equal to the QPO period roughly centered on zero lag, implying that the inter-band variability is strongly and consistently affected by these QPOs at a constant phase lag of roughly +/- pi. The synchronisation of X-ray and optical QPOs indicates that they must be produced in regions physically very close to each other; we thus propose that they can be explained by a precessing jet model, based on analogies with V404 Cyg and MAXI J1348-630.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A peculiar cyclotron line near 16 keV detected in the 2015 outburst of 4U 0115+63?

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Background Model for the High-Energy Telescope of Insight-HXMT

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Background Model for the Low-Energy Telescope of Insight-HXMT

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calibration of the Instrumental Response of Insight-HXMT/HE CsI Detectors for Gamma-Ray Monitoring

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Design and Calibration of the High Energy Particle Monitor onboard the Insight-HXMT

Three high energy particle monitors (HPMs) employed onboard the Hard X-ray Modulation Telescope Insight-HXMT) can detect the charged particles from South Atlantic Anomaly (SAA) and hence provide the alert trigger for switch-on/off of the main detectors. Here a typical design of HPM with high stability and reliability is adopted by taking a plastic scintillator coupled with a small photomultiplier tube (PMT). The window threshold of HPM is designed as 1 MeV and 20 MeV for the incident electron and proton, respectively. Before the launch of Insight-HXMT, we performed in details the ground calibration of HPM. The measured energy response and its dependence on temperature are taken as essential input of Geant4 simulation for estimating the HPM count rate given with an incident particle energy spectrum. This serves as a guidance for choosing a reasonable working range of the PMT high voltage once the real SAA count rate is measured by HPM in orbit. So far the three HPMs have been working in orbit for more than two years. Apart from providing reliable alert trigger, the HPMs data are used as well to map the SAA region.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discovery of oscillations above 200 keV in a black hole X-ray binary with Insight-HXMT

Low-frequency quasi-periodic oscillations (LFQPOs) are commonly found in black hole X-ray binaries, and their origin is still under debate. The properties of LFQPOs at high energies (above 30 keV) are closely related to the nature of the accretion flow in the innermost regions, and thus play a crucial role in critically testing various theoretical models. The Hard X-ray Modulation Telescope (Insight-HXMT) is capable of detecting emissions above 30 keV, and is therefore an ideal instrument to do so. Here we report the discovery of LFQPOs above 200 keV in the new black hole MAXI J1820+070 in the X-ray hard state, which allows us to understand the behaviours of LFQPOs at hundreds of kiloelectronvolts. The phase lag of the LFQPO is constant around zero below 30 keV, and becomes a soft lag (that is, the high-energy photons arrive first) above 30 keV. The soft lag gradually increases with energy and reaches ~0.9s in the 150-200 keV band. The detection at energies above 200 keV, the large soft lag and the energy-related behaviors of the LFQPO pose a great challenge for most currently existing models, but suggest that the LFQPO probably originates from the precession of a small-scale jet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

eXTP perspectives for the $ν$MSM sterile neutrino dark matter model

We discuss the potential of the eXTP X-ray telescope, in particular its Spectroscopic Focusing Array (SFA), Large Area Detector (LAD) and Wide Field Monitor (WFM) for the detection of a signal from keV-scale decaying dark matter. We show that the sensitivity of the eXTP is sufficient to improve existing constraints on the mixing angle of the neutrino Minimal extension of the Standard Model ($ν$MSM) by a factor of 5-10 within the dark matter mass range 2-50 keV, assuming 1% level of systematic uncertainty. We assert that the eXTP will be able to probe previously inaccessible range of $ν$MSM parameters and serve as a precursor for the Athena mission in decaying dark matter searches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

In-flight calibration of the Insight-Hard X-ray Modulation Telescope

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

In-orbit Calibration to the Point-Spread Function of Insight-HXMT

We make the in-orbit calibration to the point-spread functions (PSFs) of the collimators of the Hard X-ray Modulation Telescope with the scanning observation of the Crab. We construct the empirical adjustments to the theoretically calculated geometrical PSFs. The adjustments contain two parts: a rotating matrix to adjust the directional deviation of the collimators and a paraboloidal function to correct the inhomogeneity of the real PSFs. The parameters of the adjusting matrices and paraboloidal functions are determined by fitting the scanning data with lower scanning speed and smaller intervals during the calibration observations. After the PSF calibration, the systematic errors in source localization in the Galactic plane scanning survey are 0.010 deg, 0.015 deg, 0.113 deg for the Low-Energy Telescope (LE), the Medium-Energy telescope (ME) and the High-Energy telescope (HE), respectively; meanwhile, the systematic errors in source flux estimation are 1.8%, 1.6%, 2.7% for LE, ME and HE, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Introducing a special collection of papers in the Journal of High Energy Astrophysics on the Early Results of China's 1st X-ray Astronomy Satellite: Insight-HXMT

Insight-HXMT is the first Chinese X-ray astronomical mission, launched successfully on June 15, 2017, from China's Jiuquan Satellite Launch Center. Insight-HXMT was designed to have a broad energy coverage in X-rays, from 1-250 keV, with excellent timing and adequate energy resolution at soft X-rays, and the largest effective area at hard X-rays. Here, we present a collection of papers of the Journal of High Energy Astrophysics on the Early Results of China's 1st X-ray Astronomy Satellite Insight-HXMT. These papers cover the in-orbit performance, the background model, and all calibration results, together with several first results on observations of binaries and details on the Galactic plane survey.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Background Model of the Medium Energy X-ray telescope of Insight-HXMT

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The POLAR Gamma-Ray Burst Polarization Catalog

Despite over 50 years of research, many open questions remain about the origin and nature of GRBs. Polarization measurements of the prompt emission of these extreme phenomena have long been thought to be the key to answering a range of these questions. The POLAR detector was designed to produce the first set of detailed and reliable polarization measurements in an energy range of approximately 50-500 keV. During late 2016 and early 2017, POLAR detected a total of 55 GRBs. Analysis results of 5 of these GRBs have been reported in the past. The results were found to be consistent with a low or unpolarized flux. However, previous reports by other collaborations found high levels of polarization. We study the polarization for all the 14 GRBs observed by POLAR for which statistically robust inferences are possible. Additionally, time-resolved polarization studies are performed on GRBs with sufficient apparent flux. A publicly available polarization analysis tool, developed within the 3ML framework, was used to produce statistically robust results. The method allows to combine spectral and polarimetric data from POLAR with spectral data from the Fermi GBM and Swift-BAT to jointly model the spectral and polarimetric parameters. The time integrated analysis finds all results to be compatible with a low or zero polarization with the caveat that, when time-resolved analysis is possible within individual pulses, we observe moderate polarization with a rapidly changing polarization angle. Thus, time-integrated polarization results, while pointing to lower polarization are potentially an artifact of summing over the changing polarization signal and thus, washing out the true moderate polarization. Therefore, we caution against over interpretation of any time-integrated results and encourage one to wait for more detailed polarization measurements from forthcoming missions such as POLAR-2 and LEAP.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Exact solutions for spherically gravitational collapse around a black hole: the effect of tangential pressure

Spherically gravitational collapse towards a black hole with non-zero tangential pressure is studied. Exact solutions corresponding to different equations of state are given. We find that when taking the tangential pressure into account, the exact solutions have three qualitatively different endings. For positive tangential pressure, the shell around a black hole may eventually collapse onto the black hole, or expand to infinity, or have a static but unstable solution, depending on the combination of black hole mass, mass of the shell and the pressure parameter. For vanishing or negative pressure, the shell will collapse onto the black hole. For all eventually collapsing solutions, the shell will cross the event horizon, instead of accumulating outside the event horizon, even if clocked by a distant stationary observer.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detection of a possible X-ray Quasi-periodic Oscillation in the Active Galactic Nucleus 1H~0707-495

Quasi-periodic oscillation (QPO) detected in the X-ray radiation of black hole X-ray binaries (BHXBs) is thought to originate from dynamical processes in the close vicinity of the black holes (BHs), and thus carries important physical information therein. Such a feature is extremely rare in active galactic nuclei (AGNs) with supermassive BHs. Here we report on the detection of a possible X-ray QPO signal with a period of 3800\,s at a confidence level $>99.99\%$ in the narrow-line Seyfert 1 galaxy (NLS1) 1H~0707-495 in one data set in 0.2-10\,keV taken with {\it XMM-Newton}. The statistical significance is higher than that of most previously reported QPOs in AGNs. The QPO is highly coherent (quality factor $Q=ν/Δν\geqslant 15$) with a high rms fractional variability ($\sim15\%$). A comprehensive analysis of the optical spectra of this AGN is also performed, yielding a central BH mass $5.2\times10^6\,M_{\odot}$ from the broad emission lines based on the scaling relation. The QPO follows closely the known frequency-BH mass relation, which spans from stellar-mass to supermassive BHs. The absence of the QPO in other observations of the object suggests it a transient phenomenon. We suggest that the (high-frequency) QPOs tend to occur in highly accreting BH systems, from BHXBs to supermassive BHs. Future precise estimation of the BH mass may be used to infer the BH spin from the QPO frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Do we expect to detect electromagnetic radiation from merging stellar mass black binaries like GW150914? No

Context: The LIGO consortium announced the first direct detection of gravitation wave event GW150914 from two merging black holes; however the nature of the black holes are still not clear. Aims: We study whether electromagnetic radiation can be detected from merging stellar mass black binaries like GW150914. Methods: We briefly investigate the possible growth and merging processes of the two stellar mass black holes in the merging event of GW150914 detected by aLIGO, as clocked by a distant external observer. Our main results are: (1) The description of the black hole growth using stationary metric of a pre-existing black hole predicts strong electromagnetic radiation from merging black holes, which is inconsistent with GW150914; (2) Only gravitational wave radiation can be produced in the coalescence of two black holes such as that in the GW150914 event, if the black hole growth is described using time-dependent metric considering the influence of the in-falling matter onto a pre-existing black hole, as clocked by a distant external observer. Conclusions: Future high sensitivity detections of gravitational waves from merging black holes might be used to probe matter distribution and space-time geometry in the vicinity of the horizon. Perhaps the GW150914-like events can be identified with traditional astronomy observations only if the black holes are embedded in extremely dense medium before their final merge, when very strong electromagnetic radiation is produced and can escape from the system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of very hard electron and gamma-ray spectra of flat spectrum radio quasar in fast-cooling regime

In external Compton scenario, we investigate the formation of the very hard electron spectrum in the fast-cooling regime, using a time-dependent emission model. It is shown that a very hard electron distribution $N'_{\rm e}(\gp)\propto\gp^{-p}$ with the spectral index $p\sim1.3$ is formed below the minimum energy of injection electron when inverse Compton scattering takes place in the Klein-Nishina regime, i.e., inverse Compton scattering of relativistic electrons on broad-line region radiation in flat spectrum radio quasars. This produces a very hard gamma-ray spectrum, and can reasonably explain the very hard \emph{Fermi}-LAT spectrum of the flat spectrum radio quasar 3C 279 during the extreme gamma-ray flare in 2013 December.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Perspective of monochromatic gamma-ray line detection with the High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $γ$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $γ$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Similarity of Jet Radiation between Flat Spectrum Radio Quasars and GeV Narrow-Line Seyfert 1 Galaxies: A Universal $δ$-$L_{\rm c}$ Correlation

By modeling the broadband spectral energy distributions (SEDs) of a typical flat spectrum radio quasar (FSRQ, 3C 279) and two GeV narrow-line Seyfert 1 galaxies (NLS1s, PMN J0948+0022 and 1H 0323+342) in different flux stages with the one-zone leptonic models, we find a universal correlation between their Doppler factors ($δ$) and peak luminosities ($L_{\rm c}$) of external Compton scattering bumps. Compiling a combined sample of FSRQs and GeV NLS1s, it is found that both FSRQs and GeV NLS1s in different stages and in different sources well follow the same $δ$-$L_{\rm c}$ correlation. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. And the universal $δ$-$L_{\rm c}$ relation between FSRQs and GeV NLS1s in different stages may be further evidence that the particle acceleration and radiation mechanisms for the two kinds of sources are similar. In addition, by replacing $L_{\rm c}$ with the observed luminosity in the Fermi/LAT band ($L_{\rm LAT}$), this correlation holds, and it may serve as an empirical indicator of $δ$. We estimate the $δ$ values with $L_{\rm LAT}$ for 484 FSRQs in the Fermi/LAT Catalog and they range from 3 to 41, with a median of 16, which are statistically consistent with the values derived by other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Study of linearity of LYSO crystal for HERD calorimeter

The High Energy cosmic Radiation Detection (HERD) facility is a space mission designed for detecting cosmic ray (CR) electrons, $γ$-rays up to tens of TeV and CR nuclei from proton to iron up to several PeV. The main instrument of HERD is a 3-D imaging calorimeter (CALO) composed of nearly ten thousand cubic LYSO crystals. A large dynamic range of single HERD CALO Cell (HCC) is necessary to achieve HERD's PeV observation objectives, which means that the response of HCC should maintain a good linearity from minimum ionizing particle (MIP) calibration to PeV shower maximum. In order to study the linearity of HCC over such a large energy range, a beam test has been implemented at the E2 and E3 beam lines of BEPC. High intensity pulsed electron beam provided by E2 line are used for producing high energy density within HCC; $π^{+}$/proton provided by E3 line are used for HCC calibration. The results show that no saturation effect occurs and the linearity of HCC is better than 10% from 30 MeV (1 MIP) to 1.1$\times$10$^{3}$ TeV (energy density is 93 TeV/cm$^{3}$), which can meet the requirement mentioned above.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing Einstein's Equivalence Principle with Cosmological Fast Radio Bursts behind Clusters of Galaxies

Aims: Recently, cosmological fast radio bursts (FRBs) have been used to provide the most stringent limit up to date on Einstein's Equivalence Principle (EEP). We study how to further test EEP with FRBs. Methods: Future systematic radio surveys will certainly find abundant FRBs at cosmological distances and some of them will inevitably be located behind clusters of galaxies. Here we suggest to use those FRBs to further test EEP. Results: We find that the robustness and accuracy of testing EEP can be improved further by orders of magnitude with these FRBs. The same methodology can also be applied to any other types of fast and bright transients at cosmological distances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing Einstein's Weak Equivalence Principle With Gravitational Waves

A conservative constraint on the Einstein Weak Equivalence Principle (WEP) can be obtained under the assumption that the observed time delay between correlated particles from astronomical sources is dominated by the gravitational fields through which they move. Current limits on the WEP are mainly based on the observed time delays of photons with different energies. It is highly desirable to develop more accurate tests that include the gravitational wave (GW) sector. The detection by the advanced LIGO/VIRGO systems of gravitational waves will provide attractive candidates for constraining the WEP, extending the tests to gravitational interactions, with potentially higher accuracy. Considering the capabilities of the advanced LIGO/VIRGO network and the source direction uncertainty, we show that the joint detection of GWs and electromagnetic signals could probe the WEP to an accuracy down to $10^{-10}$, which is one order of magnitude tighter than previous limits, and seven orders of magnitude tighter than the multi-messenger (photons and neutrinos) results by supernova 1987A.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing one-zone synchrotron-self-Compton models with spectral energy distributions of Mrk 421

We test one-zone synchrotron self-Compton (SSC) models with high-quality multiwavelength spectral energy distribution (SED) data of Mrk 421. We use Markov chain Monte Carlo (MCMC) technique to fit twelve day-scale SEDs of Mrk 421 with one-zone SSC models. Three types of electron energy distribution (EED), a log-parabola (LP) EED, a power-law log-parabola (PLLP) EED and a broken power-law (BPL) EED, are assumed in fits. We find that the one-zone SSC model with the PLLP EED provides successful fits to all the twelve SEDs. However, the one-zone SSC model with the LP and BPL EEDs fail to provide acceptable fits to the highest energy X-ray data or GeV data in several states. We therefore conclude that the one-zone SSC model works well in explaining the SEDs of Mrk 421, and the PLLP EED is preferred over the LP and BPL EEDs for Mrk 421 during the flare in March 2010. We derive magnetic field $B'\sim0.01$\ G, Doppler factor $δ_{\rm D}\sim$30--50, and the curvature parameter of EED $r\sim1$--$10$ in the model with the PLLP EED. The evolutions of model parameters are explored. The physical implications of our results are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The MZ relation for local star-forming galaxies

We investigate the evolution of the mass-metallicity (MZ) relation with a large sample of 53,444 star-forming galaxies (SFGs) at $0.04<z<0.12$, selected from the catalog of MPA-JHU emission-line measurements for the Sloan Digital Sky Survey (SDSS) Data Release 7. Regarding the sample of SFGs, we correct the observational bias and raise the aperture covering fractions to check the reliability of the metallicity evolution. We show that (1) the redshift evolution of log($\rm Hα$) and log(O III) luminosities is displayed in our sample; (2) we find the metallicity evolution of $\sim 0.15$ dex at $\rm log (M_*/M_\odot)\sim9.3$ in SFGs at $0.04<z<0.12$; (3) after applying the luminosity thresholds of log$(L_{\rm H α})>41.0$ and log$(L_{\rm O\ III})>39.7$, we find that metallicity evolution is shown well, and that SFR evolution still is shown well under the latter luminosity threshold, but the evolution is not observed under the former one; (4) the evolution of the MZ relation seems to disappear at about $\rm log(M_{*}/M_\odot)>10.0$ after applying the luminosity threshold of log$(L_{\rm H α})>41.0$ or log$(L_{\rm O\ III})>39.7$; (5) we find $α=0.09$ and $α=0.07$ in the equation ($μ={\rm log}M_{*}-α\rm log(SFR)$) for log$(L_{\rm H α})>41.0$ and log$(L_{\rm O\ III})>39.7$ samples, respectively, and these imply that the evolution of the MZ relation may have a weaker dependence on SFR in our sample.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Understanding the residual patterns of timing solutions of radio pulsars with a model of magnetic field oscillation

We explain some phenomena existing generally in the timing residuals: amplitude and sign of the second derivative of a pulsar's spin-frequency ($\ddotν$), some sophisticated residual patterns, which also change with the time span of data segments. The sample is taken from Hobbs et al.\,(2010), in which the pulsar's spin-frequency and its first derivative have been subtracted from the timing solution fitting. We first classify the timing residual patterns into different types based on the sign of $\ddotν$. Then we use the magnetic field oscillation model developed in our group \citep{zhang12a} to fit successfully the different kinds of timing residuals with the Markov Chain Monte Carlo method. Finally, we simulate the spin evolution over 20 years for a pulsar with typical parameters and analyze the data with the conventional timing solution fitting. By choosing different segments of the simulated data, we find that most of the observed residual patterns can be reproduced successfully. This is the first time that the observed residual patterns are fitted by a model and reproduced by simulations with very few parameters. From the distribution of the different residual patterns in the $P-\dot P$ diagram, we argue that (1) a single magnetic field oscillation mode exists commonly in all pulsars throughout their lifetimes; (2) there may be a transition period over the lifetimes of pulsars, in which multiple magnetic field oscillation modes exist.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A self-consistent interpretation of the GeV-TeV emission from a distant blazar PKS 1424+240

We propose a scenario for self-consistent interpretation of GeV - TeV spectrum of a distant blazar PKS 1424+240. In this scenario, ultra-high energy (UHE) protons are assumed to exist in the blazar jet and produce gamma rays through synchrotron emission emitted by relativistic protons and pair cascades (resulted from $pγ$ interaction); meanwhile, some of these UHE protons may escape from the jet and are injected into intergalactic space. Therefore, we assume that UHE cosmic rays (CR) originate from relativistic protons in the jet and use energy-independent escape timescale to obtain UHECR injection spectrum. Both contributions of gamma rays injected by the source and secondary gamma rays produced in interactions of UHECRs emitted by the blazar with photon background during their propagation through intergalactic space are calculated. Our results show that this scenario is able to reproduce the GeV-TeV spectrum of PKS 1424+240 self-consistently in a broad range of redshifts $0.6<z<1.3$. The required relativistic jet power $L_p\simeq3\times10^{46}\ \rm erg\ s^{-1}$ only moderately depends on the assumed source redshift, while the proton escape timescale (and the injected UHECR luminosity) strongly depends on $z$. We compare the integral TeV fluxes predicted in our scenario with the sensitivity of the planned Cherenkov Telescope Array (CTA), and discuss the implications to future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Clumpy tori illuminated by the anisotropic radiation of accretion discs in active galactic nuclei

In this paper, we try to explain the observed correlation between the covering factor (CF) of hot dust and the properties of active galactic nuclei (AGNs), e.g., the bolometric luminosity ($L_{\rm{bol}}$) and black hole mass ($M_{\rm{BH}}$). Combining the possible dust distribution in the torus, the angular dependence of the radiation of the accretion disc, and the relation between the critical angle of torus and the Eddington ratio, there are eight possible models investigated in our work. We fit the observed CF with these models to determine the parameters of them. As a result, clumpy torus models can generally explain the observed correlations of tori, while the smooth models fail to produce the required CFs. However, there is still significant scatter even for the best-fitting model, which is the combination of a clumpy torus illuminated by the anisotropic radiation of accretion disc in an AGN. Although some of the observed scatter is due to the uncertainties in measuring $L_{\rm{bol}}$ and $M_{\rm{BH}}$, other factors are required in more realistic model. The models examined in this paper are not necessary to be the physical model of tori. However, the reasonable assumptions selected during this process should be helpful in constructing physical models of tori.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Correlation Analysis between Tibet AS-$γ$ TeV Cosmic Ray and WMAP Nine-year Data

The Wilkinson Microwave Anisotropy Probe (WMAP) team subtracted template-based foreground models to produce foreground-reduced maps, and masked point sources and uncertain sky regions directly; however, whether foreground residuals exist in the WMAP foreground-reduced maps is still an open question. Here, we use Pearson correlation coefficient (PCC) analysis with AS-$γ$ TeV cosmic ray (CR) data to probe possible foreground residuals in the WMAP nine-year data. The correlation results between the CR and foreground-contained maps (WMAP foreground-unreduced maps, WMAP template-based and MEM foreground models) suggest that: (1) CRs can trace foregrounds in the WMAP data; (2) at least some TeV CRs originate from the Milky Way; (3) foregrounds may be related to the existence of CR anisotropy (loss-cone and tail-in structures); (4) there exist differences among different types of foregrounds in the declination range of $< 15^{\circ}$. Then, we generate 10,000 mock CMB sky maps to describe the cosmic variance, which is used to measure the effect of the fluctuations of all possible CMB maps to the correlations between CR and CMB maps. Finally, we do correlation analysis between the CR and WMAP foreground-reduced maps, and find that: (1) there are significant anticorrelations; and (2) the WMAP foreground-reduced maps are credible. However, the significant anticorrelations may be accidental, and the higher signal-to-noise ratio (SNR) Planck SMICA map cannot reject the hypothesis of accidental correlations. We therefore can only conclude that the foreground residuals exist with $\sim$ 95\% probability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Correlations among Jet, Accretion Disk, and Broad Line Region of Flat Spectrum Radio Quasars

The SEDs of 18 GeV FSRQs are collected and compiled from literature, in which both the jet emission and the accretion disk radiation can be observed, in order to investigate the correlations among their jet power (P_jet), accretion disk luminosity (L_disk), and luminosity of broad line region (BLR, L_BLR). On the basis of the SED fits with the jet radiation and accretion disk radiation models, we calculate P_jet and L_disk. No correlation between P_jet with either L_disk or L_BLR is found. With a sub-sample of L_BLR for 13 GeV FSRQs, it is observed that L_BLR is strongly correlated with their L_disk. We also study the BLR covering factors of the GeV FSRQs in our sample, averagely which are smaller than that of the large samples of radio-loud and radio-quiet quasars. P_jet of some GeV FSRQs is higher than L_disk, but P_jet of all the GeV FSRQs is lower than the accretion power of black hole (BH), which is estimated by \dot{M}c^2=L_disk/0.1, indicating that the total accretion power of BH is sufficient to drive the jets in these sources; however the uncorrelation between L_disk and P_jet of the GeV FSRQs may suggest that their jets are launched by the BZ process via extracting the rotational energy of BH. Using the L_BLR-L_disk relation of the GeV FSRQs, we estimate L_disk of a BL Lac sample with their L_BLR. A comparison of L_BLR and Eddington ratio (L_disk/L_Edd) among BL Lacs, very radio-loud NLS1 galaxies, and FSRQs is also presented. It is found that along with the BL Lac-NLS1-FSRQ sequence L_BLR and L_disk/L_Edd increase, which may correspond to the change of the accretion disk structure and the transformation of the dominant mechanism for jet launching. This is also consistent with the division of their parent populations, i.e., low-excitation radio galaxies and high-excitation radio galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diagnosing the burst influence upon accretion in the clocked burster GS 1826--238

Type-I X-ray bursts on the surface of a neutron star are a unique probe to the accretion in X-ray binary systems. However, we know little about the feedback of the burst emission upon accretion. Hard X-ray shortages and enhancements of the persistent emission at soft X-rays have been observed. To put these findings in context with the aim of understanding the possible mechanism underneath, we investigated 68 bursts seen by RXTE from the clocked burster GS 1826--238. We diagnosed jointly the burst influence at both soft and hard X-rays, and found that the observations can be described as the CompTT model with variable normalization, electron temperature and optical depth. Putting these results in a scenario of coronal Compton cooling via the burst emission would lead to a shortage of the cooling power, which may suggest that additional consideration like the influence of the burst on the corona formation should be accounted for as well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamic changes of emitting electron distribution in the jet of 3C 279: signatures of acceleration and cooling

We study the dynamic changes of electron energy distribution (EED) through systematically analysing the quasi-simultaneous spectral energy distributions (SEDs) of the flat spectrum radio quasar 3C 279 in different states. With Markov chain Monte Carlo (MCMC) technique we model fourteen SEDs of 3C 279 using a leptonic model with a three-parameter log-parabola electron energy distribution (EED). The 14 SEDs can be satisfactorily fitted with the one-zone leptonic model. The observed $γ$ rays in 13 states are attributed to Compton scattering of external infrared photons from a surrounding dusty torus. The curved $γ$-ray spectrum observed during 2-8 April 2014 is well explained by the external Compton of dust radiation. It is found that there is a clear positive correlation between the curvature parameter $b$ of the EED and the electron peak energy $γ'_{\rm pk}$. No clear correlation between $b$ and the synchrotron peak frequency $ν_{\rm s}$ is found, due to the varied product of Doppler factor and fluid magnetic field from state to state. We interpret the correlation of $b-γ'_{\rm pk}$ in a stochastic acceleration scenario. This positive correlation is in agreement with the prediction in the stage when the balance between acceleration and radiative cooling of the electrons is nearly established in the case of the turbulence spectral index $q=2$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy alignment as a probe of large-scale filaments

The orientations of the red galaxies in a filament are aligned with the orientation of the filament. We thus develop a location-alignment-method (LAM) of detecting filaments around clusters of galaxies, which uses both the alignments of red galaxies and their distributions in two-dimensional images. For the first time, the orientations of red galaxies are used as probes of filaments. We apply LAM to the environment of Coma cluster, and find four filaments (two filaments are located in sheets) in two selected regions, which are compared with the filaments detected with the method of \cite{Falco14}. We find that LAM can effectively detect the filaments around a cluster, even with $3σ$ confidence level, and clearly reveal the number and overall orientations of the detected filaments. LAM is independent of the redshifts of galaxies, and thus can be applied at relatively high redshifts and to the samples of red galaxies without the information of redshifts. We also find that the images of background galaxies (interlopers) which are lensed by the gravity of foreground filaments are amplifiers to probe the filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling the evolution and distribution of the frequency's second derivative and the braking index of pulsar spin

We model the evolution of the spin frequency's second derivative $\ddotν$ and the braking index $n$ of radio pulsars with simulations within the phenomenological model of their surface magnetic field evolution, which contains a long-term power-law decay modulated by short-term oscillations. For the pulsar PSR B0329+54, a model with three oscillation components can reproduce its $\ddotν$ variation. We show that the "averaged" $n$ is different from the instantaneous $n$, and its oscillation magnitude decreases abruptly as the time span increases, due to the "averaging" effect. The simulated timing residuals agree with the main features of the reported data. Our model predicts that the averaged $\ddotν$ of PSR B0329+54 will start to decrease rapidly with newer data beyond those used in Hobbs et al.. We further perform Monte Carlo simulations for the distribution of the reported data in $|\ddotν|$ and $|n|$ versus characteristic age $τ_{\rm c}$ diagrams. It is found that the magnetic field oscillation model with decay index $α=0$ can reproduce the distributions quite well. Compared with magnetic field decay due to the ambipolar diffusion ($α=0.5$) and the Hall cascade ($α=1.0$), the model with no long term decay ($α=0$) is clearly preferred for old pulsars by the p-values of the two-dimensional Kolmogorov-Smirnov test.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-wavelength light curve evolution of Swift J1357.2-0933 during its 2011 outburst

Swift J1357.2-0933 underwent an episodic accretion in 2011 and provided very regular temporal and spectral evolution, making it an ideal source for exploring the nature of very faint X-ray transients (VFXTs). In this work, we present a detailed analysis on both X-ray and near-ultraviolet (NUV) light curves. The fluxes at all wavelengths display a near-exponential decays in the early phase and transits to a faster-decay at late times. The e-folding decay time-scales monotonically decrease with photon energies, and the derived viscous time-scale is $τ_{\rm \dot{M}} \sim 60$ days. The time-scale in the late faster-decay stage is about a few days. The high ratio of NUV luminosity to X-ray luminosity indicates that the irradiation is unimportant in this outburst, while the near-exponential decay profile and the long decay time-scales conflict with the disc thermal-viscous instability model. We thus suggest that the disc is thermally stable during the observations. Adopting the truncated disc model, we obtain a lower limit of peak accretion rate of $0.03 \dot{M}_{\rm Edd}$ and the X-ray radiative efficiency $η< 5\times10^{-4}$, which decreases as the luminosity declines. The low X-ray radiative efficiency is caused by the combined action of advection and outflows, and naturally explains that the X-ray reprocessing is overwhelmed by the viscous radiation of the outer standard disc in the NUV regime. We also propose a possibility that the outer standard disc recedes from the central black hole, resulting in the faster-decay at late times.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parameter constraints in a near-equipartition model with multi-frequency \emph{NuSTAR}, \emph{Swift} and \emph{Fermi-LAT} data from 3C 279

Precise spectra of 3C 279 in the 0.5-70 keV range, obtained during two epochs of \emph{Swift} and \emph{NuSTAR} observations, are analyzed using a near-equipartition model. We apply a one-zone leptonic model with a three-parameter log-parabola electron energy distribution (EED) to fit the \emph{Swift} and \emph{NuSTAR} X-ray data, as well as simultaneous optical and \emph{Fermi}-LAT $γ$-ray data. The Markov Chain Monte Carlo (MCMC) technique is used to search the high-dimensional parameter space and evaluate the uncertainties on model parameters. We show that the two spectra can be successfully fit in near-equipartition conditions, defined by the ratio of the energy density of relativistic electrons to magnetic field $ζ_{\rm e}$ being close to unity. In both spectra, the observed X-rays are dominated by synchrotron-self Compton photons, and the observed $γ$ rays are dominated by Compton scattering of external infrared photons from a surrounding dusty torus. Model parameters are well constrained. From the low state to the high state, both the curvature of the log-parabola width parameter and the synchrotron peak frequency significantly increase. The derived magnetic fields in the two states are nearly identical ($\sim1$\ G), but the Doppler factor in the high state is larger than that in the low state ($\sim$28 versus $\sim$18). We derive that the gamma-ray emission site takes place outside the broad-line region, at $\gtrsim$ 0.1 pc from the black hole, but within the dusty torus. Implications for 3C 279 as a source of high-energy cosmic-rays are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Primordial alignment of elliptical galaxies in intermediate redshift clusters

We measure primordial alignments for the red galaxies in the sample of eight massive galaxy clusters in the southern sky from the CLASH-VLT Large Programme, at a median redshift of 0.375. We find primordial alignment with about $3σ$ significance in the four dynamically young clusters, but null detection of primordial alignment in the four highly relaxed clusters. The observed primordial alignment is not dominated by any single one of the four dynamically young clusters, and is primarily due to a population of bright galaxies ($M_r<-20.5\ \rm{m}$) residing in the region 300 to 810 kpc from the cluster centers. For the first time, we point out that the combination of radial alignment and halo alignment can cause fake primordial alignment. Finally, we find that the detected alignment for the dynamically young clusters is real rather than fake primordial alignment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radial alignment of elliptical galaxies by the tidal force of a cluster of galaxies

Unlike the random radial orientation distribution of field elliptical galaxies, galaxies in a cluster are expected to point preferentially towards the center of the cluster, as a result of the cluster's tidal force on its member galaxies. In this work an analytic model is formulated to simulate this effect. The deformation time scale of a galaxy in a cluster is usually much shorter than the time scale of change of the tidal force; the dynamical process of the tidal interaction within the galaxy can thus be ignored. An equilibrium shape of a galaxy is then assumed to be the surface of equipotential, which is the sum of the self-gravitational potential of the galaxy and the tidal potential of the cluster at this location. We use a Monte-Carlo method to calculate the radial orientation distribution of these galaxies, by assuming the NFW mass profile of the cluster and the initial ellipticity of field galaxies. The radial angles show a single peak distribution centered at zero. The Monte-Carlo simulations also show that a shift of the reference center from the real cluster center weakens the anisotropy of the radial angle distribution. Therefore, the expected radial alignment cannot be revealed if the distribution of spatial position angle is used instead of that of radial angle. The observed radial orientations of elliptical galaxies in cluster Abell~2744 are consistent with the simulated distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super-Eddington Accretion in the Ultraluminous X-ray Source NGC1313 X-2: An Ephemeral Feast

We investigate the X-ray spectrum, variability and the surrounding ionized bubble of NGC1313 X-2 to explore the physics of super-Eddington accretion. Beyond the Eddington luminosity, the accretion disk of NGC1313 X-2 is truncated at a large radius ($\sim$ 50 times of innermost stable circular orbit), and displays the similar evolution track with both luminous Galactic black-hole and neutron star X-ray binaries. In super-critical accretion, the speed of radiatively driven outflows from the inner disk is mildly relativistic. Such ultra-fast outflows would be over ionized and might produce weak Fe K absorption lines, which may be detected by the coming X-ray mission {\it Astro-H}. If the NGC1313 X-2 is a massive stellar X-ray binary, the high luminosity indicates that an ephemeral feast is held in the source. That is, the source must be accreting at a hyper-Eddington mass rate to give the super-Eddington emission over $\sim 10^{4}-10^{5}$ yr. The expansion of the surrounding bubble nebula with a velocity of $\sim$ 100 km~s$^{-1}$ might indicate that it has existed over $\sim 10^{6}$ yr and is inflated by the radiatively driven outflows from the transient with a duty cycle of activity of $\sim$ a few percent. Alternatively, if the surrounding bubble nebula is produced by line-driven winds, less energy is required than the radiatively driven outflow scenario, and the radius of the Strömgren radius agrees with the nebula size. Our results are in favor of the line-driven winds scenario, which can avoid the conflict between the short accretion age and the apparently much longer bubble age inferred from the expansion velocity in the nebula.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super-Strong Magnetic Fields of Neutron Stars in Be X-Ray Binaries Estimated with New Torque and Magnetosphere Models

We re-estimate the surface magnetic fields of neutron stars (NSs) in Be X-ray binaries (BeXBs) with different models of torque, improved beyond Klus et al. (2014). In particular, a new torque model (Dai \& Li 2006) is applied to three models of magnetosphere radius. Unlike the previous models, the new torque model does not lead to divergent results for any fastness parameter. The inferred surface magnetic fields of these NSs for the two compressed magnetosphere models are much higher than that for the uncompressed magnetosphere model. The new torque model using the compressed-magnetosphere radius (Shi, Zhang \& Li 2014) leads to unique solutions near spin equilibrium in all cases, unlike other models that usually give two branches of solutions. Although our conclusions are still affected by the simplistic assumptions about the magnetosphere radius calculations, we show several groups of possible surface magnetic field values with our new models when the interaction between the magnetosphere and the infalling accretion plasma is considered. The estimated surface magnetic fields for NSs BeXBs in the Large Magellanic Cloud, the Small Magellanic Cloud and the Milk Way are between the quantum critical field and the maximum "virial" value by the spin equilibrium condition.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The consistency test on the cosmic evolution

We propose a new and robust method to test the consistency of the cosmic evolution given by a cosmological model. It is realized by comparing the combined quantity r_d^CMB/D_V^SN, which is derived from the comoving sound horizon r_d from cosmic microwave background (CMB) measurements and the effective distance D_V derived from low-redshift Type-Ia supernovae (SNe Ia) data, with direct and independent r_d/D_V obtained by baryon acoustic oscillation (BAO) measurements at median redshifts. We apply this test method for the Lambda-CDM and wCDM models, and investigate the consistency of the derived value of r_d/D_V from Planck 2015 and the SN Ia data sets of Union2.1 and JLA (z<1.5), and the r_d/D_V directly given by BAO data from six-degree-field galaxy survey (6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 of SDSS-III, WiggleZ and Ly-alpha forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) DR-11 over 0.1<z<2.36. We find that r_d^CMB/D_V^SN for both non-flat Lambda-CDM and flat wCDM models with Union2.1 and JLA data are well consistent with the BAO and CMB measurements within 1-sigma CL. Future surveys will further tight up the constraints significantly, and provide stronger test on the consistency.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The evolution of the magnetic inclination angle as an explanation of the long term red timing-noise of pulsars

We study the possibility that the long term red timing-noise in pulsars originates from the evolution of the magnetic inclination angle $χ$. The braking torque under consideration is a combination of the dipole radiation and the current loss. We find that the evolution of $χ$ can give rise to extra cubic and fourth-order polynomial terms in the timing residuals. These two terms are determined by the efficiency of the dipole radiation, the relative electric-current density in the pulsar tube and $χ$. The following observation facts can be explained with this model: a) young pulsars have positive $\ddotν$; b) old pulsars can have both positive and negative $\ddotν$; c) the absolute values of $\ddotν$ are proportional to $-\dotν$; d) the absolute values of the braking indices are proportional to the characteristic ages of pulsars. If the evolution of $χ$ is purely due to rotation kinematics, then it can not explain the pulsars with braking index less than 3, and thus the intrinsic change of the magnetic field is needed in this case. Comparing the model with observations, we conclude that the drift direction of $χ$ might oscillate many times during the lifetime of a pulsar. The evolution of $χ$ is not sufficient to explain the rotation behavior of the Crab pulsar, because the observed $χ$ and $\dotχ$ are inconsistent with the values indicated from the timing residuals using this model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-ray softening during the 2008 outburst of XTE J1810-189

XTE J1810-189 underwent an outburst in 2008, and was observed over $\sim 100$ d by RXTE. Performing a time-resolved spectral analysis on the photospheric radius expansion burst detected on 2008 May 4, we obtain the source distance in the range of 3.5--8.7 kpc for the first time. During its outburst, XTE J1810-189 did not enter into the high/soft state, and both the soft and hard colours decreased with decreasing flux. The fractional rms remained at high values ($\sim 30$ per cent). The RXTE/PCA spectra for 3-25 keV can be described by an absorbed power-law component with an additional Gaussian component, and the derived photon index $Γ$ increased from $1.84\pm0.01$ to $2.25\pm0.04$ when the unabsorbed X-ray luminosity in 3-25 keV dropped from $4\times10^{36}$ ergs s$^{-1}$ to $6\times10^{35}$ ergs s$^{-1}$. The relatively high flux, dense observations and broadband spectra allow us to provide strong evidence that the softening behaviour detected in the outburst of XTE J1810-189 originates from the evolution of non-thermal component rather than the thermal component (i.e. neutron star surface emission).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A state-dependent influence of the type-I bursts on the accretion in 4U 1608--52?

We investigated the possible feedback of type-I burst to the accretion process during the spectral evolution of the atoll source 4U 1608--52. By fitting the burst spectrum with a blackbody and an adjustable, persistent spectral component, we found that the latter is significant state-dependent. In the banana state the persistent flux increases along the burst evolution, while in the island state this trend holds only when the bursts are less luminous and starts to reverse at higher burst luminosities. We speculate that, by taking into account both the Poynting-Robertson drag and radiation pressure, these phenomena may arise from the interactions between the radiation field of the type-I burst and the inner region of the accretion disc.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Common Mistakes in Writing Astronomy and Physics Literature in English

This is the 3rd version with major updates and revisions to the 2nd version of the same title. In this article we have collected and corrected some common mistakes made by Chinese students in writing astronomy and physics literature in English. Brief explanations of these mistakes are given in Chinese. We plan to continue to update this collection periodically. Comments, suggestions and criticisms are welcome.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interpreting the Spin-down Evolutions of Isolated Neutron Stars with Hall Effects

The observed long-term spin-down evolution of isolated radio pulsars cannot be explained by the standard magnetic dipole radiation with a constant braking torque. However, how and why the torque varies still remains controversial, which is a major issue in understanding neutron stars. Many pulsars have been observed with significant long-term changes of their spin-down rates modulated by quasi-periodic oscillations. Applying the phenomenological model of pulsar timing noise we developed recently to the observed precise pulsar timing data of isolated neutron stars, here we show that the observed long-term evolutions of their spin-down rates and quasi-periodic modulations can be explained by Hall effects in their crusts. Therefore the evolution of their crustal magnetic fields, rather than that in their cores, dominates the observed long term spin-down evolution of these young pulsars. Understanding of the nature of pulsar timing noise not only reveals the interior physics of neutron stars, but also allows physical modeling of pulsar spin-down and thus improves the sensitivity of gravitational wave detections with pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Jet Properties of GeV-Selected Radio-Loud Narrow-Line Seyfert 1 Galaxies and Possible Connection to Their Disk and Corona

The observed spectral energy distributions of five GeV-selected narrow-line Seyfert 1 (NLS1) galaxies are fitted with a model including the radiation ingredients from the relativistic jet, the accretion disk, and the corona. We compare the properties of these GeV NLS1 galaxies with flat spectrum radio quasars (FSRQs), BL Lacertae objects (BL Lacs) and radio-quiet (RQ) Seyfert galaxies, and explore possible hints for jet-disk/corona connection. Our results show that the radiation physics and the jet properties of the GeV NLS1 galaxies resemble that of FSRQs. The luminosity variations of PMN J0948+0022 and 1H 0323+342 at GeV band is tightly correlated with the beaming factor (δ), similar to that observed in FSRQ 3C 279. The accretion disk luminosities and the jet powers of the GeV NLS1 galaxies cover both the ranges of FSRQs and BL Lacs. With detection of bright corona emission in 1H 0323+342, we show that the ratio of the corona luminosity (L_corona) to the accretion disk luminosity (L_d) is marginally within the high end of this ratio distribution for a RQ Seyfert galaxy sample, and the variation of jet luminosity may connect with Lcorona. However, it is still unclear whether a system with a high L_corona/L_d ratio prefers to power a jet.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Limits on the strength of individual gravitational wave sources using high-cadence observations of PSR B1937+21

We present the results of a search for gravitational waves (GWs) from individual sources using high cadence observations of PSR B1937+21. The data were acquired from an intensive observation campaign with the Lovell telescope at Jodrell Bank, between June 2011 and May 2013. The almost daily cadence achieved, allowed us to be sensitive to GWs with frequencies up to $4.98\times10^{-6}\,\rm {Hz}$, extending the upper bound of the typical frequency range probed by Pulsar Timing Arrays. We used observations taken at three different radio frequencies with the Westerbork Synthesis Radio Telescope in order to correct for dispersion measure effects and scattering variances. The corrected timing residuals exhibited an unmodeled periodic noise with an amplitude $~150\,\rm {ns}$ and a frequency of $3.4\rm {yr}^{-1}$. As the signal is not present in the entire data set, we attributed it to the rotational behaviour of the pulsar, ruling out the possibilities of being either due to a GW or an asteroid as the cause. After removing this noise component, we placed limits on the GW strain of individual sources equaling to $h_{\rm s}=1.53\times10^{-11}$ and $h_{\rm s}=4.99\times10^{-14}$ at $10^{-7}\,\rm {Hz}$ for random and optimal sources locations respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Possible hard X-ray shortages in bursts from KS 1731-260 and 4U 1705-44

Aims: A hard X-ray shortage, implying the cooling of the corona, was observed during bursts of IGR J17473-272, 4U 1636-536, Aql X-1, and GS 1826-238. Apart from these four sources, we investigate here an atoll sample, in which the number of bursts for each source is larger than 5, to explore the possible additional hard X-ray shortage during {\it Rossi X-ray timing explorer (RXTE)} era. Methods: According to the source catalog that shows type-I bursts, we analyzed all the available pointing observations of these sources carried out by the {\it RXTE} proportional counter array (PCA). We grouped and combined the bursts according to their outburst states and searched for the possible hard X-ray shortage while bursting. Results: We found that the island states of KS 1731-260 and 4U 1705-44 show a hard X-ray shortage at significant levels of 4.5 and 4.7 $σ$ and a systematic time lag of $0.9 \pm 2.1$ s and $2.5 \pm 2.0$ s with respect to the soft X-rays, respectively. While in their banana branches and other sources, we did not find any consistent shortage.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic Jet Properties of GeV-TeV Blazars and Possible Implications for the Jet Formation, Composition, and Cavity Kinematics

We fit the spectral energy distributions (SEDs) of a GeV-TeV FSRQ sample with the leptonic model. Their gamma_min of the relativistic electron distributions, which significantly affect the estimates of the jet properties, are constrained, with a typical value of 48. Their jet power, magnetized parameter, radiation efficiency, and jet production/radiation rates per central black hole (BH) mass are derived and compared to that of BL Lacs. We show that he FSRQ jets may be dominated by the Poynting flux and have a high radiation efficiency, whereas the BL Lac jets are likely dominated by particles and have a lower radiation efficiency than FSRQs. Being different from BL Lacs, the jet powers of FSRQs are proportional to their central BH masses. The jet production and radiation rates of the FSRQs distribute in narrow ranges and are correlated with each other, whereas no similar feature is found for the BL Lacs. We also show that the jet power is correlated with the cavity kinetic power, and the magnetic field energy in the jets may provide the cavity kinetic energy of FSRQs and the kinetic energy of cold protons in the jets may be crucial for cavity kinetic energy of BL Lacs. We suggest that the dominating formation mechanism of FSRQ jets may be the BZ process, but BL Lac jets may be produced via the BP and/or BZ processes, depending on the structures and accretion rates of accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotational Behaviors and Magnetic Field Evolution of Radio Pulsars

The observed long-term spin-down evolution of isolated radio pulsars cannot be explained by the standard magnetic dipole radiation with a constant braking torque. However how and why the torque varies still remains controversial, which is an outstanding problem in our understanding of neutron stars. We have constructed a phenomenological model of the evolution of surface magnetic fields of pulsars, which contains a long-term decay modulated by short-term oscillations; a pulsar's spin is thus modified by its magnetic field evolution. The predictions of this model agree with the precisely measured spin evolutions of several individual pulsars; the derived parameters suggest that the Hall drift and Hall waves in the NS crusts are probably responsible for the long-term change and short-term quasi-periodical oscillations, respectively. Many statistical properties of the timing noise of pulsars can be well re-produced with this model, including correlations and the distributions of the observed braking indices of the pulsars, which span over a range of more than 100 millions. We have also presented a phenomenological model for the recovery processes of classical and slow glitches, which can successfully model the observed slow and classical glitch events without biases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Magneto Hydro Dynamical Model of KHz Quasi Periodic Oscillations in Neutron Star Low Mass X-ray Binaries (II)

We study the kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low mass X-ray binaries (LMXBs) with a new magnetohydrodynamics (MHD) model, in which the compressed magnetosphere is considered. The previous MHD model (Shi \& Li 2009) is re-examined and the relation between the frequencies of the kHz QPOs and the accretion rate in LMXBs is obtained. Our result agrees with the observations of six sources (4U 0614+09, 4U 1636--53, 4U 1608--52, 4U 1915--15, 4U 1728--34, XTE 1807--294) with measured spins. In this model the kHz QPOs originate from the MHD waves in the compressed magnetosphere. The single kHz QPOs and twin kHz QPOs are produced in two different parts of the accretion disk and the boundary is close to the corotation radius. The lower QPO frequency in a frequency-accretion rate diagram is cut off at low accretion rate and the twin kHz QPOs encounter a top ceiling at high accretion rate due to the restriction of innermost stable circular orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

X-ray Scattered Halo around IGR J17544-2619

X-ray photons coming from an X-ray point source not only arrive at the detector directly, but also can be strongly forward-scattered by the interstellar dust along the line of sight (LOS), leading to a detectable diffuse halo around the X-ray point source. The geometry of small angle X-ray scattering is straightforward, namely, the scattered photons travel longer paths and thus arrive later than the unscattered ones; thus the delay time of X-ray scattered halo photons can reveal information of the distances of the interstellar dust and the point source. Here we present a study of the X-ray scattered around IGR J17544-2619, which is one of the so-called supergiant fast X-ray transients. IGR J17544-2619 underwent a striking outburst when observed with Chandra on 2004 July 3, providing a near delta-function lightcurve. We find that the X-ray scattered halo around IGR J17544-2619 is produced by two interstellar dust clouds along the LOS. The one which is closer to the observer gives the X-ray scattered at larger observational angles; whereas the farther one, which is in the vicinity of the point source, explains the halo with a smaller angular size. By comparing the observational angle of the scattered halo photons with that predicted by different dust grain models, we are able to determine the normalized dust distance. With the delay times of the scattered halo photons, we can determine the point source distance, given a dust grain model. Alternatively we can discriminate between the dust grain models, given the point source distance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black Hole Binaries and Microquasars

This is a general review on the observations and physics of black hole X-ray binaries and microquasars, with the emphasize on recent developments in the high energy regime. The focus is put on understanding the accretion flows and measuring the parameters of black holes in them. It includes mainly two parts: (1) Brief review of several recent review article on this subject; (2) Further development on several topics, including black hole spin measurements, hot accretion flows, corona formation, state transitions and thermal stability of standard think disk. This is thus not a regular bottom-up approach, which I feel not necessary at this stage. Major effort is made in making and incorporating from many sources useful plots and illustrations, in order to make this article more comprehensible to non-expert readers. In the end I attempt to make a unification scheme on the accretion-outflow (wind/jet) connections of all types of accreting BHs of all accretion rates and all BH mass scales, and finally provide a brief outlook.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Blazar Anti-Sequence of Spectral Variation within individual Blazars: Cases for Mrk 501 and 3C 279

The jet properties of Mrk 501 and 3C 279 are derived by fitting the broadband spectral energy distributions with the lepton models. The derived gamma_b (the break Lorenz factor of the electron distribution) are 10^4-10^6 for Mrk 501 and 200~600 for 3C 279 and the magnetic field strength (B) of Mrk 501 is usually one order of magnitude lower than that of 3C 279, but their Doppler factors (delta) are comparable. A spectral variation feature that the peak luminosity is correlated with the peak frequency, which is opposite to the blazar sequence, is observed in the two sources. We find that (1) the peak luminosities of the two bumps in SEDs for Mrk 501 depend on gamma_b in both the observer and co-moving frames, but they are not correlated with B and delta; (2)the luminosity variation of 3C 279 is dominated by the external Compton (EC) peak and its peak luminosity is correlated with gamma_b and delta, but anti-correlated with B. These results suggest that gamma_b may govern the spectral variation of Mrk 501 and delta and B would be responsible for the spectral variation of 3C 279. The narrow distribution of gamma_b and the correlation of gamma_b and B in 3C 279 would be due to the cooling by the EC process and the strong magnetic field. Based on our brief discussion, we propose that this spectral variation feature may be originated from the instability of the corona but not from the variation of the accretion rate as that for the blazar sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct measurement of evolving dark energy density and super-accelerating expansion of the universe

A higher value of Hubble constant has been obtained from measurements with nearby Type Ia supernovae, than that obtained at much higher redshift. With the peculiar motions of their hosts, we find that the matter content at such low redshift is only about 10% of that at much higher redshifts; such a low matter density cannot be produced from density perturbations in the background of the ΛCDM expansion. Recently the Planck team has reported a lower Hubble constant and a higher matter content. We find that the dark energy density increases with cosmic time, so that its equation-of-state parameter decreases with cosmic time and is less than -1 at low redshift. Such dark energy evolution is responsible for driving the super-accelerating expansion of the universe. In this extended ΛCDM model, the cosmological redshift represents time rather than radial coordinate, so that the universe complies to the Copernican Principle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of oxygen and nitrogen abundances and nitrogen production mechanism in massive star-forming galaxies

Utilizing the observational data of 55,318 star-forming galaxies (SFGs) selected from the catalog of MPA-JHU emission-line measurements for the SDSS DR8, we investigate the galaxy downsizing effect of their O and N enrichments, and the nitrogen production mechanism in them. We show the redshift evolution of O and N abundances and specific star formation rates for different galaxy mass ranges, demonstrating the galaxy downsizing effect caused by less massive progenitors of less massive galaxies. The O and N abundances do not remain constant for different galaxy mass ranges, and the enrichment (and hence star formation) decreases with increasing galaxy stellar mass. We find evidence of the O enrichment for galaxies with stellar masses $M_{*}>10^{11.0} $ (in units of $M_{\odot}$), i.e. $Δ({\log}({\rm O/H})) \sim 0.10$ and $Δ({\log}({\rm N/H})) \sim 0.28$ from redshift 0.023 to 0.30. Based on the evolutionary schematic model of N/O ratios in Coziol et al., who proposed the scheme that the production of nitrogen is the consequence of a sequence of bursts in SFGs, we conclude that the nitrogen production is dominated by the intermediate-mass stars, which dominate the secondary synthesis in SFGs. However, for galaxies with $M_{*}>10^{10.35} $ we find evidence of enhanced N/O abundance ratios, which are significantly above the secondary synthesis line. This suggests that outflows of massive stars, which deplete oxygen efficiently, are more important in massive galaxies. Finally we find an excellent linear relation between $M_{*}$ and log(N/O), indicating that the N/O abundance ratio is a good indicator of the stellar mass in a SFG and may be used as a standard candle for studying cosmology, if confirmed with further studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling the Evolution and Distribution of the Frequency's Second Derivative and Braking Index of Pulsar Spin with Simulations

We model the evolution of spin frequency's second derivative $\ddotν$ and braking index $n$ of radio pulsars with simulations within the phenomenological model of their surface magnetic field evolution, which contains a long-term decay modulated by short-term oscillations. For the pulsar PSR B0329+54, the model can reproduce the main characteristics of its $\ddotν$ variation with oscillation periods, predicts another $\sim 50$ yr oscillation component and another recent swing of the sign of $\ddotν$. We show that the "averaged" $n$ is different from the instantaneous $n$, and its oscillation magnitude decreases abruptly as the time span increases, due to the "averaging" effect. The simulation predicted timing residuals agree with the main features of the reported data. We further perform Monte Carlo simulations for the distribution of the reported data in $|\ddotν|$ versus characteristic age $τ_{\rm c}$ diagram. The model with a power law index $α=0.5$ can reproduce the slope of the linear fit to pulsars' distributions in the diagrams of $\log|\ddotν|-\logτ_{\rm c}$ and $\log|n|-\logτ_{\rm c}$, but the oscillations are responsible for the almost equal number of positive and negative values of $\ddotν$, in agreement with our previous analytical studies; an oscillation period of about several decades is also preferred. However the range of the oscillation amplitudes is $-11.4\lesssim\log f\lesssim-10.2$, slightly lager than the analytical prediction, $\log f\simeq-11.85$, because the "averaging" effect was not included previously.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Relaxation Behaviors of Slow and Classical Glitches: Observational Biases and Their Opposite Recovery Trends

We study the pulsar timing properties and the data analysis methods during glitch recoveries. In some cases one first fits the time-of-arrivals (TOAs) to obtain the "time-averaged" frequency $ν$ and its first derivative $\dotν$, and then fits models to them. However, our simulations show that $ν$ and $\dotν$ obtained this way are systematically biased, unless the time intervals between the nearby data points of TOAs are smaller than about $10^4$ s, which is much shorter than typical observation intervals. Alternatively, glitch parameters can be obtained by fitting the phases directly with relatively smaller biases; but the initial recovery timescale is usually chosen by eyes, which may introduce a strong bias. We also construct a phenomenological model by assuming a pulsar's spin-down law of $\dotνν^{-3} =-H_0 G(t)$ with $G(t)=1+κe^{-t/τ}$ for a glitch recovery, where $H_0$ is a constant and $κ$ and $τ$ are the glitch parameters to be found. This model can reproduce the observed data of slow glitches from B1822--09 and a giant classical glitch of B2334+61, with $κ<0$ or $κ>0$, respectively. We then use this model to simulate TOA data and test several fitting procedures for a glitch recovery. The best procedure is: 1) use a very high order polynomial (e.g. to 50th order) to precisely describe the phase; 2) then obtain $ν(t)$ and $\dotν(t)$ from the polynomial; and 3) the glitch parameters are obtained from $ν(t)$ or $\dotν(t)$. Finally, the uncertainty in the starting time $t_0$ of a classical glitch causes uncertainties to some glitch parameters, but less so to a slow glitch and $t_0$ of which can be determined from data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiation Mechanism and Jet Composition of Gamma-Ray Bursts and GeV-TeV selected Radio Loud Active Galactic Nuclei

Gamma-ray bursts (GRBs) and GeV-TeV selected radio loud Active Galactic Nuclei (AGNs) are compared based on our systematic modeling of the observed spectral energy distributions of a sample of AGNs with a single-zone leptonic model. We show that the correlation between the jet power (P_{jet}) and the prompt gamma-ray luminosity (L_{jet}) of GRBs is consistent, within the uncertainties, with the correlation between jet power and the synchrotron peak luminosity (L_{s, jet}) of flat spectrum radio quasars (FSRQs). Their radiation efficiencies (\varepsilon) are also comparable (>10% for most sources), which increase with the bolometric jet luminosity (L_{bol,jet}) for FSRQs and with the L_{jet} for GRBs with similar power-law indices. BL Lacs do not follow the P_{jet}-L_{s, jet} relation of FSRQs. They have lower \varepsilon and L_{bol, jet} values than FSRQs, and a tentative L_{bol, jet}-\varepsilon relation is also found, with a power-law index being different from that of the FSRQs. The magnetization parameters (σ) of FSRQs are averagely larger than that of BL Lacs. They are anti-correlated with $\varepsilon$ for the FSRQs, but positive correlated with \varepsilon for the BL Lacs. GeV Narrow-line Seyfert 1 galaxies potentially share similar properties with FSRQs. Based on the analogy between GRBs and FSRQs, we suggest that the prompt gamma-ray emission of GRBs is likely produced by synchrotron process in a magnetized jet with high radiation efficiency, similar to FSRQs. The jets of BL Lacs, on the other hand, are less efficient and are likely more matter dominated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

State-Dependent Orbital Modulation of X-rays in CYG X-3

We analyze all of the available {\it RXTE} observations of Cyg X-3 in order to investigate the connection between the central X-ray source and its surrounding environment. The hardness--intensity diagram of Cyg X-3 displays a "shoe" shape rather than the Q-type shape commonly seen in other black hole X-ray binary, and exhibits no apparent hysteresis effect. During the $γ$-ray outbursts, no existing data are located in the hard and intermediate states, which suggest the absence of a significant population of non-thermal electrons when the source is in these states. For the first time, we present the orbital modulation of the X-ray light curve (LC) of all five states. The different energy band LCs are in phase with each other in all five states, and the modulation amplitude of both soft and hard X-ray LCs monotonously increases with decreasing hardness from hard to soft non-thermal states. We confirm that the modulation depth decreases with increasing energy in the hard, intermediate, and very high states, as originally reported by Zdziarski et al. However, in the soft non-thermal state, the hard X-ray modulation strength significantly increases and is even larger than the soft X-ray one. Our results rule out both wind absorption and jet origins of the hard X-ray LC modulation in the soft non-thermal state, and challenge our understanding of the states of Cyg X-3.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The hard X-ray behavior of Aql X-1 during type-I bursts

We report the discovery of an anti-correlation between the soft and the hard X-ray lightcurves of the X-ray binary Aql X-1 when bursting. This behavior may indicate that the corona is cooled by the soft X-ray shower fed by the type-I X-ray bursts, and that this process happens within a few seconds. Stacking the Aql X-1 lightcurves of type-I bursts, we find a shortage in the 40--50 keV band, delayed by 4.5$\pm$1.4 s with respect to the soft X-rays. The photospheric radius expansion (PRE) bursts are different in that neither a shortage nor an excess shows up in the hard X-ray lightcurve.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Type-I X-ray bursts reveal a fast co-evolving behavior of the corona in an X-ray binary

The coronae in X-ray binaries (XRBs) still remain poorly understood, although they have been believed for a long time to play a key role in modeling the characteristic outbursts of XRBs. Type-I X-ray bursts, the thermonuclear flashes happening on the surface of a neutron star (NS), can be used as a probe to the innermost region of a NS XRB, where the corona is believed to be located very close to the NS. We report the discovery of a tiny life cycle of the corona that is promptly co-evolved with the type-I bursts superimposed on the outburst of the NS XRB IGR J17473$-$2721. This finding may serve as the first evidence of directly seeing the rapid disappearance and formation of a corona in an XRB with a cooling/heating timescale of less than a second, which can strongly constrain the accretion models in XRBs at work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wavelength Measurements of K Transitions of Oxygen, Neon, and Magnesium with X-ray Absorption Lines

Accurate atomic transition data are important in many astronomical research areas especially in line spectroscopy study. Whereas transition data of He-like and H-like ions (i.e., ions at high-charge states) are accurately calculated, that of K transitions of neutral or low-ionized metal elements are still very uncertain. Spectroscopy of absorption lines produced in the interstellar medium (ISM) has been proven to be an effective way to measure the central wavelengths of these atomic transitions. In this work we analyze 36 Chandra High Energy Transmission Grating observations and search for and measure the ISM absorption lines along sight lines of 11 low-mass X-ray binaries. We correct the Galactic rotation velocity to the rest frame for every observation and then use two different methods to merge all the corrected spectra to a co-added spectrum. However the co-added spectra obtained by these methods exhibit biases, either to the observations of high counts or high signal-to-noise ratios of the lines. We make Bayesian analysis to several significantly detected lines to obtain the systematic uncertainty and the bias correction of other lines. Compared to previous studies (e.g., Yao et al. 2009), our results improve the accuracy of wavelengths by a factor from two to five and significantly reduce the systematic uncertainties and biases. Several weak transitions (e.g., 1s--2p of Mg IV and Mg V; 1s--3p of Mg III and Mg V) are also detected for the first time, albeit with low significance; future observations with improved accuracy are required to confirm these detections.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

X-ray bursts as a probe of the corona: the case of XRB 4U 1636-536

To investigate the possible cooling of the corona by soft X-rays bursts, we have studied 114 bursts embedded in the known X-ray evolution of 4U 1636-536. We have grouped these bursts according to the ratio of the flux in the 1.5--12 keV band with respect to that in the 15--50 keV band, as monitored by RXTE/ASM and Swift/BAT, respectively. We have detected a shortage at hard X-rays while bursting. This provides hints for a corona cooling process driven by soft X-rays fed by the bursts that occurred on the surface of neutron star. The flux shortage at 30--50 keV has a time lag of 2.4$\pm$1.5 seconds with respect to that at 2--10 keV, which is comparable to that of 0.7$\pm$0.5 seconds reported in bursts of IGR 17473-2721. We comment on the possible origin of these phenomena and on the implications for the models on the location of the corona.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Blazar Anti-Sequence of Spectral Variability for Individual TeV Blazars

We compile from literature the broadband SEDs of twelve TeV blazars observed simultaneously or quasi-simultaneously with Fermi/LAT and other instruments. Two SEDs are available for each of the objects and the state is identified as a low or high state according to its flux density at GeV/TeV band. The observed SEDs of BL Lac objects (BL Lacs) are fitted well with the synchrotron + synchrotron-self-Compton (syn+SSC) model, whereas the SEDs of the two flat spectrum radio quasars (FSRQs) need to include the contributions of external Compton scattering. In this scenario, it is found that the Doppler factor delta of FSRQs is smaller than that of BL Lacs, but the magnetic field strength B of FSRQs is larger than that of BL Lacs. The increase of the peak frequency of the SEDs is accompanied with the increase of the flux for the individual sources, which seems opposite to the observational phenomena of the blazar sequence. We refer this phenomenon to blazar anti-sequence of spectral variability for individual TeV blazars. However, both the blazar sequence from FSRQs to BL Lacs and blazar anti-sequence of the spectral variability from low state to high state are accompanied by an increase of the break Lorentz factor of the electron's spectrum gamma_b and a decrease of B. We propose a model in which the mass accretion rate is the driving force behind both the blazar sequence for ensembles of blazars and the blazar anti-sequence for individual blazars. Specifically we suggest that the differences in mass accretion rate of different blazars produce the observed blazar sequence, but ΔM in each blazar results in the observed blazar anti-sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring the black hole masses in accreting X-ray binaries by detecting the Doppler orbital motion of their accretion disk wind absorption lines

So far essentially all black hole masses in X-ray binaries have been obtained by observing the companion star's velocity and light curves as functions of the orbital phase. However a major uncertainty is the estimate of the orbital inclination angle of an X-ray binary. Here we suggest to measure the black hole mass in an X-ray binary by measuring directly the black hole's orbital motion, thus obtaining the companion to black hole mass ratio. In this method we assume that accretion disk wind moves with the black hole and thus the black hole's orbital motion can be obtained from the Doppler velocity of the absorption lines produced in the accretion disk wind. We validate this method by analyzing the Chandra/HETG observations of GRO J1655-40, in which the black hole orbital motion with line of sight velocity of 90.8 (+-11.3) km/s, inferred from the Doppler velocity of disk-wind absorption lines, is consistent with the prediction from its previously measured system parameters. We obtain the black hole mass of 5.41 (+0.98, -0.57) solar masses and system inclination of 72.0 (+7.8, -7.5) degrees in GRO J1655-40. Additional observations of this source covering more orbital phases can improve estimates on its system parameters substantially. We then apply the method to the black hole X-ray binary LMC X-3 observed with HST/COS near orbital phase 0.75. We find that the disk-wind absorption lines of CIV doublet were shifted to about 50 km/s, which yields a companion-to-black-hole mass ratio of 0.6 for an assumed disk wind velocity of -400 km/s. Additional observations covering other orbital phases (0.25 in particular) are crucial to ease this assumption and then to directly constrain the mass ratio. This method in principle can also be applied to any accreting compact objects with detectable accretion disk wind absorption line features.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling pulsar time noise with long term power law decay modulated by short term oscillations of the magnetic fields of neutron stars

We model the evolution of the magnetic fields of neutron stars as consisting of a long term power-law decay modulated by short term small amplitude oscillations. Our model predictions on the timing noise $\ddotν$ of neutron stars agree well with the observed statistical properties and correlations of normal radio pulsars. Fitting the model predictions to the observed data, we found that their initial parameter implies their initial surface magnetic dipole magnetic field strength ~ 5E14 G at ~0.4 year old and that the oscillations have amplitude between E-8 to E-5 and period on the order of years. For individual pulsars our model can effectively reduce their timing residuals, thus offering the potential of more sensitive detections of gravitational waves with pulsar timing arrays. Finally our model can also re-produce their observed correlation and oscillations of the second derivative of spin frequency, as well as the "slow glitch" phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On a common misunderstanding of the Birkhoff theorem and light deflection calculation: generalized Shapiro delay and its possible laboratory test

In Newtonian gravity (NG) it is known that the gravitational field anywhere inside a spherically symmetric distribution of mass is determined only by the enclosed mass. This is also widely believed to be true in general relativity (GR), and the Birkhoff theorem is often invoked to support this analogy between NG and GR. Here we show that such an understanding of the Birkhoff theorem is incorrect and leads to erroneous calculations of light deflection and delay time through matter. The correct metric, matching continuously to the location of an external observer, is determined both by the enclosed mass and mass distribution outside. The effect of the outside mass is to make the interior clock run slower, i.e., a slower speed of light for external observer. We also discuss the relations and differences between NG and GR, in light of the results we obtained in this Lettework. Finally we discuss the Generalized Shapiro delay, caused by the outside mass, and its possible laboratory test.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Progenitors of Magnetars and Hyperaccreting Magnetized Disks

We propose that a magnetar could be formed during the core collapse of massive stars or coalescence of two normal neutron stars, through collecting and inheriting the magnetic fields magnified by hyperaccreting disk. After the magnetar is born, its dipole magnetic fields in turn have a major influence on the following accretion. The decay of its toroidal field can fuel the persistent X-ray luminosity of either an SGR or AXP; however the decay of only the poloidal field is insufficient to do so.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasars' Optical Polarization and Balmer Edge Feature Revealed by Ultra-Violet, and Polarized Visible to Near Infrared Emissions

Polarized emission from a quasar is produced by wavelength-independent electron scattering surrounding its accretion disc, and thus avoid the contamination from its host galaxy and reveal the intrinsic emission spectrum of the accretion disc. Ultra-violet (UV) emission from a quasar is normally free from the contamination from its host galaxy. Polarization fraction of the quasar's disc emission can therefore be determined by comparing total UV emission with polarized visible to near-infrared (NIR) emission; and the resulting continuum spectrum from UV to infrared can reveal the theoretically expected Balmer edge absorption feature. We fit the polarized spectra in visible and NIR bands together with the total UV spectra of two type-1 quasars (3C 95, 4C 09.72), to an extended geometrically thin and optically thick accretion disc model. In addition to the standard model, we include the Balmer edge absorption due to co-rotational neutral gas on a narrow annulus of the accretion disc. We find that the extended thin accretion disc model provides adequate description on the continuum spectra of the two quasars from UV to NIR wavelengths. A Monte-Carlo-Markov-Chain fitting to the continuum spectra is able to well constrain the true polarization fraction of the disk emission, which allows the Balmer edge feature to be completely revealed from polarized visible to UV continua. The Balmer edge feature is prominent in both quasars' spectra, and is significantly broadened due to the orbital motion of gas in the accretion disc. The broadening of the Balmer edge feature is therefore related to the quasar's inclination. This work proves the concept of determining quasar's inclination from the Balmer edge feature in their continuum spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radial Angular Momentum Transfer and Magnetic Barrier for Short-Type Gamma-Ray Burst Central Engine Activity

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiation Mechanisms and Physical Properties of GeV-TeV BL Lac Objects

BL Lac objects are the best candidates to study the jet properties since their spectral energy distributions (SEDs) are less contaminated by the emission from the accretion disk and external Compton processes. We compile the broadband SEDs observed with Fermi/LAT and other instruments from literature for 24 TeV BL Lac objects. Two SEDs, which are identified as a low or high state according to its flux density at 1 TeV, are available for each of ten objects. These SEDs can be explained well with the synchrotron+synchrotron-self-Compton model. We constrain the magnetic filed strength B and the Doppler factor delta of the radiation region by incorporating the chi^{2}-minimization technique and the gamma-ray transparency condition. Twenty-four out of the 34 SEDs can constrain B and delta in 1 sigma significance level and others only present the limits for both B and delta. The typical values of $B$ and delta are 0.1 ~ 0.6 G and 10 ~ 35, respectively. The derived values of gamma_b are significantly different among sources and even among the low and high states of a given source. Prominent flux variations with a clear spectral shift are observed and a tentative correlation between the ratio of the flux density at 1 TeV and the ratio of gamma_b in the low and high states is presented, likely indicating that the relativistic shocks for the electron acceleration may be responsible for the flux variations and the spectral shift. A weak anti-correlation between the jet power and the mass of the central black hole is observed, i.e., P_{jet} propto M^{-1}_{BH}, which disfavors the scenario of a pure accretion-driven jet. Implications for the blazar sequence and the intergalactic magnetic field from our results are also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Testing Models of Magnetic Field Evolution of Neutron Stars with the Statistical Properties of Their Spin Evolutions

We test models for the evolution of neutron star (NS) magnetic fields (B). Our model for the evolution of the NS spin is taken from an analysis of pulsar timing noise presented by Hobbs et al. (2010). We first test the standard model of a pulsar's magnetosphere in which B does not change with time and magnetic dipole radiation is assumed to dominate the pulsar's spin-down. We find this model fails to predict both the magnitudes and signs of the second derivatives of the spin frequencies ($\ddotν$). We then construct a phenomenological model of the evolution of $B$, which contains a long term decay (LTD) modulated by short term oscillations (STO); a pulsar's spin is thus modified by its B-evolution. We find that an exponential LTD is not favored by the observed statistical properties of $\ddotν$ for young pulsars and fails to explain the fact that $\ddotν$ is negative for roughly half of the old pulsars. A simple power-law LTD can explain all the observed statistical properties of $\ddotν$. Finally we discuss some physical implications of our results to models of the B-decay of NSs and suggest reliable determination of the true ages of many young NSs is needed, in order to constrain further the physical mechanisms of their B-decay. Our model can be further tested with the measured evolutions of $\dotν$ and $\ddotν$ for an individual pulsar; the decay index, oscillation amplitude and period can also be determined this way for the pulsar.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Why do the braking indices of pulsars span a range of more than 100 millions?

Here we report that the observed braking indices of the 366 pulsars in the sample of Hobbs et al. range from about $-10^8$ to about $+10^8$ and are significantly correlated with their characteristic ages. Using the model of magnetic field evolution we developed previously based on the same data, we derived an analytical expression for the braking index, which agrees with all the observed statistical properties of the braking indices of the pulsars in the sample of Hobbs et al. Our model is, however, incompatible with the previous interpretation that magnetic field growth is responsible for the small values of braking indices ($<3$) observed for "baby" pulsars with characteristic ages of less than $2\times 10^3$ yr. We find that the "instantaneous" braking index of a pulsar may be different from the "averaged" braking index obtained from fitting the data over a certain time span. The close match between our model-predicted "instantaneous" braking indices and the observed "averaged" braking indices suggests that the time spans used previously are usually smaller than or comparable to their magnetic field oscillation periods. Our model can be tested with the existing data, by calculating the braking index as a function of the time span for each pulsar. In doing so, one can obtain for each pulsar all the parameters in our magnetic field evolution model, and may be able to improve the sensitivity of using pulsars to detect gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accretion states of the Galactic microquasar GRS 1758-258

We present the results of a radio and X-ray study of the Galactic microquasar GRS 1758-258, using unpublished archival data and new observations. We focus in particular on the 2000-2002 state transitions, and on its more quiet behaviour in 2008-2009. Our spectral and timing analysis of the XMM-Newton data shows that the source was in the canonical intermediate, soft and hard states in 2000 September 19, 2001 March 22 and 2002 September 28, respectively. We estimate the disk size, luminosity and temperature, which are consistent with a black hole mass ~10 M_{sun}. There is much overlap between the range of total X-ray luminosities (on average ~0.02 L_{Edd}) in the hard and soft states, and probably between the corresponding mass accretion rates; in fact, the hard state is often more luminous. The extended radio lobes seen in 1992 and 1997 are still present in 2008-2009. The 5-GHz radio core flux density has shown variability between ~0.1-0.5 mJy over the last two decades. This firmly places GRS 1758-258 in the radio-quiet sequence of Galactic black holes, in the radio/X-ray plane. We note that this dichotomy is similar to the dichotomy between the radio/X-ray sequences of Seyfert and radio galaxies. We propose that the different radio efficiency of the two sequences is due to relativistic electron/positron jets in radio-loud black holes, and sub-relativistic, thermally dominated outflows in radio-quiet sources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Astrophysical Black Holes in the Physical Universe

In this chapter I focus on asking and answering the following questions: (1) What is a black hole? Answer: There are three types of black holes, namely mathematical black holes, physical black holes and astrophysical black holes. An astrophysical black hole, with mass distributed within its event horizon but not concentrated at the singularity point, is not a mathematical black hole. (2) Can astrophysical black holes be formed in the physical universe? Answer: Yes, at least this can be done with gravitational collapse. (3) How can we prove that what we call astrophysical black holes are really black holes? Answer: Finding direct evidence of event horizon is not the way to go. Instead I propose five criteria which meet the highest standard for recognizing new discoveries in experimental physics and observational astronomy. (4) Do we have sufficient evidence to claim the existence of astrophysical black holes in the physical universe? Answer: Yes, astrophysical black holes have been found at least in some galactic binary systems, at the center of almost every galaxy, and as the central engines of at least some long gamma-ray bursts. (5) Will all matter in the universe eventually fall into black holes? Answer: Probably "no", because "naked" compact objects, if they do exist with radii smaller than the radii of event horizons for their masses but are not enclosed by event horizons, can rescue the universe from an eternal death by re-cycling out the matter previously accreted into astrophysical black holes. Finally I also discuss briefly if we need a quantum theory of gravity in order to further understand astrophysical black holes, and what further astronomical observations and telescopes are needed to make further progress on our understanding of astrophysical black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comparing the Accretion Disk Evolution of Black Hole and Neutron Star X-ray Binaries from Low to Super-Eddington Luminosity

Low-mass X-ray binaries (LMXBs) are systems in which a low-mass companion transfers mass via Roche-lobe overflow onto a black hole (BH) or a weakly magnetized neutron star (NS). It is believed that both the solid surface and the magnetic field of an NS can affect the accretion flow and show some observable effects. Using the disk emission dominant data, we compare the disk evolution of the two types of systems from low luminosity to super-Eddington luminosity. As the luminosity decreases the disk in the NS LMXB 4U1608--522 begins to leave the innermost stable circular orbit (ISCO) at much higher luminosity ($\sim$ 0.1 $L_{\mathrm{Edd}}$), compared with BH LMXBs at much lower luminosity ($\sim$ 0.03 $L_{\mathrm{Edd}}$), due to the interaction between the NS magnetosphere and accretion flow. However, as the luminosity increases above a critical luminosity, the disks in BH and NS LMXBs trace the same evolutionary pattern, because the magnetosphere is restricted inside ISCO, and then both the NS surface emission and (dipole) magnetic field do not significantly affect the secular evolution of the accretion disk, that is driven by the increased radiation pressure in the inner region. We further suggest that the NS surface emission provides additional information of accretion disk, not available in BH systems. Through the observed NS surface emission, we argue that the disk thickness $H/R$ is less than $0.3-0.4$, and that the significant outflow from inner disk edge exists at luminosity close to Eddington luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic Field Decay Makes Neutron Stars Look Older Than They Are

It is commonly accepted that a neutron star is produced, when a massive star exhausts its nuclear fuel and ends its life in a core-collapse supernova explosion. This scenario is confirmed by the detection of pulsars, which are believed to be rapidly spinning neutron stars, in the central regions of many supernova remnants. Neutron stars and their associated supernova remnants should therefore have the same ages. As expected, the age of the Crab Pulsar, the first to be connected with a supernova remnant (the Crab Nebula), can be inferred from its current spin period and its derivative and indeed has about the same age of the supernova remnant, that was produced from a historically recorded supernova explosion in 1054. However most neutron stars appear to be much older than the ages of their associated supernova remnants, a puzzle not yet understood. Another puzzle is that so far no convincing evidence has been found in favor of magnetic field decay in neutron stars, that is predicted in most models of neutron stars. Here we show convincing evidence of magnetic field decay in some young neutron stars, and that the magnetic field decay can alter their spinning behaviors significantly such that these neutron stars appear much older than they really are.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Power-law Magnetic Field Decay and Constant Core Temperatures of Magnetars, Normal and Millisecond Pulsars

The observed correlations, between the characteristic ages and dipole surface magnetic field strengths of all pulsars, can be well explained by magnetic field decay with core temperatures of $~2\times10^{8}$ K, $\sim2\times10^{7}$ K, and $\sim10^{5}$ K, for magnetars, normal radio pulsars, and millisecond pulsars, respectively; assuming that their characteristic ages are about two orders of magnitude larger than their true ages, the required core temperatures may be reduced by about a factor of 10. The magnetic decay follows a power-law and is dominated by the solenoidal component of the ambipolar diffusion mode. In this model, all NSs are assumed to have the same initial magnetic field strength, but different core temperature which do not change as the magnetic field decays. This suggests that the key distinguishing property between magnetars and normal pulsars is that magnetars were born much hotter than normal pulsars, and thus have much longer magnetic field decay time scales, resulting in higher surface magnetic field strength even with the same ages of normal pulsars. The above conclusion agrees well with the observed correlations between the surface temperatures of magnetars and other young NSs, which do not agree with the cooling dominated evolution of neutron stars. This suggests a possible scenario that heating, perhaps due to magnetic field decay, balances neutron star cooling for observed pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation Mechanisms and Physical Properties of GeV-TeV BL Lac Objects

Broadband spectral energy distributions (SEDs) are complied from literature for 24 TeV BL Lac objects. Two SEDs, which are identified as low and high states with their flux densities at $1\ {\rm TeV}$ ($F_{1 TeV}$), are available for each of 11 objects. We fit the SEDs with the single-zone synchrotron + synchrotron-self-Compton (SSC) model and determine the physical parameters of the jets. The model well represents the SEDs, and the observed relation between two peak frequencies $ν_{\rm s}$ and $ν_{\rm c}$ also favors the model. It is found that the break of the electron spectrum ($γ_{\rm b}$) is different among sources, but magnetic field strength distributes narrowly in $0.1-0.6 $ G, indicating that the shocks in the jets are significantly different among sources and the magnetic field may be independent of the shocks. The ratio of $F_{1 TeV}$ is correlated with the ratio of the $γ_{\rm b}$ in the low and high states, indicating that the relativistic shocks in the jets may be responsible for the observed flux variations and the spectral shift. The ratio of the two peak luminosities $L_{\rm c}/L_{\rm s}$ is anti-correlated with $ν_{\rm s}$ in the co-moving frame, which is possibly due to the Klein-Nishina effect. No excess in the GeV band due to the interaction between the TeV photons and the extragalactic background light is observed, implying that the strength of the intergalactic magnetic field would be much larger than $ 10^{-16}$ G. An anti-correlation between jet power and the mass of the central black hole is observed, i.e., $P_{\rm jet}\propto M^{-1}_{\rm BH}$, disfavoring the scenario of a pure accretion-driven jet. We suggest that the spin energy extraction may be significant for powering jets in these sources, implying that the black holes with smaller mass may spin more rapidly.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The 2008 outburst of IGR J17473--2721: evidence for a disk corona?

The 2008 outburst of the atoll source IGR J17473--2721 was observed by INTEGRAL, RXTE and Swift. Tens of type-I X-ray bursts were found in this outburst. Joint observations provide sufficient data to look into the behavior of IGR J17473--2721 at the rising part of the outburst. We find that the joint energy spectra can be well fitted with a model composed of a blackbody and a cutoff power-law, with a cutoff energy decreasing from $ \sim$ 150 keV to $\sim$ 40 keV as the source leaves the quiescent state toward the low hard state. This fits into a scenario in which the corona is cooled by the soft X-rays along the outburst evolution, as observed in several other atoll sources. By using the flux measured in the 1.5--30 keV band of the type-I bursts during the outburst, we find that the linear relationship between the burst duration and the flux still holds for those bursts that occur at the decaying part of the low hard state, but with a different slope than the overall one that was estimated with the bursts happening in the whole extent of, and for the rest of the low hard state. The significance of such a dichotomy in the type-I X-ray bursts is $\sim$ 3 $σ$ under an F-test. Similar results are hinted at as well with the broader energy-band that was adopted recently. This dichotomy may be understood in a scenario where part of the accreting material forms a corona on the way of falling onto the surface of the neutron star during the decaying part of the low hard state. Based on the accretion rates of the preceding LHS, estimated from type-I X-ray bursts and from persistent emission, at least for IGR J17473-2721, most of the accretion material may fall on the neutron star (NS) surface in the LHS. Considering the burst behavior in the context of the outburst indicates a corona formed on top of the disk rather than on the NS surface.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Chinese-French SVOM mission for GRBs studies

We present the SVOM (Space-based multi-band astronomical Variable Objects Monitor) mission that the Chinese National Space Agency and the French Space Agency have decided to jointly implement. SVOM has been designed to detect all known types of gamma-ray bursts (GRBs), to provide fast and reliable GRB positions, to measure the broadband spectral shape and temporal properties of the GRB prompt emission, and to quickly identify the optical/near-infrared afterglows of detected GRBs, including high-redshift ones. Scheduled to be in orbit in the second half of the present decade, the SVOM satellite will carry a very innovative scientific payload combining for the first time a wide field X- and gamma-ray coded mask imager for GRB real-time localizations to few arcmin, a non-imaging gamma-ray monitor, and two narrow-field instruments for the study of the GRB early afterglow emission in the X-ray and visible bands. The SVOM payload is complemented by ground-based instruments including a wide-field camera to catch the GRB prompt emission in the visible band and two robotic telescopes to measure the photometric properties of the early afterglow. A particular attention is paid to the GRB follow-up in facilitating the observation of the SVOM detected GRB by the largest ground based telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

XMM-Newton Observations of NGC 247: X-ray Population and a Supersoft Ultraluminous X-ray Source

We report on a new XMM-Newton observation of NGC 247 from December 2009. The galaxy contains a supersoft, ultraluminous X-ray source (ULX) whose spectrum consists of a thermal component with a temperature about 0.1 keV and a power-law tail with a photon index around 2.5. The thermal emission is absolutely the dominant component, contributing 96% of the total luminosity in the 0.3-10 keV band. Variability is detected at timescales of 10^2 s and longer with a ν^-1 power spectrum. These properties are consistent with black hole binaries in the thermal state and suggest the presence of an intermediate mass black hole of at least 600 solar masses. However, the integrated rms power is much higher than typically found in the thermal state. An alternative explanation of the emission could be a photosphere with a radius about 10^9 cm. A possible absorption feature around 1 keV is detected, which may be due to absorption of highly ionized winds. X-ray sources within the disk of NGC 247 have a luminosity function consistent with that found in low mass X-ray binaries. We confirm previous results that X-rays from the quasar PHL 6625 may be absorbed by gas in NGC 247, mainly at energies below 0.3 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A comparison of hard X-ray photon indices and iron \ka emission lines in X-ray luminous narrow- and broad-line Seyfert 1 galaxies

We use publicly available XMM-Newton data to systematically compare the hard X-ray photon indices, $Γ_{\rm 2-10\ keV}$ and the iron K$α$ emission lines of narrow-line (NL) and broad-line Seyfert 1 (BLS1) galaxies. We compile a flux-limited ($f_{\rm 2-10\ keV} \geq 1 \times 10^{-12}$ erg s$^{-1}$ cm$^{-2}$) sample including 114 radio-quiet objects, with the 2-10 keV luminosity ranging from 10$^{41}$ to 10$^{45}$ erg s$^{-1}$. Our main results are: 1) NLS1s and BLS1s show similar luminosity distributions; 2) The weighted mean of $Γ_{\rm 2-10\ keV}$ of NLS1s, BLS1s and the total sample is $2.04\pm0.04$, $1.74\pm0.02$, $1.84\pm0.02$, respectively; a significant anti-correlation between \ga and FWHMH$β$ suggests that $Γ_{\rm 2-10\ keV} > 2.0$ may be taken to indicate X-ray luminous NLS1 type; 3) The 6.4 keV narrow iron K$α$ lines from NLS1s are generally weaker than that from BLS1s; this would indicate a smaller covering factor of the dusty tori in NLS1s, if the line emission originates from the inner boundary region of the dusty torus in an AGN; 4) all the broadened iron K$α$ lines with intrinsic width $σ>0.5$ keV correspond to FWHM\hb $\leq 4000 ~\kms$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Calibrating the correlation between black hole mass and X-ray variability amplitude: X-ray only black hole mass estimates for active galactic nuclei and ultra-luminous X-ray sources

A calibration is made for the correlation between the X-ray Variability Amplitude (XVA) and Black Hole (BH) mass. The correlation for 21 reverberation-mapped Active Galactic Nuclei (AGN) appears very tight, with an intrinsic dispersion of 0.20 dex. The intrinsic dispersion of 0.27 dex can be obtained if BH masses are estimated from the stellar velocity dispersions. We further test the uncertainties of mass estimates from XVAs for objects which have been observed multiple times with good enough data quality. The results show that the XVAs derived from multiple observations change by a factor of 3. This means that BH mass uncertainty from a single observation is slightly worse than either reverberation-mapping or stellar velocity dispersion measurements; however BH mass estimates with X-ray data only can be more accurate if the mean XVA value from more observations is used. Applying this relation, the BH mass of RE J1034+396 is found to be $4^{+3}_{-2} \times 10^6$ $M_{\odot}$. The high end of the mass range follows the relationship between the 2$f_0$ frequencies of high-frequency QPO and the BH masses derived from the Galactic X-ray binaries. We also calculate the high-frequency constant $C= 2.37 M_\odot$ Hz$^{-1}$ from 21 reverberation-mapped AGN. As suggested by Gierliński et al., $M_{\rm BH}=C/C_{\rm M}$, where $C_{\rm M}$ is the high-frequency variability derived from XVA. Given the similar shape of power-law dominated X-ray spectra in ULXs and AGN, this can be applied to BH mass estimates of ULXs. We discuss the observed QPO frequencies and BH mass estimates in the Ultra-Luminous X-ray source M82 X-1 and NGC 5408 X-1 and favor ULXs as intermediate mass BH systems (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

INTEGRAL and Swift/XRT observations of the source PKS 0208-512

The active galaxy PKS 0208-512, detected at lower energies by COMPTEL, has been claimed to be a MeV blazar from EGRET. We report on the most recent INTEGRAL observations of the blazar PKS 0208-512, which are supplemented by Swift ToO observations. The high energy X-ray and gamma-ray emission of PKS 0208-512 during August - December 2008 has been studied using 682 ks of INTEGRAL guest observer time and ~ 56 ks of Swift/XRT observations. These data were collected during the decay of a gamma-ray flare observed by Fermi/LAT. At X-ray energies (0.2 - 10 keV) PKS 0208-512 is significantly detected by Swift/XRT, showing a power-law spectrum with a photon index of ~ 1.64. Its X-ray luminosity varied by roughly 30% during one month. At hard X-/soft gamma-ray energies PKS 0208-512 shows a marginally significant (~ 3.2 sigma) emission in the 0.5-1 MeV band when combining all INTEGRAL/SPI data. Non-detections at energies below and above this band by INTEGRAL/SPI may indicate intrinsic excess emission. If this possible excess is produced by the blazar, one possible explanation could be that its jet consists of an abundant electron-positron plasma, which may lead to the emission of an annihilation radiation feature. Assuming this scenario, we estimate physical parameters of the jet of PKS 0208-512.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Low-frequency oscillations in XTE J1550-564

We present the results of timing analysis of the low-frequency Quasi-Periodic Oscillation (QPO) in the Rossi X-Ray Timing Explorer data of the black hole binary XTE J1550--564 during its 1998 outburst. The QPO frequency is observed to vary on timescales between $\sim$100 s and days, correlated with the count rate contribution from the optically thick accretion disk: we studied this correlation and discuss its influence on the QPO width. In all observations, the quality factors ($ν_0$/FWHM) of the fundamental and second harmonic peaks were observed to be consistent, suggesting that the quasi-periodic nature of the oscillation is due to frequency modulation. In addition to the QPO and its harmonic peaks, a new 1.5$ν$ component was detected in the power spectra. This component is broad, with a quality factor of $\sim$0.6. From this, we argue what the peak observed at half the QPO frequency, usually referred to as "sub-harmonic" could be the fundamental frequency, leading to the sequence 1:2:3:4. We also studied the energy dependence of the timing features and conclude that the two continuum components observed in the power spectrum, although both more intense at high energies, show a different dependence on energy. At low energies, the lowest-frequency component dominates, while at high energies the higher-frequency one has a higher fractional rms. An interplay between these two components was also observed as a function of their characteristic frequency. In this source, the transition between low/hard state and hard-intermediate state appears to be a smooth process.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the Solution to the "Frozen Star" Paradox, Nature of Astrophysical Black Holes, non-Existence of Gravitational Singularity in the Physical Universe and Applicability of the Birkhoff's Theorem

Oppenheimer and Snyder found in 1939 that gravitational collapse in vacuum produces a "frozen star", i.e., the collapsing matter only asymptotically approaches the gravitational radius (event horizon) of the mass, but never crosses it within a finite time for an external observer. Based upon our recent publication on the problem of gravitational collapse in the physical universe for an external observer, the following results are reported here: (1) Matter can indeed fall across the event horizon within a finite time and thus BHs, rather than "frozen stars", are formed in gravitational collapse in the physical universe. (2) Matter fallen into an astrophysical black hole can never arrive at the exact center; the exact interior distribution of matter depends upon the history of the collapse process. Therefore gravitational singularity does not exist in the physical universe. (3) The metric at any radius is determined by the global distribution of matter, i.e., not only by the matter inside the given radius, even in a spherically symmetric and pressureless gravitational system. This is qualitatively different from the Newtonian gravity and the common (mis)understanding of the Birkhoff's Theorem. This result does not contract the "Lemaitre-Tolman-Bondi" solution for an external observer.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations on High-z Long Gamma-Ray Burst Rate

Since the launch of Swift satellite, the detections of high-z (z>4) long gamma-ray bursts (LGRBs) have been rapidly growing, even approaching the very early Universe (the record holder currently is z=8.3). The observed high-z LGRB rate shows significant excess over that estimated from the star formation history. We investigate what may be responsible for this high productivity of GRBs at high-z through Monte Carlo simulations, with effective Swif/BAT trigger and redshift detection probabilities based on current Swift/BAT sample and CGRO/BATSE LGRB sample. We compare our simulations to the Swift observations via log N-log P, peak luminosity (L) and redshift distributions. In the case that LGRB rate is purely proportional to the star formation rate (SFR), our simulations poorly reproduce the LGRB rate at z>4, although the simulated log N-log P distribution is in good agreement with the observed one. Assuming that the excess of high-z GRB rate is due to the cosmic metallicity evolution or unknown LGRB rate increase parameterized as (1+z)^delta, we find that although the two scenarios alone can improve the consistency between our simulations and observations, incorporation of them gives much better consistency. We get 0.2<epsilon<0.6 and delta<0.6, where epsilon is the metallicity threshold for the production of LGRBs. The best consistency is obtained with a parameter set (epsilon, delta)=(~0.4, ~0.4), and BAT might trigger a few LGRBs at z~14. With increasing detections of GRBs at z>4 (~15% of GRBs in current Swift LGRB sample based on our simulations), a window for very early Universe is opening by Swift and up-coming SVOM missions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Studying Frequency Relations of kHz QPOs for 4U 1636-53 and Sco X-1: Observations Confront Theories

By fitting the frequencies of simultaneous lower and upper kilohertz Quasi-Periodic Oscillations (kHz QPOs) in two prototype neutron star QPO sources (4U~1636-536 and Sco X-1), we test the predictive power of all currently proposed QPO models. Models predict either a linear, power-law or any other relationship between the two frequencies. We found that for plausible neutron star parameters (mass and angular momentum), no model can satisfactorily reproduce the data, leading to very large chi-squared values in our fittings. Both for 4U~1636-53 and Sco X-1, this is largely due to the fact that the data significantly differ from a linear relationship. Some models perform relatively better but still have their own problems. Such a detailed comparison of data with models shall enable to identify routes for improving those models further.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Significant foreground unrelated non-acoustic anisotropy on the one degree scale in WMAP 5-year observations

The spectral variation of the cosmic microwave background (CMB) as observed by WMAP was tested using foreground reduced WMAP5 data, by producing subtraction maps at the 1$^\circ$ angular resolution between the two cosmological bands of V and W, for masked sky areas that avoid the Galactic disk. The resulting $V-W$ map revealed a non-acoustic signal over and above the WMAP5 pixel noise, with two main properties. Firstly, it possesses quadrupole power at the $\approx$ 1 $μK$ level which may be attributed to foreground residuals. Second, it fluctuates also at all values of $\ell >$ 2, especially on the $1^\circ$ scale ($200 \lesssim \ell \lesssim 300$). The behavior is {\it random and symmetrical} about zero temperature with a r.m.s. amplitude of $\approx$ 7 $μK$, or 10 % of the maximum CMB anisotropy, which would require a `cosmic conspiracy' among the foreground components if it is a consequence of their existences. Both anomalies must be properly diagnosed and corrected if `precision cosmology' is the claim. The second anomaly is, however, more interesting because it opens the question on whether the CMB anisotropy genuinely represents primordial density seeds.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Jets and environment of microquasars

Two relativistic X-ray jets have been detected with the Chandra X-ray observatory from the black hole X-ray transient XTE J1550-564. We report a full analysis of the evolution of the two jets with a gamma-ray burst external shock model. A plausible scenario suggests a cavity outside the central source and the jets first travelled with constant velocity and then are slowed down by the interactions between the jets and the interstellar medium (ISM). The best fitted radius of the cavity is $\sim$0.36 pc on the eastern side and $\sim$0.46 pc on the western side, and the densities also show asymmetry, of $\sim$0.015 cm$^{-3}$ on the east to $\sim$0.21 cm$^{-3}$ on the west. A large scale low density region is also found in another microquasar system, H 1743-322. These results are consistent with previous suggestions that the environment of microquasars should be rather vacuous, compared to the normal Galactic environment. A generic scenario for microquasar jets is proposed, classifying the observed jets into three main categories, with different jet morphologies (and sizes) corresponding to different scales of vacuous environments surrounding them.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Shot cutoff timescales in different spectral states of Cygnus X-1

We investigate the shot cutoff timescale evolutions during different spectral states in the black hole binary Cygnus X-1 with the recently proposed w spectral analysis technique. For low energy shots, their cutoff timescale decreases from the low state to intermediate state, and to high state monotonically. However in the high energy range, where shots are believed to be reprocessed in the hot corona, timescales are almost the same for different spectral states. A linear correlation is found between the energy spectrum photon index and the shot cutoff timescale at low energies. Both narrow and broad iron Kalpha lines are distinguished from timing analysis, providing dynamical evidence that broad and narrow iron lines originate respectively at small and large radii of the accretion disk. Possible mechanisms for state transitions and accretion flow geometry in accreting stellar mass black holes are discussed.

0 citations0 reviewsNo code linkedOpen access