text:
# Abundant molecular gas in the intergalactic medium of Stephan’s Quintet

## 1 Introduction

Stephan’s Quintet (Hickson Compact Group 92; hereafter SQ) is one of the best studied examples of a Hickson Compact Group. It contains four interacting galaxies, NGC 7319, NGC 7318a, NGC 7318b, and NGC 7317. A fifth galaxy, NGC 7320, happens to be a foreground object. An optical image of the group is presented in Fig. 1, (see also the excellent image in Arp & Lorre 1976). One of its most striking properties is that the major part of the gas is in the intragroup medium, most likely the result of interactions in the past and present.

A plausible scenario for the dynamical history of SQ is presented by Moles et al. (1997). They suggest that a few times 108 yr ago the group experienced a collision with NGC 7320c, a galaxy $``$4 arcmin to the east of NGC 7319 but with a very similar recession velocity (6583 km s-1, Sulentic et al. 2001) to the other galaxies in SQ. This collision removed most of the gas of NGC 7319 towards the west and east, and produced the eastern tidal tail which connects to NGC 7319. Presently, the group is experiencing another collision with the “intruder” galaxy NGC 7318b which strongly affects the interstellar medium (ISM) removed during the first collision. NGC 7318b has a recession velocity of 5765 km s-1, in contrast to the other members of the group which have velocities close to 6600 km s-1. Sulentic et al. (2001), in a multi-wavelength study of the group, confirm this scenario and suggest that the group has been visited twice by NGC 7320c. The first collision created the very faint tidal arm east of the interloper NGC 7320, whereas the second interaction produced the tidal arm which stretches from NGC 7319 eastwards.

This violent dynamical history has induced star formation at various places outside the individual galaxies. ISOCAM mid-infrared and H$`\alpha `$ observations have revealed a starburst region (object A in Xu et al. 1999, hereafter called SQ A) at the intersection of two faint optical arms stemming north from NGC 7318a/b. Several knots of star formation are visible in the tidal arm extending from NGC 7319 to the east. At this position, identified as B by Xu et al. (1999) and hereafter called SQ B, there is also mid-infrared and H$`\alpha `$ emission, although much weaker than in SQ A.

Using HST observations, Gallagher et al. (2001) found 115 candidate star clusters, most of them distributed among the tidal debris of SQ. From color-color diagrams they estimated their ages ranging from 2–3 Myr up to several Gyr. The distribution of ages sheds light on the star formation history in SQ: The youngest star clusters (with ages of less than 10 Myr) are found in SQ A and south of NGC 7318a/b, while somewhat older star clusters, with ages between 10 and 500 Myr are in the young tidal tail and around NGC 7319. This is consistent with the picture that the eastern tidal tail was produced in a previous interaction whereas the collision and star formation around SQ A is ongoing. Mendes de Oliveira et al. (2001) have carried out Fabry-Perot observations of the H$`\alpha `$ emission in SQ A and south of NGC 7318a/b, from which they determined the velocity curves for several H$`\alpha `$ emitting regions. SQ B was not covered by their observations. They identified seven Tidal Dwarf Galaxy (TDG) candidates, selected as objects exhibiting a velocity gradient compatible with rotation and possessing an $`L_\mathrm{B}/L_{\mathrm{H}\alpha }`$ ratio consistent with that of a dwarf galaxy.

Abundant atomic hydrogen is present to the east of the 3 central galaxies, around SQ A and south of NGC 7318a (Shostak et al. 1984, Williams et al. 2002, see their Fig. 5). Williams et al. (2002) suggest that the atomic gas towards the east of SQ consists of two tidal features, each connected to one of the optical tidal arms (the old and the new one, see Fig. 1). Molecular gas has been found at the position of SQ A by Gao & Xu (2000) with BIMA and by Smith & Struck (2001) with the NRAO 12m radio telescope. Both the molecular and the atomic gas in this location present two velocity components, centered at about 6000 and 6700 km s-1, which implies that they originate from different galaxies. Braine et al. (2001) in a study of the CO emission of TDGs detected 2.9$`\times `$108 M of molecular gas at SQ B. In order to follow up on these detections and to elucidate the extent, origin and fate of the CO in the intergalactic medium (IGM) in SQ, we embarked on a single dish survey, the results of which we report in this paper.

Details of our observations can be found in Section 2. In Section 3 we present the results of these observations, as well as a comparison between the molecular gas and emission at other wavelengths. We discuss our findings in Section 4 and summarize our conclusions in Section 5. Following Williams et al. (2002), who assumed a recession velocity of 6400 km s-1 and $`H_0`$=75 km s-1 Mpc-1, we adopt a distance of 85 Mpc, in which case 10′′ correspond to 4.1 kpc.

## 2 Observations

We observed the CO(1–0) and CO(2–1) lines at 115 and 230 GHz in SQ in July 2001 with the IRAM 30-meter telescope on Pico Veleta. Dual polarization receivers were used at both frequencies with the 512 $`\times `$ 1 MHz filterbanks on the CO(1–0) line and the autocorrelator or 256 $`\times `$ 4 MHz filterbanks on the CO(2–1). The observations were done in wobbler switching mode with a wobbler throw of 120′′ in azimuthal direction. Pointing was monitored on nearby quasars every 60 – 90 minutes, the rms offset being $`3`$′′. At the beginning of the observations and every time the frequency was changed, the frequency tuning was checked by observing Sagittarius B2. System temperatures were generally quite good, 150–200 K at 115 GHz and 200–300 K at 230 GHz on the $`T_\mathrm{A}^{}`$ scale. The IRAM forward efficiency, $`F_{\mathrm{eff}}`$, was 0.95 and 0.91 at 115 and 230 GHz and the beam efficiency, $`B_{\mathrm{eff}}`$, was 0.75 and 0.54, respectively. Our half-power beam size was 21′′ at 115 GHz and 11′′ at 230 GHz. All CO spectra and luminosities are presented on the main beam temperature scale ($`T_{\mathrm{mb}}`$) which is defined as $`T_{\mathrm{mb}}=(F_{\mathrm{eff}}/B_{\mathrm{eff}})\times T_\mathrm{A}^{}`$. For the data reduction, the spectra were summed over the individual positions and a constant continuum level was subtracted.

We observed three different regions of Stephan’s Quintet: (i) the region around starburst SQ A, where strong mid-infrared and H$`\alpha `$ emission had been detected (Xu et al. 1999), (ii) the area around a star forming region in the eastern tidal arm of NGC 7319 (SQ B), as well as one position (SQ B-North) at the northern tip of the HI cloud covering SQ B (Williams et al. 2002), and (iii) several positions to the south and north of NGC 7318a/b, observed by Mendes de Oliveira et al. (2001) in H$`\alpha `$. We call these positions SQ C$`n`$ where $`n`$ is the number that Mendes de Oliveira et al. (2001) gave to the corresponding region. The locations of our pointings are shown in Fig. 1.

In SQ A and SQ B the spacing of the individual pointings is 10′′ so that the CO(2–1) map is undersampled. Note that our observations are centered at the positions called SQ A and SQ B as in Xu et al. (1999) but that they comprise a much larger area.

## 3 Results

### 3.1 Molecular gas distribution

