Using electromagnetic observations to aid gravitationalwave parameter estimation of compact binaries observed with LISA
II. The effect of knowing the sky position
^{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: 17 January 2013
Accepted: 18 March 2013
In this followup paper, we continue our study of the effect of using knowledge from electromagnetic 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 the first paper, we have shown that the strong correlation between amplitude and inclination can be used for mildly inclined binaries to improve the uncertainty in amplitude, and that this correlation depends on the inclination of the system. In this paper we investigate the overall effect of the other orientation parameters, namely the sky position and the polarisation angle. We find that after the inclination, the ecliptic latitude of the source has the strongest effect in determining the GW parameter uncertainties. We ascertain that the strong correlation we found previously only depends on the inclination of the source and not on the other orientation parameters. We find that knowing the sky position of the source from electromagnetic data can reduce the GW parameter uncertainty up to a factor of ~2, depending on the inclination and the ecliptic latitude of the system. Knowing the sky position and inclination can reduce the uncertainty in amplitude by a factor larger than 40. We also find that unphysical errors in the inclinations, which we found when using the Fisher matrix, can affect the corresponding uncertainties in the amplitudes, which need to be corrected.
Key words: gravitational waves / methods: data analysis / binaries: general
© ESO, 2013
1. Introduction
Even though population synthesis studies have predicted ~10^{8} galactic whitedwarf binaries (e.g. Nelemans et al. 2001), there are only about ~50 compact binary sources that have been studied at optical, UV, and Xray wavelengths (e.g. Roelofs et al. 2010). A spacebased gravitational wave (GW) detector, such as the proposed ESA mission eLISA, is expected to observe millions of these predicted compact Galactic binaries with periods shorter than about a few hours (Nelemans 2009; AmaroSeoane et al. 2013), amongst other astrophysical sources, and resolve several thousands of these binaries (Nissanke et al. 2012). In a recent paper, we investigated the effect of using knowledge from electromagnetic (EM) observations in the gravitational wave data analysis of Galactic binaries by exploring the correlations that might exist between these parameters. We found that for binaries with relatively low inclinations (ι ≲ 45°), an EM constraint on inclination can improve the estimate of the GW amplitude by a significant factor (Shah et al. 2012, Paper I, hereafter). The improvement depends on the signaltonoise ratio (S/N) of the source and on the EM constraint on inclination itself. We also found that this correlation is a strong function of the inclination of the system, but that it is independent of other angular parameters that describe the binary, in particular, sky position, and polarisation angle. From the initial results we also found that several other correlations that were rather small as a function of inclination (e.g. Fig. 4 in Paper I) change drastically when computed as a function of sky position. We address this problem here and also investigate the general dependence of parameter uncertainties on the orientations of the binary. Additionally, we discuss the influence of knowing the sky position from EM data on the parameter uncertainties of the binary sources.
The angular resolution of cLISA^{1} for monochromatic binaries has been studied previously (Peterseim et al. 1997; Cutler 1998; Rogan & Bose 2006; Błaut 2011). Peterseim et al. (1997) have shown the uncertainty in ecliptic latitude and ecliptic longitude as a function of both the parameters for a specific case of a crosspolarised source radiating at f = 3 mHz, which would be observable by cLISA. For a source with S/N = 115, depending on the values of the ecliptic latitude and longitude, they found the uncertainties in these angles to range from 1 to 8 milliradians. However, they assumed that except for the sky position, all other parameters are known a priori. Their study was also limited by a very simplified detector response and was tailored for a system at an inclination which is favourable in extracting its parameters. Relaxing the assumption of prior knowledge of other parameters, Cutler (1998) has shown that for a binary for a few sets of sky positions, the solid angle uncertainty (σ_{Ω}) ranges from 3 to 300 square degrees, depending on the frequency of the source. In subsequent work by Rogan & Bose (2006), sky position uncertainties were estimated for signals modelled for cLISA. They included chirp as an additional parameter, and these signals were not limited to a certain polarisation. Recently, Błaut (2011) has shown sky position uncertainties for various combinations of data channels from cLISA and included detailed detector response for the two distinct cases of low and highfrequency sources. For a 10^{3} Hz source with an edgeon orientation and an S/N of 10, an observation of one year is predicted to give a solid angle resolution of ~ 0.1 steradians which is a result very similar to the study by Rogan & Bose (2006). However, none of these previous authors provided an overall effect of varying multiple angular parameters, which can have a significant effect on parameter uncertainties. In this paper we aim to generalise the uncertainties in all parameters for all possible orientations of the binary and calculate the improvement in these uncertainties that is possible thanks to prior knowledge of the sky position and/or inclination of the system. We show this for eLISA.
A monochromatic binary expected to be observed by eLISA is completely described by seven parameters: dimensionless amplitude (), frequency (f), polarisation angle (ψ), initial GW phase (φ_{0}), inclination (cosι), ecliptic latitude (sinβ), and ecliptic longitude (λ). The inclination of the binary is defined as cosι = n·L, where n is the lineofsight vector from the observer to the source and can be described by the sky position parameters β and λ. L is the orientation of the orbital plane of the binary (or angular momentum vector), which can be specified by θ_{L} and φ_{L}. One of the main findings of Paper I is that the parameter uncertainties of a monochromatic binary and the three strongest correlations between these parameters (the correlation between amplitude and inclination , the correlation between phase and polarisation c_{φψ}, and the correlation between phase and frequency c_{fφ}) are strong functions of the source’s inclination ι. Furthermore, we found that the parameter uncertainties of such a binary also change as a function of other orientation parameters. These are parameters that describe the relative geometry of the binary in relation to the detector: ψ, λ, and β. It is natural to then ask the questions, how the possible strong correlations between the parameters due to these four orientation parameters influence these uncertainties, and to which extent knowledge from the EM data can improve the remaining parameter determinations. Of the remaining three parameters, (see Eq. (3) in Paper I) is intrinsic to the source and does not depend on its relative orientation to the detector. Since it only scales the S/N of the system, varying the amplitude thus will affect parameter uncertainties as expected, where larger will give smaller uncertainties and vice versa. The initial phase φ_{0}, which does not evolve for a monochromatic binary in the eLISA band, has no effect on the signal. On the other hand, the effect of frequency on the parameter accuracies can be significant; the strength of the modulations in the signal (which is caused by the motion of the detector around the Sun and relative to the source) affects lowf sources differently than highf sources. This has been addressed in detail in previous studies, e.g. Błaut (2011). Hence, here we focus on the effects on the parameter uncertainties of a Galactic binary due to sky position and the orientation parameters.
In this study, we consider lowfrequency binaries. We start by briefly summarising the analysis methods and then describe our Monte Carlo (MC) simulations in Sect. 2. In Sect. 3, we summarise the important results and our interpretation. In Sect. 4 we show that uncertainties can be severely overestimated for parameters with clear physical bounds, such as inclination, and thus also for parameters that are tightly correlated with them. We also briefly discuss the S/N limit for which the results are valid in Sect. 4, and present our conclusions in Sect. 5.
2. Monte Carlo simulations
Our signal and noise modelling and our data analysis method have been described in Paper I and we briefly summarise them here. The signals from the binaries that are expected to be observed by eLISA are modelled as monochromatic sources characterised by seven parameters, as mentioned above. Most of these sources radiate at low frequencies (milliHertz), including some of the verification binaries like AM CVn, and at these low f, the detector response can be derived in the longwavelength approximation. The noise is composed of instrumental noise described as a Gaussian random process, and a foreground noise mostly arising from the millions of doubledetached Galactic whitedwarf binaries. After removing the binaries, which are predicted to be resolved due to their strong signals, the foreground noise is also described by a Gaussian distribution. Given that a signal has a known waveform and the noise is Gaussian, we can use the Fisher matrix to estimate the uncertainties in the parameters that describe the waveform. The inverse of the Fisher information matrix (FIM), known as the variancecovariance matrix (VCM), gives the uncertainties and the correlations between the parameters. Given a number of parameters n, one can calculate the n × n VCM. For details we refer to Paper I.
Uncertainties in parameters are proportional to the inverse of the S/N, but they are also influenced by the correlations between parameters. Generally, a system with higher S/N has smaller parameter uncertainties and vice versa. However, if there is a strong correlation between any two parameters, this (in addition to the S/N) will influence their uncertainties. As mentioned in the previous section, we find that the overall effect of the sky position and the orientation parameters of a binary has not been investigated and this is one of the goals of this study. To do this, we fixed the values of , f, and φ_{0} to those of AM CVn (see Paper I) and varied its ψ, ι, β, and λ within their allowed ranges using a Monte Carlo method. We carried out a number of of MC simulations. Each simulation corresponds to varying a different (sub)set of these four parameters simultaneously. The motivation is to understand how each of the four parameters above influence the GW data analysis individually and in combination with each other. We explain these motivations briefly below, which will become clearer in Sect. 4 where the choice of each MC simulation is explained in more detail.
First of all, we set the distance of the AM CVn system close enough so that the S/Ns are sufficient to believe the Fisher matrix results (see Sect. 4). The first set of simulations we performed for an AM CVn of a fixed , f, and φ_{0}, while its position and orientation parameters λ, β, ψ, and ι were all varied simultaneously 4000 times. Thus, from the sample of these 4000 systems with various orientation parameters we first ascertained that the dominant effect on the parameter uncertainties is indeed the inclination of the system (result from Paper I). Next, we can quantify how the other three angles affect the parameter uncertainties globally. We performed another set of calculations where we fixed inclination (to a few choices explained later), frequency, and the phase (like above), and we varied the amplitude (by changing the distance, d), the sky position, and the polarisation angle. For each choice of the inclination, MC simulations were performed for 500 systems to quantify the effect of varying the S/N (by changing the distance) and the ecliptic latitude, longitude, and the polarisation angle on the parameter uncertainties. Finally, to quantify the possible improvement in the parameter uncertainties gained by prior (EM) knowledge of the sky location, we carried out 500 MC simulations for an AM CVn at a fixed inclination and varied the rest of the orientation parameters and the distance. However, instead of a full 7 × 7 VCM, we calculated a 5 × 5 VCM, where for each system the randomly chosen sky locations were fixed to that value. Similarly, we repeated this exercise to quantify the improvement in the parameter accuracies gained by EM knowledge of inclination in addition to sky position, and thus we calculated a 4 × 4 VCM for each system. For clarity, we briefly summarise these four sets of MC simulations below:

1.
MC1: fix (i.e. d, m_{1}, and m_{2} where, m_{i} is the mass), f, and φ_{0}. Randomly pick the orientation of the binary, L(φ_{L},θ_{L}), λ, and β over their physical ranges. The parameters ψ, and ι are then given by the following equations (Cornish & Rubbo 2003):

2.
MC2: for each ι = 30°,60°,90°, randomly pick λ, β, ψ, and d.

3.
MC3: like MC2, but fix λ, and β assuming an a priori known sky position [5 × 5 FIM].

4.
MC4: like MC2, but fix λ, β, and ι assuming an a priori known sky position and inclination [4 × 4 FIM].
Each of the four parameters are randomly picked from their uniform distribution in the following range: λ ∈ [0,2π], sinβ ∈ [ − 1,1], φ_{L} ∈ [0,2π], θ_{L} ∈ [0,π], and d ∈ [0.01,0.3] kpc. The corresponding ranges in ι and ψ are [0,π]. We also define the following inclinations for presenting the results below (unless specified): faceon:= 0° ≤ ι ≤ 45° or 135° ≤ ι ≤ 180°, mildly faceon:= 45° ≤ ι ≤ 60° or 120° ≤ ι ≤ 135°, mildly edgeon:= 60° ≤ ι ≤ 80° or 100° ≤ ι ≤ 120°, and edgeon:= 80° ≤ ι ≤ 100°.
3. Results
3.1. Global dependence of the uncertainties: dominance of inclination (MC1)
We have seen in Paper I that the parameter uncertainties are a strong function of the inclination of the system, which is due to the strong dependence of the signal strength, i.e. S/N, on the inclination and to the strong correlations. In this section, we first discuss the global behaviour of the parameter uncertainties as a result of varying the sky position, inclination, and the polarisation angle at the same time. We verify the result from Paper I that inclination has a dominant effect on the parameter uncertainties. We show that after the inclination, the ecliptic latitude, β, of the system has the strongest influence on the parameter uncertainties.
Using the results from MC1, we show in Fig. 1 the S/N of binary systems as a function of the four angular parameters over which they are randomly picked. As expected, there is a strong dependence of S/N of a system on inclination, as shown in the topleft panel. Due to this dependence on the inclination, in the other panels we find that the faceon systems form a band at the top of these plots, followed by bands of mildly faceon, mildly edgeon, and edgeon systems, in that order. In the topright panel, there is a mild but noticeable increase in S/N, for the binaries at the ecliptic plane, sinβ = 0, whereas the ecliptic longitude and the polarisation angle have no significant effect on the binaries’ S/N, as shown in the bottom panels. Figure 2 is a plot of parameter uncertainties as a function of sinβ taken from MC1. In the top panels, uncertainties in , ι, φ, and ψ are smallest for edgeon systems (forming a band at the bottom) and progressively increase for the faceon systems (the top of distribution). This is discussed in Paper I and results from the correlations between these parameters. The lower inclination systems ι ≲ 45° or ι ≳ 135° have very similar signal shapes in those ranges, whereas systems with close to edgeon orientations are distinguishable by both the shape and structure for small differences in inclinations. For faceon systems a small change in inclination is indistinguishable from a small change in its intrinsic amplitude, whereas for an edgeon system, a small change in inclination produces a noticeably different signal. Thus even for their lower S/Ns, the inclinations of the edgeon systems are better determined than those of the faceon systems. A similar argument based on the strong correlation between φ and ψ for faceon systems can explain the huge uncertainties of these parameters despite their strong signals (see Paper I). Since there are no such correlations between f, λ, and β (or the solid angle Ω), their uncertainties should be a consequence of the S/N, i.e. from the top panels of Fig. 1: a higher S/N implies lower values for the uncertainties and vice versa, which indeed holds for σ_{f} and σ_{λ}. From Figs. 1 and 2 it is evident that inclination has the strongest effect on the S/N and so to the leading order this parameter affects the parameter uncertainties, which confirms our results in Paper I.
Fig. 1 Twodimensional distributions of S/N as a function of cosι, sinβ, λ, and ψ. Faceon, mildly faceon, mildly edgeon, and edgeon systems (whose exact bounds are specified in Sect. 2) can be read off from the topleft panel. In the other panels, the trend from faceon to edgeon systems follows from the top of the distribution to the bottom. 

Open with DEXTER 
Fig. 2 Twodimensional distributions of the uncertainties in the seven GW parameters of the galactic binary systems as a function of sinβ from MC1 as described in the text. Units of the angles are in radians. 

Open with DEXTER 
After inclination, the strongest effect on the S/N of the systems is produced by the ecliptic latitude. Thus, we expect the uncertainties in the parameters also to be affected by β (see topright panel in Fig. 1); sources at the ecliptic poles should have larger uncertainties than those at the ecliptic plane. In Fig. 2 we can see that for a fixed distance to the source d, most uncertainties (, σ_{ι}, σ_{φ}, σ_{ψ}, σ_{f}, and σ_{λ}) are indeed smaller towards the plane than close to the poles. This is probably largely due to the S/N as a function of sinβ. The unusual behaviour of σ_{β}, which increases at the plane compared to its value at the poles (sinβ = ± 1), has been pointed out in Cornish & Larson (2003) and is caused by the Doppler modulation Φ_{D} (Eq. (9) in Paper I). This behaviour is surprising since Φ_{D} has the largest magnitude on the plane (since it is proportional to cosβ), and so we naively expect the sky resolution to be better on the ecliptic plane. However, the resolution in β is determined by the rate of change of the Doppler modulation with respect to β (i.e. proportional to sinβ), which is the highest at the pole and lowest at the plane. This is also evident in the solid angle uncertainty σ_{Ω} in the bottomright panel in Fig. 2, where , being the unnormalised correlation between sinβ and λ (Cutler 1998). We also observe that σ_{λ} has a distinct U shape, becoming very severe at the poles. This is merely because λ becomes progressively less welldefined as a source lies closer to the pole. Note that at the pole λ is not defined, which means that it can take any value.
Fig. 3 Uncertainties in all seven parameters of the systems as a function of S/N from MC2. Open circles, squares, and filled triangles represent binary systems at fixed inclinations of 90°,60°, and 30° respectively. The size of the symbol represents the magnitude of  β  for each simulated binary. In the bottom panels the different groups of symbols overlap. 

Open with DEXTER 
Fig. 4 Twodimensional distributions of the normalised correlations as a function of sinβ from MC2. We show here only the strong (normalised) correlations, meaning  c_{ij}  > 0.5. 

Open with DEXTER 
3.2. Dependence of the uncertainties on β: S/N or correlations (MC2)?
The general behaviour of the parameter uncertainties (except σ_{β}, hence σ_{Ω}) in Fig. 2 shows a minimum at the ecliptic plane that could be attributed to the S/N, since it is the highest for sources in the plane. To investigate this change in the uncertainties as a function of sinβ, we consider the results from MC2 and the correlations between parameters from MC1. In MC2 we varied S/N (by varying the distance to the binary, d), λ, β, and ψ for a few fixed inclinations. We fixed the ι for each of the MC2 simulations because we established its dominant effect through the S/N in the previous section, therefore the goal now is to determine which contributes most to the behaviour observed in Fig. 2: the S/N (topright panel of Fig. 1), or the correlations as a function of sinβ. If the larger uncertainties at the pole compared to the plane are solely due to the corresponding S/Ns, then we should see from the MC2 simulations that uncertainties anticorrelate tightly with the S/N of the sources. In Fig. 3 we show three distinct curves in S/N vs. parameter uncertainty that roughly show monotonic functions where all the uncertainties decrease with increasing S/N. However, all uncertainties except σ_{f} have a significant horizontal scatter. In the figure, the size of the symbol represents the magnitude of the ecliptic latitude. Higher values of  β , which represent regions around both ecliptic poles, concentrate towards the righthand side of the plots (with the exception of σ_{β}), and the opposite is true for the lower values of  β , which roughly represent the ecliptic plane. In the top panels, from left to right, the curves correspond to ι = 90°,60°, and 30° respectively, and these overlap in the bottom panels. As explained above, , σ_{ι}, σ_{ψ}, and σ_{φ} are larger for faceon systems (ι = 0,30°) than for edgeon systems (ι = 90°), owing to the correlations between these parameters even though faceon systems have higher S/Ns. Since these correlations are very weak for λ, β, and ψ for all inclinations (see Paper I), their uncertainties overlap in the bottom panels of Fig. 3. This spread across the horizontal axis must be due to changes in their correlations as a function of the ecliptic latitude. This is indeed shown in Fig. 4, which is taken from MC1. In the figure we show a 2d histogram to display the global dependency of the correlations on the ecliptic latitude. We notice that the correlations involving λ, ι, ψ, φ, and (, c_{φι}, c_{φλ}, c_{ιλ}, and c_{ψλ}) are generally weaker at the ecliptic plane. However, at the poles these correlations can take any value from 0 to ± 1 and concentrate at ± 1. This is the case for all inclinations. Accordingly, close to the ecliptic poles, there are sets of λ and ψ for which the correlations can become very strong, which explains the horizontal spread of uncertainties in Fig. 2. Note that in Fig. 4, we have shown (only) the normalised correlations c_{ij}, whose absolute values can exceed 0.5. A detailed explanation of these correlations can be found in the Appendix.
3.3. Gain from a priori knowledge of sky position and inclination (MC and MC4)
From Figs. 2 and 3 we observe that uncertainties in sky location are rather large (which has also been pointed out in previous studies). For the S/Ns in this study (~ 8 − 100), the uncertainties in σ_{Ω} range from ~ 0.2 to 6 square degree, as shown in Fig. 1. If EM data can provide a more accurate sky position and possibly an inclination for AM CVn system, we can use these independent (EM) measurements to improve accuracies in other parameters. To see by which factors the uncertainties improve, we fixed the sky position of each source and calculated the corresponding 5 × 5 VCM. The ratios of uncertainties from using a 7 × 7 VCM to using a 5 × 5 VCM, i.e. the improvement when knowing the sky position (ecliptic latitude and longitude) exactly, are shown in Fig. 5, calculated from MC2 and MC3. In the figure, the size of the symbol represents the S/N of that system. The figure shows that the improvement factor depends on the inclination of the source: for edgeon systems (in open circles with ι = 90°) σ_{ψ} can improve by a factor of ~ 3, but for systems closer to a faceon configuration in squares and filled triangles the improvement can only be up to a factor of ~ 2.2. One useful result from this figure is that and σ_{ι} can improve by a factor of ~ 2 for the systems close to the ecliptic poles. This gives a motivation for independent EM observations of the sky location. However, for edgeon systems there is no improvement in even though σ_{ι} improves^{2}. This is simply due to the lack of correlation between these two parameters for ι = 90°.
Fig. 5 Uncertainty ratios of the parameters between analyses using a 7 × 7 VCM and a 5 × 5 VCM, where the sky position is fixed in the latter, as a function of sinβ (from the MC2 and MC3 simulations). Filled triangles, squares and open circles represent systems at fixed inclinations of 30°,60° and 90° respectively. The size of the symbol for each binary represents the S/N of the system. 

Open with DEXTER 
Fig. 6 Same as in Fig. 5 but now we assume that in addition to sky position we also know the inclination of the binary. The uncertainty ratios of the parameters are between 7 × 7 VCM (MC2) and 4 × 4 VCM (MC4). The improvement factors in amplitude include corrections due to overestimated uncertainties in the amplitudes of faceon systems (see Sect. 4.1). Note that the last panel is plotted in log scale in order to highlight the relatively large improvements for the systems with smaller inclinations. 

Open with DEXTER 
In Paper I we found that for the HM Cnc system, the improvement of by using the EM constraint in ι was smaller than for AM CVn, although the two have similar inclinations (thus the same (normalised) correlation between and ι). This was due to the relatively high S/N of HM Cnc compared to the latter and thus the individual improvement factors shown in Figs. 5 and 6 will also depend on the individual S/N of the systems (see Paper I).
If in addition to sky position the inclination can be constrained, as discussed previously in Paper I, the improvement factor in amplitude can be much larger because of the decoupling of the correlation between amplitude and inclination. We show the ratio of parameter uncertainties from calculating the 7 × 7 VCM to the uncertainties from a 4 × 4 VCM (where λ,β and ι are fixed) in Fig. 6. The improvements in σ_{ψ} and σ_{φ} are roughly similar to the case in Fig. 5. A remarkable difference is seen in the improvement of for binaries with ι = 30° represented in filled triangles plotted in the bottomright panel of Fig. 6. This is expected because by fixing ι, the strong correlation between amplitude and inclination is decoupled. Hence for the higher S/Ns of close to faceon systems, the uncertainties in is significantly reduce. Note that here we assumed that the inclination is known exactly.
Fig. 7 Left: 1σ uncertainty in the inclinations vs. source inclinations for two values of S/N: plusses (outer curves) correspond to S/N = 10, stars (inner curves) correspond to S/N = 53. The curves above the dotted line correspond to ι + σ_{ι}, and those below to ι − σ_{ι}. The thick horizontal dashed and dasheddotted lines are the physical limits of ι: below the line of ι ± σ_{ι} = 0° and above the line of ι ± σ_{ι} = 90°, σ_{ι} can be truncated, since that value of ι is bounded between [0°,90°]. Right: “unphysical ι” versus S/N for the systems from MC1. The thick dasheddotted line calculated from the condition ι − σ_{ι} = 0°, and the thick dashed line from the condition ι + σ_{ι} = 90° demonstrated in the left panel. At an S/N of 10, the in open circles show a lower limit and the filled circles show an upper limit at the intersections with the two curves from the conditions above. 

Open with DEXTER 
4. Discussion
4.1. Limiting inclinations to their physical range
Some of the uncertainties in ι in Fig. 2 from MC1 are noteworthy in that they exceed their physical bounds, for example σ_{ι} = 7 radians. The cause for this is that the linearised signal approximation that is internally employed in the Fisher matrix studies is not sensitive to the bounded^{3} parameters that describe the model signal, as has been pointed out in Vallisneri (2008). Thus, when an uncertainty in a (cyclicly) bounded parameter like the inclination exceeds its physically allowed range, it means the quantity cannot be determined from GW data analysis alone. For the parameter ι, we can further ask what the unphysical uncertainty in ι implies for the uncertainty in for the mildly inclined systems that has a strong correlation between the two parameters. Since the Fisher matrix estimates σ_{ι} to unphysical values while in reality it is constrained within [0°,90°], does this mean that the uncertainty in the amplitude (due to the correlation with ι) is also overestimated? Or should we not trust the uncertainty in all those parameters where the σ_{ι} becomes unphysical? Consider the error ellipse for a given signal h, with true signal values ι_{0} and , where the two parameters are strongly correlated and where σ_{ι} ≫ π. The corresponding accuracy is then also highly inflated. Now consider a second signal , inside the error ellipse, but far from the true values, and such that ι_{1} = ι_{0} + 2nπ with n > 1. From its cyclic behaviour, should hold. However, this last point would typically lie far outside the error ellipse in the plane, which would indicate that this is actually a signal that is very dissimilar to the true signal. Hence, the fact that the inclination is physically limited to the range [0,π] and that this is not taken into account in the FIM method means that we should not only limit the uncertainty in inclination to its physically allowed range, but that we should limit the uncertainty in amplitude accordingly. Specifically, the corrected amplitude uncertainty, , is given by, (3)where is the corrected inclination uncertainty.
The value of inclination where its uncertainty becomes unphysical can occur at two places, i.e. where ι + σ_{ι} > π/2^{4}, and/or ι − σ_{ι} < 0, as shown in the left panel of Fig. 7. In the panel, the 1σ uncertainty ranges of the inclination are plotted against the sources’ inclination for two cases of constant S/N of 10 and 53. The thick horizontal dashed and dasheddotted lines indicate the conditions where σ_{ι} becomes unphysical. At a fixed S/N of 10 (in plusses), these values are at ι ~ 48° and ι ~ 47°, respectively, as indicated by the open circles, whereas for a higher S/N of 53 (in stars) the above conditions for σ_{ι} occur at ι ~ 27° and ι ~ 32° indicated by filled circles. We calculated the inclinations at which these intersections occur for a range of S/Ns, 10 − 110 and show the results in the right panel of Fig. 7. These curves can be used to apply corrections to the limits of σ_{ι} in Fig. 2 from MC1, which in turn can be used to limit from Eq. (3). In Fig. 8 we compare the uncertainties in corrected in this way with those from Fig. 2. The corrected values for differ for the high S/N systems, as can been seen by comparing the two righthand panels in Fig. 8. In particular, the tail at the large present in the top panel is absent from the bottom panel.
There is a caveat to applying these corrections. In the left panel of Fig. 7 we can understand the cutoff at 0° at the bottom of the figure. However, at the top, the cutoff of 90° is not so clear. A system with a lower inclination (e.g. 20°) cannot have an error larger than ~ 50° because a system with an inclination in the range ~ 80° − 90° is distinguishable from a system with an inclination of ~ 0 − 60°, as has been pointed out in Paper I. Thus the cutoff of ι + σ_{ι} = π/2 for lower inclinations is a conservative estimate. The estimated posterior probability distribution functions (PDFs) obtained using the FIM do not predict the real PDFs for lower inclinations and lower S/Ns. Thus a full Bayesian analysis that takes into account the limits of the cyclic parameters as priors is needed to fully determine the PDF. This is beyond the scope of this paper, and we used these (probably) overestimated uncertainties to give the corresponding bounds in the amplitude.
Fig. 8 Corrections applied to the uncertainties in amplitude and inclination from Fig. 2. The row titled “new” was produced using the curves from the right panel of Fig. 7 on MC1 simulations. 

Open with DEXTER 
In the same manner that was corrected in Fig. 8, we corrected for the MC2 calculations. The predicted improvement factors of the amplitude in the bottomright panel of Fig. 6 already include the corrections mentioned above.
4.2. Use of FIM at S/N ~ 15
Another problem regarding the trustworthiness of the Fisher matrix is the limit of the S/N at which one can trust the values from VCM matrix. This has been addressed in Nicholson & Vecchio (1998), where the authors have performed Fisher matrix studies in predicting the uncertainties of the parameters of the coalescing binaries. They have compared these predictions with calculations from Bayesian uncertainties, and their main finding is that above an S/N of 15, the uncertainties from the Bayesian calculations converge to Fisher matrix predictions. At an S/N of 10, the analytic predictions of Fisher matrix deviate by about 6% and at a lower S/N of 7, they increase to about 25%. For MC1, the lowest S/N we considered is above ~ 15, as indicated in Fig. 1, and consequently the results following MC1 are reasonable. However, real observations for eLISA will of course also contain lower S/Ns.
5. Conclusion
We presented a followup study of correlations between parameters that characterise Galactic binaries that will be observed by eLISA in the frequency range 10^{4} Hz < f < 1 Hz. In Paper I we investigated the correlations between parameters by considering three of the verification binaries (which are already known optically and are thus guaranteed sources for eLISA), and we explained these correlations and parameter uncertainties as a function of the inclination of the system. In this paper, we addressed the effect of all position and orientation parameters on the uncertainties and the possible strong correlations between them. We found that the overall effect is dominated by sky position and explained that dependence. Our main conclusions are:

1.
After inclination, the ecliptic latitude of the source has thestrongest influence on its parameter uncertainties, mainly due toits dependence on the S/N and the correlations that can becomevery strong at the poles. For a given source distance, the parameteruncertainties (except for σ_{β}) are larger for a source at the ecliptic pole than for a source at the ecliptic plane. This is because a signal received from a source at the plane is stronger than that from a source at the pole.

2.
In general, the uncertainty in sky location σ_{Ω} is large, i.e. in the order of square degrees. By constraining λ and β from optical data, the uncertainties in ψ, φ, ι, and in the GW analysis can be improved by a factor of ~ 2. This factor is dependent on the inclination and the ecliptic latitude of the system.

3.
By constraining ι in addition to the sky location, the improvement in is remarkably large (factors of 10 depending on the S/N of the source). Systems with ι = 30° can have improvement factors as high as 60.

4.
The remaining angular parameters λ and ψ do not influence the correlations or uncertainties significantly.

5.
The uncertainties in amplitude for lowinclination systems for which the corresponding uncertainties in the inclinations are unphysical (i.e. beyond their physical bounds) are also exaggerated. Fisher matrix predictions for these uncertainties have to be corrected for such systems.
Note that choosing cosι instead of ι does not alleviate this problem since cosι is also a bounded quantity and there is no formalism in FIM to handle bounded quantities, whether cyclic or not. In fact, for the FIM calculations we evaluate uncertainties in cosι and sinβ, while in quoting our results they are converted to uncertainties in ι and β by scaling them as σ_{ι} = σ_{cosι}/sinι and σ_{β} = σ_{sinβ}/cosβ, since these quantities are more intuitive.
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.
References
 AmaroSeoane, P., Aoudia, S., Babak, S., et al. 2013, GW Notes, submitted [arXiv:1201.3621] (In the text)
 Błaut, A. 2011, Phys. Rev. D, 83, 3006 (In the text)
 Cornish, N. J., & Larson, S. L. 2003, Phys. Rev. D, 67, 3001 (In the text)
 Cornish, N. J., & Rubbo, L. J. 2003, Phys. Rev. D, 67, 2001 [NASA ADS] (In the text)
 Cutler, C. 1998, Phys. Rev. D, 57, 7089 [NASA ADS] [CrossRef] (In the text)
 Nelemans, G. 2009, Class. Quantum Grav., 26, 4030 [NASA ADS] [CrossRef] (In the text)
 Nelemans, G., Yungelson, L. R., Portegies Zwart, S. F., & Verbunt, F. 2001, A&A, 365, 491 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Nicholson, D., & Vecchio, A. 1998, Phys. Rev. D, 57, 4588 [NASA ADS] [CrossRef] (In the text)
 Nissanke, S., Vallisneri, M., Nelemans, G., & Prince, T. A. 2012, ApJ, 758, 131 [NASA ADS] [CrossRef] (In the text)
 Peterseim, M., Jennrich, O., Danzmann, K., & Schutz, B. F. 1997, Class. Quantum Grav., 14, 1507 [NASA ADS] [CrossRef] (In the text)
 Roelofs, G. H. A., Rau, A., Marsh, T. R., et al. 2010, ApJ, 711, L138 [NASA ADS] [CrossRef] (In the text)
 Rogan, A., & Bose, S. 2006, unpublished [arXiv:astroph/0605034] (In the text)
 Shah, S., van der Sluys, M., & Nelemans, G. 2012, A&A, 544, A153 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Vallisneri, M. 2008, Phys. Rev. D, 77, 2001 [NASA ADS] [CrossRef] (In the text)
Appendix A: Correlations as a function of β
We found that several correlations between source parameters become amplified at certain ecliptic latitudes and thus are a strong function of β, as shown in Fig. 4. Several of these correlations do not have preference on the inclination of the system, and therefore we show a 2d histogram to show the dependency on the ecliptic latitude. In the figure we show (only) the correlations that become strong, i.e. where the normalised correlation  c_{ij}  > 0.5. Some of these correlations have been explained in Paper I and are only dependent on the inclination of the system; we show them here to confirm that claim. The remaining correlations and their dependency on the latitude can be understood by considering the orientation of the binary relative to the detector. We comment on or explain each of them below:

1.
: explained in Sect. 3.1.

2.
c_{φψ} and c_{φf}: similarly to , these are strongly influenced by the source’s inclination. The correlation c_{φψ} in the bottomleft panel has a distribution that looks very similar to the case of and here too, phase and polarisation are strongly correlated for faceon and mildy faceon systems, since these two angles describe similar properties. For edgeon systems, the two angles are defined on two planes that are perpendicular to each other, making them distinct. The correlation c_{φf} in the topright panel is strong for edgeon systems because a small shift in phase is very similar to a small change in frequency. However, for faceon systems, the phase and frequency become degenerate owing to the uncertainty in φ, which is inflated due to the strong correlation with ψ. There is a slight dependency of this correlation on sinβ for the edgeon orientations since the correlation has a larger spread at the ecliptic poles. This is due to the spread of the c_{φψ} for ι > 80° systems where the spread in the phase at the poles also inflates the frequency. However, ι < 60°, c_{φψ} for has no such spread, and thus the distribution of c_{φψ} vs. sinβ for these systems lies in a narrow band throughout.

3.
c_{ψλ}: this can be understood by considering a binary with ι = 90° close to one of the poles; rotating the binary about its lineofsight is equivalent to changing its ψ and λ because they are defined on the same plane (perpendicular to the lineofsight). This also holds for a faceon binary, and thus if a source is close to the poles, ψ and λ can be highly correlated. The spread of this correlation is most likely due to the fact that λ is not well defined at the poles and thus the inflated uncertainty in the longitude implies a spread of correlations at the poles.

4.
c_{φλ}: this correlation is explained by the combination of c_{ψλ} and c_{φψ}, which are explained above. Since ψ and λ can be highly correlated at the poles and φ and ψ are strongly correlated for faceon systems, we found that φ and λ should also be highly correlated at the poles, at least for the systems with (mildly) faceon orientations. For the edgeon systems, this strong correlation can also exist because even though c_{φψ} is not as strong in this case, the values of this correlation (as pointed out above) spread at the poles for ι = 90°. This spread is responsible for the spread of c_{ψλ} around the poles for edgeon systems.

5.
c_{ιλ} and : a strong correlation between inclination and the longitude is intuitively not as clear. However, this correlation combined with for faceon systems would explain the correlation between amplitude and longitude, .

6.
c_{φι}: like above, the correlation between phase and inclination follows from the combination of c_{φλ} and c_{ιλ}, both of which are stronger close to the poles.
All Figures
Fig. 1 Twodimensional distributions of S/N as a function of cosι, sinβ, λ, and ψ. Faceon, mildly faceon, mildly edgeon, and edgeon systems (whose exact bounds are specified in Sect. 2) can be read off from the topleft panel. In the other panels, the trend from faceon to edgeon systems follows from the top of the distribution to the bottom. 

Open with DEXTER  
In the text 
Fig. 2 Twodimensional distributions of the uncertainties in the seven GW parameters of the galactic binary systems as a function of sinβ from MC1 as described in the text. Units of the angles are in radians. 

Open with DEXTER  
In the text 
Fig. 3 Uncertainties in all seven parameters of the systems as a function of S/N from MC2. Open circles, squares, and filled triangles represent binary systems at fixed inclinations of 90°,60°, and 30° respectively. The size of the symbol represents the magnitude of  β  for each simulated binary. In the bottom panels the different groups of symbols overlap. 

Open with DEXTER  
In the text 
Fig. 4 Twodimensional distributions of the normalised correlations as a function of sinβ from MC2. We show here only the strong (normalised) correlations, meaning  c_{ij}  > 0.5. 

Open with DEXTER  
In the text 
Fig. 5 Uncertainty ratios of the parameters between analyses using a 7 × 7 VCM and a 5 × 5 VCM, where the sky position is fixed in the latter, as a function of sinβ (from the MC2 and MC3 simulations). Filled triangles, squares and open circles represent systems at fixed inclinations of 30°,60° and 90° respectively. The size of the symbol for each binary represents the S/N of the system. 

Open with DEXTER  
In the text 
Fig. 6 Same as in Fig. 5 but now we assume that in addition to sky position we also know the inclination of the binary. The uncertainty ratios of the parameters are between 7 × 7 VCM (MC2) and 4 × 4 VCM (MC4). The improvement factors in amplitude include corrections due to overestimated uncertainties in the amplitudes of faceon systems (see Sect. 4.1). Note that the last panel is plotted in log scale in order to highlight the relatively large improvements for the systems with smaller inclinations. 

Open with DEXTER  
In the text 
Fig. 7 Left: 1σ uncertainty in the inclinations vs. source inclinations for two values of S/N: plusses (outer curves) correspond to S/N = 10, stars (inner curves) correspond to S/N = 53. The curves above the dotted line correspond to ι + σ_{ι}, and those below to ι − σ_{ι}. The thick horizontal dashed and dasheddotted lines are the physical limits of ι: below the line of ι ± σ_{ι} = 0° and above the line of ι ± σ_{ι} = 90°, σ_{ι} can be truncated, since that value of ι is bounded between [0°,90°]. Right: “unphysical ι” versus S/N for the systems from MC1. The thick dasheddotted line calculated from the condition ι − σ_{ι} = 0°, and the thick dashed line from the condition ι + σ_{ι} = 90° demonstrated in the left panel. At an S/N of 10, the in open circles show a lower limit and the filled circles show an upper limit at the intersections with the two curves from the conditions above. 

Open with DEXTER  
In the text 
Fig. 8 Corrections applied to the uncertainties in amplitude and inclination from Fig. 2. The row titled “new” was produced using the curves from the right panel of Fig. 7 on MC1 simulations. 

Open with DEXTER  
In the text 