We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

We report on the optical identification of the companion star to the eclipsing millisecond pulsar PSR J1824-2452H in the galactic globular cluster M28 (NGC 6626). This star is at only 0.2" from the nominal position of the pulsar and it shows optical variability (~ 0.25 mag) that nicely correlates with the pulsar orbital period. It is located on the blue side of the cluster main sequence, ~1.5 mag fainter than the turn-off point. The observed light curve shows two distinct and asymmetric minima, suggesting that the companion star is suffering tidal distortion from the pulsar. This discovery increases the number of non-degenerate MSP companions optically identified so far in globular clusters (4 out of 7), suggesting that these systems could be a common outcome of the pulsar recycling process, at least in dense environments where they can be originated by exchange interactions.

We report the discovery of five new millisecond pulsars in the globular cluster NGC 6440 and three new ones in NGC 6441; each cluster has one previously known pulsar. Four of the new pulsars are found in binary systems. One of the new pulsars, PSR J1748-2021B in NGC 6440, is notable for its eccentric (e = 0.57) and wide (P_b = 20.5 days) orbit. If the rate of advance of periastron is due solely to general relativity, we can estimate of the total mass of this binary system: 2.92 +/- 0.20 solar masses. This would imply an anomalously large mass for this pulsar, which could introduce important constraints in the study of the equation of state for cold neutron matter.

Motivated by the recent discovery of 30 new millisecond pulsars in Terzan 5, made using the Green Bank Telescope's S-band receiver and the Pulsar Spigot spectrometer, we have set out to use the same observing system in a systematic search for pulsars in other globular clusters. Here we report on the discovery of five new pulsars in NGC 6440 and three in NGC 6441; each cluster previously had one known pulsar. Using the most recent distance estimates to these clusters, we conclude that there are as many potentially observable pulsars in NGC 6440 and NGC 6441 as in Terzan 5. We present timing solutions for all of the pulsars in these globular clusters. Four of the new discoveries are in binary systems; one of them, PSR J1748-2021B (NGC 6440B), has a wide (P_b = 20.5 d) and eccentric (e = 0.57) orbit. This allowed a measurement of its rate of advance of periastron: 0.00391(18) degrees per year. If due to the effects of general relativity, the total mass of this binary system is 2.92 +/- 0.20 solar masses (1 sigma), implying a median pulsar mass of 2.74 +/- 0.21 solar masses. There is a 1 % probability that the inclination is low enough that pulsar mass is below 2 solar masses, and 0.10 % probability that it is between 1.20 and 1.44 solar masses. If confirmed, this anomalously large mass would strongly constrain the equation of state for dense matter. The other highly eccentric binary, PSR J1750-37A, has e = 0.71, and periastron advance of 0.0055(3) degrees per year, implying a total system mass of 1.97 +/-0.15 solar masses and, along with the mass function, maximum and median pulsar masses of 1.65 and 1.26 solar masses respectively.