Millisecond Radio Pulsars Research Articles

Timing of Millisecond Pulsars in NGC 6752. III. On the Presence of Nonluminous Matter in the Cluster’s Core

Millisecond pulsars are subject to accelerations in globular clusters (GCs) that manifest themselves in both the first and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of five millisecond radio pulsars in the core of the core-collapse GC NGC 6752 with the Parkes (Murriyang) and MeerKAT radio telescopes, which have allowed us to measure the proper motions, positions, and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56 × 103 M ⊙ is present in the core of NGC 6752. On the other hand, our studies highly disfavor the presence of an intermediate-mass black hole at the center of the cluster, with a mass equal to or greater than ∼3000M ⊙.

High-temporal-resolution optical spectroscopic observations of the transitional millisecond pulsar PSR J1023+0038

Transitional millisecond pulsars (tMSPs) represent a dynamic category of celestial sources that establish a crucial connection between low-mass X-ray binaries and millisecond radio pulsars. These systems exhibit transitions from rotation-powered states to accretion-powered ones and vice versa, highlighting the tight evolutionary link expected by the so-called recycling scenario. In their active phase, these sources manifest two distinct emission modes named high and low, occasionally punctuated by sporadic flares. Here, we present high-time-resolution spectroscopic observations of the binary tMSP J1023+0038, in the sub-luminous disc state. This is the first short-timescale (sim 1 min) optical spectroscopic campaign ever conducted on a tMSP. The campaign was carried out over the night of June 10, 2021 using the Gran Telescopio Canarias. The optical continuum shows erratic variability, without clear evidence of high and low modes or of orbital modulation. Besides, the analysis of these high-temporal-cadence spectroscopic observations reveals, for the first time, evidence for a significant (up to a factor of $ 2$) variability in the emission line properties (equivalent width and full width half maximum) over a timescale of minutes. Intriguingly, the variability episodes observed in the optical continuum and in the emission line properties seem uncorrelated, making their origin unclear.

A targeted radio pulsar survey of redback candidates with MeerKAT

ABSTRACT Redbacks are millisecond pulsar binaries with low-mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localization regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838−2827, J0955−3947, and J2333−5526, confirming their redback nature. PSR J0838−2827 remained undetected for 2 yr after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-yr timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333−5526 may have a neutron star mass in excess of 2 M⊙.

On the Spin Period Distribution of Millisecond Pulsars

Spin period distribution provides important clues to understand the formation of millisecond pulsars (MSPs). To uncover the intrinsic period distribution, we analyze three samples of radio MSPs in the Galactic field and globular clusters. The selection bias due to pulse broadening has been corrected but turns out to be negligible. We find that all the samples can be well described by a Weibull distribution of spin frequencies. Considering MSPs in the Galactic field or globular clusters and in isolation or binary systems, we find no significant difference in the spin distribution among these subpopulations. Based on the current known population of MSPs, we find that submillisecond pulsars are unlikely to be discovered by the Square Kilometre Array, although up to ∼10 discoveries of pulsars that spin faster than the current record holder of P = 1.4 ms are expected.

A Possible γ-Ray Pulsation from PSR J1740−5340B in the Globular Cluster NGC 6397

Recently, a new radio millisecond pulsar (MSP) J1740−5340B, hosted in the globular cluster (GC) NGC 6397, was reported with a 5.78 ms spin period in an eclipsing binary system with a 1.97 days orbital period. Based on a modified radio ephemeris updated by tool tempo2, we analyze the ∼15 yr γ-ray data obtained from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope and detect PSR J1740−5340B’s γ-ray pulsation at a confidence level of ∼4σ with a weighted H-test value of ∼26. By performing a phase-resolved analysis, the γ-ray luminosity in on-pulse interval of PSR J1740−5340B is L γ ∼ 3.8 × 1033 erg s−1 using NGC 6397's distance of 2.48 kpc. And γ-rays from the on-pulse part of PSR J1740−5340B contribute ∼90% of the total observed γ-ray emissions from NGC 6397. No significant γ-ray pulsation of another MSP J1740−5340A in the GC is detected. Considering that the previous four cases of MSPs in GCs, more data in γ-ray, X-ray, and radio are encouraged to finally confirm the γ-ray emissions from MSP J1740−5340B, especially starving for a precise ephemeris.

