Using electromagnetic observations to aid gravitationalwave parameter estimation of compact binaries observed with LISA
^{1} Department of Astrophysics/ IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands
email: s.shah@astro.ru.nl
^{2} Nikhef National Institute for Subatomic Physics, Science Park 105, 1098 XG Amsterdam, The Netherlands
^{3} Institute for Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
Received: 29 March 2012
Accepted: 17 June 2012
We present a firststage study of the effect of using knowledge from electromagnetic (EM) observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna (LISA) in the lowfrequency range, 10^{4} Hz < f < 1 Hz. In particular, we examine the extent to which the accuracy of GW parameter estimation improves if we use available information from EM data. We do this by investigating whether correlations exist between the GW parameters that describe these binaries and whether some of these parameters are also available from EM observations. We used verification binaries, which are known as the guaranteed sources for eLISA and will test the functioning of the instrument. We find that of the seven parameters that characterise such a binary, only a few are correlated. The most useful result is the strong correlation between amplitude and inclination, which can be used to constrain the parameter uncertainty in amplitude by making use of the constraint of inclination from EM measurements. The improvement can be up to a factor of ~6.5, but depends on the signaltonoise ratio of the source data. Moreover, we find that this strong correlation depends on the inclination. For mildly faceon binaries (ι ≲ 45°), EM data on inclination can improve the estimate of the GW amplitude by a significant factor. However, for edgeon binaries (ι ~ 90°), the inclination can be determined accurately from GW data alone, thus GW data can be used to select systems that will likely be eclipsing binaries for EM followup.
Key words: binaries: general / gravitational waves
© ESO, 2012
1. Introduction
The spacebased gravitational wave (GW) detector in consideration by ESA, eLISA, is expected to observe millions of compact Galactic binaries (AmaroSeoane et al. 2012; Nelemans 2009) with periods shorter than about a few hours, amongst other astrophysical sources, and resolve several thousand of these binaries (Nissanke et al. 2012). About 50 compact binary sources have been observed at optical, UV, and Xray wavelengths (e.g. Roelofs et al. 2010). The types of binaries known to us are interacting systems (AM CVn stars, ultracompact Xray binaries, and cataclysmic variables) and detached systems (double white dwarfs (WDs) and double neutron stars, Nelemans 2009, 2011). The AM CVn stars are binary systems where a WD accretes matter from a lowmass, heliumrich (hydrogendeficient) object (Solheim 2010). Their mass transfer is driven by GW radiation loss. The known ultracompact Xray binaries consist of neutron stars that are accretors whose donors are inferred to be either helium rich or carbon/oxygen rich (in’t Zand 2005). Double WDs are predicted to be the most common systems, which sometimes tend to be the outcome of many binary evolutionary paths (Webbink 1984). Of all these known systems, a handful lie in the eLISA band and will be individually detected. These are known as verification binaries since they are guaranteed sources for the detector. Parameter uncertainties in the verification binaries and Galactic binaries in general have been studied in the literature by extensively using Fisher information matrix (FIM) analyses (e.g. Cutler 1998; Stroeer & Vecchio 2006; Takahashi & Seto 2002). These studies have been done for various configurations of classic LISA (LISA Study Team 1998), which was designed to have a larger baseline of five million km with six laser links interchanging between three stations located at the vertices of a triangle that was to house two proof masses each (Vallisneri et al. 2008). Instead eLISA will have a baseline of one million km with four laser links interchanging between the proof masses. The parameter uncertainties depend intricately on the observation conditions and the geometry of the detector (Takahashi & Seto 2002) and the same is true for the (possible) correlations between the parameters of a Galactic binary. In this study, we wish to quantify whether any such correlations exist that could be useful in constraining the GW parameter estimates. Some of the GW parameters are the same as or related to electromagnetic (EM) parameters (e.g. inclination_{GW} = inclination_{EM}, f_{GW} = 2/P_{orb}, etc.). Thus, we can use an independent (EM) constraint of a (EM/GW) parameter that correlates strongly with another (GW) parameter to improve the accuracy of the latter. This work provides a first step towards developing strategic plans in observing those potentially useful EM parameters that can improve the GW parameter accuracy. In this paper, we present a FIM analysis to study whether correlations within GW parameters exist. For the useful correlations that we find, we predict quantitatively the improvement in the GW parameter when there is prior EM data in the correlated parameter. The paper is structured in the following way. In Sect. 2, we briefly summarise the signal models, both the instrumental and foreground noises, and our data analysis. We present our results and interpretations in Sect. 3. Finally, we discuss how the results differ from the old LISA detector in Sect. 4 and present our conclusions in Sect. 5.
Physical properties of verification binaries based on observations summarised in Nelemans (2011).
2. Signal modelling and data analysis
2.1. Gravitational wave signals from a Galactic binary
We consider three verification sources, AM CVn, SDSS J0651+2844 (hereafter J0651), and RX J0806.3+1527 (HM Cnc), whose physical parameters are summarised in Table 1.
AM CVn and HM Cnc are masstransferring systems known to astronomers as AM CVn binary systems, which were described in the introduction. J0651 (Brown et al. 2011) is an eclipsing detached WD binary system that was spectroscopically identified in the Sloan Digital Sky Survey (SDSS) catalogue. AM CVn and J0651 can be modelled as monochromatic GW sources, which means that they are described by seven parameters: dimensionless amplitude (), frequency (f), polarisation angle (ψ), initial GW phase (φ_{0}), inclination (cosι), ecliptic latitude (sinβ), and ecliptic longitude (λ). HM Cnc can be modelled as a mild chirper with an additional eighth parameter, the chirping frequency (ḟ). The two polarised gravitational waveforms used in the strain for slowly evolving binaries are given by (e.g. Królak et al. 2004) (1)(2)where (3)and (4)In these expressions, ℳ ≡ (m_{1}m_{2})^{3/5}/(m_{1} + m_{2})^{1/5} is the chirp mass and d is the distance to the source. The monochromatic waveforms are given by setting ḟ = 0 in the expressions above.
2.2. Detector response to the GW signals
If a GW signal is present, then the output of a detector will contain the strain, h(t), and noise, n(t). Thus, the detector registers (5)where θ is a vector characterising the seven (or eight) parameters of the binary. Since eLISA will have a motion around the Sun and a cartwheeling motion around its centre of mass, a monochromatic signal from a WD binary will be modulated in its amplitude, frequency, and phase in complicated ways. The resulting signal will be spread over a range of frequency bins of the detector (Cornish & Larson 2003). Galactic binaries typically radiate monochromatic signals at low frequencies^{1} and thus the response at the detector can be written as (Cornish & Larson 2003) (6)where (7)The functions F^{ + , × }(t) are the antenna beam patterns of the detector, and they depend on the source’s sky position (λ, β), its orientation (ψ), and the detector configuration. The phase of the signal is given by (8)where Φ_{D}(t),Φ_{P}(t) are frequency (Doppler) and phase modulations, respectively (Cornish & Larson 2003). Doppler modulation is given by (9)where f_{m} = 1/year is the modulation frequency. The phase modulation is given by (10)Different architectures of the triangular spacebased interferometer generate a number of independent data streams that provide different responses to the incoming GW signal (Vallisneri et al. 2008). For the most recent interferometer design in consideration, eLISA, the output is a single unequalarm Michelson data stream, X. This is a linear combination of phase shifts measured at the different spacecraft (by comparing the incoming light with a local reference source) shifted in time in such a way as to represent interference between two light beams travelling through the arms of the detector in opposite ways, which is a particular implementation of socalled time delay interferometry (TDI, Armstrong et al. 1999). Vallisneri (2005a) provide a detailed description of how TDI works and an explanation of why it produces an interferometry signal in which the phase shift induced by a passing GW is preserved, while the much larger shifts induced by instrumental noise are strongly suppressed. We made use of the existing numerical software Synthetic LISA (Vallisneri 2005b) to simulate accurate timedomain series of the instrumental noise and the GW signals in the form of these TDI observables. The data stream is a discrete series for a given observational time T_{obs}, where the samples are separated by Δt. The detector response to GWs and instrumental noises have been discussed extensively in the literature (e.g. Cornish & Rubbo 2003; Cutler 1998; Królak et al. 2004; Vallisneri 2005b), hence we only summarise the most essential expressions relevant to our data analysis. In Synthetic LISA, the strain at the detector, s(t), is modelled as the TDI X observable. This is the quantity we work with in our data analysis.
2.3. Noise
There are two types of noise to consider: instrumental noise and Galactic foreground noise due to unresolved compact binaries. For the particular geometry of eLISA used here, the instrumental noises (mostly from shot noise), acceleration noise, and other types of noise (e.g. optical bench noise), are characterised by their power spectral densities (PSDs^{2}) of 2.31 × 10^{38}f^{2}, 6 × 10^{48}f^{2}, and 2.76 × 10^{38}f^{2}, respectively, in units of Hz^{1}. The instrumental noise is modelled as a random, Gaussian process. We note that the sampling time, Δt, should be carefully chosen while simulating the instrumental noise in order to correctly interpret the TDI observables. The choice of Δt should correspond to a frequency that is several times higher than the highest frequency where the TDI responses have to be analysed (Vallisneri 2005b). This means that since AM CVn and J0651 have relatively low frequencies they can be analysed with samples of Δt = 64 s, whereas for HM Cnc, we need a lower sampling time of at least Δt = 16 s (see Fig. B.1 in Appendix B).
The foreground noise from the Galactic binaries is simulated using the Lisasolve (Vallisneri 2011) software where every binary (monochromatic and/or mild chirper^{3}) is modelled in the frequency domain (Cornish & Littenberg 2007). This differs from simulating signals using Synthetic LISA, where the signal is modelled accurately in the timedomain. Lisasolve instead makes use of the very slowly evolving nature of the binaries to approximately model the signals both directly and speedily in the frequency domain. We only use doubledetached WD binaries because they form the majority of the foreground noise (Nelemans et al. 2001; Petiteau 2012). We include more than 2.7 × 10^{7} detached Galactic double WDs from a simulation with the same assumptions about binary evolution and Galactic distribution as those in Nelemans et al. (2004), but with (about a factor of ten) higher intrinsic resolution. These assumptions have also been used for the simulations of the Galactic binaries on which the Mock LISA Data Challenge (MLDC) rounds are based (Littenberg 2011).
The unresolvable Galactic foreground noise is obtained by iteratively subtracting the resolvable sources (i.e. binaries expected to be detected) from the simulated population as follows:
 1.
