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. Near-infrared 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.

We report on Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory (Swift) X-ray Telescope (XRT) and Swift Burst Alert Telescope (BAT) observations of IGR J16493-4348, a wind-fed Supergiant X-ray Binary (SGXB) showing significant superorbital variability. From a discrete Fourier transform of the BAT light curve, we refine its superorbital period to be 20.058 $\pm$ 0.007 days. The BAT dynamic power spectrum and a fractional root mean square analysis both show strong variations in the amplitude of the superorbital modulation, but no observed changes in the period were found. The superorbital modulation is significantly weaker between MJD 55,700 and MJD 56,300. The joint NuSTAR and XRT observations, which were performed near the minimum and maximum of one cycle of the 20 day superorbital modulation, show that the flux increases by more than a factor of two between superorbital minimum and maximum. We find no significant changes in the 3-50 keV pulse profiles between superorbital minimum and maximum, which suggests a similar accretion regime. Modeling the pulse-phase averaged spectra we find a possible Fe K$\alpha$ emission line at 6.4 keV at superorbital maximum. The feature is not significant at superorbital minimum. While we do not observe any significant differences between the pulse-phase averaged spectral continua apart from the overall flux change, we find that the hardness ratio near the broad main peak of the pulse profile increases from superorbital minimum to maximum. This suggests the spectral shape hardens with increasing luminosity. We discuss different mechanisms that might drive the observed superorbital modulation.

Studying the physical processes occurring in the region just above the magnetic poles of strongly magnetized, accreting binary neutron stars is essential to our understanding of stellar and binary system evolution. Perhaps more importantly, it provides us with a natural laboratory for studying the physics of high temperature and high density plasmas exposed to extreme radiation, gravitational, and magnetic fields. Observations over the past decade have shed new light on the manner in which plasma falling at velocities near the speed of light onto a neutron star surface is halted. Recent advances in modeling these processes have resulted in direct measurement of the magnetic fields and plasma properties. On the other hand, numerous physical processes have been identified that challenge our current picture of how the accretion process onto neutron stars works. Observation and theory are our essential tools in this regime because the extreme conditions cannot be duplicated on Earth. This white paper gives an overview of the current theory, the outstanding theoretical and observational challenges, and the importance of addressing them in contemporary astrophysics research.

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 $\pm$ 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 $\pm$ 0.009 days from the discrete Fourier transform (DFT) of the BAT transient monitor light curve. A pulse period of 1093.1036 $\pm$ 0.0004 s is measured from a pulsar timing analysis using pointed PCA observations spanning $\sim$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_{\odot}$ and find a spectral type of B0.5 Ia for the supergiant companion through constraints on the mass and radius of the donor. Measurements of the eclipse half-angle and additional parameters describing the system geometry are provided.

We present an investigation of long-term modulation in the X-ray light curves of five little-studied candidate high-mass X-ray binaries using the Swift Burst Alert Telescope. IGR J14488-5942 and AX J1700.2-4220 show strong modulation at periods of 49.6 and 44 days, respectively, which are interpreted as orbital periods of Be star systems. For IGR J14488-5942, observations with Swift X-ray Telescope show a hint of pulsations at 33.4 s. For AX J1700.2-4220, 54 s pulsations were previously found with XMM. Swift J1816.7-1613 exhibits complicated behavior. The strongest peak in the power spectrum is at a period near 150 days, but this conflicts with a determination of a period of 118.5 days by La Parola et al. (2014). AX J1820.5-1434 has been proposed to exhibit modulation near 54 days, but the extended BAT observations suggest modulation at slightly longer than double this at approximately 111 days. There appears to be a long-term change in the shape of the modulation near 111 days, which may explain the apparent discrepancy. The X-ray pulsar XTE J1906+090, which was previously proposed to be a Be star system with an orbital period of ~30 days from pulse timing, shows peaks in the power spectrum at 81 and 173 days. The origins of these periods are unclear, although they might be the orbital period and a superorbital period respectively. For all five sources, the long-term variability, together with the combination of orbital and proposed pulse periods, suggests that the sources contain Be star mass donors.

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 has 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,-0.8), 32 +/- 2, 34 +/- 2 and 33.6 +/- 0.7 degrees for IGR J16393-4643, IGR J16418-4532, IGR J16479-4514, IGR J18027-2016 and XTE 1855-026, respectively. In wind-fed systems, the primary not exceeding the Roche-lobe size provides an upper limit on system parameters. In IGR J16393-4643, spectral types of B0 V or B0-5 III are found to be consistent with the eclipse duration and Roche-lobe, but the previously proposed donor stars in IGR J16418-4532 and IGR J16479-4514 were found to be inconsistent with the Roche-lobe size. Stars with spectral types O7.5 I and earlier are possible. For IGR J18027-2016, the mass and radius of the donor star lie between 18.6--19.4 $M_\odot$ and 17.4--19.5 $R_\odot$. We constrain the neutron star mass between 1.37--1.43 $M_\odot$. We find the mass and radius of the donor star in XTE J1855-026 to lie between 19.6--20.2 $M_\odot$ and 21.5--23.0 $R_\odot$. The neutron star mass was constrained to 1.77--1.82 $M_\odot$. Eclipse profiles are asymmetric in IGR J18027-2016 and XTE J1855-026, which we attribute to accretion wakes.

4U 1210-64 has been postulated to be a High-Mass X-ray Binary powered by the Be mechanism. X-ray observations with Suzaku, the ISS Monitor of All-sky X-ray Image (MAXI) and the Rossi X-ray Timing Explorer Proportional Counter Array (PCA) and All Sky Monitor (ASM) provide detailed temporal and spectral information on this poorly understood source. Long term ASM and MAXI observations show distinct high and low states and the presence of a 6.7101 +/- 0.0005 day modulation, interpreted as the orbital period. Folded light curves reveal a sharp dip, interpreted as an eclipse. To determine the nature of the mass donor, the predicted eclipse half-angle was calculated as a function of inclination angle for several stellar spectral types. The eclipse half-angle is not consistent with a mass donor of spectral type B5 V; however, stars with spectral types B0 V or B0-5 III are possible. The best-fit spectral model consists of a power law with index gamma = 1.85 (+0.04,-0.05) and a high-energy cutoff at 5.5 +/- 0.2 keV modified by an absorber that fully covers the source as well as partially covering absorption. Emission lines from S XVI K alpha, Fe K alpha, Fe XXV K alpha and Fe XXVI K alpha were observed in the Suzaku spectra. Out-of-eclipse, Fe K alpha line flux was strongly correlated with unabsorbed continuum flux, indicating that the Fe I emission is the result of fluorescence of cold dense material near the compact object. The Fe I feature is not detected during eclipse, further supporting an origin close to the compact object.