The first gamma ray pulsation at energies above 50 MeV was detected from the Crab pulsar (NP 0532) in 1971 by a detector consisting of a multiplate spark chamber in a balloon-borne experiment
[45][1]. These gamma ray pulsations were observed at the radio period (
๐โ 33 ms) of the source. The launch of the
Small Astronomy Satellite (SAS-2) in 1972 confirmed the pulsed gamma ray emission from the Crab and Vela pulsars
[46,47][2][3]. The
COS-B satellite, launched in 1975, confirmed and significantly improved the measurements of SAS-2 and also discovered a population of twenty unidentified gamma ray sources in the Milky Way Galaxy
[48][4]. Both SAS-2 and COS-B detected pulsed gamma ray signal at low significance level from the known radio pulsars. The number of radio pulsars increased rapidly after their first discovery, but Crab and Vela remained the only known gamma ray pulsars for more than two decades after their first detection. Subsequently, the
EGRET (Energetic Gamma Ray Experiment Telescope), onboard the Compton Gamma Ray Observatory, was launched in 1991 to search for more gamma ray pulsars apart from exploring the high-energy Universe. It detected a total of six gamma ray pulsars (Crab, B1509-58, Vela, B1706-44, B1951+32, Geminga, and B1055-52) and a few plausible candidates during its operation. Discovery of Geminga as a radio-quiet gamma ray pulsar was an important finding of EGRET, as it turned to be the prototype of young gamma ray pulsars
[49,50,51][5][6][7]. Finally, the Large Area Telescope (LAT), onboard the
Fermi satellite, has been able to discover about 340 gamma ray pulsars at energies above 100 MeV in a time span of approximately 15 years since its launch in 2008
[52][8]. It has not only confirmed the existence of radio-quiet gamma ray pulsars and MSPs, but has also established pulsars as dominant high-energy gamma ray sources in the Universe.
2. Pulsars with MACE
The gamma ray spectra of most of the pulsars observed by the
Fermi-LAT have a cut-off energy in the range 1โ10 GeV. Due to this, the gamma ray flux at higher energy end of the spectrum decreases very quickly. Therefore, ground-based telescopes with better sensitivity than the
Fermi-LAT in the energy range above 50 GeV play a crucial role in determining the pulsar spectra in the very-high-energy band. With this motivation,
wresearche
rsย have developed a formalism based on the following arguments to select a sample of the
Fermi-LAT pulsars for their observation with the MACE telescope.
In general, pulsar light curves or phasograms exhibit a two-peak structure. The first peak position at phase
๐1ย is attributed to the radio emission and is therefore referred to as the main pulse. This is assigned zero phase, and the position of the second peak at phase
๐2ย is calculated with respect to
๐1. However, there exist pulsars like Geminga, which have only one peak in their phasogram. The pulse amplitude corresponding to these peaks changes with the energy and the
๐2ย peak is observed to dominate with increasing energy. Therefore, the observation of pulsars in the very-high-energy range by ground-based telescopes like MACE is governed by the detection of pulse at
๐2 at
5๐ย confidence level in their phasogram over the huge isotropic cosmic ray background mainly due to protons. Thus, if a gamma ray pulsar is residing within the nebula, both gamma rays from the nebula and cosmic ray protons will act as the background for pulsar detection with ground-based telescopes. These background events lack in periodicity and, hence, have no phase, whereas photons originating from the pulsars have a well-defined phase associated with their rotation. Thus, the background events are uniformly distributed in the light curve or phasogram of a pulsar, and the gamma ray signal from a pulsar will pile up in a particular phase-bin. If
๐๐๐๐๐ ย is the number of bins in the phasogram of a pulsar, then the background will be reduced by a factor of
๐๐๐๐๐ . For
๐๐๐๐๐ = 30, the observation time required for detection of dominant gamma ray peak (
๐๐) in the phasogram of a pulsar at 5
๐ย statistical significance level with the MACE telescope is estimated as
where
๐ผ๐ย is the integral flux of dominating peak above 30 GeV in units of
ph cmโ2 sโ1.
WResearche
rs apply this methodology to two well-known gamma ray pulsars Crab and Geminga to estimate the expected time for MACE observations as their gamma ray spectra are available in the energy range covered by the MACE telescope.
2.1. Crab Pulsar
The Crab pulsar (PSR J0534+2200) is the most studied pulsar in all energy bands, from radio to TeV, and is also the first pulsar detected by any ground-based gamma ray telescope. It is a normal pulsar with
๐=ย 33 ms and
๐ห=4.2ร10โ13ssโ1. After many attempts, the MAGIC-I telescope detected pulsed emission above 25 GeV from Crab in 2008 and revealed the presence of high cut-off energy in the gamma ray spectrum
[77][9]. Subsequently, the VERITAS collaboration reported pulsed emission from Crab in the energy range 100โ400 GeV
[78][10]. Long-term observations of the Crab pulsar with the MAGIC telescopes resulted in the first ever detection of very high energy gamma ray emission in the energy range 70 GeVโ1.5 TeV
[79][11]. This finding strongly supported the inverse-Compton origin of the gamma ray photons from the Crab pulsar. It is among the few pulsars that have been detected across the entire electromagnetic spectrum from radio up to gamma rays of beyond 1 TeV
[79][11]. Its emission profile has three components: two dominant pulses separated by 0.4 in phase, observed from radio to TeV gamma rays, and a third bridge component. The main pulse corresponding to
๐1ย at phase 0 dominates at radio frequencies. The second pulse corresponding to
๐2ย at phase 0.4 (also referred to as interpulse) is weaker at radio frequencies and becomes dominant at very-high-energy gamma rays. The bridge component (defined as the pulse phase between the main pulse and the second pulse) is observed in gamma rays up to 150 GeV
[79][11]. The stereoscopic MAGIC telescope has detected a signal at 6.6
๐ย and 8.8
๐ย for
๐1ย and
๐2, respectively, in 152 h of crab pulsar observations at zenith angles below 35
โโ. The differential energy spectra of gamma ray photons corresponding to
๐1ย and
๐2ย are described by a power law in the range 70 GeV to 1.5 TeV. The
๐2ย spectrum is relatively hard and the corresponding integral flux is two times that of
๐1. For
๐ผ๐=๐ผ๐2=3.35ร10โ11ph cmโ2sโ1, the corresponding detection time for MACE is estimated to be
๐๐โผ61.38 h using Equation (
24).
2.2. Vela
The Vela pulsar (PSR J0835-4510) with
๐=ย 89.3 ms and
๐ห=1.2ร10โ13ssโ1ย is the second gamma ray pulsar detected, after Crab. It is located in the southern hemisphere at a distance of 280 pc from Earth. The light curve of the Vela pulsar exhibits two peaks,
๐1ย and
๐2, separated by a phase of 0.43. The two peaks are connected by a bridge emission labeled as
๐3. The largest IACT, HESS-II, made the first ground-based detection of pulsed gamma ray emission in the energy range of sub-20 GeV to 100 GeV from the Vela pulsar
[80][12]. The gamma ray signal was detected at the
๐2ย peak with a statistical significance level of more than 15
๐. The power-law model with a very soft spectral index (
ฮ๐ ย โผ4.1) very well describes the differential energy spectrum of the
๐2ย peak in the energy range 10โ100 GeV. Being located in the southern hemisphere, the Vela pulsar will be visible to the MACE telescope only at large zenith angles above 75
โ. At these zenith angles, the threshold energy for MACE may increase beyond 150 GeV. Therefore,
wresearche
rsย have not considered this source as potential candidate for MACE observations in the present work.
2.3. Geminga
The Geminga pulsar (PSR J0633+1746) is a radio-quiet gamma ray pulsar with
๐=ย 237 ms and
๐ห=1.1ร10โ14ssโ1.ย Located at a distance of 250 pc from Earth, it is one of the closest gamma ray pulsars. The MAGIC telescope has detected a pulsed gamma ray signal from Geminga at 6.3
๐ย statistical significance level in 80 h of observations
[81][13]. Only the
๐2ย peak is observed to be dominant in the phasogram in the energy range of 15โ75 GeV. The corresponding integral flux is estimated as
๐ผ๐=๐ผ๐2=2.18ร10โ11ph cmโ2sโ1. For this flux level, the total observation time required for MACE observations in the energy range 30 GeVโ5 TeV is calculated to be
๐๐โผ144 h.
2.4. Pulsar Catalog for MACE
The current population of gamma ray pulsars detected by the
Fermi-LAT is โผ340. The distribution of integral energy flux in the energy range 0.1โ100 GeV for all the
Fermi-LAT pulsars is shown in
Figure 1. It is observed that majority of the pulsars have high-energy flux โผ
1.4ร10โ11 erg cmโ2sโ1, which is typically two order less than the corresponding Crab pulsar flux. The
Fermi-LAT spectra of most of the pulsars are described by Equation (
20). The cut-off energy can be approximately calculated as
Figure 1. Histogram of energy flux of the Fermi-LAT pulsars in the energy range 0.1โ100 GeV.
The distribution of the cut-off energy (
๐ธ๐) for all the
Fermi-LAT pulsars is reported in
Figure 2. It is found that most of the pulsars have a cut-off energy close to 0.8 GeV, whereas only 33 pulsars have
๐ธ๐โฅย 2 GeV.
Figure 2. Histogram of the cut-off energy of the Fermi-LAT pulsars. The vertical dotted line indicates ๐ธ๐=2 GeV.
Applying the visibility criteria of the astrophysical sources from a given observatory, only 96
Fermi-LAT pulsars are found to be visible from the MACE site at Hanle (32.8
โ N, 78.9
o E). The next criterion is the identification of pulsars which have
๐ธ๐โฅย 2 GeV. This constraint further reduces the number of pulsars from 96 to 8. Lastly,
wresearche
rs check for the significance of pulsation in the energy range 10โ30 GeV, overlapping with MACE, to be more than 2
๐. Based on the above criterion, researchers are left with only six pulsar candidates (other than Crab and Geminga) from the 4FGL-DR3 catalog for MACE observations.