Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lav01.sinp.msu.ru/~vlk/etapiz_draft9a_PL.pdf Äàòà èçìåíåíèÿ: Fri Apr 4 13:55:54 2008 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 20:09:47 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 78

Confirmation of the 1-+ Meson Exotics in the 0 System
G. S. Adams,1 T. Adams,2 , Z. Bar-Yam,3 J. M. Bishop,2 V. A. Bodyagin,4 , D. S. Brown,5 , N. M. Cason,2 S. U. Chung,6 J. P. Cummings,1 A. I. Demianov,4 K. Danyo,6 J. P. Dowd,3 P. Eugenio,7 X. L. Fan,5 A. M. Gribushin,4 R. W. Hackenburg,6 M. Hayek,3 , § J. Hu,1 , ¶ E. I. Ivanov,8 D. Joffe,5 W. Kern,3 E. King,3 O. L. Kodolova,4 V. L. Korotkikh,4 M. A. Kostin,4 J. Kuhn,1 J. M. LoSecco,2 M. Lu,1 L. V. Malinina,4 J. J. Manak, M. Nozar,1 , C. Olchanski,6 , ¶ A. I. Ostrovidov,7 T. K. Pedlar,5 , L. I. Sarycheva,4 K. K. Seth,5 N. Shenhav,3 , § X. Shen,5, 9 , W. D. Shephard,2 N. B. Sinev,4 D. L. Stienike,2 J. S. Suh,6 , §§ S. A. Taegar,2 A. Tomaradze,5 I. N. Vardanyan,4 D. P. Weygand,9 D. B. White,1 H. J. Willutzki,6 , M. Witkowski,1 and A. A. Yershov4 (The E852 collaboration)
3

2

Department of Physics, Rensselaer Polytechnic Institute, Troy, New York 12180 Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556 Department of Physics, University of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747 4 Nuclear Physics Institute, Moscow State University, Moscow, Russian Federation 119899 5 Department of Physics, Northwestern University, Evanston, Il linois 60208 6 Physics Department, Brookhaven National Laboratory, Upton, New York 11973 7 Department of Physics, Florida State University, Tal lahassee, FL 32306 8 Department of Physics, Idaho State University, Pocatel lo, ID 83209 9 Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 (Dated: June 21, 2007)
2

1

The exclusive reaction - p 0 n, + - 0 at 18 GeV/c has been studied with a partial wave analysis on a sample of 23 492 0 n events from BNL experiment E852. A mass-dependent fit is consistent with a resonant hypothesis for the P+ wave, thus providing evidence for a neutral exotic meson with J P C = 1-+ , a mass of 1257 ± 20 ± 25 MeV/c2 , and a width of 354 ± 64 ± 60 MeV/c2 . New interpretations of the meson exotics in neutral 0 system observed in E852 and Crystal Barrel experiments are discussed.
PACS numbers: 13.20.Jx; 13.85-t; 14.40.Cs

Keywords: meson spectroscopy; exotic mesons
INTRODUCTION.

Exotic mesons with J P C = 0-- , 1-+ , 2+- , . . . do not mix with quark-antiquark mesons and thus offer a natural testing ground for QCD. Exotic mesons have been discussed [1­10] for many years but have only recently been observed experimentally. The underlying structure of the negatively charged exotic state with J P C = 1-+ observed in this experiment [11, 12] at 1400 MeV decaying into - is not yet understood. Two distinguishing characteristics of the 0 system make it an excellent one to clarify the properties of this exotic state. First, C -parity is a good quantum number in the 0 system, unlike the - system. Second, the production mechanism for the charge exchange reaction - p 0 n cannot involve the exchange of an isospin I = 0 system such as the pomeron. The Crystal Barrel experiment [13] confirmed the existence of resonant structure in the - system using stopped antiprotons in liquid deuterium in the reaction pn - 0 . Later this group analyzed data on pp annihilation at rest ¯ ¯ into 0 0 [14] and presented evidence for an exotic 1-+ resonance in the 0 system with M = (1360 ± 25) MeV/c2 and = (220 ± 90) MeV/c2 . The 0 state has been studied in the GAMS experiment [15] with - p 0 n, 2 , 0 2 at 32, 38 and 100 GeV/c. The statistics of the 38 GeV/c data was sufficient so that, using the method of Sadovsky [16], they were able to present evidence for the exotic 1 (1400). The VES experiment also observed a peak in the P+ wave of the 0 system near 1400 MeV/c2 [17]. In their most recent publication [18], using theoretical arguments the authors state that the peak can be understood without requiring an exotic meson. An analysis of E852 data with - p 0 p, 2 was recently reported [19]. A bump in the P+ wave of the 0 system was observed at M ( 0 ) = 1272 MeV/c2 with a large width ( = 660 MeV/c2 ) when fitting all the data using the method of Ref. [12]. For small t , they observe a width of 190 MeV/c2 . The authors chose not to claim evidence for exotic 1 (1400) meson production.