the S/N (see Eq. (14)) is computed for each binary against the Gaussian instrumental noise and initial Galactic foreground noise for T_{obs} of two years. The initial foreground noise is calculated using the initial catalogue of 2.7 × 10^{7} detached Galactic double WDs;
 2.
all the sources with S/N > 5 are removed from the initial catalogue/dataset. The reduced dataset is used to simulate the reduced Galactic foreground;
 3.
using the reduced Galactic foreground and the same instrumental noise, the S/N for each of the binaries in the reduced catalogue is calculated. The process is iterated with step 2.
2.4. Data analysis
For GW sources with known waveforms, one can use matched filtering methods (Cutler & Flanagan 1994; Finn 1992) to extract the signal parameters and estimate their uncertainties. Consequently, when the noise is Gaussian, the parameter uncertainties are given by their joint Gaussian probability distribution function (Cutler 1998) (11)where Γ is known as the Fisher information matrix (FIM) given by (12)The inner product (...  ...) is a generalisation of the timedomain correlation product and is conventionally defined as (13)Equation (13)^{4} holds for quasimonochromatic binaries that have an almost constant noise PSD, S_{n}(f), in the frequency region where the binary radiates (Cutler 1998). The S/N of a source is defined as the inner product of a signal with itself (14)In the limit of signals with a high signaltonoise ratio (S/N ≫ 1), the inverse of the FIM gives the variancecovariance matrix . The diagonal elements give variances (or mean square errors) in each parameter, ⟨ (σ_{i})^{2} ⟩ , and the offdiagonal elements describe the covariances (or correlations) between them. For each of our verification binaries, we calculate this matrix to investigate the correlations between the binary parameters. The derivative in Eq. (12), ∂h/∂θ_{i}, is numerically calculated in the timedomain (15)where dθ_{i} should be chosen carefully. In general, the quantity dθ_{i} should be as small as the machine accuracy allows for, but not too large to suffer from the truncation error^{5}. Thus, for wellbehaved functions, d, where ϵ ~ 10^{16} is the machine accuracy and θ_{c} is some typical value of the corresponding parameter (Press et al. 2002). To find a good choice of dθ_{i}, we compute σ_{i} for a range of dθ_{i} with logarithmic intervals and select the value for each corresponding parameter around which the standard deviations become stable. By stable, we mean that increasing or decreasing dθ_{i} by an order of magnitude should lead to values of σ_{i} that vary by no greater (smaller) factor than 1.1 (0.9). An example of the stabilisation of the variancecovariance matrix is provided in Table A.1 of Appendix A for the case of AM CVn. We perform this stability check for all our verification binaries. For instance, in the case of AM CVn, d, dθ_{f} ~ 10^{11}, dθ_{ι,ψ} ~ 10^{8} etc. All our analysis is done for T_{obs} = 2 years.
3. Results
Here we list the variancecovariance matrix for AM CVn, HM Cnc, and J0651 and discuss the strongest correlations that we found. The offdiagonal elements are specified by the normalised correlations, c_{ij}, and the diagonal elements are quoted as the square root of the variances (or standard deviations), c_{ii}, i.e. (16)Thus, c_{ij} can have values in the range [−1, +1] where, c_{ij} = +1 means maximally correlated and c_{ij} = −1 means maximally anticorrelated. We consider highly correlated parameters θ_{i},θ_{j} to be those for which c_{ij} > 0.8. In the later subsections, we show that these correlations are affected when the inclination is varied.
The correlation matrices for the three verification binaries with observed parameters as shown in Table 1 are listed below. The measured GW parameter values are shown above each matrix. Since there are no EM measurements for φ_{0} and ψ, we set them to π and π/2, respectively, for all three binaries. The choice of the φ_{0} value does not influence the results shown below, whereas the choice for the ψ has an effect that will be addressed later. In the matrices, we take the medians of values obtained using 50 different instrumentalnoise realisations and the standard deviations about these medians (σ(S/N) and σ(σ_{i})). They are listed above and below each of the corresponding matrices, respectively. The values of S/N we find here are comparable to the ones found by other groups (Petiteau 2012). The normalised correlations are accurate within the quoted precision and hence their standard deviations from the instrumental noise are not listed.
J0651, S/N = 10.72 ± 0.23. Strong correlations (as defined above) are printed in bold face.
HM Cnc, with ḟ, S/N = 39.89 ± 0.85.
Fig. 1 Left: envelopes of the GW signal (TDI X observable) for binaries with parameters of AM CVn, but for various inclinations as a function of time (innermost to outermost envelopes correspond to ι = 90°,60°,45°,40°,20°,0° respectively). The modulation of the signal is due to the complicated annual motion of eLISA around the Sun. For low inclinations (ι = 0°−45°), the signals only differ in amplitude, while the higher inclinations also differ in structure. To highlight this, we plot in the right panel the normalised ratio of the upper envelopes compared to the ι = 0° envelope using the same line styles (but now the ι = 90° is the outermost line). High inclination systems clearly have a unique structure, while low inclination envelopes only differ in amplitude, illustrating the degeneracy between ι and . 

