Globular Cluster Terzan Research Articles

Simulation-based Inference of Radio Millisecond Pulsars in Globular Clusters

Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP population in GCs through a careful study of their luminosity function. Parameters of interest are the number of MSPs hosted by a GC, as well as the mean and the width of their luminosity function, which are typically affected by large uncertainties. While, as we show, likelihood-based studies can lead to ill-behaved posteriors on the size of the MSP population, we introduce a novel, likelihood-free analysis, based on marginal neural ratio estimation, which consistently produces well-behaved posteriors. We focus on the GC Terzan 5 (or Ter 5), which currently counts 48 detected MSPs. We find that 158 −104+294 MSPs should be hosted in this GC, but the uncertainty on this number remains large. We explore the performance of our new method on simulated Terzan 5-like data sets mimicking possible future observational outcomes. We find that significant improvement on the posteriors can be obtained by adding a reliable measurement of the diffuse radio emission of the GC to the analysis or by improving the detection threshold of current radio pulsation surveys by at least a factor of 2.

Discovery of a Millisecond Pulsar Associated with Terzan 6

Observations show that globular clusters (GCs) might be among the best places to find millisecond pulsars (MSPs). However, the GC Terzan 6 seems to be an exception without any pulsar discovered, although its high stellar encounter rate suggests that it harbors dozens of them. We report the discovery of the first radio pulsar, PSR J1750–3116A, likely associated with Terzan 6 in a search of C-band (4–8 GHz) data from the Robert C. Byrd Green Bank Telescope with a spin period of 5.33 ms and dispersion measure (DM) ≃ 383 pc cm−3. The mean flux density of this pulsar is approximately 3 μJy. The DM agrees well with predictions from the Galactic free electron density model, assuming a distance of 6.7 kpc for Terzan 6. PSR J1750–3116A is likely an isolated MSP, potentially formed through dynamical interactions, considering the core-collapsed classification and the exceptionally high stellar encounter rate of Terzan 6. This is the highest radio frequency observation that has led to the discovery of a pulsar in a GC to date. While L-band (1–2 GHz) observations of this cluster are unlikely to yield significant returns due to propagation effects, we predict that further pulsar discoveries in Terzan 6 will be made by existing radio telescopes at higher frequencies.

Digging deeper into the dense Galactic globular cluster Terzan 5 with electron-multiplying CCDs

Digging deeper into the dense Galactic globular cluster Terzan 5 with electron-multiplying CCDs

Teraelectronvolt gamma-ray emission near globular cluster Terzan 5 as a probe of cosmic ray transport

Teraelectronvolt gamma-ray emission near globular cluster Terzan 5 as a probe of cosmic ray transport

Detailed chemical abundances of the globular cluster Terzan 6 in the inner bulge

We used near-infrared spectroscopy at medium-high resolution (R = 8000 − 25 000) to perform the first comprehensive chemical study of the intermediate luminosity bulge globular cluster Terzan 6. We derived detailed abundances and abundance patterns of 27 giant stars, likely members of Terzan 6, based on their accurate Hubble Space Telescope proper motions and line-of-sight radial velocities. From the spectral analysis of these stars, we determined an average heliocentric radial velocity of 143.3 ± 1.0 km s−1 with a velocity dispersion of 5.1 ± 0.7 km s−1 and an average [Fe/H] = −0.65 ± 0.01 and a low 1σ dispersion of 0.03 dex. We also measured some depletion of [Mn/Fe] with respect to the solar-scaled values and enhancement of for [Ca/Fe], [Si/Fe], [Mg/Fe], [Ti/Fe], [O/Fe], [Al/Fe], [Na/Fe], and, to a lower extent, for [K/Fe], consistent with previous measurements of other bulge globular clusters and favoring the scenario of a rapid bulge formation and chemical enrichment. Some spread in the light element abundances suggest the presence of first- and second-generation stars, typical of genuine globulars. Finally, we measured some depletion of carbon and low 12C/13C isotopic ratios, as in previous studies of field and cluster bulge giants, indicating that extra-mixing mechanisms should be at work during the post main sequence evolution in the high metallicity regime as well.

Discovery and timing of ten new millisecond pulsars in the globular cluster Terzan 5

