Subtraction-noise projection in gravitational-wave detector networks

In this paper, we present a successful implementation of a subtraction-noise projection method into a simple, simulated data analysis pipeline of a gravitational-wave search. We investigate the problem to reveal a weak stochastic background signal which is covered by a strong foreground of compact-binary coalescences. The foreground which is estimated by matched filters, has to be subtracted from the data. Even an optimal analysis of foreground signals will leave subtraction noise due to estimation errors of template parameters which may corrupt the measurement of the background signal. The subtraction noise can be removed by a noise projection. We apply our analysis pipeline to the proposed future-generation space-borne Big Bang Observer (BBO) mission which seeks for a stochastic background of primordial GWs in the frequency range $\sim 0.1-1$Hz covered by a foreground of black-hole and neutron-star binaries. Our analysis is based on a simulation code which provides a dynamical model of a time-delay interferometer (TDI) network. It generates the data as time series and incorporates the analysis pipeline together with the noise projection. Our results confirm previous ad hoc predictions which say that BBO will be sensitive to backgrounds with fractional energy densities below $\Omega=10^{-16}$Comment: 54 pages, 15 figure

BBO and the Neutron-Star-Binary Subtraction Problem

The Big Bang Observer (BBO) is a proposed space-based gravitational-wave (GW) mission designed primarily to search for an inflation-generated GW background in the frequency range 0.1-1 Hz. The major astrophysical foreground in this range is gravitational radiation from inspiraling compact binaries. This foreground is expected to be much larger than the inflation-generated background, so to accomplish its main goal, BBO must be sensitive enough to identify and subtract out practically all such binaries in the observable universe. It is somewhat subtle to decide whether BBO's current baseline design is sufficiently sensitive for this task, since, at least initially, the dominant noise source impeding identification of any one binary is confusion noise from all the others. Here we present a self-consistent scheme for deciding whether BBO's baseline design is indeed adequate for subtracting out the binary foreground. We conclude that the current baseline should be sufficient. However if BBO's instrumental sensitivity were degraded by a factor 2-4, it could no longer perform its main mission. It is impossible to perfectly subtract out each of the binary inspiral waveforms, so an important question is how to deal with the "residual" errors in the post-subtraction data stream. We sketch a strategy of "projecting out" these residual errors, at the cost of some effective bandwidth. We also provide estimates of the sizes of various post-Newtonian effects in the inspiral waveforms that must be accounted for in the BBO analysis.Comment: corrects some errors in figure captions that are present in the published versio

The confinement of phonon propagation in TiAlN/Ag multilayer coatings with anomalously low heat conductivity

TiAlN/Ag multilayer coatings with a different number of bilayers and thicknesses of individual layers were fabricated by DC magnetron co-sputtering. Thermal conductivity was measured in dependence of Ag layer thickness. It was found anomalous low thermal conductivity of silver comparing to TiAlN and Ag bulk standards and TiAlN/TiN multilayers. The physical nature of such thermal barrier properties of the multilayer coatings was explained on the basis of reflection electron energy loss spectroscopy. The analysis shows that nanostructuring of the coating decreases the density of states and velocity of acoustic phonons propagation. At the same time, multiphonon channels of heat propagation degenerate. These results demonstrate that metal-dielectric interfaces in TiAlN/Ag coatings are insurmountable obstacles for acoustic phonons propagation

Predictability of band-limited, high-frequency, and mixed processes in the presence of ideal low-pass filters

Pathwise predictability of continuous time processes is studied in deterministic setting. We discuss uniform prediction in some weak sense with respect to certain classes of inputs. More precisely, we study possibility of approximation of convolution integrals over future time by integrals over past time. We found that all band-limited processes are predictable in this sense, as well as high-frequency processes with zero energy at low frequencies. It follows that a process of mixed type still can be predicted if an ideal low-pass filter exists for this process.Comment: 10 page

Pulsar timing arrays as imaging gravitational wave telescopes: angular resolution and source (de)confusion

Pulsar timing arrays (PTAs) will be sensitive to a finite number of gravitational wave (GW) "point" sources (e.g. supermassive black hole binaries). N quiet pulsars with accurately known distances d_{pulsar} can characterize up to 2N/7 distant chirping sources per frequency bin \Delta f_{gw}=1/T, and localize them with "diffraction limited" precision \delta\theta \gtrsim (1/SNR)(\lambda_{gw}/d_{pulsar}). Even if the pulsar distances are poorly known, a PTA with F frequency bins can still characterize up to (2N/7)[1-(1/2F)] sources per bin, and the quasi-singular pattern of timing residuals in the vicinity of a GW source still allows the source to be localized quasi-topologically within roughly the smallest quadrilateral of quiet pulsars that encircles it on the sky, down to a limiting resolution \delta\theta \gtrsim (1/SNR) \sqrt{\lambda_{gw}/d_{pulsar}}. PTAs may be unconfused, even at the lowest frequencies, with matched filtering always appropriate.Comment: 7 pages, 1 figure, matches Phys.Rev.D versio

Particle Swarm Optimization and gravitational wave data analysis: Performance on a binary inspiral testbed

The detection and estimation of gravitational wave (GW) signals belonging to a parameterized family of waveforms requires, in general, the numerical maximization of a data-dependent function of the signal parameters. Due to noise in the data, the function to be maximized is often highly multi-modal with numerous local maxima. Searching for the global maximum then becomes computationally expensive, which in turn can limit the scientific scope of the search. Stochastic optimization is one possible approach to reducing computational costs in such applications. We report results from a first investigation of the Particle Swarm Optimization (PSO) method in this context. The method is applied to a testbed motivated by the problem of detection and estimation of a binary inspiral signal. Our results show that PSO works well in the presence of high multi-modality, making it a viable candidate method for further applications in GW data analysis.Comment: 13 pages, 5 figure

