Search

Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year $\textit{Fermi}$-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$\sigma$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed $\textit{in toto}$, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component--which we uncover at a significance of 8.2$\sigma$--has a power index of 2.79 $\pm$ 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation., Comment: 16+9 pages, 13+4 figures, 4+1 tables. match the published version in MNRAS

The multi-band observation, from radio to gamma-ray, plays a crucial role in constraining the emission physics and geometrical magnetospheric models of the pulsars. In this paper, We use three-dimensional magnetospheric model, with the joint of the annular gap (AG) and core gap (CG) model, to simulate the radio and gamma-ray light curves for two representative pulsars, a millisecond pulsar J2302+4442 and a young pulsar J0659+1414. It is found that for PSR J2302 + 4442, the simulated magnetic inclination angle α and viewing angle ζ are 34 ∘ and 46 ∘ , and for PSR J0659 + 1414 are 45 ∘ and 33 ∘ . This implies that the radio and gamma-ray pulse emission originates from different magnetospheric regions and the radiation geometry for millisecond and young pulsar is different. The studies show that the joint of the AG and CG model can describe the light curves not only for millisecond pulsars but also for young pulsars, and hence it can be used to explain the multi-wavelength observations of pulsars.

Accreting white dwarfs are often found in close binary systems with orbital periods ranging from tens of minutes to several hours. In most cases, the accretion process is relatively steady, with significant modulations only occurring on time-scales of ~days or longer. Here, we report the discovery of abrupt drops in the optical luminosity of the accreting white dwarf binary system TW Pictoris by factors up to 3.5 on time-scales as short as 30 minutes. The optical light curve of this binary system obtained by the Transiting Exoplanet Survey Satellite (TESS) clearly displays fast switches between two distinct intensity modes that likely track the changing mass accretion rate onto the white dwarf. In the low mode, the system also displays magnetically-gated accretion bursts, which implies that a weak magnetic field of the white dwarf truncates the inner disk at the co-rotation radius in this mode. The properties of the mode switching observed in TW Pictoris appear analogous to those observed in transitional millisecond pulsars, where similar transitions occur, although on timescales of ~tens of seconds. Our discovery establishes a previously unrecognised phenomenon in accreting white dwarfs and suggests a tight link to the physics governing magnetic accretion onto neutron stars., Comment: Submitted on 17 June 2021. Accepted for publication in Nature Astronomy on 18 August 2021

Extremely low-mass white dwarfs (ELM WDs) are the result of binary evolution in which a low-mass donor star is stripped by its companion leaving behind a helium-core white dwarf. We explore the formation of ELM WDs in binary systems considering the Convection And Rotation Boosted magnetic braking treatment. Our evolutionary sequences were calculated using the MESA code, with initial masses of 1.0 and 1.2 Msun (donor), and 1.4 (accretor), compatible with low mass X-ray binaries (LMXB) systems. We obtain ELM models in the range 0.15 to 0.27 Msun from a broad range of initial orbital periods, 1 to 25 d. The bifurcation period, where the initial period is equal to the final period, ranges from 20 to 25 days. In addition to LMXBs, we show that ultra-compact X-ray binaries (UCXB) and wide-orbit binary millisecond pulsars can also be formed. The relation between mass and orbital period obtained is compatible with the observational data from He white dwarf companions to pulsars., Comment: 16 pages, 6 figures. Accepted for publication in MNRAS

Transitional millisecond pulsars are millisecond pulsars that switch between a rotation-powered millisecond pulsar state and an accretion-powered X-ray binary state, and are thought to be an evolutionary stage between neutron star low-mass X-ray binaries and millisecond pulsars. So far, only three confirmed systems have been identified in addition to a handful of candidates. We present the results of a multi-wavelength study of the low-mass X-ray binary NGC 6652B in the globular cluster NGC 6652, including simultaneous radio and X-ray observations taken by the Karl G. Jansky Very Large Array and the Chandra X-ray Observatory, and optical spectroscopy and photometry. This source is the second brightest X-ray source in NGC 6652 ($L_{\textrm{X}}\sim1.8\times10^{34}$ erg s$^{-1}$) and is known to be variable. We observe several X-ray flares over the duration of our X-ray observations, in addition to persistent radio emission and occasional radio flares. Simultaneous radio and X-ray data show no clear evidence of anti-correlated variability. Optical spectra of NGC 6652B indicate variable, broad H $\alpha$ emission which transitions from double-peaked emission to absorption over a time-scale of hours. We consider a variety of possible explanations for the source behaviour, and conclude that based on the radio and X-ray luminosities, short time-scale variability and X-ray flaring, and optical spectra, NGC 6652B is best explained as a transitional millisecond pulsar candidate that displays prolonged X-ray flaring behaviour. However, this could only be confirmed with observations of a change to the rotation-powered millisecond pulsar state., Comment: 15 pages, 8 figures, accepted for publication in MNRAS

R-modes are the oscillations observed especially in stars rotating at high speeds, where the dominant force is the Coriolis force. Since the dominant force in these modes is the Coriolis force, the frequency of the mode is independent of the internal structure of the star but is a function of its angular velocity. Gravitational radiation drives these modes towards the zone of instability in all rotating stars, while the viscous forces inside the star play a destructive role in this oscillation. If the star is spinning slowly, the time scale for gravitational radiation is always much larger and the star is stable. In fast rotating stars, the instability of the star against these oscillations is defined as a function of temperature, since the viscous forces depend on the temperature. The X-ray spectroscopy is ideal for limiting the amplitude of these oscillations because the damping of these oscillations causes an extra temperature increase in the star. Then the r-mode instability stops and the star slowly starts to cool again. This temperature rise in the star can be examined by X-ray spectroscopy. This instability is expected to occur in fast-spinning systems, for example, newborn neutron stars, or members of low-mass X-ray binaries (LMXBs), or millisecond pulsars. r-mode instability can affect any system. According to the observations of the LMXBs, many sources are located in the theoretically expected r-mode instability window.Since the heat transport and cooling characteristics in neutron stars or LMXBs are dependent on their internal structure, they provide us with important clues about the internal structure of stars. For instance; rotation frequency – the boundary of the r-mode zone of instability in the temperature plane, and especially the minimum that determines the final rotation frequency of the star, is very different for stars with different internal compositions. In this context, potential X-ray millisecond pulsars or LMXBs with high rotation frequency from the literature will be selected and studied within the scope of the thesis. By analyzing the X-ray data taken from the archive, it was attempted to set boundaries on the surface temperatures of these sources. As a result of the analyzes to be made, it is aimed to examine the r-mode oscillations and internal structure of the studied sources.

