A Deep Pulse Search in Eleven Low Mass X-Ray Binaries

We present a systematic coherent X-ray pulsation search in eleven low mass X-ray binaries (LMXBs). We select a relatively broad variety of LMXBs, including persistent and transient sources and spanning orbital periods between 0.3 and 17 hours. We use about 3.6 Ms of data collected by the Rossi X-Ray Timing Explorer (RXTE) and XMM-Newton and apply a semi-coherent search strategy to look for weak and persistent pulses in a wide spin frequency range. We find no evidence for X-ray pulsations in these systems and consequently set upper limits on the pulsed sinusoidal semi-amplitude between 0.14% and 0.78% for ten outbursting/persistent LMXBs and 2.9% for a quiescent system. These results suggest that weak pulsations might not form in (most) non-pulsating LMXBs.Comment: submitted to ApJ, 8 page

SWIGLAL: Python and Octave interfaces to the LALSuite gravitational-wave data analysis libraries

The LALSuite data analysis libraries, written in C, implement key routines critical to the successful detection of gravitational waves, such as the template waveforms describing the merger of two black holes or two neutron stars. SWIGLAL is a component of LALSuite which provides interfaces for Python and Octave, making LALSuite routines accessible directly from scripts written in those languages. It has enabled modern gravitational-wave data analysis software, used in the first detection of gravitational waves, to be written in Python, thereby benefiting from its ease of development and rich feature set, while still having access to the computational speed and scientific trustworthiness of the routines provided by LALSuite

Using generalized PowerFlux methods to estimate the parameters of periodic gravitational waves

We investigate methods to estimate the parameters of the gravitational-wave signal from a spinning neutron star using Fourier transformed segments of the strain response from an interferometric detector. Estimating the parameters from the power, we find generalizations of the PowerFlux method. Using simulated elliptically polarized signals injected into Gaussian noise, we apply the generalized methods to estimate the squared amplitudes of the plus and cross polarizations (and, in the most general case, the polarization angle), and test the relative detection efficiencies of the various methods.Comment: 8 pages, presented at Amalid7, Sydney, Australia (July 2007), fixed minor typos and clarified discussion to match published CQG version; updated reference

Proceedings of the Fifth NASA/NSF/DOD Workshop on Aerospace Computational Control

The Fifth Annual Workshop on Aerospace Computational Control was one in a series of workshops sponsored by NASA, NSF, and the DOD. The purpose of these workshops is to address computational issues in the analysis, design, and testing of flexible multibody control systems for aerospace applications. The intention in holding these workshops is to bring together users, researchers, and developers of computational tools in aerospace systems (spacecraft, space robotics, aerospace transportation vehicles, etc.) for the purpose of exchanging ideas on the state of the art in computational tools and techniques

Implementation and characterization of BinaryWeave: A new search pipeline for continuous gravitational waves from Scorpius X-1

Scorpius X-1 (Sco X-1) has long been considered one of the most promising targets for detecting continuous gravitational waves with ground-based detectors. Observational searches for Sco X-1 have achieved substantial sensitivity improvements in recent years, to the point of starting to rule out emission at the torque-balance limit in the low-frequency range \sim 40--180 Hz. In order to further enhance the detection probability, however, there is still much ground to cover for the full range of plausible signal frequencies \sim 20--1500 Hz, as well as a wider range of uncertainties in binary orbital parameters. Motivated by this challenge, we have developed BinaryWeave, a new search pipeline for continuous waves from a neutron star in a known binary system such as Sco X-1. This pipeline employs a semi-coherent StackSlide F-statistic using efficient lattice-based metric template banks, which can cover wide ranges in frequency and unknown orbital parameters. We present a detailed timing model and extensive injection-and-recovery simulations that illustrate that the pipeline can achieve high detection sensitivities over a significant portion of the parameter space when assuming sufficiently large (but realistic) computing budgets. Our studies further underline the need for stricter constraints on the Sco X-1 orbital parameters from electromagnetic observations, in order to be able to push sensitivity below the torque-balance limit over the entire range of possible source parameters

Deep searches for X-ray pulsations from Scorpius X-1 and Cygnus X-2 in support of continuous gravitational wave searches