Parameter estimation of coalescing supermassive black hole binaries with LISA

Laser Interferometer Space Antenna (LISA) will routinely observe coalescences of supermassive black hole (BH) binaries up to very high redshifts. LISA can measure mass parameters of such coalescences to a relative accuracy of $10^{-4}-10^{-6}$, for sources at a distance of 3 Gpc. The problem of parameter estimation of massive nonspinning binary black holes using post-Newtonian (PN) phasing formula is studied in the context of LISA. Specifically, the performance of the 3.5PN templates is contrasted against its 2PN counterpart using a waveform which is averaged over the LISA pattern functions. The improvement due to the higher order corrections to the phasing formula is examined by calculating the errors in the estimation of mass parameters at each order. The estimation of the mass parameters ${\cal M}$ and $\eta$ are significantly enhanced by using the 3.5PN waveform instead of the 2PN one. For an equal mass binary of $2\times10^6M_\odot$ at a luminosity distance of 3 Gpc, the improvement in chirp mass is $\sim 11$ and that of $\eta$ is $\sim 39$. Estimation of coalescence time $t_c$ worsens by 43%. The improvement is larger for the unequal mass binary mergers. These results are compared to the ones obtained using a non-pattern averaged waveform. The errors depend very much on the location and orientation of the source and general conclusions cannot be drawn without performing Monte Carlo simulations. Finally the effect of the choice of the lower frequency cut-off for LISA on the parameter estimation is studied.Comment: 12 pages, 5 figures (eps) significant revision, accepted for publication in Phys. Rev. D. Matches with the published versio

The Effect of the LISA Response Function on Observations of Monochromatic Sources

The Laser Interferometer Space Antenna (LISA) is expected to provide the largest observational sample of binary systems of faint sub-solar mass compact objects, in particular white-dwarfs, whose radiation is monochromatic over most of the LISA observational window. Current astrophysical estimates suggest that the instrument will be able to resolve about 10000 such systems, with a large fraction of them at frequencies above 3 mHz, where the wavelength of gravitational waves becomes comparable to or shorter than the LISA arm-length. This affects the structure of the so-called LISA transfer function which cannot be treated as constant in this frequency range: it introduces characteristic phase and amplitude modulations that depend on the source location in the sky and the emission frequency. Here we investigate the effect of the LISA transfer function on detection and parameter estimation for monochromatic sources. For signal detection we show that filters constructed by approximating the transfer function as a constant (long wavelength approximation) introduce a negligible loss of signal-to-noise ratio -- the fitting factor always exceeds 0.97 -- for f below 10mHz, therefore in a frequency range where one would actually expect the approximation to fail. For parameter estimation, we conclude that in the range 3mHz to 30mHz the errors associated with parameter measurements differ from about 5% up to a factor of 10 (depending on the actual source parameters and emission frequency) with respect to those computed using the long wavelength approximation.Comment: replacement version with typos correcte

Interface-Induced Plasmon Nonhomogeneity in Nanostructured Metal-Dielectric Planar Metamaterial

Transformations of the electronic structure in thin silver layers in metal-dielectric (TiAlN/Ag) multilayer nanocomposite were investigated by a set of electron spectroscopy techniques. Localization of the electronic states in the valence band and reduction of electron concentration in the conduction band was observed. This led to decreasing metallic properties of silver in the thin films. A critical layer thickness of 23.5 nm associated with the development of quantum effects was determined by X-ray photoelectron spectroscopy. Scanning Auger electron microscopy of characteristic energy losses provided images of plasmon localization in the Ag layers. The nonuniformity of plasmon intensities distribution near the metal-nitride interfaces was assessed experimentally

Practical Methods for Continuous Gravitational Wave Detection using Pulsar Timing Data

Gravitational Waves (GWs) are tiny ripples in the fabric of space-time predicted by Einstein's General Relativity. Pulsar timing arrays (PTAs) are well poised to detect low frequency ($10^{-9}$ -- $10^{-7}$ Hz) GWs in the near future. There has been a significant amount of research into the detection of a stochastic background of GWs from supermassive black hole binaries (SMBHBs). Recent work has shown that single continuous sources standing out above the background may be detectable by PTAs operating at a sensitivity sufficient to detect the stochastic background. The most likely sources of continuous GWs in the pulsar timing frequency band are extremely massive and/or nearby SMBHBs. In this paper we present detection strategies including various forms of matched filtering and power spectral summing. We determine the efficacy and computational cost of such strategies. It is shown that it is computationally infeasible to use an optimal matched filter including the poorly constrained pulsar distances with a grid based method. We show that an Earth-term-matched filter constructed using only the correlated signal terms is both computationally viable and highly sensitive to GW signals. This technique is only a factor of two less sensitive than the computationally unrealizable optimal matched filter and a factor of two more sensitive than a power spectral summing technique. We further show that a pairwise matched filter, taking the pulsar distances into account is comparable to the optimal matched filter for the single template case and comparable to the Earth-term-matched filter for many search templates. Finally, using simulated data optimal quality, we place a theoretical minimum detectable strain amplitude of $h>2\times 10^{-15}$ from continuous GWs at frequencies on the order $\sim1/T_{\rm obs}$.Comment: submitted to Ap