We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856–1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for all ten discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748−2446ao which is an eccentric (e = 0.32) wide-orbit (orbital period Pb = 57.55 d) system. We were able to measure the rate of advance of periastron (ω̇) for this system allowing us to determine a total mass of 3.17 ± 0.02 M⊙. With a minimum companion mass (Mc) of ∼0.8 M⊙, PSR J1748−2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar (P = 2.27 ms) and the longest orbital period measured for any known DNS system. PSR J1748−2446ap has the second highest eccentricity for any recycled pulsar (e ∼ 0.905) and for this system we can measure the total mass (1.997 ± 0.006 M⊙) and estimate the pulsar and companion masses, (1.700−0.045+0.015 M⊙ and 0.294−0.014+0.046 M⊙, respectively). PSR J1748−2446ar is an eclipsing redback (minimum Mc ∼ 0.34 M⊙) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. We were also able to detect ω̇ for PSR J1748−2446au leading to a total mass estimate of 1.82 ± 0.07 M⊙ and indicating that the system is likely the result of Case A Roche lobe overflow. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies.

Discovery of a Second Eclipsing, ​​​​​Bursting Neutron Star Low-mass X-Ray Binary in the Globular Cluster Terzan 6

We have analyzed Chandra and Suzaku observations of the globular cluster Terzan 6, made when the recurrent transient GRS 1747–312 was in quiescence. Our analysis reveals the presence of a second eclipsing, bursting neutron star low-mass X-ray binary in the central regions of the cluster, in addition to GRS 1747–312. The new source, which we name Terzan 6 X2, is located only ∼0.″7 away from GRS 1747–312 in the 2021 Chandra images. The detection of a 5.14 ks long eclipse in the light curve of X2 at a time not predicted by the ephemeris of GRS 1747–312 confirms that it is an unrelated source. Using the Suzaku light curve from 2009, which, in addition to a type I X-ray burst, also showed an eclipse-like feature, we constrain the orbital period to be longer than 16.27 hr. The 0.5–10 keV luminosities of X2 vary in the range of ∼0.24–5.9 × 1034 erg s−1 on timescales of months to years. We have identified a plausible optical counterpart of X2 in Hubble Space Telescope F606W and F814W images. This star varied by 2.7 mag in V 606 between epochs separated by years. In the cluster color–magnitude diagram, the variable counterpart lies in the blue-straggler region when it was optically bright, about 1.1–1.7 mag above the main-sequence turn-off. From the orbital period–density relation of Roche-lobe filling stars, we find the mass-donor radius to be ≳0.8 R ⊙.

Algorithmic Pulsar Timer for Binaries

Pulsar timing is a powerful tool that, by accounting for every rotation of a pulsar, precisely measures the spin frequency, spin frequency derivatives, astrometric position, binary parameters when applicable, properties of the interstellar medium, and potentially general relativistic effects. Typically, this process demands fairly stringent scheduling requirements for monitoring observations as well as deep domain knowledge to “phase connect” the timing data. We present an algorithm that automates the pulsar-timing process for binary pulsars, whose timing solutions have an additional level of complexity, although the algorithm works for isolated pulsars as well. Using the statistical F-test and the quadratic dependence of the reduced χ 2 near a minimum, the global rotation count of a pulsar can be determined efficiently and systematically. We have used our algorithm to establish timing solutions for two newly discovered binary pulsars, PSRs J1748−2446aq and J1748−2446at, in the globular cluster Terzan 5, using ∼70 Green Bank Telescope observations from the last 13 yr.

A possible third body in the X-ray system GRS 1747−312 and models with higher order multiplicity

ABSTRACT GRS 1747−312 is a bright Low-Mass X-ray Binary in the globular cluster Terzan 6, located at a distance of 9.5 kpc from the Earth. It exhibits regular outbursts approximately every 4.5 months, during which periodic eclipses are known to occur. These eclipses have only been observed in the outburst phase, and are not clearly seen when the source is quiescent. Recent Chandra observations of the source were performed in 2019 June and 2021 April, June, and August. Two of these observations captured the source during its outburst, and showed clear flux decreases at the expected time of eclipse. The other two observations occurred when the source was quiescent. We present the discovery of a dip that occurred during the quiescent state. The dip is of longer duration and its time of occurrence does not fit the ephemeris of the shorter eclipses. We study the physical characteristics of the dip and determine that it has all the properties of an eclipse by an object with a well defined surface. We find that there are several possibilities for the nature of the object causing the 5.3 ks eclipse. First, GRS 1747−312 may be an X-ray triple, with an LMXB orbited by an outer third object, which could be an M-dwarf, brown dwarf, or planet. Secondly, there could be two LMXBs in close proximity to each other, potentially bound together. Whatever the true nature of the eclipser, its presence suggests that the GRS 1747−312 system is exotic.