In Fig. 2 and Fig. 3 we present a mosaic of the CO(1–0) spectra observed around SQ A and SQ B. The noise in the spectra is about 6 mK (SQ A) and 5 mK (SQ B) at a velocity resolution of 2.66 km s-1. The CO(2–1) emission was also detected at several positions where the CO(1–0) line was visible (see Fig. 4). The quality of the CO(2–1) map of SQ A is much poorer, so that no useful information can be derived from the individual spectra. Therefore we do not show it here. In Table 1 the integrated intensities and the derived molecular gas masses for the two regions are listed. All molecular gas masses in this paper are calculated using the galactic conversion factor of N(H2)/I$`{}_{CO}{}^{}=2\times 10^{20}`$ cm -2(K km s-1)-1, yielding:

$$M_{\mathrm{H}_2}[M_{}]=75I_{\mathrm{CO}}D^2\mathrm{\Omega },$$
(1)
where $`I_{\mathrm{CO}}`$ is the velocity integrated CO line intensity expressed in K km s-1, $`D`$ is the distance in Mpc and $`\mathrm{\Omega }`$ is the area covered by the observations in arcsec2 (i.e. $`\mathrm{\Omega }=1.13\mathrm{\Theta }^2`$ for a single pointing with a Gaussian beam of FWHM $`\mathrm{\Theta }`$). We direct the reader to Braine et al. (2001) for a more detailed derivation of the molecular gas mass. When referring to molecular data from other authors, we scale their values to this conversion factor throughout this paper.

The CO emission in SQ A is very extended. CO lines at a central velocity of about 6000 km s-1 are visible in individual spectra over an area of about 40$`{}_{}{}^{\prime \prime }\times `$ 50′′ (Fig. 2). The total molecular gas mass at this velocity, derived by integrating over the entire area observed, is $`2.2\times 10^9`$ M. If we were to limit the integration to the positions where the CO(1-0) line is clearly visible, we would obtain an integrated intensity of $`1.2\pm 0.06`$ K km s-1 and a total molecular gas mass of $`1.3\times 10^9`$ M. The error in these masses is about 10%, partly due to the error in the integrated intensity (5%) and partly due to the error in the total area which we estimate to be about 5 – 10%. Thus, by summing over those spectra with clearly detected signal we apparently miss a faint, extended component which accounts for $`40\pm 6`$% of the total mass present. In the integrated CO spectrum a weaker component at 6700 km s-1 is visible as well (Fig. 5). This component is not apparent in the individual spectra, and therefore we do not include this velocity range in Fig. 2. Its total molecular mass, integrated over the whole area, is $`8.6\times 10^8`$ M.

In SQ B the molecular gas has been detected over a region of about 60$`{}_{}{}^{\prime \prime }\times `$ 40′′ (Fig. 3). It is associated with the optically visible star forming region – the maximum of the CO emission coincides with it – but is more extended. The CO emission includes the region at the very tip of the optical tidal tail where weaker H$`\alpha `$ and a stellar condensation with blue colour ($`BV=0.2`$; Sulentic et al. 2001) has been found. The total molecular gas mass in this area ($`7.0\times 10^8`$ M, see Tab. 1) is more than twice as much as that found by Braine et al. (2001) who observed only the central position.

We did not detect CO at the majority of the positions of the TDG candidates observed by Mendes de Oliveira et al. (2001) south and north of NGC 7318a/b. We got a marginal detection of $`I_{\mathrm{CO}}=0.5`$ K km s-1 when averaging the positions SQ C16/17 and SQ C18/19 which we used as an upper limit. No CO was detected in the northern part of the eastern HI tail (SQ B-North). This position was chosen because it coincides with a local peak of the HI emission and shows H$`\alpha `$ emission (Sulentic et al. 2001). The rms noise of the positions with non-detections are 3–4 mK for CO(1-0) and 4–6 mK for CO(2-1), at a velocity resolution of 5 km s-1.

### 3.2 Kinematics and excitation

In Fig. 5 and Fig. 6 we show the average CO(1–0) spectra in both regions, overlaid over the HI emission (Williams et al. 2002), averaged over the same region. The line shapes in both regions are different: In SQ A the CO line is very broad, with a FWHM of $`6080`$ km s-1 and line wings extending out to more than 200 km s-1, showing that in this region large velocity gradients are present, whereas the line width in SQ B is much narrower, with a FWHM of $`3040`$ km s-1.

The HI gas of SQ A has two distinct velocity components, one at 6000 and another at 6700 km s-1. The CO and HI lines generally agree well, with the exception of a slight velocity offset between the low-velocity component of the HI line, centered at 6000 km s-1, and the CO line, centered at 6030–6040 km s-1. Nevertheless, the line width, shapes and the relative intensity of the components at 6000 and at 6700 km s-1 are the same for the HI and CO spectra. Smith & Struck (2001) find, using the NRAO 12m telescope, that both velocity components have similar strengths. The reason for why they underestimate the strength of the line at 6000 km s-1is most likely because they cover a smaller area. This can be seen when comparing the molecular gas masses that they derive for their components. Whereas their value for the component at 6700 km s-1, M$`_{\mathrm{H}_2}`$$`=7\times 10^8`$ M, agrees well with ours, their value for the component at 6000 km s-1, M$`_{\mathrm{H}_2}`$$`=7.9\times 10^8`$ M, is much lower. This can be understood if the 6700 km s-1 component is more centrally concentrated, as is the case for the HI emission (Williams et al. 2002), and as a result is covered by their beam, but the 6000 km s-1 component is not. From their Fig. 1g it is indeed obvious that their beam is smaller than the area mapped by us.

The CO emission in SQ B has a central velocity of 6625 km s-1. There is very good agreement between the CO and HI lines (from Williams et al. 2002) with respect to the line width, central positions and special features, such as the blue wing.

To better understand the dynamical behavior of the CO emitting gas we created several position-velocity diagrams using our data. Our analysis of SQ B only gave marginal indications for a weak velocity gradient of about $`1020`$ km s-1 along the east-west direction. In SQ A we did find evidence for a slope in velocity along a position angle of $`40`$ degrees for the velocity range around 6000 km s-1 (Fig. 7). Moreover, for a cut perpendicular to it, no slope was seen. This is very similar to what we found for the atomic gas using the HI data of Williams et al. (2002). A similar trend has also been seen in the H$`\alpha `$ data of Mendes de Oliveira et al. (2001) for their regions 2,3,4 and 5. The position angle of the H$`\alpha `$ velocity is the same as the one we find in CO and HI, but their observed velocity range is slightly smaller, between 5938 km s-1 and 6010 km s-1, whereas the CO data shows a gradient from about 5970 km s-1 to 6130 km s-1. This is probably due to the smaller extent of only about 30′′ of the H$`\alpha `$ observations. Normalized to the same spatial extent the H$`\alpha `$ gradient is comparable to the gradient in CO within the uncertainties.