Open with DEXTER 
Fig. 2 Left: twodimensional error ellipses of and cosι extracted from the variancecovariance matrix, , for AM CVn. The two thin black lines represent 1σ and 2σ ellipses and the black dashed lines represent the 1σ GW errors in and cosι with the cross at the true value of the parameters. The unphysical values (i.e. negative numbers along the amplitude axis) are caused by the Gaussian tails of the parameter uncertainty about their true values. The right panel is a zoom of the area that is constrained by the 1σ EM error of 4° shown in dashdotted lines. The 1σ GW standard deviation in decreases from ~1.08 × 10^{22} to ~0.165 × 10^{22}, roughly a factor of 6.5. The corresponding error ellipses for the reduced PDF are shown as thick black ellipses for 1σ and 2σ. 

Open with DEXTER 
3.1. The correlation between and cos ι
The amplitude and inclination of AM CVn and HM Cnc (when modelled without the chirp, ḟ) are highly anticorrelated, with a value of . The reason for this is that binaries with inclinations in the range from 0° to ~45° (faceon or close to faceon) have signals with the same structure (see Fig. 1), so the signals in this range can be transformed into each other by simply changing either the amplitude or the inclination. Thus, if we change the amplitude of one of these signals within its GW uncertainty, we can obtain a very similar signal by changing its inclination within its GW uncertainty as well. However, for inclinations in the range ~45°−90° (either close to edgeon or edgeon), this scaling does not apply. This is because the signals in this range not only vary in terms of the overall scale but also in their structure, when their inclination is changed. This is especially evident for a binary with ι ~ 90° (lines of ι = 60°,90° in Fig. 1). Since the inclinations of both AM CVn and HM Cnc are within the range where similar signals can be obtained by varying either the inclination or the strain amplitude, the parameters cosι and are highly correlated. Furthermore since J0561 is an edgeon system with ι ~ 90°, both cosι and are distinguishable, as can be seen in Fig. 1. This is potentially the most useful correlation from the GW analysis where an EM measurement of the inclination of a binary can constrain the error in the strain amplitude . For example for AM CVn, an EM constraint in ι with a 1σ accuracy of 4° (Roelofs et al. 2006) improves the uncertainty in from 1.084 × 10^{22} to 1.65 × 10^{23} as shown in Fig. 2. The corresponding 1σ accuracy in is estimated by selecting the twodimensional probability distribution function (2D PDF) points that lie in the range 41° < ι < 45°. This improves amplitude accuracy by a factor of 6.5 compared to the use of GW data alone.
HM Cnc’s inclination measurement derived from the EM observation is quite uncertain compared to that for AM CVn. This is mainly due to the large uncertainty in the mass ratio and the assumption about the GW radiation that goes into estimating the inclination (Roelofs et al. 2010). Hence, if we take an EM constraint in ι with a 1σ uncertainty of 7°, the amplitude improves by a factor of 1.34 as shown in Fig. 3. However, if we were to obtain a better ι constraint of 4° for HM Cnc, would improve by a factor of three, which is still a smaller improvement than that of the AM CVn’s. This is due to the higher S/N of HM Cnc: if we place AM CVn at a closer distance such that it has a higher S/N of ~40, the 1σ accuracy of 4° constrains the amplitude more tightly by a factor of 2.95, which is very similar to the case of HM Cnc. By constraining the strain amplitude, we constrain the chirp mass, the distance, or a combination of the two (see Eq. (3)). The GW frequency, the third contributor to the amplitude, is already determined very precisely. Additionally, if ḟ is measurable for a chirping binary, then from Eqs. (3) and (4) one can constrain the error in the distance (Schutz 1996).
Fig. 3 Twodimensional error ellipses and cosι extracted from the variancecovariance matrix, , for HM Cnc. See Fig. 2 for details. The 1σ EM constraint of 7° in inclination reduces the corresponding 1σ uncertainty in from ~1.22 × 10^{23} to ~0.89 × 10^{23}, roughly a factor of 1.4. 