A highly magnetized environment in a pulsar binary system.

Spider pulsars are millisecond pulsars in short-period (≲12-h) orbits with low-mass (~0.01-0.4 M⊙) companion stars. The pulsars ablate plasma from the companion star, causing time delays and eclipses of the radio emission from the pulsar. The magnetic field of the companion has been proposed to strongly influence both the evolution of the binary system1 and the eclipse properties of the pulsar emission2. Changes in the rotation measure (RM) have been seen in a spider system, implying that there is an increase in themagnetic field near the eclipse3. Here we report a diverse range of evidence for a highly magnetized environment in the spider system PSR B1744 - 24A4, located in the globular cluster Terzan 5. We observe semi-regular profile changes to the circular polarization, V, when the pulsar emission passes close to the companion. This suggests that there is Faraday conversion where the radio wave tracks a reversal in the parallel magnetic field and constrains the companion magnetic field, B (> 10 G). We also see irregular, fast changes in the RM at random orbital phases, implying that the magnetic strength of the stellar wind, B, is greater than 10 mG. There are similarities between the unusual polarization behaviour of PSR B1744 - 24A and some repeating fast radio bursts (FRBs)5-7. Together with the possible binary-produced long-term periodicity of two active repeating FRBs8,9, and the discovery of a nearby FRB in a globular cluster10, where pulsar binaries are common, these similarities suggest that a proportion of FRBs have binary companions.

A Deep Chandra X-Ray Observatory Study of the Millisecond Pulsar Population in the Globular Cluster Terzan 5

Abstract We present an analysis of 745.8 ks of archival Chandra X-Ray Observatory Advanced CCD Imaging Spectrometer data accumulated between 2000 and 2016 of the millisecond pulsar (MSP) population in the rich Galactic globular cluster Terzan 5. Eight of the 38 MSPs with precise positions are found to have plausible X-ray source matches. Despite the deep exposure, the remaining MSPs are either marginally detected or have no obvious X-ray counterparts, which can be attributed to the typically soft thermal spectra of rotation-powered MSPs, which are strongly attenuated by the high intervening absorbing column (∼1022 cm−2) toward the cluster, and in some instances to severe source crowding/blending. For the “redback” MSP binaries PSR J1748−2446P and PSR J1748−2446ad and the “black widow” binary system PSR J1748−2446O, we find clear evidence for large-amplitude X-ray variability at the orbital period consistent with an intrabinary shock origin. The third redback MSP in the cluster, PSR J1748−2446A, shows order-of-magnitude variations in flux on timescales of years, possibly due to state transitions or intense flaring episodes from a magnetically active secondary star.

Difference in Chemical Composition between the Bright and Faint Red Clump Stars in the Milky Way Bulge

Abstract The double red clump (RC) observed in color–magnitude diagrams of the Milky Way bulge is at the heart of the current debate on the structure and formation origin of the bulge. This feature can be explained by the difference between the two RCs either in distance (“X-shaped scenario”) or in chemical composition (“multiple-population scenario”). Here we report our high-resolution spectroscopy for the RC and red giant branch stars in a high-latitude field (b ∼ −8.5°) of the bulge. We find a difference in [Fe/H] between the stars in the bright and faint RC regimes, in the sense that the bright stars are enhanced in [Fe/H] with respect to the faint stars by 0.149 ± 0.036 dex. The stars on the bright RC are also enhanced in [Na/Fe] but appear to be depleted in [Al/Fe] and [O/Fe], although more observations are required to confirm the significance of these differences. Interestingly, these chemical patterns are similar to those observed among multiple stellar populations in the metal-rich bulge globular cluster Terzan 5. In addition, we find a number of Na-rich stars, which would corroborate the presence of multiple populations in the bulge. Our results support an origin of the double RC from dissolved globular clusters that harbor multiple stellar populations. Thus, our study suggests that a substantial fraction of the outer bulge stars would have originated from the assembly of such stellar systems in the early phase of the Milky Way formation.

The MAVERIC Survey: New Compact Binaries Revealed by Deep Radio Continuum Observations of the Galactic Globular Cluster Terzan 5

