Gamma-ray binaries consist of a neutron star or a black hole interacting with a normal star to produce gamma-ray emission that dominates the radiative output of the system. Only a handful of such systems have been previously discovered, all within our Galaxy. Here we report the discovery with the Fermi Large Area Telescope (LAT) of a luminous gamma-ray binary in the Large Magellanic Cloud from a search for periodic modulation in all sources in the third Fermi LAT catalog. This is the first such system to be found outside the Milky Way. The system has an orbital period of 10.3 days and is associated with a massive O5III star located in the supernova remnant DEM L241, previously identified as the candidate high-mass X-ray binary (HMXB) CXOU J053600.0-673507. X-ray and radio emission are also modulated on the 10.3 day period, but are in anti-phase with the gamma-ray modulation. Optical radial velocity measurements suggest that the system contains a neutron star. The source is significantly more luminous than similar sources in the Milky Way at radio, optical, X-ray and gamma-ray wavelengths. The detection of this extra-galactic system, but no new Galactic systems raises the possibility that the predicted number of gamma-ray binaries in our Galaxy have been overestimated, and that HMXBs may be born containing relatively slowly rotating neutron stars.

Physical parameters of both the mass donor and compact object can be constrained in X-ray binaries with well-defined eclipses, as our survey of wind-fed supergiant X-ray binaries (SGXBs) IGR J16393-4643, IGR J16418-4532, IGR J16479-4514, IGR J18027-2016 and XTE J1855-026 reveals. Using the orbital period and Kepler's third law, we express the eclipse half-angle in terms of radius, inclination angle and the sum of the masses. Pulse-timing and radial velocity curves can give masses of both the donor and compact object as in the case of the "double-lined" binaries IGR J18027-2016 and XTE J1855-026. The eclipse half angles are 15 +3 −2 , 31.7 +0.7

We report the identification from multi-wavelength observations of the Fermi Large Area Telescope (LAT) source 4FGL J1405.1-6119 (= 3FGL J1405.4-6119) as a high-mass gamma-ray binary. Observations with the LAT show that gamma-ray emission from the system is modulated at a period of 13.7135 ± 0.0019 days, with the presence of two maxima per orbit with different spectral properties. X-ray observations using the Neil Gehrels Swift Observatory X-ray Telescope (XRT) show that X-ray emission is also modulated at this period, but with a single maximum that is closer to the secondary lower-energy gamma-ray maximum. A radio source, coincident with the X-ray source is also found from Australia Telescope Compact Array (ATCA) observations, and the radio emission is modulated on the gamma-ray period with similar phasing to the X-ray emission. A large degree of interstellar obscuration severely hampers optical observations, but a near-infrared counterpart is found. Nearinfrared spectroscopy indicates an O6 III spectral classification. This is the third gamma-ray binary to be discovered with the Fermi LAT from periodic modulation of the gamma-ray emission, the other two sources also have early O star, rather than Be star, counterparts. We consider at what distances we can detect such modulated gamma-ray emission with the LAT, and examine constraints on the gamma-ray binary population of the Milky Way.

IGR J16493-4348 is an eclipsing supergiant high-mass X-ray binary (sgHMXB), where accretion onto the compact object occurs via the radially outflowing stellar wind of its early B-type companion. We present an analysis of the system's X-ray variability and periodic modulation using pointed observations (2.5-25 keV) and Galactic bulge scans (2-10 keV) from the Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA), along with Swift Burst Alert Telescope (BAT) 70-month snapshot (14-195 keV) and transient monitor (15-50 keV) observations. The orbital eclipse profiles in the PCA bulge scans and BAT light curves are modeled using asymmetric and symmetric step and ramp functions. We obtain an improved orbital period measurement of 6.7828 ± 0.0004 days from an observed minus calculated (O-C) analysis of mid-eclipse times derived from the BAT transient monitor and PCA scan data. No evidence is found for the presence of a strong photoionization or accretion wake. We refine the superorbital period to 20.067 ± 0.009 days from the discrete Fourier transform (DFT) of the BAT transient monitor light curve. A pulse period of 1093.1036 ± 0.0004 s is measured from a pulsar timing analysis using pointed PCA observations spanning ∼1.4 binary orbits. We present pulse times of arrival (ToAs), circular and eccentric timing models, and calculations of the system's Keplerian binary orbital parameters. We derive an X-ray mass function of f x (M ) = 13.2 +2.4 -2.5 M and find a spectral type of B0.5 Ia for the supergiant companion through constraints on the mass and radius of the donor.