The rotation of a magnetised accreting neutron star (NS) in a binary system is described by its spin period and two angles: spin inclination $\alpha$ with respect to the orbital momentum and magnetic angle $\chi$ between the spin and the magnetic moment. Magnetospheric accretion spins the NS up and adjusts its rotation axis, decreasing $\alpha$ to nearly perfect alignment. Its effect upon the magnetic angle is more subtle and relatively unstudied. In this work, we model the magnetic angle evolution of a rigid spherical accreting NS. We find that the torque spinning the NS up may affect the magnetic angle while both $\alpha$ and $\chi$ significantly deviate from zero, and the spin-up torque varies with the phase of the spin period. As the rotation axis of the NS is being aligned with the spin-up torque, the magnetic axis becomes misaligned with the rotation axis. Under favourable conditions, magnetic angle may increase by $\Delta \chi \sim 15^\circ-20^\circ$. This orthogonalisation may be an important factor in the evolution of millisecond pulsars, as it partially compensates the $\chi$ decrease potentially caused by pulsar torques. If the direction of the spin-up torque changes randomly with time, as in wind-fed high-mass X-ray binaries, both the rotation axis of the NS and its magnetic axis become involved in a non-linear random-walk evolution. The ultimate attractor of this process is a bimodal distribution in $\chi$ peaking at $\chi =0^\circ$ and $\chi = 90^\circ$., Comment: 12 pages, 7 figures. Resubmitted to MNRAS after minor revision. Any feedback is very welcome!

We analyse 55 ks of Chandra X-ray observations of the Galactic globular cluster M13. Using the latest radio timing positions of six known millisecond pulsars (MSPs) in M13 from Wang et al. (2020), we detect confident X-ray counterparts to five of the six MSPs at X-ray luminosities of $L_X$(0.3-8 keV)$\sim 3 \times 10^{30} - 10^{31}~{\rm erg~s^{-1}}$, including the newly discovered PSR J1641+3627F. There are limited X-ray counts at the position of PSR J1641+3627A, for which we obtain an upper limit $L_X, 9 pages, 5 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Neutron stars receive velocity kicks at birth in supernovae. Those formed in electron-capture supernovae from super asymptotic giant branch stars -- the lowest mass stars to end their lives in supernovae -- may receive significantly lower kicks than typical neutron stars. Given that many massive stars are members of wide binaries, this suggests the existence of a population of low-mass ($1.25 < M_\mathrm{psr} / $M$_\odot < 1.3$), wide ($P_\mathrm{orb} \gtrsim 10^{4}$\,day), eccentric ($e \sim 0.7$), unrecycled ($P_\mathrm{spin} \sim 1$\,s) binary pulsars. The formation rate of such binaries is sensitive to the mass range of (effectively) single stars leading to electron capture supernovae, the amount of mass lost prior to the supernova, and the magnitude of any natal kick imparted on the neutron star. We estimate that one such binary pulsar should be observable in the Milky Way for every 10,000 isolated pulsars, assuming that the width of the mass range of single stars leading to electron-capture supernovae is $\lesssim 0.2$\,M$_\odot$, and that neutron stars formed in electron-capture supernovae receive typical kicks less than 10\,km s$^{-1}$. We have searched the catalog of observed binary pulsars, but find no convincing candidates that could be formed through this channel, consistent with this low predicted rate. Future observations with the Square Kilometre Array may detect this rare sub-class of binary pulsar and provide strong constraints on the properties of electron-capture supernovae and their progenitors., Accepted to MNRAS. 7 pages, 2 figures

While the majority of massive stars have a stellar companion, most pulsars appear to be isolated. Taken at face value, this suggests that most massive binaries break apart due to strong natal kicks received in supernova explosions. However, the observed binary fraction can still be subject to strong selection effects, as monitoring of newly discovered pulsars is rarely carried out for long enough to conclusively rule out multiplicity. Here, we use the second Gaia Data Release (DR2) to search for companions to 1534 rotation-powered pulsars with positions known to better than 0.5 arcseconds. We find 22 matches to known pulsars, including one not reported elsewhere, and 8 new possible companions to young pulsars. We examine the photometric and kinematic properties of these systems and provide empirical relations for identifying Gaia sources with potential millisecond pulsar companions. Our results confirm that the observed multiplicity fraction is small. However, we show that the number of binaries below the sensitivity of Gaia and radio timing in our sample could still be significantly higher. We constrain the binary fraction of young pulsars to be $f_{\rm young}^{\rm true}\leq 5.3(8.3)\%$ under realistic(conservative) assumptions for the binary properties and current sensitivity thresholds. For massive stars ($\geq 10$ M$_{\odot}$) in particular, we find $f_{\rm OB}^{\rm true}\leq 3.7\%$ which sets a firm independent upper limit on the galactic neutron-star merger rate, $\leq 7.2\times 10^{-4}$ yr$^{-1}$. Ongoing and future projects such as the CHIME/pulsar program, MeerTime, HIRAX and ultimately the SKA, will significantly improve these constraints in the future., Comment: Accepted for publication in MNRAS; comments are welcome

