13092497
Application of asymptotic expansions for maximum likelihood estimators errors to gravitational waves from binary mergers: The single interferometer case
In this paper we apply to gravitational waves (GW) from the inspiral phase of binary systems a recently derived frequentist methodology to calculate analytically the error for a maximum likelihood estimate of physical parameters. We use expansions of the covariance and the bias of a maximum likelihood estimate in terms of inverse powers of the signal-to-noise ration (SNR)s where the square root of the first order in the covariance expansion is the Cramer Rao lower bound (CRLB). We evaluate the expansions, for the first time, for GW signals in noises of GW interferometers. The examples are limited to a single, optimally oriented, interferometer. We also compare the error estimates using the first two orders of the expansions with existing numerical Monte Carlo simulations. The first two orders of the covariance allow us to get error predictions closer to what is observed in numerical simulations than the CRLB. The methodology also predicts a necessary SNR to approximate the error with the CRLB and provides new insight on the relationship between waveform properties, SNR, dimension of the parameter space and estimation errors. For example the timing match filtering can achieve the CRLB only if the SNR is larger than the Kurtosis ofmore » the gravitational wave spectrum and the necessary SNR is much larger if other physical parameters are also unknown.« less
17415632
Approximate solutions in General Relativity via deformation of embeddings
A systematic study of deformations of four-dimensional Einsteinian space-times embedded in a pseudo-Euclidean space E N of higher dimension is presented. Infinitesimal deformations, seen as vector fields in E N , can be divided in two parts, tangent to the embedded hypersurface and orthogonal to it; only the second ones are relevant, the tangent ones being equivalent to coordinate transformations in the embedded manifold. The geometrical quantities can be then expressed in terms of embedding functions z A and their infinitesimal deformations v A z A → ˜ z A = z A + e v A . The deformations are called Einsteinian if they keep Einstein equations satisfied up to a given order in e. The system so obtained is then analyzed in particular in the case of the Schwarzschild metric taken as the starting point, and some solutions of the first-order deformation of Einstein’s equations are found. We discuss also second and third order deformations leading to wave-like solutions and to the departure from spherical symmetry towards an axial one (the approximate Kerr solution) 1
118465419
Parameter estimation from gravitational waves generated by nonspinning binary black holes with laser interferometers: Beyond the Fisher information
In this paper we apply to gravitational waves from nonspinning binary systems a recently introduced frequentist methodology to calculate analytically the error for a maximum likelihood estimate of physical parameters. While existing literature focuses on using the Cramer Rao Lower bound and Monte Carlo simulations, we use a power expansion of the bias and covariance in inverse powers of the signal to noise ratio. The use of higher order derivatives of the likelihood function in the expansions makes the prediction also sensitive to the secondary lobes of the maximum likelihood estimate probability distribution. We discuss conditions for validity of the Cramer Rao Lower bound and predict new features in regions of the parameter space currently not explored. For example, we see how the bias can become the most important contributor to the parameters' errors for high mass systems (200M(circle dot) and above).
118385730
Testing general relativity using gravitational waves from binary neutron stars: Effect of spins
We present a Bayesian data analysis pipeline for testing GR using gravitational wave signals from coalescing compact binaries, and in particular binary neutron stars. In this study, we investigate its performance when sources with spins are taken into account.
23581346
Demonstrating the feasibility of probing the neutron-star equation of state with second-generation gravitational-wave detectors.
Fisher matrix and related studies have suggested that, with second-generation gravitational-wave detectors, it may be possible to infer the equation of state of neutron stars using tidal effects in a binary inspiral. Here, we present the first fully Bayesian investigation of this problem. We simulate a realistic data analysis setting by performing a series of numerical experiments of binary neutron-star signals hidden in detector noise, assuming the projected final design sensitivity of the Advanced LIGO-Virgo network. With an astrophysical distribution of events (in particular, uniform in comoving volume), we find that only a few tens of detections will be required to arrive at strong constraints, even for some of the softest equations of state in the literature. Thus, direct gravitational-wave detection will provide a unique probe of neutron-star structure.
118457757
Towards a generic test of the strong field dynamics of general relativity using compact binary coalescence
Coalescences of binary neutron stars and/or black holes are amongst the most likely gravitational-wave signals to be observed in ground based interferometric detectors. Apart from the astrophysical importance of their detection, they will also provide us with our very first empirical access to the genuinely strong-field dynamics of General Relativity (GR). We present a new framework based on Bayesian model selection aimed at detecting deviations from GR, subject to the constraints of the Advanced Virgo and LIGO detectors. The method tests the consistency of coefficients appearing in the waveform with the predictions made by GR, without relying on any specific alternative theory of gravity. The framework is suitable for low signal-to-noise ratio events through the construction of multiple subtests, most of which involve only a limited number of coefficients. It also naturally allows for the combination of information from multiple sources to increase one's confidence in GR or a violation thereof. We expect it to be capable of finding a wide range of possible deviations from GR, including ones which in principle cannot be accommodated by the model waveforms, on condition that the induced change in phase at frequencies where the detectors are the most sensitive is comparable to the effect of a few percent change in one or more of the low-order post-Newtonian phase coefficients. In principle the framework can be used with any GR waveform approximant, with arbitrary parameterized deformations, to serve as model waveforms. In order to illustrate the workings of the method, we perform a range of numerical experiments in which simulated gravitational waves modeled in the restricted post-Newtonian, stationary phase approximation are added to Gaussian and stationary noise that follows the expected Advanced LIGO/Virgo noise curves.
119205701
Discrete derivative estimation in LISA Pathfinder data reduction
Data analysis for the LISA Technology package (LTP) experiment to be flown aboard the LISA Pathfinder mission requires the solution of the system dynamics for the calculation of the force acting on the test masses (TMs) starting from interferometer position data. The need for a solution to this problem has prompted us to implement a discrete time-domain derivative estimator suited for the LTP experiment requirements. We first report on the mathematical procedures for the definition of two methods; the first based on a parabolic fit approximation and the second based on a Taylor series expansion. These two methods are then generalized and incorporated into a more general class of five-point discrete derivative estimators. The same procedure employed for the second derivative can be applied to the estimation of the first derivative and of a data smoother allowing defining a class of simple five-point estimators for both. The performances of three particular realizations of the five-point second-derivative estimator are analyzed with simulated noisy data. This analysis pointed out that those estimators introducing large amount of high-frequency noise can determine systematic errors in the estimation of low-frequency noise levels.
9497849
Interaction between stray electrostatic fields and a charged free-falling test mass.
We present an experimental analysis of force noise caused by stray electrostatic fields acting on a charged test mass inside a conducting enclosure, a key problem for precise gravitational experiments. Measurement of the average field that couples to the test mass charge, and its fluctuations, is performed with two independent torsion pendulum techniques, including direct measurement of the forces caused by a change in electrostatic charge. We analyze the problem with an improved electrostatic model that, coupled with the experimental data, also indicates how to correctly measure and null the stray field that interacts with the test mass charge. Our measurements allow a conservative upper limit on acceleration noise, of 2  (fm/s2)/Hz(1/2) for frequencies above 0.1 mHz, for the interaction between stray fields and charge in the LISA gravitational wave mission.
123492651
Axion to magnon conversion. A scheme for the detection of galactic axions
Abstract We consider the possibility of detecting axions in the halo of our Galaxy by means of converting them into collective excitations of a ferromagnet. We calculate that, under optimal conditions, one can convert about 6000 axions per day in a 1 m 3 ferromagnet. The signal has an intrinsic diurnal modulation of order unit due to the combined effect of Earth rotation and motion of the solar system with respect to the galactic center. We discuss the ultimate limitations arising from thermal fluctuations of the magnetization.
118459662
Bayesian parameter estimation in the second LISA Pathfinder mock data challenge
A main scientific output of the LISA Pathfinder mission is to provide a noise model that can be extended to the future gravitational wave observatory, LISA. The success of the mission depends thus upon a deep understanding of the instrument, especially the ability to correctly determine the parameters of the underlying noise model. In this work we estimate the parameters of a simplified model of the LISA Technology Package instrument. We describe the LISA Technology Package by means of a closed-loop model that is used to generate the data, both injected signals and noise. Then, parameters are estimated using a Bayesian framework, and it is shown that this method reaches the optimal attainable error, the Cramer-Rao bound. We also address an important issue for the mission: how to efficiently combine the results of different experiments to obtain a unique set of parameters describing the instrument.
117779864
Ferromagnetic core SQUIDs: perspectives as detectors of gravitomagnetic fields and exotics.