While kilohertz quasi-periodic oscillations (kHz QPOs) have been well studied for decades since their initial discovery, the cause of these signals remains unknown, as no model has been able to accurately predict all of their spectral and timing properties. Separately, X-ray reverberation lags have been detected in AGN and stellar-mass black hole binaries, and reverberation may be expected to occur in neutron star systems as well, producing lags of the same amplitude as the lags measured of the kHz QPOs. Furthermore, the detection of a relativistically reflected Fe K line in the time-averaged spectra of many neutron star systems provides an additional motivation for testing reverberation. While it has been shown that the lag-energy properties of the lower kHz QPOs are unlikely to be produced by X-ray reverberation, the upper kHz QPOs have not yet been explored. We therefore model the upper kHz QPO lag-energy spectra using relativistic ray-tracing functions and apply them to archival RXTE data on 4U 1728-34 where upper kHz QPOs have been detected. By modeling the time-averaged spectra in which upper kHz QPOs had been significantly detected, we determine the reflected flux fraction across all energies and produce a model for the lag-energy spectra from X-ray reverberation. We explore the dependence of the modeled lag properties on several different types of reflection models, but are unable to successfully reproduce the measured lags of 4U 1728-34. We conclude that reverberation alone does not explain the measured time lags detected in upper kHz QPOs.

Atoll sources are accreting neutron star (NS) low-mass X-ray binaries. We present a spectral analysis of four persistent atoll sources (GX 3+1, 4U 1702$-$429, 4U 0614+091, and 4U 1746$-$371) observed for $\sim20$ ks each with NuSTAR to determine the extent of the inner accretion disk. These sources range from an apparent luminosity of $0.006-0.11$ of the Eddington limit (assuming the empirical limit of $3.8\times10^{38}$ ergs s$^{-1}$). Broad Fe emission features shaped by Doppler and relativistic effects close to the NS were firmly detected in three of these sources. The position of the disk appears to be close to the innermost stable circular orbit (ISCO) in each case. For GX 3+1, we determine $R_{in}=1.8^{+0.2}_{-0.6}\ R_{\mathrm{ISCO}}$ (90% confidence level) and an inclination of $27^{\circ}-31^{\circ}$. For 4U 1702$-$429, we find a $R_{in}=1.5_{-0.4}^{+1.6}\ R_{\mathrm{ISCO}}$ and inclination of $53^{\circ}-64^{\circ}$. For 4U 0614+091, the disk has a position of $R_{in}=1.3_{-0.2}^{+5.4}\ R_{\mathrm{ISCO}}$ and inclination of $50^{\circ}-62^{\circ}$. If the disk does not extend to the innermost stable circular orbit, we can place conservative limits on the magnetic field strength in these systems in the event that the disk is truncated at the Alfv\'{e}n radius. This provides the limit at the poles of $B\leq6.7\times10^{8}$ G, $3.3\times10^{8}$ G, and $14.5\times10^{8}$ G for GX 3+1, 4U 1702$-$429, and 4U 0614+091, respectively. For 4U 1746$-$371, we argue that the most plausible explanation for the lack of reflection features is a combination of source geometry and strong Comptonization. We place these sources among the larger sample of NSs that have been observed with NuSTAR.

Although the most luminous class of neutron star low mass X-ray binaries, known as Z sources, have been well studied, their behavior is not fully understood. In particular, what causes these sources to trace out the characteristic Z-shaped pattern on color-color or hardness-intensity diagrams is not well known. By studying the physical properties of the different spectral states of these sources, we may better understand such variability. With that goal in mind, we present a recent NuSTAR observation of the Z source GX 349+2, which spans approximately 2 days, and covers all its spectral states. By creating a hardness-intensity diagram we were able to extract five spectra and trace the change in spectral parameters throughout the Z-track. GX 349+2 shows a strong, broad Fe K$\alpha$ line in all states, regardless of the continuum model used. Through modeling of the reflection spectrum and Fe K$\alpha$ line we find that in most states the inner disk radius is consistent with remaining unchanged at an average radius of 17.5 $R_g$ or 36.4 km for a canonical 1.4 $M_\odot$ neutron star. During the brightest flaring branch, however, the inner disk radius from reflection is not well constrained.