Open with DEXTER 
3.2. Influence of the inclination, ι
We have seen above that depends on the inclination of the binary. For the special case of exactly faceon binaries (ι = 0,π), the Fisher matrix is illdefined owing to the symmetry in the signal, namely, h(t,ι + dι) = h(t,ι − dι). The derivatives are zero, so the terms including ι in the Fisher matrix are also zero, resulting in an illdefined matrix that cannot be inverted. In general, it is hard to track the behaviour of the variancecovariance matrix from the FIM alone. Hence, to determine the influence of inclination, we fix all parameters except , ι, and λ, and calculate the matrices Γ and for AM CVn as a function of ι. This can be justified, as most correlations with and ι are not very strong; the largest correlation (of ~ 0.6) is for λ as shown in Fig. 4. This means that for all inclinations
thus we can focus on , cos ι, and λ independently and reduce the 7 × 7 matrix to a 3 × 3 matrix, which is easier to interpret^{6}. In Fig. 5, we show the S/N, (relative) parameter uncertainties and the normalised correlation for these three parameters as a function of inclination. It is obvious that for a source at a given distance, the S/N for an edgeon binary is lower than that for a faceon binary owing to the less favourable orientation of the latter. Thus, the S/N is a function of cosι and has a minimum at ι = 90° (cosι = 0). In general, we expect the parameter uncertainties to vary accordingly, i.e., a higher S/N should correspond to a higher precision and viceversa. The longitude follows this behaviour in the bottomleft panel. However, for and cosι this is not the case as the detector cannot accurately determine the inclination in the range of 0.7 ≲  cosι  ≤ 1, as shown in the topright and middleleft panels. For ι ≤ 10° and ι ≥ 170° in particular, the uncertainties in ι and become very large. However, as shown in the middleright panel, the correlations are maximal for these almost faceon binaries, thus by constraining one of the parameters from EM (or other) data, we constrain the range of values that the other parameter can take. This method cannot be applied to edgeon binaries as the correlation is very weak in these cases. If we take a value of dι for which the computation of the 3 × 3 variancecovariance matrix is numerically stable for all ι and use it to produce the full 7 × 7 matrix, the parameter uncertainties behave as shown in Fig. 6. The uncertainties σ_{φ0} and σ_{ψ} behave similarly to and σ_{ι} since φ_{0}, and ψ are also degenerate for nearly faceon systems. In contrast, the uncertainties σ_{f},σ_{λ}, and σ_{sinβ} behave as may be expected from the S/N since they are not strongly tied to any of the four parameters with strong correlations.
Fig. 4 Correlations of the remaining parameters with and ι for a binary with AM CVn’s parameter values as a function of cosι. Most correlations are weak for all inclinations except . 