The wealth of detections of millisecond pulsars (MSPs) in $\gamma$-rays by {\em Fermi} has spurred searches for these objects among the several unidentified $\gamma$-ray sources. Interesting targets are a sub-class of binary MSPs, dubbed "Black Widows" (BWs) and "Redbacks" (RBs), which are in orbit with low-mass non-degenerate companions fully or partially ablated by irradiation from the MSP wind. These systems can be easily missed in radio pulsar surveys owing to the eclipse of the radio signal by the intra-binary plasma from the ablated companion star photosphere, making them better targets for multi-wavelength observations. We used optical and X-ray data from public databases to carry out a systematic investigation of all the unidentified $\gamma$-ray sources from the Fermi Large Area Telescope (LAT) Third Source Catalog (3FGL), which have been pre-selected as likely MSP candidates according to a machine-learning technique analysis. We tested our procedure by recovering known binary BW/RB identifications and searched for new ones, finding possible candidates. At the same time, we investigated previously proposed BW/RB identifications and we ruled out one of them based upon the updated $\gamma$-ray source coordinates., Comment: 22 pages, 14 figures, 21 pages - This manuscript is based on the MSc Thesis work of C. Braglia, defended on April 2020 at the University of Milan, accepted for publication in MNRAS

The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fast Folding Algorithm searches. We present 19 new discoveries, bringing the total to 27 discoveries in the programme, and we report the results of follow-up timing observations at Parkes for 26 of these pulsars, including the millisecond pulsar PSR J1421-4409; the faint, highly-modulated, slow pulsar PSR J1646-1910; and the nulling pulsar PSR J1337-4441. We present new timing solutions for 23 pulsars, and we report flux densities, modulation indices, and polarization properties., Comment: Accepted to MNRAS; data available at https://doi.org/10.5281/zenodo.3900980

In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion., Comment: 15 pages, 7 figures; accepted for publication in MNRAS

Recent observational studies of $\gamma$-ray emission from massive globular clusters (GCs) have revealed possible evidence of dark matter (DM) annihilation within GCs. It is, however, still controversial whether the emission comes from DM or from milli-second pulsars. We here present the new results of numerical simulations, which demonstrate that GCs with DM can originate from nucleated dwarfs orbiting the ancient MW. The simulated stripped nuclei (i.e., GCs) have the central DM densities ranging from 0.1 to several ${\rm M_\odot pc^{-3}}$, depending on the orbits and the masses of the host dwarf galaxies. However, GCs born outside the central regions of their hosts can have no/little DM after their hosts are destroyed and the GCs become the Galactic halo GCs. These results suggest that only GCs originating from stellar nuclei of dwarfs can possibly have DM. We further calculate the expected $\gamma$-ray emission from these simulated GCs and compare them to observations of $\omega$ Cen. Given the large range of DM densities in the simulated GCs, we suggest that the recent possible detection of DM annihilation from GCs should be more carefully interpreted., Comment: 5 pages, 4 figures, to be published in MNRAS

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

Binary millisecond pulsars (MSPs) provide several opportunities for research of fundamental physics. However, finding them can be challenging. Several subdwarf B (sdB) binary systems with possible neutron star companions have been identified, allowing us to perform a targeted search for MSPs within these systems. Six sdBs with companions in the neutron star mass range, as determined from their optical light curves, were observed with the Green Bank and Westerbork radio telescopes. The data were searched for periodic signals as well as single pulses. No radio pulsations from sdB systems were detected, down to an average sensitivity limit of 0.11 mJy. We did, however, discover a pulsar in the field of sdB HE0532-4503. Follow-up observations with the Giant Metrewave Radio Telescope showed that this pulsar, J0533-4524, is not spatially coincident with the sdB system. The pulsar has a relatively low magnetic field but still emits giant pulses. We place an upper limit of three to the number of radio pulsars in the six sdB systems. The non-detections may be explained by a combination of the MSP beaming fraction, luminosity, and a recycling fraction, 12 pages, 8 figures. Accepted for publication in MNRAS

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.

This paper presents the first detection of Shapiro delay from the binary millisecond pulsar PSR J1811–2405. We report a 11σ measurement of the orthometric amplitude, h3 = 6.8(6) × 10−7, and a 16σ measurement of the orthometric ratio, ς = 0.81(5). Given the relatively high orbital inclination, i = 79(2)°, of this binary system, we obtain constraints on the companion mass of $m_{\rm {c}}=0.31^{+0.08 }_{ -0.06}\, \mathrm{M}_{\odot }$. The pulsar mass is currently less well constrained, with a value of $2.0^{+0.8 }_{ -0.5}\, \mathrm{M}_{\odot }$. The companion mass and the orbital period are in agreement with the prediction made by previous numerical calculations of the evolution of compact binary systems. From a study of the polarization, we find that the orbital inclination angle is ∼100° and that PSR J1811–2405 is an orthogonal rotator. In addition, the μs-level timing precision together with its narrow profile makes PSR J1811–2405 a good candidate for inclusion in the pulsar timing arrays being used to detect nHz gravitational waves.

Young isolated neutron stars (INS) most commonly manifest themselves as rotationally powered pulsars (RPPs) which involve conventional radio pulsars as well as gamma-ray pulsars (GRPs) and rotating radio transients (RRATs). Some other young INS families manifest themselves as anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) which are commonly accepted as magnetars, i.e. magnetically powered neutron stars with decaying superstrong fields. Yet some other young INS are identified as central compact objects (CCOs) and X-ray dim isolated neutron stars (XDINSs) which are cooling objects powered by their thermal energy. Older pulsars, as a result of a previous long episode of accretion from a companion, manifest themselves as millisecond pulsars and more commonly appear in binary systems. We use Dirichlet process Gaussian mixture model (DPGMM), an unsupervised machine learning algorithm, for analyzing the distribution of these pulsar families in the parameter space of period and period derivative. We compare the average values of the characteristic age, magnetic dipole field strength, surface temperature and transverse velocity of all discovered clusters. We verify that DPGMM is robust and provides hints for inferring relations between different classes of pulsars. We discuss the implications of our findings for the magneto-thermal spin evolution models and fallback discs., Accepted by MNRAS (14 January 2020). 11 pages, 4 figures, 7 tables