Properties of PSR J1023+0038 Based on Kepler, TESS, and FAST

The pulsar system PSR J1023+0038 is a very interesting binary with radio millisecond pulsar and low-mass X-ray binary transitions, which depend on the stellar accretion process. We have used Kepler light curves with 1 minute cadence and TESS light curves with 2 minute cadence to study the flare properties of PSR J1023+0038. We detected 516 flare events in the K2 data and 84 flare events in the TESS survey data. We determined the flare rise and delay times, the flare durations and amplitudes, and flare energies. We obtained a value of the power-law index of 0.359 ± 0.023 between the flare energies and decay times, which is similar to the theoretically predicted value of 1/3. This in turn indicates that PSR J1023+0038 might be exhibiting physical mechanisms described by magnetic reconnection theory. Using the maximum likelihood estimation method, we also calculated the power-law index of the cumulative flare frequency distribution, finding 1.87 ± 0.27 for the Kepler data and 1.74 ± 0.29 for the TESS data. Meanwhile, we reduced six FAST observations of PSR J1023+0038 from 2019–2021 using a standard pulsar search procedure. We report that we detected a prompt signal with a period compatible with that found in previous observations of PSR J1023+0038. However, there are several reasons (statistical, excess DM, pulse shape, etc.) as to why the association is not likely.

Observational Evidence for a Spin-up Line in the P– Diagram of Millisecond Pulsars

It is believed that millisecond pulsars attain their fast spins by accreting matter and angular momentum from companion stars. Theoretical modeling of the accretion process suggests a spin-up line in the period–period derivative (P-) diagram of millisecond pulsars, which plays an important role in population studies of radio millisecond pulsars and accreting neutron stars in X-ray binaries. Here we present observational evidence for such a spin-up line using a sample of 143 radio pulsars with P < 30 ms. We also find that PSRs J1823−3021A and J1824−2452A, located near the classic spin-up line, are consistent with the broad population of millisecond pulsars. Finally, we show that our approach of Bayesian inference can probe accretion physics, allowing constraints to be placed on the accretion rate and the disk–magnetosphere interaction.

Millisecond pulsar kicks cause difficulties in explaining the Galactic Centre gamma-ray excess

ABSTRACT The unexplained excess gamma-ray emission from the Milky Way’s Galactic Centre has puzzled astronomers for nearly a decade. Two theories strive to explain the origin of this excess: self-annihilating dark matter particles or an unresolved population of radio millisecond pulsars. We examine the plausibility of a pulsar origin for the GeV excess using N-body simulations. We simulated millisecond pulsars in a realistic dynamical environment: (i) pulsars were born from the known stellar mass components of our Galaxy; (ii) pulsars were given natal velocity kicks as empirically observed from two different studies (or, for comparison, without kicks); (iii) pulsars were evolved in a Galactic gravitational potential consistent with observations. Multiple populations of pulsars (with different velocity kicks) were simulated over 1 Gyr. With final spatial distributions of pulsars, we constructed synthetic gamma-ray surface brightness profiles. From comparisons with published Fermi-LAT surface brightness profiles, our pulsar simulations cannot reproduce the concentrated emission in the central degrees of the Bulge, though models without natal velocity kicks approach the data. We considered additive combinations of our (primordial MSP) simulations with models where pulsars are deposited from destroyed globular clusters in the Bulge, and a simple model for pulsars produced in the nuclear star cluster. We can reasonably reproduce the measured central gamma-ray surface brightness distribution of Horiuchi and collaborators using several combinations of these models, but we cannot reproduce the measured distribution of Di Mauro with any combination of models. Our fits provide constraints on potential pathways to explain the gamma-ray excess using MSPs.