In Fig. 8 and 9 we show the CO(2–1) lines, averaged over the same area as the corresponding CO(1–0) lines. In order to compute the CO(2–1)/CO(1–0) integrated intensity ratio in SQ A and SQ B, the CO(2–1) spectra have been convolved to the spatial resolution of the CO(1–0) beam. To minimize the error of this ratio, we have included only spectra at positions with a S/N $`>3\sigma `$ of the CO(1-0) line. In SQ A, the CO(2–1)/CO(1–0) ratio is $`0.69\pm 0.16`$ and in SQ B $`0.56\pm 0.13`$. Within the errors, there is no significant difference in the line ratios. They are slightly lower than the average found for a mixed sample of nearby galaxies (CO(2–1)/CO(1–0)$`=0.89\pm 0.06`$, Braine & Combes 1992). Better CO(2–1) data are necessary in order to judge whether this difference is significant.

### 3.3 Comparison of the CO and HI distribution

In spite of the practically identical kinematics, the atomic and molecular gas distributions show differences both in SQ A and in SQ B. In SQ A, the CO lines are wide and have low maxima. As a result, it is difficult to display the velocity integrated intensities of the CO emission because of their large uncertainties. We therefore show in Fig. 10 a comparison of the HI emission around 6000 km s-1 with the observed positions in CO. The positions where we clearly ($`>3\sigma `$) detected a line are indicated by a cross and the observed positions as circles. Note that, as mentioned earlier, also at the positions without a clear CO detection some faint emission must be present. From Fig. 10 it is clear that the overall distribution of the CO emission is different from the HI emission. Even though some CO marginally coincides with the HI peak, most of it is located to the south-east along a faint HI ridge.

In Fig. 11 we display the contours of the integrated CO intensity of SQ B overlayed onto contours of an HI integrated surface brightness map. The maxima of the CO and HI emission are close, but do not exactly coincide. Instead, the CO peaks slightly (10 – 15 ′′) west from the HI ridge, towards a steep HI gradient, suggestive of compression in the gas which implies that this over-density in HI may have provoked the formation of molecular gas. The molecular-to-atomic gas ratio is highest at the peak of the CO emission and west of it. This means, under the assumption that the molecular gas is being formed in situ, as has been suggested for TDGs (Braine et al. 2001), that this compressed region is the place of the most efficient transformation of HI into H2.

The HI mass in SQ A, integrated over the area mapped in CO, is $`1.6\times 10^9`$ M($`9\times 10^8`$ M) for the velocity component at 6000 km s-1 (6700 km s-1) (see Tab. 1). This means that the molecular gas mass is comparable to the atomic gas mass with a ratio of M$`_{\mathrm{H}_2}`$/MHI= 1.4 (1.0). The HI mass of the entire area around SQ A is $`3.1\times 10^9`$ M, which is of the same order as the total mass of the gas in molecular form. The atomic gas mass in SQ B, integrated over the region where CO was detected, is $`1.4\times 10^9`$ M. This yields a molecular-to-atomic gas mass ratio of 0.5, higher than the values of less than about 0.2 found for other TDGs (Braine et al. 2001).

### 3.4 Comparison of the CO and H$`\alpha `$ distribution

The H$`\alpha `$ emission in SQ A consists, as do the HI and CO data, of two velocity components, centered around 6000 km s-1 and 6700 km s-1, respectively (Sulentic et al. 2001, Iglesias-Páramo & Vílchez 2001, Xu et al. 1999). The distribution of both velocity components is similar in the regions that our CO observations covered. Using the data of Xu et al. (1999), taken with narrowband filters, we present in Fig. 10 the spatial distribution of the H$`\alpha `$ emission at 6000 km s-1 and our CO detections. We observe that the CO gas follows closely the H$`\alpha `$ emission in the north-south direction, but surprisingly little CO is seen close to the brightest H$`\alpha `$ peak in the starburst region at SQ A.

Similarly, we present in Fig. 11 the H$`\alpha `$ image of SQ B overlayed with the contours of the molecular gas. The H$`\alpha `$ emission is situated 5′′ north of the maximum of the CO emission and approximately 15′′ north-west from the maximum of the HI emission. We conclude that the H$`\alpha `$ and CO emission coincide within the limited spatial resolution (21′′) of our data. Also, the velocities agree, with the velocity of the CO line being 6625 km s-1 and that of the H$`\alpha `$ emission 6617 km s-1 (Sulentic et al. 2001). Altogether, this is pointing to a causal connection between molecular gas and recent star formation.

In Tab. 1 we list the molecular gas mass divided by the SFR (derived from the H$`\alpha `$ luminosity) which gives the gas consumption time (i.e. the inverse of the star formation efficiency). The values we derive for SQ A and SQ B are typical to those of spiral galaxies, which have been estimated to lie between 4 and 40 $`\times 10^8`$ yr (Kennicutt 1998).

### 3.5 Molecular gas fraction and star formation in the non-detections

We observed various positions (see Fig. 1) south of NGC 7318a/b and north of SQ B in CO without getting a detection. The results are summarized in Tab. 1.

The position north of SQ B (SQ B-North) coincides with the north tip of the HI cloud and was chosen because Sulentic et al. (2001) detected H$`\alpha `$ emission at the same velocity as the HI and suggested that an extragalactic HII region may be forming. Gallagher et al. (2001) identified two starclusters (30 and 31) at the same position. The clusters have blue colours and are therefore young, less than 10 – 500 Myr, and must have formed in situ. At this position we derived an upper limit of M$`_{\mathrm{H}_2}`$/MHI$`0.08`$. This value is much lower than at the other places in SQ where we detected molecular gas, but it is not unusual compared to typical values for TDGs (Braine et al. 2001). The upper limit of M$`_{\mathrm{H}_2}`$/SFR$`<4\times 10^8`$ yr is at the low end of the typical range for spiral galaxies.

Around NGC 7318a/b we observed several of the positions where Mendes de Oliveira et al. (2001) have found H$`\alpha `$ emission (their positions 8, 16/17, 18/19, 20, 21, 22 and 23; note that positions 16 and 17 as well as 18 and 19 fall within the IRAM beam). The authors suggest that some of them (8, 20, 21, 22 and 23) are TDG candidates, mainly based on the detection of a velocity gradient in the H$`\alpha `$ data which they interpret as rotation. We observed in addition to these TDG candidates positions 16/17 and 18/19 because of their relatively high HI column density. The upper limits of the molecular gas fraction, M$`_{\mathrm{H}_2}`$/MHI, are low for some objects, but not unusual compared to spiral galaxies and TDGs. The upper limits for the gas consumption time, M$`_{\mathrm{H}_2}`$/SFR, on the other hand, are, for most objects, much lower than the values found for TDGs (see Braine et al. 2001) or spiral galaxies. This is particularly surprising, as the large size of the CO(1-0) beam of 21′′, corresponding to 8.6 kpc, should comprise the molecular gas associated with the star formation even if the H$`\alpha `$ and CO emission do not coincide exactly, because, e.g. the ongoing star formation has destroyed the molecular gas locally. A low metallicity could explain the non-detections but seems unlikely as we did detect CO at other places in SQ. We checked for the presence of star clusters from Gallagher et al. (2001) at the observed positions in order to find out whether star clusters have been formed since the beginning of the collision with NGC 7318b of less than $`10^8`$ yr ago (Sulentic et al. 2001). Region C23 coincides with an old star-cluster (number 117), with an age of several Gyr, which was therefore born before the start of the latest interaction. Regions C8, C16/17 and C18/19 coincide with star clusters 133, 134, 135, 137, 122, 128, and 130, some of which are younger than 10 Myr. At the other observed positions no star clusters were found within a radius of 10′′. Thus, young star clusters formed after the beginning of the collision have only been found close to the parent galaxy NGC 7318a and not in the TDG candidates C20 – C23.

