Next Contents Previous

5.1. Sky Maps of VELMOD Residuals

VELMOD does not assign galaxies a unique distance (Section 2.2.1). Thus, there is no unique measure of the amount by which their observed and predicted peculiar velocities differ. However, there is a well-defined, expected apparent magnitude for each object,

Equation 21 (21)

where P(m|eta, cz) is given by equation (11). Similarly, the rms dispersion about this expected value is

Equation 22 (22)

Note that Deltam is not equal to the TF scatter; it also includes the combined effects of velocity noise, peculiar velocity gradients, and density changes along the line of sight. At large distances, Deltam tends toward sigmaTF (although dispersion bias can make it smaller; cf. Willick 1994). With the above definitions, one can define a normalized magnitude residual for each galaxy,

Equation 23 (23)

The normalized residual has the virtue of having unit variance for all objects. In contrast, the variance of the unnormalized magnitude residual m - E(meta, cz) depends on distance (velocity noise is more important for nearby objects), while the variance of a peculiar velocity residual formed from the magnitude residual (eq. [24]) grows with distance. The normalized magnitude residual deltam is a measure of the correctness of the IRAS velocity model. If, in a given region, deltam > 0 in the mean, galaxies in that part of space must be more distant than IRAS predicts them to be, i.e., they have negative radial peculiar velocity relative to the IRAS prediction. Regions in which deltam < 0 in the mean have positive radial peculiar velocities relative to IRAS.

We will use the normalized magnitude residual below in our quantitative analysis of residuals, but let us first visualize this residual field on the sky, by converting deltam into the corresponding radial peculiar velocity residual deltau. If we were to do this for each galaxy individually, the ~ 20% distance errors due to TF scatter would completely hide the systematic departures from the IRAS model. Instead, we will compute smoothed velocity residuals. This procedure is the most well behaved if we first smooth deltam and then convert the result into deltau. We first place each galaxy at the distance d assigned it (13) by the IRAS velocity model; this is a redshift-space distance, so our calculation is unaffected by Malmquist bias. Then, for each galaxy i, we compute a smoothed residual deltam,is as the weighted sum of the residuals deltam of itself and its neighbors j, where the weights are wij = exp (-dij2 / 2Si2), and dij is the IRAS-predicted distance between galaxies i and j. We take the smoothing length Si to be Si = di / 5. The smoothed residual deltam,is is converted into a smoothed velocity residual according to

Equation 24 (24)

where Deltami is given by equation (22). The quantity fi is given by exp (-Deltai2 / 2), where Deltai = 0.46 Deltami / (sumj wij)1/2; it guarantees that deltauis, which is lognormally distributed, has an expectation value of zero if Deltami (which is normally distributed) does also (cf. Willick 1991, Section 6.3, for details).

In Figures 10, 11, and 12, we plot VELMOD velocity residuals on the sky for the redshift ranges 0-1000 km s-1, 1000-2000 km s-1, and 2000-3000 km s-1, respectively. In each figure, the top panel shows residuals from the betaI = 0.6 (no-quadrupole) fit, and the bottom panel shows residuals from the betaI = 0.5 (quadrupole-modeled) fit, the VELMOD runs closest to the maximum likelihood value of betaI for each case. The plots reveal why the addition of the quadrupole results in a large increase of likelihood. In each redshift range, the no-quadrupole fits show coherent negative velocity residuals both in the Ursa Major region (l appeq 150°, b appeq 65°) and at b appeq -60°, l ltapprox 30°, and l gtapprox 330°. In both of these regions, the addition of the quadrupole reduces the amplitude of the residuals significantly. In other parts of the sky, smaller but still significant coherent residuals are reduced with the addition of the quadrupole. This shows that the pattern of departure from the pure IRAS velocity field is well modeled by a quadrupolar flow of modest amplitude and therefore has the simple physical interpretation we discussed in Section 4.4.

Figure 10

Figure 10. VELMOD velocity residuals plotted on the sky in Galactic coordinates, for objects with 0 < czLG leq 1000 km s-1. The top panel is for the betaI = 0.6 run without the quadrupole. The bottom panel is for the betaI = 0.5 run with the quadrupole modeled. The open circles indicate objects that are inflowing relative to the IRAS prediction; the stars indicate outflowing objects.

Figure 11

Figure 11. Same as Fig. 10, but for objects with 1000 < czLG leq 2000 km s-1.

Figure 12

Figure 12. Same as Fig. 10, but for objects with 2000 < czLG leq 3000 km s-1.

