www.elsevier.com/locate/physletb
Instructions to authors
Aims and Scope Physics Letters B ensures the rapid publication of letter-type communications in the fields of Nuclear Physics, Particle Physics and Astrophysics. Articles should influence the physics community significantly. Submission Electronic submission is strongly encouraged. The electronic file, accompanied by a covering message, should be e-mailed to one of the Editors indicated below. Easy submission via the LANL-preprint server is certainly possible; please visit http://www.elsevier.com/locate/plbsubmission. If electronic submission is not feasible, submission in print is possible, but it will delay publication. In the latter case manuscripts (one original + two copies), accompanied by a covering letter, should be sent to one of the following Editors: L. Alvarez-Gaumé, Theory Division, CERN, CH-1211 Geneva 23, Switzerland, E-mail address: Luis.Alvarez-Gaume@CERN. CH Theoretical High Energy Physics (General Theory) J.-P. Blaizot, ECT*, Strada delle Tabarelle, 266, I-38050 Villazzano (Trento), Italy, E-mail address: [email protected]. FR Theoretical Nuclear Physics M. Cvetic,ˇ David Rittenhouse Laboratory, Department of Physics, University of Pennsylvania, 209 S, 33rd Street, Philadelphia, PA 19104-6396, USA, E-mail address: [email protected] Theoretical High Energy Physics M. Doser, EP Division, CERN, CH-1211 Geneva 23, Switzerland, E-mail address: [email protected] Experimental High Energy Physics D.F. Geesaman, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, IL 60439, USA, E-mail address: [email protected] Experimental Nuclear Physics H. Georgi, Department of Physics, Harvard University, Cambridge, MA 02138, USA, E-mail address: Georgi@PHYSICS. HARVARD.EDU Theoretical High Energy Physics G.F. Giudice, CERN, CH-1211 Geneva 23, Switzerland, E-mail address: [email protected] Theoretical High Energy Physics N. Glover, Institute for Particle Physics Phenomenology, Department of Physics, Science Laboratories, University of Durham, South Road, Durham DH1 3LE, UK, E-mail address: [email protected] Theoretical High Energy Physics W. Haxton, Institute for Nuclear Theory, Box 351550, University of Washington, Seattle, WA 98195-1550, USA, E-mail address: [email protected] Theoretical Nuclear Physics and Nuclear Astrophysics V. Metag, II. Physikalisches Institut, Universität Giessen, Heinrich-Buff-Ring 16, D-35392 Giessen, Germany, E-mail address: [email protected]. UNI-GIESSEN.DE Experimental Nuclear Physics L. Rolandi, EP Division, CERN, CH-1211 Geneva 23, Switzerland, E-mail address: [email protected] Experimental High Energy Physics W.-D. Schlatter, CERN, CH-1211 Geneva 23, Switzerland, E-mail address: [email protected] Experimental High Energy Physics
0370-2693/2005 Published by Elsevier B.V. doi:10.1016/S0370-2693(05)00566-6 vi Instructions to authors
H. Weerts, 3247 Biomedical and Physical Sciences Building, Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824-1111, USA, E-mail address: [email protected] Experimental High Energy Physics T. Yanagida, Department of Physics, Faculty of Science, University of Tokyo, Tokyo 113-0033, Japan, E-mail address: [email protected] Theoretical High Energy Physics The authors should indicate in which of the following four sections they would like to see their article published: Astrophysics & Cosmology, covered by L. Alvarez-Gaumé, M. Cveticˇ and W. Haxton. Experiments, covered by M. Doser, D.F. Geesaman, V. Metag, L. Rolandi, W.-D. Schlatter and H. Weerts. Phenomenology, covered by J.-P. Blaizot, M. Cvetic,ˇ H. Georgi, G.F. Giudice, N. Glover, W. Haxton and T. Yanagida. Theory, covered by L. Alvarez-Gaumé, J.-P. Blaizot, M. Cvetic,ˇ H. Georgi, G.F. Giudice, N. Glover, W. Haxton and T. Yanagida. For submissions in High Energy Physics authors are encouraged to provide the number of their Los Alamos preprint to the Physics Letters B Editor at the moment of submission. Original material. By submitting a paper for publication in Physics Letters B the authors imply that the material has not been published previously nor has been submitted for publication elsewhere and that the authors have obtained the necessary authority for publication. Refereeing. Submitted papers will be refereed and, if necessary, authors may be invited to revise their manuscript. If a submitted paper relies heavily on unpublished material, it would be helpful to have a copy of that material for the use of the referee. Publication speed The Editors and Publisher cooperate closely to ensure minimal publication delays. All proofreading will be done by the Publisher and proofs are not sent to the author(s). In order to keep delays to a minimum it is of utter importance for the author(s) closely to observe the guidelines given in the “Preparation of Manuscripts” below before submission of the manuscript and to supply an E-mail address and/or telefax number of the corresponding author. Preparation of manuscripts The following requirements as regards presentation of the manuscript should be met: (1) Manuscript: All manuscripts should be written in good English. The original typescript should be typed on one side of the paper, with double spacing and a wide margin. Instead of the original typescript, a copy of high quality (not a carbon copy) is also acceptable. Please adhere to the following order of presentation: Article title, Author(s), Affiliation(s), Abstract, Classification codes and keywords, Main text, Acknowledgements, Appendices, References, Figure captions, Tables. The name, complete postal address, telephone and fax numbers and the E-mail address of at least the corresponding author should be given on the first page of the manuscript. The title page should contain title, author(s), address(es) and abstract. The main text should start on a new page. All pages should be numbered. (2) Length: The total length of the paper should preferably not exceed six journal pages, equivalent to ten typewritten pages with double spacing, including the list of authors, abstract, references, figure captions and three figures. In the case that more figures are required, the text should be shortened accordingly. (3) Title: The title should be brief and such that it conveys to the informed reader the particular nature of the contents of the paper. (4) Address: The name(s) of the author(s) and the name and address of the institute where the research work was done should be indicated on the manuscript. The name of the author to whom correspondence is to be addressed should be underlined and an E-mail address and/or a telefax number supplied. (5) Abstract: An abstract of less than 60 words is required. It should contain the keywords of the paper as well as the essence of the results achieved. (6) Classification codes and keywords: Supply one to four classification codes (PACS and/or MSC) and up to six keywords of your own choice that describe the content of your article in more detail. (7) Formulae: Displayed formulae should be numbered and typed or clearly and unambiguously written by hand. Symbols requiring bold-face type, like vectors, etc., should be identified properly in the margin of the manuscript. Long equa- tions should be avoided as much as possible by introduction of suitable abbreviations of component expressions. Special attention should be paid to symbols that can easily be misread, such as i (lower case), I (cap.), 1 (el), 1 (one), (prime), ◦ o (lower case), O (cap.), 0 (zero), (degree), u,v (vee), ν (Greek nu), V (cap.), x, ×,X,z,Z,p,P,ρ (Greek rho), etc. (8) Footnotes: The footnotes may be typed at the end of the page on which they are alluded to, or at the end of the paper, or on a separate sheet. Please do not mix footnotes and references. Instructions to authors vii
(9) References: In the text, reference to other parts of the paper should be made by section (or equation) number, not by page number. References to other papers should be consecutively numbered in the text using square brackets and should be listed by number on a separate sheet at the end of the paper. Please do not combine multiple references to different papers into one numbered reference. The references should be as complete as possible and be presented as follows: For a book: B. de Wit and J. Smith, Field theory in particle physics, Vol. 1 (North-Holland, Amsterdam, 1986). For a paper in a journal: UAl Collab., G. Arnison et al., Phys. Lett. B 177 (1986) 244. For a paper in a contributed volume: R. Jackiw, in: Progress in quantum field theory, eds. H. Ezawa and S. Kamefuchi (North-Holland, Amsterdam, 1986) p. 83. For an unpublished paper: J.F. Gunion and H.E. Haber, UCD Report 86-12 (1986), unpublished. For a preprint: A. Lahanas and D.V. Nanopoulos, CERN preprint CERN-TH 4400/86 (1986). For a conference report: M.B. Green, Superstrings and the unification of forces and particles, in: Proc. fourth Marcel Gross- mann Meeting on General relativity (Rome, June 1985), Vol. 1, ed. R. Ruffini (North-Holland, Amsterdam, 1986) p. 203. (10) Figures: Each figure should also be submitted in triplicate: one master figure and two copies, the figure must be referred to in the text, be numbered and have a caption. The captions should be collected on a separate sheet. The appropriate place of each figure should be indicated in the margin. Axes of figures must be labelled properly. The (line) drawings for the figures must be submitted on separate sheets, drawn in black India ink and carefully lettered (with the use of stencils). The lettering as well as the essential details should have proportionate dimensions so as not to become illegible or unclear after the usual reduction by the printers (ideal lettering size after reduction of the drawing to one-column width is 1.8 mm). The drawings should preferably be of the same size as the typescript and designed for a reduction factor in print of two to three. The photographs should be originals, with somewhat more contrast than is required in the printed version. They should be unmounted unless part of a composite figure. Any scale markers should be inserted on the photograph, not drawn below it. The figures should be identified by the name of the first author, the journal name and the figure number. Instead of original drawings, sharp and contrasty glossy prints of about typescript size or high quality laserprints are also acceptable. If requested, original drawings will be returned to the author(s) upon publication of the paper. For detailed instructions on the preparation of electronic artwork, consult the Author Gateway from Elsevier at http:// authors.elsevier.com. (11) Colour illustrations: Illustrations in colour will be accepted in cases when the use of colours is judged by the Editor to be essential for the presentation. The Publisher and the author will each bear part of the extra costs involved. The costs charged to the authors of articles containing colour figures will be € 635 (approximately US$ 760) for the first page containing colour and € 318 for each additional page with colour, independent of the number of colour figures on each page. These prices are exclusive of Value Added Tax (VAT). Authors will be billed in Euros; the dollar price is for guidance only. The author receives 200 reprints of an article with colour illustrations free of charge. More reprints can be ordered at the usual rates for (black and white) reprints, there will be no additional charge for reprints containing colour illustrations. Colour illustrations should be submitted in the form of good quality colour photographs, transparencies, colour printer output, or 35 mm slides. Polaroid colour prints should be avoided. When supplying colour photographs and transparencies they should be close to the final size expected for publication. It should be noted that, in some cases, printing requirements will prevent figures from being located in the most preferred position in the text. (12) Tables: Tables should be typed on separate sheets and each table should have a number and a title. The appropriate places for the insertion of the tables should be indicated in the margin.
After acceptance Notification. You will be notified by the Editor of the journal of the acceptance of your article and invited to supply an electronic version of the accepted text, if this is not already available. Copyright transfer. You will be asked to transfer the copyright of the article to the Publisher. This transfer will ensure the widest possible dissemination of information. No proofs. In order to speed up publication, all proofreading will be done by the Publisher and proofs are not sent to the author(s). Enquiries. Visit the Author Gateway from Elsevier (http://authors.elsevier.com) for the facility to track accepted articles and set up e-mail alerts to inform you of when an article’s status has changed. The Author Gateway also provides detailed artwork guidelines, copyright information, frequently asked questions and more. Contact details for questions arising after acceptance of an article, especially those relating to proofs, are provided when an article is accepted for publication. viii Instructions to authors
Instructions for LATEX manuscripts The Publisher welcomes the receipt of an electronic version of your accepted manuscript (preferably encoded in LATEX). If you have not already supplied the final, accepted version of your article to the journal Editor, you are requested herewith to send a file with the text of the accepted manuscript directly to the Publisher by E-mail to the address given below. If the electronic file is suitable for processing by the Publisher, the article will be published without rekeying the full text. The article should be encoded in LATEX, preferably using the Elsevier document class ‘elsart’, or alternatively the standard document class ‘article’ or the document style ‘revtex’. The Elsevier LATEX package (including detailed instructions for LATEX preparation) can be obtained from Elsevier’s web site: www.elsevier.com/locate/latex, or from the Comprehensive TEX Archive Network (CTAN). The Elsevier package consists of the files: elsart.cls (use this file if you are using LaTeX2e, the current version of LATEX), elsart.sty and elsart12.sty (use these two files if you are using LaTeX2.09, the previous version of LATEX), instraut.dvi and/or instraut.ps and/or instraut.pdf (instruction booklet), model-harv.tex or model-num.tex (model files with instructions), template-harv.tex or template-num.tex (template files). To obtain the package from CTAN, use direct access via FTP at ftp.dante.de (Germany), ftp.tex.ac.uk (UK), or ctan.tug.org (Massachussets, USA) and go to the directory /tex-archive/macros/latex/contrib/supported/elsevier, or search for Elsevier with one of the CTAN search engines (http://ctan.tug.org/CTANfind.html, http://www.tex.ac.uk/CTANfind.html or http://www.dante.de/cgi-bin/ctan-index). CTAN is a mirrored network of ftp.tex.ac.uk, ftp.dante.de and ctan.tug.org, which are widely mirrored (see ftp://ctan.tug.org/ tex-archive/README.mirrors) and hold up-to-date copies of all the public-domain versions of TEX, LATEX, Metafont and ancillary programs. Questions concerning the LATEX author-prepared article and requests for the booklet with instructions to authors should be directed to the address given below. No changes from the version accepted by the Editor of the journal are permissible, without the prior and explicit approval by the Editor. Such changes should be clearly indicated on an accompanying printout of the file. The Publisher reserves the right to decide whether to process a manuscript from authors’ files or not. Articles coded in a simple manner with no user-defined macros are most likely to be handled this way. If sent via electronic mail, files should be accompanied by a clear identification of the article (name of journal, Editor’s reference number) in the “subject field” of your electronic-mail message. Authors should include an ASCII table (available from the Publisher) in their files, to enable any transmission errors to be detected. Addresses Ð Mail: Drs. E.J. van der Wolk, Physics Letters B, Editorial-Production Journals Department, Elsevier B.V., P.O. Box 2759, 1000 CT Amsterdam, The Netherlands. E-mail: [email protected]. Author benefits No page charges. Publishing in Physics Letters B is free. Free offprints. The corresponding author will receive 25 offprints free of charge. An offprint order form will be supplied by the Publisher for ordering any additional paid offprints. Discount. Contributors to Elsevier journals are entitled to a 30% discount on all Elsevier books. Contents Alert. Physics Letters B is included in Elsevier’s pre-publication service Contents Alert (for information, please contact: [email protected]). ContentsDirect. Physics Letters B is included in Elseviers free E-mail service ContentsDirect. Please register at http://www. elsevier.com/locate/ContentsDirect. Publication scheme Physics Letters B will be published weekly. All correspondence with the Editors or Publisher should contain full reference to the paper concerned, i.e., the names of all the contributors, the full title, as well as the reference number. Contributors are reminded that, once their contribution has been accepted for publication, all further correspondence should be sent directly to the Publisher and not to the Editors (Drs. E.J. van der Wolk, Physics Letters B, Editorial-Production Journals Department, Elsevier B.V., P.O. Box 2759, 1000 CT Amsterdam, The Netherlands; telephone number +31 20 485 2634; telefax number +31 20 485 2431; electronic mail: [email protected] (the receipt of an E-mail message will be acknowledged; in the absence of this acknowledgement, the mailer should assume that his message was never received by the Publisher)). A prompt reply may be expected whenever appropriate. If no reaction is received within three weeks, authors are kindly requested to seek contact again. PHYSICS LETTERS B
EDITORS
L. ALVAREZ-GAUMÉ J.-P. BLAIZOT M. CVETICˇ GENEVA VILLAZZANO (TRENTO) PHILADELPHIA, PA
M. DOSER D.F. GEESAMAN H. GEORGI GENEVA ARGONNE, IL CAMBRIDGE, MA
G.F. GIUDICE N. GLOVER W. HAXTON GENEVA DURHAM SEATTLE, WA
V. METAG L. ROLANDI W.-D. SCHLATTER GIESSEN GENEVA GENEVA
H. WEERTS T. YANAGIDA EAST LANSING, MI TOKYO
VOLUME 615, 2005
Amsterdam – Boston – Jena – London – New York – Oxford Paris – Philadelphia – San Diego – St. Louis Physics Letters B 615 (2005) 1–13 www.elsevier.com/locate/physletb
Nutty dyons
Yves Brihaye a, Eugen Radu b
a Physique-Mathématique, Universite de Mons-Hainaut, Mons, Belgium b Department of Mathematical Physics, National University of Ireland Maynooth, Maynooth, Ireland Received 25 February 2005; received in revised form 6 April 2005; accepted 7 April 2005 Available online 15 April 2005 Editor: N. Glover
Abstract We argue that the Einstein–Yang–Mills–Higgs theory presents nontrivial solutions with a NUT charge. These solutions ap- proach asymptotically the Taub–NUT spacetime and generalize the known dyon black hole configurations. The main properties of the solutions and the differences with respect to the asymptotically flat case are discussed. We find that a non-Abelian mag- netic monopole placed in the field of gravitational dyon necessarily acquires an electric field, while the magnetic charge may take arbitrary values. 2005 Elsevier B.V. All rights reserved.
1. Introduction
A feature of certain gauge theories is that they admit classical solutions which are interpreted as represent- ing magnetic monopoles. For non-Abelian gauge fields interacting with a Higgs scalar, there exist even regular configurations with a finite mass, as proven by the famous ’t Hooft–Polyakov solution [1]. Typically, the mag- netic monopoles admit also electrically charged generalizations—so-called dyons, the Julia–Zee solution [2] of the SU(2)-Higgs theory possibly being the best known case. These solutions admits also gravitating generalizations, both regular and black hole solutions being considered in the literature (see [3] for a general review of this top- ics). In SU(2)-Einstein–Yang–Mills–Higgs (EYMH) theory, a branch of globally regular gravitating dyons emerges smoothly from the corresponding flat space solutions. The non-Abelian black hole solutions emerge from the glob- ally regular configurations, when a finite regular event horizon radius is imposed [4,5]. These solutions cease to exist beyond some maximal value of the coupling constant α (which is proportional to the ratio of the vector meson mass and Planck mass).
E-mail address: [email protected] (Y. Brihaye).
0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.04.016 2 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13
It has been speculated that such configurations might have played an important role in the early stages of the evolution of the Universe. Also, various analyses indicate that the monopole and dyon solutions are important in quantum theories. Since general relativity shares many similarities with gauge theories, one may ask whether Einstein’s equations present solutions that would be the gravitational analogous of the magnetic monopoles and dyons. The first example of such a solution was found in 1963 by Newman, Unti and Tamburino (NUT) [6,7]. This metric has become renowned for being “a counterexample to almost anything” [8] and represents a generalization of the Schwarzschild vacuum solution [9] (see [10] for a simple derivation of this metric and historical review). It is usually interpreted as describing a gravitational dyon with both ordinary and magnetic mass.1 The NUT charge which plays a dual role to ordinary mass, in the same way that electric and magnetic charges are dual within Maxwell theory [11].By continuing the NUT solution through its horizon one arrives in the Taub universe [7], which may be interpreted as a homogeneous, nonisotropic cosmology with the spatial topology S3. As discussed by many authors (see, e.g., [13,14]), the presence of magnetic-type mass (the NUT parameter n) introduces a “Dirac-string singularity” in the metric (but no curvature singularity). This can be removed by appropriate identifications and changes in the topology of the spacetime manifold, which imply a periodic time coordinate. Moreover, the metric is not asymptotically flat in the usual sense although it does obey the required fall-off conditions. A large number of papers have been written investigating the properties of the gravitational analogs of magnetic monopoles [15,16], the vacuum Taub–NUT solution being generalized in different directions. The corresponding configuration in the Einstein–Maxwell theory has been found in 1964 by Brill [17]. This Abelian solution has been generalized for the matter content of the low-energy string theory, a number of NUT-charged configurations being exhibited in the literature (see, e.g., [18] for a recent example and a large set of references). A discussion of the non-Abelian counterparts of the Brill solution is presented in [19]. These configurations generalize the well-known SU(2)-Einstein–Yang–Mills hairy black hole solutions [20], presenting, as a new feature, a nontrivial electric potential. However, the “no global non-Abelian charges” results found for asymptotically flat EYM static configurations [21] are still valid in this case, too. Here we present arguments for the existence of NUT-charged generalizations of the known EYMH black hole solutions [4,5]. Apart from the interesting question of finding the properties of a Yang–Mills–Higgs dyon in the field of a gravitational dyon, there are a number of other reasons to consider this type of solutions. In some super- symmetric theories, closure under duality forces us to consider NUT-charged solutions. Furthermore, dual mass solutions play an important role in Euclidean quantum gravity [22] and therefore cannot be discarded in spite of their causal pathologies. Also, by considering this type of asymptotics, one may hope to attain more general features of gravitating non-Abelian dyons. The Letter is structured as follows: in the next section we present the general framework and analyse the field equations and boundary conditions. In Section 3 we present our numerical results. We conclude with Section 4, where our results are summarized.
2. General framework and equations of motion
2.1. Action principle
The action for a gravitating non-Abelian SU(2) gauge field coupled to a triplet Higgs field with vanishing Higgs self-coupling is
1 Note that the Taub–NUT spacetime plays also an important role outside general relativity. For example, the asymptotic motion of monopoles in (super-)Yang–Mills theories corresponds to the geodesic motion in a Euclideanized Taub-NUT background [12].However,these developments are outside the interest of this work. Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 3 √ R 1 1 S = −gd4x − Tr F F µν − Tr D ΦDµΦ , (1) 16πG 2 µν 4 µ with Newton’s constant G. The field strength tensor is given by Fµν = ∂µAν − ∂νAµ − ie[Aµ,Aν], with Dµ = ∂µ − ie[Aµ, ] being the covariant derivative and e the Yang–Mills coupling constant. µν Varying the action (1) with respect to g , Aµ and Φ we have the field equations
1 1 √ 1 1 √ R − g R = 8πGT , √ D −gF µν = ie Φ,DνΦ , √ D −gDµΦ = 0, µν 2 µν µν −g µ 4 −g µ (2) where the stress-energy tensor is 1 1 1 T = 2Tr F F gαβ − g F F αβ + Tr D ΦD Φ − g D ΦDαΦ . (3) µν µα νβ 4 µν αβ 2 µ ν 4 µν α
2.2. Metric ansatz and symmetries
We consider NUT-charged spacetimes whose metric can be written locally in the form dr2 θ 2 ds2 = + P 2(r) dθ2 + sin2 θdϕ2 − N(r)σ2(r) dt + 4n sin2 dϕ , (4) N(r) 2 the NUT parameter n being defined as usually in terms of the coefficient appearing in the differential dt + 4n sin2(θ/2)dϕ.Hereθ and ϕ are the standard angles parametrizing an S2 with ranges 0 θ π,0 ϕ 2π. Apart from the Killing vector K0 = ∂t , this line element possesses three more Killing vectors characterizing the NUT symmetries θ K = sin ϕ∂ + cos ϕ cot θ∂ + 2n cos ϕ tan ∂ , 1 θ ϕ 2 t θ K = cos ϕ∂ − sin ϕ cot θ∂ − 2n sin ϕ tan ∂ , 2 θ ϕ 2 t K3 = ∂ϕ − 2n∂t . (5) These Killing vectors form a subgroup with the same structure constants that are obeyed by spherically symmetric solutions [Ka,Kb]= abcKc. The n sin2(θ/2) term in the metric means that a small loop around the z-axis does not shrink to zero at θ = π. This singularity can be regarded as the analogue of a Dirac string in electrodynamics and is not related to the usual degeneracies of spherical coordinates on the two-sphere. This problem was first encountered in the vacuum NUT metric. One way to deal with this singularity has been proposed by Misner [8]. His argument holds also independently of the precise functional form of N and σ . In this construction, one considers one coordinate patch in which the string runs off to infinity along the north axis. A new coordinate system can then be found with the string running off to infinity along the south axis with t = t + 4nϕ, the string becoming an artifact resulting from a poor choice of coordinates. It is clear that the t coordinate is also periodic with period 8πnand essentially becomes an Euler angle coordinate on S3. Thus an observer with (r,θ,ϕ)= const follows a closed timelike curve. These lines cannot be removed by going to a covering space and there are no reasonable spacelike surface. One finds also that surfaces of constant radius have the topology of a three-sphere, in which there is a Hopf fibration of the S1 of time over the spatial S2 [8]. Therefore for n different from zero, the metric structure (4) generically shares the same troubles exhibited by the vacuum Taub–NUT gravitational field [23], and the solutions cannot be interpreted properly as black holes. 4 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13
2.3. Matter fields ansatz
While the Higgs field is given by the usual form
Φ = φτ3, (6) the computation of the appropriate SU(2) connection compatible with the Killing symmetries (5) isamorein- volved task. This can be done by applying the standard rule for calculating the gauge potentials for any spacetime group [24,25]. According to Forgacs and Manton, a gauge field admit a spacetime symmetry if the spacetime trans- L = L formation of the potential can be compensated by a gauge transformation [24] Ki Aµ DµWi , where stands for the Lie derivative. Taking into account the symmetries of the line element (4) we find the general form 1 2 θ A = dt + 4n sin dϕ u(r)τ3 + ν(r)τ3 dr + ω(r)τ1 +˜ω(r)τ2 dθ 2e 2 + cos θτ3 + ω(r)τ2 −˜ω(r)τ1 sin θ dϕ . (7)
This gauge connection remains invariant under a residual U(1) gauge symmetry which can be used to set ν = 0. Also, because the variables ω and ω˜ appear completely symmetrically in the EYMH system, the two amplitudes must be proportional and we can always set ω˜ = 0 (after a suitable gauge transformation). Thus, similar to the n = 0 case, the gauge potential is described by two functions ω(r) and u(r) which we shall refer to as magnetic and electric potential, respectively.
2.4. Field equations and known solutions
Within the above ansatz, the classical equations of motion can be derived from the following reduced action 2 3 1 n σ N S = dr dt σ 1 − NP 2 − PP N + 2P (σ NP ) + 8πG P 2 2 2 2 2 2 2 1 σ(ω − 1 + 2nu) P u ω u 1 − Nσω 2 + − − + σNP2φ 2 + σω2φ2 , (8) e2 2P 2 2σ 2 σN 2 where the prime denotes the derivative with respect to the radial√ coordinate r. At this point, we fix the metric gauge by choosing P(r)= r2 + n2, which allows a straightforward analysis of the relation with the Abelian configurations. Dimensionless quantities are obtained by considering the rescalings r → r/(ηe), φ → φη, n → n/(ηe), u → ηeu (where η is the asymptotic√ magnitude of the Higgs field). As a result, the field equations depend only on the coupling constant α = 4πGη. The EYMH equations reduce to the following system of five nonlinear differential equations 2 n N rN = − N + σ 2 − 1 2 3 1 P 2 2 2 2 1 P u ω u 1 − 2α2 Nω 2 + ω2 − 1 + 2nu 2 + + + NP2φ 2 + ω2φ2 , 2P 2 2σ 2 σ 2N 2 2 2 2 2 2 n σ(1 − σ ) α σ 2ω u σ = + P 2φ 2 + 2ω 2 + , rP2 r σ 2N 2 2 2 (ω − 1 + 2nu) u (Nσ ω ) = σω + φ2 − , P 2 σ 2N Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 5 NσP2φ = 2σω2φ, P 2u 2ω2u 2nσ = − ω2 − 1 + 2nu . (9) σ σN P 2 Two explicit solution of the above equations are well known. The vacuum Taub–NUT one corresponds to 2(Mr + n2) ω(r) =±1,u(r)= 0,σ(r)= 1,φ(r)= 1,N(r)= 1 − . (10) r2 + n2 The U(1) Brill solution [17] has the form nQ − Q r ω(r) = 0,u(r)= u + m e ,σ(r)= 1,φ(r)= 1, 0 r2 + n2 2(Mr + n2) α2(Q2 + Q2 ) N(r)= 1 − + e m , (11) r2 + n2 4(r2 + n2) and describes a gravitating Abelian dyon with a mass M, electric charge Qe and magnetic charge Qm ≡ 1 − 2u0n, u0 being an arbitrary constant, corresponding to the asymptotic value of the electric potential. It can be stressed that the Brill solution possesses two, one or zero horizons, according to the values of the free parameters Qe, M, u(∞). In the same way as in the case of Reissner–Nordström solutions, the extremal Brill solution can be defined as the solutions with a degenerate horizon at r = r0. This gives the following conditions, fixing M and r0 α2 r = M, M2 + n2 − Q2 + Q2 = 0. (12) 0 4 e m As we will see later, it is convenient to further specify the arbitrary constant u(∞) in such a way the u(r0) = 0, this implying 1 − Q nQ − MQ m + m e = 0, (13) 2n M2 + n2 which fixes Qm and leaves Qe as the only remaining free parameter. In the following we will refer to this solution as to the extremal Brill solution. As far as we could see, it is not possible to express M and Qm in a closed form depending on (α,n,Qe), but the solution can be constructed numerically.
2.5. Boundary conditions
We want the metric (4) to describe a nonsingular, asymptotically NUT spacetime outside an horizon located at r = rh.HereN(rh) = 0 is only a coordinate singularity where all curvature invariants are finite. A nonsingular extension across this null surface can be found just as at the event horizon of a black hole. If the time is chosen to be periodic, as discussed above, this surface would not be a global event horizon, although it would still be an apparent horizon. The regularity assumption implies that all curvature invariants at r = rh are finite. The corresponding expansion as r → rh is 2 2 N(r)= N1(r − rh) + O(r − rh) ,σ(r)= σh + σ1(r − rh) + O(r − rh) , 2 2 3 ω(r) = ωh + ω1(r − rh) + O(r − rh) ,u(r)= u1(r − rh) + u2(r − rh) + O(r − rh) , 2 φ(r)= φh + φ1(r − rh) + O(r − rh) , (14) 2 = 2 + 2 where Ph rh n and 1 (ω2 − 1)2 1 u2P 2 N = − α2 h + 1 h + ω φ2 , 1 1 2 2 2 h h rh 2Ph 2 σh 6 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 2 2 2 2 2 n σh(1 − σ ) α σ 2ω u σ = h + h P 2φ2 + ω2 + 1 1 , 1 2 h 1 2 1 2 2 rhP rh σ N h h 1 2 2 2 2 2 ω ω − 1 σ u nσ (ω − 1) u1ω u r 2ω φh ω = h h + φ2 ,u= 1 1 − h h + h − 1 h ,φ= h , 1 2 h 2 4 2 2 1 2 (15) N1 Ph 2σh Ph N1Ph Ph N1Ph σh, u1, ωh, φh being arbitrary parameters. The analysis of the field equations as r →∞gives the following expression in terms of the constants c, u0, Qe, ˜ φ1, M 2M 2n2 − α2(φ˜2 + (1 − 2nu2)2 + Q2) M(2n2 + α2φ˜2) N(r)∼ 1 − − 1 0 e + 1 +···, r r2 r3
2 ˜2 2 ˜2 α φ 4α φ M − 1−u2r σ ∼ 1 − 1 − 1 +···,ω(r)∼ ce 0 +···, 2r2 3r3 ˜ ˜ 2 2 ˜2 φ φ M Q n(1 − 2nu ) Qe(6n + α φ ) φ ∼ 1 − 1 + 1 +···,u(r)∼ u − e + 0 − 1 +···. (16) r r2 0 r r2 6r3 Note that similar to the n = 0 asymptotically flat case, the magnitude of the electric potential at infinity cannot 2 exceed that of the Higgs field, |u0| < 1. The constant M appearing in the asymptotic expansion of the metric function N(r) can be interpreted as the total mass of solutions (this can be proven rigorously by applying the general formalism proposed in [28]). Note that M and n are unrelated on a classical level. Also, no purely monopole solution can exist for a nonvanishing NUT charge (i.e., one cannot consistently set u = 0 unless ω =±1, in which case the vacuum Taub–NUT solution is recovered). Thus, a non-Abelian magnetic monopole placed in the field of gravitational dyon necessarily acquires an electric field. We close this section by remarking that the definition of the non-Abelian charges is less clear for n = 0. Although we may still define a ’t Hooft field strength tensor, in the absence of a nontrivial two-sphere at infinity on which to integrate, the only reasonable definition the non-Abelian magnetic and electric charges is in terms of the asymptotic (3) 2 behavior of the gauge field. By analogy to the asymptotically flat case, Qe and Qm are defined from Ftr Qe/r (3) = − and Fθφ Qm sin θ (a similar problem occurs for an U(1) field [27]). Thus, since Qm 1 2nu0, the usual quantization relation for the magnetic charge is lost for n = 0, which is a consequence of the pathological large scale structure of a NUT-charged spacetime.
3. Numerical results
Although an analytic or approximate solution appears to be intractable, we present in this section numerical arguments that the known EYMH black hole solutions can be extended to include a NUT parameter. The equations of motion (9) have been solved for a large set of the parameters (α, n, Qe,rh), looking for solu- tions interpolating between the asymptotics (14) and (16). NUT-charged solutions are found for any n = 0EYMH dyonic black hole configuration by slowly increasing the parameter n (since the transformation n →−n leaves the field equations unchanged except for the sign of the electric potential, we consider here only positive values of n). As expected, these configurations have many features in common with the n = 0 solutions discussed in [4]; they also present new features that we will pointed out in the discussion. Typical profiles for the metric functions N(r) and σ(r) and for the electric potential u(r) are presented in Fig. 1, for a dyonic black hole solution as well as for two NUT-charged solutions. The gauge function ω(r) and the Higgs scalar φ(r) interpolates monotonically between some constant values on the event horizon and zero respectively one at infinity, without presenting any local extremum (see Fig. 3).
2 This depends on the asymptotic structure of the spacetime. For example, in an anti-de Sitter spacetime, u0 may take arbitrary values [26]. Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 7
Fig. 1. The functions N(r), σ(r) and u(r) are plotted for three typical solutions for the same values of (rh,Qe,α).
The domain of existence of the non-Abelian nutty dyons can be determined in the space of parameters. If we fix the electric charge Qe of the solution, then there likely exist a volume VQ in the parameter space of (α,n,rh) inside which non-Abelian solutions exist and on the side of which they become singular and/or bifurcate into Abelian solution of the type of the Brill solution. For n = 0 the domain of the (α, rh) plane where non-Abelian solutions exist was determined in [5] for Qe = 0 and in [4] for Qe = 0. The determination of VQ is of course a huge task. In this Letter, we will not attempt to determine the shape of VQ accurately but rather attempt to sketch it by analyzing the pattern of solutions on some generic lines in the space of parameters. For definiteness we set Qe = 0.2 in our numerical analysis, although nontrivial solutions have been found also for other values of the electric charge.
3.1. n varying
First, we have integrated the system of equations (9) with fixed values for α, rh and Qe and increased the NUT charge n. Our values here are α = 1.0, rh = 0.2 and Qe = 0.2 corresponding to a generic values for the parame- ters (the corresponding n = 0 gravitating dyon was constructed in [4]). As far as the function u(r) is concerned, there exists a main difference between the case n = 0 and n = 0. Indeed in the case n = 0 this function behave 4 asymptotically like u(r) ∼ u0 + Qe/r + O(1/r ) while in the presence of a NUT charge the behaviour is instead 2 u(r) ∼ u0 + Qe/r + K/r , where, as seen from (16), the constant K increases with n. Thus, when n becomes large, it becomes more difficult to construct numerical solutions with a good enough accuracy,3 for a given value of the electric charge Qe. The effect of increasing n apparently depend strongly of the value α.Forα small (typically α 1) the pattern can be summarized by the following points: (i) No local extrema of N(r) are found for small enough values of n. For larger n, the function N(r) develops a local maximum and also a local minimum, say NM and Nm at some intermediate, n-depending values of r.Forn large enough, we have Nmax > 1. No local minimum of N persist for large enough n, the minimum of N(r) (Nm = 1) being attained as r →∞. (ii) The second metric function σ(r) still remains monotonically increasing but the value σ(rh) diminishes when n increases. (iii) The asymptotic
3 To integrate the equations, we used the differential equation solver COLSYS which involves a Newton–Raphson method [29]. 8 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13
Fig. 2. The values of the parameters M, σ(rh), Nm, NM , φ(rh), ω(rh) and u(∞) = u0 are shown as a function of n for solutions with rh = 0.2, Qe = 0.2 and two different values of the coupling constant α. value u(∞) also decreases for increasing n. With the values choosen, we have u(∞) ≈ 0.189 for n = 0; we find. u(∞) = 0forn ≈ 0.25 and negative values for larger n. These effects are illustrated on Fig. 2a. On this figure we have set 0 Fig. 3. The metric functions N(r), σ(r)and the matter functions ω(r), φ(r) and u(r) are shown as a function of r for fixed values of (rh,Qe,n) and three different values of α. The functions N(r) and u(r) of the corresponding extremal Brill solution (with ω(r) = 0,h(r)= 1) are also exhibited. Nevertheless, it seems that there are two possible patterns for n →∞: for values of α smaller than a critical value αˆ , solutions with large values of n seem to occur, while for α>αˆ , the solutions bifurcate into a Brill solution (for Qe = 0.2 we find αˆ ∼ 1.5). The occurence of these two patterns is reminiscent to the case of n = 0 gravitating dyons. Note also that, as shown in these plots, the mass parameter M takes negative values for large enough values of n. This is not a surprise, since something similar happens already for the U(1) Brill solution (11). 3.2. α varying We now discuss the behaviour of the solutions for a varying α and the other parameters fixed. In absence of ∈[ + 2 ] a NUT charge it is know [4,5] that non-Abelian dyonic black hole exist for rh 0, 3 4Qe/2 . For fixed Qe and rh and increasing α they bifurcates into an extremal Reissner–Nordström solution at α ∼ αc.Thevalueαc depends of course on rh and Qe.For rh 1thevalueαc ≈ 1.4 is found numerically and depends weakly on Qe. ∼ + 2 ≈ + 2 + 2 For rh 3 4Qe/2wehaveαc (3 4Qe)/(1 Qe)/2. For n>0 we see (e.g., on Fig. 3) that the local maximum characterizing the function N(r) of a nutty solution (at least for large enough values of the NUT charge n) progressively disappears in favor of a local minimum when α increases. This minimum appears far outside the event horizon r = rh and becomes deeper. In fact, the minimal value Nm approaches zero when α tends to a critical value, say αc(n, Qe,rh). If we denote by rm the value of the radial variable where N(rm) = 0 (with rm >rh) our numerical results strongly indicate that the non-Abelian solution converges into an extremal Brill solution on the interval r ∈[rm, ∞) for α → αc. Indeed, the matter functions’ profiles u, w, φ and the metric functions σ,N all approaches the profiles of the corresponding extremal Brill solution with the same αc,Qe,n. This result is illustrated on Fig. 3 for n = 1, rh = 0.2 and Qe = 0.2; in this case, we find αc ≈ 2.35 but we believe that the result holds for generic values of (n, rh,Qe). The determination of the critical value αc(n, rh,Qe) is not aimed in this Letter. However, it seems that the value αc depends weakly of n, for example, we find αc ≈ 2.22 for n ∈[3, 4]. Nevertheless we can conclude that nutty dyons exist on a finite interval of α and bifurcate into extremal Brill solutions for α = αc. 10 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 Fig. 4. The profiles of the functions N(r), σ(r) and u(r) are represented for Qe = 0.2, n = 0.5, α = 1 and three different values of rh. 3.3. rh varying In the case of gravitating dyonic black holes with event horizon rh, the solutions approach the corresponding regular gravitating solution on the interval ]rh, ∞[ when the limit rh → 0 is considered. It is therefore a natural question to investigate how nutty-dyons behave in the same limit. Considering this problem for a few generic values of (α, n) we reach the conclusion that, in the limit rh → 0, the nutty dyon becomes singular at r = 0 because the value σ(rh) tends to zero. This situation is illustrated on Fig. 4 where the functions N(r), σ(r) and u(r) are plotted for three different values of rh and α = 1, n = 0.5. Remarkably, this figure reveals that the functions σ(r) and u(r) are rather independant of rh (it is also true for w(r), φ(r) which are not represented)√ while the function N(r) indeed involves nontrivially with rh.Notealso that for the metric gauge choice P(r)= r2 + n2, the area of two-sphere dΩ2 = P 2(r)(dθ 2 + sin2 θdϕ2) does not vanish at r = 0. However, by choosing a Schwarzschild gauge choice P(r)= r, a straightforward analysis of the corresponding field equations (which can easily be derived from (8)) implies that it is not possible to take a consistent set of boundary conditions at r = 0 without introducing a curvature singularity at that point. Therefore, no globally regular EYMH solutions are found for n = 0. The determination of the domain of nutty dyons for fixed (α,n,Qe) and increasing the horizon radius rh is very likely an involved problem. Already in the case n = 0, discussed in [4] the numerical analysis reveals several (up to three) branches of solutions on some definite intervals of the parameter rh. We believe that similar patterns could occur for n>0 but their analysis is out of the scope of this Letter. 4. Further remarks In this work we have analysed the basic properties of gravitating YMH system in the presence of a NUT charge. We have found that despite the existence of a number of similarities to the n = 0 case (for example, the presence of a maximal value of the coupling constant α), the NUT-charged solutions exhibits some new qualitative features. The static nature of a n = 0 spherically symmetric gravitating non-Abelian solution implies that it can only produce a “gravitoelectric” field. There both non-Abelian monopole and dyon black hole solutions are possible to exist, with a well defined zero event horizon radius limit. For a nozero NUT charge, the existence of the cross Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 11 metric term gϕt shows that the solutions have also a “gravitomagnetic” field. The term gϕt does not produce an ergoregion but it will induce an effect similar to the dragging of inertial frames [30]. In this case we have found that only non-Abelian dyons are possible to exist and the usual magnetic charge quantization relation is lost. The total mass of these solutions may be negative and the configurations do not survive in the limit of zero event horizon radius. A discussion of possible generalizations of this work should start with the radially excited nutty dyons, for which the gauge function ω(r) possesses nodes. These configurations are very likely to exist, continuing for n>0 the excited configurations discussed in [4]. Also, in our analysis, to simplify the general picture, we set the Higgs potential V(φ)to zero. We expect to find the same qualitative results for a nonvanishing scalar potential (at least if the parameters are not to large). It would be a challenge to construct axially symmetric NUT-charged dyons (the corresponding n = 0 monopoles are discussed in [31]). Such dyon solutions would present a nonvanishing angu- lar momentum, generalizing the Abelian Kerr-Newman–NUT configurations (a set of asymptotically flat rotating solutions have been considered recently in [32]). Similar to the case n = 0, the solutions discussed in this work can also be generalized by including a more general matter content. However, we expect that these more general configurations will present the same generic properties discussed in this work. This may be important, since there are many indications that the NUT charge is an important ingredient in low energy string theory [27], conclusion enhanced by the discovery of “duality” transformations which relate superficially very different configurations. In many situations, if the NUT charge is not included in the study, some symmetries of the system remain unnoticed (see, e.g., [33] for such an example). Therefore, we may expect the NUT charge to play a crucial role in the duality properties of a (supersymmetric-) theory presenting gravitating non-Abelian dyons. Unfortunately, the pathology of closed timelike curves is not special to the vacuum Taub–NUT solution but afflicts all solutions of Einstein equations solutions with “dual” mass in general [28]. This condition emerges only from the asymptotic form of the fields, and is completely insensitive to the precise details of the nature of the source, or the precise nature of the theory of gravity at short distances where general relativity may be expected to break down [23]. This a causal behavior precludes the nutty dyons solutions discussed in this Letter having a role classically and implies a number of pathological properties of these configurations. Nevertheless, there are various features suggesting that the Euclidean version of NUT-charged solutions play an important role in quantum gravity [22]. For example, the entropy of such solutions generically do not obey the simple “quarter-area law”. As usual, a positive-definite metric is found by considering in (4) the analytical continuation t → it, n → in, which gives P 2(r) = r2 − n2. In this case, the absence of conical singularities at the root rh of the function N(r) imposes a periodicity in the Euclidean time coordinate = 4π β , (17) N (rh)σ (rh) which should be equal with the one to remove the Dirac string β = 8πn. In the usual approach, the solution’s parameters must be restricted such that the fixed point set of the Killing vector ∂t is regular at the radial position r = rh. We find in this way two types of regular solutions, “bolts” (with arbitrary rh = rb >n) or “nuts” (rh = n), depending on whether the fixed point set is of dimension two or zero (see [13] for a discussion of these solutions in the vacuum case and [34] for a recent generalization with anti-de Sitter asymptotics). We expect that the Euclidean nutty dyons will present some new features as compared to the Lorentzian coun- terparts. For example, globally regular solutions may exist in this case, since r = rh corresponds to the origin of the coordinate system (note also that the SU(2) Yang–Mills system is known to present self-dual solutions in the background of a vacuum Taub–NUT instanton [35]). In the absence of closed form solutions, the properties of these non-self dual EYMH solutions cannot be predicted directly from those of the Lorentzian configurations. However, similar to the Lorentzian case, they can be studied in a systematic way, by using both analytical and numerical arguments. For example, the magnitude of the electric potential at infinity of the Euclidean solutions, is 12 Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 not restricted. Also, the condition β = 8πn implies N (rh)σ (rh) = 2n and introduces a supplementary constraint on the matter functions as r → rh. A study of such solutions may be important in a quantum gravity context. Acknowledgements E.R. thanks D.H. Tchrakian for useful discussions. Y.B. is grateful to the Belgian FNRS for financial support. The work of E.R. is carried out in the framework of Enterprise-Ireland Basic Science Research Project SC/2003/390 of Enterprise-Ireland. References [1] G. ’t Hooft, Nucl. Phys. B 79 (1974) 276; A.M. Polyakov, JETP Lett. 20 (1974) 194, Pis’ma Zh. Eksp. Teor. Fiz. 20 (1974) 430 (in Russian). [2] B. Julia, A. Zee, Phys. Rev. D 11 (1975) 2227. [3] M.S. Volkov, D.V. Gal’tsov, Phys. Rep. 319 (1999) 1; D.V. Gal’tsov, hep-th/0112038. [4] Y. Brihaye, B. Hartmann, J. Kunz, N. Tell, Phys. Rev. D 60 (1999) 104016, hep-th/9904065. [5] P. Breitenlohner, P. Forgacs, D. Maison, Nucl. Phys. B 383 (1992) 357. [6] E.T. Newman, L. Tamburino, T. Unti, J. Math. Phys. 4 (1963) 915. [7] C.W. Misner, J. Math. Phys. 4 (1963) 924; C.W. Misner, A.H. Taub, Sov. Phys. JETP 28 (1969) 122. [8] C.W. Misner, in: J. Ehlers (Ed.), Relativity Theory and Astrophysics I: Relativity and Cosmology, in: Lectures in Applied Mathematics, vol. 8, American Mathematical Society, Providence, RI, 1967, p. 160. [9] S.W. Hawking, G.F.R. Ellis, The Large Structure of Space–Time, Cambridge Univ. Press, Cambridge, 1973. [10] D. Lynden-Bell, M. Nouri-Zonoz, Rev. Mod. Phys. 70 (1998) 427, gr-qc/9612049. [11] M. Damianski, E.T. Newman, Bull. Acad. Pol. Sci. 14 (1966) 653; J.S. Dowker, Gen. Relativ. Gravit. 5 (1974) 603. [12] G.W. Gibbons, N.S. Manton, Nucl. Phys. B 274 (1986) 183. [13] S.W. Hawking, C.J. Hunter, Phys. Rev. D 59 (1999) 044025, hep-th/9808085. [14] R.B. Mann, Phys. Rev. D 60 (1999) 104047, hep-th/9903229. [15] D.J. Gross, M.J. Perry, Nucl. Phys. B 226 (1983) 29. [16] R.D. Sorkin, Phys. Rev. Lett. 51 (1983) 87. [17] D.R. Brill, Phys. Rev. 133 (1964) B845. [18] C.V. Johnson, H.G. Svendsen, Phys. Rev. D 70 (2004) 126011, hep-th/0405141. [19] E. Radu, Phys. Rev. D 67 (2003) 084030, hep-th/0211120. [20] M.S. Volkov, D.V. Galtsov, JETP Lett. 50 (1989) 346; H.P. Kuenzle, A.K. Masood-ul-Alam, J. Math. Phys. 31 (1990) 928; P. Bizon, Phys. Rev. Lett. 64 (1990) 2844. [21] P. Bizon, O.T. Popp, Class. Quantum Grav. 9 (1992) 193; A.A. Ershov, D.V. Galtsov, Phys. Lett. A 150 (1990) 159; D.V. Galtsov, A.A. Ershov, Phys. Lett. A 138 (1989) 160. [22] S.W. Hawking, in: S.W. Hawking, W. Israel (Eds.), General Relativity. An Einstein Centenary Survey, Cambridge Univ. Press, Cambridge, 1979, p. 746. [23] M. Mueller, M.J. Perry, Class. Quantum Grav. 3 (1986) 65. [24] P. Forgacs, N.S. Manton, Commun. Math. Phys. 72 (1980) 15. [25] P.G. Bergmann, E.J. Flaherty, J. Math. Phys. 19 (1978) 212. [26] J.J. van der Bij, E. Radu, Int. J. Mod. Phys. A 18 (2003) 2379, hep-th/0210185. [27] C.V. Johnson, R.C. Myers, Phys. Rev. D 50 (1994) 6512, hep-th/9406069. [28] S. Ramaswamy, A. Sen, J. Math. Phys. 22 (1981) 2612; A. Magnon, J. Math. Phys. 27 (1986) 1066; A. Magnon, J. Math. Phys. 28 (1987) 2149. [29] U. Ascher, J. Christiansen, R.D. Russell, Math. Comp. 33 (1979) 659; U. Ascher, J. Christiansen, R.D. Russell, ACM Trans. 7 (1981) 209. Y. Brihaye, E. Radu / Physics Letters B 615 (2005) 1–13 13 [30] R.L. Zimmerman, B.Y. Shahir, Gen. Relativ. Gravit. 21 (1989) 821. [31] B. Hartmann, B. Kleihaus, J. Kunz, Phys. Rev. D 65 (2002) 024027, hep-th/0108129. [32] B. Kleihaus, J. Kunz, F. Navarro-Lerida, Phys. Lett. B 599 (2004) 294, gr-qc/0406094. [33] N. Alonso-Alberca, P. Meessen, T. Ortin, Class. Quantum Grav. 17 (2000) 2783, hep-th/0003071. [34] D. Astefanesei, R.B. Mann, E. Radu, JHEP 0501 (2005) 049, hep-th/0407110. [35] C.N. Pope, A.L. Yuille, Phys. Lett. B 78 (1978) 424. Physics Letters B 615 (2005) 14–18 www.elsevier.com/locate/physletb Search for a Lorentz invariance violation contribution in atmospheric neutrino oscillations using MACRO data G. Battistoni a, Y. Becherini b, S. Cecchini b,c, M. Cozzi b,H.Dekhissib,d, L.S. Esposito b,G.Giacomellib, M. Giorgini b,G.Mandriolib,S.Manzoorb,e, A. Margiotta b, L. Patrizii b,V.Popab,f, M. Sioli b,∗,G.Sirrib, M. Spurio b, V. Togo b a INFN Sezione di Milano, I-20133 Milano, Italy b Dipartimento di Fisica dell’Università di Bologna and INFN, I-40127 Bologna, Italy c INAF-IASF Sezione di Bologna, I-40129 Bologna, Italy d Facultè des Sciences, Universitè Mohamed Ier, Oujda, Morocco e RPD, PINSTECH, P.O. Nilore, Islamabad, Pakistan f Institute for Space Sciences, R77125 Bucharest, Romania Received 7 March 2005; received in revised form 5 April 2005; accepted 6 April 2005 Available online 13 April 2005 Editor: L. Rolandi Abstract The energy spectrum of neutrino-induced upward-going muons in MACRO has been analysed in terms of relativity prin- ciples violating effects, keeping standard mass-induced atmospheric neutrino oscillations as the dominant source of νµ → ντ transitions. The data disfavor these exotic possibilities even at a subdominant level, and stringent 90% C.L. limits are placed on −24 −26 the Lorentz invariance violation parameter |v| < 6 × 10 at sin 2θv = 0and|v| < 2.5–5 × 10 at sin 2θv =±1. These limits can also be re-interpreted as upper bounds on the parameters describing violation of the equivalence principle. 2005 Elsevier B.V. All rights reserved. PACS: 14.60.Pq; 14.60.Lm; 11.30.Cp Keywords: Neutrino mass and mixing; Lorentz and Poincaré invariance Neutrino mass-induced oscillations are the best ex- Two flavor νµ → ντ oscillations are strongly favored planation of the atmospheric neutrino problem [1–4]. over a wide range of alternative solutions such as νµ → νsterile oscillations [5,6], νµ → νe oscillations [3,4] or other exotic possibilities [7,8]. * Corresponding author. E-mail addresses: [email protected], In this Letter, we assume standard mass-induced [email protected] (M. Sioli). neutrino oscillations as the leading mechanism for fla- 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.04.010 G. Battistoni et al. / Physics Letters B 615 (2005) 14–18 15 vor transitions and we treat Lorentz invariance flavor In this case, the νµ survival probability is transitions as a subdominant effect [9]. More specifi- P(ν → ν ) cally, we constrain mass-induced neutrino oscillation µ µ 2 2 18 parameters to the ones obtained with a global fit of = 1 − sin 2θv sin 2.54 × 10 vLEν , (4) all MACRO neutrino data. Then we study the en- where v = (vνv − vνv ) is the neutrino MAV differ- ergy spectrum of a limited sample of neutrino-induced 3 2 ence in units of c. Notice that neutrino flavor oscil- upward-going muons and its compatibility with the lations induced by VLI are characterized by an LE inclusion of competitive oscillation scenarios. In the ν dependence of the oscillation probability (Eq. (4)), to literature, neutrino oscillations induced by violation of be compared with the L/E behavior of mass-induced (CPT-conserving) Lorentz invariance (VLI) and viola- ν oscillations (Eq. (2)). tion of the equivalence principle (VEP) are described When both mass-induced transitions and VLI- within the same formalism. In the following we will induced transitions are considered simultaneously, the mention only VLI for simplicity. muon neutrino survival probability can be expressed In this scenario, neutrinos can be described in terms as [9–11] of three distinct bases: flavor eigenstates, mass eigen- 2 2 states and velocity eigenstates, the latter being char- P(νµ → νµ) = 1 − sin 2Θ sin Ω, (5) acterized by different maximum attainable velocities (MAVs) in the limit of infinite momentum. where the global mixing angle Θ and the term Ω are Both mass-induced oscillations and VLI transitions given by: are treated in the two-family approximation and we 2Θ = atan(a1/a2), assume that mass and velocity mixings occur inside | | = 2 + 2 the same families (e.g., ν2 and ν3 ). Ω a1 a2 (6) The usual interpretation of the atmospheric neu- → with trino oscillations is νµ ντ induced by the mixing | m | m 2 of the two mass eigenstates ν2 and ν3 , and two a1 = 1.27 m sin 2θmL/Eν | | weak eigenstates νµ and ντ , i.e., + × 18 iη 2 10 v sin 2θvLEνe , m m m m |ν = ν cos θ + ν sin θ , = 2 µ 2 23 3 23 a2 1.27 m cos 2θmL/Eν m m m m |ν =−ν sin θ + ν cos θ , (1) 18 τ 2 23 3 23 + 2 × 10 v cos 2θvLEν . (7) m where θ (≡ θm) is the flavor-mass mixing angle. The 2 23 Here m , L and Eν are expressed, as in Eqs. (2) survival probability of muon neutrinos at a distance L and (4),ineV2, km and GeV, respectively. The addi- from production is tional factor eiη connects the mass and velocity eigen- states, and for the moment it is assumed to be real P(νµ → νµ) (η = 0orπ). Note that formulae (2) and (5) do not = − 2 2 2 1 sin 2θm sin 1.27m L/Eν , (2) depend on the sign of the mixing angle and/or on 2 2 = 2 − 2 2 the sign of the v and m parameters; this is not where m (mνm mνm ) is expressed in eV , L in 3 2 so in the case of mixed oscillations, where the rela- km and the neutrino energy Eν in GeV. Notice the de- 2 tive sign between the mass-induced and VLI-induced pendence on L/Eν in the argument of the second sin oscillation terms is important. The whole domain of term. variability of the parameters can be accessed with the In the VLI case, the two flavor eigenstates |νµ, |ντ 2 requirements m 0, 0 θm π/2, v 0 and and the two velocity eigenstates |νv, |νv are con- 2 3 −π/4 θv π/4. θ v ≡ θ nected through the mixing angle 23 ( v) in analogy The same formalism also applies to violation of the with mass-induced oscillations: equivalence principle, after substituting v/2 with the v v v v adimensional product |φ|γ ; γ is the difference of |νµ= ν cos θ + ν sin θ , 2 23 3 23 the coupling constants for neutrinos of different types | =− v v + v v ντ ν2 r sin θ23 ν3 cos θ23. (3) to the gravitational potential φ [12]. 16 G. Battistoni et al. / Physics Letters B 615 (2005) 14–18 Fig. 1. Energy dependence of the νµ → νµ survival probability for mass-induced oscillations alone (continuous line) and mass-induced + VLI −25 oscillations for v = 2 × 10 and sin 2θv = 0, ±0.3, ±0.7and±1 (dashed lines for positive values, dotted lines for negative values). 2 2 The neutrino pathlength was fixed at L = 10000 km and we assumed m = 0.0023 eV , θm = π/4. As shown in [10], in a footnoted comment of [13] current interactions νµ + N → µ + X; upgoing muons and more recently in [11], the most sensitive tests of were identified with the streamer tube system (for VLI can be made by analysing the high energy tail tracking) and the scintillator system (for time-of-flight of atmospheric neutrinos at large pathlength values. measurement). Early results concerning atmospheric As an example, Fig. 1 shows the energy dependence neutrinos were published in [15] and in [1] for the up- of the νµ → νµ survival probability as a function of throughgoing muon sample and in [16] for the low the neutrino energy, for neutrino mass-induced oscil- energy semi-contained and upgoing-stopping muon lations alone and for both mass and VLI-induced os- events. Matter effects in the νµ → νsterile channel were cillations for v = 2 × 10−25 and different values of presented in [5] and a global analysis of all MACRO sin 2θv parameter. Note the large sensitivity for large neutrino data in [2]. neutrino energies and large mixing angles. Given the In order to analyse the MACRO data in terms of very small neutrino mass (mν 1 eV), neutrinos with VLI, we used a subsample of 300 upthroughgoing energies larger than 100 GeV are extremely relativis- muons whose energies were estimated via multiple tic, with Lorentz γ factors larger than 1011. Coulomb scattering in the 7 horizontal rock absorbers MACRO [14] was a multipurpose large area detec- in the lower apparatus [17,18]. The energy estimate tor (∼10000 m2 sr acceptance for an isotropic flux) was obtained using the streamer tubes in drift mode, located in the Gran Sasso underground Lab, shielded which allowed to considerably improve the spatial res- by a minimum rock overburden of 3150 hg/cm2.The olution of the detector (∼3 mm). The overall neutrino detector had global dimensions of 76.6 × 12 × 9.3m3 energy resolution was of the order of 100%, mainly and used limited streamer tubes and scintillation coun- dominated by muon energy losses in the rock be- ters to detect muons. νµ’s were detected via charged low the detector (note that Eµ0.4Eν). Upgoing G. Battistoni et al. / Physics Letters B 615 (2005) 14–18 17 muon neutrinos of this sample have large zenith angles (> 120◦) and the median value of neutrino pathlengths is slightly larger than 10000 km. Following the analysis in Ref. [18], we selected a low and a high energy sample by requiring that the rec rec reconstructed neutrino energy Eν should be Eν < rec 30 GeV and Eν > 130 GeV. The number of events surviving these cuts is Nlow = 49 and Nhigh = 58, re- spectively; their median energies, estimated via Monte Carlo, are 13 GeV and 204 GeV (assuming mass- induced oscillations). The analysis then proceeds by fixing the neutrino mass oscillation parameters at the values obtained with the global analysis of all MACRO neutrino data [2]: Fig. 2. 90% C.L. upper limits on the Lorentz invariance viola- 2 2 2 ¯ tion parameter v versus sin 2θv. Standard mass induced oscilla- m = 0.0023 eV ,sin 2θm = 1. Then, we scanned tions are assumed in the two-flavor νµ → ντ approximation, with 2 2 the plane of the two free parameters (v, θv)using m = 0.0023 eV and θm = π/4. The dashed line shows the limit the function obtained with the same selection criteria of Ref. [18] to define the low and high energy samples; the continuous line is the final result obtained with the selection criteria optimized for the present analy- high MC 2 ¯ 2 Ni − αN (v, θv; m , θm) sis (see text). χ2 = i , (8) σi i=low The energy cuts described above (the same used in Ref. [18]), were optimized for mass-induced neutrino where N MC is the number of events predicted by i oscillations. In order to maximize the sensitivity of the Monte Carlo, α is a constant which normalizes the analysis for VLI induced oscillations, we performed number of Monte Carlo events to the number of ob- a blind analysis, based only on Monte Carlo events, served events and σ is the overall error comprehensive i to determine the energy cuts which yield the best per- of statistical and systematic uncertainties. formances. The results of this study suggest the cuts We used the Monte Carlo simulation described Erec < 28 GeV and Erec > 142 GeV; with these cuts in [18] with different neutrino fluxes in input [19–22]. ν ν the number of events in the real data are N = 44 The largest relative difference of the extreme values low events and N = 35 events. The limits obtained with of the MC expected ratio N /N is 13%. How- high low high this selection are shown in Fig. 2 by the continuous ever, in the evaluation of the systematic error, the main line. As expected, the limits are now more stringent sources of uncertainties for this ratio (namely the pri- than for the previous choice. mary cosmic ray spectral index and neutrino cross In order to understand the dependence of this re- sections) have been separately estimated and their ef- 2 fects added in quadrature (see [18] for details): in this sult with respect to the choice of the m parameter, 2 work, we use a conservative 16% theoretical system- we varied the m values around the best-fit point. 2 ± atic error on the ratio Nlow/Nhigh. The experimental We found that a variation of m of 30% moves systematic error on the ratio was estimated to be 6%. up/down the upper limit of VLI parameters by at most In the following, we show the results obtained with the a factor 2. computation in [22]. Finally, we computed the limit on v marginal- The inclusion of the VLI effect does not improve ized with respect to all the other parameters left free 2 2 2 the χ in any point of the (v, θv) plane, compared to to variate inside the intervals: m = m ± 30%, ¯ mass-induced oscillations stand-alone, and proper up- θm = θm ± 20%, −π/4 θv π/4 and any value of per limits on VLI parameters were obtained. The 90% the phase η. We obtained the 90% C.L. upper limit −25 C.L. limits on v and θv, computed with the Feld- |v| < 3 × 10 . man and Cousins prescription [23],areshownbythe An independent and complementary analysis was dashed line in Fig. 2. performed on a sample of events with a reconstructed 18 G. Battistoni et al. / Physics Letters B 615 (2005) 14–18 rec neutrino energy 25 GeV Measurement of exclusive ρ+ρ− production in mid-virtuality two-photon interactions and study of the γγ∗ → ρρ process at LEP L3 Collaboration P. Achard t, O. Adriani q, M. Aguilar-Benitez y, J. Alcaraz y, G. Alemanni w, J. Allaby r, A. Aloisio ac, M.G. Alviggi ac, H. Anderhub aw, V.P. Andreev f,ah,F.Anselmoh, A. Arefiev ab, T. Azemoon c, T. Aziz i,P.Bagnaiaam,A.Bajoy,G.Baksayz,L.Baksayz, S.V. Baldew b,S.Banerjeei,Sw.Banerjeed, A. Barczyk aw,au, R. Barillère r, P. Bartalini w, M. Basile h,N.Batalovaat, R. Battiston ag,A.Bayw, F. Becattini q, U. Becker m,F.Behneraw, L. Bellucci q, R. Berbeco c, J. Berdugo y,P.Bergesm, B. Bertucci ag,B.L.Betevaw,M.Biasiniag, M. Biglietti ac,A.Bilandaw, J.J. Blaising d, S.C. Blyth ai, G.J. Bobbink b,A.Böhma, L. Boldizsar l,B.Borgiaam,S.Bottaiq, D. Bourilkov aw, M. Bourquin t, S. Braccini t,J.G.Bransonao,F.Brochud,J.D.Burgerm, W.J. Burger ag,X.D.Caim, M. Capell m, G. Cara Romeo h, G. Carlino ac, A. Cartacci q, J. Casaus y, F. Cavallari am, N. Cavallo aj, C. Cecchi ag, M. Cerrada y,M.Chamizot, Y.H. Chang ar,M.Chemarinx,A.Chenar,G.Cheng,G.M.Cheng,H.F.Chenv, H.S. Chen g, G. Chiefari ac, L. Cifarelli an, F. Cindolo h,I.Clarem,R.Clareal, G. Coignet d,N.Colinoy, S. Costantini am,B.delaCruzy, S. Cucciarelli ag, R. de Asmundis ac, P. Déglon t, J. Debreczeni l,A.Degréd,K.Dehmeltz, K. Deiters au, D. della Volpe ac, E. Delmeire t, P. Denes ak, F. DeNotaristefani am,A.DeSalvoaw, M. Diemoz am, M. Dierckxsens b, C. Dionisi am, M. Dittmar aw,A.Doriaac,M.T.Dovaj,5, D. Duchesneau d, M. Duda a, B. Echenard t,A.Eliner,A.ElHagea, H. El Mamouni x, A. Engler ai,F.J.Epplingm, P. Extermann t, M.A. Falagan y, S. Falciano am,A.Favaraaf, J. Fay x,O.Fedinah, M. Felcini aw, T. Ferguson ai, H. Fesefeldt a, E. Fiandrini ag, J.H. Field t, F. Filthaut ae,P.H.Fisherm, W. Fisher ak, I. Fisk ao, G. Forconi m, K. Freudenreich aw,C.Furettaaa, Yu. Galaktionov ab,m, S.N. Ganguli i, P. Garcia-Abia y, M. Gataullin af, S. Gentile am,S.Giaguam, Z.F. Gong v,G.Grenierx,O.Grimmaw, M.W. Gruenewald p, M. Guida an, V.K. Gupta ak,A.Gurtui,L.J.Gutayat, D. Haas e, D. Hatzifotiadou h, T. Hebbeker a,A.Hervér,J.Hirschfelderai,H.Hoferaw, 0370-2693/$ – see front matter 2005 Published by Elsevier B.V. doi:10.1016/j.physletb.2005.04.011 20 L3 Collaboration / Physics Letters B 615 (2005) 19–30 M. Hohlmann z,G.Holzneraw,S.R.Houar,B.N.Jing, P. Jindal n, L.W. Jones c, P. de Jong b, I. Josa-Mutuberría y, M. Kaur n, M.N. Kienzle-Focacci t,J.K.Kimaq, J. Kirkby r, W. Kittel ae, A. Klimentov m,ab, A.C. König ae, M. Kopal at, V. Koutsenko m,ab, M. Kräber aw, R.W. Kraemer ai,A.Krügerav, A. Kunin m, P. Ladron de Guevara y, I. Laktineh x, G. Landi q, M. Lebeau r,A.Lebedevm,P.Lebrunx,P.Lecomteaw, P. Lecoq r,P.LeCoultreaw,J.M.LeGoffr,R.Leisteav,M.Levtchenkoaa, P. Levtchenko ah,C.Liv, S. Likhoded av,C.H.Linar,W.T.Linar,F.L.Lindeb, L. Lista ac, Z.A. Liu g, W. Lohmann av, E. Longo am,Y.S.Lug,C.Luciam, L. Luminari am, W. Lustermann aw,W.G.Mav,L.Malgerir, A. Malinin ab,C.Mañay,J.Mansak, J.P. Martin x, F. Marzano am, K. Mazumdar i, R.R. McNeil f,S.Meler,ac,L.Merolaac, M. Meschini q, W.J. Metzger ae,A.Mihulk,H.Milcentr, G. Mirabelli am,J.Mnicha, G.B. Mohanty i,G.S.Muanzax, A.J.M. Muijs b,B.Musicarao,M.Musyam, S. Nagy o, S. Natale t, M. Napolitano ac, F. Nessi-Tedaldi aw,H.Newmanaf,A.Nisatiam, T. Novak ae,H.Nowakav, R. Ofierzynski aw,G.Organtiniam,I.Palat,C.Palomaresy, P. Paolucci ac, R. Paramatti am,G.Passalevaq,S.Patricelliac,T.Paulj, M. Pauluzzi ag, C. Paus m,F.Paussaw, M. Pedace am,S.Pensottiaa, D. Perret-Gallix d,D.Piccoloac, F. Pierella h, M. Pioppi ag,P.A.Pirouéak, E. Pistolesi aa, V. Plyaskin ab, M. Pohl t, V. Pojidaev q, J. Pothier r,D.Prokofievah, G. Rahal-Callot aw, M.A. Rahaman i, P. Raics o,N.Rajai, R. Ramelli aw,P.G.Rancoitaaa,R.Ranieriq, A. Raspereza av, P. Razis ad,D.Renaw, M. Rescigno am, S. Reucroft j, S. Riemann av,K.Rilesc,B.P.Roec, L. Romero y,A.Roscaav, C. Rosemann a, C. Rosenbleck a,S.Rosier-Leesd,S.Rotha, J.A. Rubio r, G. Ruggiero q, H. Rykaczewski aw,A.Sakharovaw, S. Saremi f,S.Sarkaram, J. Salicio r, E. Sanchez y, C. Schäfer r, V. Schegelsky ah, H. Schopper u, D.J. Schotanus ae, C. Sciacca ac,L.Servoliag, S. Shevchenko af,N.Shivarovap, V. Shoutko m, E. Shumilov ab, A. Shvorob af,D.Sonaq,C.Sougax, P. Spillantini q, M. Steuer m, D.P. Stickland ak,B.Stoyanovap,A.Straessnert, K. Sudhakar i, G. Sultanov ap, L.Z. Sun v,S.Sushkova,H.Suteraw,J.D.Swainj, Z. Szillasi z,3,X.W.Tangg,P.Tarjano, L. Tauscher e, L. Taylor j, B. Tellili x, D. Teyssier x,C.Timmermansae, Samuel C.C. Ting m, S.M. Ting m, S.C. Tonwar i,J.Tóthl,C.Tullyak,K.L.Tungg, J. Ulbricht aw, E. Valente am,R.T.VandeWalleae, R. Vasquez at, V. Veszpremi z, G. Vesztergombi l, I. Vetlitsky ab,G.Viertelaw, S. Villa al, M. Vivargent d,S.Vlachose, I. Vodopianov z,H.Vogelai,H.Vogtav, I. Vorobiev ai,ab, A.A. Vorobyov ah, M. Wadhwa e, Q. Wang ae,X.L.Wangv,Z.M.Wangv, M. Weber r, S. Wynhoff ak,L.Xiaaf,Z.Z.Xuv, J. Yamamoto c,B.Z.Yangv,C.G.Yangg,H.J.Yangc,M.Yangg,S.C.Yehas, An. Zalite ah, Yu. Zalite ah, Z.P. Zhang v,J.Zhaov,G.Y.Zhug,R.Y.Zhuaf, H.L. Zhuang g,A.Zichichih,r,s,B.Zimmermannaw, M. Zöller a a III Physikalisches Institut, RWTH, D-52056 Aachen, Germany 1 b National Institute for High Energy Physics, NIKHEF, and University of Amsterdam, NL-1009 DB Amsterdam, The Netherlands L3 Collaboration / Physics Letters B 615 (2005) 19–30 21 c University of Michigan, Ann Arbor, MI 48109, USA d Laboratoire d’Annecy-le-Vieux de Physique des Particules, LAPP, IN2P3-CNRS, BP 110, F-74941 Annecy-le-Vieux cedex, France e Institute of Physics, University of Basel, CH-4056 Basel, Switzerland f Louisiana State University, Baton Rouge, LA 70803, USA g Institute of High Energy Physics, IHEP, 100039 Beijing, China 6 h University of Bologna, and INFN, Sezione di Bologna, I-40126 Bologna, Italy i Tata Institute of Fundamental Research, Mumbai (Bombay) 400 005, India j Northeastern University, Boston, MA 02115, USA k Institute of Atomic Physics and University of Bucharest, R-76900 Bucharest, Romania l Central Research Institute for Physics of the Hungarian Academy of Sciences, H-1525 Budapest 114, Hungary 2 m Massachusetts Institute of Technology, Cambridge, MA 02139, USA n Panjab University, Chandigarh 160 014, India o KLTE-ATOMKI, H-4010 Debrecen, Hungary 3 p Department of Experimental Physics, University College Dublin, Belfield, Dublin 4, Ireland q INFN, Sezione di Firenze, and University of Florence, I-50125 Florence, Italy r European Laboratory for Particle Physics, CERN, CH-1211 Geneva 23, Switzerland s World Laboratory, FBLJA Project, CH-1211 Geneva 23, Switzerland t University of Geneva, CH-1211 Geneva 4, Switzerland u University of Hamburg, D-22761 Hamburg, Germany v Chinese University of Science and Technology, USTC, Hefei, Anhui 230 029, China 6 w University of Lausanne, CH-1015 Lausanne, Switzerland x Institut de Physique Nucléaire de Lyon, IN2P3-CNRS, Université Claude Bernard, F-69622 Villeurbanne, France y Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas, CIEMAT, E-28040 Madrid, Spain 4 z Florida Institute of Technology, Melbourne, FL 32901, USA aa INFN, Sezione di Milano, I-20133 Milan, Italy ab Institute of Theoretical and Experimental Physics, ITEP, Moscow, Russia ac INFN, Sezione di Napoli, and University of Naples, I-80125 Naples, Italy ad Department of Physics, University of Cyprus, Nicosia, Cyprus ae Radboud University and NIKHEF, NL-6525 ED Nijmegen, The Netherlands af California Institute of Technology, Pasadena, CA 91125, USA ag INFN, Sezione di Perugia, and Università Degli Studi di Perugia, I-06100 Perugia, Italy ah Nuclear Physics Institute, St. Petersburg, Russia ai Carnegie Mellon University, Pittsburgh, PA 15213, USA aj INFN, Sezione di Napoli, and University of Potenza, I-85100 Potenza, Italy ak Princeton University, Princeton, NJ 08544, USA al University of California, Riverside, CA 92521, USA am INFN, Sezione di Roma, and University of Rome “La Sapienza”, I-00185 Rome, Italy an University and INFN, Salerno, I-84100 Salerno, Italy ao University of California, San Diego, CA 92093, USA ap Bulgarian Academy of Sciences, Central Laboratory of Mechatronics and Instrumentation, BU-1113 Sofia, Bulgaria aq The Center for High Energy Physics, Kyungpook National University, 702-701 Taegu, South Korea ar National Central University, Chung-Li, Taiwan as Department of Physics, National Tsing Hua University, Taiwan at Purdue University, West Lafayette, IN 47907, USA au Paul Scherrer Institut, PSI, CH-5232 Villigen, Switzerland av DESY, D-15738 Zeuthen, Germany aw Eidgenössische Technische Hochschule, ETH Zürich, CH-8093 Zürich, Switzerland Received 18 March 2005; received in revised form 5 April 2005; accepted 6 April 2005 Available online 14 April 2005 Editor: L. Rolandi 22 L3 Collaboration / Physics Letters B 615 (2005) 19–30 Abstract + − ∗ Exclusive ρ ρ production in two-photon collisions between a quasi-real√ photon, γ , and a mid-virtuality photon, γ ,is studied with data collected at LEP at centre-of-mass energies 183 GeV s 209 GeV with a total integrated luminosity of − ∗ + − 684.8pb 1. The cross section of the γγ → ρ ρ process is determined as a function of the photon virtuality, Q2,andthe 2 2 2 two-photon centre-of-mass energy, Wγγ, in the kinematic region: 0.2GeV Q 0.85 GeV and 1.1GeV Wγγ 3GeV. + − ∗ These results, together with previous L3 measurements of ρ0ρ0 and ρ ρ production, allow a study of the γγ → ρρ process over the Q2-region 0.2GeV2 Q2 30 GeV2. 2005 Published by Elsevier B.V. 1. Introduction duction were found to have a similar dependence on Wγγ and to be of similar magnitude. However, the The L3 Collaboration has recently measured the ex- ρ+ρ− cross section is systematically higher than the clusive production of ρ0ρ0 [1,2] and ρ+ρ− [3] pairs ρ0ρ0 one. This is in contrast with the suppression + − in the two-photon fusion process: and different Wγγ dependence of ρ ρ production 0 0 + − + − ∗ + − [6] with respect to ρ ρ [7] observed in data with e e → e e γγ → e e ρρ, (1) 2 Q ≈ 0 and Wγγ 2 GeV. We note that despite the where the beam electrons7 radiate virtual photons wide range of theoretical models [8,9], ρ-pair produc- which interact and produce a hadronic final state. One tion at Q2 ≈ 0 is still not well understood. Therefore of the photons, γ , is quasi-real, characterised by a the experimental study of the Q2-evolution of ρ-pair 2 = small value of its squared four momentum, Pγ production is important to understand vector meson m2 ≈ 0, whereas the other one, γ ∗, has a significant pair-production in two-photon interactions. γ 0 0 2 =− 2 =− 2 − 2 Previously, we performed a measurement of ρ ρ virtuality, Q Pγ ∗ mγ ∗ mγ . Our mea- surements cover the two-photon centre-of-mass en- production [2] for intermediate virtualities: ergy 0.2GeV2 Q2 0.85 GeV2. (3) 1.1GeV W 3GeV. (2) γγ In this Letter, we complement that study with the first The two measurements [1,3] done at large virtual- measurement of the process ities, 1.2GeV2 Q2 30 GeV2, provide a testing + − + − ∗ + − + − ground for a recently-developed QCD-based model e e → e e γγ → e e ρ ρ (4) [4]. This model describes well the Q2-dependence of in the kinematic region (2) and (3). These data allow the ρ0ρ0 production at large momentum transfer [5]. + − to follow the Q2-evolution of the ρρ-production over The measured cross sections for ρ0ρ0 and ρ ρ pro- two orders of magnitude in this variable. The analysis techniques employed in this study are 1 Supported by the German Bundesministerium für Bildung, similar to those of our previous measurements [2,3]. Wissenschaft, Forschung und Technologie. The data used, corresponding to an integrated lumi- 2 Supported by the Hungarian OTKA fund under contract num- nosity of 684.8 pb−1,arethesameasinRef.[2] bers T019181, F023259 and T037350. 3 Also supported by the Hungarian OTKA fund under contract and were collected by the L3 detector√[10] at LEP at number T026178. centre-of-mass energies 183 GeV s 209 GeV. 4 Supported also by the Comisión Interministerial de Ciencia y Scattered beam electrons which have radiated photons Tecnología. with virtualities in the range (3) can be “tagged” by 5 Also supported by CONICET and Universidad Nacional de La the Very Small Angle Tagger (VSAT) [11].TheVSAT Plata, CC 67, 1900 La Plata, Argentina. is an electromagnetic calorimeter, constructed with 6 Supported by the National Natural Science Foundation of China. BGO crystals, with a geometrical acceptance covering 7 Throughout this Letter, the term “electron” denotes both elec- the polar angle range 5 mrad θ 10 mrad, for az- tron and positron. imuthal angles in the ranges −1.25 rad φ 1.25 rad L3 Collaboration / Physics Letters B 615 (2005) 19–30 23 and π − 1.25 rad φ π + 1.25 rad. When the elec- the fits is taken. To make the selection robust against tron with the largest scattering angle is detected in instrumental noise and backgrounds and to reduce the the VSAT, the virtuality of the photon it radiated is, sensitivity to the Monte Carlo simulation of fake pho- within 1% precision, equal to the transverse momen- tons, we retain events with one additional photon, not 2 0 0 tum squared, pt , of the final state hadron system: used in the π π pair, if the photon energy is less than 300 MeV and does not exceed 10% of the energy of 2 = − ≈ 2 ≈ 2 Q 2EbEs(1 cos θs) EbEsθs pt , (5) the π 0π 0 pair. 2 where Eb is the beam energy, and Es and θs are the The transverse momentum squared, pt , of the four- energy and the scattering angle of the tagged electron. pion system is used to measure the Q2 of the event Therefore the VSAT is not used to directly measure and is required to be in the range 0.2–0.85 GeV2.For Q2, but rather to select exclusive final states by corre- selection of an exclusive final state, the acoplanarity lating the direction of the transverse momentum vector angle, φaco, calculated from the difference between the of the tagged electron with the detected hadron sys- azimuthal angle of the tagged electron, φtag,shown tem. in Fig. 1(b), and the azimuthal angle of the four-pion system, is required to be less than 0.4 rad, as shown in Fig. 1(c). The data contain a contribution from η 2. Event selection production, as visible in the π +π −π 0 mass spectrum, shown in Fig. 1(d). This background is removed by re- The reaction e+e− → e+e−ρ+ρ− contributing to quiring M(π+π −π 0)>0.65 GeV. the process After all cuts, 414 events are retained. Their four- + − → + − + − 0 0 pion mass spectrum is shown in Fig. 2(a). The region e e e etagπ π π π (6) 1.1GeV Wγγ 3 GeV is populated by 387 events, is identified by one and only one scattered electron, which are used for the cross section determination. A ± etag, detected in the VSAT, two charged pions mea- strong signal from ρ production is observed in the sured in the tracking chamber, and energy clusters M(π±π 0) spectrum, shown in Fig. 2(b). The cluster- from the two-photon decays of the π 0’s, deposited ing of entries at the crossing of the ρ± mass bands in in the BGO electromagnetic calorimeter. These events the correlation plot of the masses of the π ±π 0 com- are collected by two independent track-triggers [12]. binations, shown in Fig. 2(c), gives evidence for a The trigger efficiency, as determined from the data it- signal from ρ+ρ− intermediate states. No structure is self, is (60 ± 3)%. observed in the correlation plot of the masses of the Single-tagged events are selected by requiring just π +π − and π 0π 0 combinations, shown in Fig. 2(d). one electromagnetic cluster with energy greater then We also inspected the two- and three-pion mass distri- 50% of the beam energy reconstructed in the VSAT. butions, shown in Fig. 3, for production of higher-mass The event candidates must have exactly two tracks resonances. The only statistically-significant signal is ± ± 0 0 with zero total charge. The tracks must come from the from the a2 (1320) state in the π π π mass spec- interaction vertex, have transverse momentum greater trum, as seen in Fig. 3(f). than 100 MeV and an energy loss in the tracking chamber compatible with the pion hypothesis. The se- lected events should contain a π 0π 0 pair, therefore we 3. Data analysis consider event candidates that have four or five pho- tons, identified as isolated clusters in the electromag- 3.1. Monte Carlo modelling netic calorimeter. Photons having energies greater than 60 MeV are paired to reconstruct neutral pions, which To estimate the number of ρ+ρ− events in the are required to be in the mass window 100 MeV selected four-pion data sample, we consider non- M(γγ) 170 MeV, as shown in Fig. 1(a). The mass interfering contributions from the processes: of a π 0 candidate is constrained to the nominal value ∗ → + − by a 1-C kinematic fit. If more than one π 0π 0 com- γγ ρ ρ , ∗ ± ∓ bination exists, the one with the smallest χ2 sum of γγ → ρ π π 0, 24 L3 Collaboration / Physics Letters B 615 (2005) 19–30 + − Fig. 1. Distributions for π π π0π0 candidates. (a) Two-photon invariant mass for the selected π0’s (two entries per event); (b) azimuthal angle, φtag, of the tagged electron for tags in the inner side of the LEP ring (in) and, folded over it, for tags in the outer side of the LEP ring + − 0 0 + − 0 (out); (c) acoplanarity angle, φaco, between the tagged electron and the π π π π system and (d) mass of the π π π system (two entries per event). The data are compared to the four-pion Monte Carlo. The estimated background is indicated by the hatched histograms. The arrows indicate the selection cuts. ∗ → ± ∓ γγ a2 (1320)π , 3.2. Background estimation ∗ + − γγ → π π π 0π 0, non-resonant. (7) The contribution to the selected events from e+e− About 40 million Monte Carlo events of the processes annihilation and from the process e+e− → e+e−τ +τ − (7) are generated with the EGPC [13] program, which is negligible. Random coincidences with off-momen- uses the luminosity function from Ref. [14]. Particle tum beam electrons, which give signals in the VSAT, production and decay is uniform in phase-space. The are a source of background. The flux of these par- generated events are passed through the full L3 detec- ticles is dominantly on the outer side of the LEP tor simulation using the GEANT [15] and GHEISHA ring. Therefore, this background would cause an ex- [16] programs and processed in the same way as the cess in the number of events having a tag on the data, reproducing the detector behaviour as monitored outer side of the accelerator ring, Nout, with respect in the different data-taking periods. to the inner side, Nin. In the selected data, the ratio For acceptance calculations, Monte Carlo events Nout/Nin = 1.04 ± 0.10 is close to unity, indicating are assigned a Q2-dependent weight, evaluated us- that this background is small. This conclusion is cor- ing the GVDM form-factor [17] for both interacting roborated by the good agreement observed between photons. The detection efficiencies of the process (4) the φtag distribution of the selected data and Monte 2 are listed in Tables 1 and 2 for bins in Q and Wγγ. Carlo event samples, shown in Fig. 1(b). The efficiencies for the four-pion final states of all the Two sources of background remain. The first is processes (7) are of similar magnitude. partially-reconstructed events from two-photon inter- L3 Collaboration / Physics Letters B 615 (2005) 19–30 25 ± 0 Fig. 2. Mass distributions for the selected events: (a) the four-pion system, Wγγ;(b)theπ π combinations (four entries per event); (c) corre- − + + − lation between the π π0 and π π0 pairs (two entries per event) and (d) correlation between the π π and π0π0 pairs. The two-dimensional distributions have a bin width of 55 × 55 MeV2, the size of the boxes is proportional to the number of entries and both plots have the same vertical scale. actions with higher particle multiplicities, when tracks passing the selection, are combined with the distribu- or photons escape detection. The second is signal tion of selected π +π −π 0π 0 Monte Carlo events so as events with one or more photons substituted by pho- to reproduce the φaco distribution observed in the data, ton candidates due to noise. To estimate the accepted as shown in Fig. 1(c). The estimated background lev- background we use background-like event samples ex- els are listed in Tables 1 and 2. As data samples are tracted from the experimental data. The first back- used in the background estimation, they contain also ground is modelled with selected π ±π ±π 0π 0 events, a fraction of events with fake tags and thus take into in which at least two charged particles have not been account the effect of this background. detected and by π +π −π 0π 0π 0 events in which one π 0 is excluded from consideration. An event-mixing 3.3. Fit method technique is employed in order to reproduce events from the second background: one or two photons In order to determine the differential ρ+ρ− pro- forming a π 0 are excluded from a selected event and duction rate, a maximum likelihood fit of the data replaced by photons from another data event. The to a sum of Monte Carlo samples of the processes 2 φaco distributions of the background-like data samples, (7) is performed in intervals of Q and Wγγ using 26 L3 Collaboration / Physics Letters B 615 (2005) 19–30 ± Fig. 3. (a), (c), (e), (g) Mass distributions of the π π0 combinations (four entries per event) in four Q2-intervals. Distributions for the entire 2 2 2 + − 0 0 kinematic region 1.1GeV Wγγ 3 GeV and 0.2GeV Q 0.85 GeV of: (b) The sum of the π π and π π mass spectra (two entries per event). (d) The neutral three-pion combinations (two entries per event). (f) The charged three-pion combinations (two entries per + − ± event). (h) The sum of the π π π0 and π π0π0 mass spectra (four entries per event). The points represent the data, the hatched areas + − show the ρ ρ component and the open areas show the sum of the other contributing processes. The fraction of the different components are determined by the fit and the total normalisation is to the number of the events. a box method [1–3,18]. The inputs to the fit are the six The analysis procedure is optimised for deriving two-pion masses in an event, namely, the four combi- the ρ+ρ− contribution and only the ρ+ρ− content and nations π ±π 0 and the two combinations π +π − and the sum of the rest of the contributing processes, de- π 0π 0. They provide a complete description of a four- noted as “other 4π”, are considered for cross section pion event in our model of isotropic production and measurements. To check the quality of the fit, the two- phase space decay. and three-pion mass distributions of the data are com- L3 Collaboration / Physics Letters B 615 (2005) 19–30 27 Table 1 + − + − + − ∗ + − Detection efficiencies, ε, background fractions, Bg, and cross sections of the reactions e e → e e ρ ρ , γγ → ρ ρ and of the sum 2 of the rest of the contributing processes, “other 4π”, as a function of Q for 1.1GeV Wγγ 3 GeV. The values of the differential cross sections are corrected to the centre of each bin. The first uncertainties are statistical, the second systematic. An overall normalization uncertainty of 5% for the trigger is not included 2 2 2 Q range ε Bg σee [pb] dσee/dQ [pb/GeV ] σγγ [nb] σγγ [nb] + − + − + − [GeV2] [%] [%] ρ ρ ρ ρ ρ ρ other 4π 0.20–0.28 0.8 14 7.4 ± 2.4 ± 1.992± 29 ± 23 5.7 ± 1.8 ± 1.410.9 ± 2.2 ± 1.5 0.28–0.40 1.2 14 5.7 ± 1.8 ± 1.347± 15 ± 10 4.3 ± 1.4 ± 1.012.2 ± 1.8 ± 1.4 0.40–0.55 1.1 15 5.6 ± 1.6 ± 1.137± 11 ± 7.34.9 ± 1.4 ± 1.013.3 ± 2.0 ± 1.8 0.55–0.85 0.7 18 7.7 ± 2.5 ± 2.025± 8.2 ± 6.55.3 ± 1.7 ± 1.412.1 ± 2.2 ± 1.9 Table 2 + − + − + − ∗ + − Detection efficiencies, ε, background fractions, Bg, and cross sections of the reactions e e → e e ρ ρ , γγ → ρ ρ andofthesumof 2 2 2 the rest of the contributing processes, “other 4π”, as a function of Wγγ for 0.2GeV Q 0.85 GeV . The first uncertainties are statistical, the second systematic. An overall normalization uncertainty of 5% for the trigger is not included Wγγ-range ε Bg σee [pb] σγγ [nb] σγγ [nb] + − + − [GeV] [%] [%] ρ ρ ρ ρ other 4π 1.10–1.40 0.6 25 4.9 ± 1.8 ± 1.33.9 ± 1.5 ± 1.19.0±2.4±1.7 1.40–1.65 0.9 18 6.7 ± 1.6 ± 1.37.6 ± 1.9 ± 1.514.8±2.7±2.5 1.65–1.85 1.1 15 5.1 ± 1.5 ± 0.98.4 ± 2.4 ± 1.615.8±3.1±2.3 1.85–2.10 1.1 13 3.9 ± 1.4 ± 0.85.9 ± 2.0 ± 1.218.3±3.0±2.7 2.10–2.40 1.2 10 2.2 ± 1.0 ± 0.53.2 ± 1.4 ± 0.811.5±2.1±1.8 2.40–3.00 1.2 10 2.2 ± 1.0 ± 0.51.9 ± 0.9 ± 0.58.5±1.5±1.5 pared in Fig. 3 with those of a mixture of Monte Carlo and Wγγ bin using the program GALUGA [20], which event samples from the processes (7), in proportions performs O(α4) QED calculations. The same proce- determined by the fit. The observed experimental dis- dure was used in our previous studies [1–3]. The cross tributions are reasonably well described by the Monte section σγγ is derived from the measured cross sec- Carlo model. tion using the relation σγγ = σee/LTT. Thus, σγγ represents an effective cross section containing contri- butions from both transverse and longitudinal photon 4. Results polarisations. The cross section of the process γγ∗ → ρ+ρ− is listed in Table 1 as a function of Q2 and in The cross sections of the process Table 2 as a function of Wγγ. The sum of the cross sec- tions of the other contributing processes is also given + − → + − + − e e e e ρ ρ in Tables 1 and 2. 2 in bins of Q and Wγγ, σee, are listed in Tables 1 Several sources of systematic uncertainty are con- and 2. The statistical uncertainties, also listed in Ta- sidered. The contribution of the selection procedure bles 1 and 2, are those of the fit. The differential cross is in the range 12–18%; Monte Carlo statistics in the 2 section, dσee/dQ , derived from σee, is listed in Ta- range 1.3–2.1%; the fit procedure in the range 11– ble 1. When evaluating the differential cross section, a 20%. Half of the changes of the acceptance when no correction based on the Q2-dependence of the ρ+ρ− form factor re-weighting of the Monte Carlo events is Monte Carlo sample is applied, so as to assign the performed is considered as model uncertainty. It is in cross section value to the centre of the corresponding the range 0.5–5%. The background correction proce- Q2-bin [19]. dure introduces systematic uncertainties in the range To evaluate the cross section, σγγ, of the process 2–6%. All contributions are added in quadrature to ob- γγ∗ → ρ+ρ−, the integral of the transverse photon tain the systematic uncertainties, quoted in Tables 1 2 luminosity function, LTT, is computed for each Q and 2. Finally, a normalization uncertainty of 5% ac- 28 L3 Collaboration / Physics Letters B 615 (2005) 19–30 counts for the uncertainty of the trigger efficiency de- 0.85 GeV2 there is a clear enhancement of ρ0ρ0 pro- termination. duction relative to ρ+ρ−.Thisissimilartowhatwas observed at Q2 ≈ 0 [6,7], but in contrast with the high Q2-region, where both cross sections have similar de- ∗ + − 5. Study of γγ → ρρ process pendence on Wγγ and the ρ ρ is systematically higher than the ρ0ρ0. These differences are clearly Combining the present results with the L3 data seen in the ratio on ρρ production from Refs. [1–3], we compare the + − + − 0 0 σ (ρ ρ ) ρ ρ to the ρ ρ channels and their evolution as R = ee , 0 0 afunctionofQ2. The cross section of the process σee(ρ ρ ) γγ∗ → ρρ is plotted in Fig. 4 as a function of where the sum is for the region 1.1GeV Wγγ W .ForW 2.1 GeV and 0.2GeV2 Q2 γγ γγ 2.1 GeV. In the domain 0.20 GeV2 Q2 0.85 GeV2 we measure R = 0.62 ± 0.10(stat.) ± 0.09(syst.),a value that can only be explained by the presence of an isospin I = 2 intermediate state or by a mix- ture of different states [8,9]. The value of this ra- tio for 1.2GeV2 Q2 8.5GeV2 is R = 1.81 ± 0.47(stat.) ± 0.22(syst.) [3], close to the factor 2, ex- pected for an isospin I = 0 state. Such variation sug- gests different ρ-pair production mechanisms at low and high Q2. 2 The differential cross section dσee/dQ of the re- action e+e− → e+e−ρρ is shown in Fig. 5(a). The L3 measurements span a Q2-region of two orders of mag- nitude, over which the differential cross sections show a monotonic fall of more than four orders of magni- tude. The ρρ data are fitted to a form [21] expected from QCD-based calculations [22]: dσ 1 ee ∼ , 2 n 2 2 2 (8) dQ Q (Q +Wγγ ) where n is a constant and Wγγ is the average Wγγ value, 1.9 GeV for this measurement. Although this 2 formula is expected to be valid only for Q Wγγ, we find it provides a good parametrisation of the Q2- evolution of the ρρ data. A fit to the ρ+ρ− data finds an exponent n = 2.3 ± 0.2 with χ2/d.o.f. = 1.4/7. A value n = 2.9 ± 0.1 was found for ρ0ρ0 with χ2/d.o.f. = 6.9/10 [2]. Only the statistical uncertain- ties are considered in the fits. The results of the fits are shown in Fig. 5(a). The fits indicate a cross-over of the differential cross sections of ρ+ρ− and ρ0ρ0 produc- ∗ Fig. 4. Cross section of the process γγ → ρρ as a function of tion in the vicinity of Q2 ≈ 1GeV2. 2 Wγγ in three Q regions. The results from this measurement, full The measured cross section of the process γγ∗ → points in (a), are compared to previous L3 measurements of the ρρ 2 production [1–3]. The bars show the statistical uncertainties. Some ρρ as a function of Q is shown in Fig. 5(a). The + − 0 0 points from the previous measurements are displaced horizontally change of the relative magnitude of ρ ρ and ρ ρ for better readability. production is clearly visible when comparing the low- L3 Collaboration / Physics Letters B 615 (2005) 19–30 29 2 Fig. 5. The ρρ production cross section as a function of Q ,for1.1GeV Wγγ 3 GeV: (a) differential cross section of the process + − + − ∗ e e → e e ρρ and (b) cross section of the process γγ → ρρ. The results from this measurement, full points in the region Q2 < 1GeV2, are presented together with previous L3 measurements of the ρρ production [1–3]. The bars indicate the statistical uncertainties. Some points from the previous measurements are displaced horizontally for better readability. The lines in (a) represent the results of fits using the QCD-inspired form of Eq. (9). The lines in (b) represent the results of a fit to the ρ0ρ0 data based on the GVDM model [17] and of a fit basedonaρ-pole parametrisation. and the high-Q2 regions. A parametrisation, based on tion and allow to follow the evolution of ρρ cross the GVDM model [17]: sections over a Q2-region of two orders of magnitude. 2 2 A QCD-based form, derived for the description of σγγ Wγγ,Q = σγγ(Wγγ)F Q , 2 the differential cross section dσee/dQ of the process + 2 2 + − → + − 2 2 1 Q /4mV 0.22 e e e e ρρ at high Q , is found to provide a F Q = rV + , 2 (1 + Q2/m2 )2 1 + Q2/m2 good parametrisation of the Q -evolution of the ρρ V =ρ,ω,φ V 0 2 2 2 (9) data in the entire interval 0.2GeV Q 30 GeV , over which the differential cross sections show a with r = 0.65, r = 0.08, r = 0.05 and m = ρ ω φ 0 monotonic decrease of more than four orders of mag- 1.4 GeV reproduces well the Q2-dependence of the nitude, for 1.1GeV W 3GeV. ρ0ρ0 data, as shown in Ref. [2] and indicated by γγ The Q2-dependence of the cross section of the the line in Fig. 5(b). The fit finds a cross section of ∗ process γγ → ρ0ρ0 is well reproduced by a para- 13.6 ± 0.7nbfortheWγγ region 1.1GeV Wγγ 2 2 + − metrisation based on the GVDM model over the entire 3 GeV at Q = 0. The Q -evolution of ρ ρ data + − Q2-region. On the other hand, the ρ ρ data cannot cannot be satisfactorily described by this form. In ad- 0 0 be satisfactorily described by such a parametrisation. dition, as shown in Fig. 5(b), the ρ ρ data cannot 2 2 A ρ-pole description of the Q -dependence for both be described by the much steeper Q -fall of a ρ-pole 0 0 + − + − ρ ρ and ρ ρ data is excluded. parametrisation [2]; the same is true for the ρ ρ + − The relative magnitude of ρ ρ and ρ0ρ0 produc- cross section since it is decreasing with Q2 less steeply tion changes in the vicinity of Q2 ≈ 1GeV2, suggest- than the ρ0ρ0 one. ing different ρ-pair production mechanisms at low and high Q2. 6. Conclusions We have performed the first measurement of exclu- References sive ρ+ρ− production in mid-virtuality two-photon collisions. These results complement previous L3 [1] L3 Collaboration, P. Achard, et al., Phys. Lett. B 568 (2003) + − measurements of exclusive ρ ρ and ρ0ρ0 produc- 11. 30 L3 Collaboration / Physics Letters B 615 (2005) 19–30 [2] L3 Collaboration, P. Achard, et al., Phys. Lett. B 604 (2004) G. Alexander, et al., Z. Phys. C 30 (1986) 65; 48. B.A. Li, K.F. Liu, Phys. Lett. B 118 (1982) 435; [3] L3 Collaboration, P. Achard, et al., Phys. Lett. B 597 (2004) B.A. Li, K.F. Liu, Phys. Lett. B 124 (1983) 550; 26. B.A. Li, K.F. Liu, Phys. Rev. D 30 (1984) 613; [4] M. Diehl, et al., Phys. Rev. Lett. 81 (1998) 1782; B.A. Li, K.F. Liu, Phys. Rev. Lett. 58 (1987) 2288; N. Kivel, L. Mankiewicz, M.V. Polyakov, Phys. Lett. B 467 S.J. Brodsky, G. Köpp, P.M. Zerwas, Phys. Rev. Lett. 58 (1987) (1999) 263; 443. A. Freund, Phys. Rev. D 61 (2000) 074010. [9] J.L. Rosner, Phys. Rev. D 70 (2004) 034028. [5] I.V. Anikin, B. Pire, O.V. Teryaev, Phys. Rev. D 69 (2004) [10] L3 Collaboration, B. Adeva, et al., Nucl. Instrum. Methods 014018. A 289 (1990) 35; [6] ARGUS Collaboration, H. Albrecht, et al., Phys. Lett. B 217 M. Chemarin, et al., Nucl. Instrum. Methods A 349 (1994) 345; (1989) 205; M. Acciarri, et al., Nucl. Instrum. Methods A 351 (1994) 300; ARGUS Collaboration, H. Albrecht, et al., Phys. Lett. B 267 I.C. Brock, et al., Nucl. Instrum. Methods A 381 (1996) 236; (1991) 535; A. Adam, et al., Nucl. Instrum. Methods A 383 (1996) 342. CELLO Collaboration, H.-J. Behrend, et al., Phys. Lett. B 218 [11] T. van Rhee, Ph.D. thesis, University of Utrecht, 1999. (1989) 493. [12] P. Béné, et al., Nucl. Instrum. Methods A 306 (1991) 150; [7] TASSO Collaboration, R. Brandelik, et al., Phys. Lett. B 97 D. Haas, et al., Nucl. Instrum. Methods A 420 (1999) 101. (1980) 448; [13] F.L. Linde, Ph.D. thesis, Rijksuniversiteit Leiden, 1988. M. Althoff, et al., Z. Phys. C 16 (1982) 13; [14] V.M. Budnev, et al., Phys. Rep. C 15 (1975) 181. MARK II Collaboration, D.L. Burke, et al., Phys. Lett. B 103 [15] GEANT version 3.21 is used; (1981) 153; R. Brun, et al., CERN report CERN DD/EE/84-1 (1984), re- CELLO Collaboration, H.-J. Behrend, et al., Z. Phys. C 21 vised 1987. (1984) 205; [16] GHEISHA Collaboration, H. Fesefeldt, et al., RWTH Aachen TPC/Two-Gamma Collaboration, H. Aihara, et al., Phys. Rev. report PITHA 85/2 (1985). D 37 (1988) 28; [17] J.J. Sakurai, D. Schildknecht, Phys. Lett. B 40 (1972) 121; PLUTO Collaboration, Ch. Berger, et al., Z. Phys. C 38 (1988) I.F. Ginzburg, V.G. Serbo, Phys. Lett. B 109 (1982) 231. 521; [18] D.M. Schmidt, R.J. Morrison, M.S. Witherell, Nucl. Instrum. ARGUS Collaboration, H. Albrecht, et al., Z. Phys. C 50 Methods A 328 (1993) 547. (1991) 1. [19] G.D. Lafferty, T.R. Wyatt, Nucl. Instrum. Methods A 355 [8] N.N. Achasov, et al., Phys. Lett. B 108 (1982) 134; (1995) 541. N.N. Achasov, et al., Z. Phys. C 16 (1982) 55; [20] G.A. Schuler, Comput. Phys. Commun. 108 (1998) 279. N.N. Achasov, et al., Phys. Lett. B 203 (1988) 309; [21] M. Diehl, private communication. G. Alexander, et al., Phys. Rev. D 26 (1982) 1198; [22] M. Diehl, T. Gousset, B. Pire, Phys. Rev. D 62 (2000) 073014. Physics Letters B 615 (2005) 31–38 www.elsevier.com/locate/physletb + − + − Measurement of the KL → e e e e decay rate A. Lai, D. Marras Dipartimento di Fisica dell’Università e Sezione dell’INFN di Cagliari, I-09100 Cagliari, Italy A. Bevan, R.S. Dosanjh, T.J. Gershon, B. Hay 1, G.E. Kalmus, C. Lazzeroni, D.J. Munday, M.D. Needham 2, E. Olaiya, M.A. Parker, T.O. White, S.A. Wotton Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK 3 G. Barr, G. Bocquet, A. Ceccucci, T. Cuhadar, D. Cundy, G. D’Agostini, N. Doble, V. Falaleev, L. Gatignon, A. Gonidec, B. Gorini, G. Govi, P. Grafström, W. Kubischta, A. Lacourt, A. Norton, S. Palestini, B. Panzer-Steindel, G. Tatishvili 4,H.Taureg, M. Velasco, H. Wahl 5 CERN, CH-1211 Genève 23, Switzerland C. Cheshkov, P. Hristov, V. Kekelidze, D. Madigojine, N. Molokanova, Yu. Potrebenikov, A. Zinchenko Joint Institute for Nuclear Research, Dubna, Russian Federation I. Knowles, V. Martin, R. Sacco, A. Walker Department of Physics and Astronomy, University of Edinburgh, JCMB King’s Buildings, Mayfield Road, Edinburgh EH9 3JZ, UK M. Contalbrigo, P. Dalpiaz, J. Duclos, P.L. Frabetti, A. Gianoli, M. Martini, F. Petrucci, M. Savrié Dipartimento di Fisica dell’Università e Sezione dell’INFN di Ferrara, I-44100 Ferrara, Italy A. Bizzeti 6, M. Calvetti, G. Collazuol 7, G. Graziani, E. Iacopini, M. Lenti, G. Ruggiero Dipartimento di Fisica dell’Università e Sezione dell’INFN di Firenze, I-50125 Firenze, Italy 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.03.078 32 A. Lai et al. / Physics Letters B 615 (2005) 31–38 H.G. Becker, M. Eppard, H. Fox 8, K. Eppard, A. Kalter, K. Kleinknecht, U. Koch, L. Köpke, P. Lopes da Silva, P. Marouelli, I. Melzer-Pellmann, A. Peters, B. Renk, S.A. Schmidt, V. Schönharting, Y. Schué, R. Wanke, A. Winhart, M. Wittgen Institut für Physik, Universität Mainz, D-55099 Mainz, Germany 9 J.C. Chollet, L. Fayard, L. Iconomidou-Fayard, J. Ocariz, G. Unal, I. Wingerter-Seez Laboratoire de l’Accélératur Linéaire, IN2P3-CNRS, Université de Paris-Sud, 91406 Orsay, France 10 G. Anzivino, P. Cenci, E. Imbergamo, G. Lamanna 11, P. Lubrano, A. Mestvirishvili, A. Nappi, M. Pepe, M. Piccini Dipartimento di Fisica dell’Università e Sezione dell’INFN di Perugia, I-06100 Perugia, Italy R. Casali, C. Cerri, M. Cirilli 12, F. Costantini, R. Fantechi, L. Fiorini, S. Giudici, I. Mannelli, G. Pierazzini, M. Sozzi Dipartimento di Fisica, Scuola Normale Superiore e Sezione INFN di Pisa, I-56100 Pisa, Italy J.B. Cheze, J. Cogan 13, M. De Beer, P. Debu, A. Formica, R. Granier de Cassagnac 14, E. Mazzucato, B. Peyaud, R. Turlay, B. Vallage DSM/DAPNIA, CEA Saclay, F-91191 Gif-sur-Yvette, France M. Holder, A. Maier, M. Ziolkowski Fachbereich Physik, Universität Siegen, D-57068 Siegen, Germany 15 R. Arcidiacono, C. Biino, N. Cartiglia, F. Marchetto, E. Menichetti, N. Pastrone Dipartimento di Fisica Sperimentale dell’Università e Sezione dell’INFN di Torino, I-10125 Torino, Italy J. Nassalski, E. Rondio, M. Szleper, W. Wislicki, S. Wronka Soltan Institute for Nuclear Studies, Laboratory for High Energy Physics, PL-00-681 Warsaw, Poland 16 H. Dibon, G. Fischer, M. Jeitler, M. Markytan, I. Mikulec, G. Neuhofer, M. Pernicka, A. Taurok, L. Widhalm Österreichische Akademie der Wissenschaften, Institut für Hochenergiephysik, A-1050 Wien, Austria 17 Received 16 March 2005; accepted 31 March 2005 Available online 9 April 2005 Editor: W.-D. Schlatter A. Lai et al. / Physics Letters B 615 (2005) 31–38 33 Abstract + − + − The decay rate of the long-lived neutral K mesonintothee e e e final state has been measured with the NA48 detector at the CERN SPS. Using data collected in 1998 and 1999, a total of 200 events has been observed with negligible background. + − + − −8 This observation corresponds to a branching ratio of Br(KL → e e e e ) = (3.30 ± 0.24stat ± 0.23syst ± 0.10norm) × 10 . 2005 Elsevier B.V. All rights reserved. 1. Introduction [1,2] and thus depends on the structure of the KL → γ ∗γ ∗ vertex. Phenomenological models include vec- tor meson dominance of the photon propagator [3], + − + − The KL → e e e e decay is expected to pro- QCD inspired models [4], intermediate pseudo-scalar ∗ ∗ ceed mainly via the intermediate state KL → γ γ and vector mesons [5] and models based on chi- ral perturbation theory [6]. The probability for both virtual photons to convert into e+e− pairs is calcu- lated to be in the range (5.89–6.50) × 10−5 [2,7]. E-mail address: [email protected] (M. Lenti). 1 The chiral model prediction of [7] corresponds to Present address: EP Division, CERN, 1211 Genève 23, → + − + − = × −8 Switzerland. Br(KL e e e e ) 3.85 10 , including the 2 Present address: NIKHEF, PO Box 41882, 1009 DB Amster- effect of a form factor, which increases the width by dam, The Netherlands. 4%. The interference term due to the identity of parti- 3 Funded by the UK Particle Physics and Astronomy Research cles has been calculated to change the branching ratio Council. by only 0.5%. 4 On leave from Joint Institute for Nuclear Research, Dubna The decay was first observed by the CERN NA31 141980, Russian Federation. 5 Also at Dipartimento di Fisica dell’Università e Sezione experiment [8] based on 2 observed events and has dell’INFN di Ferrara, I-44100 Ferrara, Italy. been confirmed by later measurements [9].Herewe 6 Dipartimento di Fisica dell’Università di Modena e Reggio report the result obtained from the 1998 and 1999 data Emilia, via G. Campi 213/A, I-41100 Modena, Italy. taking periods by the NA48 experiment at the CERN 7 Present address: Scuola Normale Superiore e Sezione INFN di SPS. Pisa, I-56100 Pisa, Italy. 8 Present address: Physikalisches Institut, Universität Freiburg, D-79104 Freiburg, Germany. 9 Funded by the German Federal Minister for Research and 2. Experimental setup and data taking Technology (BMBF) under contract 7MZ18P(4)-TP2. 10 Funded by Institut National de Physique des Particules et de The NA48 experiment is designed specifically to Physique Nucléaire (IN2P3), France. measure the direct CP violation parameter Re(/) 11 Present address: Dipartimento di Fisica dell’Università di Pisa e Sezione INFN di Pisa. using simultaneous beams of KL and KS. To produce 12 Present address: CERN, CH-1211 Genève 23, Switzerland. the KL beam, 450 GeV/c protons are extracted from 13 Present address: Centre de Physique des Particules de Mar- the accelerator during 2.4 s every 14.4 s and 1.1×1012 seille, Université de la Méditerranée, IN2P3-CNRS, F-13288 Mar- of these are delivered to a beryllium target. Using di- seille, France. 14 pole magnets to sweep away charged particles and Present address: Laboratoire Leprince-Ringuet, Ecole Poly- collimators to define a narrow beam, a neutral beam technique, IN2P3-CNRS, F-91128 Palaiseau, France. × 7 ± 15 Funded by the German Federal Minister for Research and of 2 10 KL per burst and divergence 0.15 mrad Technology (BMBF) under contract 056SI74. enters the decay region. The fiducial volume begins 16 Supported by the Committee for Scientific Research grant 126 m downstream of the target and is contained in 2P03B07615 and using computing resources of the Interdisciplinary an evacuated cylindrical steel vessel 89 m long and Center for Mathematical and Computational Modelling of the Uni- 2.4 m in maximum diameter. The vessel is terminated versity of Warsaw. 17 Funded by the Austrian Ministery for Traffic and Research un- at the downstream end by a Kevlar-fiber composite −3 der the contract GZ 616.360/2-IV GZ 616.363/2-VIII, and by the window of a thickness corresponding to 3 × 10 Fonds für Wissenschaft und Forschung FWF No. P08929-PHY. radiation length and is followed immediately by the 34 A. Lai et al. / Physics Letters B 615 (2005) 31–38 main NA48 detector. The sub-detectors which are used More details on the 4-track trigger can be found in + − + − in the KL → e e e e analysis are described be- [12]. low. A magnetic spectrometer consisting of a dipole magnet is preceded and followed by two sets of drift 3. Data analysis chambers. The drift chambers are each comprised of eight planes of sense wires, two horizontal, two The decay vertex of candidate events was recon- vertical and two along each of the 45◦ directions. structed from the 4-track barycenter position in the Only the vertical and horizontal planes are instru- transverse direction, calculated as a function of the mented in the third chamber. The volume between the vertex longitudinal position; each track is weighted by chambers is filled with helium at atmospheric pres- its momentum to take into account the multiple scat- sure. The momentum resolution is p/p = 0.65% at tering effect. The 4-track vertex longitudinal position 45 GeV/c. was calculated by minimizing the sum of the squared Two segmented plastic scintillator hodoscope planes transverse distance of each track from the transverse are placed after the helium tank and provide signals for vertex position; the closest distance of approach of the the trigger. 4-track vertex is defined as the square root of this sum A liquid krypton filled calorimeter (LKr) is used at its minimum. for measuring the energy, position and time of elec- Events were preselected by requiring two positive tromagnetic showers. Space and time resolutions of and two negative tracks; each couple of tracks must better than 1.3 mm and 300 ps, respectively, have been form a 2-track vertex with distance of closest approach achieved for energies above 20 GeV. The energy res- smaller than 10 cm and an axial position of the ver- olution was determined to be σ(E) = 0√.032 ⊕ 0.090 ⊕ tex less than 210 m downstream of the target; each E E E 0.0042, with E measured in GeV. track must be compatible in time with any other within A hadron calorimeter composed of 48 steel plates, 8ns. each 24 mm thick, interleaved with scintillator is The tracks extrapolated at the LKr were required to used in trigger formation and for particle identifica- be in a fiducial area given by an octagon about 5 cm tion. smaller than the outside perimeter of the calorimeter A detailed description of the detector can be found and an inner radius of 15 cm; the distance to any dead in [10]. calorimeter cell (about 80 out of 13500) had to exceed The data used in this analysis were recorded in the 2 cm. The separation between each track extrapolated 1998–1999 data taking period. Candidate events were at LKr entry face was required to be greater than 5 cm. selected by a two-stage trigger. At the first level, a trig- The momentum of each track, measured by the mag- ger requiring adjacent hits in the hodoscope is put in netic spectrometer, was required to exceed 2 GeV, well coincidence with a total energy condition ( 35 GeV), above the detector noise of 100 MeV per cluster in the defined by adding the energy deposited in the hadronic LKr. calorimeter with that seen by the trigger in the LKr Electron candidates were identified by requiring calorimeter. The second level trigger uses information that cluster centers in the LKr be within 1.5 cm of the from the drift chambers to reconstruct tracks and in- extrapolation of each track (the rms width of electro- variant masses. For the 4-track part of the trigger, the magnetic showers in LKr is 2.2 cm). To reject pion number of clustered hits in each of the first, second, showers, the ratio of cluster energy to track momen- and fourth drift chamber had to be between 3 and 7. tum E/p was required to be greater than 0.9. Track- All possible 2-track vertices were calculated online. associated clusters with E/p < 0.8 were classified as At least two vertices within 6 m of each other in the pions. axial direction had to be found. In order to determine The fiducial volume was defined by the axial po- the efficiency of the 4-track trigger, downscaled events sition of the vertex being between 127.5 and 210 m that passed a trigger based just on the total energy con- downstream of the target. Within this volume, 4-track dition were recorded (the downscaling factor changed vertices were determined with a typical longitudinal from 100 to 60 depending on the data taking period). resolution of 0.5 m, as estimated by the Monte Carlo A. Lai et al. / Physics Letters B 615 (2005) 31–38 35 simulation. The total energy had to be within 50 and 200 GeV. + − + − 3.1. Selection of KL → e e e e candidates and background rejection + − + − Candidate events for the decay KL → e e e e with all tracks identified as electrons were selected. The following four classes of background sources were relevant: • Events with two decays KL → πeν occurring at the same time and for which the pions were misidentified as electrons. Being due to two co- incident kaon decays the invariant mass of the system could be around and above the nominal KL mass. These events were largely rejected by requiring a good vertex quality: the 4-track ver- + − + − tex closest distance of approach (defined above) Fig. 1. Correlation of e e e e invariant mass with the squared 2 had to be smaller than 5 cm. Because of miss- transverse momentum pt of the reconstructed kaon. The box is the ing transverse momentum in this and most other signal region. background decays, we required the square of 2 + − the transverse momentum pt of the reconstructed ton conversion to a e e pair. Each pair of op- kaon with respect to the line joining the de- positely charged tracks was therefore required to cay vertex and the KL target to be less than be separated by more than 2 cm in the first drift 0.0005 (GeV/c)2. We chose not to cut harder chamber. According to the Monte Carlo simu- in order to include most signal events with final lation, there was no remaining background with state radiation. The position of the cut is indicated converted photons. + − + − in Fig. 1. The Monte Carlo simulation indicates • Events KL → π π e e [11,12], with the pions that 8.3% of the signal events were lost by the misidentified as electrons. Due to the misidenti- 2 requirement on pt . In addition, as already men- fication probability of 0.5% [13] this background tioned in the preselection, it was required that each was found to be negligible. track be compatible in time with any other within + − + − 8 ns. A study of sidebands in this time distribution The invariant e e e e mass distribution result- shows that the background from this source was ing from this selection is shown in Fig. 2. Note the negligible. slightly asymmetric shape of the KL mass peak, which 0 0 0 0 0 0 • Events KL → π π ,π π π , where the π s un- is due to photons radiated off the electrons in the final dergo single or double Dalitz decays or photons state. convert in the material of the detector, so that 2 Finally, a mass window of 475 MeV/c2 < + − + − positive and 2 negative electrons are detected. Due m(e e e e )<515 MeV/c2 was set to define the to the missing photons, the invariant mass of the final sample. In total, 200 candidate events were se- + − + − e e e e system is below the nominal KL mass. lected, 62 from the 1998 data period and 138 from the + − • Events KL → γγ and KL → e e γ , with con- 1999 one. version of the photons in the material upstream The Monte Carlo simulation of the background → + − 0 of the spectrometer also yield invariant masses shows that the contribution from KL π π πDalitz, → 0 0 0 → 0 0 around the nominal KL mass. The conversion KL π πDalitzπDalitz and KL πDalitzπDalitz decays probability in the material of the NA48 detector in the signal region was negligible (less than 0.1% at is of similar magnitude as that for internal pho- 90% C.L.). As a cross-check in the data, we defined 36 A. Lai et al. / Physics Letters B 615 (2005) 31–38 In addition, at least one extra cluster in the calorime- ter with energy larger than 2 GeV, separated by more than 15 cm from each extrapolated track was required, with a time compatible within 3 ns with the average track time. The invariant mass of the e+e−γ system was re- quired to be in the range of 120–140 MeV/c2 and the + − 0 invariant mass of the π π πDalitz system had to be in the range of 475–515 MeV/c2. Monte Carlo studies + − 0 showed that background from KL → π π π with one of the external photons converting in the mater- ial of the detector was completely eliminated by the requirement on the minimum distance of the electron tracks in the first drift chamber being larger than 2 cm. All other backgrounds have been estimated to be neg- ligible. After applying all selection criteria, a total of × 6 → + − 0 2.988 10 KL π π πDalitz decays were found (0.822 × 106 in the 1998 data sample and 2.166 × 106 + − + − Fig. 2. Invariant mass of the e e e e system. The data are shown in the 1999 one). as dots with error bars, the Monte Carlo prediction for the signal and backgrounds, normalized to the data, is shown as histogram. The position of the mass cut is also indicated. 3.3. Acceptance determination and systematic uncertainty 2 2 a control region 0.0010 The difference between signal and normalization Table 1 → + − + − trigger efficiency should cancel when computing the Systematic uncertainty contributions to Br(KL e e e e ) decay rate. We cross-checked this assumption by com- Source paring the 4-track trigger efficiency as a function of Trigger efficiency ±2.0% 4-track total energy (instead of the kaon energy) with Background estimation ±1.0% the signal spectrum. We also used a partially biased not E/p cut efficiency ±2.0% ± downscaled trigger based only on LKr information to Detector acceptance 3.0% Radiative corrections ±5.6% cross-check directly the signal trigger efficiency, ob- ± taining results in agreement with the previous method. Total 7.0% We assigned a systematic uncertainty of 2% to the 4-track trigger efficiency. the momentum dependence of the E/p cut efficiency In order to be insensitive to the real kaon energy were taken into account. The correction induced by the spectrum, the acceptance correction was applied to E/p cut decreases the branching ratio by a factor of the data in bins of kaon energy (5 GeV wide) and 0.944 ± 0.020 where the error is systematic. + − + − the corresponding average KL → e e e e /KL → In Table 1 the different contributions to the system- + − π π π 0 acceptance ratio was evaluated to be atic uncertainty are listed. Dalitz + − + − 2.89. The average acceptance for KL → e e e e → + − 0 was 5.67% and 2.03% for KL π π πDalitz decays, for events generated in the range 50 GeV 0.03syst and γCP = 0.13 ± 0.10stat ± 0.03syst. Clearly, the obtained precision on the parameters βCP and γCP is limited by the event statistics. The main sources for the systematic uncertainties are the φ dependence of the detector acceptance and the effect of wrong pair- ings. By imposing γCP = 0, the fitted value of βCP is −0.13 ± 0.10stat. Our result is consistent with the hy- pothesis of a CP =−1 amplitude as expected from a KL decay. Acknowledgements It is a pleasure to thank the technical staff of the par- ticipating laboratories, universities and affiliated com- puting centers for their efforts in the construction of the NA48 apparatus, in operation of the experiment, Fig. 3. The acceptance corrected φ distribution of the data (dots with and in the processing of the data. error bars) in arbitrary units. The result of the fit with βCP and γCP as free parameters is shown as a solid line; the fit with γCP ≡ 0is also shown (dashed line). References no direct CP violation, the angular distribution in φ is [1] Z.E.S. Uy, Phys. Rev. D 43 (1991) 802; given by [17,18] Z.E.S. Uy, Phys. Rev. D 43 (1991) 1572; dn T. Miyasaki, Nuovo Cimento 5 (1972) 125. ∝ 1 + βCP cos 2φ + γCP sin 2φ, [2] T. Miyazaki, E. Taksugi, Phys. Rev. D 8 (1973) 2051. dφ [3] J.J. Sakurai, Phys. Rev. 156 (1967) 1508; with M. Moshe, P. Singer, Phys. Rev. D 6 (1972) 1379. [4] M.A. Shifman, et al., Nucl. Phys. B 120 (1977) 316. 1 −|r|2 2Re(r) [5] L. Bergström, et al., Phys. Lett. B 131 (1983) 229. β = B, γ = C, CP 1 +|r|2 CP 1 +|r|2 [6] G. D’Ambrosio, J. Portoles, J. Nucl. Phys. B 492 (1997) 417; J.L. Goity, L. Zhang, Phys. Lett. B 398 (1997) 3874. and the parameter of indirect CP violation and the [7] L. Zhang, J.L. Goity, Phys. Rev. D 57 (1998) 7031. ratio r of the amplitudes for K1 and K2 decaying into [8] G.D. Barr, et al., Phys. Lett. B 259 (1991) 389. the e+e−e+e− final state. In the limit of no indirect [9] M.R. Vagins, et al., Phys. Rev. Lett. 71 (1993) 35; CP violation, the term in sin 2φ vanishes and the clas- T. Akagi, et al., Phys. Rev. D 47 (1993) R2644; P. Gu, et al., Phys. Rev. Lett. 72 (1994) 3000; sical formula by Kroll and Wada is obtained, with the G.D. Barr, et al., Z. Phys. C 65 (1995) 361; constant B =±0.20 for CP =±1 coming from inte- T. Akagi, et al., Phys. Rev. D 51 (1995) 2061; gration over phase space. At the moment no calcula- A. Alavi-Harati, et al., Phys. Rev. Lett. 86 (2001) 5425. tion exists for C. [10] A. Lai, et al., Eur. Phys. J. C 22 (2001) 231. For the measurement of the angle φ, the ambigu- [11] A. Alavi-Harati, et al., Phys. Rev. Lett. 84 (2000) 408. + − [12] A. Lai, et al., Eur. Phys. J. C 30 (2003) 33. ity in the e e pairing was resolved by choosing the [13] A. Lai, et al., Phys. Lett. B 602 (2004) 41. combination that minimized the product of invariant [14] Particle Data Group, Phys. Lett. B 592 (2004) 1. masses of the two pairs. Monte Carlo studies showed [15] E. Barbario, B. van Eijk, Z. Was, Comput. Phys. Commun. 66 that in 98% of the cases we obtain the correct combi- (1991) 115. nation by this method. The remaining wrong pairings [16] GEANT Detector Description and Simulation Tool, CERN Program Library Long Write-up W5013, 1994. are uniformly distributed in φ. The acceptance cor- [17] N.M. Kroll, W. Wada, Phys. Rev. 98 (1955) 1355. rected φ distribution is shown in Fig. 3. By fitting this [18] P. Gu, et al., Phys. Rev. Lett. 72 (1994) 3000. distribution we measured βCP =−0.13 ± 0.10stat ± Physics Letters B 615 (2005) 39–49 www.elsevier.com/locate/physletb Measurement of the γγ → π +π − and γγ → K+K− processes at energies of 2.4–4.1 GeV Belle Collaboration H. Nakazawa j, S. Uehara j,K.Abej,K.Abeat,H.Aiharaav,M.Akatsux,Y.Asanoba, V. Aulchenko b,T.Aushevn, S. Bahinipati f,A.M.Bakichaq,Y.Banaj, S. Banerjee ar, I. Bedny b,U.Bitenco, I. Bizjak o, S. Blyth ac, A. Bondar b, A. Bozek ad,M.Brackoˇ j,v,o, J. Brodzicka ad,A.Chenz,B.G.Cheond,R.Chistovn, Y. Choi ap, Y.K. Choi ap, A. Chuvikov ak,J.Dalsenow, M. Danilov n, M. Dash bc,A.Drutskoyf, S. Eidelman b, Y. E n ar i x,D.Epifanovb,S.Fratinao, N. Gabyshev b,A.Garmashak, T. Gershon j, G. Gokhroo ar, K. Hayasaka x, H. Hayashii y, Y. Hoshi at,S.Houz,W.-S.Houac, T. Iijima x, A. Imoto y, K. Inami x,A.Ishikawaj,R.Itohj,M.Iwasakiav,Y.Iwasakij, J.H. Kang bd, J.S. Kang q, S.U. Kataoka y, N. Katayama j,H.Kawaic, T. Kawasaki af, H.R. Khan aw,H.Kichimij,H.J.Kims,J.H.Kimap,S.K.Kimao,S.M.Kimap, S. Korpar v,o, P. Križan u,o,P.Krokovnyb, R. Kulasiri f,C.C.Kuoz,A.Kuzminb, Y. - J . Kwo n bd,G.Lederm,S.E.Leeao,T.Lesiakad,J.Lian,S.-W.Linac,D.Liventsevn, G. Majumder ar,F.Mandlm, T. Matsumoto ax,A.Matyjaad,W.Mitaroffm,H.Miyakeah, H. Miyata af,R.Mizukn, D. Mohapatra bc,T.Moriaw, T. Nagamine au, Y. Nagasaka k, E. Nakano ag, M. Nakao j, Z. Natkaniec ad,S.Nishidaj,O.Nitohay, T. Nozaki j, S. Ogawa as, T. Ohshima x, T. Okabe x, S. Okuno p,S.L.Olseni, W. Ostrowicz ad, P. Pakhlov n,H.Palkaad,C.W.Parkap,H.Parks,N.Parslowaq, L.S. Peak aq, R. Pestotnik o, L.E. Piilonen bc,H.Sagawaj,Y.Sakaij,N.Satox, T. Schietinger t, O. Schneider t, J. Schümann ac,K.Senyox, H. Shibuya as,B.Shwartzb, J.B. Singh ai, A. Somov f,N.Soniai, R. Stamen j,S.Stanicˇ ba,1,M.Staricˇ o,K.Sumisawaah, T. Sumiyoshi ax, S.Y. Suzuki j, O. Tajima j, F. Takasaki j,N.Tamuraaf, M. Tanaka j, Y. Teramoto ag, X.C. Tian aj, T. Tsuboyama j, T. Tsukamoto j,T.Uglovn,S.Unoj, G. Varner i, S. Villa t,C.C.Wangac,C.H.Wangab, M. Watanabe af, Y. Watanabe aw, B.D. Yabsley bc, A. Yamaguchi au,Y.Yamashitaae, M. Yamauchi j, J. Ying aj,Y.Yusaau, L.M. Zhang an,Z.P.Zhangan, V. Zhilich b,D.Žontaru,o 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.03.067 40 Belle Collaboration / Physics Letters B 615 (2005) 39–49 a Aomori University, Aomori, Japan b Budker Institute of Nuclear Physics, Novosibirsk, Russia c Chiba University, Chiba, Japan d Chonnam National University, Kwangju, South Korea e Chuo University, Tokyo, Japan f University of Cincinnati, Cincinnati, OH, USA g University of Frankfurt, Frankfurt, Germany h Gyeongsang National University, Chinju, South Korea i University of Hawaii, Honolulu, HI, USA j High Energy Accelerator Research Organization (KEK), Tsukuba, Japan k Hiroshima Institute of Technology, Hiroshima, Japan l Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, PR China m Institute of High Energy Physics, Vienna, Austria n Institute for Theoretical and Experimental Physics, Moscow, Russia o J. Stefan Institute, Ljubljana, Slovenia p Kanagawa University, Yokohama, Japan q Korea University, Seoul, South Korea r Kyoto University, Kyoto, Japan s Kyungpook National University, Taegu, South Korea t Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne, Switzerland u University of Ljubljana, Ljubljana, Slovenia v University of Maribor, Maribor, Slovenia w University of Melbourne, Victoria, Australia x Nagoya University, Nagoya, Japan y Nara Women’s University, Nara, Japan z National Central University, Chung-li, Taiwan aa National Kaohsiung Normal University, Kaohsiung, Taiwan ab National United University, Miao Li, Taiwan ac Department of Physics, National Taiwan University, Taipei, Taiwan ad H. Niewodniczanski Institute of Nuclear Physics, Krakow, Poland ae Nihon Dental College, Niigata, Japan af Niigata University, Niigata, Japan ag Osaka City University, Osaka, Japan ah Osaka University, Osaka, Japan ai Panjab University, Chandigarh, India aj Peking University, Beijing, PR China ak Princeton University, Princeton, NJ, USA al RIKEN BNL Research Center, Brookhaven, NY, USA am Saga University, Saga, Japan an University of Science and Technology of China, Hefei, PR China ao Seoul National University, Seoul, South Korea ap Sungkyunkwan University, Suwon, South Korea aq University of Sydney, Sydney, NSW, Australia ar Tata Institute of Fundamental Research, Bombay, India as Toho University, Funabashi, Japan at Tohoku Gakuin University, Tagajo, Japan au Tohoku University, Sendai, Japan av Department of Physics, University of Tokyo, Tokyo, Japan aw Tokyo Institute of Technology, Tokyo, Japan ax Tokyo Metropolitan University, Tokyo, Japan ay Tokyo University of Agriculture and Technology, Tokyo, Japan az Toyama National College of Maritime Technology, Toyama, Japan ba University of Tsukuba, Tsukuba, Japan bb Utkal University, Bhubaneswer, India bc Virginia Polytechnic Institute and State University, Blacksburg, VA, USA bd Yonsei University, Seoul, South Korea Belle Collaboration / Physics Letters B 615 (2005) 39–49 41 Received 22 December 2004; accepted 29 March 2005 Available online 7 April 2005 Editor: M. Doser Abstract + − + − − We have measured π π and K K production in two-photon collisions using 87.7fb 1 of data collected with the Belle + − detector at the asymmetric energy e e collider KEKB. The cross sections are measured to high precision in the two-photon ∗ center-of-mass energy (W) range between 2.4GeV PACS: 12.38.Qk; 13.25.Gv; 13.66.Bc; 13.85.Lg Keywords: Two-photon collisions; Mesons; QCD; Charmonium 1. Introduction pion or kaon pairs at large energy and momentum transfers, in which the amplitude is expressed by a Exclusive processes with hadronic final states test hard γγ → qq¯ subprocess and a form factor describ- various model calculations motivated by perturbative ing the soft transition from qq¯ to the meson pair. DKV, and non-perturbative QCD. Two-photon production of as well as BL, predict the sin−4 θ ∗ dependence of the exclusive hadronic final states is particularly attractive angular differential cross section, which is an impor- due to the absence of strong interactions in the initial tant test of these approaches. It is interesting to in- state and the possibility of calculating γγ → qq¯ am- vestigate experimentally an energy scale where those plitudes. The perturbative QCD calculation by Brod- theoretical predictions become valid. The recent mea- sky and Lepage (BL) [1] is based on factorization surements of γγ → π +π − and K+K− performed by of the amplitude into a hard scattering amplitude for ALEPH [5] are consistent, within their errors, with the γγ → qqq¯ q¯ and a single-meson distribution ampli- BL’s prediction of the energy dependence, but not the tude. Their prediction gives the dependence√ on the normalization. However, their dataset is not sufficient center-of-mass (c.m.) energy W(≡ s) and scatter- to conclusively test the W and sin−4 θ ∗ dependences. ing angle θ ∗ for γγ → M+M− processes In this report, we measure γγ → π +π − and → + − 2 2 γγ K K processes with high precision, and dσ + − 16πα |F (s)| γγ → M M ≈ M , (1) make quantitative comparisons with QCD predictions. | ∗| 4 ∗ − d cos θ s sin θ This analysis is based on an 87.7fb 1 data sample col- where M represents a meson and FM denotes its elec- lected at or near the Υ(4S) resonance energy, accumu- tromagnetic form factor. Vogt [2], based on the pertur- lated with the Belle detector [6] located at KEKB [7]. bative approach, claimed a need for soft contributions, as his result for the hard contribution was well below the experimental cross section obtained by CLEO [3]. 2. KEKB accelerator and Belle detector Diehl, Kroll and Vogt (DKV) proposed [4] the soft handbag contribution to two-photon annihilation into KEKB is a colliding beam accelerator of 8 GeV electrons and 3.5 GeV positrons designed to produce ¯ E-mail address: [email protected] (H. Nakazawa). copious BB meson pairs to observe CP violation. 1 On leave from Nova Gorica Polytechnic, Nova Gorica, Slove- The Belle detector, with a 1.5 T solenoidal mag- nia. netic field, surrounds the interaction point and sub- 42 Belle Collaboration / Physics Letters B 615 (2005) 39–49 ◦ ◦ tends the polar angle range 17 <θlab < 150 , mea- ened in the next section), the invariant mass of these sured from the z axis, which is aligned opposite the two tracks be below 4.5GeV/c2, and that the squared positron beam. It is described in detail in Ref. [6]. missing mass of the event be above 2 GeV2/c4. Here, Briefly, charged track momenta and their decay points the two tracks are assumed to be massless parti- are measured by the central drift chamber (CDC) and cles. The latter two requirements eliminate radiative silicon vertex detector (SVD). The hadron identity of Bhabha and initial state radiation events. The remain- these charged tracks is determined using information ing events consist of two-photon production of e+e−, from the time-of-flight counters (TOF), the aerogel µ+µ−, π +π −, K+K−, and pp¯ final states as well as threshold Cerenkovˇ counters (ACC), and the specific unvetoed e+e− → τ +τ − events according to a Monte ionization in the CDC. Hadron/electron discrimina- Carlo (MC) study [8]. tion is performed using the above information as well The predicted versus measured range and trans- as the energy deposition and shower profile in the verse deviation of hits in the KLM are used to con- segmented CsI(Tl)-crystal electromagnetic calorime- struct a normalized likelihood Rµ that a track extrap- ter (ECL). Hadron/muon discrimination is achieved olated from the CDC is a muon rather than a pion using information from the neutral-kaon and muon de- or kaon. Here, we classify an event as arising from + − tector (KLM), which consists of glass resistive plate γγ → µ µ if either track has Rµ > 0.66. Simi- counters embedded in the solenoid’s iron flux return. larly, the TOF, ACC, CDC, and ECL information is used to construct a normalized likelihood Re that a re- 3. Event selection constructed track is an electron rather than a hadron. We classify an event as arising from γγ → e+e− if The signal events are collected predominantly by either track has Re > 0.66. For above two separations a trigger that requires two charged tracks penetrating 93% of signal events survive for both modes accord- through the CDC and TOF, with an opening angle in ing to the MC study described later. The TOF, ACC the rϕ plane (perpendicular to the z axis) of at least and CDC information is used to construct another nor- ◦ 135 . malized likelihood Rp that a reconstructed track is We select signal candidates according to the fol- a proton rather than a kaon, with a high value corre- lowing criteria. There must be exactly two oppositely- sponding to a proton-like track. The scatterplot of this charged reconstructed tracks satisfying the following quantity for the negative track vs. that for the positive conditions: −0.47 cos θlab 0.82 for the polar an- track in each event is shown in Fig. 1(a). Note that the + − + − gle θlab of each track; pt > 0.8GeV/c for the mo- peak near the origin contains both K K and π π mentum component in the rϕ plane of each track; candidates. Events above the hyperbolic curve dr 1 cm and |dz| < 2 cm for the origin of each − + + − R − . R − . = . , track relative to the nominal e e collision point; and p 1 01 p 1 01 0 0101 |dz1 − dz2| 1 cm for the two tracks’ origin differ- shown in the inset of Fig. 1(a), are deemed to arise ence along the z axis, where the origin is defined by the from γγ → pp¯. closest approach of the track to the nominal collision After removing events that appear to arise from point in the rϕ plane. The event is vetoed if it contains two-photon production of µ+µ−, e+e−, and pp¯ ac- any other reconstructed charged track with transverse cording to the above criteria, the remaining sample momentum above 0.1GeV/c. consists of two-photon production of K+K−, π +π −, Cosmic rays are suppressed by demanding that and residual µ+µ−,aswellase+e− → τ +τ − pro- the opening angle α between the two tracks satisfy duction where each τ lepton decays to a single pion cos α −0.997. The signal is enriched relative to or muon. Information from the TOF, ACC, and CDC other backgrounds by requiring that the scalar sum of is used to form a normalized likelihood ratio RK that the momentum of the two tracks be below 6 GeV/c, a reconstructed track is a kaon rather than a pion (or the total energy deposited in the ECL be below 6 GeV, muon), with a high value corresponding to a kaon- the magnitude of the net transverse momentum of the like track. The scatter plot of this quantity for the above-selected two charged tracks in the e+e− c.m. negative track vs. that for the positive track is shown frame be below 0.2GeV/c (this condition is tight- in Fig. 1(b). Note that the peak near the origin con- Belle Collaboration / Physics Letters B 615 (2005) 39–49 43 (a) (b) Fig. 1. Two-dimensional plots of likelihood ratios for hadron identification: (a) Rp and (b) RK . The cut boundaries are shown in the top-view insets. tains π +π −, τ +τ −, and residual µ+µ− events. Events 4. Background rejection − + above the diagonal line R + R = 1.2, shown in K K + − the inset of Fig. 1(b), are classified as K K candi- The spectrum of the residual γγ → µ+µ− back- − + dates, while events below the line R + R = 0.8are ground within the π +π − sample is obtained from a + − K K + − classified as π π candidates (including τ τ and MC simulation program AAFH [9], based on a full + − residual µ µ backgrounds). Events in the diagonal O(α4) QED calculation, with a data sample corre- band between these two lines are discarded. sponding to an integrated luminosity of 174.2fb−1 + − The π π sample is somewhat contaminated by that is processed by the full detector simulation pro- + − non-exclusive two-photon background γγ → π π X gram and then subjected to trigger simulation and the + − + − as well as the e e → τ τ process, in roughly equal above event selection criteria. After calibration of the proportion. We note that these backgrounds appear at muon identification efficiency to match that in the high values of the magnitude of the net transverse mo- data using identified γγ → µ+µ− events, the resid- | + + −| + − + − mentum pt pt in the e e c.m. frame, and are ual µ µ background is scaled by the integrated lu- often accompanied by photons from the prompt decay minosity ratio and then subtracted from the π +π − of a neutral pion in the final state. Therefore, we reject sample. (The contamination amounts to 19–42%, de- + − events in the π π sample that contain a photon with pending on W .) energy above 400 MeV (Eγ -veto). The distributions The excess in the Eγ -vetoed histogram of Fig. 2(a) | + + −| + − of pt pt for the π π candidates before and af- above the smooth curve from the signal MC, de- ter application of this veto are shown as the histograms scribed in more detail below, is attributed to non- in Fig. 2(a). exclusive γγ → π +π −X events that are not rejected + − + − + − + − The yields of the π π and K K events are by the Eγ -veto; most of the e e → τ τ events expressed as functions of three variables: W derived are rejected by this veto. A similar excess appears ∗ from the invariant mass of the two mesons, | cos θ | in Fig. 2(b) for the γγ → K+K− process. Assum- | + + −| and pt pt . Eighty-five 20 MeV wide bins in W ing that this remaining background is proportional to times six bins in the cosine of the γγ c.m. scattering net transverse momentum, we determine the slope us- ∗ angle θ times twenty bins in net transverse momen- ing the difference between data and MC in the range | + + −| tum are used in the ranges 2.4GeV + − + − + − + − Fig. 2. |p + p | distribution for π π (a) and K K (b) candidates. The dashed and solid histograms in π π indicate the distribution t t + − of events before and after Eγ -veto (which is not applied to the K K candidates), respectively. The arrows indicate the upper boundaries of | + + −| + − pt pt for the signal. The residual muon background has been subtracted from the π π distribution. The curves show the signal MC distribution which is normalized to the signal candidates at the leftmost bin. ∗ + − + − Fig. 3. Number of events (| cos θ | < 0.6) obtained for the γγ → π π (solid) and γγ → K K (dashed) samples after the background subtraction. there is no dependence on the scattering angle θ ∗.Us- γγ → π +π − (γ γ → K+K−). Thus, we obtain 6919 ing the smoothed slope, the estimated non-exclusive and 6234 signal events for π +π − and K+K−, respec- background, which amounts to at most 4% (9%) for tively. The background-subtracted yields, integrated π +π − (K+K−) below 3.1 (3.3) GeV but as much as over net transverse momentum and scattering angle, 28% (26%) in the highest W bin, is subtracted from are shown as a function of W in Fig. 3. each bin. Finally, we restrict our signal region to net The signatures of the χc0(1P) and χc2(1P) reso- transverse momentum below 0.05 (0.10) GeV/c for nances are observed in both π +π − and K+K− chan- Belle Collaboration / Physics Letters B 615 (2005) 39–49 45 nels. By fitting each W distribution outside the range value of the trigger efficiency is ∼ 93% for events in 3.3–3.7 GeV to a cubic polynomial, we see an ex- the acceptance. cess of 129 (153) events in the π +π − (K+K−) chan- The efficiency-corrected measured differential + − + − nel in the χc0 range of 3.34–3.44 GeV and a corre- cross sections for γγ → π π and γγ → K K , ∗ sponding excess of 54 (33) events in the χc2 range normalized to the partial cross section σ0 for | cos θ | < of 3.54–3.58 GeV. We obtain consistent results from 0.6, are shown in Fig. 4 for each 100 MeV wide W a fit of each distribution to a cubic polynomial plus bin. The partial cross sections σ0 for both processes, a Breit–Wigner (χc0) or a Gaussian (χc2) peak. The integrated over the above scattering angle range, are + − χc0 statistical significance is 6.2σ(8.2σ)in the π π shown in Fig. 5 (along with their ratio) and itemized (K+K−) channel, where σ is the standard deviation. in Table 1 as a function of W . The χc2 statistical significance is 4.8σ (3.7σ )inthe π +π − (K+K−) channel. The significances are taken from the square root of the difference of the goodness- 6. Systematic errors of-fit values from the two fits where the peak term in the above fit function is included or excluded. As- 4 ∗ suming a flat (sin θ ) shape for the χc0 (χc2) reso- The dominant systematic errors are listed in Ta- nance [10], we subtract the above excesses bin by bin ble 2. The error in track finding efficiency is estimated from each angular distribution in the above W ranges. by comparing η → π +π −π 0 and η → γγ samples in data and MC. The uncertainty due to trigger efficiency is estimated by comparing the yields of γγ → µ+µ− 5. Derivation of the cross section in real and simulated data [9] after accounting for the background from e+e− → µ+µ− nγ events (varying The differential cross section for a two-photon with W from 0.5–4.6%), which have the same topol- process to a two-body final state arising from an ogy [13]. To estimate the systematic error in the K/π + − + − electron–positron collision is given by separation for the π π (K K ) mode, which is dominated by the contribution from the fake rate for + − + − dσ | ∗|; → K K (π π ) events, we have compared the two- ∗ W, cos θ γγ X d| cos θ | dimensional likelihood ratio distribution of the experi- | ∗|; + − → + − mental data in Fig. 1(b) to that from the corresponding = N(W, cos θ e e e e X) ∗ ∗ , signal MC, normalized appropriately. We assign the Lγγ(W) W | cos θ | (W,| cos θ |) L dt (2) so-estimated background contribution to the system- atic error from this source. The uncertainty in the rel- where N and denote the number of the signal events ative muon identification efficiency between real and and a product of detection and trigger efficiencies, simulated data is used to determine the error asso- respectively; L dt is the integrated luminosity, and ciated with the residual µ+µ− subtraction from the L is the luminosity function, defined as γγ π +π − sample. We use an error of 100% of the sub- dσ ; + − → + − tracted value for the non-exclusive background sub- dW (W e e e e X) Lγγ(W) = . traction. We allow the number of χ events to fluctu- σ(W; γγ → X) cJ ate by up to 20% of the measured excess to estimate ∗ + − The efficiencies (W,| cos θ |) for γγ → π π the error due to the χc subtraction that is applied for and γγ → K+K− are obtained from a full Monte the energy bins in the range 3.3GeV −1 | ∗| + − + − Fig. 4. Angular dependence of the cross section, σ0 dσ/d cos θ ,fortheπ π (closed circles) and K K (open circles) processes. The − ∗ curves are 1.227 sin 4 θ . The errors are statistical only. + − + − ∗ − Fig. 5. Cross section for (a) γγ → π π ,(b)γγ → K K in the c.m. angular region | cos θ | < 0.6 together with a W 6 dependence line derived from the fit of s|RM |. (c) Cross section ratio. The solid line is the result of the fit for the data above 3 GeV. The errors indicated by short ticks are statistical only. Belle Collaboration / Physics Letters B 615 (2005) 39–49 47 Table 1 + − + − ∗ Cross sections and errors for the γγ → π π and γγ → K K processes, in the angular range | cos θ | < 0.6 γγ → π+π− γγ → K+K− W, GeV σ0, nb stat. err., nb syst. err., nb σ0, nb stat. err., nb syst. err., nb 2.4–2.5 0.832 0.026 0.083 0.595 0.019 0.062 2.5–2.6 0.625 0.024 0.070 0.549 0.018 0.059 2.6–2.7 0.636 0.025 0.067 0.497 0.018 0.054 2.7–2.8 0.530 0.023 0.059 0.377 0.016 0.043 2.8–2.9 0.448 0.022 0.048 0.364 0.016 0.043 2.9–3.0 0.407 0.021 0.042 0.313 0.015 0.038 3.0–3.1 0.302 0.019 0.032 0.266 0.014 0.034 3.1–3.2 0.247 0.017 0.026 0.205 0.013 0.027 3.2–3.3 0.202 0.016 0.022 0.199 0.013 0.027 3.3–3.4 0.153 0.016 0.023 0.130 0.012 0.022 3.4–3.5 0.103 0.014 0.018 0.093 0.011 0.018 3.5–3.6 0.093 0.014 0.017 0.093 0.011 0.016 3.6–3.7 0.079 0.011 0.013 0.070 0.009 0.012 3.7–3.8 0.068 0.011 0.014 0.062 0.008 0.011 3.8–3.9 0.043 0.009 0.010 0.048 0.007 0.009 3.9–4.0 0.035 0.008 0.010 0.039 0.007 0.008 4.0–4.1 0.041 0.009 0.013 0.039 0.007 0.008 Table 2 where M denotes either a pion or a kaon. Here, they in- Contributions to the systematic errors. A range is shown when the troduce the “annihilation form factor” R2M (s), which error has a W dependence can be determined experimentally. + − + − Source π π K K Fig. 4 shows normalized angular distributions as + − + − Tracking efficiency 4% 4% a function of W ,forπ π and K K modes. The Trigger efficiency 4% 4% solid curves indicate the expectations from a sin−4 θ ∗ K/π separation 0–1% 2–4% behavior predicted by BL and DKV models µµ background subtraction 5–17% 0% Non-exclusive background subtraction 4–27% 7–20% 1 dσ 1.227 = . (4) Luminosity function 5% 5% σ d| cos θ ∗| 4 ∗ Integrated luminosity 1% 1% 0 sin θ χc subtraction (3.3GeV Table 3 + − Results for the product of the two-photon decay width and the branching fraction, Γγγ(χcJ )B(χcJ → M M ). The second column gives the observed χcJ yields in the W region of 3.34–3.44 GeV (3.54–3.58 GeV) for χc0 (χc2). The first and second errors for ΓγγB are statistical and systematic, respectively + − Number of events Γγγ(χcJ )B(χcJ → M M ) [eV] Significance + − γγ → χc0 → π π 129 ± 18 15.1 ± 2.1 ± 2.36.2σ + − γγ → χc0 → K K 153 ± 17 14.3 ± 1.6 ± 2.38.2σ + − γγ → χc2 → π π 54 ± 10 0.76±0.14±0.11 4.8σ + − γγ → χc2 → K K 33 ± 80.44±0.11±0.07 3.7σ 2 2 0.05 GeV and s|R2K |=0.68 ± 0.01 ± 0.05 GeV rules. This compensates for the partial account of the by fits to the data in this W range. Fig. 5 shows the SU(3) breaking by BL who used the same wave func- observed cross sections in comparison with the recent tions for pions and kaons so that the cross section ratio ALEPH data [5]. The lines indicate expectations with is equal to the fourth power of the ratio of the kaon the s|R2M (s)| values from the fit above. Our data are and pion decay constants. consistent with W −6 behavior predicted by BL [1] and BF [14] models. We can also directly obtain the power n of the W - 8. Two-photon decay width of χ cJ resonances n dependence (σ0 ∝ W ) from the data. We find n = −7.9 ± 0.4 ± 1.5forπ +π − and −7.3 ± 0.3 ± 1.5 + − The measured yields of χc0 and χc2 events can be for K K ,for3.0GeV 9. Conclusion ence and Technology of India; the BK21 program of the Ministry of Education of Korea and the CHEP − Using 87.7fb 1 of data collected with the Belle SRC program of the Korea Science and Engineering detector at KEKB, we have measured with high pre- Foundation; the Polish State Committee for Scientific + − cision the cross sections for the γγ → π π and Research under contract No. 2P03B 01324; the Min- + − γγ → K K processes in the two-photon c.m. en- istry of Science and Technology of the Russian Fed- ergy range 2.4GeV Transition from baryonic to mesonic freeze-out J. Cleymans a, H. Oeschler b,K.Redlichc, S. Wheaton a a UCT-CERN Research Centre and Department of Physics, University of Cape Town, Rondebosch 7701, South Africa b Institut für Kernphysik, Technische Universität Darmstadt, D-64289 Darmstadt, Germany c Institute of Theoretical Physics, University of Wrocław, Pl-45204 Wrocław, Poland Received 6 December 2004; received in revised form 24 February 2005; accepted 22 March 2005 Available online 9 April 2005 Editor: V. Metag Abstract + + The recently discovered sharp peak in the K /π ratio in relativistic heavy-ion collisions is discussed in the framework of the statistical model. In this model a rapid change is expected as the hadronic gas undergoes a transition from a baryon- = dominated to a meson-dominated gas. The transition occurs√ at a temperature T 140 MeV and baryon chemical potential µB = 410 MeV corresponding to an incident energy of sNN = 8.2 GeV. The maximum in the Λ/π ratio is well reproduced + + by the statistical model, but the change in the K /π ratio is much less pronounced than the one observed by the NA49 − − − + − + Collaboration. The calculated smooth increase of the K /π ratio and the shape of the Ξ /π and Ω /π ratios exhibiting maxima at different incident energies is consistent with the presently available experimental data. We conclude that the measured particle ratios with 20–30% deviations agree with a hadronic freeze-out scenario. These deviations seem to occur just in the transition from baryon-dominated to meson-dominated freeze-out. 2005 Elsevier B.V. All rights reserved. The NA49 Collaboration has recently performed difference in heavy-ion collisions. This transition has a series of measurements of Pb–Pb collisions at 20, been referred to in Ref. [1] as the “horn”. A strong 30, 40, 80 and 158 A GeV beam energies [1–3]. variation with energy of the Λ/π ratio has been pre- When these results are combined with measurements dicted on the basis of arguments put forward in [10]. at lower beam energies from the AGS [4–9] they re- It has also been suggested recently in Ref. [11] that veal an unusually sharp variation with beam energy in this is a signal of the special critical point of the QCD the Λ/π, with π≡3/2(π + + π −), and K+/π + phase diagram [12–14] at high baryon density. In this ratios. Such a strong variation with energy does not Letter we explore another, less spectacular, possibility occur in pp collisions and therefore indicates a major for the origin of the sharp maximum, namely as be- ing due to the transition from a baryon-dominated to a meson-dominated hadronic gas. The distinction be- E-mail address: [email protected] (J. Cleymans). ing based on whether the entropy of the hadronic gas 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.03.074 J. Cleymans et al. / Physics Letters B 615 (2005) 50–54 51 √ Fig. 1. The net baryon density as a function of sNN calcu- Fig. 2. The entropy density normalised to T 3 as a function of the lated along the chemical freeze-out curve [15] corresponding to beam energy as calculated in the statistical model using THERMUS E/N = 1GeV. [21]. is dominated by baryons or by mesons. For this pur- pose we study various quantities√ along the freeze-out Above this value the entropy is carried mainly by curve [15] as a function of sNN. mesonic degrees of freedom. It is remarkable that the In the statistical model a steep rise at low ener- entropy density divided by T 3 is constant over the en- gies and a subsequent flattening off leading to a mild tire freeze-out curve, except for the low-energy region maximum in the K+/π + ratio, was predicted many corresponding to the SIS energy region. The line de- years ago [16–18]. The sharpness of the observed peak noting the transition from a baryon-dominated to a therefore comes as a surprise. On the other hand, a meson-dominated hadron gas is shown in Fig. 3.This sharp peak in the Λ/π ratio was predicted by the line crosses the freeze-out curve at a temperature of statistical model [18] and is in good agreement with T = 140 MeV, when the baryon chemical potential the data. While the statistical model cannot explain equals µB = 410 MeV. The corresponding invariant + + √ the sharpness of the peak in the K /π ratio, there energy is sNN = 8.2GeV. are nevertheless several phenomena, giving rise to The strong decrease in the net baryon density seen the rapid change, which warrant a closer look at the in Fig. 1 is due to the fact that low energies are char- model. In Fig. 1 we show the net baryon density cal- acterized by a very low multiplicity of mesons and, culated along√ the chemical freeze-out curve [15] as a correspondingly, a very large baryon-to-meson ratio. function of sNN. This curve shows a clear maximum As a consequence, the baryon chemical potential is with the net baryon density decreasing rapidly towards also very large. As the beam energy is increased, me- higher energies. son production increases and the baryon chemical po- To get a better estimate of the statistical parameters tential decreases. The number of strange baryons pro- in the transition region we show in Fig. 2 the entropy duced in heavy-ion collisions at different collision en- density as a function of beam energy following the ergies will follow the net baryon density since a large freeze-out curve given in [15]. The separate contri- baryon chemical potential will also enhance the num- bution of mesons and of baryons to the total entropy ber of hyperons. + + is also shown in this figure. There√ is a clear change The corresponding K /π ratio is shown in Fig. 4. of baryon to meson dominance at sNN = 8.2GeV. As is well known [19,20], the statistical model de- 52 J. Cleymans et al. / Physics Letters B 615 (2005) 50–54 Fig. 3. The chemical freeze-out curve [15] together with the regions where baryonic (or mesonic) contributions to the entropy density dominate, separated by the dashed line. Calculated in the statistical model using THERMUS [21]. scription leads to a mild maximum in this ratio which does not reproduce the so-called “horn” observed by the NA49 Collaboration [1]. The observed deviations at the highest SPS energy have been interpreted as a lack of full chemical equilibrium in the strangeness Fig. 4. (a) The Λ/π ratio as a function of beam energy. (b) The + + − − sector, leading to a strangeness suppression factor, γs , K /π and K /π ratios as a function of energy. The solid and deviating from its equilibrium value by about thirty dashed lines are the predictions of the statistical model calculated percent. Detailed fits using the statistical model in using THERMUS [21]. The data points are from Refs. [1–4,7,9]. the region of the “horn” show rapid variations in γs [20] which do not lend themselves to any interpreta- tion. There is no corresponding peak in the K−/π − combination of the facts that strangeness has to be bal- ratio because the production of K− is not tied to anced, the baryon chemical potential decreases rapidly that of baryons. As the relative number of baryons with energy and the multi-strange baryons have suc- decreases with increasing energy, there is no corre- cessively higher thresholds. The values are listed in sponding decrease in the number of K− as is the Table 1. case with K+ as these must be balanced by strange It is to be expected that if these maxima do not all baryons. occur at the same temperature, i.e., at the same beam It is worth noting that the maxima in the ratios energy, then the case for a phase transition is not very for multi-strange baryons occur at ever higher beam strong. The observed behavior seems to be governed energies. This can be seen clearly in Fig. 5 for the by properties of the hadron gas. More detailed experi- Ξ −/π + ratio which peaks at a higher value of the mental studies of multi-strange hadrons will allow the beam energy. The ratio Ω−/π + also shows a (very verification or disproval of the trends shown in this weak) maximum, as can be seen in Fig. 5. The higher Letter. It should be clear that the Ω−/π + ratio is very the strangeness content of the baryon, the higher in broad and shallow and it will be difficult to find a max- energy is the maximum. This behavior is due to a imum experimentally. J. Cleymans et al. / Physics Letters B 615 (2005) 50–54 53 If the change in properties of the above excitation functions were associated with a genuine deconfine- ment phase transition one would expect these changes to occur at the same beam energy. It is clear that more data are needed to clarify the precise nature of the sharp variation observed by the NA49 Collabora- tion. Acknowledgements We thank C. Blume for his help with the NA49 data. We acknowledge the support of the German Bun- desministerium für Bildung und Forschung (BMBF), the Polish State Committee for Scientific Research (KBN) grant 2P03 (06925), the National Research Foundation (NRF, Pretoria) and the URC of the Uni- versity of Cape Town. − + − + Fig. 5. The Ξ /π (full line) and Ω /π (dashed line) ratios as a function of beam energy calculated using THERMUS [21].The data points are from Refs. [22–24]. The square points correspond to + − + References the (Ω + Ω)/π¯ ratio while the round points are for the Ξ /π ratio. [1] M. Ga´zdzicki, NA49 Collaboration, J. Phys. G: Nucl. Part. Phys. 30 (2004) S701. [2] S.V. Afanasiev, et al., NA49 Collaboration, Phys. Rev. C 66 Table 1 (2002) 054902. Maxima in particle ratios [3] T. Anticic, et al., Phys. Rev. Lett. 93 (2004) 022302; √ Ratio Maximum at sNN (GeV) Maximum value C. Blume, et al., SQM04 Proceedings, J. Phys. G, in press. [4] L. Ahle, et al., E802 Collaboration, Phys. Rev. C 57 (1998) Λ/π 5.10.052 − + 466. Ξ /π 10.20.011 + + [5] L. Ahle, et al., E802 Collaboration, Phys. Rev. C 60 (1999) K /π 10.80.22 − + 044904; Ω /π 27.00.0012 L. Ahle, et al., E802 Collaboration, Phys. Rev. C 60 (1999) 064901. [6] L. Ahle, et al., E866/E917 Collaboration, Phys. Lett. B 490 In conclusion, while the statistical model cannot (2000) 53. + + [7] S. Albergo, et al., Phys. Rev. Lett. 88 (2002) 062301. explain the sharpness of the peak in the K /π [8] S. Ahmad, et al., Phys. Lett. B 381 (1996) 3. ratio, its position corresponds precisely to a transi- [9] J. Klay, et al., E895 Collaboration, Phys. Rev. C 68 (2003) tion from a baryon-dominated to a meson-dominated 054905. hadronic gas. This transition occurs at a temperature [10] M. Ga´zdzicki, M.I. Gorenstein, Acta Phys. Pol. B 30 (1999) = = 2705. T 140 MeV, a baryon√ chemical potential µB = [11] R. Stock, J. Phys. G: Nucl. Part. Phys. 30 (2004) S633. 410 MeV and an energy sNN 8.2 GeV. In the [12] M.G. Alford, K. Rajagopal, F. Wilczek, Phys. Lett. B 442 statistical model this transition leads to a sharp peak (1998) 247; in the Λ/π ratio, and to moderate peaks in the Y.H. Atta, T. Ikeda, Phys. Rev. D 67 (2003) 014028. K+/π +, Ξ −/π + and Ω−/π + ratios. Furthermore, [13] Z. Fodor, S.D. Katz, J. High Energy Phys. 0203 (2002) 014; these peaks are at different energies in the statistical Z. Fodor, S.D. Katz, hep-lat/0402006. [14] F. Karsch, K. Redlich, A. Tawfik, Phys. Lett. B 571 (2003) 67. model. The statistical model predicts that the maxima − + − + [15] J. Cleymans, K. Redlich, Phys. Rev. Lett. 81 (1998) 5284. in the Λ/π, Ξ /π and Ω /π occur at increas- [16] J. Cleymans, H. Oeschler, K. Redlich, Phys. Rev. C 59 (1999) ing beam energies. 1663. 54 J. Cleymans et al. / Physics Letters B 615 (2005) 50–54 [17] F. Becattini, J. Cleymans, A. Keränen, E. Suhonen, K. Redlich, [20] F. Becattini, M. Ga´zdzicki, A. Keränen, J. Manninen, R. Stock, Phys. Rev. C 64 (2001) 024901. Phys. Rev. C 69 (2004) 024905. [18] P. Braun-Munzinger, J. Cleymans, H. Oeschler, K. Redlich, [21] S. Wheaton, J. Cleymans, hep-ph/0407174. Nucl. Phys. A 697 (2002) 902. [22] C. Meurer, et al., J. Phys. G 30 (2004) S14; [19] P. Braun-Munzinger, K. Redlich, J. Stachel, in: R. Hwa, X.N. S.V. Afanasiev, et al., Phys. Lett. B 358 (2002) 275. Wang (Eds.), Quark–Gluon Plasma 3, World Scientific, Singa- [23] P. Chung, et al., Phys. Rev. Lett. 91 (2003) 202301. pore, 2003, nucl-th/0304013. [24] C. Alt, et al., nucl-ex/0409004. Physics Letters B 615 (2005) 55–60 www.elsevier.com/locate/physletb Absence of structure in the 20,22Ne + 118Sn quasi-elastic barrier distribution E. Piasecki a,Ł.Swiderski´ a,P.Czosnykaa, M. Kowalczyk a,K.Piaseckia, M. Witecki a, T. Czosnyka b,J.Jastrz˛ebski b, A. Kordyasz b, M. Kisielinski´ b, T. Krogulski c, M. Mutterer d, S. Khlebnikov e, W.H. Trzaska f, K. Hagino g,N.Rowleyh a Institute of Experimental Physics, Warsaw University, Poland b Heavy Ion Laboratory, Warsaw University, Poland c Institute of Experimental Physics, University in Bialystok, Poland d Institut für Kernphysik, Technische Universität, Darmstadt, Germany e Radium Institute, St. Petersburg, Russia f University of Jyväskylä, Finland g Tohoku University, Sendai, Japan h Institut de Recherches Subatomiques/Universite Louis Pasteur (UMR 7500), Strasbourg, France Received 15 October 2004; received in revised form 4 February 2005; accepted 24 March 2005 Available online 7 April 2005 Editor: V. Metag Abstract Motivated by the extreme deformation parameters of the projectile, we have measured quasi-elastic scattering for 20Ne + 118Sn. In contrast to calculations based on known collective states, the experimental barrier distribution is structureless. A com- parison with the system 22Ne + 118Sn shows that this smoothing is unlikely to be due to nucleon- or α-transfer channels, and is more likely to be due to coupling to many other weak channels. 2005 Elsevier B.V. All rights reserved. PACS: 25.70.Bc; 25.70.Hi; 25.70.Jj Keywords: Coupled channels; Coulomb barrier distribution; 20Ne; 22Ne; Quasi-elastic scattering Our understanding of the interplay of nuclear reac- of barrier distributions, which show the remarkable ef- tion channels has been greatly enhanced by the study fects of strong coupling to collective excitations. Such work has mainly taken place in the context of fusion E-mail address: [email protected] (E. Piasecki). [1], where the structures are observed in the “experi- 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.03.070 56 E. Piasecki et al. / Physics Letters B 615 (2005) 55–60 mental barrier distribution” [2] coupling to collective states gives rise to a struc- 2 ture in the barriers that must be crossed if the d (E σfus) Dfus(E) = . (1) target and projectile nuclei are to come into close dE2 contact. However, after crossing the barriers, the Here σfus is the fusion cross section and E is the inci- system may quasi-fission (binary separation be- dent centre-of-mass energy. In many cases, this func- fore forming an equilibrated compact CN) or form tion gives a clear fingerprint of the couplings involved a compound nucleus that may then yield evapora- (rotational or vibrational states in the target and/or pro- tion residues or may itself fission. This is partic- jectile in varied degrees of complexity). In some cases ularly true of very heavy systems leading to the multi-nucleon transfer channels may also play an im- creation of superheavy elements. Thus three cross portant role [3–5]. sections must be measured, two of which (fusion– One can also define barrier distributions in terms of fission and quasi-fission) yield angular distribu- the quasi-elastic cross section, that is the sum of elas- tions which must be integrated. Thus a complex tic, inelastic and transfer channels [6]. This has several experiment is necessary if one wishes to exploit advantages: Eq. (1). However, quasi-elastic scattering is the complementary part of the barrier-crossing flux • The quasi-elastic barrier distribution is defined in (and related to it by unitarity) and in principle terms of a first derivative rather than a second there are no additional complications in the quasi- derivative: elastic distribution for such systems, though there is of course the practical consideration of discrim- d(σqe/σR) D =− . (2) inating against fission fragments in the detectors. qe dE This leads to a smaller statistical error in Dqe com- Because of the above potentially interesting appli- pared with Dfus for a given statistical error in the cations it is important to understand as well as possible corresponding cross section. This could make the the relationship between Dqe and Dfus. use of Dqe interesting for low-intensity, for exam- Despite the above-mentioned advantages, there is ple radioactive, beams. a possible disadvantage to Dqe. While theoretical cal- • Eq. (2) is most readily derived and understood for culations generally give very similar results for Dqe = a scattering angle θ π, though we can apply the and Dfus [6,7], it has been noted that the experimental formula for any large scattering angle by intro- quasielastic distribution, while having the same over- ducing an effective energy into the cross section: all shape as the fusion distribution, may possess less σqe(Eeff) ≈ σqe(E, θ) with [6] detailed structure. This disagreement was observed in 16 + 144 2 Ref. [6] for the system O Sm. Here, a small sec- E = E. (3) ondary peak (arising essentially from coupling to the eff 1 + cosec(θ/2) lowest octupole phonon state in the target) is present This corrects for centrifugal effects and means in Dfus but not in Dqe. On the other hand, for the sys- 40 90 that σqe can be measured simultaneously at several tem Ca + Zr, similar phonon structures are present effective energies with a single beam energy by in both distributions [6]. The reason for this difference using detectors disposed at different scattering an- is not clear. gles. Along with the use of degraders, this greatly The systems 16O + 116,119Sn have already been facilitates studies with a cyclotron accelerator [8]. measured at the Warsaw cyclotron [8]. For those sys- • The measurement of the total quasi-elastic cross tems, σfus has also been measured [9] and neither Dqe section means that relatively simple detectors can nor Dfus display any significant structure, either theo- be exploited since there is no need for a great res- retically or experimentally. olution either in mass, charge or energy. Our aim in the present work is to try to elucidate • For systems leading to a heavy compound nu- the above problem by taking a system where the pro- cleus (CN), Eq. (1) should be written for the total jectile has a known large deformation that leads to a “barrier-passing” cross section. In other words the prediction of very marked structure in both Dqe and E. Piasecki et al. / Physics Letters B 615 (2005) 55–60 57 Dfus and to measure its quasi-elastic scattering. We have chosen the system 20Ne+ 118Sn for the very large deformation of the projectile: β2 = 0.46, β4 = 0.27. These values seem to be representative of results of many scattering experiments with projectiles like elec- trons, protons, pions and α-particles [10]. The experimental method was very similar to that described in Ref. [8]. That is, we measured quasielas- tic large-angle scattering excitation functions using 30 10×10 mm PIN diodes placed at 130◦, 140◦ and 150◦ in the laboratory system and two “Rutherford” semi- conductor detectors (of 6 mm diameter) placed at 35◦ with respect to the beam. The experiment consisted of measuring the backward/forward counting ratio as a function of the projectile energy. The geometry of the Fig. 1. Examples of the Q-value spectra measured for 20,22Ne scat- ◦ ◦ PIN diode set-up (including the detector–target sepa- tered at 140 and 35 on 118Sn for energy corresponding to max- rations) was the same as in Ref. [8]. In addition, 4 tele- imum of the barriers. The calculated distribution is presented by a scopes at backward angles (two each at 120◦ and 160◦) dashed line. For easier comparison the spectra have been normalized to a common peak value. were used to measure the relative intensity of light charged particle transfers. The light charged particle transfer was not measured at forward angles, however, tion, the stability of the gain, the offset and resolution the CC calculations indicate that at these energies it is of the electronics and detectors as well as the beam- strongly backward peaked. energy were continuously monitored using a precise The 20Ne beam with intensity of 1–40 pnA (de- pulse generator and an α-particle source. pending on the energy) from the Warsaw Cyclotron Data analysis was performed as in Ref. [8]: from impinged on the 130 µg/cm2 target of 118Sn (enriched the kinetic energy spectra, assuming two-body kine- 2 matics, we calculate the Q-value spectrum (see Fig. 1). to 96.6%) on a 65 µg/cm Al2O3 backing, which was taken into account in data analysis. Then the number of counts in the integration limits − The use of nickel foils as energy degraders and the ( 3, 15) MeV for the forward detectors and − different detectors angles enable measurements at ef- ( 3, 22) MeV for the backward ones are obtained, E is calculated and the σ /σ ratios constructed. fective energies Eeff with sufficiently small steps. En- eff qe R ergy calibration was performed using a precise pulse Data were binned over 0.75 MeV intervals and a generator and Bi–Po α-particle sources (the estimated Savitzky–Golay filter (cubic, 5 points) was applied. pulse height defect of the semiconductor detectors was Then, using Eq. (2), the barrier distributions were de- negligibly small). termined. The energy resolution of the experiment is of crit- The telescopes enabled us to check that the projec- ical importance for assessing the meaning of the re- tile-like fragments following elastic, inelastic, p- and sults. It was continuously monitored during the exper- α-transfer reactions were contained within the integra- iment and turned out to vary in time within the limits tion limits, while protons and α-particles (for example, 0.8–1.6 MeV (FWHM). Apart from the contributions evaporated from fusion products) were rejected. Nei- coming from straggling in the degrader and the energy ther fission nor quasi-fission fragments were observed. loss in the target (≈ 0.3 MeV each), the main contri- The possible products of n-transfer reaction, in our bution to the width comes from characteristics of the measurement not experimentally distinguishable from beam. Precise information on the beam-energy spread the inelastic excitations, are in the same energy range was determined from the energy spectra registered in and are also necessarily included in the quasi-elastic the forward detectors. Other effects, for example, the events. geometry of detection, had very little influence on the The experimental results for σqe/σR and Dqe are energy resolution. In addition to checking the resolu- shown in Fig. 2(a) and (b) (filled circles). Also shown 58 E. Piasecki et al. / Physics Letters B 615 (2005) 55–60 20 118 Fig. 3. (a) Experimental Dqe for Ne + Sn compared with the theoretical distribution (solid line). The dotted line shows the no-coupling result. See text. Both calculated distributions are folded with the experimental resolution 1.6 MeV. The dashed line shows 22 118 the results without folding. (b) Experimental Dqe for Ne + Sn compared with the theoretical calculation folded with a 2.0 MeV resolution (solid line), and ignoring the resolution (dashed line). Fig. 2. (a) Experimental values of σqe/σR as a function of Eeff (see text) for the systems 20Ne + 118Sn (filled circles) and 22Ne + 118Sn (open circles). (b) The derived barrier distributions Dqe. essentially shifts the entire distribution to lower ener- gies (effectively a renormalisation of the central po- are the results for a 22Ne beam (open circles) which tential). Here the shift is 0.9 MeV. In addition, we will be discussed below. We note that there is no sig- included couplings to the lowest vibrational state in ∗ = = nificant structure in either of the Dqe distributions. the target (E 1.41 MeV; β2 0.11). Also shown in This is in evident disagreement with theoretical cal- Fig. 3(a) is the no-coupling calculation (dotted curve). culations which we shall present below. We see that the coupled channels calculations result in In Fig. 3(a), the solid line is the result of a coupled- a Dqe distribution which possesses two distinct peaks channels (CC) calculation performed using the pro- in contrast to the experimental curve which is rela- gram CCQUEL [11]. (Unless otherwise stated, the tively smooth and centered around 62 MeV. theoretical curves take into account the experimental One should emphasize that this structure cannot be + resolution.) This calculation includes the states 0 , removed by changing β2 and β4 values within the ob- 2+,4+,6+ in the 20Ne ground state rotational band. served scatter quoted in Ref. [10]. Calculations were The results converge rapidly as the number of states performed with inclusion of reorientation terms but is increased. It was verified that truncation of calcula- their omission affects only the details of the distrib- tions at the 6+ level is entirely sufficient for our pur- ution and not the presence of the structure. The two- poses. In addition, this calculation takes account of the peaked structure arising from the 20Ne rotational band strong octupole phonon state in the projectile (E∗ = differs from the smooth barrier distribution obtained 154 5.62 MeV; β3 = 0.39). As observed in other cases for heavy deformed nuclei (for example, Sm [13]) [12], however, the presence of a high-lying phonon due to the large excitation energies in the present case E. Piasecki et al. / Physics Letters B 615 (2005) 55–60 59 (E(2+) = 1634 keV as opposed to 82 keV). These tems 40Ca + 90,96Zr are very different from each other reduce considerably the influence of the higher-lying and it has been suggested [5,6,15] that this is due to the members of the rotational band and force a rapid presence of strong multi-nucleon transfer channels [4] convergence to a barrier distribution possessing pro- for the 96Zr case (6 neutrons outside the N = 50 closed nounced structure. We note, however, that our calcu- shell with favorable Q values). (Though note that this lations seem to give the correct overall width of the explanation is not universally accepted [16].) Unfor- barrier distribution. tunately it is very difficult to include transfer channels It is important to mention that the correctness of in the CC calculations and we chose, therefore, to do the experimental procedure and of the resolution deter- an additional experiment using a 22Ne beam. Here mination was confirmed a posteriori through the bar- we have much smaller effective Q values [17] than rier distribution measurements for 20Ne+ natNi, which for 20Ne (see Table 1) and this should give rise to a were performed simultaneously with the present mea- reduced probability of α transfer. Indeed, our experi- surements. Despite the fact that the Ni target was a nat- ment confirmed that while in the case of 20Ne beam ural mixture of isotopes and despite its smaller atomic at the barrier energy the α-transfer contributes 3% of number (which would normally be expected to give the total scattering cross section at backward angles rise to a weaker barrier structure), a clear structure, in (120◦–160◦) and much less below the barrier, in the agreement with calculations, was found [14]. case of 22Ne it is weaker by a factor of about 6. The Discussing the possible reasons of disagreement contribution of proton transfer is still smaller. between experimental and calculated barrier distribu- The calculated curve in Fig. 3(b) takes the values tions, we would like to draw attention to Fig. 1.We β2 = 0.47 and β4 = 0.10 [10]. Although the defor- note there that while the measured large-angle nega- mation parameter β4 is considerably smaller than that tive Q-value spectrum extends to quite high energies, for 20Ne, the calculation predicts a similar though less the calculated spectrum fails to reproduce anything pronounced structure in the Dqe distribution. However, like the observed intensities in that region. This indi- the experimental results for Dqe (Fig. 1(b)) show that cates the importance of high-lying excitations in ei- for 22Ne this function is still very smooth, suggesting ther the target or projectile, which may in some way that the absence of structure in the 20Ne case is not wash out the predicted structures. However, cutting out due to strong α transfers. This is suggested also by the these high Q-value contributions to the quasi-elastic similarity of the Q spectra measured with both neon cross section has rather little effects on the distribution isotopes (Fig. 1). of barriers. It seems, therefore, that it is not simply Can the barrier distribution be smoothed out by the presence of extra contributions to the cross section neutron transfer? Indeed, in the case of 16O + 144Sm, which destroy the structure in the barrier distribution, the disappearance of the high-energy peak [6] was at- but rather the way they interfere with (damp) the cou- tributed to the dominance of the contribution of neu- plings to the collective states. tron transfer channels to σqe at that energy (where Can these high-lying excitations be caused by σqe/σR is small). In our case, however, the structures 20 transfer reactions? Since Ne is a good α-cluster lie in an energy range where σqe/σR is still of the order nucleus, it could be that α-transfers affect the above of 0.5–0.8. We have no measure of the neutron-transfer results by providing more channels, more barriers and cross sections, but judging from the negative Q values thus a greater overlap and smoothing of the structures. (see Table 1) they are not expected to be important ei- In fact it has been observed that the Dfus for the sys- ther. Table 1 118 * Effective Q values for stripping in the Ne + Sn collisions (in parentheses the Qgg values are given) 1n 2n 1p 2p α 20Ne −9.6 (−10.4) −11.2 (−12.9) −2.1 (−7.7) +2.8 (−8.5) +9.5 (−3.6) 22Ne −3.2 (−3.9) 0.0 (−1.5) −4.4 (−10.2) −2.8 (−14.1) +4.2 (−8.6) * For pick-up reactions the effective Q values are all negative, with the exception of 2n-transfer, where Q =+0.8 MeV. In this experiment pick-up of charged particles was never observed. 60 E. Piasecki et al. / Physics Letters B 615 (2005) 55–60 Thus, for the moment we do not consider the pos- fact that present results confront us with clear limi- sible influence of transfer reactions on the Dqe distri- tations of the standard CC approach. Experimentally, bution as very likely. investigations of this could consist from measurements It is appropriate to mention that since the high en- of Dfus for present systems and looking for other ones. ergy excitations were observed also for the 16O + ASn Theoretically, one should look for possibilities of re- systems [8], while it was much weaker for the 20Ne + laxing the present limitations on the model space. natNi system, where the predicted structure was ob- served [14], it is likely that this high-lying intensity corresponds mainly to excitations of the 118Sn tar- Acknowledgements get rather than to excitations of 20Ne projectile. We have, however, already noted that high-lying collec- We are grateful to Marie-Antoinette Saettel for her tive states simply shift the entire barrier distribution to help in providing the targets for this experiment. The lower energies, and so it might be expected that this work was funded in part by Grant No. 2 P03B 02624 of intensity corresponds to weak coupling to many non- the KBN and also supported by the cooperation agree- collective states. ment between the IN2P3 (France) and Polish Labora- tories. In this context it would be valuable to measure σfus 20,22 118 for Ne + Sn and to find out if Dfus is similar to Dqe or whether it still possess the predicted structure. If the predicted structures were observed, this could References again point to the influence of many weak channels, [1] M. Dasgupta, D.J. Hinde, N. Rowley, A.M. Stefanini, Annu. which may contribute differently in fusion and scatter- Rev. Nucl. Part. Sci. 48 (1998) 401. ing. This would mean that principal assumption of the [2] N. Rowley, G.R. Satchler, P.H. Stelson, Phys. Lett. B 254 standard CC model about truncation of available chan- (1991) 25. nels to just a few, strongly coupled collective ones, is [3] H. Timmers, et al., Phys. Lett. B 399 (1997) 35. too radical for the description of scattering of 20,22Ne [4] G. Montagnoli, et al., J. Phys. G: Nucl. Part. Phys. 23 (1997) A 1439. on Sn nuclei. If they were not, this would mean that [5] N. Rowley, in: Yu.Ts. Oganessian, V. Zagrebaev (Eds.), Proc. the model space is too limited even for description of of International Conference on Fusion Dynamics at the Ex- fusion. tremes, Dubna, 25–27 May, 2000, World Scientific, Singapore, This rises a number of interesting questions. Can 2001, p. 296. [6] H. Timmers, et al., Nucl. Phys. A 584 (1995) 190. we simulate the smoothing action of the weak chan- [7] K. Hagino, N. Rowley, Phys. Rev. C 69 (2004) 054610. nels by schematically including some extra couplings, [8] E. Piasecki, et al., Phys. Rev. C 65 (2002) 054611. recovering the missing strength in high Q-value bins? [9] V. Tripathi, et al., Phys. Rev. C 65 (2001) 014614. Is it possible to do it in more realistic way? Is the im- [10] G.S. Blanpied, et al., Phys. Rev. C 38 (1988) 2180. portance of weak channels a reason for similar effects [11] K. Hagino, N. Rowley, unpublished. [12] K. Hagino, et al., Phys. Rev. Lett. 79 (1997) 2014. observed with other systems [6]? How to recognize [13] J.X. Wei, et al., Phys. Rev. Lett. 67 (1991) 3368. systems not posing such problems? [14] L. Swiderski,´ et al., Int. J. Mod. Phys., in press. Concluding, the usefulness of the quasi-elastic scat- [15] V.I. Zagrebaev, Phys. Rev. C 67 (2003) 061601. tering for the barrier distribution investigations de- [16] G. Montagnoli, et al., Eur. Phys. J. A 15 (2002) 351. pends on a deeper understanding of some of these [17] R. Bass, Nuclear Reactions with Heavy Ions, Springer-Verlag, Berlin, 1980, 152. points. However, perhaps even more important is the Physics Letters B 615 (2005) 61–67 www.elsevier.com/locate/physletb Influence of NN-rescattering effect on the photon asymmetry of d(γ,π −)pp reaction Eed M. Darwish Physics Department, Faculty of Science, South Valley University, Sohag 82524, Egypt Received 14 December 2004; received in revised form 28 January 2005; accepted 6 April 2005 Available online 15 April 2005 Editor: J.-P. Blaizot Abstract − The influence of final-state NN-rescattering on the beam asymmetry Σ for linearly polarized photons in π photoproduc- tion on the deuteron in the energy range from π-threshold through the ∆(1232)-resonance has been investigated. Numerical results for this spin observable are predicted and compared with recent experimental data from the LEGS Spin Collaboration. Final-state NN-rescattering is found to be quite important and leads to a better agreement with existing experimental data. Furthermore, the differences with other theoretical models have been discussed. 2005 Elsevier B.V. All rights reserved. PACS: 24.70.+s; 14.20.-c; 29.27.Hj; 25.30.Fj Keywords: Polarization phenomena in reactions; Spin observables; Polarized beams; Final-state interactions 1. Introduction fact that in contrast to the differential cross section, which is a sum of the absolute squares of the t-matrix It is a well-known fact that polarization observables elements, these polarization observables contain inter- allow a further and much more detailed analysis of ference terms of the various reaction amplitudes in the process under study compared to the differential different combinations and, therefore, may be more cross section alone. Because polarization observables sensitive to small amplitudes and to small contribu- contain a much richer information on the dynamics of tions of interesting dynamical effects. the system than attainable without beam and/or target In recent years a great effort, both from theoret- polarization and without polarization analysis of the ical [1–7] and experimental [8–11] points of view, particles in the final state. The reason for this is the has been devoted to the analysis of single-pion photo- production with polarized beams and/or polarized tar- gets. In [1], π − photoproduction on the deuteron has E-mail address: [email protected] (E.M. Darwish). been studied within a diagrammatic approach includ- 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.04.018 62 E.M. Darwish / Physics Letters B 615 (2005) 61–67 ing nucleon–nucleon (NN) and pion–nucleon (πN) can only be considered as a first step towards a more rescattering in the final state. Special emphasize was realistic description of spin observables. given for the analyzing powers connected to beam In this Letter we investigate, therefore, the in- and target polarization, and to polarization of one fluence of final-state NN interaction on the photon of the final protons. First preliminary model calcula- asymmetry for the reaction d(γ,π −)pp in the en- tions for the photon asymmetry Σ of the γd → π −pp ergy region from π-threshold through the ∆(1232)- reaction have been given in the pure impulse ap- resonance. To our knowledge, the influence of NN- proximation (IA) [2]. The comparison between these FSI effect on this spin observable has never been predictions for Σ and the preliminary experimen- studied before. The πN-rescattering contribution has tal data from LEGS Spin Collaboration [12] gives a been considered as negligible in the region of the clear indication that the effects of final-state inter- ∆(1232)-resonance [14,15] and thus it is not consid- action (FSI) may be important. The deuteron tensor ered in the present work. Our main goal is to analyze analyzing powers of the reaction d(γ,π −)pp have the recent experimental data from LEGS [12]. Further- been studied in the IA [3] without inclusion of any more, it was an open question whether the inclusion FSI or two-body exchange current contributions. In of rescattering contributions would lead to a good de- our previous papers [4–6], various polarization ob- scription of the available data. servables in inclusive single-pion photoproduction on In the next section we will define the photon asym- the deuteron using a polarized photon beam and/or metry Σ in terms of the transition matrix amplitude. an oriented deuteron target have been investigated in In Section 3 we will present and discuss the numeri- the pure IA only, i.e., by neglecting any FSI effects cal results of our calculations and compare them with and possible two-body contributions to the produc- the experimental data and other predictions. Finally, tion operator. In particular, a complete survey on all we summarize our conclusions in Section 4. single- and double-polarization observables like beam and target asymmetries was given. In [7] the influence of final-state NN-rescattering on the helicity struc- 2. Linear photon asymmetry ture of the inclusive reaction γ d → π −pp has been investigated. The differential polarized cross-section difference for the parallel and antiparallel helicity The beam asymmetry Σ for linearly polarized pho- states has been predicted and compared with recent tons is defined in analogy to deuteron photodisinte- experimental data from MAMI (Mainz/Pavia) [13].It gration [16] writing the differential cross section for has been shown that the effect of NN-rescattering linearly polarized photons and unpolarized deuterons is much less important in the polarized differen- in the form tial cross-section difference than in the unpolarized dσ dσ0 one. (θ ,φ ) = (θ ) dΩ π π dΩ π The photon asymmetry Σ is very sensitive to the in- π π × + γ ternal mechanisms of the reaction and, therefore, can 1 P Σ(θπ ) cos 2φπ , (1) be a very useful test to impose constraints on the theo- retical models. The work has been partly motivated by where dσ0/dΩπ denotes the semi-inclusive unpolar- preliminary experimental results, for the γd → π −pp ized differential cross section of incoherent pion pho- toproduction on the deuteron, where only the final channel, with the LEGS Brookhaven National Labo- γ ratory [12] which shows strong and not trivial angu- pion is detected without analyzing its energy [14], P lar dependences of this observable. In agreement with is the degree of linearly polarized photons [16], θπ and these preliminary data, one can see in [2,6] that the φπ represent the polar and azimuthal pion angles and predictions in the pure IA can hardly provide a reason- Σ is the photon asymmetry for linearly polarized pho- able description of the data since major discrepancies tons. Then one has [4,16] are found. As already noted in [2,6], the effect of NN- rescattering is quite important. This means in particu- dσ0 Σ =−W00, (2) lar that the calculation in the spectator nucleon model dΩπ E.M. Darwish / Physics Letters B 615 (2005) 61–67 63 with for NN-rescattering in the final state. Further details with respect to the matrix elements are not discussed 1 1−m 110 W00 = √ (−) d C md −m 0 here and can be found in [14]. 2 3 d smt,mγ md md qmax × dq dΩ ρ M(tµ) 3. Results and discussion pNN s sm,mγ md 0 Here we present and discuss our results for the pho- (tµ) × M (3) ton asymmetry Σ for linearly polarized photons of s−m,mγ −m d π − photoproduction on the deuteron with inclusion j1j2j denoting with Cm1m2m a Clebsch–Gordan coefficient, of NN-rescattering in the final state. These results mγ the photon polarization, md the spin projection are also compared to the preliminary experimental of the deuteron, s and m total spin and its projection data of LEGS [12] and the preliminary IA calculations of the two outgoing nucleons, respectively, t their to- of Lee [2]. The results presented here are calculated tal isospin, µ the isospin projection of the pion, qmax by using the elementary photoproduction operator of the maximum value of pion momentum, ΩpNN the Schmidt et al. [17] and the deuteron wave function solid angle of the relative momentum p NN of the final of Paris potential [18]. For the half-off-shell NN- NN system and ρs the phase space factor. For further scattering amplitude, the separable representation [19] details with respect to the kinematical variables and of the realistic Paris potential has been used. All par- quantum numbers we refer to our previous work [14]. tial waves with total angular momentum J 3are For the transition M-matrix we include, in this included. work, besides the pure IA, the driving term from NN- We start with presenting our results for the linear rescattering, so that the total transition matrix reads photon asymmetry Σ at different photon lab-energies = M(tµ) = M(tµ) IA + M(tµ) NN ωγ 200, 270, 330, 370, 420 and 500 MeV in Fig. 1 smm m smm m smm m , (4) γ d γ d γ d as a function of emission pion angle θπ in the lab- where the first term represents the transition amplitude oratory frame. The solid curves show the results of in the pure IA and the second is the corresponding one the full calculation, i.e., when NN-rescattering is in- − Fig. 1. Linear photon asymmetry of the differential cross section for linearly polarized photons for π photoproduction on the deuteron as a function of emission pion angle in the laboratory frame at fixed values of photon lab-energies. Notation of curves: dashed—IA; solid—IA + NN-rescattering. 64 E.M. Darwish / Physics Letters B 615 (2005) 61–67 + Fig. 2. The ratio of the linear photon asymmetry with NN-rescattering ΣIA NN to the one in the IA ΣIA as a function of emission pion angle at the same photon lab-energies as in Fig. 1. cluded, while the dashed curves show the contribu- Fig. 3 shows the sensitivity of our results for the tion of the IA alone in order to clarify the importance linear photon asymmetry Σ to the photon lab-energy ◦ ◦ of NN-FSI effect. In order to show in greater detail ωγ at three fixed values of pion angle θπ = 0 ,90 the relative influence of NN-rescattering effect on the and 180◦ for photon lab-energies between 200 and linear photon asymmetry, we show in Fig. 2 the ef- 500 MeV. In order to show the relative influence of fect of NN-rescattering relative to the IA by the ratio NN-rescattering effect on the linear photon asymme- ΣIA+NN/ΣIA, where ΣIA denotes the photon asym- try, we show in the bottom panels of Fig. 2 the ef- metryintheIA and ΣIA+NN the one including the fect of NN-rescattering relative to the IA by the ratio contribution of NN-rescattering. ΣIA+NN/ΣIA, where ΣIA denotes the photon asym- In the photon energy domain of this work, the mag- metry in the IA and ΣIA+NN the one including the netic multipoles dominate over the electric ones, due contribution of NN-rescattering. In view of these re- to the excitation of the ∆-resonance. This is clear sults, one notes that NN-rescattering—the difference from the dominantly negative values of Σ as shown in between the solid and the dashed curves—is quite Fig. 1. On the contrary, the left-top and right-bottom small, almost completely negligible at extreme for- panels in Fig. 1 show that small positive values are ward and backward angles. found at ωγ = 200 and 500 MeV. We see also that the Fig. 4 shows a comparison of our numerical re- asymmetry Σ is sensitive to the energy of the incom- sults for the linear photon asymmetry Σ in the pure ing photon. It is noticeable, that the photon asymmetry IA (dashed curves) and with NN-rescattering (solid ◦ ◦ Σ vanish at θπ = 0 which is not the case at 180 .At curves) with experimental data. In view of the fact extreme forward and backward emission pion angles that experimental data for this spin observable are one sees, that the photon asymmetry is relatively small not available in a final form, we compare our predic- in comparison to the results when θπ changes from tions with the preliminary experimental data from the about 30◦ to 120◦. One notices also, that the contribu- LEGS Spin Collaboration [12] as depicted in Fig. 4. tion from NN-rescattering is much important in this We see that the general feature of the data is repro- region, in particular in the peak position. For lower and duced. However, the discrepancy is rather significant higher photon energies, one finds the strongest effect in the region where the photon energy close to the by NN-rescattering. ∆-resonance. This could be due to the higher-order E.M. Darwish / Physics Letters B 615 (2005) 61–67 65 − Fig. 3. Linear photon asymmetry Σ for d(γ,π )pp as a function of the photon lab-energy ωγ at fixed values of pion angle θπ in the laboratory frame with NN-rescattering (top panels) and their ratios with respect to the pure IA (bottom panels). Notation as in Fig. 1. − Fig. 4. Photon asymmetry Σ for the reaction d(γ,π )pp at ωγ = 270 and 330 MeV photon lab-energy in comparison with the preliminary data from LEGS [12]. Notation as in Fig. 1. rescattering mechanisms which are neglected in this ferences between theory and experiment are still ev- work. In the same figure, we also show the results from ident. It is appear that our model is still not capable the IA only (dashed curves). It is seen that the NN- of describing the measured photon asymmetry, even rescattering yields an about 10% effect in the region if NN-FSI is included. Future efforts must be made of the peak position. We found that this is mainly due to remove the remaining discrepancies such as a com- to the interference between the IA amplitude and the plete three-body treatment of the final πNN system. NN-FSI amplitude. In agreement with our previous Now we compare our results for the linear pho- results [6], one notes that the pure IA (dashed curves ton asymmetry with the preliminary model predic- in Fig. 4) cannot describe the experimental data. The tions of Lee [2] as shown in Fig. 5. The solid curves inclusion of NN-FSI leads at ωγ = 270 MeV to a show the results of the present calculations when NN- quite satisfactory description of the data, whereas at rescattering is included and the dashed ones show our 330 MeV NN-FSI effect is small and therefore dif- results in the IA alone. The preliminary IA results 66 E.M. Darwish / Physics Letters B 615 (2005) 61–67 − Fig. 5. Photon polarization asymmetry Σ for d(γ,π )pp reaction at two photon lab-energies. Dashed curves: IA of present calculations; solid curves: IA plus NN-FSI of present calculations; dotted curves: preliminary IA calculations of Lee [2]. of Lee [2] are represented in this figure by the dot- well as with the ‘totally neglected’ rescattering mech- ted curves. It is very clear that the results for the IA anisms in [2] would be very interesting to understand of the present work showed certain significant differ- the origin of this discrepancy. ences to the preliminary IA results [2] which cannot be attributed to the use of different elementary pion photoproduction operators and/or from different real- 4. Summary istic NN potential models used for the deuteron wave function. It is clear from the right panel of Fig. 5 that In this Letter we have investigated the influence of the discrepancy is rather significant in the region of the final-state NN-rescattering effect on the linear pho- ∆(1232)-resonance, in particular, in the peak position. ton asymmetry Σ for the γd → π −pp reaction in the As already mentioned in the beginning of this sec- photon energy range from π-threshold to 500 MeV. tion, our results are calculated using the effective La- We have found that the effect due to the final two- grangian model developed by Schmidt et al. [17]. nucleon interaction to be small, but it can has signif- The main advantage of this model is that it has been icant contribution to the photon asymmetry through constructed to give a realistic description of the ∆- its interference with the dominant term from the im- resonance region. It is also given in an arbitrary frame pulse approximation. Furthermore, the linear photon of reference and allows a well defined off-shell contin- asymmetry is found to be sensitive to the energy of the uation as required for studying pion production on nu- incident photon. In comparison with the preliminary clei. As shown in Figs. 1–3 in [14], the results for this experimental data from LEGS [12], the inclusion of model are in good agreement with recent experimental NN-FSI effect leads to a better agreement with exper- data as well as with other theoretical predictions. On imental data. With respect to the comparison with the the other hand, the well-known dynamical model of preliminary results of [2] in the IA, we found a large Sato and Lee [20] has been used in [2]. This model has difference between both calculations in the peak posi- given also a successful description of the pion photo- tion. The origin of this difference is still not clear. production data. Therefore, the big difference between We would like to conclude that the results presented both predictions in the IA results cannot be attributed here for linear photon asymmetry of d(γ,π −)pp can to the use of different elementary operators. This can be used as a basis for the simulation of the behav- be interpreted as lack of understanding of the nonres- iour of this asymmetry and for an optimal planning of onant background, which in dynamical models is re- new experiments of this reaction with polarized pho- lated to the pion cloud. It seems that pion cloud effects ton beams. An experimental check of these predictions are not yet consistently included in dynamical mod- for the linear photon asymmetry covering a large range els. An independent evaluation for Σ in the pure IA as for the pion angle would provide an additional sig- E.M. Darwish / Physics Letters B 615 (2005) 61–67 67 nificant test of our present understanding of this spin [4] E.M. Darwish, Nucl. Phys. A 735 (2004) 200. observable. Furthermore, an independent evaluation in [5] E.M. Darwish, Int. J. Mod. Phys. E 13 (2004) 1191. the framework of effective field theory would be very [6] E.M. Darwish, J. Phys. G: Nucl. Part. Phys. 31 (2005) 105. [7] E.M. Darwish, Nucl. Phys. A 748 (2005) 596. interesting. Future improvements should include fur- [8] St. Goertz, W. Meyer, G. Reicherz, Prog. Part. Nucl. Phys. 49 ther investigations including FSI as well as two-body (2002) 403; effects. St. Goertz, W. Meyer, G. Reicherz, Prog. Part. Nucl. Phys. 51 (2003) 309, Erratum. [9] See, for example, in: M. Anghinolfi, M. Battaglieri, R. De Acknowledgements Vita (Eds.), Proceedings of the 2nd International Symposium on the Gerasimov–Drell–Hearn 2002, Sum Rule and the Spin Structure of the Nucleon, Genova, Italy, 3–6 July 2002, World This work is supported in part by the Biblio- Scientific, Singapore, 2003. theca Alexandrina—Center for Special Studies and [10] B. Krusche, S. Schadmand, Prog. Part. Nucl. Phys. 51 (2003) Programs—under grant No. 2602314 Sohag 2nd- 399. Sohag. I am indebted to Prof. H. Arenhövel as well [11] V. Burkert, T.-S.H. Lee, Int. J. Mod. Phys. E 13 (2004) 1035. [12] A. Sandorfi, M. Lucas, private communication; as the members of his work group for fruitful discus- M. Lucas, in: LOWq Workshop on Electromagnetic Nuclear sions and valuable comments. I would like to thank Reactions at Low Momentum Transfer, 23–25 August 2001, Profs. T.-S. Harry Lee, T. Sato for useful discussions Halifax, Nova Scotia, Canada, 2001; and Prof. S.A.E. Khallaf for a careful reading of this A. Sandorfi in [9]. Letter. [13] P. Pedroni, private communication; C.A. Rovelli, Diploma Thesis, University of Pavia, Italy, 2002. [14] E.M. Darwish, H. Arenhövel, M. Schwamb, Eur. Phys. J. A 16 (2003) 111. References [15] M.I. Levchuk, M. Schumacher, F. Wissmann, nucl-th/0011041. [16] H. Arenhövel, Few-Body Syst. 4 (1988) 55. [1] A.Yu. Loginov, A.A. Sidorov, V.N. Stibunov, Phys. At. [17] R. Schmidt, H. Arenhövel, P. Wilhelm, Z. Phys. A 355 (1996) Nucl. 63 (2000) 391. 421. [2] T.-S.H. Lee, private communication; [18] M. Lacombe, et al., Phys. Lett. B 101 (1981) 139. T.-S.H. Lee, in: LOWq Workshop on Electromagnetic Nuclear [19] J. Haidenbauer, W. Plessas, Phys. Rev. C 30 (1984) 1822; Reactions at Low Momentum Transfer, 23–25 August 2001, J. Haidenbauer, W. Plessas, Phys. Rev. C 32 (1985) 1424. Halifax, Nova Scotia, Canada, 2001. [20] T. Sato, T.-S.H. Lee, Phys. Rev. C 54 (1996) 2660. [3] A.Yu. Loginov, A.V. Osipov, A.A. Sidorov, V.N. Stibunov, nucl-th/0407045. Physics Letters B 615 (2005) 68–78 www.elsevier.com/locate/physletb Top quark pair production and decay at a polarized photon collider A. Brandenburg a,1,Z.G.Sib,2 a DESY-Theorie, 22603 Hamburg, Germany b Department of Physics, Shandong University, Jinan, Shandong 250100, China Received 17 March 2005; received in revised form 2 April 2005; accepted 2 April 2005 Available online 12 April 2005 Editor: N. Glover Abstract Top quark pair production by (polarized) γγ collisions offers an interesting testing ground of the Standard Model and its extensions. In this Letter we present results for differential cross sections of top quark pair production and decay including QCD radiative corrections. We take into account the full dependence on the top quark spins. We give analytic and numerical results for single and double differential angular distributions of tt¯ decay products which are due to top quark polarizations and spin correlations in the intermediate state. 2005 Elsevier B.V. All rights reserved. Keywords: Photon collider; Top quarks; QCD corrections; Polarization; Spin correlations 1. Introduction At a future linear lepton collider, backscattered laser light may provide very high-energy photons [1], which would allow for a very interesting physics program [2,3]. In particular, top quark pair production is possible with large rates in (un)polarized photon–photon fusion. The measurement of the process γγ → ttX¯ is an important test of the Standard Model (SM). The first order QCD corrections to this process have already been calculated and found to be large [4–7]. The electroweak virtual plus soft-photonic O(α) corrections are also known [8].This process will also provide information on possible anomalous γtt¯ couplings [9,10] without contributions from Ztt¯ couplings present in e+e− collisions. Once the Higgs boson is discovered, it will be of primary importance to determine whether its properties are as predicted within the SM. In this respect, the process γγ → ttX¯ may play an important role. For example, E-mail address: [email protected] (Z.G. Si). 1 Work supported by a Heisenberg fellowship of DFG. 2 Work supported in part by NSFC and NCET. 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.04.003 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 69 heavy quark production in polarized γγ collisions will help to determine the parity of the Higgs boson produced as a resonance and decaying into top quark pairs [6,7]. In particular, if a Higgs boson is no CP eigenstate, spin correlations of the top quark pairs will help to probe the scalar and pseudoscalar couplings of the Higgs boson to the top quark [11,12]. For this kind of studies, predictions for top quark pair production and decay at a photon collider must be as precise as possible within the SM. In particular, the spin state of the intermediate tt¯ pair must be taken into account. (The role of the top quark polarization in probing the tt¯ threshold dynamics in γγ collisions was discussed in [13].) The purpose of this Letter is therefore to study the processes γγ → ttX¯ → + + X, + jets + X, all jets, (1.1) where stands for a charged lepton, with polarized photons from backscattered laser beams. We include QCD radiative corrections and take into account polarization and spin correlation effects of the intermediate tt¯ pairs. Leading order results and QCD corrections for the cross section and for top quark spin observables in the process ¯ 2 γγ → ttX are summarized in Sections 2 and 3. Numerical results to order α αs for the effective lepton collider cross section and for several decay distributions are given in Section 4. 2. Kinematics and leading order results 2 0 The production of top quark pairs by photon scattering at leading order α αs is described by the reaction ¯ γ(p1,λ1) + γ(p2,λ2) → t(k1,st ) + t(k2,st¯). (2.1) Here, p1,p2,k1 and k2 denote the momenta of the particles, λ1 and λ2 are the helicities of the photons, and the vectors st and st¯ describe the spins of top quark and antiquark. These fulfil the relations 2 = 2 =− · = · = st st¯ 1 and k1 st kt¯ st¯ 0. (2.2) In the (anti)top rest frame the spin of the (anti)top is described by a unit vector sˆt (sˆt¯). We choose the specific rest frames that are obtained by a rotation-free Lorentz boost from the zero momentum frame of the tt¯ quarks (tt¯-ZMF). Both the tt¯-ZMF and the t and t¯ rest frames will be used to construct spin observables from the final state momenta of the tt¯ decay products. We use the tt¯-ZMF rather than the c.m. frame of the colliding high-energy photons, since the latter system is probably more difficult to reconstruct experimentally. For the 2 → 2 process of Eq. (2.1), the two frames coincide. The differential cross section for the process of Eq. (2.1) can be written as follows: N 2 dσ(λ1,λ2,st ,st¯) = |M0| dΓ2, (2.3) 2sγγ 2 where the two-particle phase space measure is denoted by dΓ2, sγγ = (p1 + p2) and N = 3 is the number of colours. A simple calculation gives: 2 4 2 2 16α Qt π |M | = A + B p · (s + s¯) + B | ↔ p · (s + s¯) + C (s · s¯) + D (p · s )(p · s¯) 0 − 2 2 2 0 0 1 t t 0 λ1 λ2 2 t t 0 t t 0 1 t 2 t (1 β z ) + D0|z→−z(p2 · st )(p1 · st¯) , (2.4) with 2 2 4 2 2 2 2 4 2 2 A0 = 1 + 2β 1 − z − β 1 + 1 − z + λ1λ2 1 − 2β 1 − z − β z 2 − z , (2.5) 4m 2 2 2 2 2 B0 = λ1 1 − 2β + β z + λ2 1 − β z , (2.6) sγγ 70 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 2 4 2 2 2 4 2 2 C0 = 1 − 2β + β 1 + 1 − z + λ1λ2 1 − 2β + β z 2 − z , (2.7) 2 2 4(1 + βz)(1 − z )(1 − λ1λ2)β D0 =− . (2.8) sγγ Here, Qt = 2/3, m is the top quark mass, 4m2 β = 1 − , (2.9) s ¯ ˆ ˆ and z is the cosine of the scattering angle in the tt-ZMF, i.e., z = pˆ γ · k, where pˆ γ (k) is the direction of one of the photons (of the top quark) in that frame. 3. NLO results for γγ → ttX¯ In this section we present results for the inclusive reaction γγ → ttX¯ (3.1) 2 to order α αs . Apart from the cross section we study observables that depend on the spins of the top quark and antiquark. For polarized photons, observables of the form s O = 2St · aˆ (3.2) can have non-zero expectation values. Here, aˆ is an arbitrary reference direction and St is the top quark spin operator. The expectation value of Os is related to a single spin asymmetry: σ(↑) − σ(↓) Os = , (3.3) σ(↑) + σ(↓) where the arrows on the right-hand side refer to the spin state of the top quark with respect to the quantization axis aˆ. We will consider here two choices for aˆ, aˆ = kˆ (helicity basis), aˆ = pˆ (beam basis), (3.4) where kˆ denotes the direction of the top quark in the tt¯-ZMF and pˆ is the direction of the lepton beam coming from the left in that frame, which coincides to good approximation with the direction of one of the high-energy photons. Top quark polarization perpendicular to the plane spanned by pˆ and kˆ is induced by absorptive parts in the one-loop amplitude. This effect is, however, quite small (∼ a few percent) [14]. Analogous observables may of course be defined for the top antiquark. Apart from the above single spin observables, we also consider observables of the form d ˆ O = 4(St · aˆ)(St¯ · b). (3.5) ˆ Here, aˆ and b are arbitrary reference directions and St¯ is the top antiquark spin operator. The expectation value of Od is related to a double spin asymmetry: σ(↑↑) + σ(↓↓) − σ(↑↓) − σ(↓↑) Od = . (3.6) σ(↑↑) + σ(↓↓) + σ(↑↓) + σ(↓↑) For the reference directions we will consider here aˆ =−bˆ = kˆ, aˆ = bˆ = pˆ. (3.7) A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 71 Finally, we also present results for the observable ˜ d 4 O = S · S¯. (3.8) 3 t t The above double spin asymmetries Od and O˜ d have also proved useful for an analysis of spin correlations of top quark pairs in hadronic collisions [15]. The NLO cross section for the reaction γγ → ttX¯ may be written in terms of two scaling functions: α2Q4 σ(s,m,λˆ ,λ ) = t c(0)(ρ, λ ,λ ) + 4πα c(1)(ρ, λ ,λ ) . (3.9) 1 2 m2 1 2 s 1 2 Likewise, the unnormalized expectation values of the above spin observables are of the form α2Q4 σ O = t d(0)(ρ, λ ,λ ) + 4πα d(1)(ρ, λ ,λ ) , (3.10) a m2 a 1 2 s a 1 2 where a = 1 corresponds to the observable O˜ d defined in Eq. (3.8), a = 2(3) to the observable Od defined in Eq. (3.5) in the helicity (beam) basis, and a = 4(5) corresponds to the single spin observable Os defined in Eq. (3.2) in the helicity (beam) basis. The variable ρ is defined as 4m2 ρ = . (3.11) sγγ (0) (0) The lowest order scaling functions c √and da can be computed analytically. We use the following auxiliary functions, which vanish in the limit β = 1 − ρ → 0: 1 1 2 1 2 2 = ln(x) + 2β ,= ln(x) + 2β + β3 ,= ln(x) + 2β + β3 + β5 , 1 β 2 β3 3 3 β5 3 5 (3.12) where x = (1 − β)/(1 + β). We then obtain: ρ2 c(0)(ρ, λ ,λ ) = Nπβρ 1 + ρ − ρ2 + λ λ − 1 + ρ − − λ λ , 1 2 1 2 2 1 2 1 Nπβρ ρ2 d(0)(ρ, λ ,λ ) =− 1 + ρ + ρ2 + (1 + 2ρ)λ λ + 1 − − (1 + ρ)λ λ , 1 1 2 3 1 2 2 1 2 1 1 + 8ρ − 7ρ2 + ρ3 + (5 − 3ρ + ρ2)λ λ d(0)(ρ, λ ,λ ) = Nπβρ 1 2 2 1 2 3 ρ3 ρ2 − −1 + ρ − 2ρ2 + + 1 + λ λ , 2 2 1 2 2 √ 9 − 20 ρ − 6ρ + 14ρ3/2 + 29ρ2 + 6ρ5/2 − 20ρ3 + 3ρ4 d(0)(ρ, λ ,λ ) = Nπβρ − 3 1 2 15 √ −21 − 20 ρ + 25ρ + 14ρ3/2 − 16ρ2 + 6ρ5/2 − 3ρ3 + λ1λ2 15 √ ρ3 + −1 − 4ρ − ρ2 − 2 ρ + − 2ρ3/2 + ρ5/2 2 √ ρ2 √ + 1 + 2 ρ + 3ρ + λ λ (1 − ρ)2 , 2 1 2 3 72 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 Fig. 1. Left: scaling functions c(0)(ρ, 1, 1) (dotted), c(0)(ρ, 1, −1) (dash-dotted), c(1)(ρ, 1, 1) (full), and c(1)(ρ, 1, −1) (dashed). Right: scaling (0) (0) − (1) (1) − functions d1 (ρ, 1, 1) (dotted), d1 (ρ, 1, 1) (dash-dotted), d1 (ρ, 1, 1) (full), and d1 (ρ, 1, 1) (dashed). (0) (0) − (1) (1) − Fig. 2. Left: scaling functions d2 (ρ, 1, 1) (dotted), d2 (ρ, 1, 1) (dash-dotted), d2 (ρ, 1, 1) (full), and d2 (ρ, 1, 1) (dashed). Right: (0) (0) − (1) (1) − scaling functions d3 (ρ, 1, 1) (dotted), d3 (ρ, 1, 1) (dash-dotted), d3 (ρ, 1, 1) (full), and d3 (ρ, 1, 1) (dashed). ρ d(0)(ρ, λ ,λ ) = Nπβρ 1 − ρ − ln(x) (λ + λ ), 4 1 2 2 1 2 √ 1 − ρ √ 1 − ρ √ d(0)(ρ, λ ,λ ) =−Nπβρ 1 − 6 ρ + ρ3/2 + 2 + 2 ρ + 3ρ − ρ3/2 − ρ2 (λ − λ ). 5 1 2 3 2 2 1 2 (3.13) (1) (1) The functions c and da are obtained by a numerical integration. The scaling functions for the spin-averaged cross section and all spin observables are plotted in Figs. 1–3 for different choices of the photon helicities as a function of 1 η = − 1. (3.14) ρ The result for unpolarized photons can be inferred from 1 c(0),(1)(ρ, 0, 0) = c(0),(1)(ρ, 1, 1) + c(0),(1)(ρ, 1, −1) , (3.15) 2 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 73 (0) (0) − (1) (1) − Fig. 3. Left: scaling functions d4 (ρ, 1, 1) (dotted), d4 (ρ, 1, 1) (dash-dotted), d4 (ρ, 1, 1) (full), and d4 (ρ, 1, 1) (dashed). Right: (0) (0) − (1) (1) − scaling functions d5 (ρ, 1, 1) (dotted), d5 (ρ, 1, 1) (dash-dotted), d5 (ρ, 1, 1) (full), and d5 (ρ, 1, 1) (dashed). 1 d(0),(1)(ρ, 0, 0) = d(0),(1)(ρ, 1, 1) + d(0),(1)(ρ, 1, −1) , (3.16) 1,2,3 2 1,2,3 1,2,3 (0),(1) = d4,5 (ρ, 0, 0) 0. (3.17) As a check we compared our result for the functions c(0),(1)(ρ, 0, 0) with the results given in Fig. 2 of Ref. [4] and found perfect agreement. We further compared the functions c(0),(1)(ρ, 1, ±1) to the results given in Table 1 of Ref. [6]. After a trivial rescaling to account for the different conventions used in the definition of the scaling functions, we also found agreement. 4. Effective cross sections and spin observables 4.1. The effective cross section for γγ → ttX¯ The total tt¯ cross section at a photon collider may be written at NLO QCD as (cf., e.g., [5]) ymax ymax 2 4 α Qt e e (0) (1) σ ¯ = dy dy f (y ,P ,P )f (y ,P ,P ) c + 4πα c . (4.1) tt m2 1 2 γ 1 e L γ 2 e L s 0 0 e The function fγ (y1,Pe,PL) is the normalized energy spectrum of the photons resulting from Compton backscat- tering of laser light off the high energy electron beam. It is explicitly given by: − 1 f e(y, P ,P ) = N 1 − y + (2r − 1)2 − P P xr(2r − 1)(2 − y) . (4.2) γ e L 1 − y e L Here, Pe (PL) is the polarization of the electron (laser) beam, and y is the fraction of the electron energy in the c.m. frame transferred to the photon. It takes values in the range x 0 y ≡ y , (4.3) x + 1 max with 4E E x = L e , 2 (4.4) me 74 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 Table 1 √ Results for the effective cross section at see = 500 GeV (P ,P ; P ,P )σLO (fb) σ NLO (fb) K = σ NLO/σ LO e1 e2 L1 L2 tt¯ tt¯ tt¯ tt¯ (0, 0; 0, 0) 49.81 76.44 1.53 (0.85, 0.85;−1, −1) 175.86 260.77 1.48 (0.85, 0.85;+1, +1) 15.96 26.89 1.68 (0.85, −0.85;−1, +1) 48.99 71.93 1.47 where EL (Ee) is the energy of the laser (electron) beam and me is the electron mass. In order to avoid the creation of an e+e− pair from the backscattered laser beam and the low energy laser beam, the maximal value for x is √ xmax = 2 1 + 2 . (4.5) For a beam energy Ee = 250 GeV, this leads to an optimal laser energy EL ≈ 1.26 eV, (4.6) which will be used in the following numerical results. Finally, y r = . (4.7) x(1 − y) The normalization factor N in Eq. (4.2) is determined by ymax e = fγ (y, Pe,PL)dy 1. (4.8) 0 (0),(1) 2 The scaling functions c have to be evaluated at ρ = 4m /(y1y2see) and for polarizations = (i) (i) = λi Pγ yi,Pe ,PL ,i1, 2. (4.9) The function Pγ (y, Pe,PL) describes the degree of polarization of photons scattered with energy fraction y, which is given by 1 1 P (y, P ,P ) = xrP 1 + (1 − y)(2r − 1)2 − (2r − 1)P + 1 − y . (4.10) γ e L N e e L − fγ (y, Pe,PL) 1 y √ Numerical results for σtt¯ are given in Table 1 for see = 500 GeV and different polarizations of the laser and = = = = electron beam. We use the values m√t 178 GeV, α 1/128 and αs(µ mt ) 0.1. The QCD corrections to σtt¯ are quite large. This is because for see = 500 GeV most of the top quark pairs are produced close to threshold where the Coulombic β−1 singularity from soft gluons is important. 4.2. Spin observables The spin observables for γγ → ttX¯ discussed in Section 3 translate into observables built from the momenta of the tt¯ decay products. Thesinglespinasymmetries(3.2) cause a non-trivial one-particle inclusive decay distribution of the form 1 dσ(γγ → a1 + X) 1 = (1 + Bi cos θ1). (4.11) σ d cos θ1 2 Here, θ1 is the angle between the direction of a top quark decay product a1 measured in the top quark rest frame and one of the reference directions aˆ defined in Eq. (3.4). The coefficient Bi , with i = heli, beam for the helicity A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 75 and beam bases, is determined by the top quark spin asymmetry (3.2) and by the so-called spin analysing power of the decay product a1, which will be discussed below. The double spin asymmetries (3.5) lead to a two-particle inclusive decay distribution of the following form: → + 1 dσ(γγ a1a2 X) 1 ¯ = (1 + Bi cos θ1 + Bi cos θ2 − Ci cos θ1 cos θ2), (4.12) σ d cos θ1d cos θ2 4 where θ1 is defined as above and θ2 is analogously the angle between one of the top antiquark decay products and ˆ one of the reference directions b defined in Eq. (3.7). The coefficients Ci are determined by the double spin asym- metries (3.5) and the spin analysing powers of the two decay products a1 and a2. Finally, a non-zero expectation value of the observable defined in Eq. (3.8) leads to a distribution of the form 1 dσ(γγ → a a + X) 1 1 2 = (1 − D cos ϕ), (4.13) σ d cos ϕ 2 where ϕ is the angle between the direction of flight of the top decay product a1 and the antitop decay product a2 defined in the t and t¯ rest frames, respectively. We recall that these rest frames have to be obtained by a rotation-free boost from the tt¯-ZMF. The spin analysing power of the t and t¯ decay products is encoded in the one-particle inclusive angular distrib- utions dΓ/d cos θ for the decays ¯ t(st ) → a1(q1) + X1, t(st¯) → a2(q2) + X2. (4.14) Here q1 and q2 are the momenta of a1 and a2, respectively, defined in the rest frame of the (anti)top quark and θ is the angle between the polarization vector of the (anti)top quark and the direction of flight of a1(a2). For a fully polarized ensemble of top quarks (antiquarks) these distributions are of the form dΓ (1,2) Γ (1,2) = 1 ± κ(1,2) cos θ , (4.15) d cos θ 2 where Γ (1,2) is the partial width of the respective decay channel. The quantity κ(1,2) is the (anti)top-spin analysing power of a1,2. For the case of the standard (V − A) charged current interactions these distributions were computed to order αs for the semileptonic and non-leptonic channels in Refs. [16,17], respectively. (2) (1) As we work to lowest order in the electroweak couplings, Γ = Γ and κ2 = κ1 to all orders in αs ,ifthe channel a2 + X2 is the charge-conjugate of a1 + X1. + For semileptonic top decays t → b ν(g), the charged lepton is the most efficient analyser of the spin of the top quark. In the case of non-leptonic decays t → bqq¯(g) a good top spin analyser that can be identified easily is the least-energetic light quark jet. In Ref. [17] the coefficients κ(f ) were given to NLO accuracy for different choices of the spin analyser. To compute the coefficients Bi ,Ci and D we need the partial widths (sl,h) = (sl,h) + (sl,h) Γ a0 4παsa1 , (4.16) where the indices sl and h stand for semileptonic and hadronic decay modes. Further, we need the dimensionful coefficients (sl,h) (,j) = (,j) + (,j) Γ κ b0 4παsb1 , (4.17) where (j) refers to using the charged lepton (least-energetic light quark jet) as spin analyser. For the determination −5 −2 of these coefficients we use the Fermi constant GF = 1.16639 × 10 GeV , m = 178 GeV, mW = 80.42 GeV, 76 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 ΓW = 2.12 GeV, mb = 4.75 GeV, and all other quark and lepton masses are put to zero. (We do not use the narrow width approximation for the intermediate W boson.) We obtain, putting the CKM matrix elements |Vtb|= |Vqq |=1: ah ah = 0.52221 GeV,asl = 0 ,ah =−0.01968(15) GeV,asl =−0.01097(5) GeV. (4.18) 0 0 N 1 1 For the relevant coefficients b0,1 we obtain: = sl j = =− j =− b0 a0 ,b0 0.26950 GeV,b1 0.01118(8) GeV,b1 0.02375(26) GeV. (4.19) j The Durham algorithm was used as jet clustering scheme to obtain the four parton contribution to b1 . Within the leading pole approximation for the intermediate top quarks and antiquarks, the coefficients of the single and double differential distributions (4.11)–(4.13) are obtained in terms of the following quantities: ymax ymax α2Q4 1 σ = t dy dy f e(y ,P ,P )f e(y ,P ,P ) c(0)a(1) + πα c(1)a(1) + c(0)a(1) , s 2 1 2 γ 1 e L γ 2 e L 0 4 s 0 1 (4.20) m Γt 0 0 ymax ymax α2Q4 1 σ = t dy dy f e(y ,P ,P )f e(y ,P ,P ) d 2 2 1 2 γ 1 e L γ 2 e L m Γt 0 0 × (0) (1) (2) + (1) (1) (2) + (0) (1) (2) + (0) (1) (2) c a0 a0 4παs c a0 a0 c a1 a0 c a0 a1 , (4.21) ymax ymax α2Q4 1 N s = t dy dy f e(y ,P ,P )f e(y ,P ,P ) d(0)b(1) + πα d(1)b(1) + d(0)b(1) , r 2 1 2 γ 1 e L γ 2 e L r 0 4 s r 0 r 1 m Γt 0 0 (4.22) ymax ymax α2Q4 1 N d = t dy dy f e(y ,P ,P )f e(y ,P ,P ) r 2 2 1 2 γ 1 e L γ 2 e L m Γt 0 0 × (0) (1) (2) + (1) (1) (2) + (0) (1) (2) + (0) (1) (2) dr b0 b0 4παs dr b0 b0 dr b1 b0 dr b0 b1 . (4.23) We then get to NLO in αs : d d d s s N1 N2 N3 N4 N5 D = , Cheli = , Cbeam = , Bheli = , Bbeam = . (4.24) σd σd σd σs σs The LO and NLO results for these quantities are shown in Table 2 for favorable and realistic choices of electron and laser polarizations, using the same parameters as in Table 1. In most cases the QCD corrections are of the order of a few percent and thus much smaller than the corrections to the total tt¯ cross section. This was to be expected, since the bulk of the corrections is due to soft gluons which do not affect the tt¯ spin state. The biggest correction (∼ 11%) occurs for the coefficient Bbeam if the least energetic light jet is used as spin analyser. Photon polarization is an important asset: it is necessary to obtain polarized top quarks and thus non-zero coefficients Bheli and Bbeam. Further, the choice (Pe1,Pe2; PL1,PL2) = (0.85, 0.85;−1, −1), which increases the total yield of tt¯ pairs by more than a factor of 3 (see Table 1), in addition leads to larger tt¯ spin correlations. In particular, in the helicity basis the correlation is then almost 100% in the dilepton channel. A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 77 Table 2 √ Results for double and single spin asymmetries at see = 500 GeV (Pe1,Pe2; PL1,PL2) Dilepton Lepton–jet Jet–jet LO NLO LO NLO LO NLO D (0, 0; 0, 0) −0.670 −0.686 −0.346 −0.338 −0.178 −0.167 (0.85, 0.85;−1, −1) −0.806 −0.801 −0.416 −0.394 −0.215 −0.194 Cheli (0, 0; 0, 0) 0.811 0.826 0.418 0.408 0.216 0.201 (0.85, 0.85;−1, −1) 0.985 0.981 0.508 0.483 0.262 0.238 Cbeam (0, 0; 0, 0) −0.580 −0.606 −0.299 −0.299 −0.154 −0.148 (0.85, 0.85;−1, −1) −0.808 −0.804 −0.417 −0.396 −0.215 −0.195 (Pe1,Pe2; PL1,PL2) Lepton + XJet+ X LO NLO LO NLO Bheli (0.85, 0.85;−1, −1) 0.658 0.655 0.340 0.323 Bbeam (0.85, −0.85;−1, 1) −0.684 −0.637 −0.353 −0.314 So-called non-factorizable corrections do neither contribute at NLO QCD to σtt¯ nor to the angular correlations considered above. A proof of this statement is given in [15]. 5. Conclusions ¯ 2 We have computed a variety of spin observables for the process γγ → ttX up to order α αs . Together with the differential rates of polarized top and antitop quark decays at order αs , we have obtained the NLO QCD contributions to the fully differential cross section with intermediate top quark pair production at a photon collider.√ We have applied the above results to tt¯ production and decay at a future linear collider operating at s = 500 GeV. We have shown that for an appropriate choice of the polarizations of the laser and electron beam, the cross section and the double/single spin asymmetries can be quite large. While the QCD corrections to the cross section can be very large, most of the double/single spin asymmetries are affected at the level of only a few percent. The observables considered here will provide useful tools to analyse in detail the top quark pair production and decay dynamics. In particular, their precise measurement will test whether the top quark truly behaves as a quasi- free fermion as predicted in the Standard Model. Acknowledgements We wish to thank W. Bernreuther and P.M. Zerwas for discussions. A.B. was supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft. Z.G. Si wishes to thank DFG and MoE of China for financial support of his short visit in Germany, and also thanks DESY Theory Group for its hospitality. References [1] I.F. Ginzburg, G.L. Kotkin, S.L. Panfil, V.G. Serbo, V.I. Telnov, Nucl. Instrum. Methods A 219 (1984) 5. [2] S.J. Brodsky, P.M. Zerwas, Nucl. Instrum. Methods A 355 (1995) 19. [3] B. Badelek, et al., ECFA/DESY Photon Collider Working Group, Int. J. Mod. Phys. A 19 (2004) 5097, hep-ex/0108012. [4] M. Drees, M. Kramer, J. Zunft, P.M. Zerwas, Phys. Lett. B 306 (1993) 371. [5] J.H. Kuhn, E. Mirkes, J. Steegborn, Z. Phys. C 57 (1993) 615. [6] G. Jikia, A. Tkabladze, Phys. Rev. D 54 (1996) 2030, hep-ph/9601384. [7] G. Jikia, A. Tkabladze, Phys. Rev. D 63 (2001) 074502, hep-ph/0004068. 78 A. Brandenburg, Z.G. Si / Physics Letters B 615 (2005) 68–78 [8] A. Denner, S. Dittmaier, M. Strobel, Phys. Rev. D 53 (1996) 44, hep-ph/9507372. [9] S.Y. Choi, K. Hagiwara, Phys. Lett. B 359 (1995) 369, hep-ph/9506430. [10] B. Grzadkowski, Z. Hioki, K. Ohkuma, J. Wudka, Nucl. Phys. B 689 (2004) 108, hep-ph/0310159; B. Grzadkowski, Z. Hioki, K. Ohkuma, J. Wudka, Phys. Lett. B 593 (2004) 189, hep-ph/0403174. [11] H. Anlauf, W. Bernreuther, A. Brandenburg, Phys. Rev. D 52 (1995) 3803, hep-ph/9504424; H. Anlauf, W. Bernreuther, A. Brandenburg, Phys. Rev. D 53 (1996) 1725, Erratum. [12] W. Bernreuther, A. Brandenburg, M. Flesch, Phys. Rev. D 56 (1997) 90, hep-ph/9701347. [13] V.S. Fadin, V.A. Khoze, M.I. Kotsky, Z. Phys. C 64 (1994) 45, hep-ph/9403246. [14] W. Bernreuther, J.P. Ma, B.H.J. McKellar, Phys. Rev. D 51 (1995) 2475, hep-ph/9404235. [15] W. Bernreuther, A. Brandenburg, Z.G. Si, P. Uwer, Nucl. Phys. B 690 (2004) 81, hep-ph/0403035. [16] A. Czarnecki, M. Jezabek, J.H. Kühn, Nucl. Phys. B 351 (1991) 70. [17] A. Brandenburg, Z.G. Si, P. Uwer, Phys. Lett. B 539 (2002) 235, hep-ph/0205023. Physics Letters B 615 (2005) 79–86 www.elsevier.com/locate/physletb Decay constants of the pseudoscalar charmonium and bottomonium Zhi-Gang Wang a, Wei-Min Yang b, Shao-Long Wan b a Department of Physics, North China Electric Power University, Baoding 071003, PR China b Department of Modern Physics, University of Science and Technology of China, Hefei 230026, PR China Received 15 November 2004; received in revised form 22 February 2005; accepted 23 March 2005 Available online 14 April 2005 Editor: W. Haxton Abstract In this Letter, we investigate the structures of the pseudoscalar charmonium and bottomonium in the framework of the coupled rainbow Schwinger–Dyson equation and ladder Bethe–Salpeter equation with the confining effective potential (infrared modified flat bottom potential). As the current masses are very large, the dressing or renormalization for the c and b quarks are tender, however, mass poles in the timelike region are absent. The Euclidean time Fourier transformed quark propagator has no mass poles in the timelike region which naturally implements confinement. The Bethe–Salpeter wavefunctions for those mesons have the same type (Gaussian type) momentum dependence and center around zero momentum with spatial extension to about q2 = 1GeV2 which happen to be the energy scale for chiral symmetry breaking, the strong interactions in the infrared region result in bound states. The decay constants for those pseudoscalar heavy quarkonia are compatible with the values of experimental extractions and theoretical calculations. 2005 Elsevier B.V. All rights reserved. PACS: 14.40.-n; 11.10.Gh; 11.10.St; 12.40.Qq Keywords: Schwinger–Dyson equation; Bethe–Salpeter equation; Decay constant; Confinement 1. Introduction the heavy quarks), the soft scale (the relative momen- tum of the heavy quark–antiquark |p|) and the ultrasoft Heavy quarkonium, bound state of the heavy quark scale (the typical kinetic energy of the heavy quark– and antiquark, characterized by at least three widely antiquark E), plays a special role in probing the strong separated energy scales: the hard scale (the mass m of interactions in both the perturbative and nonperturba- tive regions. By definition of the heavy quark, m is large in comparison with the typical hadronic scale E-mail address: [email protected] (Z.-G. Wang). ΛQCD, the corresponding processes can be success- 0370-2693/$ – see front matter 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.physletb.2005.03.080 80 Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 fully described in perturbative quantum chromody- well as the meson structures, such as electromagnetic namics (QCD) due to the asymptotic freedom. How- form factors, radius, decay constants [18,22,23]. ever, the lower scales |p| and E, which are responsible During the past two years, the experiments have for the binding, cannot be accessible by perturbation discovered a number of new states, for example, the ηc − + theory. The appearance of multiscales in the dynamics in exclusive B → KKSK π decays by Belle [24], of the heavy quarkonium makes its quantitative study the narrow DsJ states by BaBar, CLEO and Belle [25], extremely difficult, the properties of the bound states evidence for the Θ+(1540) with quantum numbers and their decays can provide powerful test for QCD in of K+n [26], and the X(3872) through decay to both the perturbative and nonperturbative regions. π +π −J/ψ by Belle [27]. New experimental results The physicists propose many original approaches call for interpretations, offer opportunities to extend to deal with the long distance properties of QCD, such our knowledge about hadron spectrum and challenge as chiral perturbation theory [1], heavy quark effective our understanding of the strong interaction; further- theory [2], QCD sum rules [3], lattice QCD [4], pertur- more, they revitalize the study of heavy quarkonia and bative QCD [5], coupled Schwinger–Dyson equation stimulate a lot of theoretical analysis through the char- (SDE) and Bethe–Salpeter equation (BSE) method monia and bottomonia have been thoroughly investi- [6], nonrelativistic QCD [7], potential nonrelativistic gated. QCD [8], etc. All of those approaches have both out- The decay constants of the pseudoscalar charmo- standing advantages and obvious shortcomings in one nium and bottomonium (ηc and ηb) mesons play an or other ways. The coupled rainbow SDE and ladder important role in modern physics with the assumption BSE have given a lot of successful descriptions of of current-meson duality. The precise knowledge of the long distance properties of the low energy QCD the those values fηc and fηb will provide great im- and the QCD vacuum (for example, Refs. [9–12],for provements in our understanding of various processes recent reviews one can see Refs. [13,14]). The SDE convolving the ηc and ηb mesons, for example, the can naturally embody the dynamical symmetry break- process B → ηcK, where the mismatches between ing and confinement which are two crucial features of the theoretical and experimental values are large [28]. QCD, although they correspond to two very different The ηc meson is already observed experimentally, the energy scales [15,16]. On the other hand, the BSE is current experimental situation with the ηb meson is a conventional approach in dealing with the two body rather uncertain, yet the discovery of the ηb meson relativistic bound state problems [17].Fromthesolu- is one of the primary goals of the CLEO-c research tions of the BSE, we can obtain useful information program [29]; furthermore, the ηb meson may be ob- about the under-structure of the mesons and obtain served in run II at the Tevatron through the decay powerful tests for the quark theory. However, the obvi- modes into charmed states D∗D(∗) [30]. It is inter- ously drawback may be the model dependent kernels esting to combine those successful potential mod- for the gluon two-point Green’s function and the trun- els within the framework of coupled SDE and BSE cations for the coupled divergent SDE and BSE series to calculate the decay constants of the pseudoscalar in one or the other ways [18]. Many analytical and nu- heavy quarkonia such as ηc and ηb. For previous stud- merical calculations indicate that the coupled rainbow ies about the electroweak decays of the pseudoscalar SDE and ladder BSE with phenomenological potential mesons with the SDE and BSE, one can consult models can give model independent results and satis- Refs. [6,9–14]. In this Letter, we use an infrared mod- factory values [6,9–14]. The usually used effective po- ified flat-bottom potential (IMFBP) which takes the tential models are confining Dirac δ function potential, advantages of both the Gaussian distribution potential Gaussian distribution potential and flat bottom poten- and the FBP to calculate the decay constants of those tial (FBP) [13,14,19–21]. The FBP is a sum of Yukawa pseudoscalar heavy quarkonia. potentials, which not only satisfies chiral invariance The Letter is arranged as follows: we introduce the and fully relativistic covariance, but also suppresses IMFBP in Section 2; in Section 3–5,wesolvethe the singular point that the Yukawa potential has. It rainbow SDE and ladder BSE, explore the analytic- works well in understanding the dynamical chiral sym- ity of the quark propagators, investigate the dynam- metry breaking, confinement and the QCD vacuum as ical dressing and confinement, finally obtain the de- Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 81 dV(0) d2V(0) dnV(0) cay constants for those pseudoscalar heavy quarkonia; = =···= = 0. (4) Section 6 is reserved for conclusion. dr dr2 drn The aj can be determined by solve the equations in- ferred from the flat bottom condition Eq. (4).Asin 2. Infrared modified flat bottom potential previous literature [18,21–23], n issettobe9.The phenomenological effective potential (IMFBP) can be The present techniques in QCD calculation cannot taken as give satisfactory large r behavior for the gluon two- 2 = 2 + 2 point Green’s function to implement the linear poten- G k G1 k G2 k . (5) tial confinement mechanism, in practical calculation, the phenomenological effective potential models al- ways do the work. As in our previous work [18],we 3. SchwingerÐDyson equation use a Gaussian distribution function to represent the infrared behavior of the gluon two-point Green’s func- The SDE can provide a natural framework for in- tion, vestigating the nonperturbative properties of the quark and gluon Green’s functions. By studying the evo- 2 2 2 2 − k lution behavior and analytic structure of the dressed 4πG1 k = 3π e ∆ , (1) ∆2 quark propagators, we can obtain valuable information which determines the quark–antiquark interaction about the dynamical dressing phenomenon and con- through a strength parameter and a range parame- finement. In the following, we write down the rainbow ter ∆. This form is inspired by the δ function poten- SDE for the quark propagator, tial (in other words the infrared dominated potential) − S 1(p) = iγ · p +ˆm used in Refs. [19,20], which it approaches in the limit c,b ∆ → 0. For the intermediate momentum, we take the d4k λa λa + 4π γµ S(k)γν Gµν(k − p), FBP as the best approximation and neglect the large (2π)4 2 2 momentum contributions from the perturbative QCD (6) calculations as the coupling constant at high energy where is very small. The FBP is a sum of Yukawa poten- − S 1(p) = iA p2 γ · p + B p2 tials which is an analogy to the exchange of a series of particles and ghosts with different masses (Euclidean ≡ A p2 iγ · p + m p2 , (7) form), k k G (k) = δ − µ ν G k2 , (8) n µν µν 2 2 aj k G2 k = , (2) k2 + (N + jρ)2 ˆ j=0 and mc,b stands for the current quark mass that explic- itly breaks chiral symmetry. where N stands for the minimum value of the masses, The full SDE for the quark propagator is a divergent ρ is their mass difference, and aj is their relative cou- series of coupled nonlinear integral equations for the pling constant. Due to the particular condition we take propagators and vertexes, we have to make truncations for the FBP, there is no divergence in solving the SDE. in one or other ways. The rainbow SDE has given a lot In its three-dimensional form, the FBP takes the fol- of successful descriptions of the QCD vacuum and low lowing form: energy hadron phenomena [6,13–16], in this Letter, we n e−(N+jρ)r take the rainbow SDE. If we go beyond the rainbow V(r)=− a . (3) λa j approximation, the bare vertex γµ 2 has to be substi- = r a j 0 tuted by the full quark–gluon vertex Γµ (qqg), which In order to suppress the singular point at r = 0, we satisfies the Slavnov–Tayler identity. In the weak cou- 2 → take the following conditions: pling limit, g 0, two Feynman diagrams contribute a to the vertex Γµ (qqg) at one-loop level due to the V(0) = const, non-Abelian nature of QCD, i.e., the self-interaction 82 Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 of gluons [31]. If we neglect the contributions from the BSE for the pseudoscalar quarkonia, a three-gluon vertex Γµ (ggg) and retain an Abelian ver- a −1 P −1 P sion, the vertex Γ a(qqg) can be taken as λ Γ (qqp), S+ q + χ(q,P)S− q − µ 2 µ 2 2 where the vertex Γµ(qqp) is the quark–photon vertex 16π d4k which satisfies the Ward–Takahashi identity. In practi- = γ χ(k,P)γ G (q − k), 4 µ ν µν (9) cal calculation, we can take the vertex Γµ(qqp) to be 3 (2π) the Ball–Chiu and Curtis–Pennington vertex [32,33] where S(q) is the quark propagator, Gµν(k) is the so as to avoid solving the coupled SDE for the vertex gluon propagator, Pµ is the four-momentum of the Γµ(qqp). However, the nonperturbative properties of center of mass of the pseudoscalar quarkonia, qµ is QCD at the low energy region suggest that the SDEs the relative four-momentum between the quark and are strongly coupled nonlinear integral equations, no antiquark, γµ is the bare quark–gluon vertex, and theoretical work has ever proven that the contributions χ(q,P)is the Bethe–Salpeter wavefunction (BSW) of a from the vertex Γµ (ggg) can be safely neglected due the bound state. to the complex Dirac and tensor structures. The one We can perform the Wick rotation analytically and Feynman diagram contributions version of the ver- continue q and k into the Euclidean region.1 In the a tex Γµ (qqg), i.e., neglecting the contributions from lowest order approximation, the BSW χ(q,P) can be a a the vertex Γµ (ggg) in dressing the vertex Γµ (qqg) written as is inconsistent with the Slavnov–Tayler identity [31]. = 0 + · 0 If we take the assumption that the contributions from χ(q,P) γ5 iF1 (q, P ) γ PF2 (q, P ) the vertex Γ a(ggg) are not different greatly from the + · · 1 µ γ qq PF3 (q, P ) vertex Γ a(qqg), we can multiply the contributions µ + [ · · ] 0 a i γ q,γ P F4 (q, P ) . (10) from the vertex Γµ (qqg) by some parameters which effectively embody the contributions from the vertex The ladder BSE can be projected into the following a Γµ (ggg) [34]. four coupled integral equations, In this Letter, we assume that a Wick rotation to ∞ π Euclidean variables is allowed, and perform a rotation 0,1 = 3 2 analytically continuing p and k into the Euclidean re- H(i,j)Fj (q, P ) k dk sin θK(i, j), gion. The Euclidean rainbow SDE can be projected j j 0 0 (11) into two coupled integral equations for A(p2) and B(p2). Alternatively, one can derive the SDE from the expressions of the H(i,j)and K(i,j) are cumber- the Euclidean path-integral formulation of the theory, some and neglected here. thus avoiding possible difficulties in performing the Here we will give some explanations for the expres- Wick rotation [35]. As far as only numerical results sions of H(i,j).TheH(i,j)’s are functions of the 2 + are concerned, the two procedures are equal. In fact, quark’s Schwinger–Dyson functions (SDF) A(q 2 + · 2 + 2 + · 2 + 2 − the analytical structures of quark propagators have in- P /4 q P), B(q P /4 q P), A(q P /4 · 2 + 2 − · teresting information about confinement, we will make q P) and B(q P /4 q P). The relative four- detailed discussion about the c and b quarks propaga- momentum q is a quantity in the Euclidean space- tors respectively in Section 5. time while the center of mass four-momentum P must be continued to the Minkowski spacetime, i.e., P 2 = −m2 q · P ηc,ηb , this results in the varying throughout a complex domain. It is inconvenient to solve the SDE 4. BetheÐSalpeter equation at the resulting complex values of the quark momen- tum, especially for the heavy quarks. As the dressing The BSE is a conventional approach in dealing with effect is minor, we can expand A and B in terms of the two body relativistic bound state problems [17]. The precise knowledge about the quark structures of 1 To avoid possible difficulties in performing the Wick rotation, the mesons will result in better understanding of their one can derive the BSE from the Euclidean path-integral formula- properties. In the following, we write down the ladder tion of the theory. Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 83 Taylor series of q · P , for example, 5. Coupled rainbow SDE and ladder BSE and the decay constants A q2 + P 2/4 + q · P In this section, we explore the coupled equations of = A q2 + P 2/4 + A q2 + P 2/4 q · P +···. the rainbow SDE and ladder BSE for the pseudoscalar heavy quarkonia numerically, the final results for the The other problem is that we cannot solve the SDE SDFs and BSWs can be plotted as functions of the in the timelike region as the two point gluon Green’s square momentum q2. function cannot be exactly inferred from the SU(3) In order to demonstrate the confinement of quarks, color gauge theory even in the low energy spacelike re- we have to study the analyticity of SDFs for the c and gion. In practical calculations, we can extrapolate the b quarks, and prove that there no mass poles on the real values of A and B from the spacelike region smoothly timelike q2 axial. In the following, we take the Fourier to the timelike region with suitable polynomial func- transform with respect to the Euclidean time T for the tions. To avoid possible violation with confinement in scalar part (Ss ) of the quark propagator [6,13,37], sense of the appearance of pole masses q2 =−m2(q2) +∞ in the timelike region, we must be care in choos- 2 ∗ dq4 B(q ) S (T ) = eiq4T , ing the polynomial functions [20].Fortheηc meson, s 2 2 2 2 2 2π q A (q ) + B (q ) = the mass is about 3.0 GeV, the extrapolation to the −∞ q 0 − 2 (14) timelike region with the quantity mηc /4 can be per- formed easily, however, the large mass of the ηb meson where the 3-vector part of q is set to zero. If S(q) has a makes the extrapolation into the deep timelike region mass pole at q2 =−m2(q2) in the real timelike region, troublesome. Although the η meson has not been ∗ −mT b the Fourier transformed Ss (T ) would fall off as e observed experimentally yet, the theoretical calcula- ∗ =− for large T or log Ss mT . tions indicate that its mass is about 9.4GeV[36].As In our numerical calculations, for small T ,theval- the dressed quark propagators comprise the notation ues of S∗ are positive and decrease rapidly to zero and 2 s of constituent quarks by providing a mass m(q ) = beyond with the increase of T , which are compatible 2 2 B(q )/A(q ), which corresponding to the dynamical with the result (curve tendency with respect to T ) from symmetry breaking phenomena for the light quarks. ∗ lattice simulations [38];forlargeT , the values of Ss We can simplify the calculation greatly and avoid the are negative, except occasionally a very small fraction problems concerning the extrapolations in solving the ∗ of positive values. The negative values for Ss indicate BSE by take the following propagator for the c and b an explicit violation of the axiom of reflection posi- quarks, tivity [39], in other words, the quarks are not physical observable, i.e., confinement. −1 2 S q = iγ · q + Mc,b, (12) For the c and b quarks, the current masses are very large, the dressing or renormalization is tender and the where the Mc,b is the Euclidean constituent quark curves are not steep which in contrast to the dynam- 2 = 2 2 = 2 mass with Mc,b mc,b(q ) q obtained from the ical chiral symmetry breaking phenomenon for the solution of the SDE (6). light quarks, mc(0)/mˆ c 1.5 and mb(0)/mˆ b 1.1, Finally we write down the normalization condition however, mass poles in the timelike region are ab- for the BSW, sent. At zero momentum, mc(0) = 1937 MeV and m (0) = 5105 MeV, while the Euclidean constituent b 4 −1 quark masses M = 1908 MeV and M = 5096 MeV, d q ∂S+ −1 c b N χ¯ χ(q,P)S 2 2 2 2 c 4 Tr − which defined by M = m (q ) = q , are compati- (2π) ∂Pµ ble with the constituent quark masses in the litera- −1 −1 ∂S− ture. From the plotted BSWs (see Fig. 1 as an ex- +¯χS+ χ(q,P) = 2Pµ, (13) ∂Pµ ample), we can see that the BSWs for pseudoscalar mesons have the same type (Gaussian type) momen- + where χ¯ = γ4χ γ4. tum dependence while the quantitative values are dif- 84 Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 b quarks are taken in solving the BSE as the dress- ing is tender. We borrow some idea from the fact that the simple phenomenological model of Cornell po- tential (Coulomb potential plus linear potential) with constituent quark masses can give satisfactory mass spectrum for the heavy quarkonia3 and take larger values for the strength parameter and range para- meter ∆, i.e., = 2.2 GeV and ∆ = 2.9GeV2,in the infrared region comparing with the correspond- ing ones used in Ref. [18]. Furthermore the masses of the pseudoscalar mesons are taken as input parame- ters. If we take the Euclidean constituent quark masses Mc = mc(0) and Mb = mb(0), the decay constants for = the ηc and ηb mesons change slightly, fηc 357 MeV Fig. 1. BSWs for charmonium. = and fηb 289 MeV. ferent from each other. Just like the lighter qq¯ and qQ¯ pseudoscalar mesons [18], the Gaussian type BSWs 6. Conclusion center around zero momentum with spatial extension to about q2 = 1GeV2 which happen to be the energy In this Letter, we investigate the under-structures scale for chiral symmetry breaking, the strong inter- of the pseudoscalar heavy quarkonia η and η in the actions in the infrared region result in bound states. c b framework of the coupled rainbow SDE and ladder Finally we obtain the values for the decay constants of BSE with the confining effective potential (IMFBP). those pseudoscalar mesons which are defined by After we solve the coupled rainbow SDE and ladder ifπ Pµ =0|¯qγµγ5q π(P) , BSE numerically, we obtain the SDFs and BSWs for the pseudoscalar heavy quarkonia η and η .Asthe d4k c b = N γ γ χ(k,P) , current masses of the c and b quarks are very large, c Tr µ 5 4 (15) (2π) the dressing or renormalization for the SDFs is ten- hereweuseπ to represent the pseudoscalar mesons,2 der and the curves are not steep which in contrast to the explicitly dynamical chiral symmetry break- = = fηc 349 MeV,fηb 287 MeV, (16) ing phenomenon for the light quarks, however, mass which are compatible with the results from the ex- poles in the timelike region are absent. We can sim- perimental extractions and theoretical calculations, plify the calculation greatly and avoid the problems = ± ≈ concerning the extrapolations in solving the BSE by fηc 335 75 MeV (exp) [40]; fηc 400 MeV = ± = ± making the substitution B(q2) → M and A(q2) → 1. (exp) [41]; fηc 420 52 MeV, fηb 705 27 MeV = ± The BSWs for the pseudoscalar heavy quarkonia have (theor) [42]; fηc 292 25 MeV (theor) [43]; ≈ = ± the same type (Gaussian type) momentum depen- fηc 350 MeV (theor) [44]; fηc 300 50 MeV (theor) [45]. In calculation, the values of mˆ c and mˆ b dence while the quantitative values are different from are taken as the current quark masses, mˆ c = 1250 MeV each other. The Gaussian type BSWs center around 2 = and mˆ b = 4700 MeV; the input parameters for the zero momentum with spatial extension to about q 2 FBP are N = 1.0Λ, V(0) =−11.0Λ, ρ = 5.0Λ and 1GeV which happen to be the energy scale for chi- Λ = 200 MeV, which are determined in study of the ral symmetry breaking, the strong interactions in the qq¯ and qQ¯ pseudoscalar mesons [18]. In this Letter, infrared region result in bound states. Our numerical the Euclidean constituent quark masses for the c and results for the values of the decay constants of the 2 3 Here we write down the Nc explicitly according to the normal- For an excellent review of the potential models, one can consult ization condition (13). Ref. [46]. Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 85 pseudoscalar heavy quarkonia are compatible with the M.A. Ivanov, Yu.L. Kalinovsky, P. Maris, C.D. Roberts, Phys. corresponding ones obtained from the experimental Lett. B 416 (1998) 29; extractions and theoretical calculations. Once the sat- M.A. Ivanov, Yu.L. Kalinovsky, P. Maris, C.D. Roberts, Phys. Rev. C 57 (1998) 1991. isfactory SDFs and BSWs for the pseudoscalar heavy [13] C.D. Roberts, S.M. Schmidt, Prog. Part. Nucl. Phys. 45 (2000) quarkonia are known, we can use them to investigate a S1. lot of important quantities involving the B, ηc and ηb [14] C.D. Roberts, nucl-th/0304050; mesons. P. Maris, C.D. Roberts, Int. J. Mod. Phys. E 12 (2003) 297. [15] V.A. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories, Word Scientific, Singapore, 1993. [16] R. Alkofer, L.V. Smekal, hep-ph/0007355; Acknowledgements C.S. Fischer, R. Alkofer, hep-ph/0301094. [17] E.E. Salpeter, H.A. Bethe, Phys. Rev. 84 (1951) 1232; N. Nakanishi, Suppl. Prog. Theor. Phys. 43 (1969) 1. This work is supported by National Natural Science [18] Z.G. Wang, W.M. Yang, S.L. Wan, Phys. Lett. B 584 (2004) Foundation, Grant No. 10405009, and Key Program 71; Foundation of NCEPU. The authors are indebted to Z.G. Wang, W.M. Yang, S.L. Wan, Nucl. Phys. A 744 (2004) Dr. J. He (IHEP) for numerous help, without him, 156. [19] H.J. Munczek, A.M. Nemirovsky, Phys. Rev. D 28 (1983) 181. the work would not be finished. One of the authors [20] H.J. Munczek, P. Jain, Phys. Rev. D 46 (1991) 438; (Z.-G. Wang) would like to thank Dr. Gogohia for P. Jain, H.J. Munczek, Phys. Rev. D 44 (1991) 1873; helpful discussion. P. Jain, H.J. Munczek, Phys. Rev. D 48 (1993) 5403. [21] K.L. Wang, S.L. Wan, Phys. Rev. D 47 (1993) 2098. [22] Z.G. Wang, S.L. Wan, K.L. Wang, Phys. Lett. B 498 (2001) 195; References Z.G. Wang, S.L. Wan, Phys. Lett. B 536 (2001) 241; Z.G. Wang, J. Phys. G 28 (2002) 3007; [1] J. Gasser, H. Leutwyler, Ann. Phys. 158 (1984) 142; Z.G. Wang, S.L. Wan, K.L. Wang, Commun. Theor. Phys. 35 (2001) 697; J. Gasser, H. Leutwyler, Nucl. Phys. B 250 (1985) 465. [2] M. Neubert, Phys. Rep. 245 (1994) 259. Z.G. Wang, S.L. Wan, K.L. Wang, Chin. Phys. 10 (2001) 497; Z.G. Wang, S.L. Wan, K.L. Wang, Chin. Phys. Lett. 18 (2001) [3] M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. B 147 (1979) 385; 743. [23] S.L. Wan, K.L. Wang, J. Phys. G 22 (1996) 1287. M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. B 147 (1979) 448. [24] Belle Collaboration, S.-K. Choi, et al., Phys. Rev. Lett. 89 [4] R. Gupta, hep-lat/9807028, for an elementary introduction. (2002) 102001; [5] S.J. Brodsky, G.P. Lepage, Phys. Lett. B 87 (1979) 359; Belle Collaboration, S.-K. Choi, et al., Phys. Rev. Lett. 89 S.J. Brodsky, G.P. Lepage, Phys. Rev. D 22 (1980) 2157. (2002) 129901, Erratum. [6] C.D. Roberts, A.G. Williams, Prog. Part. Nucl. Phys. 33 (1994) [25] BaBar Collaboration, B. Aubert, et al., Phys. Rev. Lett. 90 477. (2003) 242001; CLEO Collaboration, D. Besson, et al., Phys. Rev. D 68 (2003) [7] G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. D 51 (1995) 1125; 032002; Belle Collaboration, P. Krokovny, et al., Phys. Rev. Lett. 91 G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. D 55 (1997) 5853, Erratum. (2003) 262002. [26] LEPS Collaboration, Phys. Rev. Lett. 91 (2003) 012002; [8] A. Pineda, J. Soto, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 428; N. Brambilla, A. Pineda, J. Soto, A. Vairo, Nucl. Phys. B 566 V.V. Barmin, et al., Phys. At. Nucl. 66 (2003) 1715; CLAS Collaboration, Phys. Rev. Lett. 91 (2003) 252001; (2000) 275; N. Brambilla, A. Pineda, J. Soto, A. Vairo, hep-ph/0410047. CLAS Collaboration, Phys. Rev. Lett. 92 (2004) 032001; CLAS Collaboration, Phys. Rev. Lett. 92 (2004) 049902, Erra- [9] Y.B. Dai, C.S. Huang, D.S. Liu, Phys. Rev. D 43 (1991) 1717; Y.B. Dai, Y.B. Ding, C.S. Huang, C.L. Wang, Commun. Theor. tum; SAPHIR Collaboration, Phys. Lett. B 572 (2003) 127; Phys. 18 (1992) 313. [10] P. Maris, C.D. Roberts, Phys. Rev. C 56 (1997) 3369; A.E. Asratyan, et al., hep-ex/0309042; HERMES Collaboration, Phys. Lett. B 585 (2004) 213; P. Maris, C.D. Roberts, P.C. Tandy, Phys. Lett. B 420 (1998) 267. SVD Collaboration, hep-ex/0401024; M. Abdel-Bary, et al., hep-ex/0403011; [11] P. Maris, P.C. Tandy, Phys. Rev. C 60 (1999) 055214; P. Maris, Nucl. Phys. A 663 (2000) 621. ZEUS Collaboration, hep-ex/0403051. [27] Belle Collaboration, S.-K. Choi, et al., Phys. Rev. Lett. 91 [12] M.A. Ivanov, Yu.L. Kalinovsky, C.D. Roberts, Phys. Rev. D 60 (1999) 034018; (2003) 262001. 86 Z.-G. Wang et al. / Physics Letters B 615 (2005) 79–86 [28] L. Li, Z.G. Wang, T. Huang, Phys. Rev. D 70 (2004) 074006. ALEPH Collaboration, A. Heister, et al., Phys. Lett. B 530 [29] H. Stock, hep-ex/0310021. (2002) 56, a review of the mass of the ηb meson from theo- [30] F. Maltoni, A.D. Polosa, hep-ph/0405082. retical calculations. [31] A.I. Davydychev, P. Osland, L. Saks, Phys. Rev. D 63 (2001) [37] P. Maris, Phys. Rev. D 52 (1995) 6087. 014022; [38] M.S. Bhagwat, M.A. Pichowsky, C.D. Roberts, P.C. Tandy, A.I. Davydychev, P. Osland, L. Saks, Nucl. Phys. B (Proc. Phys. Rev. C 68 (2003) 015203. Suppl.) 89 (2000) 277; [39] J. Glimm, A. Jaffee, Quantum Physics. A Functional Point of V. Gogohia, hep-th/0406064. View, Springer-Verlag, New York, 1981. [32] J.S. Ball, T.W. Chiu, Phys. Rev. D 22 (1980) 2542. [40] CLEO Collaboration, K.W. Edwards, et al., Phys. Rev. Lett. 86 [33] D.C. Curtis, M.R. Pennington, Phys. Rev. D 42 (1990) 4165; (2001) 30. D.C. Curtis, M.R. Pennington, Phys. Rev. D 44 (1991) 536; [41] Particle Data Group, Phys. Rev. D 54 (1996) 1. D.C. Curtis, M.R. Pennington, Phys. Rev. D 46 (1992) 2663. [42] D.S. Hwang, G.H. Kim, Z. Phys. C 76 (1997) 107. [34] M.S. Bhagwat, A. Holl, A. Krassnigg, C.D. Roberts, P.C. [43] G. Cvetic, C.S. Kim, G.L. Wang, W. Namgung, Phys. Lett. Tandy, Phys. Rev. C 70 (2004) 035205; B 596 (2004) 84. M.S. Bhagwat, P.C. Tandy, Phys. Rev. D 70 (2004) 094039. [44] V.A. Novikov, L.B. Okun, M.A. Shifman, A.I. Vainshtein, V.I. [35] S.T. Stainsby, R.T. Cahill, Phys. Lett. A 146 (1990) 467. Zakharov, Phys. Rep. 41 (1978) 1. [36] B.A. Kniehl, A.A. Penin, A. Pineda, V.A. Smirnov, M. Stein- [45] N.G. Deshpande, J. Trampetic, Phys. Lett. B 339 (1994) 270. hauser, Phys. Rev. Lett. 92 (2004) 242001; [46] C. Quigg, J.L. Rosner, Phys. Rep. 56 (1979) 167. Physics Letters B 615 (2005) 87–92 www.elsevier.com/locate/physletb Enlarged Galilean symmetry of anyons and the Hall effect P.A. Horváthy a, L. Martina b,P.C.Stichelc a Laboratoire de Mathématiques et de Physique Théorique, Université de Tours, Parc de Grandmont, F-37 200 Tours, France b Dipartimento di Fisica dell’Università, and Sezione INFN di Lecce, Via Arnesano, CP 193, I-73 100 Lecce, Italy c An der Krebskuhle 21, D-33 619 Bielefeld, Germany Received 23 March 2005; accepted 2 April 2005 Available online 12 April 2005 Editor: N. Glover Abstract Enlarged planar Galilean symmetry, built of both space–time and field variables and also incorporating the “exotic” central extension is introduced. It is used to describe non-relativistic anyons coupled to an electromagnetic field. Our theory exhibits an anomalous velocity relation of the type used to explain the anomalous Hall effect. The Hall motions, characterized by a Casimir of the enlarged algebra, become mandatory for some critical value(s) of the magnetic field. The extension of our scheme yields the semiclassical effective model of the Bloch electron. 2005 Elsevier B.V. All rights reserved. 1. Introduction field yields the first-order phase space Lagrangian