### 3.6 Comparison of the CO distribution to other tracers

In Fig. 12, 13 and 14 we compare the CO distribution to that of the 15 $`\mu `$m emission observed by ISOCAM (Xu et al. 1999), the X–ray emission from the ROSAT archive, and the 20cm radio continuum emission (Williams et al. 2002), respectively.

The overall distribution of the dust traced by the ISOCAM 15 $`\mu `$m image is similar to the CO emission in both SQ A and SQ B (Fig. 12). In SQ A the 15 $`\mu `$m emission coincides with the maximum in the H$`\alpha `$ and indicates the position of the intergalactic starburst. The CO is slightly offset towards the south and west of the 15 $`\mu `$m peak. The CO emission coincides with faint 15 $`\mu `$m emission at all our observed positions, showing that the presence of molecular gas is accompanied by the presence of dust. In SQ B the 15 $`\mu `$m peak is marginally offset (5-10′′) from the CO maximum and is located towards the peak of the HI surface density.

In the intergalactic medium around SQ A, radio continuum (Williams et al. 2002) and X–ray emission (Pietsch et al. 1997) are present. As previously noted (e.g. Sulentic et al. 2001), there is good agreement between the radio continuum, the X–ray and H$`\alpha `$ emission in the elongated north-south feature north of NGC 7318a/b. All along this ridge we find CO emission (Fig. 13 and 14). CO is also present in a region without X–ray emission towards the north which coincides with the edge of the ISOCAM peak. In SQ B no X–ray emission is detected, but weak ($`2\sigma `$) radio continuum emission can be seen, coinciding with the peak of the CO emission.

The radio continuum emission in both regions is mostly non-thermal (i.e. synchrotron emission) as can be deduced by calculating the expected thermal radio emission from the H$`\alpha `$ emission. The thermal radio emission derived from the H$`\alpha `$ flux yields (following eq. (3) and (4) in Condon 1992), a few percent of the radio continuum flux in SQ A and about 10% in SQ B. This thermal fraction at 1.49 GHz is typical for star-forming galaxies (Condon 1992).

## 4 Discussion

### 4.1 Molecular gas, shocks and star formation in SQ A

The enormous amount of molecular gas that we have found and its large extent are surprising. In the following we will discuss how this molecular gas fits into the current picture of the IGM in SQ A.

The IGM around SQ A shows two velocity components. Sulentic et al. (2001) suggest that the gas at 6700 km s-1 was stripped during the earlier visit of NGC 7320c and part of it has become shock-ionized due to the present interaction with NGC 7318b. In this view, the H$`\alpha `$ emission at 6700 km s-1 traces the shock front. Strong support for this shock hypothesis is the fact that X–ray (tracing hot gas) and radio continuum emission (tracing cosmic ray electrons produced in shocks) are coincident with the H$`\alpha `$ emission. The gas at velocities around 6000 km s-1 originates from NGC 7318b and is being removed from this galaxy in the ongoing collision. The H$`\alpha `$ emissions at 6000 and 6700 km s-1 overlap to a large extent in the region where we observed CO. In HI, the two velocity components have their maxima at the same position, but the component at 6700 km s-1 is spatially much less extended than that at 6000 km s-1 (Williams et al. 2002). Sulentic et al. (2001) suggest that the apparent spatial coincidence is a projection and that NGC 7318b, and thus the shock caused by it, are moving along the line of sight.

The presence of molecular gas indicates a rather cool environment. It is therefore surprising to find large amounts and extended CO emission at positions where shocks are present and are heating and ionizing the gas. J-shocks with velocities above about 50 km s-1 are able to effectively destroy molecular gas (Hollenbach & McKee 1980; 1989). Shocks with these velocities are plausibly present given the large velocity difference (600 – 1000 km s-1) between the intruder NGC 7318b and the rest of the SQ galaxies.

This strongly suggests that shocks do not affect the gas at 6000 km s-1, the velocity at which we found the major part of the CO emission. This means that the X–ray and radio continuum emission tracing shocks are only associated with the gas at 6700 km s-1. In addition to the presence of extended molecular gas, support for this hypothesis comes from the atomic gas. Similarly to the CO, the HI component at 6000 km s-1 is very extended, indicating that shock heating is not affecting this area so that cool gas can exist. Contrary to this, the HI component at 6700 km s-1 is spatially much more concentrated. Furthermore, the spatial coincidence between the 15 $`\mu `$m and and CO emission support the absence of shocks because dust is very easily destroyed by them (e.g. Borkowski & Dwek 1995, Jones et al. 1996).

The molecular gas mass in the starburst region111 I.e. adding the 4 pointings around the position (-5′′, 5′′), where the maximum of the H$`\alpha `$ emission is found, and summing both velocity components.  is $`7.0\pm 1.0\times 10^8`$ M. The surprising result of our observations is that the molecular gas content of the starburst region is only a small part (about one fourth) of the total molecular gas content and that the starburst does not mark the maximum of the CO emission. Our molecular gas mass in the starburst region is only slightly higher than the mass estimated from interferometric observations by Gao & Xu (2000) ($`2.45.3\times 10^8`$ M). The observations of Gao & Xu (2000) showed CO emission only at the position of the starburst. This could indicate that the extended emission that we have found is more smoothly distributed and therefore not detected by interferometric observations.

The origin of the molecular gas in SQ A is unclear. It could have been stripped together with the atomic gas from NGC 7318b or it could have formed in situ. In the center of NGC 7318b, a molecular gas mass of $`1.2\times 10^9`$ M has been found (Smith & Struck 2001). This yields, together with its blue luminosity of $`4.6\times 10^{10}`$ M, a ratio of $`L_\mathrm{B}/M_{\mathrm{H}_2}=38`$, considerably higher than the average value of $`L_\mathrm{B}/M_{\mathrm{H}_2}6`$ (with a scatter of a factor 3) that has been found for isolated galaxies of this luminosity (Perea et al. 1997). Thus, the molecular gas in SQ A could have its origin in NGC 7318b. On the other hand, if the molecular gas has formed in situ, a large scale process and not gravitational collapse must be responsible for its formation over a scale of about 25 kpc. The velocity of this gas, 6000 km s-1 (which is higher than the central velocity of NGC 7318b of 5765 km s-1) could be an indication that it has already been accelerated by the collision with the intergalactic gas and the slight offset between the molecular (6030 km s-1) and atomic gas (6000 km s-1) central velocities might show that the molecular gas is forming in the most accelerated parts.

### 4.2 Origin of the molecular gas in SQ B