2 The present analysis studies - p 0 n, + - 0 at 18 GeV/c in E852. The advantage of this mode over the all-neutral final state is that the production vertex is defined by charged tracks. This improves the mass resolution as well as the ability to require that the interaction took place within the target.
EXPERIMENTAL SETUP AND DATA SELECTION

The data for this analysis were obtained at the Alternating Gradient Synchrotron (BNL USA) in 1995. Using an 18 GeV/c - beam interacting in a liquid hydrogen target, a total of 750 million triggers were acquired of which 108 million were of a type designed to enrich the exclusive final state - p + - 4 n. A total of 6 million events of this type were fully reconstructed. The data were kinematically fit [20] to select events consistent with the - + 0 0 n hypothesis (with a confidence level of at least 0.01%) yielding about 4 million events. Of those, 85 228 events passed a mass cut enhancing mesons, m( - + 0 ) < 0.65 GeV/c2 , and 74 549 passed a cut to remove events passing through a low-efficiency region in the drift chambers. A final kinematic fit selected 23 492 events for the partial wave analysis (PWA), which were consistent with the 0 n, + - 0 hypothesis at a minimum confidence level of 1%. A strong meson signal is observed in this final data sample (Fig. 1a) with a mass of 539.2 ± 0.3 MeV/c2 and a (resolution-dominated) width of 23.7 ± 0.2 MeV/c2 . The filled regions in the figure indicate the signal region and the side-band regions used in the analysis. In the signal region, the signal-to-background ratio is about 4 to 1 for all 0 masses and 5 to 1 for m( 0 ) > 1.1 GeV/c2 . The 0 mass spectrum shown in Fig.1b has two clear peaks: the a0 (980) and the a0 (1320). 0 2 The non- background was estimated as a function of 0 mass using the side-band and signal regions. The background fraction varies between 24% and 14% going from lower to higher mass in the region 0.78 < m( 0 ) < 1.74 GeV/c2 . The experimental acceptance was determined using a Monte Carlo event sample generated with isotropic angular distributions in the Gottfried-Jackson frame. The detector simulation was based on the E852 detector simulation package SAGEN [12]. The experimental acceptance was incorporated into the PWA by means of Monte Carlo normalization integrals [12]. The acceptance as a function of mass and as a function of t is flat.
PARTIAL WAVE ANALYSIS

The partial-wave analysis (PWA) method described in [12] (see also [21, 22]) was used to study the spin-parity structure of the 0 system in this data set. The PWA was carried out using the extended maximum likelihood method separately in each mass bin in the mass region between 0.78 and 1.74 GeV in mass bins of 0.04 GeV for 0 < |t | < 1.0 (GeV/c)2 using the likelihood function
n

lnL
i

lnI (i ) -

d () I ().

(1)

Here I () is the predicted angular distribution, () is the angular acceptance, and the sum is over the event sample. The partial waves are parameterized by a set of five numbers: J P C m , where J is the angular momentum, P and C are the parity and the C-parity of the 0 system, m is the absolute value of the angular momentum pro jection and is the reflectivity. We use a simplified notation where each partial wave is denoted by a letter indicating the 0 system's angular momentum in standard spectroscopic notation, and a subscript which can take the values 0, +, or -, for m = 0- , 1+ , or 1- respectively. We assume that the contribution from partial waves with m > 1 is small and can be neglected [12]. The amplitudes used are the unnatural parity-exchange waves (UNPW) S0 , P0 , P- , D0 , D- , and the natural parity-exchange waves (NPW) P+ , D+ . The NPW waves interfere between themselves as do the UNPW waves but the NPW waves do not interfere with the UNPW waves. The P+ wave would be an exotic J P C = 1-+ (denoted by 0 1 ) if the wave is resonant. For each partial wave the complex production amplitudes were determined from an extended maximum likelihood fit [22]. The spin 1/2 nature of the target proton leads to spin-flip and spin-non-flip amplitudes and thus to a production spin-density matrix with maximal rank two. The PWA fit presented in this paper was carried out with the assumption that a spin-density matrix of rank one was sufficient [12]. An isotropic incoherent background was included. The magnitude of the background was fixed as determined from the side bands. We investigate the quality of the fits by comparing the moments of the decay angular distributions H (LM ), L 4 [12, 22], of the data with those predicted


