V A C U U M S T R U C T U R E O F Q C D IN A N E F F E C T I V E L A G R A N G I A N A P P R O A C H By Todd Darwin Fugleberg B. Sc., The University of Saskatchewan, 1993 M. Sc., The University of British Columbia, 1996 A THESIS THE SUBMITTED IN PARTIAL REQUIREMENTS DOCTOR FULFILLMENT OF FOR T H E DEGREE OF OF PHILOSOPHY in THE FACULTY OF GRADUATE DEPARTMENT OF STUDIES PHYSICS We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA August 2000 © Todd Darwin Fugleberg, 2000 In presenting this thesis i n partial fulfilment of the requirements for an advanced degree at the University of B r i t i s h Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Physics T h e University of B r i t i s h C o l u m b i a 6224 A g r i c u l t u r a l R o a d Vancouver, B . C . , C a n a d a V 6 T 1Z1 Date: Abstract The vacuum structure of QCD is studied using an anomalous effective Lagrangian approach. This approach makes it possible to determine how physical observables depend on the strong CP violation parameter, 6. The ^-dependence of QCD and the phenomenology of the light pseudoscalar mesons in this theory are illustrated. The vacuum structure of QCD is shown to be quite complex with the prediction of a number of different types of nontrivial vacuum states. Two specific examples of nontrivial vacuum states are analysed in more detail. The decay rate of a metastable vacuum state for the 6 = 0 case is nonperturbatively calculated in both the zero and high temperature limits. The formation of a nontrivial vacuum state in heavy ion collisions is predicted in a simplified numerical model. These results have implications for the study of the evolution of the early universe near the QCD phase transition and may be tested experimentally very soon in heavy ion collision experiments. ii Table of Contents Abstract ii Table of Contents iii L i s t of Figures v Acknowledgement 1 2 Introduction 1 1.1 Review of Quantum Chromodynamics and the ^-parameter 2 1.2 Motivation 9 1.3 Overview 12 T h e A n o m a l o u s Effective L a g r a n g i a n and Implications 15 2.1 The Anomalous Effective Lagrangian 15 2.2 Theta Dependence 20 2.3 Vacuum Structure 21 2.3.1 q=l 23 2.3.2 q^l 25 2.4 3 vi Phenomenology of the Pseudoscalar Mesons for Non-Zero 9 25 False V a c u u m Decay 35 3.1 The Domain Wall Solution 38 3.2 Semiclassical Theory 40 iii 3.3 3.4 3.5 4 5 Quantum Corrections at Zero Temperature 42 3.3.1 Positive Eigenvalues 43 3.3.2 Zero and Negative Eigenvalues 53 3.3.3 Decay Rate for Zero Temperature 54 Quantum Corrections for High Temperature 55 3.4.1 Positive Eigenvalues 55 3.4.2 Zero and Negative Eigenvalues 58 3.4.3 Decay Rate for High Temperature 59 Summary 59 P r o d u c t i o n of N o n t r i v i a l T h e t a V a c u a 61 4.1 DCC 62 4.2 Nontrivial #-vacua 64 4.3 Numerical Evolution 67 Conclusions 74 5.1 Results 74 5.2 Future Research 76 Bibliography 79 Appendices 84 A D e r i v a t i o n of the Effective P o t e n t i a l for Q C D 84 B A . l Anomalous Effective Lagrangian for QCD 84 A.2 Determining the values of p and q 91 Hyperspherical Harmonics i n Four Dimensions 94 iv List of Figures 2.1 Branches of the effective potential before the infinite volume limit is taken. 19 2.2 Piecewise smooth effective potential after the infinite volume limit is taken. 20 2.3 Two flavour effective potential 24 2.4 Vacuum energy as a function of vacuum angle 9 24 2.5 Effective potential for q=2 and N = 2 26 2.6 7T° mass as a function of 9 31 f 2.7 TJ mass as a function of 9 2.8 31 Tj' mass as a function of 9 31 2.9 The Euler 9 angle as a function of 9 33 2.10 The Euler <p angle as a function of 9 33 2.11 The Euler ip angle as a function of 9 . . . 33 3.12 Effective potential at 9 = 0 for equal chiral phases rp = <f> = <j> = <f> . . . u d 39 s 3.13 The domain wall solution 40 3.14 Expansion of the partition function 43 3.15 The "potential" U"{<f> ( )) r 4 6 b 3.16 Comparison of the perturbation V (r) pert 3.17 The phase shift <5O(^A) and the approximate perturbation. 48 shown as a function of UJ\ for various values of n. . 49 4.18 Time evolution of \4>k\ shown for various \k\ 4.19 The Fourier distribution of the field is shown at three different times. . . 69 70 4.20 The zero mode and a non-zero mode are shown as a function of time for three different volumes 71 v Acknowledgement I would like to thank my supervisor, D r . A r i e l Zhitnitsky, for his help and guidance. I have learned a great deal through working w i t h h i m and I am grateful to h i m for pushing me to make this thesis the very best it could be. I would also like to thank the other members of my supervisory committee: D r . Janis M c K e n n a , D r . John N g and D r . G o r d o n Semenoff for their invaluable help and comments. I would particularly like to thank D r . M c K e n n a for her advice on other career related matters. I would like to thank L o r i Paniak, Sebastian Jaimungal, K i r k Buckley, M a r k Laidlaw, and Igor H a l p e r i n for their helpful advice and for many stimulating and enjoyable conversations. I would also like to thank all of my other friends here at U . B . C . for making my time i n Vancouver very pleasant. In particular I would like to thank Jaret, Sofia, Greg, Mateo, A l i s o n , K e v i n and all the members of "Where's M y H e a d " . Thanks are due to the entire faculty and staff of the Physics and A s t r o n o m y Department at U . B . C . for providing a stimulating and positive work environment. I would like to acknowledge that this work was funded i n part by the N a t u r a l Sciences and Engineering Research Council of Canada. Words cannot express my gratitude to my family for the constant support and encouragement they have shown me throughout my studies and, indeed, throughout my life. I would especially like to thank my wonderful wife L i s a and my fantastic son N a t h a n for keeping me sane during this demanding process. T h i s thesis is a testament to their love and support. vi Chapter 1 Introduction It is generally accepted that Q u a n t u m Chromodynamics ( Q C D ) is the theory that describes the strong nuclear force. In this theory a l l hadrons are composite structures of particles called quarks, and quarks interact w i t h other quarks v i a exchange of gluons . 1 Quarks have a flavour, an electric charge and a colour charge, while gluons only carry colour charge. The existence of quarks is supported by various experimental results. A s well, a l l the observed hadronic particles can be described as colour singlet combinations of three quarks or quark-antiquark pairs. Q C D can be used to calculate observables i n the high energy (weak coupling) limit, such as in deep inelastic lepton-hadron scattering, and the theory agrees well w i t h the experiments in this regime. Unfortunately, at low energies Q C D exhibits essentially non-perturbative behaviour such as confinement and dynamical breaking of chiral symmetry. A s well, Q C D is asymptotically free, which means that the perturbative weak coupling limit applies i n the high energy regime while the realm of ordinary nuclear matter is the nonperturbative strong coupling l i m i t . Therefore, while Q C D is an excellent candidate for a theory of the strong nuclear force, it is difficult to use it to obtain any quantitative low energy results for ordinary nuclear matter. Anomalous effective Lagrangians are a useful tool for analysing Q C D i n the low energy limit. The anomalous effective Lagrangian used in this thesis makes it possible to study the vacuum structure of Q C D , which could have important implications for study of the As well as by the exchange of the electroweak force carriers x 1 Chapter Introduction 1. 2 evolution of the early universe. As well, it allows us to make predictions that can be tested in heavy ion collisions. Probing the physics of the early universe directly in heavy ion collisions is certainly an exciting concept. We will return to the discussion of effective Lagrangians after a review of Quantum Chromodynamics. 1.1 R e v i e w of Q u a n t u m C h r o m o d y n a m i c s a n d the ^-parameter QCD is an SU(N ) gauge theory defined by the QCD Lagrangian: 2 C 1 N f C = --G%G*r (i-fDf + - m ^) 1 % + G^G?, (1.1) where G^ is the field strength tensor of QCD, W are the quark fields, D^f is the covariant u Q derivative and a is the QCD strong coupling parameter. The dual field strength tensor s is defined as G^" — Sap^G^ where e appbV is the totally antisymmetric four index tensor. The index i labels the quark flavour and the sum on i is over the number of quark flavours, Nf. The Greek indices a and j3 label the quark colour and run over the number of colours, N . The Greek indices p and v are Lorentz indices and run over the spacetime c dimensions. The Roman index a on the field strength tensor is the colour index for the gluon fields and runs over the range: a = 1 . . . N% — 1. The theory of QCD has N = 3 c and this value is used in this thesis even though we leave the symbol N in the formulas c so that N dependence is explicit. We will be concerned with low energy physics so we c will only include the light quarks. In Chapter 2 and Chapter 3 we will include the up, down, and strange quark flavours (Nf = 3) and in Chapter 4 our numerical results only include the up and down flavours (Nf = 2). The field strength tensor in terms of the gluon (gauge) fields is: G% = d,Al SU(N ) 2 C - d Al v - QsUcA^Al, is the group of N x N unitary matrices with unit determinant c c (1.2) Chapter 1. Introduction where g = A-Ka and f 2 s s abc 3 are the structure constants of the gauge group SU(N ). C The quark covariant derivative is: (1.3) where A"^ are the generators of QCD is a SU(N ) C coupled SU(N ) C SU(N ). C gauge theory which means that the Lagrangian is invariant under transformations of the quark and gluon fields similar to electromag- netism. We will not comment further on this fact other than to comment that physical observables must be invariant under these gauge transformations and are said to be gauge invariant. The "extra" term in (1.1) involving the ^-parameter is a total derivative and thus is irrelevant to perturbation theory. However, nonperturbative effects due to instantons indicate that this term should be included in general. We will say more about this ^-parameter after discussing the symmetries of the Lagrangian. In the chiral limit, which is defined by = 0, this Lagrangian is invariant under the transformations: V -+ e x p ^ A M j a (1.4) % -> exp(zAA^ ) \F 5 7 Q where A is an arbitrary real constant and A^- are the generators of the group of Nf x Nf unitary matrices, U(Nf). The transformations correspond to the multiplication of a vec- tor of quarks in flavour space by a member of U(Nf). The invariance under these two types of transformations means that the Lagrangian has an U (Nf)y ector x U (iV/) Axial sym- metry. The U(Nf) matrices can be written as the product of the matrices corresponding to subgroups of U(Nf): U(N ) = SU(N ) ® U(l) f f (1.5) Chapter 1. Introduction 4 where SU(Nf) has already been defined and U(l) is simply an overall phase which also fixes the determinant of the U(Nf) matrix. The symmetry group of the massless QCD Lagrangian is then : 3 Sl^JV/) Vector ® tf(l)vector ® SU(N ) f AxM <g> E/(l) Axial- (1-6) Symmetries of the Lagrangian can be shown to imply the existence of conserved currents. For example, the vector current defined in terms of the quark fields is conserved: <V" = \ ( * 7 * ) = 0. M (1.7) The conserved currents lead to conserved charges: As well, current conservation leads to certain exact identities, called Ward Identities, which are very important and useful. Symmetries of a field theory at the classical level may not be true symmetries of the quantised field theory. Quantising the fields and applying perturbation theory may reveal that the symmetry is not respected by the quantum theory. In this case the symmetry is said to have an anomaly. The associated currents and charges are no longer conserved, but the Ward Identities that come from the original symmetry are not lost, and are simply altered to include the effect of the anomaly. They are then referred to as anomalous Ward Identities. The t/(l) Axial symmetry of massless QCD is anomalous and this anomaly is taken into account in anomalous effective Lagrangians, which we discuss in the next section. As well, some of the symmetries of the Lagrangian may not be respected by the vacuum state. This is called spontaneous symmetry breaking and is indicated by the 3 There is also the SU(N )c [or c 0 gauge symmetry and a scaling symmetry that we ignore here. Chapter 1. Introduction 5 existence of a non-zero vacuum expectation value of some combination of fields, which is called a condensate. Another result is the appearance of a massless particle, called a Goldstone boson, for every broken symmetry. In the QCD vacuum, the condensate: (1.9) spontaneously breaks the ( / ( A / ) i a i subgroup of the symmetry group (1.6). This corr Ax responds to nine broken symmetries and should indicate the presence of nine Goldstone bosons. These symmetries, however, are explicitly broken by the mass terms which we will discuss next. For non-zero quark masses, ^ 0, the QCD Lagrangian is no longer invariant under the transformations (1.4) and we say that the symmetry is explicitly broken . 4 Fortunately the Ward Identities can again be simply altered to include the effects of the explicit breaking of the symmetry. Further, it can be shown that the Goldstone bosons will acquire small masses proportional to the square root of the quark masses and are then referred to as pseudo-Goldstones. The pions, kaons and the 77-meson are identified as the 8 pseudo-Goldstones of the spontaneously broken SU(A^/) Axial symmetries. The 77' meson, however, which should be the remaining pseudo-Goldstone of the broken £7(1) Axial symmetry, is heavier than it should be. This is the famous [7(1) problem. The last term in the QCD Lagrangian (1.1) involving the ^-parameter is related to the vacuum structure of QCD and the axial anomaly. This term is important for this thesis so we give a brief explanation of its origin (for a detailed discussion see [1]). It can be shown that all possible configurations of the gauge field, A^, can be all be classified by an integer, n, called the winding number, which is determined by: ,3 (1.10) For equal masses of the quarks there is a remnant SU(Nf) flavour symmetry that will remain unbroken. This point will be relevant to the meson mixing angles in Chapter 2. 4 Chapter 1. Introduction 6 where Tr denotes the matrix trace over the generating matrices, A of SU(N ). In particua C lar the vacuum configurations of the fields fall into classes labelled by the winding number and these are referred to as winding vacua. However, there exist operators corresponding to physical observables that transform between winding vacua: 0\n) = \n + l), (1.11) which implies that the QCD vacuum should involve contributions from all winding vacua. The 0-vacuum states, which are coherent superpositions of winding vacua: |0) = 5 > " " > > , (1.12) n are the vacuum states that allow no transitions between different vacuum states: (9\0\9') ~<J(0-0')- (1.13) Using these 9-vacuum states in QCD has a nontrivial effect on the Lagrangian which we illustrate next. The vacuum-vacuum transition amplitude in terms of the \0 = 0}-vacua is defined to be: (6 = 0\O\9 = 0) = J VA^V^Vy O exp if d xC 4 , (1.14) where C is the QCD Lagrangian (1.1) without the 9 term. This amplitude can be written in terms of the winding vacua as: (9 = 0\O\9 = 0) = ^2(m\0\n). (1.15) m,n The presence of non-zero 9 leads to an extra phase: (9\0\9) = 2 e - {m\0\n). m,n l{m n)e (1.16) Chapter 1. Introduction 7 It can be shown that: i / ^ G ^ m - n , (1.17) which describes a gauge field configuration which changes winding number over time. Combining equations (1.14), (1.15), (1.16) and (1.17) we obtain the result: (6\0\d) = where we have used J VA.V^Vty g = 47ra! . 2 s s O exp [i j d x {£ + ^ G ^ G T } ] 4 , (1.18) The integrand in the exponent is exactly the QCD Lagrangian we quoted in (1.1). This shows that the general QCD Lagrangian should include the 0-term. The 9 term is connected with the axial anomaly. The axial anomaly can be restated as the fact that the axial current: = E * W * \ (1-19) i=i is not conserved. It can be shown that the divergence of the axial current in the quantum theory is given by: « = -krG^G't. (1.20) As was mentioned earlier, G^G^ is a total divergence, which means it can be written: G% G^ = d^. / (1.21) is defined in terms of the gauge fields, A^, and is not gauge invariant. This means that we can define a new conserved current: J5, = J 5 , - ^ K „ (1.22) with a corresponding conserved charge: Q = Jd?x J . 5 5 0 (1.23) Chapter 1. Introduction 8 It can be shown that this charge generates changes of the 9 parameter: e \9) = iAQs (1.24) \e-2N A), f for arbitrary real constant A. This does not contradict our earlier comments because Q 5 is not gauge invariant and therefore does not correspond to a physical observable. If all the quarks are massless then the conservation of the anomalous axial current means that the ^-dependence of the physical theory can be removed by a Axial rotation. The 9 term in ( 1 . 1 ) breaks some discrete symmetries of the rest of the Lagrangian. Without this term, the Lagrangian is invariant under the charge*conjugation (C) transformation, where particle and anti-particles are switched (^ o ty), and the Parity (P) transformation where all spatial vectors are reflected through the origin (x —> —x). The term involving the ^-parameter explicitly breaks P and CP invariance and we can say that the 9 parameterises the amount of CP violation in the strong interaction. Stringent limits have been placed on the value of the ^-parameter in the universe today by measurements of the neutron electric dipole moment[2] which predict a value of 6 < 1 0 ~ 1 0 [3]. The problems of the ^-parameter are not eliminated, however, as the ^-parameter in the full theory will contain, for instance, in addition to any fundamental parameter, contributions from the determinant of the quark mass matrix(see [ 1 ] ) : 9 = 9 QCD + arg(det M). (1.25) The strong CP problem can then be stated in the following way: "Why do completely unrelated contributions to the ^-parameter cancel each other out to such a high degree of accuracy?". A solution to this problem can be found by promoting the ^-parameter to the status of a dynamical field called the axion. This field can be shown to dynamically relax to zero, providing a solution to the strong CP problem[4]. The axion is an extremely Chapter 1. Introduction 9 weakly interacting particle and has so far evaded detection despite the best efforts of many physicists. It is important to note, however, that this research does not require the existence of the axion. In fact, even if there is no dynamical axion to solve the strong CP problem, we will argue in Chapter 4 that it is possible to create regions of space where dynamical fields will behave as if the ^-parameter were non-zero. This "induced" ^-parameter could differ significantly from the ^-parameter in the universe as a whole. This type of non-zero ^-parameter is much more likely to be observed in heavy ion collisions and is perhaps a more important motivation for considering non-zero values of 9. We will not distinguish between the two types of 9 for most of the results of this thesis. 1.2 Motivation As was mentioned previously, it is difficult to use the full theory of QCD to calculate observables in the low energy regime. However, in this regime we are able to use effective Lagrangians which have proven to be a very powerful tool in Quantum Field Theory. These effective Lagrangians are not meant to be a complete theory, but are used to study the behaviour of the low energy degrees of freedom of QCD in order to understand the theory in this regime. There are two main types of effective Lagrangians in use in Quantum Field Theory. The first type is the Wilsonian effective action which describes the low energy dynamics of the lightest particles in the theory. The heavier particles are in some sense "integrated out" as can actually be done in some quantum field theories. However, in many cases (such as QCD), it is not known how to do this and one simply writes down an effective Lagrangian based on constraints such as symmetries. Effective Chiral Lagrangians are of this type and are very useful in their domain of application. Chapter 1. Introduction 10 Another type of effective Lagrangian is defined as the Legendre transformation of the generating functional for the connected Greens functions. This formulation implements (at the Lagrangian level) anomalous Ward identities relating vacuum condensates of the fields and is referred to as the anomalous effective Lagrangian. This type of approach is very useful in studying vacuum properties of the theory and is the approach used in this thesis. The utility of the second type of approach was recognised long ago for supersymmetric gauge theories for both the pure gauge case [5] and for full supersymmetric QCD[6]. More recently it has been applied to non-supersymmetric Yang Mills theory[7] and generalised to QCD[8, 9]. The Witten-Veneziano-Di Vecchia effective chiral Lagrangian[10], which was studied even earlier, represents a hybrid of the two types of effective Lagrangians. It is based on symmetry properties, as are effective chiral Lagrangians, but it also represents the first attempt to include the effects of the anomaly which leads to the second type of approach. It should be stressed that the form of the anomalous effective Lagrangian used in this thesis is only one possible candidate for a low energy effective Lagrangian for QCD. It is certainly a promising candidate, but we can not say with certainty that it is the only one. Study of this anomalous effective Lagrangian approach is interesting for four reasons. First, this effective Lagrangian provides a generalisation of the large N Di Vecchiac Veneziano-Witten effective chiral Lagrangian ( W W ECL) for arbitrary 7V . The W W C ECL was determined purely from symmetry considerations and inclusion of the anomaly. The Lagrangian that we use is a generalisation in the sense that starting from this Lagrangian, integrating out the heavy degrees of freedom, and then passing to the large N limit, we reproduce the W W ECL. Furthermore, this approach to the derivation of c the W W ECL fixes all dimensional parameters in terms of the measurable quark and Chapter 1. Introduction 11 gluon condensates. I cannot stress enough that we are not saying that the W W ECL is wrong. The effective Lagrangian we use is simply a possible generalisation of the W W effective chiral Lagrangian to finite values of N . c Second, this approach makes it possible to address the problem of the ^-dependence of QCD. The ^-parameter is the strong CP violation parameter which was discussed in the previous section. The form of the effective Lagrangian that we use solves a number of puzzles relating to the ^-dependence of physical observables in general and the energy dependence of a nontrivial ^-vacuum. The ^-dependence of QCD is directly related to the axion potential since the axion arises by promoting the ^-parameter to the status of a dynamical field. Third, this effective Lagrangian allows one to address the phenomenology of the light pseudoscalar mesons and the n' meson without any further phenomenological input. We obtain important results concerning the rj' mass and we are able to calculate the masses of the light neutral pseudoscalar mesons (IT and rj). We obtain values for the r\ — rj 0 mixing angle for unbroken SU(2) flavour symmetry and 7r° — 77 — 77' mixing angles for fully broken flavour symmetry. As well, we generalise the above results to include the effect of a non-zero ^-parameter. Finally, we show that the mesons cease to be pure pseudoscalars in the presence of a non-zero 0-parameter. The ^-dependence in the phenomenology of the pseudoscalar mesons is particularly important for the detection of some of the nontrivial vacuum states that we will discuss next. The final, and perhaps the most important reason, for studying this effective Lagrangian, is that it predicts the possibility of a large variety of metastable vacua in certain regions of parameter space. Nontrivial vacua have been shown to exist in Yang Mills theories in the large N limit using the AdS/CFT correspondence[ll]. The same phenomenon c was observed in the analysis of soft breaking of supersymmetric models[12]. The existence of nontrivial vacua appears to be a general phenomenon for gauge theories in the strong Chapter 1. Introduction 12 coupling limit. The metastable vacua could have played an important role in the evolution of the early universe at around the time of the QCD phase transition. The possible effects include dark matter (axions [13] and MACHOs [14]), baryogenesis[15, 16], the formation of large scale magneticfields[17]and many others[18]. Therefore, understanding the physics of these nontrivial vacuum states could be important in understanding the early evolution of the universe. What is even more exciting is that most of the results of this thesis are potentially testable in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratories. One of the results presented in this thesis suggests that a nontrivial 9-vacuum state is potentially observable at RHIC. The ^-dependence of the phenomenology of the pseudoscalar mesons discussed above provides a number of ways to look for this ^-vacuum state. Probing the physics of the early universe directly in heavy ion collisions is certainly an exciting possibility. The reasons for believing that this effective Lagrangian might be physically meaningful are also fourfold. First, by construction, it reproduces the chiral and conformal Ward identities of QCD. Second, it reduces to the W W ECL in the appropriate limits as was mentioned above. Third, it reproduces the known dependence of physical observables in 9 for small values of 9, but also leads to 2TT periodicity of physical observables in 9. Fourth, the related effective Lagrangian for gluodynamics [7] has the desirable property: d E (9) d9 2k vac 2k 2k (1.26) which was argued for by Veneziano in order for the 17(1) problem to be resolved [19]. 1.3 Overview The vacuum properties of QCD are studied starting from an anomalous effective Lagrangian. First, the form of the anomalous effective Lagrangian is discussed. Chapter 1. Introduction 13 Second, the ^-dependence of the vacuum of QCD is studied in this approach and shown that it has the correct ^-dependence for small 9 and that it also produces 2ix periodicity in 9 of physical observables. Third, the phenomenology of the pseudoscalar mesons is studied. We obtain a mass formula for the rf meson which is a generalisation of the mass formula from [8] to non-zero values of 9. We obtain values for the masses of the neutral pseudoscalars for arbitrary values of 9. We also obtain a formula for the r\ — rf mixing angles and a method of calculating the ir — r\ — rf mixing angles for non-zero 9. Finally we can show that for 0 non-zero values of 9 the pseudoscalar mesons cease to be pure pseudoscalars, but acquire scalar components. Fourth, the decay rate is calculated for a false vacuum at 9 = 0 consisting of a chiral fields at a nontrivial minimum of the effective potential. We nonperturbatively calculate the decay rate per unit volume, T/V, for both the zero and high temperature theories. In addition, the method used in these calculations represents a new approach that could be useful in other calculations of this type. Finally, a numerical study is presented of the evolution of the light degrees of freedom in regions of non-zero 9 in order to show that it is possible to form a nontrivial 9-vacuum state in relativistic heavy ion collisions. The goal of this thesis is develop the understanding of the vacuum structure of QCD in this particular effective Lagrangian approach and make predictions about observable consequences that could be tested in heavy ion collisions. In Chapter 2 the anomalous effective Lagrangian is discussed and the ^-dependence of the vacuum, the vacuum structure and the phenomenology of the pseudoscalar mesons is studied in this theory. In Chapter 3 we calculate the decay rate of a chiral false vacuum for the low temperature theory and the high temperature theory. In Chapter 4 we discuss a numerical calculation performed to demonstrate the possibility of producing Chapter 1. Introduction 14 a nontrivial 0-vacuum state in a heavy ion collision. Chapter 5 discusses conclusions and areas for future research. Chapter 2 T h e A n o m a l o u s Effective L a g r a n g i a n a n d Implications The purpose of this chapter is to introduce the anomalous effective Lagrangian for QCD and discuss some of its implications . In the first section we discuss the form of the 1 anomalous effective Lagrangian. In the second section we discuss the ^-dependence of the vacuum energy. In the third section we discuss the structure of the QCD vacuum. In the fourth section we discuss the phenomenology of the light neutral pseudoscalar mesons and the rj' meson, including the masses and mixing angles and the appearance of a scalar component of the mesons for non-zero 6. 2.1 T h e A n o m a l o u s Effective L a g r a n g i a n The construction of the effective Lagrangian used in these studies was already accomplished before this research. The details of the effective Lagrangian necessary for the understanding of this work are presented here. For completeness, a review of the derivation by Zhitnitsky and Halperin[7, 8] is included in Appendix A. A useful tool for studying the properties of the QCD vacuum in the presence of the ^-parameter is the anomalous effective Lagrangian for QCD. The anomalous effective Lagrangian is defined as the Legendre transform of the generating functional for zero momentum correlation functions. In QCD the relevant operators are G^G^ (topological density operator), G = G^G^ (gluon scalar current) andtyty(quark condensate). The 2 correlation functions of both gluon operators are fixed in terms of the gluon condensate: 1 Many of the results of this chapter were published i n 15 [9]. Chapter 2. The Anomalous Effective Lagrangian and Implications 16 [20, 21]. Actually, only the potential part of this Lagrangian can be constructed (a /irG ) 2 s as this is the only part that can be fixed by low energy theorems relating the zero momentum correlation functions. The potential part determines the whole Lagrangian, however, once the standard kinetic term is inserted, and so we will use the terms effective potential and effective Lagrangian interchangeably. The most important point is that the effective potential constructed in this way will contain all the information provided in known low energy theorems of QCD. In the effective Lagrangian approach the light matter fields are described by the 2 unitary matrix, Uij, corresponding to the phases of the chiral condensate: (2.27) with 3. ir \ a a U = exp tV2 — r z- h fn UU* = (2.28) 1 , ir is the octet of pseudoscalar fields SU(Nf) singlet pseudoscalar field (meson). SU(Nf), (pions, kaons and the eta meson) and rf is the We use the values f rf yJNf frf J where A are the Gell-Mann matrices of a 2 a = 132 MeV[22] and fa = 86 MeV[20] for the meson decay n constants . The indices, % and j, of U^ run over the quark flavour index, i,j = 1 4 Nf, which correspond to the quark flavours up(u), down(d) and strange(s), respectively. The effective potential in terms of these fields is given by: W (U, eff U) + = - Jim ^°° v i v log I £ VEcos - - ( 0 - « l o g d e t U) + —l exp V \l=0 + V ^VTr(MU + MW) j , (2.29) Note that the rf is not really very light, but it enters the theory in this way and the mass of the rf is one of the topics that we address in this thesis. Note that mixing of the flavour eigenstates is ignored at this level. See Section 2.4 for more details. It should be noted that what we refer to as f /V% is sometimes denoted f in the literature. As well, what we refer to as f >, is in the notation of [20]. 2 3 4 n v w Chapter 2. The Anomalous Effective Lagrangian and Implications 17 where V is the total four volume, M = d i a g ( m j | | ) and E = ba /(32ix)(G ) 2 5 s , b = y i V — |7V/ is the first coefficient of the Gell-Mann - Low /^-function of QCD (see c Appendix A) and Tr denotes the matrix trace. The sum over the integer I arises in the full derivation (see Appendix A) and is a sum over different branches of the potential 6 (see below). It is necessary in order that the full effective potential, including the gluon degrees of freedom, be single valued and bounded from below. The sum remains after the heavy gluon degrees of freedom are integrated out. The prescription of summing over all branches was first suggested in [23] and applied to the anomalous effective Lagrangian of [5] to cure problems similar to this case. As well, this prescription is consistent with the results of [11] where the problem of ^-dependence is studied using the AdS/CFT correspondence. The parameters p and q are relatively prime integers of the order: p = 0(N ) and C q = 0(N C °) = 0 ( 1 ) . Different approaches to fixing these parameters give different results. For completeness, the details of these approaches are discussed in Appendix A. Here we simply give the two sets of parameters that have been derived: • q= 1 and p — N, • q= 8 and p = 3b c QCD = UN - 27V), C where bq.cD is the first coefficient of the beta function of QCD. The first set of values is perhaps more natural and is in agreement with supersymmetric theories, but there are some difficulties with the arguments for transferring these results unchanged to nonsupersymmetric theories . The second set of values has been derived from a number of 7 We use the values (a /-n G ) = 0.012 GeV , (#^$ ) = (# * ) = -(240 MeV) and ( * * ) = -0.8(240MeV) , except in Chapter 3 where we use ( $ * ) = -(240MeV) . If a single value for I was used in this effective potential, instead of summing over I, then the effective potential would simply be the argument of the exponential for a specific value of I. This is what is meant by the branch of the effective potential and is labelled by the integer I. See Appendix A for brief details and references. 5 2 4 u d d 3 s 3 6 7 s s 3 s s Chapter 2. The Anomalous Effective Lagrangian and Implications 18 different approaches and cannot be ruled out. The question of which set of parameters is correct is very significant as different values of these parameters can lead to some different qualitative results. Most of the results of this thesis do not depend on specific values of these parameters. When results do depend on specific values, this fact will be pointed out and discussed. The effective potential (2.29) may not be represented by a single analytic function in the V —> oo limit. The thermodynamic limit selects, at each value of 6 — i logdet U, one particular branch (ie. a particular value of the integer 7) and cusp singularities occur where the branches coincide. The different branches of the effective potential before the infinite volume limit is taken are shown in Figure 2.1. The infinite volume limit picks the lowest energy branch out from all the other branches of the effective potential (see Figure 2.2). The cusps are present because level crossing occurs and the branch which is the lowest energy changes due to this crossing. The present form of the effective potential (2.29) is most conducive to discussing the topic of chiral rotations. Notice that one can perform the global transformation: U -> diag(e^) U, (2.30) to transform to a Lagrangian of the same form since the kinetic term of the form, T r ^ f / c ^ C / ) , will be unaffected. This fact was used in [10] to introduce the ^-parameter -1 into the effective potential. This fact was also used in [10] to show that if any of the quark masses vanish, the ^-parameter can have no effect on the physics. This result will be briefly touched on later in this chapter. Finally, we will use a chiral rotation of this type in Chapter 4 to introduce the concept of the induced ^-parameter. Now we return to the effective potential for the light matter fields which, when we take U — diag(e ) l(t>i 8 8 , becomes a piecewise smooth potential for the phases of the chiral T h i s is n o t a r e s t r i c t i o n since the q u a r k mass m a t r i x c a n always be d i a g o n a l i s e d . Chapter 2. The Anomalous Effective Lagrangian and Implications 19 w,eff -71 / \ / \ \ / / \ / \ 4> K Figure 2.1: Branches of the effective potential before the infinite volume limit is taken for 9 = 0, Nf = 1, q — 5 and p = 3b cDQ condensate, with cusp singularities: V cos</>i , / = 0,1,... ,p - 1, P P P (2.31) J if ( 2 / - i ) ^ <e-Y <f> <(2i i i + i)^ . (2.32) The chiral phases of the quark condensate are denoted by (pi where the index ' i ' runs over the quark flavour index (ie. u, d, s). They are also referred to as the light chiral fields and are related to the pseudoscalar meson fields by recalling Eq. (2.28). At this point it is useful to give some indication of what the parameters p and q mean in terms of the effective potential. The parameter 'p' is the number of branches of the effective potential. It also corresponds to the number of distinct vacua before the infinite volume limit is taken. The parameter V determines the number of distinct vacua with different energies that exist before the infinite volume limit is taken. For q even there are (q + 2)/2 different energy vacua; for q odd there are (q + l)/2 different energy vacua. Chapter 2. The Anomalous Effective Lagrangian and Implications w 20 eff Figure 2.2: Piecewise smooth effective potential after the infinite volume limit is taken for 9 = 0, Nf = 1, q = 5 and p = 3&QC£>- For q — 1 all 'p' vacua are degenerate in energy. The parameter '<?' is also the number of physically distinct vacua that exist at 9 = 0 in the infinite volume limit. This form of the effective potential (2.31) for the light chiral fields is the starting point of all the research of this thesis. 2.2 T h e t a Dependence i The effective potential (2.31) for the phases of the chiral condensate described in the last section is a piecewise smooth potential with cusp singularities. The form of this potential is important because it satisfies two seemingly contradictory requirements. The first requirement is that for pure gluodynamics the vacuum energy should depend on 9 only in the specific combination 9/N in order that the U(l) problem be resolved. With c p = 0(N ) the effective potential appears to have the correct form and in fact has C exactly the desired form if we use the beta function of pure gluodynamics. The second Chapter 2. The Anomalous Effective Lagrangian and Implications 21 requirement is that since 9 is an angular variable, all physical observables should be periodic in 9 with period 27r. This requirement seems at first glance to be impossible to meet when the first requirement is met. This potential does satisfy both these properties, however, and the piecewise smooth structure turns out to be crucial (see Figure 2.2). The different branches of the effective potential each locally depend only on the combination 0/N , but the fact that in the infinite volume limit only the lowest energy state survives c restores the 2-K periodicity of physical observables. 2.3 V a c u u m Structure From the effective potential above (2.31) one can derive the equations of minimisation for each of the chiral fields in the I — 0 branch: sin ( ^ - ^ . ) = ^ f s i „ ( « . (2-33) They coincide with those of [10] to lowest order in 1/N and can be used to find the C solutions for the vacuum values of the chiral fields as a function of 9. For general values of Si = an analytic solution is not possible. For three light quarks pMi/qE the realistic case, e ,Ed <C l,s u s (Nf = 3) in ~ 1, approximate solutions can be found as solutions of the equations: 0f 4 +4 = 0 ) = 0, 0 ) = = 0 ) e sin^ (2.34) = e sin0j / = 0 ) u =o) d , where 0(e ,Ed) terms have been neglected in (2.33). The solution to these equations is: u • j,(i=o) _ S m 0 M = m sm9 d [ml + ml + 2m m u • JI('=O) d m sin0 u cos 9]V* + ^ £ d ) ' Chapter 2. The Anomalous Effective Lagrangian and Implications sinrf=°> = 22 0(s ,e ). u d This solution for the I — 0 branch agrees with the solution in [10]. Notice that for 9 = 0 the solutions are fa = 0. Finally notice that for m-i = m <C AQ D C the lowest energy state is: (2.36) which, when substituted back into the effective potential for I = 0, produces a formula for the vacuum energy as a function of 9: (2.37) This formula is the vacuum energy of the true vacuum state in a universe with a given value of 9 and it is independent of the values of q and p. This relation actually has a very deep significance. It shows that the energy difference between the 9 = 0 vacuum and the 9^0 vacuum is proportional to the parameter Mi, in agreement with a general theorem that physical ^-dependence appears only in combination with m and goes away in the q chiral limit. It is also important for the results in Chapter 4 that the 0-vacuum energy 9 is not proportional to the parameter E, as suggested by the 0-dependence of (2.31), but is proportional to the smaller parameter Mj. This means that the formation of nontrivial 0-vacuum may be possible in heavy ion collisions. It must now be stated that the exact form of the effective potential is not crucial to the 0-vacuum energy dependence. The crucial feature is that the ^-parameter occurs only in the very special form: (9 — Y, <t>)- This feature is also present in [10] and in both cases we see that the large parameter in front of the term containing this combination does not prevent the formation of a state with nontrivial 0-value. Instead it means that the state will have nontrivial values of the chiral phases which will depend on 9 and the vacuum energy will depend on Mi oc m . q 9 T h i s was s h o w n i n [10]: p h y s i c a l ^-dependence vanishes i f any q u a r k mass goes t o zero. Chapter 2. The Anomalous Effective Lagrangian and Implications 23 For further discussion of the vacuum structure we will distinguish between two situations: q = 1 and q 2.3.1 1. q=l In the case where q = 1 we find that the solution of the corresponding equation of minimisation for Z = 1 does not lead to a minimum of the effective potential but instead leads to a saddle point. This is illustrated in Figure 2.3 for 9 = 0 and Nf = 2. The solutions for all other values of I coincide with one of these solutions (modulo 27r). Therefore there exists a unique vacuum. We have checked numerically, using physical values of the quark masses in the case Nf = 3, that for all values of 9 there is only one vacuum, in agreement with [10]. In particular, at 9 = TT no extra metastable vacua appear. When the masses satisfy Dashen's constraint[24]: mm u d > m \m s d - m \, u (2.38) the vacuum doubling phenomenon—the appearance of a metastable vacuum state— occurs for values of 9 near TT. This metastable state is a "double" of the original vacua because at the critical point, 9 = TT, the two vacua are exactly degenerate. These results agree with the results of [10]. For three (Nf = 3) light equal-mass quarks we find that the metastable state exists for 9 values from 7r/2 to 37r/2. This is in agreement with the results of [25]. The vacuum energy is shown as a function of 9 in this case in Figure 2.4. None of the results of this section are new, but they show that the form of the effective potential that we use reproduces previous results [10, 25] for the case q = 1. Chapter 2. The Anomalous Effective Lagrangian and Implications 24 Figure 2.4: Vacuum energy as a function of vacuum angle 9 for the case q = 1 with Nf = 3 light flavours of equal masses. The solid and dashed lines represent minima and the dotted line is the intermediate saddle point. Chapter 2. The Anomalous Effective Lagrangian and Implications 2.3.2 25 q^l Interesting vacuum structure arises when we allow q ^ 1. Notice that the effective potential (2.31) has the property that <> / —)• <j> + 2TT simply corresponds to a relabelling of the branches 1—^1 + q. Therefore for 9 = 0 we have q physically distinct vacuum states corresponding to the minima of the branches I = 0,..., q — 1. For illustrative purposes, consider the specific case q = 2, Nf = 2 and 9 = 0. The effective potential for this case is illustrated in Figure 2.5. Notice that there are q = 2 minima of the effective potential: one stable vacuum solution and one metastable. This behaviour carries over to all values of q, with q — 1 metastable vacua and 1 stable vacuum. This is a new feature of the theory and leads to the phenomenon of false vacuum decay which is discussed in the next chapter. One desirable feature of this potential is that the vacuum doubling phenomenon, in which an extra metastable vacuum state appears and becomes degenerate with the highest energy metastable vacuum state at the critical values: 6 = (2k + l)^, fc = 0 , l , . . . , p - l , (2.39) occurs for all values of the light quark masses. The extreme sensitivity of the theory to the values of the light quark masses is avoided when q ^ 1. 2.4 P h e n o m e n o l o g y of the Pseudoscalar M e s o n s for N o n - Z e r o 9 In this section we specialise to the case of Nf = 3 in order to calculate the ^-dependence of the 7T°, rj and r( masses and the r] — rf and 7r° — 77 — 77' mixing angles. As well we show that the pseudoscalar mesons, which are pure pseudoscalars for 9 = 0, acquire scalar components for non-zero values of 9. Chapter 2. The Anomalous Effective Lagrangian and Implications 26 Figure 2.5: Effective potential for q=2 and Nf = 2. The point indicated as a) is the global minimum. The two local minima indicated as b) are identified. Chapter 2. The Anomalous Effective Lagrangian and Implications 27 T h e chiral condensate is parameterised in the following way: r9 _ h I- U = UQ exp iV2 — fix ' (2.40) yjNf fv' where U is the solution to the minimisation equations, A are the G e l l - M a n n matrices a 0 of SU(Nf), 7 r is the pseudoscalar octet, 7r is the singlet field that is identified w i t h a 9 the rf i n the absence of mixing. The 7r and 7r fields have vanishing expectation values. a 9 Inserting this U into E q . (2.29) and expanding this term to quadratic order in the fields we find that the linear terms vanish and we obtain the mass matrix from the terms of the form: 1 £ M ^ a,6=3,8,9 (2.41) 7r 7r . a f l 6 6 Note that the indices do not run over the quark flavour group but over the U(3) indices associated w i t h the G e l l - M a n n matrices . T h e 7r , TT and 7r fields are associated w i t h 10 3 8 9 the G e l l - M a n n matrices: 1 ( A, = 0 0^ 0 - 1 0 0 0 0 1 0 0 1 0 0 \ ^ 1 0 An = which are matrices w i t h respect to the quark flavour basis. (UM — dd)/y/2, (uu + dd — 2ss)/y/6 and (uu + dd + ss)/\/3, 0^ 0 1 0 0 0 (2.42) 1 These fields behave like: respectively. These are the flavour eigenstates usually associated with the neutral mesons 7r°, n and n', respectively, in the absence of mixing. The components of the mass m a t r i x are: m 3 3 = 2 f 2 2 (M cos (j> + M u u cos <f> ), d d J 7T 2 (M cos 4> + M cos <f> + 4M cos 4> ), u u d d S s 3/? We include in this term the identity matrix, denoted by A , which is not properly a Gell-Mann matrix. 10 9 Chapter 2. The Anomalous Effective Lagrangian and Implications 28 m;99 3S = m ™39 = m m 2 83= 89 = m v3# 2 = d u 9 3 = \ITT ~T7~ ( 98 (2.43) M V Nf m ( u COS (j) - M COS (j) ) , MU d ~ C0S MD ' C0S Ujrf \ r^r~F~r 2 cos ( Mu -0« + d M C fa - °S 2M S COS 0), S V 37V/ J^Jjji where the fa are the vacuum values of the particular vacuum about which we are expanding. Notice that the 77199 component depends on 9 explicitly and a l l of the components depend on 9 i m p l i c i t l y through the vacuum values of fa. Neglecting mixing effects, mg coincides w i t h the physical mass of the rj. 9 In fact for 9 = 0 and particular values, p = 36QCD and q = 8, we reproduce exactly the relation derived in [8]: i f m * k ^ Nf< = ° 2 ~ w ^ ) / m { q > i % ) + ° { m ^ • ( 2 - 4 4 ) u,d,s T h i s relation shows that the 77' is heavier than the other pseudoscalar mesons because of the contribution from the gluon condensate. contributes 20-40% of the mass of the rf. of 0(m ) 2 m It also predicts that the strange quark This suggests that perturbative corrections (see below) could be quite significant. Taking the formal limit N c —y 0 leads to the mass relation in [21] where m , scales like Nf/N , 2 q c w i t h [26]. If we instead take the limit N c m , = 0(m ), -> 0 0 and i n agreement —> 0 0 w i t h fixed quark mass m q ^ 0, then as for ordinary pseudo-Goldstone bosons. T h e mass relation (2.44) gives 2 q the numerical result for m g 9 ( ~ m^/)=900 M e V . T h i s result is encouraging as it is of the right order of magnitude, (see below) and is the first connection between the unrelated quantities, m i and ( G ) . 2 1 1 v The fact that the mass m a t r i x M ab defined by (2.41) has off diagonal terms, as seen in (2.43), means that the physical fields are not the TT fields, but linear combinations of 0, The gluon condensate is measured in the charmonium system which is totally unrelated to the rj'. Chapter 2. The Anomalous Effective Lagrangian and Implications 29 them determined by the eigenvectors of the mass matrix. The eigenvalues of the mass matrix give the masses of the physical mass eigenstates corresponding to the 7r°, 77 and 77' mesons. For 9 = 0, p = iV , q = 1, and the standard values m = 4 MeV, m = 7 MeV c u d and m = 150 MeV we obtain the values: s m,o = 131 MeV, m, = 412 MeV, m , = 1015 MeV. v (2.45) Comparing these values to the measured masses of 135 MeV, 547 MeV and 958 MeV respectively, we see that the mass values are all of the right order of magnitude and in approximately the right ratios to one another, which is very encouraging. These values were obtained using specific values for gluon and strange quark condensates, the values of which are fairly uncertain. They could be off by 20% or more and this would have a significant effect on the masses. A change to either of the condensate values would have a similar relative effect on each of the 77 and rj' masses, but little effect on the 7r° mass. In addition these values could be altered significantly by higher order perturbative corrections. An analytical form for the masses of the neutral mesons as a function of 9 would be extremely messy, but is obtained very easily numerically, and the results are shown in Figures 2.6, 2.7 and 2.8 (Note the different scales on the vertical axes.). We can see that there is a mild ^-dependence in the masses of the rj and 77' and a strong ^-dependence in the mass of the 7r°. This can be easily understood from the components of the mass matrix (2.43). The masses of the 77 and 77' are predominantly determined by the diagonal terms m | and m ^ . The 7 n | term is dominated by the M term, because the strange 8 8 s quark mass is so much larger than the up and down quark masses. This term is very weakly dependent on 9 through the cos((f> ) term as can be seen by recalling the vacuum s solution (2.35). The variation of the ml term with respect to 9 is similarly determined 9 by the M term because the term involving the larger parameter, E, is constant since s Chapter 2. The Anomalous Effective Lagrangian and Implications 30 the vacuum solutions have (Yi 4>i — 9) = 0. T h e 7T° mass, however, is predominantly determined by the m | term which is determined solely by the up and down quark terms 3 which are strongly dependent on the ^-parameter as illustrated by the vacuum solutions (2.35) for (j) and (j) . T h i s 0-dependence of the masses is one of the possible signals of u d the formation of a nontrivial 0-vacuum discussed i n Chapter 4. In the l i m i t of pure SU(2) flavour symmetry (m u = m d = m, which also means 4>u = <t>d) we can consider the simplest case of 77 — rj mixing. For 9 = 0 the m i x i n g matrix is: ( —^(2rn (ss)+rn(uu)) m 77-77- V ^(m (ss) s - m(uu)) s \ g j y ; (m (ss) — m(uu)) s - ^ ^ ( ^ )G + 7 ' which can easily be shown to coincide with an accuracy 0(m ) 2 , (2.46) 2 J w i t h the m a t r i x given by Veneziano[26]: { r)-r)' — \(\m\-m\) (2.47) m V 2V2, where we replace the topological susceptibility, x, evaluated i n pure gluodynamics, w i t h the value i n Q C D proportional to the gluon condensate. It should be stressed that the fact that both sides of the relation (2.44) depend on quantities evaluated i n Q C D represents a significant improvement i n our understanding. A Q C D analog of the Witten-Veneziano formula[26] can be written down for the particular values p = 36QCD and q = 8: ,OL 36 = -f* K, + m - 2774) + 2 0(m\ (2.48) T h i s equation is a generalisation of the rj mass relation (2.44) to include the effects of 77 — rj mixing. The 77 — 77' m i x i n g angle as a function of 9 is: Chapter 2. The Anomalous Effective Lagrangian and Implications TC° m a s s (MeV) 0 0 -n/2 -Tt n/2 TC Figure 2.6: 7r° mass as a function of 9 r\ m a s s (MeV) 4JJ. 410 \ \ 40 8 406 / \ 402 / -TC \ 404 e n/2 -n/2 Figure 2.7: r\ mass as a function of 9 r\' m a s s (MeV) < y'l014 1013 1012 \ \ \ 1011 1010 \ -TC -TC/2 Tt/2 Figure 2.8: 77' mass as a function of 9 31 Chapter 2. The Anomalous Effective Lagrangian and Implications where m | is given above i n (2.43) and, ra / and m , are eigenvalues of the mass m a t r i x 2 9 32 2 for a given value of 9. For 9 = 0 and fa given by the vacuum solutions i n this case, we obtain an r\ — 77' m i x i n g angle of about —22°. T h i s value is reasonable as phenomenological estimates from experimental data range from —10° to —23°[27]. A s well, we can study the full case of TT° — TJ — n' mixing angles when the SU(2) flavor symmetry is broken. The mass eigenvector coordinate system, TC — i) — fj' is defined by 0 the eigenvectors of the mass matrix. T h e orientation of the mass eigenvector coordinate system w i t h respect to the SU(3) flavor eigenstate system, #3, 7r and 7r , can be specified 8 by Euler angles. For the specific case where 9 = 0 and m u 9 = 4 MeV, m d = 7 M e V and m = 150 M e V we obtain the Euler angles[28]: s 9 uler E = l9°, (f>Euler = -0.83°, ^ Euler = 2.5°. (2.50) These values show that at least for 9 = 0 the m i x i n g of the TT° meson w i t h the TJ and 77' is fairly small and breaking of the SU(2) flavor symmetry does not have a large effect on the 77 — 77' mixing. In general, the Euler angles can be found numerically and they are shown as a function of 9 in Figures 2.9, 2.10 and 2.11. We can see that there is a significant dependence on 9 in the case of the full 7r° — 77 — 77' mixing. A s i n the case of the masses, the m i x i n g angles would likely be significantly affected by changes to the values used for the gluon and strange quark condensates and by the inclusion of perturbative corrections. The final topic of this section is the behaviour of the pseudoscalar mesons under the P a r i t y Operation (P) for non-zero values of 9. For 9 = 0 they are eigenfunctions of Parity w i t h eigenvalue —1, which means they are pure pseudoscalars. For 9^0 they cease to be eigenstates of P a r i t y and become mixtures of scalar and pseudoscalar states . 12 Recalling the definition (2.27) of the chiral condensate matrix, U, and the diagonal 1 2 T h i s fact was p r e v i o u s l y noted i n [29]. Chapter 2. The Anomalous Effective Lagrangian and Implications 33 Euler Angle 9 (°) Figure 2.9: T h e Euler 9 angle for the mass eigenvector coordinate system w i t h respect to the flavor eigenstate coordinate system as a function of 9. Euler Angle cp (°) -l •* \ -1.5 '• -n/2 e n/2 -2.5 Figure 2.10: T h e Euler (f) angle for the mass eigenvector coordinate system w i t h respect to the flavor eigenstate coordinate system as a function of 9. Euler Angle \|/ (°) 9 g • • 7 • 6 5 \ -71 4 y ^71/2 7C Figure 2.11: T h e Euler ip angle for the mass eigenvector coordinate system with respect to the flavor eigenstate coordinate system as a function of 9. Chapter 2. The Anomalous Effective Lagrangian and Implications 34 form of the vacuum solution, the quark condensates in the 9-vacuum: (* *' )e = (*Li*Rj)o=o e*', Li (2.51) Rj can be represented as a chiral rotation of the usual vacuum: $' = U $Ul = (ty ty ) =o , R Ll , Rj e (2.52) with: {U ) R = lk foe** (Ul) = 8 k3 k i e& . (2.53) The quark fields transform under this rotation in the following way: = (UR)* * R , < K 3 = *L {U ) ] k] • (2.54) In the chirally rotated basis, the spin content of the pseudoscalars is the same as in the standard basis. However, in terms of the unrotated basis the spin content will generally have a very different form. Consider for example the 7r° field in the 9-vacuum: |7T°) ~ \u'i^ u' — d'i~f d') = cos^ {ui^u) — s i n 0 \uu) — (u 4-> d). 5 5 u u (2.55) The second term in this equation, which vanishes for 9 = 0, is the additional scalar component which the 7r° field acquires for non-zero values of 9. The other pseudoscalar particles undergo a similar transformation. The scalar component of the pseudoscalar mesons for non-zero values of 9 could have extremely important implications. It may be possible to create regions with non-zero c^-parameter in heavy ion collisions. If pseudoscalar particles originating in this region manage to escape from the reaction volume, observation of their scalar admixture would provide very convincing proof of the existence of a region of non-zero ^-parameter. In this chapter we have introduced the anomalous effective Lagrangian and discussed some of the interesting implications of this anomalous effective Lagrangian. In the following chapters we expand on some of these implications. Chapter 3 False Vacuum Decay In this chapter we discuss the decay of the lowest metastable vacuum for 0 = 0, described in Chapter 2, in both the zero and high temperature limits . The existence of this 1 specific metastable vacuum is somewhat controversial, as it requires the counterintuitive parameter value, q ^ 1, as we will discuss further below. However, nontrivial vacua have been shown to exist in Yang-Mills theories in the large -/V limit using the AdS/CFT c correspondence and the same phenomenon was observed in the analysis of soft breaking of supersymmetric models[12]. It appears that the existence of nontrivial vacua is a general phenomenon for gauge theories in the strong coupling regime. Therefore this calculation can be considered as one example of a general class of problems. Therefore, in addition to the specific results we obtain for this particular calculation, another important goal of this chapter is to develop techniques which could be useful in the study of this type of phenomenon. The decay rate of the false vacuum in large N Yang Mills theory was estimated in c [31]. The calculation of the decay rate of the false vacuum presented here differs from [31] because it is valid for finite iV and because the heavy glueball degrees of freedom 2 c have been integrated out in our approach, while their calculation derives entirely from gluodynamics. For more comments on the elimination of heavy degrees of freedom see below. The decay rate of a false vacuum was estimated in [25] in the case that 9 w ir and N is finite. Our calculation differs from [25] in the inclusion of the singlet rf field and c 1 T h e results of this chapter are contained in [30] 2 W e use this term to refer to gluonic degrees of freedom which take into account colour confinement. 35 Chapter 3. False Vacuum Decay 36 choice of parameters q = 8 and p = l l i V — 2Nj. As well, both of these estimates only c use the semiclassical approximation, while we determine the effect of the first quantum corrections. The form of the semiclassical decay rate in our case is identical to that of [31, 25]. We make no numerical comparisons as the false vacua involved are all different. These false vacua and domain walls could lead to many interesting consequences in the study of the evolution of the early universe at around the time of the QCD phase transition. One example of this is related to baryogenesis and dark matter, and is described in [15]. The zero temperature decay rate of the false vacua calculated in this chapter is relevant to this particular application. As well, these metastable states hopefully can be experimentally studied at RHIC and the high temperature decay rate calculated in this chapter would be relevant to this research. Bubbles of this false vacuum would display CP odd signatures such as those described in [32] where the large iV limit c was assumed. Even without consideration of these interesting applications, our method of calculation of the determinantal prefactor is useful as an alternative to previous methods. Previous calculations [33] [34] [35] [36] of bubble nucleation rates use a particular method for calculating the determinant ratio of operators of the form: M = -V + 2 UJ + 2 cxV(r), (3.56) involving a theorem from [33]. We prefer to use a more direct approach. Our procedure provides a method for obtaining both analytical approximations and exact numerical calculations for this determinantal prefactor. Our numerical approximation involves only numerical integration, unlike the method of [37] which also involves numerical solution of differential equations and might not be reliable in some cases involving non-smooth perturbation potentials. Our methods are more along the lines of [38] but are sufficiently different to constitute independent results. As well, we use the same methods to calculate Chapter 3. False Vacuum Decay 37 the decay rate in the zero temperature theory which has not been done before. We calculate the decay rate of the lowest metastable vacuum for 8 = 0 in the simplified setting where = 4MeV and the chiral phases are equal. This amounts to studying only radial motion in </>-space. We analyse the problem in the spirit of [39] and only consider transitions between the lowest energy metastable state and the physical vacuum. The results should be easily generalisable to other transitions. This calculation is an approximation for a number of reasons. Most obviously because we make several approximations in order to obtain an analytical answer, but we believe it should give the dominant contribution within a factor of the order of unity. It should be noted that this approximation is in some sense nonperturbative, since it should contain contributions from perturbation diagrams of all orders. The second reason this calculation is an approximation is that we have ignored some complications so that we are able to carry out the calculation. The inclusion of fermions (nucleons) could drastically change the results, but they also drastically increase the difficulty of the calculation, so they are left out in this first approximation. As well, there is the consideration that heavy intrinsic degrees of freedom, such as glueballs, might play a role. It has been suggested[40] that the cusps in the effective potential are due to the integration out of heavy degrees of freedom, that the presence of cusps invalidates the construction of the domain wall from the effective Lagrangian and that a complete calculation must include the heavy degrees of freedom. The decay rate in both these cases can only be decreased by the inclusion of these heavy degrees of freedom. In order to calculate the decay rate, however, we must ignore these effects and content ourselves with calculating an upper bound for the decay rate. The results of this chapter are the only results which depend on specific values for the integer parameters p and q. The very existence of the metastable vacuum state at 8 = 0 requires q ^ 1. The parameters p and q are left in the formulas and so are Chapter 3. False Vacuum Decay 38 otherwise arbitrary, but we must keep this restriction in mind throughout the chapter. When numerical results are given, the values p = 3bQ D = 11AT — 2Nf and q = 8 are C C used. 3.1 The Domain Wall Solution An approximate solution for the field configuration of the domain wall between the metastable and true vacua is needed to calculate the decay rate. For ease of calculation we re-scale and shift the chiral field 4> -> (2/fn<jN~f)<t> ~ ^/(qNf) in order to have the standard normalisation of the kinetic term and a symmetrised form of the potential. The effective potential for 9 = 0 becomes: W(<P) E 1 — cos E 1 — cos Mf(<p) pf* r if <T5 > 0 M/(</>) if 0 < 0 p TT cos JnJNf / — cos QNf 27T \ (3.57) The effective potential (3.57) has a global minimum at 4> = (TTf )/'(2q^Nf) + local minimum at and a n = ~(Trf )/(2q,jN~f), with a cusp singularity between them (see v Figure 3.12). The minima are interpreted as two vacua separated by a high potential barrier (~ G ) which is fairly wide, while the energy splitting, AE, between the states 2 is fairly small in comparison: 2TT AE = m N (tyty) q cos qN f ' + 0K) (3.58) f/ Therefore we can use the thin wall approximation [39] in our calculations. The domain wall solution in this approximation corresponding to the effective potential (3.57) is: <f>d.w.( ) = Pfir x 2qJN f •— + 4arctan < tan f — ] exp\a(x P \ PJ I 4 — xo)l if X < X, 0 Chapter 3. False Vacuum Decay 39 w eff Figure 3.12: Effective potential at 6 = 0 for equal chiral phases 4> = <j) = (f> u = d ^^Z— — — 4 arctan itan ( — | exp — [ u(x — x )] if x > x ,(3.59) 2q^N f [P I 0 \4pJ 0 where x is the position of the centre of the domain wall and: 0 lcPW eff 2q^N~ VE p= f (3.60) PU is the inverse width of the wall. The solution (3.59) is shown as a function of x — x 0 in Figure (3.13). Its first derivative is continuous at x = XQ, but the second derivative exhibits a discontinuity. The energy density of the wall is given by: a = / . \ / < | ^ } ( l - cos ^ + 0(m fl). q (3.61) In what follows we use this domain wall solution in our calculation of the decay rate per unit volume of the false vacuum for both zero and high temperature field theory. The next section describes the semiclassical calculation. Chapter 3. False Vacuum Decay 40 <r> f -R R i X - X, J <!>- • Figure 3.13: The domain wall solution. 3.2 Semiclassical Theory The decay of a false vacuum was discussed in [39] and the formula for the decay rate per unit volume, T/V, was shown to be: T/V = Ae- [l B/h + 0{h)}. (3.62) The semiclassical approximation at zero temperature tells us that B is given by the Euclidean action of <fi: I drd?x B = S 4 (3.63) where U(<f>) is the potential, W ff, for the chiral phases described in the previous section e neglecting the energy difference, AE, between the two vacua: E i - c o s t e ^ - z ) • El-cos if ^ > 0 if0<O. (3.64) Chapter 3. False Vacuum Decay 41 In order for this to be finite we must have lim^oo <p(r) = = —7rf /2qJNf where 7T r — y/x + y + z + t' . The solution of this problem is the four dimensional equivalent 2 2 2 2 to the solution which Coleman calls "the bounce" [39]. It describes a bubble of the true vacuum in the false vacuum at the origin which forms, grows to a maximum size and then shrinks to nothing again (ie. an 0(4) invariant bubble). In the thin wall approximation the bounce solution describing a false vacuum bubble with radius, R, is: (j) for r <C R + Mr) = | (f> . .(R - r) for r « R , (3.65) d w <p_ for r 3> R and the value of B is calculated to obtain[39]: ~2(K)J^j , (3-66) where a is the domain wall energy and: . 3 • 2? • T r y ( l ~ cos ^ 3 4 ^ " J ^ M» ( l _ c o s ^_) _ 27 A 3 - 256 2 * ; j ^ G 2 ) m3|(M>| ' 3 2 ( j is the Euclidean action of the 4D bounce solution. The thin wall approximation is valid because the radius of the 4D bubble, which is found to be R = 3a/AE by minimising the value S^ffo] , is much larger than the width of the domain wall, 1/p. For finite temperature QCD, the semiclassical approximation is slightly different. At sufficiently high temperature the bubble solution becomes a stable 0(3) invariant bubble with radius[41]: R(T) = (3.68) where the temperature dependence of the domain wall energy, S\(T), is not known. In Chapter 3. False Vacuum Decay 42 this case the calculation of B gives[41]: T/V oc e x p ( - S M / T ) -.exp ( - ^ ^ ) • 3 (3-69) It should be noted that this decay rate is for the ground state of the metastable well. The decay rate for excited energy states above the metastable vacuum via thermally activated transitions, while similar in form, is not the same[42]. In the next section we calculate the quantum corrections and determine the decay rate for the zero temperature theory. The section following that presents the same calculations for the high temperature theory. 3.3 Q u a n t u m C o r r e c t i o n s at Zero T e m p e r a t u r e The quantum corrections at zero temperature correspond to the coefficient A in Equation(3.62)[43]: det[-d^ + U"(<f> )} (3.70) b det[-d d» ll + U"(<j>-)] This is the determinantal prefactor mentioned at the beginning of this chapter. The spectrum of the operator in the numerator consists of a discrete spectrum with zero and negative eigenvalues and a continuous positive eigenvalue spectrum starting at LU = 2 U"((/>-). These two parts of the spectrum must be analysed separately and it can be shown that this factor separates into three parts: M /<MI) ) wi 2 d e t ' [ - < V + U"((j> ) (w- ) d e t [ - c V + to b 2 5 2 (3.71) The first term comes from the zero eigenvalues. The second term comes from the negative eigenvalue. The third term is the determinant of the continuous positive eigenvalue spectrum where det' signifies that the zero and negative eigenvalues are to be omitted. With these eigenvalues omitted, the perturbed operator in the numerator has five less Chapter 3. False Vacuum Decay 43 Figure 3.14: Expansion of the partition function. eigenvalues in the spectrum because five of the eigenvalues of the unperturbed operator in the denominator have become part of the discrete spectrum of the perturbed operator. Assuming these eigenvalues have originated from the bottom of the unperturbed continuous spectrum we divide the third term by a factor of co for each omitted eigenvalue. 2 The contributions from the zero and negative eigenvalues are normalised with a factor of co keeping each of the three terms dimensionless. 3.3.1 P o s i t i v e Eigenvalues The contribution of the positive eigenvalues requires evaluation of the determinant ratio: det't-fyd" + U" {<!>„)] (UJ- ) 2 5 However, since the (3.72) + to } 2 in the denominator corresponds to a set of measure zero in (co~ ) 2 det[-d„d» 5 the continuous eigenvalue spectrum, we can omit it and the notation det' which indicates omission of a discrete set of eigenvalues. This determinant ratio is infinite, and to obtain a finite answer we must divide by an infinite factor. The determinant in the numerator can be expanded as: det | - + U"(<j> )\ = exp [Tr log {-d^ b = det| — exp Tr 2 + UJ + V (r)] 2 pert I -UJ X Vyert( ) r -d^d* + CO 2 1 _pert(r)_ 2 d^ N + co , 2 ' 1 +/ + 77 3 Vpert(r) \ [-8^+CO ) 2 (3.73) Chapter 3. False Vacuum Decay 44 It should be noted that the determinant is equal to the partition function of the self interacting massive scalar particle and that the second factor in the last line is expanded up to one loop contributions involving three interactions with the effective external potential (see Figure 3.14). Tracing over a Cartesian basis, we can see that the one and two interaction contributions are divergent but the three interaction contribution is finite: Tr Vpertjr) (3-74) - c V + UJ 2 Tr V ert(r) j -d^ + LO p V (r) Tr where VI _ pert -d^d" + V (k) pert J J4T, Vpert(k) >4 = 2 f d 4 p Vpert(-k) (k + co ) ((k + p) 2 P d 4 q 2 Vpertjk) 2 (3.75) + UJ ) 2 V (-k V (p) pert pert - p) 43.76) (k + UJ ) ( p 2 + UJ2) ((k+p + q) + uj f 2 2 LO 2 2 2 is the Fourier transform of the V (r). We only get a finite answer if we pert divide through by the infinite factors. The actual calculation is most easily done using hyperspherical coordinates in four dimensions, (r, 6, cf>, ip), and expanding the eigenfunctions x( , <f>, ip) in terms of the 4D r hyperspherical harmonics: oo n \l\ I \ (3.77) n=0 l=-nm=-\l\ which are discussed in Appendix B. In this situation, the Laplacian, when operating on each term of (3.77) with quantum number, n, becomes: An + 8n-l-3' 4r 2 ^ - > 1372 dr 2 r 3/2 = J_T> 3 / 2 r 2 (3.78) There are 2n(n + 2) + 1 such terms corresponding to different values of I and m but with the same eigenvalue. Therefore: d e t [ - ^ ^ + U"} = (det[-2? + n u"]) ^ . 2n +1 (3.79) Chapter 3. False Vacuum Decay 45 First consider the denominator of (3.72), det[—V + UJ ]. In order to calculate this 2 n determinant we need to solve the eigenvalue equation: -co [D n 2 + X]u{r) = 0. (3.80) The solutions to this differential equation that are well behaved at r = 0 and r = oo are: u(r) where UJ = X y/X — LO . j 2 n = {uj rf' j {uj r) n (3.81) = ^J^-J {to r), 2 x x n+1 x are the 4D analogs of the spherical Bessel functions and J i(x) n+ are Bessel functions of the first kind. These solutions become purely oscillatory as r —> oo. Notice that this solution is only well defined for X > uo and indeed there are no solutions 2 for values of A < to which are well behaved at r = 0 and r = oo. Therefore the continuous 2 spectrum of eigenvalues can be written as A = UJ + co and: 2 2 det[-V + co ] = Y[{UJ + uo ). 2 2 n The numerator of (3.72) involves the "potential" (3.82) 2 x For the symmetrised ef- U"(fa(r)). fective potential (3.64), the "potential" for this problem is approximately constant except for a small perturbation in a small region near the bubble wall (see Figure 3.15). The constant potential has been analysed above. The solutions with or without the perturbations are identical for r <C R and almost identical for r 3> R. The perturbed solution differs from the unperturbed solution in this latter region at most by normalisation and a phase shift, cd r = uj r + 5(uj ). x x x In this situation, for each value of n, we can obtain the ratio of determinants for a discrete spectrum by [44]: n 4±4 ^^UJ Z + UJ{ = exp (E In ^ ± 4 ) y— * exp uo rUJ ) lJ 1 l x fc y— 1 ( f* " uo + uj l l x *') • J (3.83) which becomes in the continuum: exp 1 -I /•«> TX JO 2LO 5{UJ ) x UJ 1 X + uoi (3.84) Chapter 3. False Vacuum Decay 46 204000 203000 0 R Figure 3.15: The "potential" U"(<f> (r)). b This formula gives the determinant ratio for a particular value of n if we know the phase shifts, 8{LO\). In order to calculate the phase shifts, however, we must make further approximations. Consider both the perturbed and the unperturbed equations: An + 8n d 2 2 dr 2 d 2 An 2 AT + 8n + 3 dr 2 +3 + (A - UJ ) 2 2 Ar 2 (A - UJ + 2 u{r) = 0, (3.85) 0. (3.86) V {r)) v(r) = pert Multiply (3.85) by v(r) and (3.86) by u(r), subtract, and integrate over r from 0 to oo to obtain: roo / (u"(r)v(r) Jo - v"(r)u(r)) dr = (u'(r)v(r) - v'(r)u(r)) poo |g° = / Jo V (r)u(r)v(r)dr. pert (3.87) Using the fact that both solutions vanish at r = 0 and using the asymptotic form of the Chapter 3. False Vacuum Decay 47 Bessel functions we obtain an exact formula for the phase shift, sin5(ujx) - 1 r°° — / LU Jo S(u)\): (3.88) V (r)u(r)v(r)dr. pert X Since the perturbed differential equation is extremely difficult to solve, we use perturbation theory to obtain: 8{u ) x ~ 1 6{u ) x 0 = — / f JO ld\ 7T 0 0 V {r)u(r) dr = - 2 pert / /*°° (3.89) V (r)rJ (uu r) dr, pert Z JO n+l 2 x for small phase shifts. Using this result for the phase shift in Equation (3.84), we obtain the ratio of determinants: det[-2> + £/"(&)] det[-X\ + LU ' n 2 roo exp Jo re dto\dr JO u) r (3.90) V ert{r)J +l{LO rf X p UU +LU 2 n X 2 for each value of n. The complete determinant ratio then becomes: d e t f - d ^ + LU ' 2 exp £ ( 2 n ( n + 2) + l) n c °° dLU dr rCO n=0 X " o , U x r <--^ry 2 Pert{r)Jn+l{LU 2 V X uc ~r LUX LU + JO 0 Z x (3.91) for an arbitrary perturbation potential. The next step in our analysis is to take into account the specific potential in our problem. It should be noted that this approximation, (3.91), gives the same result as the exact answer expanded to one loop (see Eq.3.73 and 3.74) where the trace is taken instead over \rnlm), \knlm) bases. This formula is suitable for numerical evaluation, but in order to obtain an analytical answer we must approximate the perturbation in the potential U" by: -0 R 0 2 ii 1 0 > r > n otherwise , (3.92) Chapter 3. False Vacuum Decay 48 Figure 3.16: Comparison of the perturbation V (r) and the approximate perturbation. The perturbation V (r) is the dotted line and the approximate perturbation V (r) is the solid line. pert pert pert Chapter 3. False Vacuum Decay 49 0.03 0.01 0.025 0.008 0.02 0.015 >0 0 . 0 0 6 0.01 0 . 004 0.005 0.002 400 800 400 1200 800 1200 CO, 1200 CO, n=5 n=0 0 . 006 0 .005 0 . 006 0.004 5 0 >0 0.003 0.004 0 . 002 0 . 002 0.001 400 800 co, 1200 400 800 n=15 n=10 Figure 3.17: The phase shift 8 (u ) Notice the changing vertical scale. Q where B= (£/"(</>(0)) — shown as a function of to for various values of n. x x Figure 3.16). With this approximate potential we U"(4>(R)))(see find the phase shift via (3.89): 8 (UJ ) = 0 X -200-^ \^r R 2 {Jn i(ruj f 2 + - x J (ruj )J (ru )) n -4-2n^4 r (ruj ) |2+2n , / 4 x 4)!(„ + T ) T | (n + 2 ( n + 2)(n 4) + -r Lj\ 11 2- - r R(ru ) 2+2n 2 2n 3 x 10 (n + f)[(n + l)!] 2 ( 3 X n+2 x 3 (n + F 2 x \ n + 4, 2n + 3; -rW ) x (n + |, n + 5, 2n + 4; - r ^ ) 2 2 5 7 \ n + -, n + -, n + 2, n + -, 2n + 3; -r u J Using this result we can verify that 2 5(to ) x <C J -i 11/lOfl . (3.93) 2 X 1 which means that the approximations of (3.89) are justified. The phase shift as a function of u> is shown for a few values of x n in Figure (3.17). However, (3.93) is very difficult to work with so we make another approximation. For co large enough, we can use the asymptotic approximation for the x Chapter 3. False Vacuum Decay 50 Bessel function in (3.89), J +i{x) = \ — \p i(x) COS X n n+ V 7TX I -{ D n+ £]-Gn i(*)sin [ x - ( n + | ) | ] } , + (3.94) where: 1 _ (p-l)(p-9) ( M - 1 ) ( M - 9)(/i - 25) (ja - + 2! (8xy 4! (8x) n ^ - k^ ) (M-1)0"-9)(//-25) , Vrc+lW — „1!, /(8a;) „ x „ . (8xY / „ x , 13! /i = 4(n + l ) . 49) 4 , (3.95) 1 (3.96) (3.97) 2 In this way we find: <*oM This estimate for 307ru; TT A ui\ (3.98) x > has an infrared cutoff in u\, below which it is not a good 8 (LJ\) 0 estimate, and this will be true no matter how many terms in the approximation of the Bessel function we keep. Therefore, as a first approximation we keep only the first term and cut off the approximation for co\ 5 (id\) 0 integral at the position of the first peak in 5 (id\). 0 Using this we obtain a first estimate of the ratio of determinants for each value of the quantum number n given by: det[-2? + U"{<j> )] IRQ exp det[-V + £/"(</>_)] 30™ arctan w h id (3.99) n where: u)(n) «33.16n + 72 = an + b, is the infrared cutoff given by the location of the first peak in (3.100) S (id ). 0 X Using the result (3.99) we can calculate a first approximation for the complete determinant ratio. For large values of n the terms of the sum in the exponent of (3.91) approach: P 4R 30TTU; 2 n \an) (3.101) Chapter 3. False Vacuum Decay 51 and therefore the sum is infinite. The first estimate for the complete determinant ratio based on (3.91) is infinite. This was to be expected based on comments at the beginning of this section and, since this approximation exactly coincides with (3.74), we subtract this term from the exponent of (3.91) as our normalisation prescription for the complete determinant ratio. Note that we have used the approximation of (3.98) to calculate the infinite factor, but the equality with (3.74) holds before this approximation was made. In order to find a finite result for the determinant ratio we adjust our earlier approximation for the phase shift: 1 = arcsin(<5 (u; )) = 6(LO ) 0 X A S (LO ) 0 X + 3 -5 (WA) 0 3 + — S (LO ) 5 0 X + (3.102) which would add correction terms to the exponent in (3.99). The second term in (3.102) leads to a factor of the form: TT 2 exp 24/0 f°° \3 LO 5 (LO ) , d U _ 3 X X - co 2 0 + (3.103) X LO 2 X for each value of n. This term is also divergent but does not exactly coincide with the next term in the expansion (3.73). However, the divergent contribution in each must be the same. Therefore dividing through by the previously obtained infinite factor (ie. the right hand side of (3.99)) leads to the renormalised value: TT {2RP\ f LO det[-V +U"((f> ) LO — — — — <, , \ — arctan exp det[-V + [/"(</>_)] 24 V 30vro; LO(n) M ). 2 n 3 b n n (3.104) For large values of n the terms of the sum in (3.91) approach: (3.105) which also leads to a divergent sum and therefore an infinite value for the determinant ratio. The divergence here must be contained in the divergence in the two interaction Chapter 3. False Vacuum Decay 52 term of the expansion (3.73), as are other logarithmic divergences obtained from keeping more terms in the approximation of the Bessel function (3.94). However, subtracting the two interaction term will almost certainly leave a finite contribution at the next order in n. We assume this is the case, but since we are doing an approximate calculation, we will not calculate the finite contribution from the two interaction term or any correction terms to our approximation, since they will have the same physical dependence as the approximation we will give. We should further note that correction terms are probably not calculable analytically and are not likely to be a problem for our results. They are obtained from integrating oscillatory functions (see second term in (3.98)) over u>\ in (3.103) and should not contribute very much. The renormalised approximation to the complete determinant ratio (3.91) that we obtain from (3.104) is given by: 7T / 2R8 \ 24 3 ^ J 2 6 X P X P r £ ( 2 n ( n + 2) + l) Ln=0 TT (2RR X 24 [SOTTUJ J 2 6 3 LO ~£2n 71=1 100 y"(2n(n + 2) + 1) < 3 LO 4 w 2 ,3 LO 3a n 3 3 3.106) LO - arctan 3 4W , . . LO" - j N ( 2 n ( n + 2) + l) 71=1 arctan Lo(n) 100 LO LO" ^ (3.107) I where £ is the Riemann Zeta function: oo -i (3.108) 71=1 Evaluating gives: TT [2RPY exp[-2 x 10~ . 2 exp -5300 5 24 \30TVLOJ (3.109) Because of the small value of the exponent we can see that the exponential is extremely well approximated by: exp -53007T 24 2 (2Rp\ ' 3 \ 307TLOJ _ 53007T 24 2 (2RpV [SOITLOJ (3.110) Chapter 3. False Vacuum Decay 53 which is our result for the determinant ratio (3.72). This is a nonperturbative calculation because Eq.(3.84) is a nonperturbative resummation of the perturbation expansion (3.73) of the determinant. While we make an approximation through Eq.(3.102), this expansion is nonperturbative since the terms of the expansions do not coincide. Our approximate calculation should therefore have contributions from all perturbative diagrams. This nonperturbative resummation of diagrams in the determinantal prefactor is very similar in spirit to that of [45]. This result constitutes the contribution of the positive eigenvalues to the quantum corrections to the decay rate in the zero temperature theory. In the next section we consider the zero and negative eigenvalues. 3.3.2 Zero a n d N e g a t i v e Eigenvalues The zero eigenvalues contribute ^JB/2TV per collective coordinate[39]. The action of the 4D bubble is independent of the centre of the bubble which means there are 4 collective coordinates leading to the first factor in Eq.(3.71). The eigenfunctions of zero eigenvalue are: d Xo>(x, V, z, t) = ^Mr) dv = d x^ d Mr) = Tr T ~^Mr), (3.111) where r = y/x + y + z + t and: 2 2 2 2 t= r cos(V>) ~ Y , (3.112) z= r cos(0) sin(V>) ~ Ino, (3.113) 100 y = r cos(0) sin(0) sin(^) ~ Y 1U x= + Y„ n u r sin(0) sin(0) sin(t/>) ~ Y - r _ i . in n (3.114) (3.115) Therefore these eigenfunctions all correspond to n = 1 and since there are no radial nodes we can be sure that there are no negative eigenvalues with n ^ O : Chapter 3. False Vacuum Decay 54 In [46] Coleman argued that there is only a single negative eigenvalue for an 0(4) invariant bounce. We assume this and determine its value. We use the method of 3 Coleman in [39] with a slight modification. As was argued by Coleman, the only possible eigenfunctions of negative eigenvalue are those that are bound to the bubble wall. For such eigenfunctions we can approximate the centrifugal potential in (3.78) by a constant determined by its value at the bubble wall (r = R): An 2 ^p ~r~ ^pn + 8n + 3 (3.116) AR 2 where X is a number independent of n. We know that for n = 1 the lowest eigenvalue is p zero: A01 An + 15 (3.117) 0. AR 2 Therefore we can obtain the lowest eigenvalue for n = 0: A_ = A,oo A + 0 3 15 AR AR? 2 3 + _3_ AR 2 (3.118) R' 2 This value is different from the value obtained by Coleman, but only by a factor of 2. We cannot explain this discrepancy but can only stand by our calculation. 3.3.3 Decay R a t e for Zero T e m p e r a t u r e Therefore combining the zero temperature semiclassical result (3.66) with the quantum corrections obtained in the last section, we estimate the decay rate per unit volume of the false vacuum in the zero temperature theory to be: T/V = det'i-d^ + U"(<f> )] LU exp (-S ) ( u; ( ^ 4 {LU- ) 2 5 / 2 7 7 r V \ (27ir o Lu \ 2 A 2 2 RLU det[-0„d" + 4 5300TT 2 1 1.55 x 10" MeV * 3 x l O ^ f n T V . 18 (3.119) b 4 1 24 LU } 2 2RpY 307m; j -1/2 (3.120) (3.121) This is most likely a good assumption but not proven due to the cusp in our potential. For more details see [43],[46] and [47]. 3 Chapter 3. False Vacuum Decay 55 It should be noted that the quantum corrections are negligible so the decay rate is basically determined by the semiclassical result. We believe that this observation would be unchanged by the inclusion of corrections that we have ignored in this calculation. While the quantum corrections did not turn out be significant in this case, the fact that they are not is relevant to the baryogenesis mechanism mentioned at the start of this chapter. As well, the techniques applied to the problem may prove useful in other calculations of this type. In the next section we perform the same calculation in the high temperature limit. 3.4 Q u a n t u m C o r r e c t i o n s for H i g h T e m p e r a t u r e The quantum corrections to the decay rate at high temperature are[41]: 1/2 det'[-$# + £/"(&)] (3.122) which again factors into three parts corresponding to zero, negative and positive eigenvalues respectively. 3.4.1 P o s i t i v e Eigenvalues The calculation in three dimensions is extremely similar. The expansion of the determinant ratio is exactly the same as in 4D (3.73). Tracing over a Cartesian basis, we can see that the one loop term is divergent while the two loop term is finite: (3.123) Tr (k + UJ ) ((k + p) + UJ )' 2 2 2 2 (3.124) Chapter 3. False Vacuum Decay 56 The calculation is done using spherical coordinates and the eigenfunctions are expanded in terms of the spherical harmonics: 0 0 7 1 x(r,M) = £ 1 11 1 £ E n=0 l = - n G (r) (3.125) ^ Y U M ) . m m=-\l\ T The Laplacian, when operating on each term of (3.125) with quantum number I, becomes: 1(1 + 1) d 2 d d -> r dr l l 1 (3.126) r = -T>ir. r 2 2 r Therefore: + U"} = TJ(det[-2?i + U"}) . 1=0 deti-did 1 (3.127) 2l+l The solutions to this differential equation that are well behaved at r = 0 and r = 0 0 are: u(r) = u rji(uj r) x ^^J {oo r), = x l+l/2 (3.128) x where j are the usual 3D spherical Bessel functions and LO = y/\ — LO . We obtain: 2 t X S(LO ) V {r)u{r) dr W— / X L0 X =- / 2 pert JO V (r)rJ (co r) dr, (3.129) 2 pert l+1/2 x Z JO thus giving : 45 \det[-d d + U"(<l> )]\ ROO RE D\r = e x p ^ ( 2 / + 1) roo /-oo ;V (r)J (oj rY * •<•> ^ deti-did + a* Jo Jo duoxdrLO +L0l ^ 1=0 i i b A 1 pert l+1/2 x 2 2 (3.130) As in the zero temperature case we can find the exact expression for the phase shift with the approximate potential but it is too difficult to work with. Instead we approximate the Bessel functions as in (3.94) where (n + 1) —> (I + 1/2). In this way we find: 5(LO ) X W S (UX) 0 2Rf3 + 30TTLO 0(l/u ). 2 (3.131) X Note that we have again dropped the factor of (to ) and the "det"' notation, as the omitted eigenvalues correspond to a set of measure zero. The summand in this expression appears in [33], but is only used for the high energy modes. The non-divergent term in what follows does not appear in [33]. 4 5 2 3 Chapter 3. False Vacuum Decay 57 A first estimate of the ratio of determinants for each value of the quantum number, 1, is given by: det[-V + U"(fa)} 2R(3 t exp det[-2?, + U"{4>-)} where 0(1) LO arctan 307TCJ (3.132) 4 9 . 5 / + 8 2 = cl + d is the infrared cutoff given by the location of the first peak in 5 (LO\) as a function of LO\. For large values of I the terms of the sum in 0 (3.130) approach: , 307ro; (%) . \cL M (3.133) and the sum is infinite and therefore the complete determinant ratio is also infinite. Again, this was to be expected and, since this approximation exactly coincides with the one loop term, our renormalisation prescription is to remove the one loop term from the exponent of (3.130). Again we adjust our earlier approximation for the phase shift: 1 S(u ) = arcsin(<5 (a; )) = ^o(w ) + 0 x A A and the correction is implemented as in 3 + -S (LO ) 3 0 X —S (LO ) 5 0 + X • (3.134) Evaluating the contribution to the (3.103). determinant ratio leads to: det[-Vi + U"((f> )] exp det[-V + f/"(0_)] TT2 b 3 (2RP" For large values of I the terms of the sum in 2 (0 LO 2A I 30TTO; t TT LO (3.130) ( 2Rp X — arctan ui MO, (3.135) approach: LO (3.136) The sum is finite so we do not need to further normalise and the result for the complete determinant ratio (3.130) is: Tr (2RP\ \ r ° 2 exp 2 4 Zrl0 r 0 0 LO \30TTU 1 0 0 f , ,3 1 O, ,3 ( LO W). , ,3 (3.137) Chapter 3. False Vacuum Decay 58 Evaluating gives: TT 2 exp 82— 24 (2Rfi\ J exp[3 x 10~ ;. (3.138) 7 Because of the small value of the exponent, we can see that the exponential is extremely well approximated by: 82TT (2R(3Y 24 \307ruj) (2Rj3\ 2 827T : exp ~2A V307ra;) 2 (3.139) This result constitutes the contribution of the positive eigenvalues to the quantum corrections to the decay rate in the high temperature theory. In the next section we consider the zero and negative eigenvalues. 3.4.2 Zero a n d Negative Eigenvalues In the 0(3) invariant bubble there are 3 collective coordinates leading to the first term in (3.122). The eigenfunctions of zero eigenvalue are: d dr d dx^ dr Mr) x» d r dr M )> r (3.140) where r = y/x + y + z and: 2 2 2 z = r cos(0) ~ Yin, (3.141) y = r cos(</>) sin(0) ~ Y + (3.142) x = r sin(^) sin(0) ~ Y — l i - i - (3.143) n u Therefore these eigenfunctions all correspond to I = 1 and since there are no radial nodes we can be sure that there are no negative eigenvalues with I ^ 0. We assume, as in the previous section, that there is only a single negative eigenvalue and concentrate on obtaining an approximation. The only possible eigenfunctions of Chapter 3. False Vacuum Decay 59 negative eigenvalue are those that are bound to the bubble wall. For such eigenfunctions we can approximate the centrifugal potential in (3.126) by a constant: A„ = A P 1(1 + 1) + p ( 3 . 1 4 4 ) where A is a number independent of /. We know that for / = 1 the lowest eigenvalue is p zero: Aoi = A + = 0. 0 (3.145) Therefore we can obtain the lowest eigenvalue for I = 0: A_ = A = A = - - ^ . 00 3.4.3 (3.146) 0 Decay Rate for High Temperature Therefore combining the high temperature semiclassical result (3.69) with the quantum corrections obtained in the last section we estimate the decay rate per unit volume of the false vacuum in the high temperature theory to be: Rw R / V = ^[ U 3 /l67rg (T) y ^ 1 [HAEW) 3 / / 167^ (T) \ ) rw ] 3 exp 7 827T 2 4 2 (2RB\ 31 "1/2 U O T T U J (3.147) We cannot obtain a numerical estimate for the decay rate at this time since the temperature dependence of the semiclassical result is not currently known. This result could be important to the study of false vacuum states of this type in heavy ion collisions. 3.5 Summary In this chapter we nonperturbatively calculated the decay rate of a particular metastable vacuum state for 9 = 0 predicted by the form of the effective potential for the light chiral fields. We calculated the decay rate in both the zero and high temperature theories. The Chapter 3. False Vacuum Decay 60 low temperature result is significant to the baryogenesis mechanism described in [15]. The high temperature results could be important for testing the theory in heavy ion collisions. Both results could have important implications for the study of the development of the early universe near the time of the Q C D phase transition. In addition the method used in these calculations represents a new approach that could be useful in other calculations of this type. Chapter 4 Production of Nontrivial Theta Vacua In this chapter we show in a simplified numerical model that nontrivial 0-vacua can be realized . This nontrivial vacuum is very different from the metastable vacuum discussed 1 in the previous chapter, in which we had 9 = 0 throughout space and the metastable vacuum was definitely metastable. The nontrivial (9-vacuum that we discuss in this chapter has 9 ^ 0 by definition and is the true vacuum state in a universe with 9 ^ 0 everywhere. The energy of this vacuum state is given by Equation (2.37). As we have mentioned previously, we live in a universe where it is known that 9 < 10~ [3] so the 10 true vacuum state of our universe corresponds very closely to the 9 = 0 vacuum state. However, it is possible that nontrivial 9-vacua can be produced inside a high temperature fireball created in heavy ion collisions. While this is a simplified model which differs somewhat from the real physics associated with high energy ion collisions, there is no reason to believe that the nontrivial 9-vacua cannot be realized at RHIC. Our results also give a rough estimate of the time it takes for these nontrivial 0-vacua to be formed. In order to have a hope of observing them, they must form within the time that the central region of the fireball is isolated from the true vacuum. The real importance of observing effects of non-zero 9 vacua in relativistic heavy ion collisions is that it provides us with an opportunity to study the physics of the early universe at around the time of the QCD phase transition. In the very early universe the conditions would have been similar to the inside of the heavy ion collisions and the 1 T h e results of this chapter were p u b l i s h e d i n [48] 61 Chapter 4. Production of Nontrivial Theta Vacua 62 existence of a ^-parameter could have had significant implications to the evolution of the early universe. Note that the formulas in this chapter are valid for arbitrary values of the parameters p and q. The numerical results were obtained using the specific values p = N and q = 1 c In the next section we will discuss the Disoriented Chiral Condensate (DCC) which is closely related to the 9-vacuum state and is fairly well understood. In the section after that we will show how to extend the idea of the DCC to obtain the 0-vacuum state. 4.1 D C C The creation of a 9-vacuum state is very similar to the creation of the Disoriented Chiral Condensate (DCC) in heavy ion collisions, which has been a subject of interest for some time [49, 50, 51] (among many others, see also a nice review [52] for a discussion of DCC as an example of an out of equilibrium phase transition). Therefore, we discuss the DCC first and then develop the idea of a 0-vacuum state afterwards. DCC refers to regions of space in which the chiral condensate points in a different direction from that of the ground state. The DCC could exist inside a hot shell of debris produced in heavy ion collisions and be protected from the true vacuum by this shell. For both the DCC and the 0-vacuum state, the difference in energy between the created state and the true vacuum state is proportional to the small parameter m and negligible at high temperature. q The chiral condensate, (qq), is the vacuum expectation value of the composite operator qq made up of quark(g) and antiquark(g) fields. Mesons are bound states of quarks and anti-quarks, and so the chiral condensates are just vacuum expectation values of the meson fields. The low energy dynamics of pions was modelled in [50] by a linear sigma model involving the four-component field (a,if). The potential of the four component field is a steep 4D Mexican hat potential with a slight extra slope (~ m ) in the a q Chapter 4. Production of Nontrivial Theta Vacua 6 3 direction. In the region exterior to the hot shell of debris created in a heavy ion collision, the fields have the vacuum values (v,0). In the interior region, which has recently been in contact with the extremely high energy decay products, the pion fields can initially take non-zero values and the (CT, 7?) vector can point in an arbitrary direction in the four dimensional configuration space. In particular the vacuum expectation value of the fields differs from the vacuum value: (v,0). Non-zero vacuum expectation values are defined as condensates and the four vector of field condensates differs from the vacuum direction or vacuum orientation. This is what is meant by a Disoriented Chiral Condensate. The condensates of course eventually relax to the vacuum values. The absolute value of the chiral condensate right after the phase transition is expected to be close to its final (zero temperature) magnitude because of the steepness of the Mexican hat potential. However, the vacuum orientation of the formed condensate takes a longer time to relax due to the small free energy difference ~ m between the DCC q state and the vacuum state. The high energy products of the collision expand outwards at relativistic speeds. Inside this shell is a region of DCC isolated from the true vacuum and at effectively zero temperature. If the cooling process of the interior is very rapid the system is initially out of equilibrium, but its evolution can be determined using the zero temperature equations of motion. It has been shown[50] that in the quenched approximation, where the temperature goes from ~ 200 MeV to zero instantaneously, the subsequent evolution exhibits a temporary growth of long wavelength spatial modes of the pion field corresponding to domains where the chiral condensate is approximately correlated, but oriented differently from the vacuum value. It is useful to discuss the DCC in an effective Lagrangian approach to emphasise the similarity with a nontrivial 0-vacuum state. Due to the fact that the absolute value of the chiral condensate right after the phase transition is expected to be close to the Chapter 4. Production of Nontrivial Theta Vacua 64 zero temperature magnitude, we can parameterise Goldstone fields by the same unitary matrix, U, introduced in Chapter 2: U = e^ \ f (¥ V ) L = - | ( * A > | U^. R (4.148) The energy density of the DCC is determined by the mass term in an effective Lagrangian: = -\Tr{MU where we put m u = = m + M^U ) = -2m|(tftf)| ] cos(^), for simplicity. Eq.(4.149) implies that any not a stable vacuum state because -g^l^o (4.149) (f> ^ 0(mod 2ir) is 0; d therefore the vacuum is misaligned. a n On the other hand, the energy difference between the misaligned state and true vacuum (<f> = 0) is small and proportional to m . Therefore, the probability to create a state q with an arbitrary <j> at high temperature T ~ T c is proportional to exp[—V(E$ — Eo)/T], where V is 3D volume, and depends on (j) only very weakly, i.e. (p is a quasi-flat direction. Right after the phase transition when (tftf) becomes non-zero, the pion field begins to roll down the potential slope toward <> / = 0, and of course overshoots it. Thereafter, the 4> field oscillates. One should expect coherent oscillations of the ix meson field which would correspond to a zero-momentum condensate of pions. This is exactly what was found in[50]. Eventually these classical oscillations produce real TT mesons which hopefully can be observed. 4.2 Nontrivial 0-vacua Now we turn to the nontrivial 0-vacuum state when the 9 term is also non-zero and therefore the 0-vacuum state can be formed. In the evolution of the early universe this non-zero ^-parameter could be a fundamental 9 associated with a dynamical axion. In the case of heavy ion collisions, however, it is much more likely that we might observe what is called an "induced" ^-parameter. Chapter 4. Production of Nontrivial Theta Vacua 65 The induced ^-parameter could arise if a region is created in a heavy ion collision over which the chiral phases are significantly correlated about a non-zero value. If this occurs then we perform a chiral rotation of the type discussed in Chapter 2: U -> e-^f^U (4.150) = e, i0 where —fais the non-zero value about which the chiral phases are correlated and we obtain a new theory with an effective value of 9 induced by the initial conditions. The fields behave as if the ^-parameter were non-zero. As we said above, we are far more likely to observe an induced 9 vacuum in heavy ion collisions, so we have this type of 9 parameter in mind throughout this chapter. The production of nontrivial 9-vacua occurs in much the same way as for the DCC discussed above. The new element is that, in addition to chiral fields differing from their true vacuum values, the ^-parameter of QCD, which is zero in the real world, becomes effectively non-vanishing in a large region of space. The starting point is again the effective potential for the chiral fields (2.31). The energy density of the misaligned vacuum is given by the potential: V(fa,9) = -Ecos ,5>-* ]+-i / = 0,l,...,p-l, (4.151) for (21 - 1 ) - < 9 - Y] fa < (21 + 1)- . q (4.152) q As discussed previously in Chapter 3, and stress again here, the crucial point is that the ^-parameter appears only in the combination, J2fa - 9. The most important difference between Eqs. (4.149) and (4.151) is the presence of the parametrically large term ~ E ^> m \(tyty)\ q in the expression for energy (4.151). This term does not go away in the chiral limit and provides a non-zero mass for the rf meson which is expressed Chapter 4. Production of Nontrivial Theta Vacua 66 in terms of the parameter, E. It was for exactly this reason that it was thought until recently[53] that the nontrivial 0-vacua would involve too large an energy cost to be produced because of the large parameter associated with the 0-parameter. The key point is that, for arbitrary phases fa, the energy of a misaligned state with non-zero 9 differs by a huge amount (~ E) from the vacuum energy. However, when the relevant combination fa —9) from Eq.(4.151) is close by an amount of the order 0(m ) g to its vacuum value, a Boltzmann suppression due to the term proportional to E is absent, and an arbitrary misaligned 0-state can be formed. Indeed, in the limit Mj <C E, because of the large parameter in the first term, the vacuum state that is favoured for non-zero 9 is the solution fa « 9/Nf. Substituting this back into the potential we reproduce the well-known result[10] for the vacuum energy as a function of 9: V{9) xi-E-YtMi cos(9/N ). (4.153) f This shows that the energy cost of creating a nontrivial 9-vacuum is proportional to the much smaller parameter, Mi oc m , as in the case of the DCC. This is a very important g result because it means that it may be possible to create a nontrivial (9-vacuum state in heavy ion collisions. At this point we can apply the same philosophy as for the DCC. The chiral fields , 2 fa, are allowed to take random values and after the phase transition roll toward the true solution fa « 9/Nf and of course overshoot it. The situation is very similar to what was described for the DCC with the only difference that, in general, we expect an arbitrary ^-disoriented state to be created in heavy ion collisions, not necessarily the \9 = 0) state. The difference in energy between these states is proportional to m from (4.153). q Now we study the evolution of the chiral fields with time. For this study we assume A s was discussed in Chapter 2, for 6 ^ 0 the meson fields are not pure pseudoscalar fields, but rather have a scalar component as well; the mixing angle between the singlet and octet combinations also depends on 9. 2 Chapter 4. Production of Nontrivial Theta Vacua 67 the following situation. The rapid expansion of the high energy shell leaves behind an effectively zero temperature region in the interior which is isolated from the true vacuum. The high temperature non-equilibrium evolution is very suddenly stopped, or "quenched", leaving the interior region in a non-equilibrium initial state that then begins to evolve according to (almost) zero temperature Lagrangian dynamics. Starting from an initial non-equilibrium state we can study the behaviour of the chiral fields using the zero temperature equations of motion. The equations of motion are non-linear and cannot be solved analytically but we can solve them numerically in order to determine the behaviour of the fields. The production of a nontrivial ^-vacuum is indicated by the fact that the chiralfieldsrelax to constant and equal non-zero values on a time scale over which spatial oscillations of the fields vanish. This indicates the presence of a condensate parameterised by the ^-parameter. The formation of a non-perturbative condensate is also supported by observation of the phenomenon of coarsening (see below) and by a test of volume-independence of our results. 4.3 Numerical Evolution The equations of motion for the phases of the chiral condensate with two quark flavours consist of two coupled, second-order, nonlinear partial differential equations: fa- where V 2 V fa 2 + yfa + ^-V(fa,e) =0, i = l,2 (4.154) is the three dimensional spatial Laplacian and the potential is given in (4.151). Emission of pions and expansion of the domain contribute to the damping, 7 , as might other processes. We do not know exactly how they would contribute, but we simulate these unknown effects by including a damping term with a natural value for the damping constant, 7 ~ AQ D C ~ 200 MeV. Chapter 4. Production of Nontrivial Theta Vacua 68 In this section we present the results of a numerical solution of the equations of motion. The initial data for each of the chiral fields fa was chosen on a 3D grid of 16 3 points. The initial data consisted of random values of fa and fa = 0. The initial data was evolved in time steps using a Two-Step Adams-Bashforth-Moulton Predictor-Corrector method[54] for each grid point with the spatial Laplacian approximated at each grid point using a finite difference method. We used periodic spatial boundary conditions. The grid spacing was determined by the length of a side of the spatial grid which was varied in order to vary the volume. The size of the time step between successive spatial grids was much smaller than the spatial grid spacing and was fixed at about l0~ MeV~ . b l We evolved the data for 8000 time steps and then applied a Fast Fourier Transform [55] to the spatial data at evenly spaced time steps. We then binned the data in small increments of the magnitude of the wave vector and averaged over the values in each bin in order to obtain the angular averaged power spectrum. This procedure was carried out for different volumes. In all cases the results were qualitatively the same. We saw an initial growth of long wavelength modes as in [50] and subsequent damped oscillation of all modes. The \k\ = 0 modes oscillate and approach the non-zero equilibrium values of the fields. They exhibit this behaviour in the same time frame in which the Fourier coefficients of the modes with non-zero wave vectors fall to a tiny fraction of the zero mode coefficient. This qualitative behaviour occurs for different total volumes and grid sizes suggesting that this behaviour is not due to finite size effects. The Fourier zero mode corresponds to a spatially constant field value. Therefore the fact that the zero-mode goes to non-zero constant value means that the field approaches a classical spatially constant field configuration. This clearly corresponds to a nonvanishing vacuum expectation value of the field which is the definition of a condensate. Therefore we can say that a condensate has been formed. This state is the ground state Chapter 4. Production of Nontrivial Theta Vacua 69 Time (fm/c) Figure 4.18: Time evolution of \fa\ shown for various \k\. Notice that the zero mode settles down to fa 9/Nf. within the region where 0 / 0 . We refer to this condensate as the nontrivial 0-vacua. We should note that the \k\ = 0 mode is really only a quasi-zero mode as it is obtained in a finite spatial volume with periodic boundary conditions. However, the quasi-zero mode approaches the same value irrespective of the total spatial volume indicating that this really is a condensate. If it were not, we would expect the value of the coefficient to decrease when the volume of the system increases. The evolution of the Fourier modes of the <f>fieldsis shown in Figure 4.18 for the specific case of 9 = 27r/16 and a spatial grid of 10 fm on a side. The initial values of the <t>fieldswere randomly chosen within the range —77r/16 to 77r/16. The zero mode clearly settles down to a non-zero value. All higher momentum modes vanish extremely rapidly and are negligible long before the zero mode settles down to its equilibrium value. The instantaneous distribution of Fourier modes for the evolution above is shown in Figure 4.19 at a few different times. This graph clearly shows the amplification of the zero mode as time increases. This phenomenon of coarsening and the formation of a Chapter 4. Production of Nontrivial Theta Vacua 70 0.15 0.1 t =2t 2 0.05 t =1.6x10- MeV" 3 1 o 10 5 1 15 20 25 Magnitude of k (arbitrary units) Figure 4.19: The Fourier distribution of the field is shown at three times within the first 1000 time steps of the evolution. The amplification of the zero mode demonstrates that the system exhibits the coarsening phenomenon. nonperturbative condensate is very similar to earlier discussions in [52]. In Figure 4.20, \(j)k \ is shown as a function of time for three different volumes. We chose 9 = 7r/16 and the volumes (8fm) , (16fm) , and (32fm) . For each volume, we plotted 3 3 3 the zero mode and the same non-zero mode. Notice that the zero mode is independent of the volume of the system, while the magnitude of the non-zero mode decreases with increasing volume. This is the signature that a real condensate has been formed. For a total volume of (10 fm) and 9 = 2TT/16 the time for relaxation from the initial 3 non-equilibrium state following the quench to the nontrivial ^-vacuum is approximately 10 fm/c~ 4 x 10~ s. This is of the same order of magnitude as the time which we might 23 expect the central region of the fireball in a heavy ion collision to be isolated from the usual vacuum. In order to be able to observe a nontrivial 9-vacuum in a relativistic heavy ion collision, it must form before the protective shell breaks down. As we have stated from the beginning, this is a simplified model. We have not Chapter 4. Production of Nontrivial Theta Vacua 71 Time (fm/c) Figure 4.20: The zero mode and a non-zero mode are shown as a function of time for three different volumes. The heavily dashed line represents the smallest volume, the medium dashed line the middle volume, and the solid line represents the largest volume. This graph illustrates the volume independence of the zero mode. attempted to correctly account for the collision geometry in heavy ion collisions since we have used a cubic lattice for our simulations. This fact is probably not relevant since the qualitative features of the behaviour are very robust under changes to the periodic lattice. This implies the behaviour is independent of the boundaries and as such would occur in any geometry. The volume we have used is just at the upper limit of what we would expect in heavy ion collisions at RHIC. Finally, we only included a phenomenological damping term. A complete calculation would have to include a better understanding of the effects of the expansion of the domain and emission of pions. However, our simplified calculation suggests the possibility of producing nontrivial 0-vacua in heavy ion collisions at RHIC. The most natural question to ask at this point is: "If a nontrivial 6-vacuum state can be formed in heavy ion collisions on a time scale over which they might be observable, Chapter 4. Production of Nontrivial Theta Vacua 72 how would we detect them?". One of the first suggestions of how to detect them was that, due to the fact that the 0-vacua are odd under charge conjugation times parity (CP), their decay must produce some CP odd correlations suggested in [32]. However, one should expect that the signal will be considerably (if not completely) washed out by the re-scattering of the pions and their interactions in the final states, which mimic true CP-odd effects. In practice it is quite difficult to overcome the problem of separating a true CP violation from its simulation due to the final state interactions[56]. A more promising direction is to look at 77(77') —» TTTT decays[32, 57] which are strongly forbidden in our world, but nevertheless, will be of order one if 9 / 0 due to the presence of a scalar component in the pseudoscalar mesons as discussed in Chapter 2. Calculations in [57] indicate that this effect could be quite noticeable and provide a definitive signature of a nontrivial 0-vacuum state. One could also look for electromagnetic decays of the pseudoscalar mesons which show the effects of shifted meson masses and lack of definite parity. This signature[57] might be easier to detect, as photons have a much better chance of penetrating the protective shell since they do not interact strongly. Finally, one could look for signatures related to the decay of the 9-vacuum state itself[57, 58]. These signatures range from production of axions to the enhancement of low energy photons and dileptons. The photons and dileptons could be produced indirectly, from the decay of the pseudoscalar mesons which are produced in predominantly low momentum modes or, directly, from the decay of an electromagnetic condensate coupled to the QCD condensate. In a recent paper[58] we obtained simple estimates for the TT , 77 and 77' 0 production in the decay of a 9-vacuum. These results show an amplification of production in the 10 MeV momentum range with many more 77' mesons produced than either ir° or 77 mesons. This would result in the enhancement of low energy dileptons which could possibly provide an explanation for the problem of the excess of low energy dileptons seen at CERN[59]. Indeed, it was shown in [60] that a sufficiently large enhancement of the Chapter 4. Production of Nontrivial Theta Vacua * 73 ?/-production could easily explain the excess of dileptons seen at CERN. In summary, 3 there are a number of different ways that nontrivial ^-vacuum states could be detected. Therefore we have shown, through a numerical calculation of the zero temperature equations of motion with a non-zero induced 0 parameter, that the chiral fields, 0, after a quench from high temperature go to a spatially constant non-zero value related to the 0-parameter. This occurs on a time scale of the order of I O -23 seconds. The fact that all other non-zero modes fall to negligible values long before, indicates that we have formed a condensate or a nontrivial 0-vacuum state. This prediction is testable in heavy ion collisions at RHIC. 3 Note that the rj' enhancement mechanism in [60] was very different from that described here. Chapter 5 Conclusions 5.1 Results The overall result of this thesis is a better understanding of the complexity of the QCD vacuum. The ^-dependence of the vacuum energy has been clarified by analysis in the anomalous effective Lagrangian approach. The form of the effective potential in this approach simultaneously satisfies the naively contradictory requirements that the vacuum energy should depend only on the combination 9/N c and should also be periodic in 9 with period 27T. Furthermore, we have discussed the fact that the 9-vacuum energy is determined by a small parameter proportional to the quark mass m , which is similar to the disoriented q chiral condensate. This is despite the fact that 9 appears in the effective Lagrangian associated with a large parameter independent of m . q This feature actually does not require the use of the anomalous effective Lagrangian we use, and is true for any effective chiral Lagrangian which respects the anomalous Ward Identities. This vacuum energy dependence is actually very significant to the possibility of production of nontrivial 9vacua in heavy ion collisions. We have presented results related to the phenomenology of the pseudoscalar mesons in a nontrivial 0-vacuum. We obtained a generalisation of the 77' mass formula expressed in terms of (G ) and 2 to non-zero 9. We obtained values for the masses in this model for the light neutral pseudoscalars for arbitrary values of 9. We obtained 77 — 77' and 7T — 77 — 77' mixing angles for non-zero 9. Finally we showed 0 74 Chapter 5. Conclusions 75 that for non-zero values of 9 the pseudoscalar mesons cease to be pure pseudoscalars, but acquire scalar components. This ^-dependence of the phenomenology of the pseudoscalar mesons is particularly important to the detection of nontrivial 0-vacua. Analysis of the vacuum structure of QCD in the anomalous effective Lagrangian approach leads to the realization that the vacuum structure may be very complex. The analysis shows the existence of metastable vacua (possibly a large number of them) for non-zero values of 9 and the possible existence of metastable vacuum states even for 9 = 0. This is in agreement with the predictions of supersymmetric models. There is even the possibility of stable solitonic configurations between all the different vacuum states. A l l of these different vacuum configurations are potentially producible in relativistic heavy ion collisions and the increased understanding of the vacuum structure of QCD from this research will be crucial to their observation. The decay rate for a single metastable vacuum for 9 = 0 to the true vacuum was calculated in both the low and high temperature limits. The results we obtained are given by Eqs. (3.119) and (3.147). These results are important for testing the theory in heavy ion collisions. In addition, the method used in these calculations represents a new approach that could be useful in other calculations of this type. The final results we have obtained concern the formation of nontrivial 9-vacua in heavy ion collisions. Our numerical results in a simplified model suggest that a nontrivial 9vacuum will likely be formed in heavy ion collisions. We also calculated a formation time of about 10~ seconds, which is short enough that this nontrivial ^-vacuum state could 23 be formed inside the shell of hot debris before the shell breaks down, while the interior region is still isolated from the exterior vacuum. This means that it may be possible to detect the nontrivial 0-vacua in heavy ion collusions. As has been emphasised so far in this conclusion, all of the results of this thesis could be testable in and have important implications for the physics of heavy ion collisions. Chapter 5. Conclusions 76 The even more exciting implication of this fact is that we may have the opportunity to study the physics of the early universe at the time of the QCD phase transition in experiments at RHIC in the near future. The presence of nontrivial vacuum structure could have important implications for this epoch in the evolution of the early universe[13, 14, 15, 16, 17, 18]. 5.2 F u t u r e Research There are many possible areas of future research which expand upon the results of this thesis. There are many methods by which the masses and the mixing angles of the pseudoscalar mesons we obtained could be altered. For instance, perturbative corrections to the mass matrix could have quite a significant effect on these values. The inclusion in the effective Lagrangian of terms of higher order in the quark mass would also lead to corrections. These corrections are expected to be quite small and are less important than perturbative corrections. It might be useful to determine how these corrections would affect these values and if they could bring about better agreement with measured masses. As was mentioned above, the method used in the calculation of the decay rate of the false vacuum in Chapter 3 represents a new approach. The method could be applied to other decay problems already analysed using different techniques to test their agreement. As well, there is currently interest in other types of domain walls in low energy QCD for which these techniques might be well suited. The inclusion of heavy degrees of freedom in the calculation of the false vacuum decay would be highly desirable for a complete understanding of the decay process. The inclusion of nucleons or heavy intrinsic degrees of freedom (glueballs) could have a significant effect on the decay rate of the false vacuum. Unfortunately, it is not clear Chapter 5. Conclusions 77 at present how to include these effects and whether they would complicate the problem beyond the hope of a solution. The chiral false vacua at 9 = 0 studied in Chapter 3 was only the most simple of the possible metastable vacua of this type. For q = 8 there are three other metastable vacua that one could study. As well, there is the added possibility of metastable solitonic configurations that interpolate between different vacua. All of these different vacuum states are potentially producible in heavy ion collisions and could have implications for the study of the development of the early universe. As was stated in Chapter 4, the numerical calculation presented there was done in a simplified model system. While the results are certainly important, it would be more satisfying to do the calculation in a more realistic model. The cubic lattice that we studied is highly unrealistic as the reaction volume in a heavy ion collision is shaped more like an elongated ellipsoid. The qualitative result of the formation of the condensate is almost certainly unaffected by the shape of the lattice. It is unlikely that the formation time of the condensate would be affected by the shape of the lattice but definitely possible. It would be useful to attempt to repeat this calculation taking into account the collision geometry and compare the results with those shown in this thesis. A potentially important improvement to this simple model would be developing a more physically motivated understanding of the way in which the system is damped. We have included the effects of expansion of the domain and emission of pions in a very simple phenomenological damping term. Better understanding of how these processes would affect the dynamics will almost certainly be important in determining the formation time of, and thus the likelihood of observing, nontrivial 0-vacua in heavy ion collisions. Another important area where more research is required is in the area of determining experimental signatures of the nontrivial vacuum structure. There have been a number of ideas put forward about how we might detect nontrivial #-vacua[32, 57, 58]. These ideas Chapter 5. Conclusions 78 must be studied in more detail and used to derive testable predictions for experimentalists at RHIC to measure. As well, it has been suggested that some of these signatures[57, 58] may already be present in existing data for lower energy heavy ion collisions performed at CERN[59]. Demonstration of this claim would certainly provide further impetus for research in this area. In particular, in a recent paper[58], we have obtained simple estimates for the 7r°, r\ and if production in the decay of a 9-vacuum state which show enhancement in the 10 MeV momentum range with many more 77" s produced than 7r°'s or 77's. This would result in an enhancement of low energy dileptons would could explain the excess of low energy dileptons seen at CERN[59]. In fact, it has already been shown in [60] that a sufficiently large enhancement of the rj' production could easily easily 1 explain the excess of dileptons seen at CERN. 1 Note that the 77' enhancement mechanism in [60] was very different from that described here. Bibliography [1] J. Donoghue, E. Golowich and B. Holstein, Dynamics of the Standard Model Cambridge University Press (1992). [2] P. Harris, et al, Phys. Rev. Lett. 82 (1999) 904. K. Smith, et al, Phys. Lett. B234 (1990) 191. [3] B. Borasoy, Phys. Rev. D61 (2000) 114017, hep-ph/0004011. R. Crewther, P. Di Vecchia, G. Veneziano and E. Witten, Phys. Lett. B88 (1979) 123. [4] R. Peccei and H. Quinn, Phys. Rev. Lett. 38 (1977) 1440; Phys. Rev. D16 (1977) 1791. S. Weinberg, Phys. Rev. Lett. 40 (1978) 223. F. Wilczek, Phys. Rev. Lett. 40 (1978) 279. [5] G. Veneziano and S. Yankielowicz, Phys. Lett. B113 (1982) 231. [6] T. Taylor, G. Veneziano and S. Yankielowicz, Nucl. Phys. B218 (1983) 439. [7] I. Halperin and A. Zhitnitsky, Phys. Rev. D58 (1998) 054016. [8] I. Halperin and A. Zhitnitsky, Phys. Rev. Lett. 81 (1998) 4071, hep-ph/9803301. [9] T. Fugleberg, I. Halperin, and A. Zhitnitsky, Phys. Rev. D59 (1999) 074023. [10] E. Witten, Ann. Phys. 128 (1980) 363. P. Vecchia and G. Veneziano, Nucl. Phys. B171 (1980) 253. [11] E. Witten, Phys. Rev. Lett. 81 (1998) 2862. [12] M . Shifman, Prog. Part. Nucl. Phys. 39 (1997) 1, hep-th/9704114. [13] J. E. Kim, Phys. Rep. 150 (1987) 1. P. Sikivie, "Dark matter axions '96" (1996) hep-ph/9611339. E. Shellard and R. Battye, Phys. Rept. 307 (1998) 227, astro-ph/9808220. E. Shellard and R. Battye, "Cosmic Axions" (1998) astro-ph/9802216. 79 Bibliography 80 [14] S. Banerjee, et al, "The Cosmological Quark-Hadron Transition and Massive Compact Halo Objects" (2000) astro-ph/0002007. N. Borghini, W. Cottingham and R. Vinh Mau, J . Phys. G 2 6 (2000) 771, hepph/0001284. A. Bhattacharyya, et al, Nucl. Phys. A 6 6 1 (1999) 629, hep-ph/9907262. H. Kim, B.H. Lee and C. H. Lee, "Relics of Cosmological Quark-Hadron Phase Transition" (1999) astro-ph/9901286.. E. Witten, Phys. Rev. D 3 0 (1984) 272. and references therein. [15] R. Brandenberger, I. Halperin and A. Zhitnitsky, "Baryogenesis at the QCD Scale", hep-ph/9903318, based on a plenary talk by R.B. at SEWM-98 and an invited talk by A.Z. at COSMO-98, published in the proceedings of SEWM-98. [16] S. Nayak and U. Yajnik, "Baryogenesis through Axion Domain Wall" (2000), hepph/0002219. [17] M. Forbes and A. Zhitnitsky, "Primordial galactic magnetic fields from axion domain walls at the QCD phase transition", hep-ph/0004051. [18] G. Fuller, G. Mathews and C. Alcock, Phys. Rev. D 3 7 (1988) 1380. D. Enstrom, "Gamma-Ray Bursts and Dark Matter - A Joint Origin?" (1998) hepph/9810335. B. Kampfer, Annalen Phys. 9 (2000) 1. A. Coley and T. Trappenberg, Phys. Rev. D 5 0 (1994) 4881. M. Hindmarsh, Phys. Rev. D 4 5 (1992) 1130. [19] G. Veneziano, Nucl. Phys. B159 (1979) 213. [20] V.A. Novikov, M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 1 9 1 (1981) 301. [21] I. Halperin and A. Zhitnitsky, Mod. Phys. Lett. A 1 3 (1998) 1955, hep-ph/9707286. [22] D.E. Groom, et al, The European Physical Journal C 1 5 (2000) 1. [23] A. Kovner and M. Shifman, Phys. Rev. D 5 6 (1997) 23961, hep-th/9702174. [24] R. Dashen, Phys. Rev. D 3 (1971) 1879. [25] A. Smilga, Phys. Rev. D 5 9 (1999) 114021, hep-ph/9805214. [26] E. Witten, Nucl. Phys. B 1 5 6 (1979) 269. G. Veneziano, Nucl. Phys. B159 (1979) 213. Bibliography 81 [27] F-G. Cao and A. I. Signal, Phys. Rev. D 6 0 (1999) 114012, hep-ph/9908481. L. Burakovsky and T. Goldman, Phys. Lett. B 4 2 7 (1998) 361, hep-ph/9802404. T. Feldmann and P. Kroll, Phys. Rev. D 5 8 (1998) 114006, hep-ph/9802409. E.P. Venugopal and B. R. Holstein, Phys. Rev. D 5 7 (1998) 4397, hep-ph/9710382. and references therein. [28] H. Goldstein, Classical Mechanics, Addison-Wesley Publishing Co. (1980) 143. [29] Z. Huang, K.S. Viswanathan and D.-d. Wu, Mod. Phys. Lett. A 7 (1992) 3147. [30] T. Fugleberg, "False Vacuum Decay in QCD within an Effective Lagrangian Approach", hep-ph/9906430. [31] M . Shifman, Phys. Rev. D 5 9 (1999) 021501. [32] D. Kharzeev, R.D. Pisarski, and M.Tytgat, Phys. Rev. Lett. 81 (1998) 512. D. Kharzeev, R.D. Pisarski, Phys. Rev. D 6 1 (2000) 111901, hep-ph/9906401. [33] W. Cottingham, D. Kalafatis and R. Vinh Mau, Phys. Rev B 4 8 (1993) 6788. [34] J. Baacke and V.G. Kiselev, Phys. Rev. D 4 8 (1993) 5648. [35] J. Baacke, Phys. Rev. D 5 2 (1995) 6760. [36] A. Strumia and N. Tetradis, Nucl. Phys. B542 (1999) 719. [37] J. Baacke and A. Surig, Z. Phys. C 7 3 (1997) 369. [38] L. Carson, X. Li, L. McLerran and R.T. Wang, Phys. Rev. D 4 2 (1990) 2127. [39] M.B. Voloshin, I. Yu. Kobzarev and L.B. Okun, Sov. J. Nucl. Phys. 20 (1975) 644. S. Coleman, Phys. Rev. D 1 5 (1977) 2929. [40] I. Kogan, A. Kovner and M. Shifman, Phys. Rev. D 5 7 (1998) 5195. [41] A. Linde, Phys. Lett. B 1 0 0 (1981) 37. A. Linde, Nuc. Phys. B 2 1 6 (1983) 421. [42] I. Affleck, Phys. Rev. Lett. 46 (1981) 388. [43] C. Callan and S. Coleman, Phys. Rev. D 1 6 (1977) 1762. [44] A. Vainshtein, V. Zakharov, V. Novikov and M . Shifman, Sov. Phys. Usp. 25(4) (1982) 195. [45] K. Kiers and M. Tytgat, Phys. Rev. D 5 7 (1998) 5970, hep-ph/9807412. Bibliography 82 [46] S. Coleman, The Uses of Instantons in Aspects of Symmetry, Cambridge University Press (1985). [47] S. Coleman, V. Glaser and A. Martin, Comm. Math. Phys. 58 (1978) 211. [48] K. Buckley, T. Fugleberg and A. Zhitnitsky, Phys. Rev. Lett. 84 (2000) 4814, hepph/9910229. [49] A.A. Anselm and M.G. Ryskin, Phys. Lett. B266 (1991) 482. J.-P. Blaizot, A. Krzywicki, Phys. Rev. D46 (1992) 246. J. Bjorken, Int. J. Mod. Phys. A 7 (1992) 561. [50] K. Rajagopal and F. Wilczek, Nucl. Phys. B404 (1993) 577. [51] K. Rajagopal, "The Chiral Phase Transition in QCD: Critical Phenomena and Long Wavelength Pion Oscillations", in Quark Gluon Plasma 2 edited by R. Hwa, World Scientific (1995), hep-ph/9504310. [52] D. Boyanovsky and H.J. de Vega, "Dynamics of Symmetry Breaking Out of Equilibrium: From Condensed Matter to QCD and the Early Universe" in ANSI-2000 published by the National Academy of Sciences of India, hep-ph/9909372. [53] I. Halperin and A. Zhitnitsky, Phys. Lett. B440 (1998) 77. [54] R. Burden and J. Faires, Numerical Analysis, PWS-KENT Publishing Co. (1989) 257. [55] W. Press, S. Teukolsky, W. Vetterling and B. Flannery, Numerical Recipes in C, Cambridge University Press (1992). [56] R. Peccei, "CP Violation: A Theoretical Review", hep-ph/9508389; "Overview of Kaon Physics" published in the Proceedings of Kaon99, hep-ph/9909236. [57] A. Zhitnitsky, "Signatures of the Induced O vacuum state in Heavy Ion Collisions", hep-ph/0003191. [58] K. Buckley, T. Fugleberg and A. Zhitnitsky, "Induced theta-vacuum states in heavy ion collisions: A possible signature", hep-ph/0006057. [59] G. Agakichiev et al, CERES collaboration, Phys. Rev. Lett. 75 (1995) 1272. N. Masera for the HELIOS-3 collaboration, Nucl. Phys. A 590 (1995) 93c. [60] J. Kapusta, D. Kharzeev and L. McLerran, Phys. Rev. D53 (1996) 5028, hepph/9507343. Bibliography [61] C. Callan, R. Dashen and D. Gross, Phys. Lett. B63 (1976) 172. R. Jackiw and C. Rebbi, Phys. Rev. Lett. 37 (1976) 172. [62] F. Sannino and J. Schechter, Phys. Rev. D57 (1998) 170, hep-th/9708113. [63] J. H. Kuhn and V.I. Zakharov, Phys. Lett. B252 (1990) 615. [64] I. Halperin and A. Zhitnitsky, Nucl. Phys. B539 (1999) 166, hep-th/9802095. [65] K. Intriligator, R.G. Leigh, and N. Seiberg, Phys. Rev. D50 (1994) 1092. K. Intriligator, Phys. Lett. B336 (1994) 409. [66] K. Intriligator and N . Seiberg, Nucl. Phys. Proc. Suppl. 45BC (1996) 1, th/9509066. Appendix A D e r i v a t i o n of the Effective P o t e n t i a l for Q C D The purpose of this appendix is to review the details of the derivation of the Effective Potential for the light fields in QCD. The first section reviews the derivation. The second section discusses some methods of determining two unknown parameters, p and q, of the Lagrangian. A.l A n o m a l o u s Effective L a g r a n g i a n for Q C D The derivation of the anomalous effective Lagrangian was originally presented in [8]. It relies heavily on the derivation of the anomalous effective Lagrangian for Yang Mills theory [7]. The Anomalous Effective Lagrangian is defined as the Legendre transform of the generating functional for zero momentum correlation functions. In Yang Mills theory the relevant operators are the G ^ G ^ (topological density operator) and G = G ^ G ^ 2 (gluon scalar current). The correlation functions of these operators are fixed in terms of the gluon condensate:(O. /TT G ) [20, 21]. Actually only the potential part of this S 2 Lagrangian can be constructed as this is the only part that can be fixed by the Low Energy Theorems. This effective Lagrangian, while not useful for calculating the S- matrix, is just what we need to study the vacuum properties. Construction of the effective Lagrangian is considerably simplified by defining new 84 Appendix A. Derivation of the Effective Potential for QCD 85 complex fields, H and H, which are linear combinations of the relevant operators: H = UiM * G where (3(a ) s i +l i> A | f i = G = -b o?J{2Tx)+0{al) i / ^ ) » _ JL ^ A , G . ) G j (A 155 is the Gell-Mann - Low /^-function for gluodynamics 1 YM and the first coefficient of the /3-function is given by b YM = (11/3)AT . £ YM C is a gener- ally unknown parameter which, parameterises the correlation function of the topological density: lim i Jdxe^(0\T^GG(x) ^ G G ( O ) } \0) YM = i YM ( ^ 2 G ) Y M . (A.156) The reason for using these composite fields has to do with the holomorphic structure of zero momentum correlation functions: lim i J dxe {0\T{H{x) if(0)}|0) = -1(H) , lim i = -4(F), iqx J dxe (0\T{H{x) iqx iJ(0)}|0) lim % J dxe (0\T{H(x) #(0)}|0) (A.157) = 0 , iqx where the dimensional transmutation formula: /—^±Q \ — 2 N 8TT 1 t cons 8^ M exp ( - , R 2 bgl which comes from renormalisation group arguments, has been used. (A.158) Holomorphy is analyticity of functions of complex variables. Holomorphy requires that a function of a complex variable, z, have independent functional dependence on z and z*. We can immediately see that the fields H and H decouple in the low energy theorems. This approach is motivated by fact that holomorphy is an exact property of the effective superpotential in the supersymmetric case. These details are mentioned here for completeness I n [9] t h e one-loop ^ - f u n c t i o n was used b u t m o s t o f t h e discussion is f o r m a l l y correct keeping t h e full /3-function. 1 Appendix A. Derivation of the Effective Potential for QCD 86 and have some relevance to details discussed in the next section but are not essential for the understanding. It can be seen that the n-point zero momentum correlation function of the operator H equals (—4) (.r7). Multi-point correlation functions of the operator H are analogously n-1 expressed in terms of its vacuum expectation value (H). At the same time, it is easy to check that the decoupling of the fields H and H holds for arbitrary n-point functions of H and H. The goal of the construction is an effective Lagrangian that reproduces at tree level all the Low Energy Theorems for the composite fields H and H. Consider the generating functional of connected Green functions with sources J and J : 2 exp [iW(JJ)] = £ [VA exp 'ez J n .9 + 27vn' / dxG + i * ™ f dxGG J 327T J 2 +iJ J dxH + iJ J dxH (A.159) where A are the gauge fields in terms of which G is defined. Effective zero momentum fields, h and h, are defined as: h " - ™ - h =ir- (jdxH) , Jdxh = - h (A ieo) which satisfy the equation: jdxh = (JdxH) . (A.161) The effective action F(h, h) is defined as the Legendre transformation of the generating functional W(J, J): T(h, h) = -W(J, J) + J dxJh + j dxJh. (A.162) Note the summation over the integer, n', which automatically ensures the 2TT periodicity in 6 and quantisation of topological charge. It is completely equivalent to the way that the vacuum angle, 6, originally appeared in Y M theory[61]. 2 Appendix A. Derivation of the Effective Potential for QCD 87 The definition of generating functional (A. 159), the Low Energy Theorems (A. 157), and their extensions to arbitrary n-point functions can be combined to obtain the following set of equations for W(J, J): Qn+l 8J n+1 Qn+l dJ n+1 W\ j = i jdxdx ...dx (T{H(x )...H(x )H(0)}) = (-A) J W\ = i Jdxdx ...dx (T{H(x )...H(x )H(0)}) = J= n =0 1 n 1 n J=J=0 1 n 1 dx(H), n n n (-4) Jdx(H), n Qk+l dJ dJ k l = 0. W\ J=J =0 (A.163) The solution to these equations is: j dx(H)e~ - ^ J dx(H)e- , W(J, J) = iJ (A.164) iJ which substituted into (A. 160) and solving for J, J gives: J =-HM) • =- Mw>J l - (A i65) These expressions can then be substituted back into equation (A.164) to obtain W as a function of h, h. The definition of T (A. 162) now becomes a differential equation for the effective potential, U(h, h) = -(l/V)T(h, h): ,8U TT u 1,, dU T - Th- m -l h h = ( h T x - - » + h) (A 166 One solution to this equation is: U^K'h) = \h\og^ + Ifciogi + fc), (A.167) where C, C and D are arbitrary constants. This function is not bounded from below and for complex argument is not single valued because of the branches of the logarithm. When dealing with the Veneziano-Yankielowicz effective Lagrangian [5] this problem was dealt with by summing over the branches of the logarithm [23]. In this case there exists yet another solution: U (h,h) = ~h\og(^\ 2 + - L M o g ( i ) +D(h-h), (A.168) Appendix A. Derivation of the Effective Potential for QCD 88 where a is an arbitrary real number. We consider the specific case: U {h, h) = \h L o g ^ + ^h L o g | + | 3 (k + n^j (h - h), (A.169) where Log stands for the principal branch of the logarithm, p and q (p/q = a) are relatively prime integers and n and k are integers in the ranges n = 0, ± 1 . . . and k = 0,1,. .., q — 1. This is the form of the function that will allow both the 9/N dependence and 27r periodicity discussed in Chapter 2. Therefore this is a natural place to discuss the introduction of the ^-parameter into the effective Lagrangian. The form of the generating functional (A. 159) contains an infinite sum over an integer parameter, n', which enters in the combination 6 + 2-Kn'. The only term in Eq. (A.169) which can accommodate the ^-parameter in this form is the last since n runs over all integers and the ^-parameter is introduced in the combination: 9 + 27rn. The final answer for the improved effective potential W(h, h) (here "improved" refers to the necessity of summation over the integers n,k in Eq.(A.170), see below) reads[7]: -ivwM) e = g . where the constants CYM,CYM g T 1* L L o g h + FI L O G J L | q 9 + 2nn\ h-h\ /, ,. can be taken to be real and expressed in terms of the vac- uum energy in Y M theory at 9 = 0, C Y M = C = -2eE^ \Q) M Y M = -2e(-b a /(32Tr)G ), 2 YM s and V is the 4-volume. The integer numbers p and q are relatively prime and related to the parameter f introduced in Eq.(A.156): q/p = 2£. Thus, we expect that the parameter £ defined in Equations(A.155), (A.156) is a rational number. This expectation is motivated by the fact that it turns out to be the case in all existing proposals to fix the value of £, to be discussed in the next section, and by experience with supersymmetric models. (In all likelihood, irrational values of £ would produce a non-differentiable 9 dependence for Y M theory.) On general grounds, it follows that p = 0(N ) and q = 0{N °) = 0(1). C C Appendix A. Derivation of the Effective Potential for QCD 89 Now this result needs to be generalised to the case of full Q C D w i t h Nf light flavours and N colours. In the effective Lagrangian approach the light matter fields are described 3 c by the unitary matrix, U^, corresponding to the phases of the chiral condensate: {¥ V ) L = -\{ty V )\U , R L R (A.171) lj with : 4 r- TT \ — 2 a a U = exp i\/2 rf , UU = 1 , M^JNf fv' /TT where A are the G e l l - M a n n matrices of SU(Nf), i\ a (A.172) + is the octet of pseudoscalar fields a (pions, kaons and the eta meson), n' is the SU(Nf) singlet pseudoscalar field (meson) Nf is the number of quark flavours. We use the values /„. = 132 MeV[22] and and f > = 86 M e V [ 2 0 ] for the meson decay constants . The indices, i and j , of U^ run over 5 v the quark flavor index, i,j = 1 Nf, which correspond to the quark flavours up(u), down(d) and strange(s), respectively. A s is well known, the effective potential for the U field is uniquely determined by the chiral anomaly and amounts to the substitution: 9^9-iTr (logU), (A.173) where T r refers to the trace of the matrix. In the sense of anomalous conformal W a r d Identities [20], Q C D reduces to pure Y a n g - M i l l s theory with the substitutions {G )YM 2 and b = bQCD —> b YM (G )QCD 2 —> Analogously, an effective Lagrangian for Q C D should transform to that of pure Y a n g - M i l l s theory when the chiral fields, U , are "frozen". Therefore the effective potential for Q C D is given by: e -ivmH,v) f + = g e x p (_ iV L L . __ (, q 9 - ilogdet U+ 2irn\ h - h + ITTV f k + 1 2 _ j _ _ o g h + £ L o h g \ i , , J + - y T r ( M U + h.c.) \ , TrrT /lt rrr (A.174) Note that the T?' is not really very light, but it enters the theory in this way. Note that mixing of the flavor eigenstates is ignored at this level. It should be noted that what we refer to as f /V2 is sometimes denoted f in the literature. As well, what we refer to as f i, is fn'/VH in the notation of [20]. 3 4 5 w n n Appendix A. Derivation of the Effective Potential for QCD where M = diac/(m |(tf tf )|), E = (ba /(32ir)G ), i i V is the 4-volume and the inte- 2 i 90 s gers p and q are relatively prime and related to the parameter £ from (A.155) (A.156) by q/p = 2£. T h e physical inputs to this equation are the values of the vacuum condensates {a /n G ) = 0.012 G e V , (tf"tf ) = ( t f t f ) 2 4 u d d s = -(240MeV) 3 and ( t f t f ) s 5 = —0.8 (240 M e V ) , and the quark masses (m , md and m ). 3 u s A s mentioned above, the sum over k prevents an ambiguity due to the multi-valued nature of the log function by summing over a l l branches of the effective potential[23][8]. The sum over n is the summation over a l l integer topological charges. T h e values of the parameters p and q are not definitely known as different proposals for their determining lead t o different values which we w i l l discuss i n the next section. Wherever possible we w i l l leave the values of p and q arbitrary but we should note that specific results i n Chapter 3 use the the values p = l l i V — 2Nf and q = 8 and the whole basis of the false c vacuum decay i n that chapter requires q / 1. T h e numerical results i n Chapter 4 use the values p = N and q = 1 but the theory of that chapter is valid for arbitrary p and q. c Integrating out the heavy "glueball" fields, h and h, by minimising the effective potential w i t h respect to these variables one obtains the effective potential for the light chiral fields U and U : 6 exp VE cos — -{8-i\ogDet U) + —l P P which is Eq.(2.29) and the starting point of all the research presented i n this thesis. One should note that most of the results of Chapters 2 and 4 are not sensitive to the specific functional form of the cos [—q/p (0 — i log Det U)\ term i n the effective potential. 6 Note that the summation variable k in (A. 174) has been renamed / in this formula. Appendix A. Derivation of the Effective Potential for QCD 91 However, the requirement of Veneziano [19] for the resolution of the U(l) problem: I d E (9) 2k vac d9 j 2k 9=0 / 2fc " i=l / dxi(GG(x{)... GG{x )) -! 2fc x ~ (—) 2k (A.176) *' X1 c supports the cos(9/N ) behaviour of the potential. c The next section w i l l discuss the values of the integer parameters p and q. A.2 D e t e r m i n i n g t h e values o f p a n d q In the previous section we saw the anomalous effective potential depends on the integer parameters p = O(N ) and q = 0(N 0 C °) = 0(1). There are different approaches to the fixing of these parameters that lead to very different vacuum structure. I w i l l discuss the physical implications of these parameters and mention some of the different proposals for fixing their values. The parameter p' is the number of branches of the effective potential. It also corl responds to the number of distinct vacua before the infinite volume limit is taken. T h e parameter q' determines the number of distinct vacua w i t h different energies that exist l before the infinite volume limit is taken. For q even there are (q + 2)/2 different energy vacua; for q odd there are (q + l)/2 different energy vacua. For q = 1 all of the p vacua are degenerate i n energy. T h e parameter ' g ' is also the number of physically distinct vacua that exist at 9 = 0 i n the infinite volume limit. It can be shown that the vacuum energy of the 0-vacuum has the following form: hrv E v { 9 ) = { 6 l ~3^ hcv G 2 l 9 ) = ( Q | " 327 ° 2 l 0 ) C O S { m = ~ m C ° s { m > (A - 177 ) where the £ is a parameter from the following low energy theorem: lim z Jdxe^(0\T{^GG(x) ^GG(0)} |0> = £ < ^ G 2 2 > • (A.178) Appendix A. Derivation of the Effective Potential for QCD 92 C o m p a r i n g equation (A.177) w i t h (2.31) for fa = 0 and 1 = 0 we identify 2£ = q/p. Fur- ther, E q u a t i o n (A. 177) can only have 27T periodicity if 9 is a rational number. Therefore restriction to q and p relatively prime integers is physically motivated. T h e first method of fixing the values of p and q comes from a self duality hypothesis i n [20] where it was suggested that gluodynamics should be a holomorphic i n the composite operators: G ± iGG. C o m p a r i n g w i t h (A.155) and using fi(a ) = -b MOs /(27r), 2 the 2 Y s G e l l - M a n n L o w /^-function for gluodynamics w i t h b YM = ( l l / 3 ) i V , we.find that this c corresponds to q/p = 2£ = 12/11TV . C T h e second approach is based on the analysis of softly broken supersymmetric theories by counting the number of degenerate vacua. A detailed discussion i n the case of pure gluodynamics is given i n [12]. T h e gluino condensate is given by: (AA) = exp(i9/N c + 2-Kk/N ), k = 0 , . . . , N — 1; this corresponds to N degenerate vacua. This sets the value c c c of 'p' at N . The vacua are a l l degenerate which sets the value of V at 1. The problem c w i t h this approach is that it is only valid when the gluon mass, m g <C A YM S where A.SYM is the dynamical mass scale; pure gluodynamics corresponds to the opposite l i m i t , m g 3> A.SYM- We should also mention that a non-standard non-soft supersymmetric breaking approach [62] leads to the value p = 11 N . c Another approach was motivated by a suggestion i n [63] that in Q C D w i t h massless quarks, nonperturbative m a t r i x elements should be holomorphic in the P a u l i V i l l a r s fermion mass, M . A s s u m i n g this holomorphy, they relate the proton m a t r i x element R of the topological density (p\GG\p) w i t h the corresponding m a t r i x element for the gluon scalar current (p\G \p). 2 A similar approach was used in [21] for pure gluodynamics to relate the zero momentum two-point functions of GG and G . 2 Pure gluodynamics was considered as the low energy l i m i t of a theory including a heavy quark and holomorphy i n the physical quark mass was argued to hold i n the m q the values p = 3b M = H i V Y c and q = 8. —> oo l i m i t . Holomorphy requires Appendix A. Derivation of the Effective Potential for QCD A n inverse route was used i n [64] to determine 'p' and q\ l 93 T h i s approach again relies on the fact that holomorphy of pure gluodynamics can be clarified by coupling it to a heavy fermion as in the previous method. T h e difference is that instead of starting from the coupled theory, in this approach the starting point is the holomorphic effective potent i a l of pure gluodynamics and a heavy fermion is introduced at the effective Lagrangian level by the "integrating in" procedure familiar in the context of supersymmetric theories [65, 66]. T h e consistency condition used i n [64], that the holomorphic structure of the pure gluodynamics should arise from the holomorphic structure of gluodynamics plus a heavy fermion, leads to the values p = 3b YM = H-ZVc, and q = 8. T h e agreement of this method w i t h the previous approach is encouraging. We do not pretend that the above discussion fully describes the methods or even lists all the possible ways of setting the values of the parameters p' and q'. l l O u r goal was simply to illustrate the uncertainty about the values and give the two m a i n possibilities. The results can be summarised as follows. The self duality hypothesis, non-standard S U S Y breaking and the consistency conditions of theories w i t h and without a heavy fermion lead to the values p = 3b = HiV YM values are then p = 3bQ D C = HN C c and q = 8 or 12. For Q C D the corresponding — 2Nj and q = 8 or 12. T h e soft S U S Y breaking approach leads to the values p = N c and q = 1. A s was stated at the beginning of this section, different values of these parameters can lead to some different qualitative results and are particularly important for the results of Chapter 3. Most of the results of Chapters 2 and 4 do not depend on specific values of these parameters. Appendix B Hyperspherical Harmonics i n Four Dimensions Hyperspherical coordinates i n 4 D dimensions are related to the Cartesian coordinates by: x = r sinO/0 sin(0) sin(<£), (B.179) y = r s i n ( ^ ) sin(0) cos(<£), (B.180) z = r sin(^)cos(0), (B.181) w = r cos(V>). (B.182) T h e L a p l a c i a n i n 4 D i n hyperspherical coordinates is: 1 1 —d (r d ) + . r r sin^ 3 3 r v r 7 2 1 1 (sin ^ ) + . dg (sin0<9<,) + . , ' r sin ?/> sin 0 r sin^sm 2 2 v 2 z 2 2 2 2 0 9 (B.183) Assuming separable solutions and treating 0 and <f> coordinates i n exactly the same way as i n three dimensions we obtain the differential equation: V"{ib) + 2 cotO/O^'O/O - 1(1 + 1) csc 0/0$(VO 2 = A * ( B . 1 8 4 ) W i t h the substitution u = cos(ip) this becomes: (1 - u )U"(u) - 3uU'(u) - ; }u(u) = BU(u). l 1 — u 2 l( l + l (B.185) If B = n(n + 2) and 1(1 + 1) = /'(/' + 1) this can be identified as the differential equation satisfied by the associated type II Chebyshev functions. These can be obtained from the 94 Appendix B. Hyperspherical Harmonics in Four Dimensions 95 type I I Chebyshev differential equation in exactly the same way as associated Legendre functions are obtained from the Legendre differential equation. As an aside notice that the type I I Chebyshev equation is a special case of the Geigenbauer (Ultraspherical) equation: (1 - x )-^Ci \x) 2 - (2a + l ) ^ ) ( i ) a - n(n + 2a)C£\x) = 0, (B.186) with a = 1. The Legendre polynomial equation corresponds to a = 1/2. Hyperspherical coordinates in all higher dimensions will lead to associated Geigenbauer equations with integer or half integer a. The hyperspherical harmonics in four dimensions are given by: Y {9A) Y im(9, <j>, ip) = A(n, I, m) I 0<l<n U (cos(xp)) l lm n , n [ Y\i\ {6,<f>) U - {cos{ibr)) n m where (B.187) - n < l < - l l l are the usual 3D spherical harmonics and U are associated Chebyshev Yi (6,(p) l n M type I I functions defined by: U (x) n l = (1 - x *y/^ (x) Un = (1 - x )'/ ;-l)"(n + l ) V i d 2"+ (n +1/2)! dx l 2 1 for I < -*V -^ {<i-* >" dx 2 1 1 /, !, +,/s } 0.188) > 0 and by: = - [ ( i - * 2 > n + I i • ^ for I < —2. These hyperspherical harmonics form a complete orthogonal basis for the functions of the angular variables in four dimensions: E E E* ^ r , M , V 0 W ^ y ^ ) = n=0 l=—n m=—l ROC JO R-K RN Jo JO JO ^^^'P^-^), sin sin tp 6 (B.190) r2-K ^rni ^i t ^) n'Vm'{ ^A,^) Y r ( ) Y r = S >S 5 . nn W W (B.191)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Vacuum structure of QCD in an effective Langrangian...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Vacuum structure of QCD in an effective Langrangian approach Fugleberg, Todd Darwin 2000-12-31
pdf
Page Metadata
Item Metadata
Title | Vacuum structure of QCD in an effective Langrangian approach |
Creator |
Fugleberg, Todd Darwin |
Date | 2000 |
Date Issued | 2009-07-23T17:37:16Z |
Description | The vacuum structure of QCD is studied using an anomalous effective Lagrangian approach. This approach makes it possible to determine how physical observables depend on the strong CP violation parameter, θ. The θ -dependence of QCD and the phenomenology of the light pseudoscalar mesons in this theory are illustrated. The vacuum structure of QCD is shown to be quite complex with the prediction of a number of different types of nontrivial vacuum states. Two specific examples of nontrivial vacuum states are analysed in more detail. The decay rate of a metastable vacuum state for the θ = 0 case is nonperturbatively calculated in both the zero and high temperature limits. The formation of a nontrivial vacuum state in heavy ion collisions is predicted in a simplified numerical model. These results have implications for the study of the evolution of the early universe near the QCD phase transition and may be tested experimentally very soon in heavy ion collision experiments. |
Extent | 4490681 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2009-07-23 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0085708 |
URI | http://hdl.handle.net/2429/11141 |
Degree |
Doctor of Philosophy - PhD |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2000-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2000-565467.pdf [ 4.28MB ]
- Metadata
- JSON: 831-1.0085708.json
- JSON-LD: 831-1.0085708-ld.json
- RDF/XML (Pretty): 831-1.0085708-rdf.xml
- RDF/JSON: 831-1.0085708-rdf.json
- Turtle: 831-1.0085708-turtle.txt
- N-Triples: 831-1.0085708-rdf-ntriples.txt
- Original Record: 831-1.0085708-source.json
- Full Text
- 831-1.0085708-fulltext.txt
- Citation
- 831-1.0085708.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0085708/manifest