The molecular gas in SQ B could have either been stripped from the parent galaxy together with the atomic gas or it could have formed in situ. Studying a sample of TDGs, Braine et al. (2001) argue in favor of the in-situ formation. Their main argument is that the CO distribution in spiral galaxies is normally concentrated towards their central parts whereas the atomic gas is much more widely distributed and therefore more easily stripped. In SQ the situation is more complicated because more than two galaxies are interacting. However, the arguments of Braine et al. (2001) seem to be applicable to NGC 7319, the galaxy thought to be the progenitor of the eastern HI cloud and the young optical tidal tail. All of the HI from this galaxy has been stripped, whereas CO is still found within its optical disk (Yun et al. 1997, Smith & Struck 2001). Its molecular gas content of $`4.8\times 10^9`$ M (Smith & Struck 2001) and its blue luminosity ($`L_\mathrm{B}=5.5\times 10^{10}`$ L; Verdes-Montenegro et al. 1998) yield $`L_\mathrm{B}/M_{\mathrm{H}_2}=11`$ which is normal compared to other galaxies (see above in Sect. 4.1). This supports the idea that in NGC 7319, as in other interacting galaxies, the HI was more easily stripped whereas CO stayed within the optical disk.

In SQ B the molecular gas is found in close association with the atomic gas and with indicators of star formation, as described in detail in the previous sections. Altogether, this spatial and (in the case of HI also kinematical) coincidence suggests that we are seeing in-situ molecular gas formation possibly provoked by compression of the atomic gas, and the subsequent star formation traced by H$`\alpha `$ emission and other tracers. Further support for this comes from the distribution of M$`_{\mathrm{H}_2}`$/MHI. If the molecular gas had been stripped together with the atomic gas, the molecular-to-atomic gas mass ratio would have been preserved, if no molecular gas formation or destruction took place after the stripping. In this case it would be an unlikely coincidence that we find the highest value of M$`_{\mathrm{H}_2}`$/MHI close to the HI peak.

### 4.3 Tidal Dwarf Galaxies in Stephan’s Quintet?

We observed very different regions in SQ and our CO observations have shown that the properties of the ISM are very diverse. Some of our observed regions have been claimed to be TDG candidates (e.g. Mendes de Oliveira et al. 2001, Iglesias-Páramo & Vílchez 2001). Duc et al. (2000) defined a TDG as a self-gravitating entity formed out of tidal debris. Although it is, without a detailed kinematical study, difficult to decide whether a given TDG candidate is indeed self-gravitating, Braine et al. (2000) have argued that the detection of molecular gas forming in situ gives strong support to the hypothesis of self-gravitation. Most objects observed in SQ are, however, different in one way or another from the sample of TDGs observed by Braine et al. (2000; 2001).

SQ B is the object that most closely resembles TDGs. Its has a clear tidal origin and we have summarized evidence that the molecular gas is forming in situ. The ratio of molecular gas mass to the H$`\alpha `$ luminosity is normal for star forming galaxies. The star formation rates that can be derived from different tracers such as H$`\alpha `$ emission, 15 $`\mu `$m emission (following the estimate in Xu et al. 1999) or radio continuum (applying eq. 21 in Condon (1992), extrapolated with a Salpeter initial mass function down to stars with masses above 0.1 M) give similar values of about 0.2 Myr-1. These values are in the same range as those found for TDGs ($`0.010.3`$ Myr-1) and several nearby dwarf galaxies ($`0.050.9`$ Myr-1) (Braine et al. 2001) showing that the observations of SQ B are consistent with the formation of a dwarf galaxy in the tidal tail. However, the large extent of the molecular gas covering the entire tip of the optical tail is surprising and different from other TDGs. One reason is probably the fact that star formation happens at two positions along the optical tidal tail: apart from SQ B there is a weaker star-forming knot showing a blue colour and H$`\alpha `$ emission at the very end of the optical tidal tail (at $`\alpha _{2000}=22^h36^m12.2^s`$ and $`\delta _{2000}=33^{}57^{}30`$′′). Molecular gas was detected at this position. Higher resolution observations are necessary in order to find out whether the molecular gas is concentrated towards the two star forming regions as would be expected if TDGs are forming.

The situation in SQ A is entirely different. Whereas the star formation is very concentrated at the central position of SQ A, the distribution of the molecular gas is much more extended with no maximum at the starburst. This lack of association of the molecular gas with the starburst region, together with the unusual origin of the ISM in this region brings up strong doubts whether SQ A is a TDG, i.e. an object that is bound by gravitation. These doubts are reinforced by the large dynamical mass that is deduced from the CO velocity gradient in SQ A (see Sect. 3.2). Under the assumption that the region around SQ A is gravitationally bound and that the velocity gradient is due to rotation we can calculate the dynamical mass ($`M_{\mathrm{dyn}}=R\mathrm{}V^2/G`$, where $`R`$ is the radius of the region, $`\mathrm{}V`$ is the velocity gradient along $`R`$ and $`G`$ is the gravitational constant). With $`R=30`$′′ and $`\mathrm{}V=80`$ km s-1 we derive $`M_{\mathrm{dyn}}=1.9\times 10^{10}`$ M. This is much larger than the mass of the molecular and atomic gas at this velocity ($`3.8\times 10^9`$ M). The stellar mass in the starburst region has been estimated by Xu et al. (1999) from a K-band image to be $`0.81.6\times 10^7`$ M, much smaller than the gas mass. Thus, the dynamical mass is much higher than the visible mass meaning that this region would need a high fraction of dark matter in order to be gravitationally bound. This is in contrast to what has been found for a sample of TDGs (Braine et al. 2001) where the visible and dynamical masses were found to be very close obviating the need for dark matter. We interpret the high dynamical mass as evidence that this region is not gravitationally bound on the size scales observable at the current resolution. The extended and homogeneous distribution of the molecular gas in this region makes it unlikely that a gravitational collapse is responsible for its formation. The star formation in the starburst region SQ A has rather been triggered by local processes, as e.g. by the pressure of the hot surrounding gas, as suggested by (Xu et al. 1999).

Most of the objects in SQ C possess a gradient in the H$`\alpha `$ emission indicative of rotation (Mendes de Oliveira et al. 2001) and might therefore be self-gravitating objects. However, the lack of gas, both atomic (except for C16/17 and C18/19 the atomic gas densities are $`10^{20}`$ cm-2) as well as molecular, makes them different from other TDGs. For this reason it is unlikely that they will become like TDGs that resemble dwarf irregulars or blue compact dwarfs as they do not have the gas reservoir to continue forming stars for a long time.

## 5 Conclusions

We have observed CO in several regions in SQ. Our main results are:

1. In both SQ A and SQ B we detected large amounts of molecular gas of $`3.1\times 10^9`$ M and $`7.0\times 10^8`$ M, respectively. The molecular-to-atomic gas mass ratios are high (1.2, respectively 0.5), much larger than what has been found in typical TDGs (Braine et al. 2001). The molecular gas is very extended, over regions of about 15 – 25 kpc.
2. The CO spectrum in SQ A consists of lines at different velocities, centered at 6030 and 6700 km s-1, with most of the emission coming from the low velocity component, whereas in SQ B the CO is at 6625 km s-1. The kinematical agreement (velocity, width and spectral shape) between the CO and HI lines is very good in both regions, except for a small offset in the low-velocity components of the CO and HI lines in SQ A (the HI line peaks at 6000 km s-1).
3. In SQ A, the CO emission at 6000 km s-1 is found south of the starburst region, covering the southern part of the HI distribution at this velocity and coinciding with X–ray and radio continuum emission which are indicative of shocks. The much weaker CO emission at 6700 km s-1 is more spatially concentrated. We conclude that the presence of large amounts of extended molecular gas at 6000 km s-1 implies that the shocks traced by X–ray and radio continuum emission only affect the gas at 6700 km s-1 and that the apparent spatial coincidence is a projection.
4. In SQ B there is good spatial agreement between the tracers of star formation (H$`\alpha `$, 15 $`\mu `$m and radio continuum) and the CO emission. The CO peak is slightly offset from the HI peak towards a steep HI gradient. We interpret this as indicating that the molecular gas is forming in-situ, possibly in a region of compressed HI and that the star formation is taking place subsequently. This, together with the fact that the emission at different wavelengths is similar to that of normal star forming galaxies and that SQ B has a clear tidal origin, makes it the object in SQ that most resembles a TDG.
5. We have searched without success for CO in several objects around NGC 7318a/b in which Mendes de Oliveira et al. (2001) found H$`\alpha `$ emission. The upper limits for the ratios between molecular gas mass and SFR (derived from the H$`\alpha `$ emission) are lower than those found in spiral galaxies and TDGs. The molecular-to-atomic gas mass ratio is also low, though not unusual. The low atomic and molecular gas content of these regions makes them different from TDGs and will not allow them to sustain star formation for a long time.

###### Acknowledgements.

We would like to thank J. Iglesias-Páramo and J. Vílchez for making available their H$`\alpha `$ data to us, L. Verdes-Montenegro for the HI data, C. Xu for the ISOCAM and H$`\alpha `$ data and S. Gallagher for useful information about the star clusters in SQ. Thanks are also due to the referee, G. Petitpas, for the careful reading of the manuscript and detailed suggestions for improvements. VC would like to acknowledge the partial support of JPL contract 960803 and EB acknowledges financial support from CONACyT via project 27607-E. We made use of data from the Canadian Astronomy Data Center, which is operated by the Dominion Astrophysical Observatory for the National Research Council of Canada’s Herzberg Institute of Astrophysics.


format: text

text:
# 1 Introduction

## 1 Introduction

The BNL $`g_\mu 2`$ Collaboration has announced an improved result for $`a_\mu =(g_\mu 2)/2`$ with twice the precision of their 2001 analysis. Here we interpret these results in the context of supersymmetry. It has been recognized for some time that the muon anomalous moment can act as an important probe of physics beyond the standard model especially of supersymmetry. This is so because the new physics contribution to the leptonic anomalous magnetic moment $`a_l`$ scales as $`m_l^2/\mathrm{\Lambda }^2`$ and is proportional to the square of the lepton mass. Thus $`a_\mu `$ provides a more sensitive probe of new physics than $`a_e`$ even though $`a_e`$ is more sensitively determined while the $`a_\tau `$ determination is far less sensitive to be competitive. Because of the above any improvement in the determination of $`a_\mu `$ has important implications for new physics beyond the standard model. Thus the improved measurement of BNL experiment in 2001 which gave $`a_\mu ^{exp}=11659203(15)\times 10^{10}`$ resulted in a large amount of theoretical activity to explore the implications of the results. At the time of the BNL2001 result the standard model prediction consisting of the qed, electro-weak and hadronic corrections was estimated to be $`a_\mu ^{SM}=11659159.7(6.7)\times 10^{10}`$ which gave $`a_\mu ^{exp}a_\mu ^{SM}=43(16)\times 10^{10}`$ implying a $`2.6\sigma `$ difference between experiment and theory. The above difference was based on a light by light (LbL) hadronic correction of  $`a_\mu ^{had}(LbL)=8.5(2.5)\times 10^{10}`$ which was later found to be in error. Thus the analysis of Knecht et.al. gives  $`a_\mu ^{had}(LbL)=8.3(1.2)\times 10^{10}`$ while the revised analysis of Hayakawa and Kinoshita gives $`a_\mu ^{had}(LbL)=8.9(1.5)\times 10^{10}`$. A partial analysis of $`a_\mu ^{had}(LbL)`$ has also been given in Refs.. Thus the partial analysis of Ref. finds $`a_\mu ^{had}(LbL)=8.3(3.2)\times 10^{10}`$ and the work of Ref. which computed the pion pole part finds $`a_\mu ^{had}(LbL:\pi ^0pole)=5.6\times 10^{10}`$. Again the sign of these corrections agree with the sign of the reevaluation of this quantity in Refs.. Using the average of the first two which are the more complete calculations and using the hadronic vacuum polarization correction of Ref. one finds $`a_\mu ^{SM}=11659176.8(6.7)\times 10^{10}`$ and the BNL2001 result gives $`a_\mu ^{exp}a_\mu ^{SM}=26(16)\times 10^{10}`$. This difference corresponds only to a $`1.6\sigma `$ deviation between theory and experiment.

Before proceeding further we wish to discuss a bit further the issue of errors in the hadronic corrections specifically the LbL correction and the $`\alpha ^2`$ vacuum polarization correction. Regarding the LbL hadronic correction, in addition to the analyses mentioned above, i.e., Refs. there is also the analysis of Ref. based on chiral perturbation theory. This analysis finds $`a_\mu ^{had}(LbL)=(1.3_6^{+5}+3.1\stackrel{~}{C})\times 10^{10}`$ where $`\stackrel{~}{C}`$ is a correction arising from the sub leading contributions. These contributions are either of $`O(\alpha ^3p^2/\mathrm{\Lambda }^2)`$, where p is a mass of order $`m_\mu `$ or $`m_\pi `$ and $`\mathrm{\Lambda }`$ is a hadronic scale $`1`$ GeV, which are not enhanced by a factor of $`N_C`$ (the number of quark colors) or of $`O(N_C\alpha ^3p^2/\mathrm{\Lambda }^2)`$ but are not enhanced by large logarithms. The coefficients $`\stackrel{~}{C}`$ in this analysis is an unknown low energy constant (LEC) which contains the contributions of nonperturbative physics at short distance. The authors of Ref. view the parameter $`\stackrel{~}{C}`$ as basically unconstrained except through the measurement of the anomalous moment itself. In our analysis we assume the validity of the analyses of Refs.. Aside from the LbL contribution, the other source of error in the hadronic correction is the contribution from the $`\alpha ^2`$ hadronic vacuum polarization correction and this error is an important component in extracting the deviation between experiment and the standard model.

The new BNL result gives the world average on $`a_\mu `$ as .

$$a_\mu ^{exp}=(11659203)(8)\times 10^{10}$$
(1)
With the LbL correction of Refs. and assuming the leading hadronic correction to lie in the range $`692(6)\times 10^{10}`$ to $`702(8)\times 10^{10}`$ the standard model prediction lies in the range $`a_\mu ^{SM}=11659177(7)\times 10^{10}`$ to $`a_\mu ^{SM}=11659186(8)\times 10^{10}`$ which leads to

$`a_\mu ^{exp}a_\mu ^{SM}=(26)(10)\times 10^{10}`$
$`to(17)(11)\times 10^{10}`$ (2)
which corresponds to a difference between experiment and theory of about $`1.6\sigma `$ to $`2.6\sigma `$.