On the recurrence times of neutron star X-ray binary transients and the nature of the Galactic Centre quiescent X-ray binaries

ABSTRACT The presence of some X-ray sources in the Galactic Centre region which show variability, but do not show outbursts in over a decade of monitoring has been used to argue for the presence of a large population of stellar mass black holes in this region. A core element of the arguments that these objects are accreting black holes is the claim that neutron stars (NSs) in low-mass X-ray binaries (LMXBs) do not have long transient recurrence times. We demonstrate in this paper that about half of the known transient LMXBs with clear signatures for NS primaries have recurrence times in excess of a decade for outbursts at the sensitivity of MAXI. We furthermore show that, in order to reconcile the expected total population of NS LMXBs with the observed one and with the millisecond radio pulsar (MSRP) population of the Galaxy, systems with recurrence times well in excess of a century for outbursts detectable by instruments like MAXI must be the dominant population of NS LMXBs, and that few of these systems have yet been discovered.

Gamma-Ray Spectral Properties of the Galactic Globular Clusters: Constraint on the Number of Millisecond Pulsars

Abstract We study the γ-ray spectra of 30 globular clusters (GCs) thus far detected with the Fermi Gamma-ray Space Telescope. Presuming that γ-ray emission of a GC comes from millisecond pulsars (MSPs) contained within the GC, a model that generates spectra for the GCs is built based on the γ-ray properties of the detected MSP sample. We fit the GCs’ spectra with the model, and for 27 of them, their emission can be explained as arising from MSPs. The spectra of the other three, NGC 7078, 2MS-GC01, and Terzan 1, cannot be fit with our model, indicating that MSPs’ emission should not be the dominant one in the first two and the third one has a unique hard spectrum. We also investigate six nearby GCs that have relatively high encounter rates compared to the comparison cases. The candidate spectrum of NGC 6656 can be fit with that of one MSP, supporting its possible association with the γ-ray source at its position. The five others do not have detectable γ-ray emission. Their spectral upper limits set limits of ≤1 MSPs in them, consistent with the numbers of radio MSPs found in them. The estimated numbers of MSPs in the γ-ray GCs generally match well those reported for radio pulsars. Our studies of the γ-ray GCs and the comparison nearby GCs indicate that the encounter rate should not be the only factor determining the number of γ-ray MSPs a GC contains.

Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: inferences in the stochastic gravitational-wave background search

ABSTRACT We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 yr. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, ENTERPRISE, and TEMPONEST+FORTYTWO, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $\gamma = 3.78_{-0.59}^{+0.69}$ (95 per cent credible regions). Fixing the spectral index to γ = 13/3 as expected from the GWB by circular, inspiralling supermassive black hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS, and EPHEMGP, to address possible Solar system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension.

Combined analysis of neutron star natal kicks using proper motions and parallax measurements for radio pulsars and Be X-ray binaries