Binaries harbouring millisecond pulsars enable a unique path to determine neutron star masses: radio pulsations reveal the motion of the neutron star, while that of the companion can be characterised through studies in the optical range. PSR J1012+5307 is a millisecond pulsar in a 14.5-h orbit with a helium-core white dwarf companion. In this work we present the analysis of an optical spectroscopic campaign, where the companion star absorption features reveal one of the lightest known white dwarfs. We determine a white dwarf radial velocity semi-amplitude of K_2 = 218.9 +- 2.2 km/s, which combined with that of the pulsar derived from the precise radio timing, yields a mass ratio of q=10.44+- 0.11. We also attempt to infer the white dwarf mass from observational constraints using new binary evolution models for extremely low-mass white dwarfs, but find that they cannot reproduce all observed parameters simultaneously. In particular, we cannot reconcile the radius predicted from binary evolution with the measurement from the photometric analysis (R_WD=0.047+-0.003 Rsun). Our limited understanding of extremely low-mass white dwarf evolution, which results from binary interaction, therefore comes as the main factor limiting the precision with which we can measure the mass of the white dwarf in this system. Our conservative white dwarf mass estimate of M_WD = 0.165 +- 0.015 Msun, along with the mass ratio enables us to infer a pulsar mass of M_NS = 1.72 +- 0.16 Msun. This value is clearly above the canonical 1.4 Msun, therefore adding PSR J1012+5307 to the growing list of massive millisecond pulsars., 13 pages, 7 figures. Accepted for publication in MNRAS

Past studies of the $\textit{Fermi}$ Galactic Center Excess (GCE) have found evidence for a ``high-energy tail'' in the GCE spectrum. Such a signature could be the result of inverse-Compton (IC) gamma rays produced by the injection of electrons/positrons into the interstellar medium by the putative population of millisecond pulsar (MSPs) responsible for the GCE. For this ICRC 2021 contribution, we present the results of an exhaustive study on simulated data, in which we analyze the detection potential of the forthcoming Cherenkov Telescope Array (CTA) to the high-energy tail of the $\textit{Fermi}$ GeV excess. In particular, we find that CTA will have sufficient sensitivity to detect this signal for physically reasonable electron/positron acceleration efficiencies, provided that the Galactic diffuse emission model (GDE) is well understood. Furthermore, we discuss the necessary conditions for a reliable CTA discovery in the case of a high degree of uncertainties in the GDE model. We also show that in the event that CTA observes an excess of diffuse gamma rays in the Galactic bulge, it will be able to discriminate between the dark matter and MSPs hypotheses, based on their distinct spatial morphologies.

Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes-curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars-has long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year Fermi-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$\sigma$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed in toto, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component-which we uncover at a significance of 8.2$\sigma$-is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this inverse Compton component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation.

The Fermi Large Area Telescope has observed a gamma-ray excess toward the center of the Galaxy at ~ GeV energies. The spectrum and intensity of the excess are consistent with the annihilation of dark matter with a mass of ~100 GeV and a velocity-averaged cross section of ~ 1e-26 cubic centimeter per second. In the meantime, a population of unresolved millisecond pulsars (MSPs) in the Galactic center remains a possible source of the excess. Furthermore, recent analyses have shown that the excess prefers the spatial morphology of the stellar bulge distribution in the Galactic center, supporting a MSP origin. The new discovery makes it imperative to further study the signals from MSPs. This dissertation studies the gamma-ray emission from Galactic millisecond pulsars to provide new insights into the origin of the Galactic center excess. Using the GALPROP code, we simulate the propagation of e± injected by the putative MSPs in the Galactic bulge and calculate the inverse Compton (IC) emission caused by the e± losing energy in the interstellar radiation field. We find recognizable features in the spatial maps of the IC. Above TeV energies, the IC morphology tends to follow the distribution of the injected e±. Then, we study the Cherenkov Telescope Array (CTA) sensitivity to the IC signal from MSPs. We find that the CTA has the potential to robustly discover the IC signature when the MSP e± injection efficiencies are in the range ≈ 2.9-74.1%. The CTA can also discriminate between an MSP and a dark matter origin for the radiating e± based on their different spatial maps. Next, we analyze the Fermi data from directions of Galactic globular clusters. The globular clusters are thought to be shining in gamma rays because of the MSP population they host. By analyzing their gamma-ray spectra, we reveal evidence for an IC component in the high-energy tail of Fermi data. Based on the IC component in the globular cluster spectra, the e± injection efficiency of millisecond pulsars is estimated to be slightly smaller than 10%. Finally, we study the spatial morphology of the 511 keV signal toward the Galactic center using data from INTEGRAL/SPI. We confirm that the 511 keV signal also traces the old stellar population in the Galactic bulge, which is similar to the Fermi GeV excess. Using a 3D smoothing kernel, we find that the signal is smeared out over a characteristic length scale of 150 ± 50 pc. We show that positron propagation prior to annihilation can explain the overall phenomenology of the 511 keV signal. Doctor of Philosophy Dark matter means matter that does not interact with light; therefore, they are invisible to traditional observations. We know that dark matter exists based on plenty of gravitational evidence: the motions of stars in galaxies, the large-scale structure of the Universe, the temperature fluctuations in the cosmic microwave background. However, we still know very little about the particle nature of dark matter. Detecting dark matter is one of the most extensive missions of modern physics. In indirect detection, the dark matter particles are expected to annihilate or decay in the cosmos, producing messenger particles that include gamma rays, cosmic rays, and neutrinos. Astronomical observations could detect those signals and confirm the nature of dark matter. However, understanding the astrophysical sources is essential for indirect detection of dark matter as they may emit similar signals. For a recent example, the Fermi Large Area Telescope launched by NASA is the most sensitive gamma-ray telescope in the energy range of ~ 100 MeV to ~ 100 GeV. It has detected an excess of gamma-ray signals toward the Galactic center consistent with what we expect from dark matter annihilation. However, millisecond pulsars, a type of fast rotating neutron stars, may also generate similar gamma-ray signals. Therefore, the origin of the signal remains unsettled. In this dissertation, we study different prospective of the gamma-ray emission from the millisecond pulsars in the Milky Way. We first study the inverse Compton signal from the millisecond pulsars in the Galactic bulge, caused by the relativistic e± injected by the millisecond pulsars. We find that the signal traces the original distribution of the e± above TeV energies. Next generation ground-based gamma-ray observatories like the Cherenkov Telescope Array (CTA) could be used to detect the signal. We study the CTA sensitivity to such an inverse Compton signal. We find that CTA can detect the inverse Compton signal from millisecond pulsars and discriminate it from a dark matter signal. We also study the gamma-ray emission from globular clusters in the Milky Way. They are dense collections of old stars orbiting our Galaxy, and they are known for hosting many millisecond pulsars. We reveal evidence for inverse Compton emission from the gamma-ray data of globular clusters. Our discovery helps us better understand the high-energy property of millisecond pulsars. Last, we study the morphology of the Galactic 511 keV signal caused by positron annihilation. Compact objects including millisecond pulsars are potential sources of the positrons. We find that the old stellar distribution with a smearing scale of ~ 150 pc best describes the 511 keV signal. Positron propagation from their sources prior to annihilation could explain the measured smearing scale.