Neutron stars in low mass X-ray binaries are hypothesised to emit continuous gravitational waves that may be detectable by ground-based observatories. The torque balance model predicts that a higher accretion rate produces larger-amplitude gravitational waves, hence low mass X-ray binaries with high X-ray flux are promising targets for gravitational wave searches. The detection of X-ray pulsations would identify the spin frequency of these neutron stars, and thereby improve the sensitivity of continuous gravitational-wave searches by reducing the volume of the search parameter space. We perform a semi-coherent search for pulsations in the two low mass X-ray binaries Scorpius X-1 and Cygnus X-2 using X-ray data from the \textit{ Rossi X-ray Timing Explorer} Proportional Counter Array. We find no clear evidence for pulsations, and obtain upper limits (at $90\%$ confidence) on the fractional pulse amplitude, with the most stringent being $0.034\%$ for Scorpius X-1 and $0.23\%$ for Cygnus X-2. These upper limits improve upon those of Vaughan et al. (1994) by factors of $\sim 8.2$ and $\sim 1.6$ respectively.Comment: 10 pages, 11 figure

Valency of rare earths in RIn3 and RSn3: Ab initio analysis of electric-field gradients

In RIn3 and RSn3 the rare earth (R) is trivalent, except for Eu and Yb, which are divalent. This was experimentally determined in 1977 by perturbed angular correlation measurements of the electric-field gradient on a 111Cd impurity. At that time, the data were interpreted using a point charge model, which is now known to be unphysical and unreliable. This makes the valency determination potentially questionable. We revisit these data, and analyze them using ab initio calculations of the electric-field gradient. From these calculations, the physical mechanism that is responsible for the influence of the valency on the electric-field gradient is derived. A generally applicable scheme to interpret electric-field gradients is used, which in a transparent way correlates the size of the field gradient with chemical properties of the system.Comment: 10 page

Template Lattices for a Cross-Correlation Search for Gravitational Waves from Scorpius X-1

We describe the application of the lattice covering problem to the placement of templates in a search for continuous gravitational waves from the low-mass X-Ray binary Scorpius X-1. Efficient placement of templates to cover the parameter space at a given maximum mismatch is an application of the sphere covering problem, for which an implementation is available in the LatticeTiling software library. In the case of Sco X-1, potential correlations, in both the prior uncertainty and the mismatch metric, between the orbital period and orbital phase, lead to complications in the efficient construction of the lattice. We define a shearing coordinate transformation which simultaneously minimizes both of these sources of correlation, and allows us to take advantage of the small prior orbital period uncertainty. The resulting lattices have a factor of about 3 fewer templates than the corresponding parameter space grids constructed by the prior straightforward method, allowing a more sensitive search at the same computing cost and maximum mismatch.Comment: 21 pages, 8 figure

N-dimensional electron in a spherical potential: the large-N limit

We show that the energy levels predicted by a 1/N-expansion method for an N-dimensional Hydrogen atom in a spherical potential are always lower than the exact energy levels but monotonically converge towards their exact eigenstates for higher ordered corrections. The technique allows a systematic approach for quantum many body problems in a confined potential and explains the remarkable agreement of such approximate theories when compared to the exact numerical spectrum.Comment: 8 pages, 1 figur

Searching for Galactic White Dwarf Binaries in Mock LISA Data using an F-Statistic Template Bank

We describe an F-statistic search for continuous gravitational waves from galactic white-dwarf binaries in simulated LISA Data. Our search method employs a hierarchical template-grid based exploration of the parameter space. In the first stage, candidate sources are identified in searches using different simulated laser signal combinations (known as TDI variables). Since each source generates a primary maximum near its true "Doppler parameters" (intrinsic frequency and sky position) as well as numerous secondary maxima of the F-statistic in Doppler parameter space, a search for multiple sources needs to distinguish between true signals and secondary maxima associated with other, "louder" signals. Our method does this by applying a coincidence test to reject candidates which are not found at nearby parameter space positions in searches using each of the three TDI variables. For signals surviving the coincidence test, we perform a fully coherent search over a refined parameter grid to provide an accurate parameter estimation for the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters with errors comparable to statistical expectations, although there is still some systematic excess with respect to statistical errors expected from Gaussian noise. An experimental iterative pipeline with seven rounds of signal subtraction and re-analysis of the residuals allows us to increase the number of signals recovered to a total of 3419 with 29 false alarms.Comment: 29 pages, 11 figures; submitted to Classical and Quantum Gravit