Open with DEXTER 
3.3. The correlation between φ_{0} and ψ
The orientation parameters φ_{0} and ψ are also highly anticorrelated with c_{φ0ψ} ≈ − 0.99 for AM CVn and HM Cnc. This can also be explained by examining the geometry of the physical binary. The orbital phase φ is measured in the plane of the binary (perpendicular to the binary’s angular momentum vector, ), while ψ is measured in the plane perpendicular to the line of sight, . Both AM CVn and HM Cnc have an inclination in the range where the system appears to be faceon to the detector (see Sect. 3.1). Thus, φ_{0} and ψ are practically degenerate for these binaries. In contrast, for J0651 and are almost perpendicular, and the effects of φ_{0} and ψ can be more easily distinguished by eLISA. This correlation is not very noteworthy, because the parameters involved are of little astrophysical interest.
3.4. The correlation between f and φ_{0}
For J0651, there is a strong anticorrelation between f and φ_{0}, with a value of − 0.89. This is because if the frequency of a signal is changed by a small fraction, δf, the signal accumulates a phase shift towards the end of the wave. This signal can be deshifted in phase to obtain another signal with a similar match to the original, in turn making the two slightly dissimilar signals indistinguishable. Hence, the two parameters f and φ_{0} are degenerate. This degeneracy is not present for AM CVn and HM Cnc owing to the strong anticorrelation between φ_{0} and ψ and large errors in these two parameters.
Fig. 5 Signaltonoise ratio, , σ_{cosι}, σ_{λ}, , , and c_{cosιλ} as a function of cosι, keeping the other parameters fixed to the values of AM CVn. The vertical dashed line is the measured value of ι for this binary, 43°. The best known value for the amplitude is . As expected, the S/N is higher for faceon than edgeon orientations. The relative uncertainties in the amplitude and inclination are very large for inclinations up to 10°. The bottomright panel shows the normalised correlations as the dotted line and c_{cosιλ} as the solid line. 

Open with DEXTER 
3.5. The correlation between f and ḟ
If we model HM Cnc with an additional eighth parameter, ḟ, we find a strong anticorrelation between f and ḟ. This can be explained in a very similar way to the case of c_{fφ0} above. Consider a signal which is obtained by changing the chirp by δḟ, within its uncertainty. This will slightly change the frequency towards the end of the signal. By additionally changing f within its uncertainty, we can obtain a very similar signal to the one for which we changed only its chirp by δḟ. Thus, the two parameters are degenerate for a mildly chirping binary. This could be a very useful correlation for constraining the ḟ of the binary by using an EM constraint on f. However, the GW data analysis already provides a relative accuracy of ~10^{6} for the frequency, which is hard to improve with EM observations.
After providing a potential measurement of ḟ for the J0651 source (M. Kilic, priv. comm.), the correlation between f and ḟ could be used to test one of the predictions of General Relativity in Eq. (4), which holds for a binary that evolves only because of GW radiation loss. Given that the individual masses and their f are measured, Eq. (4) can be used to test the prediction of ḟ.
Fig. 6 Signaltonoise ratio and (relative) parameter uncertainties for AM CVn as a function of cosι. The general behaviour for , σ_{cosι} and σ_{λ} is consistent with Fig. 5 above. The uncertainties in σ_{sinβ},σ_{λ}, and σ_{f} behave as may be expected from the S/N, where parameters are more accurately determined for higher S/N. 