ABSTRACT Supernova explosion and the associated neutron star (NS) natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary, or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain NS natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass, and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries (BeXs). Our best model for natal kicks is consistent with both measurements and contains a fraction of w = 0.2 ± 0.1 weak natal kicks with $\sigma _1 = 45^{+25}_{-15}$ km s−1, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: σ2 = 336 km s−1. We found that currently used models for natal kicks of NSs produced by electron capture supernova (ecSN; combination of Maxwellian σ = 265 km s−1 and σ = 30 km s−1 for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and BeXs.

Formation of millisecond pulsars with helium white dwarfs, ultra-compact X-ray binaries, and gravitational wave sources

ABSTRACT Close-orbit low-mass X-ray binaries (LMXBs), radio binary millisecond pulsars (BMSPs) with extremely low-mass helium white dwarfs (ELM He WDs) and ultra-compact X-ray binaries (UCXBs) are all part of the same evolutionary sequence. It is therefore of uttermost importance to understand how these populations evolve from one specie to another. Moreover, UCXBs are important gravitational wave (GW) sources and can be detected by future space-borne GW observatories. However, the formation and evolutionary link between these three different populations of neutron star (NS) binaries are not fully understood. In particular, a peculiar fine-tuning problem has previously been demonstrated for the formation of these systems. In this investigation, we test a newly suggested magnetic braking prescription and model the formation and evolution of LMXBs. We compute a grid of binary evolution models and present the initial parameter space of the progenitor binaries which successfully evolve all the way to produce UCXBs. We find that the initial orbital period range of LMXBs, which evolve into detached NS + ELM He WD binaries and later UCXBs, becomes significantly wider compared to evolution with a standard magnetic braking prescription, and thus helps to relieve the fine-tuning problem. However, we also find that formation of wide-orbit BMSPs is prohibited for strong versions of this new magnetic braking prescription, which therefore calls for a revision of the prescription. Finally, we present examples of the properties of UCXBs as Galactic GW sources and discuss their detection by the LISA, TianQin, and Taiji observatories.

Constraining the ellipticity of millisecond pulsars with observed spin-down rates

A spinning neutron star (NS) that is asymmetric with respect to its spin axis can emit continuous gravitational wave (GW) signals. The spin frequencies and their distribution of radio millisecond pulsars (MSPs) and accreting MSPs provide some evidences of GW radiation, and MSPs are ideal probes detecting high frequency GW signals. It is generally thought that MSPs originate from the recycled process, in which the NS accretes the material and angular momentum from the donor star. The accreted matter would be confined at the polar cap zone by an equatorial belt of compressed magnetic field fixed in the deep crust of the NS, and yields "magnetic mountain". Based on an assumption that the spin-down rates of three transitional MSPs including PSR J1023+0038 are the combinational contribution of the accretion torque, the propeller torque, and the GW radiation torque, in this work we attempt to constrain the ellipticities of MSPs with observed spin-down rates. Assuming some canonical parameters of NSs, the ellipticities of three transitional MSPs and ten redbacks are estimated to be $\epsilon=(0.9-23.4)\times 10^{-9}$. The electrical resistivities of three transitional MSPs are also derived to be in the range $\eta=(1.2-15.3)\times 10^{-31}~\rm s$, which display an ideal power law relation with the accretion rate. The characteristic strains ($h_{\rm c}=(0.6-2.5)\times10^{-27}$) of GW signals emitting by these sources are obviously beyond the sensitivity scope of the aLIGO. We expect that the third-generation GW detectors like the Einstein Telescope can seize the GW signals from these sources in the future.

A method of ground target positioning by observing radio pulsars

We use radio millisecond pulsars (MSPs) to determine the position of the Parkes telescope to explore the feasibility of pulsar navigation for terrestrial application, as well to investigate the relations between pulsar observations and the positional errors for potential X-ray pulsar navigation (XNAV). Different from the analytical algorithm (Han et al. Astrophys. Space Sci. 364, 3, 2018) which derives the relations between observatory coordinates and pulsar timing residuals, a grid search method is used in our work. The method divides the space adjacent to the radio telescope into small grids, from which we form the rms timing residuals and fit pulsar parameters to find the position with minimal χ2 value. Six pulsars from the Parkes Pulsar Timing Array (PPTA) are selected to perform position determinations. The positional errors are analyzed together with the pulsar reference phase, timing residuals, and observation numbers. Using a hybrid database of pulsar timing and cartographic data, we obtained best positional accuracy of around 10 m by PSR J0437 − 4715. To investigate the feasibility of pulsar terrestrial navigation, only few observations from PSR J0437 − 4715 are used for position determinations. The results show that the precision of about 100 m is achievable with only 3 or 4 observations. In this paper, we introduce our method in detail and discuss the positioning performance in different scenarios.

Radio and X-ray monitoring of the accreting millisecond X-ray pulsar IGR J17591−2342 in outburst

ABSTRACT IGR J17591−2342 is a new accreting millisecond X-ray pulsar that was recently discovered in outburst in 2018. Early observations revealed that the source’s radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS–LMXB) at comparable X-ray luminosity, and assuming its likely ≳6 kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591−2342. In total, we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift–XRT) observations, which make IGR J17591−2342 one of the best-sampled NS–LMXBs. We use these to fit a power-law correlation index $\beta = 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities (LR ∝ LXβ). However, our monitoring revealed a large scatter in IGR J17591−2342’s radio luminosity (at a similar X-ray luminosity, LX ∼1036 erg s−1, and spectral state), with LR ∼ 4 × 1029 erg s−1 during the first three reported observations, and up to a factor of 4 lower LR during later radio observations. None the less, the average radio luminosity of IGR J17591−2342 is still one of the highest among NS–LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591−2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. None the less, we cannot rule out that IGR J17591−2342 becomes a radio millisecond pulsar during quiescence.