3

1000 500

Events/0.002 GeV

400

a)
Events/0.02 GeV

800

b)

600

300

200

400

100

200

0 0.45

0 0.5 0.55
0+-

0.6
2

0.65

0.8

1

1.2

1.4

1.6

1.8
2

2

2.2

2.4

M( ),GeV/c

M( 0 ),GeV/c

FIG. 1: (a) Fit of the + - 0 mass distribution in the mass region. There are two entries per event: one for each way to assign a 0 to the decay. The signal region and the side-band regions are shown shaded. (b) The uncorrected 0 effective-mass distribution for events consistent with the reaction - p 0 n .

by Monte Carlo events generated with the fit amplitudes. We also directly compare the angular distributions for cos(GJ ) and T J between the data and those Monte Carlo events. The quality of the fits is good. Since natural-parity exchange (NPE) and unnatural-parity exchange (UNPE) amplitudes have different |t | dependences, a fit to the |t | distribution using a function of the form N (t ) = n1 |t |e-b1 |t | + n2 e-b2 |t | was carried out to determine the relative contributions of the two exchanges. The fitted parameters are b1 = (7.41 ± 0.08)(GeV/c)2 , b2 = (2.68 ± 0.07)(GeV/c)2 , and thus the ratio of UNPE and NPE contributions is equal to 0.71 ± 0.03. A value of about 70% for the ratio of UNPE to NPE at 18 GeV/c is expected based on interpolation of experimental - p charge exchange data between 3 and 40 GeV/c in the Regge model [16]. There is mathematical ambiguity in the description of a system of two pseudo-scalar mesons [23]. For our set of amplitudes there are eight ambiguous solutions, each of which leads to identical angular distributions. These solutions were found analytically starting from one solution by means of the Barrelet zeros method [22]. The eight solutions in each mass bin are shown in Fig. 2 as a point for every ambiguous solution. (In some cases the solutions are too close together to be visible as separate solutions. This is particularly true in the mass bins at M = 1.56 GeV/c2 and at M = 1.72 GeV/c2 .) The spread between the various ambiguous solutions dominates the systematic uncertainty in the resonance parameter determination (see below).
MASS DEPENDENT FITS

To study the resonant structure in the partial waves, we used three different procedures. The first method utilizes a Mass Dependent Fit (MDF) of the average solutions in the NPW sector. The PWA results in each mass bin were averaged between ambiguous solutions [12]. The mass dependence of the P+ and D+ intensities as well as their relative phase difference were then fit by relativistic Breit-Wigner (BW) functions (in both the P+ and D+ waves) with mass-dependent widths and Blatt-Weisskopf barrier factors[12]. The mass and width of the a0 are well known 2 and were fixed using values[24]: M = 1320 MeV/c2 , and = 120 MeV/c2 . The width of the a2 (1320) includes the experimental mass resolution. There are three free parameters in the fit of the D+ intensity, |D+ |2 : one for the magnitude and two parameterizing the smooth background for the D+ wave, as was done in [12]. In the MDF of the |P+ |2 distribution and the relative phase (D+ - P+ ) there are four free parameters: three from the BW function and one for the production phase (assumed constant). The fit was carried out in the mass interval 1.1 - 1.74 GeV/c2 . The resonant hypothesis for D+ and P+ waves with a mass-independent production phase gives a 2 /DoF=1.14 for 28 degrees of freedom. The non-resonant hypothesis (no phase variation for the P+ wave) gives 2 /DoF=3.02. It is clear from Fig. 2c that a single resonant phase for the a2 (1320) (dotted line) with a constant (non-resonant) P+ wave is not satisfactory. The contrast between the D+ - P+ phase variation and the D+ phase variation with mass (Fig. 2c, solid line) shows clearly that the observed phase variation is consistent with interference between two resonant waves. It's worth pointing out that the mass dependence of the D+ - P+ relative phase would be nearly flat if each wave were resonant with very similar masses and widths. The P+ resonant parameters from the fit with the average solutions and the average error matrix [12] are: M =


4
5000 4000 3000 2000 1000 0
2

)

Events / (0.04 GeV/c

a

|D +|

2

0.8

1

1.2

1.4

1.6

2

2000 1800 1600 1400 1200 1000 800 600 400 200 0

)