The 559 Hz black-widow pulsar PSR J1555-2908, originally discovered in radio, is also a bright gamma-ray pulsar. Timing its pulsations using 12 yr of Fermi-LAT gamma-ray data reveals long-term variations in its spin frequency that are much larger than is observed from other millisecond pulsars. While this variability in the pulsar rotation rate could be intrinsic "timing noise", here we consider an alternative explanation: the variations arise from the presence of a very-low-mass third object in a wide multi-year orbit around the neutron star and its low-mass companion. With current data, this hierarchical-triple-system model describes the pulsar's rotation slightly more accurately than the best-fitting timing-noise model. Future observations will show if this alternative explanation is correct., 6 pages, 2 figures, submitted to ApJL

We report the discovery of two Black Widow millisecond pulsars in the globular cluster M28 with the MeerKAT telescope. PSR J1824$-$2452M (M28M) is a 4.78-ms pulsar in a $5.82\,$hour orbit and PSR J1824$-$2452N (M28N) is a 3.35-ms pulsar in a $4.76\,$hour orbit. Both pulsars have dispersion measures near $119.30\,$pc$\,$cm$^{-3}$ and have low mass companion stars ($\sim$$0.01-0.03\,$M$_\odot$), which do not cause strong radio eclipses or orbital variations. Including these systems, there are now five known black widow pulsars in M28. The pulsar searches were conducted as a part of an initial phase of MeerKAT's globular cluster census (within the TRAPUM Large Survey Project). These faint discoveries demonstrate the advantages of MeerKAT's survey sensitivity over previous searches and we expect to find additional pulsars in continued searches of this cluster., Comment: Accepted to ApJ on 12/28/2021 2/17/2022: Added Acknowledgements, Software, and Facilities sections. Fixed missing 'Total Proper Motion' units in table 2

We present a pulse timing analysis of NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658 during the outburst that started on 2022 August 19. Similar to previous outbursts, after decaying from a peak luminosity of $\simeq 1\times10^{36} \, \mathrm{erg \, s^{-1}}$ in about a week, the pulsar entered in a $\sim 1$ month-long reflaring stage. Comparison of the average pulsar spin frequency during the outburst with those previously measured confirmed the long-term spin derivative of $\dot{\nu}_{\textrm{SD}}=-(1.15\pm0.06)\times 10^{-15} \, \mathrm{Hz\,s^{-1}}$, compatible with the spin-down torque of a $\approx 10^{26} \, \mathrm{G \, cm^3}$ rotating magnetic dipole. For the first time in the last twenty years, the orbital phase evolution shows evidence for a decrease of the orbital period. The long-term behaviour of the orbit is dominated by a $\sim 11 \, \mathrm{s}$ modulation of the orbital phase epoch consistent with a $\sim 21 \, \mathrm{yr}$ period. We discuss the observed evolution in terms of a coupling between the orbit and variations in the mass quadrupole of the companion star., Comment: 11 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

The electric field accelerating electrons in the Timokhin-Arons polar-cap model as applied to millisecond pulsars is so high that the electron Lorentz factors are limited either by radiation reaction or by the Breit-Wheeler process. In the former case, it is possible to obtain an upper limit for curvature radiation momentum components perpendicular to the local magnetic field which is independent of the flux-line radius of curvature. The threshold value for single-photon conversion to pairs is high but is possibly reached in J0030+0451. However, owing to the high polar-cap temperature reported, direct pair production by the Breit-Wheeler process is probably more important. If the existence of coherent radio emission in millisecond pulsars is assumed to need a high-multiplicity pair plasma, then it follows that the primary electrons also produce gamma-rays in the Fermi-LAT energy band for which the magnetosphere is completely transparent. The absence of these, in phase with the radio emission, would be an immediate indication that ultra-high energy electrons and an active Timokhin-Arons polar cap are not present in J0030+0451., Comment: 4 pages; to be published in MNRAS Letters

4FGL J0935.5+0901, a $\gamma$-ray source recently identified as a candidate redback-type millisecond pulsar binary (MSP), shows an interesting feature of having double-peaked emission lines in its optical spectrum. The feature would further suggest the source as a transitional MSP system in the sub-luminous disk state. We have observed the source with \xmm\ and Five-hundred-meter Aperture Spherical radio Telescope (FAST) at X-ray and radio frequencies respectivelyfor further studies. From the X-ray observation, a bimodal count-rate distribution, which is a distinctive feature of the transitional MSP systems, is not detected, while the properties of X-ray variability and power-law spectrum are determined for the source. These results help establish the consistency of it being a redback in the radio pulsar state. However no radio pulsation signals are found in the FAST observation, resulting an upper limit on the flux density of $\sim 4\,\mu$Jy. Implications of these results are discussed., Comment: 10 pages, 5 figures

We present the results of simultaneous observations of the transitional millisecond pulsar (tMSP) candidate CXOU J110926.4-650224 with the XMM-Newton satellite and the MeerKAT telescope. The source was found at an average X-ray luminosity of $L_{\rm X}\simeq7\times10^{33}$ erg s$^{-1}$ over the 0.3-10 keV band (assuming a distance of 4 kpc) and displayed a peculiar variability pattern in the X-ray emission, switching between high, low and flaring modes on timescales of tens of seconds. A radio counterpart was detected at a significance of 7.9$\sigma$ with an average flux density of $\simeq$33$\mu$Jy at 1.28 GHz. It showed variability over the course of hours and emitted a $\simeq$10-min long flare just a few minutes after a brief sequence of multiple X-ray flares. No clear evidence for a significant correlated or anticorrelated variability pattern was found between the X-ray and radio emissions over timescales of tens of minutes and longer. CXOU J110926.4-650224 was undetected at higher radio frequencies in subsequent observations performed with the Australia Telescope Compact Array, when the source was still in the same X-ray sub-luminous state observed before, down to a flux density upper limit of 15$\mu$Jy at 7.25 GHz (at 3$\sigma$). We compare the radio emission properties of CXOU J110926.4-650224 with those observed in known and candidate tMSPs and discuss physical scenarios that may account for its persistent and flaring radio emissions., Comment: 10 pages, 4 figures. Accepted for publication on A&A