Open with DEXTER 
4. Discussion
The S/N of the data for the verification binaries listed in Table 1 are relatively low. In particular, J0651 has S/N ~ 10. This may not satisfy S/N ≫ 1, one of the criteria under which FIM studies are valid. To exclude the effects of S/N, we place the verification binaries at a closer distance so that they all have S/N of ~100. We find that all of the parameter uncertainties decrease as expected, except for the uncertainty in . However, the strong correlations we found for , c_{φ0ψ}, c_{fφ0} and c_{fḟ} listed above in the matrices remain the same. Furthermore, these correlations do not change when we consider the sixlink configuration of the classic LISA with a five million km baseline and slightly different instrumental noise. From the sixlink interferometer, three independent data streams can be synthesised, of which A, E, and T is one of the combinations. Naturally, this means that the S/N of the source will be higher, and that there may be additional information from these data streams. However, we found that when we consider these optimal TDIs, the strong correlations we found are unaffected. Additionally, the weak correlations between the rest of the parameters are also unaffected.
It is wellknown that in the lowfrequency regime there are many overlapping sources, which can significantly degrade the parameter uncertainty. Crowder & Cornish (2004) performed a thorough FIM analysis and show that the degradation of parameter uncertainties depends on the frequency, the number of binaries per frequency bin, and the observation time. However, the degradation is uniform for all the seven parameters, which means that each parameter is affected in the same way by the overlapping sources. Hence, the presence of overlapping sources affects the uncertainties in the parameters, but not the normalised correlation between the parameters. In addition the potential biases from imperfect subtraction of resolvable sources (Sect. 2.3) can affect the results, although this mostly affects the S/N degradation of the source and thus not the normalised correlations.
Our study should be used to implement EM priors in an MCMC data analysis based on a Bayesian framework. While the latter leads to more accurate results, the FIM tools provides useful predictions (Vallisneri 2008). In the case of Galactic binaries, it has been shown that the onedimensional marginalised posterior PDFs are wellmatched by the FIM predictions when the data analysis of a single source is considered (Cornish & Crowder 2005).
In a forthcoming paper, we will discuss the dependence of the aforementioned correlations and accuracies on the sky position and polarisation angle. Our test calculations show that the strong correlations , c_{φ0ψ}, c_{fφ0}, and c_{fḟ} found above do not vary as a function of these parameters.
Finally, we observe that the uncertainty in the inclination of a binary indicates whether the system is eclipsing. In Sect. 3.1, we found that for J0651, the uncertainties in ι,φ_{0}, and ψ are small compared to those of HM Cnc, even though its S/N is a factor of three lower than that of HM Cnc. An uncertainty of σ_{cosι} = 0.043 translates to σ_{ι} ~ 2.5°. Hence, for J0651like systems with S/N of ≳ 10, the small uncertainty inclination inferred from a GW analysis may be used to find candidate eclipsers.
5. Conclusion
We have performed Fishermatrix studies to investigate whether there are correlations between the parameters that characterise Galactic binaries with short periods, between about six minutes and a few hours, which lie in the eLISA frequency band, 10^{4} Hz − 1 Hz. We focused on three verification binaries, AM CVn, SDSS J0651+2844, and the mildly chirping HM Cnc. Our main findings are:

1.
There are strong correlations between the strain amplitude and inclination, and between the phase and polarisation angle for AM CVn and HM Cnc. For the latter, there is an additional strong anticorrelation between f and ḟ when it is modelled as a mildly chirping source.

2.
The remaining parameters are not or only weakly correlated for the binaries considered.

3.
The correlation between strain amplitude and inclination can be very useful in constraining the amplitude of the binary, which is a function of the two masses, the distance to the binary, and its frequency. Since the frequency is very accurately determined (of the order of 10^{9}Hz) by the GW data analysis, a combination of the masses and the distance can be constrained by using an EM constraint on the inclination.

4.
These correlations depend strongly on the inclination of the system: approximately faceon (ι ≲ 45°,ι ≳ 135°) binaries have strongly correlated parameters, whereas for edgeon binaries, the correlations become very weak.

5.
For binaries with very low inclinations (ι < 10° and ι > 170°), the uncertainties in amplitude, inclination, phase, and polarisation angle derived from GW data are very large.

6.
We have found that the influence of S/N on the correlations is not significant; placing the sources at shorter distances to give them higher S/N yields the same normalised correlations.

7.
The strong correlation between amplitude and inclination is the most useful correlation found, since inclinations are typically measured to a higher accuracy by EM observations than with GW measurements alone. Hence, this correlation can be used to constrain the amplitude by a factor of six for systems with orientations and S/N similar to those of AM CVn.

8.
We have found that some correlations also depend on sky position and polarisation angle, which we will discuss in a forthcoming paper.