Abstract Owing to its massive, dense core, Terzan 5 has the richest population of millisecond pulsars (MSPs) known among Galactic globular clusters. Here we report new deep 2–8 GHz radio continuum observations of Terzan 5 obtained with the Karl G. Jansky Very Large Array. We have identified a total of 24 sources within the cluster half-light radius, including 17 within the core radius. Nineteen are associated with previously studied MSPs and X-ray binaries. Three of the new radio sources have steep radio spectra and are located within the cluster core, as expected for MSPs. These three sources have hard X-ray photon indices (Γ = 1.3–1.5) and highly variable X-ray emission, suggesting they are binary MSPs belonging to the spider class. For the most X-ray luminous of these sources, the redback spider classification is confirmed by its X-ray light curve, which shows an orbital period of 12.32 hr and double-peaked structure around X-ray maximum. The likely discovery of bright binary MSPs in a well-studied cluster like Terzan 5 highlights how deep radio continuum imaging can complement pulsar search and timing observations in finding probable eclipsing systems. The other new radio source in the core has a flat radio spectrum and is X-ray faint ( erg s−1) with a photon index Γ = 2.1 ± 0.5, consistent with the properties expected for a quiescent stellar-mass black hole X-ray binary.

The MeerKAT telescope as a pulsar facility: System verification and early science results from MeerTime

AbstractWe describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (${\sim}2.8\,\mbox{K Jy}^{-1}$) low-system temperature (${\sim}18\,\mbox{K at }20\,\mbox{cm}$) radio array that currently operates at 580–1 670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar$\mbox{J}0737{-}3039\mbox{A}$, pulse profiles from 34 millisecond pulsars (MSPs) from a single 2.5-h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR$\mbox{J}0540{-}6919$, and nulling identified in the slow pulsar PSR J0633–2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright MSPs confirm that MeerKAT delivers exceptional timing. PSR$\mbox{J}2241{-}5236$exhibits a jitter limit of$<4\,\mbox{ns h}^{-1}$whilst timing of PSR$\mbox{J}1909{-}3744$over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1 000 pulsars per day and the future deployment of S-band (1 750–3 500 MHz) receivers will further enhance its capabilities.

A MUSE study of the inner bulge globular cluster Terzan 9: a fossil record in the Galaxy

Context. Moderately metal-poor inner bulge globular clusters are relics of a generation of long-lived stars that formed in the early Galaxy. Terzan 9, projected at 4°.12 from the Galactic center, is among the most central globular clusters in the Milky Way, showing an orbit which remains confined to the inner 1 kpc. Aims. Our aim is the derivation of the cluster’s metallicity, together with an accurate measurement of the mean radial velocity. In the literature, metallicities in the range between −2.0 < [Fe/H] < −1.0 have been estimated for Terzan 9 based on color-magnitude diagrams and CaII triplet (CaT) lines. Methods. Given its compactness, Terzan 9 was observed using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope. The extraction of spectra from several hundreds of individual stars allowed us to derive their radial velocities, metallicities, and [Mg/Fe]. The spectra obtained with MUSE were analysed through full spectrum fitting using the ETOILE code. Results. We obtained a mean metallicity of [Fe/H] ≈ −1.10 ±0.15, a heliocentric radial velocity of vhr = 58.1 ± 1.1 km s−1, and a magnesium-to-iron [Mg/Fe] = 0.27 ± 0.03. The metallicity-derived character of Terzan 9 sets it among the family of the moderately metal-poor Blue Horizontal Branch clusters HP 1, NGC 6558, and NGC 6522.

Exploring the Mass Segregation Effect of X-Ray Sources in Globular Clusters. II. The Case of Terzan 5

Abstract Using archival Chandra observations with a total effective exposure of 734 ks, we derive an updated catalog of point sources in the massive globular cluster (GC) Terzan 5. Our catalog covers an area of 58.1 arcmin2 (R ≤ 4.′3) with 489 X-ray sources, and more than 75% of these sources are first detected in this cluster. We find significant dips in the radial distribution profiles of X-ray sources in Terzan 5, with the projected distance and width of the distribution dips for bright (L X ≳ 9.5 × 1030 erg s−1) X-ray sources larger than those of the faint (L X ≲ 9.5 × 1030 erg s−1) sources. By fitting the radial distribution of the X-ray sources with a “generalized King model,” we estimated an average mass of 1.48 ± 0.11 and 1.27 ± 0.13 M ⊙ for the bright and faint X-ray sources, respectively. These results are in agreement with that observed in 47 Tuc, which may suggest a universal mass segregation effect for X-ray sources in GCs. Compared with 47 Tuc, we show that the two-body relaxation timescale of Terzan 5 is much smaller, but its dynamical age is significantly younger than 47 Tuc. These features suggest that the evolution of Terzan 5 is not purely driven by two-body relaxation, and the tidal stripping effect also plays an important role in accelerating the dynamical evolution of this cluster.