## 2 Interpreting the New BNL Result for Supersymmetry

Next we discuss the implications of these results for supersymmetry. The supersymmetric electroweak contribution to $`a_\mu `$ arises from the chargino and neutralino exchange corrections. For the CP conserving case the chargino ($`\chi ^\pm `$) and the sneutrino ($`\stackrel{~}{\nu }`$) exchange contributions are typically the larger contributions and here one finds

$$a_\mu ^{\stackrel{~}{\chi }^\pm }=\underset{a=1,2}{}(\frac{m_\mu ^2}{48\pi ^2}\frac{A_{R}^{(a)}{}_{}{}^{2}}{m_{\stackrel{~}{\chi }_a^\pm }^2}F_1(\left(\frac{m_{\stackrel{~}{\nu }}}{m_{\stackrel{~}{\chi }_a^\pm }}\right)^2)+\frac{m_\mu }{8\pi ^2}\frac{A_R^{(a)}A_L^{(a)}}{m_{\stackrel{~}{\chi }_a^\pm }}F_2(\left(\frac{m_{\stackrel{~}{\nu }}}{m_{\stackrel{~}{\chi }_a^\pm }}\right)^2))$$
(3)
where $`A_L(A_R)`$ are the left(right) chiral amplitudes and $`F_1,F_2`$ are form factors as defined in Ref.. The accuracy of the result of Ref. was further tested by taking the supersymmetric limit of the results of Ref. in Ref. where in addition the effects of CP violation on $`g_\mu 2`$ were also investigated. The amplitude of Eq.(3) is dominated by the chiral interference term proportional to $`A_LA_R`$. Since $`A_L1/\mathrm{cos}\beta `$ ($`\mathrm{tan}\beta =<H_2>/<H_1>`$ where $`<H_2>`$ gives mass to the up quark and $`<H_1>`$ gives mass to the down quark and the lepton) one finds that $`a_\mu ^{SUSY}\mathrm{tan}\beta `$ for large $`\mathrm{tan}\beta `$. Further, one finds that $`A_L`$ depends on the sign of $`\mu \stackrel{~}{m}_2`$, where $`\stackrel{~}{m}_2`$ is $`SU(2)`$ gaugino mass and $`\mu `$ is the higgs mixing parameter both taken at the electroweak scale, and thus the sign of $`a_\mu ^{SUSY}`$ is determined by the sign of $`\mu \stackrel{~}{m}_2`$. The above leads one to the result that over a large part of the parameter space unless $`\mathrm{tan}\beta 1`$ one finds

$$a_\mu ^{SUSY}>0,\stackrel{~}{m}_2\mu >0;a_\mu ^{SUSY}<0,\stackrel{~}{m}_2\mu <0$$
(4)
where we have used the sign convention on $`\mu `$ of Ref.. These results have been confirmed numerically in a wide class of models. The implications of the BNL result of 2001 has been analyzed extensively in the literature within the framework of mSUGRA as well as in a variety of other models.

We analyze now the implications of the new result within mSUGRA. We note in passing that for a class of models with large extra dimensions the correction to g-2 from Kaluza-Klein modes is rather small once the constraints arising from Kaluza-Klein corrections to $`G_F`$ are taken account of. Thus the extra dimensions do not interfere in the extraction of implications of the new result for supersymmetry. In our analysis we will use the $`2.6\sigma `$ deviation in Eq.(2) as the default value for the deviation but we will take a $`2\sigma `$ error corridor in the analysis. This error corridor includes the case corresponding to the smaller value of the deviation in Eq.(2), i.e., $`1.6\sigma `$ and taking a $`1\sigma `$ error corridor in the analysis. Our assumption thus gives us the constraint

$$6\times 10^{10}(a_\mu ^{exp}a_\mu ^{SM})46\times 10^{10}$$
(5)
Attributing the entire difference $`a_\mu ^{exp}a_\mu ^{SM}`$ to supersymmetry and using Eqs.(2) and (4) we find that the sign of $`\mu `$ is positive confirming a similar result arrived at in analyses based on the BNL2001 data. The sign of $`\mu `$ is of great importance for dark matter and for $`b\tau `$ unification. As is well known $`bs+\gamma `$ imposes an important constraint on the parameter space of supersymmetric models. The constraint is very stringent for negative $`\mu `$ eliminating most of the parameter space of models. For positive $`\mu `$ the constraint is not very strong for small $`\mathrm{tan}\beta `$ but becomes a strong constraint for large values of $`\mathrm{tan}\beta `$. The imposition of this constraint involves the standard model value and the experimental value both of which have significant errors. The standard model prediction of this decay including the next to the leading order correction has been given by several authors. Further, there are several recent experimental determinations of this decay. In our analysis we take a $`2\sigma `$ error corridor between experiment and the prediction of the standard model and thus allow the decay branching ratio to lie in the range $`2\times 10^4<B(bs+\gamma )<4.5\times 10^4`$. Another important aspect of mSUGRA is that it leads to the lightest neutralino to be the lightest supersymmetric particle (LSP) using the renormalization group analysis (see, e.g., Ref.) and hence with R parity invariance a candidate for cold dark matter. The current astrophysical data allows the amount of cold dark matter to lie in the range $`0.1<\mathrm{\Omega }_{CDM}h^2<0.3`$ where $`\mathrm{\Omega }_{CDM}=\rho _{CDM}/\rho _c`$ where $`\rho _{CDM}`$ is the matter density due to cold dark matter and $`\rho _c`$ is the critical matter density needed to close the universe, and $`h`$ is the Hubble parameter in units of 100km/sMpc. In our analysis we assume the entire cold dark matter as arising from the lightest neutralino and thus impose the constraint on the neutralino relic density so that $`0.1<\mathrm{\Omega }_{\chi ^0}h^2<0.3`$.