The origins of the diffuse flux of cosmogenic PeV neutrinos detected in the IceCube experiment during 2010 - 2017 remain unidentified. A population of extragalactic newly born fast spinning pulsars are investigated as possible candidates for generating the PeV energy scale neutrinos. A two-step mechanism of particle acceleration is adopted for transferring energy from star rotation to high energy electrons. Electrons might be boosted up to $\approx 0.01$ EeV energies and above, and produce ultra-high-energy (UHE) neutrinos and gamma rays when these electrons interact with low energy positrons and soft radiations in the acceleration zone. The theoretically derived extragalactic diffuse muon neutrino flux in the energy range [1-10] PeV is found consistent with the IceCube level if only a fraction ($\eta_{k} \approx 0.31\%$) of the total bolometric luminosity of pulsars are transferred to power the PeV neutrinos. Using the above value of the cosmic ray electron loading parameter $\eta_{k}$, the diffuse gamma ray flux from inverse Compton scattering off UHE electrons in the soft radiation field might be predicted from the model., Comment: 10 pages, 0 figures

We report the results of a deep search for faint Chandra X-ray sources, radio sources, and optical counterparts in the nearby, core-collapsed globular cluster, NGC 6752. We combined new and archival Chandra imaging to detect 51 X-ray sources (12 of which are new) within the 1.9 arcmin half-light radius. Three radio sources in deep ATCA 5 and 9 GHz radio images match with Chandra sources. We have searched for optical identifications for the expanded Chandra source list using deep Hubble Space Telescope photometry in B, R, H-alpha, UV, and U. Among the entire sample of 51 Chandra sources, we identify 18 cataclysmic variables (CVs), 9 chromospherically active binaries (ABs), 3 red giants (RGs), 3 galaxies (GLXs), and 6 active galactic nuclei (AGNs). Three of the sources are associated with millisecond pulsars (MSPs). As in our previous study of NGC 6752, we find that the brightest CVs appear to be more centrally concentrated than are the faintest CVs, although the effect is no longer statistically significant as a consequence of the inclusion in the faint group of two intermediate brightness CVs. This possible difference in the radial distributions of the bright and faint CV groups appears to indicate that mass segregation has separated them. We note that photometric incompleteness in the crowded central region of the cluster may also play a role. Both groups of CVs have an inferred mass above that of the main-sequence turnoff stars. We discuss the implications for the masses of the CV components., 35 pages, 30 figures

The Variables and Slow Transients Survey (VAST) on the Australian Square Kilometre Array Pathfinder (ASKAP) is designed to detect highly variable and transient radio sources on timescales from 5 s to$\sim\!5$yr. In this paper, we present the survey description, observation strategy and initial results from the VAST Phase I Pilot Survey. This pilot survey consists of$\sim\!162$h of observations conducted at a central frequency of 888 MHz between 2019 August and 2020 August, with a typical rms sensitivity of$0.24\ \mathrm{mJy\ beam}^{-1}$and angular resolution of$12-20$arcseconds. There are 113 fields, each of which was observed for 12 min integration time, with between 5 and 13 repeats, with cadences between 1 day and 8 months. The total area of the pilot survey footprint is 5 131 square degrees, covering six distinct regions of the sky. An initial search of two of these regions, totalling 1 646 square degrees, revealed 28 highly variable and/or transient sources. Seven of these are known pulsars, including the millisecond pulsar J2039–5617. Another seven are stars, four of which have no previously reported radio detection (SCR J0533–4257, LEHPM 2-783, UCAC3 89–412162 and 2MASS J22414436–6119311). Of the remaining 14 sources, two are active galactic nuclei, six are associated with galaxies and the other six have no multi-wavelength counterparts and are yet to be identified.

High sensitivity radio searches of unassociated $\gamma$-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of \textit{Fermi} Large Area Telescope (LAT) $\gamma$-ray sources. On Feb. 27$^{th}$, 2018 we discovered an isolated millisecond pulsar (MSP), PSR J0318+0253, coincident with the unassociated $\gamma$-ray source 3FGL J0318.1+0252. PSR J0318+0253 has a spin period of $5.19$ milliseconds, a dispersion measure (DM) of $26$ pc cm$^{-3}$ corresponding to a DM distance of about $1.3$ kpc, and a period-averaged flux density of $\sim$11 $\pm$ 2 $\mu$Jy at L-band (1.05-1.45 GHz). Among all high energy MSPs, PSR J0318+0253 is the faintest ever detected in radio bands, by a factor of at least $\sim$4 in terms of L-band fluxes. With the aid of the radio ephemeris, an analysis of 9.6 years of \textit{Fermi}-LAT data revealed that PSR J0318+0253 also displays strong $\gamma$-ray pulsations. Follow-up observations carried out by both Arecibo and FAST suggest a likely spectral turn-over around 350 MHz. This is the first result from the collaboration between FAST and the \textit{Fermi}-LAT teams as well as the first confirmed new MSP discovery by FAST, raising hopes for the detection of many more MSPs. Such discoveries will make a significant contribution to our understanding of the neutron star zoo while potentially contributing to the future detection of gravitational waves, via pulsar timing array (PTA) experiments., Comment: 7 pages, 2 figures, 1 table

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semi-analytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling shown the capability of coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including $\mathcal{O}(10^2)$ pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also high beneficial to include DNest in current toolboxes of PTA analysis., 20 pages, 10 figures, submitted to ApJ