Continued cooling of the accretion-heated neutron star crust in the X-ray transient IGR J17480–2446 located in the globular cluster Terzan 5

We present a new Chandra observation (performed in July 2016) of the neutron star X-ray transient IGR J17480-2446, located in the globular cluster Terzan 5. We study the continued cooling of the neutron star crust in this system that was heated during the 2010 outburst of the source. This new observation was performed two years after the last observation of IGR J17480-2446, hence, significantly extending the cooling baseline. We reanalysed all available Chandra observations of the source (but excluding observations during which one of the known transients in Terzan 5 was in outburst) and fitted the obtained cooling curve with our cooling code NSCool, which allows for much improved modelling than what was previously performed for the source. The data and our fit models indicate that the crust was still cooling ~5.5 years after the outburst ended. The neutron star crust has likely not reached crust-core thermal equilibrium yet, and further cooling is predicted (which can be confirmed with additional Chandra observations in >5 years). Intriguingly, we find indications that the thermal conductivity might be relatively low in part of the crust compared to what has been inferred for other crust-cooling sources and tentatively suggest that this layer might be located around the neutron drip. The reason for this difference is unclear, but might be related to the fact that IGR J17480-2446 harbours a relatively slowly rotating neutron star (with a spin of 11 Hz) that has a relatively strong inferred surface magnetic field ($10^{9-10}$ Gauss) compared to what is known or typically assumed for other cooling sources.

Loosely coherent search in LIGO O1 data for continuous gravitational waves from Terzan 5 and the Galactic Center

We report results of a search for continuous gravitational waves from a region covering the globular cluster Terzan 5 and the galactic center. Continuous gravitational waves are expected from fast-spinning, slightly non-axisymmetric isolated neutron stars as well as more exotic objects. The regions that we target are believed to be unusually abundant in neutron stars. We use a new loosely coherent search method that allows to reach unprecedented levels of sensitivity for this type of search. The search covers the frequency band 475-1500 Hz and frequency time derivatives in the range of [-3e-8, +1e-9] Hz/s, which is a parameter range not explored before with the depth reached by this search. As to be expected with only a few months of data from the same observing run, it is very difficult to make a confident detection of a continuous signal over such a large parameter space. A list of parameter space points that passed all the thresholds of this search is provided. We follow-up the most significant outlier on the newly released O2 data and cannot confirm it. We provide upper limits on the gravitational wave strength of signals as a function of signal frequency.

Halo intruders in the Galactic bulge revealed by HST and Gaia: the globular clusters Terzan 10 and Djorgovski 1

Context. The low-latitude globular clusters Terzan 10 and Djorgovski 1 are projected in the Galactic bulge, in a Galactic region highly affected by extinction. A discrepancy of a factor of ∼2 exists in the literature in regards to the distance determination of these clusters. Aims. We revisit the colour-magnitude diagrams (CMDs) of these two globular clusters with the purpose of disentangling their distance determination ambiguity and, for the first time, of determining their orbits to identify whether or not they are part of the bulge/bar region. Methods. We use Hubble Space Telescope CMDs, with the filters F606W from ACS and F160W from WFC3 for Terzan 10, and F606W and F814W from ACS for Djorgosvski 1, and combine them with the proper motions from Gaia Data Release 2. For the orbit integrations, we employed a steady Galactic model with bar. Results. For the first time the blue horizontal branch of these clusters is clearly resolved. We obtain reliable distances of d⊙ = 10.3 ± 1.0 kpc and 9.3 ± 0.5 kpc for Terzan 10, and Djorgovski 1 respectively, indicating that they are both currently located in the bulge volume. From Gaia DR2 proper motions, together with our new distance determination and recent literature radial velocities, we are able to show that the two sample clusters have typical halo orbits that are passing by the bulge/bar region, but that they are not part of this component. For the first time, halo intruders are identified in the bulge.

Polarimetry of the Eclipsing Pulsar PSR J1748–2446A

Abstract Observations with the Parkes radio telescope of the eclipsing millisecond binary pulsar PSR J1748−2446A, which is in the globular cluster Terzan 5, are presented. These include the first observations of this pulsar in the 3 GHz frequency band, along with simultaneous observations in the 700 MHz band, and new observations around 1400 MHz. We show that the pulsar signal is not eclipsed in the 3 GHz band and observe the known eclipses in the lower-frequency bands. We find that the observed pulse signal becomes depolarized during particular orbital phases and postulate that this depolarization occurs because of rotation-measure fluctuations resulting from turbulence in the stellar wind responsible for the eclipses.