The result of the analysis is exhibited in Fig.1 in the $`m_0M_{\frac{1}{2}}`$ plane, where $`m_0`$ is the universal scalar mass and $`M_{\frac{1}{2}}`$ is the universal gaugino mass, for $`\mu `$ positive and for values of $`\mathrm{tan}\beta `$ of 5,10,30 and 50. Thus yellow region-I at the top right hand side is disallowed either due to the absence of EWSB or due to the lighter chargino mass lying below the current experimental lower limit. Further, yellow region-II at the bottom is disallowed due to stau becoming the LSP. The red area on the boundary between the yellow and the white region is the zone satisfying the relic density constraint $`0.1<\mathrm{\Omega }_{\chi ^0}h^2<0.3`$ (the blue filled circles are the additional allowed points if the relic density constraint includes the lower region $`0.02<\mathrm{\Omega }_{\chi ^0}h^2<0.1`$). As discussed above $`bs+\gamma `$ is not a strong constraint for positive $`\mu `$ and small $`\mathrm{tan}\beta `$. As a consequence the $`bs+\gamma `$ constraint does not eliminate any parameter space for the cases $`\mathrm{tan}\beta =5`$ and $`\mathrm{tan}\beta =10`$ as can be seen from Figs.(1a) and (1b). However, for large values of $`\mathrm{tan}\beta `$ the $`bs+\gamma `$ constraint becomes important even for positive $`\mu `$ as can be seen from Figs.(1c) and (1d). We turn now to the imposition of the $`a_\mu ^{exp}a_\mu ^{SM}`$ constraint of Eq.(2). In the theoretical computation of the supersymmetric electroweak correction we take into account the full one loop supersymmetric contribution including the chargino-sneutrino and the neutralino-smuon exchanges. The result of imposing $`2\sigma `$ variation around the central value gives the corridor shown in Figs.1 by the upper ($`2\sigma `$) and the lower ($`+2\sigma `$) solid lines. Specifically the ranges of sparticle masses, Higgs boson masses and the $`\mu `$ parameter corresponding to $`1\sigma `$ and $`2\sigma `$ constraints are exhibited in Fig.2. The spectra span a wide range and with the lower end of the scale corresponding to a low values of fine tuning and the high end corresponding to high values of fine tuning (see, e.g., Ref.). The spectrum of Fig.2 leads us to the conclusion that sparticles satisfying $`1\sigma `$ constraint are well within reach of the LHC and some of the sparticles could also show up at RUNII of the Tevatron. The ranges allowed by the $`2\sigma `$ constraint constraint have upper limits which in some cases reach as high as 3 TeV and may exceed the reach of the LHC. However, Fig. 2 does not contain the constraints of $`bs+\gamma `$ or of relic density. With these additional constraints the upper limits would be substantially smaller.

We discuss now the direct detection of dark matter under the BNL2002 $`g_\mu 2`$ constraint. The analysis of Fig.1 shows that there are substantial regions of the parameter space where all the constraints, i.e., the g-2 constraint the $`bs+\gamma `$ constraint and the relic density constraints are simultaneously satisfied. We discuss now the direct detection of dark matter under these constraints. There has been considerable work on the theoretical predictions of the neutralino-proton cross section $`\sigma _{\chi p}`$ which enters in the direct detection of dark matter and some recent literature can be found in Ref.. Here we give the analysis of these cross sections in mSUGRA for positive $`\mu `$ under the new $`g_\mu 2`$ constraint. Results are exhibited in Figure 3 for values of $`\mathrm{tan}\beta `$ of 5,10,30 and 50. We find that the neutralino-proton cross sections consistent with all the constraints lie in the range $`10^{46}cm^2\sigma _{\chi p}2\times 10^{43}cm^2`$. A significant part of this range will be accessible at future dark matter detectors GENIUS and ZEPLIN.

## 3 Conclusion

The new Brookhaven measurement of $`g_\mu 2`$ has confirmed the earlier measurement with twice the precision. However, interpretation of what this implies for physics beyond the standard model is very sensitive to the errors in the standard model prediction. The main surprise at the end of last year was the switch in sign of the LbL contribution. There is now a general agreement that the sign of this contribution is positive and also an agreement between two independent evaluations on its size. An important exception to this is the result of chiral perturbation theory, which although a more fundamental approach is beset by the appearance of unknown low energy constants which require a nonperturbative approach such as lattice gauge calculation for its evaluation. In the present analysis we have assumed the validity of the analysis on LbL of Knecht and Nyffeler and the revised analysis of Hayakawa and Kinoshita. Taking account of all the errors the difference between experiment and the standard model is now predicted to lie between $`1.6\sigma `$ and $`2.6\sigma `$. We have revisited the implications of the $`g_\mu 2`$ constraint for supersymmetry in view of the new result. Using the $`2.6\sigma `$ difference and a $`2\sigma `$ error corridor our analysis points to significant regions in the mSUGRA parameter space consistent with the $`bs+\gamma `$, relic density and $`g_\mu 2`$ constraints. Within $`1\sigma `$ error corridor all of the sparticles are accessible at the Large Hadron collider. Further, for low values of $`\mathrm{tan}\beta `$ sparticles may also be accessible at Fermilab Tevatron. We also carried out an analysis of the neutralino-proton cross section $`\sigma _{\chi ^0p}`$ and find that a significant part of the parameter space would be explored by the future dark matter detectors such as CDMS (Soudan), GENIUS and ZEPLIN. Finally, the $`g_\mu 2`$ experiment could lead to an even more stringent constraint if there was a reduction in the error associated with the standard model prediction. This could come about by eliminating the ambiguity on the LbL correction through a lattice gauge calculation analysis, and through improved low energy data on $`e^+e^{}hadrons`$ in the analysis of hadronic vacuum polarization correction.

We thank William Marciano for an interesting discussion. This work is supported in part by NSF grant PHY-9901057.

Figure Captions
Figure 1: Analysis in mSUGRA exhibiting the allowed and the disallowed regions in the $`m_0M_{\frac{1}{2}}`$ plane due to the radiative breaking of the electro-weak symmetry, the $`bs+\gamma `$ constraint and the neutralino relic density constaint as discussed in the text for values of $`\mathrm{tan}\beta `$ of 5,10,30 and 50. The area to the right of the dotted blue line in each case is the allowed region by the $`bs+\gamma `$ constraint and the (red) dark area satisfy the relic density constraint $`0.1<\mathrm{\Omega }_{\chi ^0}h^2<0.3`$ (while the (blue) dark filled circles satisfy $`0.02<\mathrm{\Omega }_{\chi ^0}h^2<0.1`$). The region I is discarded because of the absence of the electroweak radiative symmetry breaking or via the lower limit of chargino mass constraint. The region II is eliminated because of stau becoming the LSP. The solid black lines exhibit the lower limit ($`a_\mu (2\sigma )`$) and the upper ($`a_\mu (2\sigma )`$) limit corresponding to the $`\pm 2\sigma `$ variation around the central value given by the first line of Eq.(2). The LEP ’higgs signal’ is also indicated in the figures.

Figure 2: Mass ranges of sparticles and Higgs bosons as well as the $`\mu `$ parameter in mSUGRA for $`\mu >0`$ when the $`1\sigma `$ and $`2\sigma `$ $`g_\mu 2`$ constraints are imposed and $`\mathrm{tan}\beta `$ ranges from 2-60 and $`A_0=0`$. The $`bs+\gamma `$ constraint and the relic density constraint are not imposed. For each (s)particle the vertical line to the left is the range in mass allowed by a $`2\sigma `$ $`g_\mu 2`$ constraint while the vertical line to the right is the mass range allowed by a $`1\sigma `$ $`g_\mu 2`$ constraint.

Figure 3: A plot of the scalar proton-LSP cross-section $`\sigma _{\chi p}`$ vs LSP mass $`m_\chi `$ which enters in the direct detection rates. The blue dots are the points allowed in mSUGRA and the red filled circles satisfy the constraint $`0.1<\mathrm{\Omega }h^2<0.3`$. Black squares which form the black region satisfy the $`2\sigma `$ constraints on $`bs+\gamma `$ and $`g_\mu 2`$. The enclosed region on the top left hand side is the region where the DAMA collaboration claims a signal. The dot-dashed line on top is the upper limit from CDMS while the solid line is the latest upper limit from EDELWEISS experiment. The lower dashed line is the sensitivity that will be reached by the CDMS experiment at Soudan mine while the dotted line is the sensitivity that will be achieved by the GENIUS detector.