Millisecond Pulsars and Black Holes in Globular Clusters

Abstract Over 100 millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic N-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses ≈2 × 105 M ⊙ can produce up to 10–20 MSPs, while a very massive GC model with mass ≈106 M ⊙ can produce close to 100. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated toward the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) are less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS–NS and NS–BH binaries in GCs, although the predicted numbers are very small.

Modeling Pulsars in dense star clusters

AbstractOver a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic N-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses ≍ 2 × 105M⊙ can produce up to 10 – 20 MSPs, while a very massive GC model with mass ≍ 106M⊙ can produce close to 100. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations.

Investigation of the Millisecond Pulsar Origins by their Spin Periods at the Wavebands of Radio, X-Ray, and γ-Ray

To track the formation and evolution links of the millisecond pulsars (MSPs) powered by accretion and rotation in the Galactic field, we investigate the spin period (P) and spin-down power () distributions of the MSPs observed at the wavebands of radio, X-ray, and γ-ray. We find that all but one (119/120) of the γ-ray MSPs have been detected with the radio signals (radio+γ MSPs); on the contrary, nearly half of the radio MSPs (118/237) have not been detected with γ-rays (radio-only MSPs). In addition, the radio+γ MSPs are shown to be the relative faster and more energetic objects ( and ) compared with the radio-only MSPs ( and ), while the spin periods of these two MSP populations are compatible with the log-normal distributions by the statistical tests. Most rotation-powered MSPs (RMSPs) with the radio eclipsing (31/34) exhibit the radio+γ signals, which share the faster spin () and larger spin-down power () distributions than the non-eclipsing ones (, ), implying the radio+γ MSPs to be younger than the radio-only MSPs. It is noticed that the spin distribution of the accretion-powered X-ray MSPs shows a clustering phenomenon around ∼1.6–2.0 ms, which is not observed in RMSPs, hinting that the RMSPs may experience the multiple possible origins. Particularly, all the three super-fast spinning RMSPs with P ∼ 1.4–1.6 ms exhibit the non-eclipsing, and we argue that they may be the distinctive sources formed by the accretion induced collapse (AIC) of white dwarfs.

Optical monitoring of a transitional millisecond pulsar: PSR J1023+0038

PSR J10203+0038 is a transitional millisecond pulsar (tMSP) in an eclipsing binary system, which has been observed to switch between the radio loud millisecond pulsar (MSP) and low-mass-X-ray binary (LMXB) states. This behavior offers a great opportunity to study the origin of MSPs and confirming the ‘recycling’ scenario, a theoretical model of MSP’s origin. We develop an automated pipeline to monitor the system using Python programming language and Source-Extractor software for detecting the objects and measuring its magnitude. We obtain a series of observations with the 0.6m PROMPT-8 telescope at Cerro Tololo in Chile. The magnitude threshold for alert has been set 16.884 mag in R filter. When the magnitude of the system increases over the limit, 16.884 mag in filter R, the pipeline will alert us about the next possible switching of this system. The pipeline has been running on server at National Astronomical Research of Thailand (NARIT) since January 2018. We have found that, during January and February 2018, the system still remains in LMXB state.