X-ray pulsar-based navigation (XNAV) is a promising autonomous navigation method, and the pulse phase is the basic measurement of XNAV. However, the current methods for estimating the pulse phase for orbiting spacecraft have a high computational cost. This paper proposes a stellar angle measurement-aided pulse phase estimation method for high Earth orbit (HEO) spacecraft, with the aim of reducing the computational cost of pulse phase estimation in XNAV. In this pulse phase estimation method, the effect caused by the orbital motion of the spacecraft is roughly removed by stellar angle measurement. Furthermore, a deeply integrated navigation method using the X-ray pulsar and the stellar angle is proposed. The performances of the stellar angle measurement-aided pulse phase estimation method and the integrated navigation method were verified by simulation. The simulation results show that the proposed pulse phase estimation method can handle the signals of millisecond pulsars and achieve pulse phase estimation with lower computational cost than the current methods. In addition, for HEO spacecraft, the position error of the proposed integrated navigation method is lower than that of the stellar angle navigation method.

We report observed and derived timing parameters for three millisecond pulsars (MSPs) from observations collected with the Parkes 64-m telescope, Murriyang. The pulsars were found during re-processing of archival survey data by Mickaliger et al. One of the new pulsars (PSR J1546-5925) has a spin period $P=7.8$ ms and is isolated. The other two (PSR J0921-5202 with $P=9.7$ ms and PSR J1146-6610 with $P=3.7$ ms) are in binary systems around low-mass ($>0.2 M_{\odot}$) companions. Their respective orbital periods are $38$.2 d and $62.8$ d. While PSR J0921-5202 has a low orbital eccentricity $e=1.3 \times 10^{-5}$, in keeping with many other Galactic MSPs, PSR J1146-6610 has a significantly larger eccentricity, $e = 7.4 \times 10^{-3}$. This makes it a likely member of a group of eccentric MSP-He white dwarf binary systems in the Galactic disk whose formation is poorly understood. Two of the pulsars are co-located with previously unidentified point sources discovered with the Fermi satellite's Large Area Telescope, but no $\gamma$-ray pulsations have been detected, likely due to their low spin-down powers. We also show that, particularly in terms of orbital diversity, the current sample of MSPs is far from complete and is subject to a number of selection biases., Comment: 8 pages, 2 figures, accepted for publication in MNRAS

Pulsar winds interacting with sources of external pressure are well-established as efficient and prolific TeV accelerators in our Galaxy. Yet, enabled by observations from Fermi-LAT, a growing class of non-accreting pulsars in binaries has emerged and these are likely to become apparent as TeV emitters in the CTA era. This class consists of the black widows and redbacks, binaries in which a millisecond pulsar interacts with its low-mass companion. In such systems, an intrabinary shock can form as a site of particle acceleration and associated nonthermal emission. We motivate why these sources are particularly interesting for understanding pulsar winds. We also describe our new multizone code which models the X-ray and gamma-ray synchrotron and inverse Compton spectral components for select spider binaries, including diffusion, convection, and radiative energy losses in an axially symmetric, steady-state approach. This new multizone code simultaneously yields energy-dependent light curves and orbital-phase-resolved spectra. It also better constrains the multiplicity of electron/positron pairs that have been accelerated up to TeV energies and are necessary to power orbitally-modulated synchrotron emission components between the X-rays and MeV/GeV bands potentially observed in some systems. This affords a more robust prediction of the expected high-energy and VHE gamma-ray flux. Nearby MSPs with hot or flaring companions may be promising targets for CTA, and it is possible that spider binaries could contribute to the observed AMS-02 energetic positron excess., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC 2021)

If the mysterious Fermi-LAT GeV gamma-ray excess is due to an unresolved population of millisecond pulsars (MSP) in the Galactic bulge, one expects this very same population to shine in X rays. For the first time, we address the question of what is the sensitivity of current X-ray telescopes to an MSP population in the Galactic bulge. To this end, we create a synthetic population of Galactic MSPs, building on an empirical connection between gamma- and X-ray MSP emission based on observed source properties. We compare our model with compact sources in the latest Chandra source catalog, applying selections based on spectral observables and optical astrometry with Gaia. We find a significant number of Chandra sources in the region of interest to be consistent with being bulge MSPs that are as yet unidentified. This motivates dedicated multi-wavelength searches for bulge MSPs: Some promising directions are briefly discussed., 14 pages, 7 figures. Matches version published in PRD

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background., 10 pages, 5 figures

Ground-based Air-Cherenkov telescopes have detected pulsations at energies above 50 GeV from a growing number of {\sl Fermi} pulsars. These include the Crab, Vela, PSR~B1706-44 and Geminga pulsar, with the first two having pulsed detections above 1 TeV. In some cases, there appears to be very-high-energy (VHE) emission that is an extension of the {\sl Fermi} spectra to higher energies, while in other cases, additional higher-energy spectral components that require a separate emission mechanism may be present. We will present results of broadband spectral modeling using global magnetosphere fields and multiple emission mechanisms that include synchro-curvature (SC) and inverse Compton scattered (ICS) radiation from accelerated particles (primaries) and synchrotron self-Compton (SSC) emission from lower-energy pairs. Our models predict three distinct VHE components: SC from primaries whose high-energy tail can extend to 100 GeV, SSC from pairs that can extend to several TeV and ICS from primary particles accelerated in the current sheet, scattering pair synchrotron radiation, that appears beyond 10 TeV. Model spectra show a wide range of VHE flux, with detectable pair SSC and primary ICS components expected for Crab-like pulsars and some millisecond pulsars but only a primary ICS component for Vela. Our models suggest that H.E.S.S.-II and MAGIC have detected the high-energy tail of the primary SC component that produces the {\sl Fermi} spectrum in Vela, Geminga and PSR B1706$-$44. We argue that the ICS component peaking above 10 TeV from Vela has been seen by H.E.S.S. Detection of this emission component from the Crab and other pulsars is possible with HAWC and CTA, and directly measures the maximum particle energy in pulsars.