Events / (0.04 GeV/c

b

|P +|

2

0.8

1

1.2

1.4

1.6

0 Phase Di erence (rad) -0.5 -1 -1.5 -2 -2.5 -3 0.8

c
( D + - P +)
1 1.2 1.4 2 1.6

M( ) (GeV/c )

0

FIG. 2: The Partial Wave Analysis (PWA) and Mass-Dependent Fit (MDF) results. The points shown in each mass bin are the eight ambiguous PWA solutions. a) the D+ wave intensity; b) the P+ wave intensity; and c) the relative phase between the P+ and D+ waves. The lines show the MDF results (method 2). The average PWA solution in each mass bin is plotted using grey (yellow) points. The dotted line in (c) is the phase difference if the P+ phase is constant.

1265 ± 21 MeV/c2 and = 411 ± 64 MeV/c2 . The second method was similar to the first except that instead of fitting the average solutions, a large number ( 103 ) of randomly chosen combinations of ambiguous solutions in each mass bin were used as input to the massdependent fit. The obtained distributions of the mass and width of the P+ resonance for those fits with acceptable values of 2 (that is, those with 2 /DoF < 2) were then fitted by a Gaussian. The mean values of these distributions are: M = 1257 ± 20 MeV/c2 and = 354 ± 64 MeV/c2 . The curves shown in Fig. 2 are drawn using these mean values. The RMS values of these distributions are M = 25 MeV/c2 and = 58 MeV/c2 . The systematic errors are obtained using these RMS values. Our analysis shows that this spread due to the different ambiguous solutions dominates the systematic error. The third method used was to carry out a Mass-Dependent Partial Wave Analysis (MDPWA) [12]. In this procedure, an extended maximum likelihood function is generalized to include not only the angular distribution, but also the 0 mass distribution for each wave. This analysis is free from the problem of ambiguous solutions but it is necessary to parameterize the mass dependence of every partial wave (including the UNPWs) and all relative phases. We use the same parametrization for the D+ and P+ waves as in the first two methods. The mass dependence of the UNPW waves were chosen to be polynomials of second order with constant phases except for the S0 wave. The S0 wave was fitted with a BW function using the a0 (980) resonance parameters. The MDPWA results for the P+ wave are M = 1256 ± 10 MeV/c2 and = 319 ± 34 MeV/c2 , consistent with the results of methods 1 and 2. In [12] it was shown that a pure D+ wave can artificially induce a P+ wave due to the effects of finite acceptance and resolution. This "leakage" leads to a P+ wave that mimics the D+ intensity and phase. In our case, the P+ intensity would therefore have an intensity with the shape of the a2 (1320) and a D+ - P+ phase difference which doesn't depend on mass. These features allowed us to include in the MDPWA fit a leakage term with these features. We observed that the leakage contribution to the P+ wave from the D+ wave is negligible. Evidence for a resonance interpretation for the P+ wave is primarily the behavior of the D+ - P+ relative phase (Fig. 2c). Since the D+ phase variation is well known because of the a2 (1320) production, it is clear that the P+ phase cannot be constant (see dotted curve in Fig. 2c) and it is well-described by a BW phase variation. The ratio of the P+ and D+ intensities in the range 1.24 < M ( 0 ) < 1.34 GeV is equal to |P+ |2 /|D+ |2 = 0.43 ± 0.10.


5 This ratio is larger than that for the - system, as reported in Ref. [12]. The mass of the neutral exotic 1-+ state, decaying into 0 , observed here (M = 1257 ± 20 ± 25 MeV/c2 ) is lower than the mass observed in the Crystal Barrel experiment [14] (M = 1360 ± 25 MeV/c2 ) by about 100 MeV although the results are barely consistent within errors. The width measured here ( = 354 ± 64 ± 58 MeV/c2 ) is also consistent with that from the Crystal Barrel measurement ( = 220 ± 90 MeV/c2 ). It should also be noted that our result is similar to those obtained from the low-t fits in Ref. [19]. They obtained M = 1301±14 MeV/c2 and = 190±32 MeV/c2 in one analysis, and M = 1386±32 MeV/c2 and = 363±81 MeV/c2 in a fit to the experimental moments. No systematic errors were given for either method. The lower mass found in our analysis may be a consequence of interference between the resonant state and background in the 0 system, some of which may be from rescattering between the and the 0 . It is also possible that two or more resonant 1-+ states may be present in the 0 decay channel in the mass interval mass between 1200 and 1400 MeV as might be expected if the exotic state is a four quark state. - 0 The mass for the 1 is some 100 MeV below the negatively charged 1 measured by Crystal Barrel [13] and BNLE852 [11, 12]. Of course, this does not constitute compelling evidence for more than one exotic state in this mass region. However, if we assume that there are two states at 1280 and at 1380 MeV respectively, then one expects the lower-mass state to have a substantial branching ratio into (the reaction - + p 0 + n is mediated by exchange), whereas the higher-mass state would couple to the f2 (1270) channel or perhaps to the P omeron + channel (the reaction - + p - + p is produced in part by the exchange of the Pomeron).
CONCLUSIONS