In the bottom panels, it is difficult to find any well-sampled region within 2000 km s-1 where deltau gtapprox 100 km s-1. This is all the more remarkable because the TF errors themselves are of order 300 km s-1 per galaxy at a distance of 1500 km s-1. Figure 12 does show several high-amplitude residuals. However, at 2500 km s-1, the TF residual for a single object is 500 km s-1, so when the effective number of galaxies per smoothing length is only a few, velocity residuals of several hundred km s-1 are expected from TF scatter only. In well-sampled regions, one sees that in general deltau ltapprox 150 km s-1, the only exception being a patch of large (gtapprox 250 km s-1) positive residuals at l appeq 330°, b appeq -20°. In the b > 0° part of the Great Attractor region at l appeq 300°, the residuals are less than 100 km s-1, even in this highest redshift shell. This is significant, given the often heard claims that the IRAS model cannot fit the observed flow into the Great Attractor. Although some residual coherence is apparent, we will demonstrate in Section 5.2 that this is largely because of the smoothing.

In Figures 13, 14, and 15, we again plot VELMOD residuals on the sky for the three redshift ranges, now for the two values of betaI most strongly disfavored by the likelihood statistic in the range studied, betaI = 0.1 (top panels) and betaI = 1.0 (bottom panels). In each plot, the quadrupole of Figure 4 has been included. These plots, which should be compared with the bottom panels of Figures 10, 11, and 12, demonstrate why very low and high betaI do not fit the TF data well. In each redshift range, these models exhibit large coherent residuals. For betaI = 0.1, we see large negative peculiar velocities relative to IRAS in the Ursa Major region at cz leq 2000 km s-1. Indeed, the residual plot for betaI = 0.1 (with quadrupole included) shows many of the same features as the no-quadrupole model with betaI = 0.6, because the IRAS field itself contributes some of the needed quadrupole. However, the IRAS contribution scales with betaI and is thus inadequate at low betaI. At betaI = 1.0, many of the systematic residuals associated with the quadrupole are gone, especially in Ursa Major. However, other regions show highly significant residuals: at l appeq 150°, b appeq -20°, and cz leq 1000 km s-1, for example, one sees negative peculiar velocity residuals of amplitude gtapprox 200 km s-1, which is significant at such small distances. In the same redshift range, at l = 270°-360°, b < 0°, there are positive velocity residuals of amplitude gtapprox 150 km s-1. These regions exhibit much smaller residuals in the betaI = 0.5 model.

Figure 13

Figure 13. Same as Fig. 10, except now the results for betaI = 0.1 and betaI = 1.0 are shown. In each case, the quadrupole is the same as it was for the best-fit model (betaI = 0.5).

Figure 14

Figure 14. Same as Fig. 13, but for objects with 1000 < czLG leq 2000 km s-1.

Figure 15

Figure 15. Same as Fig. 14, but for objects with 2000 < czLG leq 3000 km s-1.

In the higher redshift shells, the poor fit of the betaI = 1.0 model is evidenced chiefly in the direction of the Great Attractor (GA) (l appeq 300°, b appeq 20°). For 1000 < cz leq 2000 km s-1, this model predicts much too large positive peculiar velocities, so that the data exhibit inflow relative to the model. In the highest redshift bin, the betaI = 1.0 model exhibits both positive and negative velocity residuals of high amplitude in the GA direction; residuals of both signs are seen in this region for betaI = 0.5 as well, but they are of much smaller amplitude (bottom panel of Fig. 12). The betaI = 0.1 model, on the other hand, predicts too small positive peculiar velocities in the GA direction at the highest redshifts. Indeed, note that in the 2000 < czLG leq 3000 km s-1 shell, nearly all data points exhibit outflow relative to the betaI = 0.1 IRAS predictions, whereas at lower velocities the residuals typically indicate inflow. This global mismatch is more general than the insufficient quadrupole mentioned in the previous paragraph, showing that low betaI could not yield a good fit even if we were to give VELMOD full freedom in fitting the quadrupole at all betaI.

Although sky plots of residuals argue in favor of the betaI = 0.5 plus quadrupole model, the residuals from that model are not manifestly negligible. We will address this issue quantitatively below. For now, however, we can demonstrate qualitatively that the residuals seen in the betaI = 0.5 plus quadrupole model are not unexpected by comparing them with the mock catalogs, for which the IRAS velocity predictions are known to be a good fit. Figure 16 plots VELMOD velocity residuals with respect to betaI = 1 (the correct value) for a single mock catalog. The same three redshift ranges used for the real data are shown. The mock catalog residuals are comparable in amplitude and apparent coherence to the real data. Generally speaking, velocity residuals in well-sampled regions are ltapprox100 km s-1 within 1000 km s-1 and are ltapprox200 km s-1 at larger distances. One also sees apparent coherence in the mock catalog residual map, as was the case with the real data. The similar amount of apparent coherence in the real and mock data indicates that the former is not the result of a poor fit. The apparent coherence in the residual sky maps is an artifact of the smoothing used to generate them, as we show in the next section.

Figure 16

Figure 16. VELMOD velocity residuals for a single mock catalog run using betaI = 1.0, the true value for the mock catalog. (The particular simulation used had a maximum likelihood value of betaI = 0.963.) The three panels show residuals for the three redshift ranges used in analyzing the real data.

13 We take d to be the "crossing point distance" w defined in Section 2.2.2. In the case of triple-valued zones, we take the central distance. Back.

Next Contents Previous