The repeating fast radio burst (FRB) localized to a globular cluster in M81 challenges our understanding of FRB models. In this Letter, we explore dynamical formation scenarios for objects in old globular clusters that may plausibly power FRBs. Using N-body simulations, we demonstrate that young neutron stars may form in globular clusters at a rate of up to $\sim50\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ through a combination of binary white dwarf mergers, white dwarf--neutron star mergers, binary neutron star mergers, and accretion induced collapse of massive white dwarfs in binary systems. We consider two FRB emission mechanisms: First, we show that a magnetically-powered source (e.g., a magnetar with field strength $\gtrsim10^{14}\,$G) is viable for radio emission efficiencies $\gtrsim10^{-4}$. This would require magnetic activity lifetimes longer than the associated spin-down timescales and longer than empirically-constrained lifetimes of Galactic magnetars. Alternatively, if these dynamical formation channels produce young rotation-powered neutron stars with spin periods of $\sim10\,$ms and magnetic fields of $\sim10^{11}\,$G (corresponding to spin-down lifetimes of $\gtrsim10^5\,$yr), the inferred event rate and energetics can be reasonably reproduced for order unity duty cycles. Additionally, we show that recycled millisecond pulsars or low-mass X-ray binaries similar to those well-observed in Galactic globular clusters may also be plausible channels, but only if their duty cycle for producing bursts similar to the M81 FRB is small., 11 pages, 2 tables, 1 figure. Accepted for publication in ApJ Letters. Comments welcome

A numerical rotating neutron star solver is used to study the temporal evolution of accreting neutron stars using a multi-polytrope model for the nuclear equation of state named ACB5. The solver is based on a quadrupole expansion of the metric, but confirms the results of previous works, revealing the possibility of an abrupt transition of a neutron star from a purely hadronic branch to a third-family branch of stable hybrid stars, passing through an unstable intermediate branch. The accretion is described through a sequence of stationary rotating {stellar} configurations which lose angular momentum through magnetic dipole emission while, at the same time, gaining angular momentum through mass accretion. The model has several free parameters which are inferred from observations. The mass accretion scenario is studied in dependence on the effectiveness of angular momentum transfer which determines at which spin frequency the neutron star will become unstable against gravitational collapse to the corresponding hybrid star on the stable third-family branch. It is conceivable that the neutrino burst which accompanies the deconfinement transition may trigger a pulsar kick which results in the eccentric orbit. A consequence of the present model is the prediction of a correlation between the spin frequency of the millisecond pulsar in the eccentric orbit and its mass at birth., Comment: 7 pages, 4 figures, 1 table, accepted for publication in Astronomische Nachrichten

Pulsed γ-ray properties of 22 millisecond pulsars (MSPs) with double peaks are studied in the revised versions of the outer gap (OG) and the stripped wind (SW) models. The major differences between these two models are magnetic field structures, γ-ray production sites, and radiation mechanisms. In the models, γ-ray light curves of these MSPs are calculated through a Markov chain Monte Carlo (MCMC) method to pick best-fitting model parameters. Our results indicate that (1) both models can reproduce observed double-peak structures of the MSPs well, for most MSPs, a relatively large magnetic inclination angle (α ∼ 50°–90°) and small viewing angle (ζ ∼ 30°–90°) are obtained in the OG model, but a relatively small α ∈ (20°, 60°) and large ζ ∈ (70°, 150°) in the SW model; (2) phase-averaged spectra, cut-off energy, and γ-ray luminosity calculated in both models are consistent with observed those by Fermi-Large Area Telescope (LAT). Therefore, it may be concluded that the OG and the SW models for describing high-energy emissions of MSPs cannot be ruled out at present, and further studies are required.

We investigate the spin dynamics of a millisecond pulsar (MSP) in a tightly bounded orbit around a massive black hole. These binaries are progenitors of the extreme-mass-ratio-inspirals (EMRIs) and intermediate-mass-ratio-inspirals (IMRIs) gravitational wave events. The Mathisson-Papapetrou-Dixon (MPD) formulation is used to determine the orbital motion and spin modulation and evolution. We show that the MSP will not be confined in a planar Keplerian orbit and its spin will exhibit precession and nutation induced by spin-orbit coupling and spin-curvature interaction. These spin and orbital behaviours will manifest observationally in the temporal variations in the MSP's pulsed emission and, with certain geometries, in the self-occultation of the pulsar's emitting poles. Radio pulsar timing observations will be able to detect such signatures. These extreme-mass-ratio binaries (EMRBs) and intermediate-mass-ratio binaries (IMRBs) are also strong gravitational wave sources. Combining radio pulsar timing and gravitational wave observations will allow us to determine the dynamics of these systems in high precision and hence the subtle behaviours of spinning masses in strong gravity.

We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distribution with power-slope $-1.75$, fully consistent with the prediction of Bahcall-Wolf pioneering work. Using the vast amount of data available, we infer the rate for tidal disruption events, being $4 \times 10^{-6}$ per yr, and estimate the number of objects that emit gravitational waves during the phases preceding the accretion onto the super-massive black hole, $\sim 270$ per Gyr. We show that some of these sources could form extreme mass-ratio inspirals. We follow the evolution of binary stars population, showing that the initial binary fraction of $5\%$ drops down to $2.5\%$ inside the inner parsec. Also, we explored the possible formation of binary systems containing a compact object, discussing the implications for millisecond pulsars formation and the development of Ia Supernovae., Comment: Accepted by MNRAS after a major revision

We report the discovery of a new low-mass X-ray binary near the center of the unassociated Fermi GeV gamma-ray source 4FGL J0540.0-7552. The source shows the persistent presence of an optical accretion disk and exhibits extreme X-ray and optical variability. It also has an X-ray spectrum well-fit by a hard power law with a Gamma = 1.8 and a high ratio of X-ray to gamma-ray flux. Together, these properties are consistent with the classification of the binary as a transitional millisecond pulsar (tMSP) in the sub-luminous disk state. Uniquely among the candidate tMSPs, 4FGL J0540.0-7552 shows consistent optical, X-ray, and gamma-ray evidence for having undergone a state change, becoming substantially brighter in the optical and X-rays and fainter in GeV gamma-rays sometime in mid-2013. In its current sub-luminous disk state, and like one other candidate tMSP in the Galactic field, 4FGL J0540.0-7552 appears to always be in an X-ray "flare mode", indicating that this could be common phenomenology for tMSPs., ApJ in press