Mass dependent fits of the D+ and P+ amplitudes and their relative phase using three different methods (described above) all lead to the conclusion that the P+ wave is well-described by a resonance hypothesis and is inconsistent with 0 having a constant phase. The resonance parameters for the observed 1 are given by M = 1257 ± 20 ± 25 MeV/c2 2 and = 354 ± 64 ± 58 MeV/c . Here the first error is statistical and the second is systematic. (We have chosen to take the resonant parameters and errors from method 2.) This result, together with the previous results from Crystal Barrel [14] and E852 [19] provide strong indications for one or more spin-exotic mesons near 1400 MeV/c2 decaying to 0 . We would like to express our appreciation to the members of the MPS group, and to the staffs of the AGS, BNL, and the collaborating institutions for their efforts. This research was supported in part by the National Science Foundation, the US Department of Energy, and the Russian Ministry for Education and Science. The Southeastern Universities Research Association (SURA) operates the Thomas Jefferson National Accelerator Facility for the United States Department of Energy under contract DE-AC05-84ER40150.

§ ¶ §§

[1] [2] [3] [4] [5] [6] [7] [8] [9] [10]

Present address: Department of Physics, Florida State University, Tallahassee, FL 32306 Deceased Present address: Department of Physics, University of Maryland, College Park, MD 20742 Permanent address: Rafael, Haifa, Israel Present address: TRIUMF, Vancouver, B.C., V6T 2A3, Canada Present address: Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 Present address: Laboratory for Nuclear Studies, Cornell University, Ithaca, NY 14853 Permanent address: Institute of High Energy Physics, Bejing, China Present address: Department of Physics, Kyungpook National University, Daegu, Korea R.L. Jaffe and K. Johnson, Phys. Lett. 60B, 201 (1976); R.L. Jaffe, Phys. Rev. D15, 267 (1977). T. Barnes and F.E. Close, Phys. Lett. B 116, 365 (1982). T. Barnes et al., Nucl. Phys. B 224, 241 (1983). N. Isgur and J. Paton, Phys. Rev. D 31, 2910 (1985). I. Balitsky et al., Z. Phys. C 33, 265 (1986); J.I. Latorre, P. Pascual, and S. Narison, Z. Phys, C 34, 347 (1987); J. Govaerts et al., Nucl. Phys. B 284, 674 (1987). T. Barnes et al., Phys. Rev. D 52, 5242 (1995). F.E. Close and P.R. Page, Nucl. Phys. B 443, 233 (1995). Y. Uehara et al., Nucl. Phys. A 606, 357 (1996). P. Lacock, C. Michael, P. Boyle and P. Rowland, Phys. Rev. D 54, 6997 (1996); C. Bernard et al., Nucl. Phys. B (Proc. Suppl.) 53, 228 (1997). T. Barnes, Acta Phys. Polon. B31, 2545 (2000).


6
[11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] D.R. Thompson et al., Phys. Rev. Lett. 79, 1630 (1997). S.U. Chung et al., Phys. Rev. D 60, 092001 (1999). A. Abele et al., Phys. Lett. B 423, 175 (1998). A. Abele et al., Phys. Lett. B 446, 349 (1999). D. Alde et al. Phys. Atom. Nucl., 62, 421 (1999). S.A. Sadovsky, Nucl. Phys. A 655, 131c (1999). G. M. Beladidze et al., Phys. Lett. B 313, 276 (1993). D.V. Amelin et al., Phys.Atom.Nucl. 68, 359 (2005). A.R. Dzierba et al., Phys. Rev. D 67, 094015 (2003). O.I. Dahl et al., "SQUAW kinematic fitting program", Univ. of California, Note P-126, unpublished (1968). J. P. Cummings and D. P. Weygand, "The New BNL PWA Programs", BNL-64637, unpublished (1997). S.U. Chung, Phys. Rev. D 56, 7299 (1997). S.A. Sadovsky, Soviet Physics Doklady 36, 537 (1991). W. -M. Yao et al., J. of Phys. G33, 1 (2006).