9.
Even for signals with an S/N of ~10, we have been able to reliably determine the inclination when the system is edgeon. This will enable efficient searches for eclipsing systems with EM instruments.
Another reason to consider a smaller matrix is that it is difficult to stabilise for low inclinations ι ≤ 10°. We find that the optimal choice of dθ_{i} for which the parameter uncertainties σ_{i} are stable varies as a function of inclination. For low values of ι, the optimal dθ_{i} is orders of magnitude larger than for high values of ι. Note that the FIM (as opposed to ) is stable for the same choice of dθ_{i}. In addition, the FIM choice of dθ_{i} that is made for the inclination affects only matrix elements that contain terms related to ι. However, this influence spreads to all matrix elements of , because of the matrix inversion.
Acknowledgments
This work was supported by funding from FOM. We are very grateful to Michele Vallisneri for providing support with the Synthetic LISA and Lisasolve softwares. We thank the anonymous referee for the helpful suggestions to improve the paper.
References
 AmaroSeoane, P., Aoudia, S., Babak, S., et al. 2012, GW Notes, submitted [arXiv:1201.3621] [Google Scholar]
 Armstrong, J. W., Estabrook, F. B., & Tinto, M. 1999, ApJ, 527, 814 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, W. R., Kilic, M., Hermes, J. J., et al. 2011, ApJ, 737, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Cornish, N. J., & Crowder, J. 2005, Phys. Rev. D, 72, 043005 [NASA ADS] [CrossRef] [Google Scholar]
 Cornish, N. J., & Larson, S. L. 2003, Class. Quant. Grav., 20, 163 [NASA ADS] [CrossRef] [Google Scholar]
 Cornish, N. J., & Littenberg, T. B. 2007, Phys. Rev. D, 76, 083006 [NASA ADS] [CrossRef] [Google Scholar]
 Cornish, N. J., & Rubbo, L. J. 2003, Phys. Rev. D, 67, 022001 [NASA ADS] [CrossRef] [Google Scholar]
 Cutler, C. 1998, Phys. Rev. D, 57, 7089 [NASA ADS] [CrossRef] [Google Scholar]
 Cutler, C., & Flanagan, É. E. 1994, Phys. Rev. D, 49, 2658 [NASA ADS] [CrossRef] [Google Scholar]
 Finn, L. S. 1992, Phys. Rev. D, 46, 5236 [NASA ADS] [CrossRef] [Google Scholar]
 in’t Zand, J. J. M. 2005, A&A, 441, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Królak, A., Tinto, M., & Vallisneri, M. 2004, Phys. Rev. D, 70, 022003 [NASA ADS] [CrossRef] [Google Scholar]
 LISA Study Team. 1998, LISA: Laser Interferometer Space Antenna for the Detection and Observation of Gravitational Waves, PrePhase A Report, 2nd edn. (Garching, Germany: Max Planck Institut fur Quantenoptik) [Google Scholar]
 Littenberg, T. B. 2011, Phys. Rev. D, 84, 063009 [NASA ADS] [CrossRef] [Google Scholar]
 Nelemans, G. 2009, Class. Quant. Grav., 26, 094030 [NASA ADS] [CrossRef] [Google Scholar]
 Nelemans, G. 2011, Wiki Verification Sources [Google Scholar]
 Nelemans, G., Yungelson, L. R., & Portegies Zwart, S. F. 2001, A&A, 375, 890 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nelemans, G., Yungelson, L. R., & Portegies Zwart, S. F. 2004, MNRAS, 349, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Nissanke, S., Vallisneri, M., Nelemans, G., & Prince, T. A. 2012, ApJ, submitted [arXiv:1201.4613] [Google Scholar]
 Petiteau, A. 2012, Configurations of the eLISA (NGO) mission [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2002, Numerical recipes in C++: the art of scientific computing, eds. W. H. Press, S. A. Teukolsky, W. T. Vetterling, & B. P. Flannery [Google Scholar]
 Roelofs, G. H. A., Groot, P. J., Nelemans, G., Marsh, T. R., & Steeghs, D. 2006, MNRAS, 371, 1231 [NASA ADS] [CrossRef] [Google Scholar]
 Roelofs, G. H. A., Groot, P. J., Benedict, G. F., et al. 2007, ApJ, 666, 1174 [NASA ADS] [CrossRef] [Google Scholar]
 Roelofs, G. H. A., Rau, A., Marsh, T. R., et al. 2010, ApJ, 711, L138 [NASA ADS] [CrossRef] [Google Scholar]
 Schutz, B. F. 1996, Class. Quant. Grav., 13, A219 [NASA ADS] [CrossRef] [Google Scholar]
 Skillman, D. R., Patterson, J., Kemp, J., et al. 1999, PASP, 111, 1281 [NASA ADS] [CrossRef] [Google Scholar]
 Solheim, J.E. 2010, PASP, 122, 1133 [NASA ADS] [CrossRef] [Google Scholar]
 Stroeer, A., & Vecchio, A. 2006, Class. Quant. Grav., 23, 809 [NASA ADS] [CrossRef] [Google Scholar]
 Strohmayer, T. E. 2005, ApJ, 627, 920 [NASA ADS] [CrossRef] [Google Scholar]
 Takahashi, R., & Seto, N. 2002, ApJ, 575, 1030 [NASA ADS] [CrossRef] [Google Scholar]
 Timpano, S. E., Rubbo, L. J., & Cornish, N. J. 2006, Phys. Rev. D, 73, 122001 [NASA ADS] [CrossRef] [Google Scholar]
 Vallisneri, M. 2005a, Phys. Rev. D, 72, 042003 [NASA ADS] [CrossRef] [Google Scholar]
 Vallisneri, M. 2005b, Phys. Rev. D, 71, 022001 [NASA ADS] [CrossRef] [Google Scholar]
 Vallisneri, M. 2008, Phys. Rev. D, 77, 042001 [NASA ADS] [CrossRef] [Google Scholar]
 Vallisneri, M. 2011, Lisasolve fastbinary package [Google Scholar]
 Vallisneri, M., Crowder, J., & Tinto, M. 2008, Class. Quant. Grav., 25, 065005 [NASA ADS] [CrossRef] [Google Scholar]
 Webbink, R. F. 1984, ApJ, 277, 355 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Stabilisation of the variancecovariance matrix
Here we give an example of the Stabilisation of the variancecovariance matrix for AM CVn. The dθ_{i} are varied for each parameter over seven orders of magnitude as described in Sect. 2.4. In this example, d, dθ_{f}, dθ_{ι}, dθ_{ψ} are varied from (10^{4},...,10^{2}) × × θ_{i}, where, θ_{i} are the true values for AM CVn. The values of dθ_{φ}, dθ_{β}, dθ_{λ} are varied from 10^{6},...,10^{0}, more than for parameters above, because the numerical derivatives of φ, β, and λ are unstable for lower values. The first row of the uncertainties in the parameters correspond to the first set of dθ_{i}, et cetera. In principle, the uncertainties σ_{i} should level off around the third row. The uncertainties for all our verification binaries are calculated in this way.
Parameter uncertainties, σ_{i}, in a range of dθ_{i} for the case of AM CVn, with S/N ~ 11 for T_{obs} = 2 years.
Appendix B: Sampling
Fig. B.1 Average PSDs of instrumental noise for eLISA with varying sampling time, Δt. For lowerfrequency sources such as AM CVn and J0651, a sampling as fine as 64 s can be used. However, for HM Cnc, whose frequency is higher, a high Δt would overestimate its S/N, since the level of the noise is too low for the longest sampling times. 

Open with DEXTER 
In Fig. B.1, we show the effect of the sampling time Δt on the noise PSD of the instrumental noise, which is the average spectrum of the TDI X observable. At f < 5 × 10^{3} Hz, the
sampling time does not affect the level of the PSD. Above these frequencies, the PSD with a relatively long sampling time of 64 s underestimates its true level.
All Tables
Physical properties of verification binaries based on observations summarised in Nelemans (2011).
Parameter uncertainties, σ_{i}, in a range of dθ_{i} for the case of AM CVn, with S/N ~ 11 for T_{obs} = 2 years.
All Figures
Fig. 1 Left: envelopes of the GW signal (TDI X observable) for binaries with parameters of AM CVn, but for various inclinations as a function of time (innermost to outermost envelopes correspond to ι = 90°,60°,45°,40°,20°,0° respectively). The modulation of the signal is due to the complicated annual motion of eLISA around the Sun. For low inclinations (ι = 0°−45°), the signals only differ in amplitude, while the higher inclinations also differ in structure. To highlight this, we plot in the right panel the normalised ratio of the upper envelopes compared to the ι = 0° envelope using the same line styles (but now the ι = 90° is the outermost line). High inclination systems clearly have a unique structure, while low inclination envelopes only differ in amplitude, illustrating the degeneracy between ι and . 

Open with DEXTER  
In the text 
Fig. 2 Left: twodimensional error ellipses of and cosι extracted from the variancecovariance matrix, , for AM CVn. The two thin black lines represent 1σ and 2σ ellipses and the black dashed lines represent the 1σ GW errors in and cosι with the cross at the true value of the parameters. The unphysical values (i.e. negative numbers along the amplitude axis) are caused by the Gaussian tails of the parameter uncertainty about their true values. The right panel is a zoom of the area that is constrained by the 1σ EM error of 4° shown in dashdotted lines. The 1σ GW standard deviation in decreases from ~1.08 × 10^{22} to ~0.165 × 10^{22}, roughly a factor of 6.5. The corresponding error ellipses for the reduced PDF are shown as thick black ellipses for 1σ and 2σ. 

Open with DEXTER  
In the text 
Fig. 3 Twodimensional error ellipses and cosι extracted from the variancecovariance matrix, , for HM Cnc. See Fig. 2 for details. The 1σ EM constraint of 7° in inclination reduces the corresponding 1σ uncertainty in from ~1.22 × 10^{23} to ~0.89 × 10^{23}, roughly a factor of 1.4. 

Open with DEXTER  
In the text 
Fig. 4 Correlations of the remaining parameters with and ι for a binary with AM CVn’s parameter values as a function of cosι. Most correlations are weak for all inclinations except . 

Open with DEXTER  
In the text 
Fig. 5 Signaltonoise ratio, , σ_{cosι}, σ_{λ}, , , and c_{cosιλ} as a function of cosι, keeping the other parameters fixed to the values of AM CVn. The vertical dashed line is the measured value of ι for this binary, 43°. The best known value for the amplitude is . As expected, the S/N is higher for faceon than edgeon orientations. The relative uncertainties in the amplitude and inclination are very large for inclinations up to 10°. The bottomright panel shows the normalised correlations as the dotted line and c_{cosιλ} as the solid line. 

Open with DEXTER  
In the text 
Fig. 6 Signaltonoise ratio and (relative) parameter uncertainties for AM CVn as a function of cosι. The general behaviour for , σ_{cosι} and σ_{λ} is consistent with Fig. 5 above. The uncertainties in σ_{sinβ},σ_{λ}, and σ_{f} behave as may be expected from the S/N, where parameters are more accurately determined for higher S/N. 

Open with DEXTER  
In the text 
Fig. B.1 Average PSDs of instrumental noise for eLISA with varying sampling time, Δt. For lowerfrequency sources such as AM CVn and J0651, a sampling as fine as 64 s can be used. However, for HM Cnc, whose frequency is higher, a high Δt would overestimate its S/N, since the level of the noise is too low for the longest sampling times. 

Open with DEXTER  
In the text 