V A C U U M S T R U C T U R E O F Q C D I N 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 T H E S I S S U B M I T T E D I N P A R T I A L F U L F I L L M E N T O F T H E R E Q U I R E M E N T S F O R T H E D E G R E E O F D O C T O R O F P H I L O S O P H Y in T H E F A C U L T Y O F G R A D U A T E S T U D I E S D E P A R T M E N T O F P H Y S I C S We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y O F B R I T I S H C O L U M B I A August 2000 © Todd Darwin Fugleberg, 2000 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of Br i t i sh 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 The University of Br i t i sh Columbia 6224 Agricul tural Road Vancouver, B . C . , Canada V 6 T 1Z1 Date: Abstract The vacuum structure of QCD is studied using an anomalous effective Lagrangian ap-proach. This approach makes it possible to determine how physical observables depend on the strong CP violation parameter, 6. The ^ -dependence of QCD and the phenomenol-ogy 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 anal-ysed 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 for-mation 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 Abs t rac t i i Table of Contents i i i L is t of Figures v Acknowledgement v i 1 In t roduct ion 1 1.1 Review of Quantum Chromodynamics and the ^-parameter 2 1.2 Motivation 9 1.3 Overview 12 2 The Anomalous Effective Lagrangian 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 Phenomenology of the Pseudoscalar Mesons for Non-Zero 9 25 3 False V a c u u m Decay 35 3.1 The Domain Wall Solution 38 3.2 Semiclassical Theory 40 iii 3.3 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 3.4 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 3.5 Summary 59 4 P roduc t ion of N o n t r i v i a l The ta Vacua 61 4.1 DCC 62 4.2 Nontrivial #-vacua 64 4.3 Numerical Evolution 67 5 Conclusions 74 5.1 Results 74 5.2 Future Research 76 Bibl iography 79 Appendices 84 A Der iva t ion of the Effective Poten t ia l for Q C D 84 A . l Anomalous Effective Lagrangian for QCD 84 A.2 Determining the values of p and q 91 B Hyperspher ica l Harmonics i n Four Dimensions 94 iv Lis t 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 Nf = 2 26 2.6 7T° mass as a function of 9 31 2.7 TJ mass as a function of 9 31 2.8 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>u = <j>d = <f>s. . . 39 3.13 The domain wall solution 40 3.14 Expansion of the partition function 43 3.15 The "potential" U"{<f>b(r)) 4 6 3.16 Comparison of the perturbation Vpert(r) and the approximate perturbation. 48 3.17 The phase shift <5O(^ A) shown as a function of UJ\ for various values of n. . 49 4.18 Time evolution of \4>k\ shown for various \k\ 69 4.19 The Fourier distribution of the field is shown at three different times. . . 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, Dr . Ar i e l Zhitnitsky, for his help and guidance. I have learned a great deal through working wi th h im and I am grateful to h im 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: Dr . Janis McKenna , Dr . John N g and Dr . Gordon Semenoff for their invaluable help and comments. I would particularly like to thank Dr . M c K e n n a for her advice on other career related matters. I would like to thank Lor i Paniak, Sebastian Jaimungal, K i r k Buckley, Mark Laidlaw, and Igor Halperin for their helpful advice and for many stimulating and enjoyable con-versations. I would also like to thank all of my other friends here at U . B . C . for making my time in Vancouver very pleasant. In particular I would like to thank Jaret, Sofia, Greg, Mateo, Al ison , Kev in and all the members of "Where's M y Head". Thanks are due to the entire faculty and staff of the Physics and Astronomy Depart-ment at U . B . C . for providing a stimulating and positive work environment. I would like to acknowledge that this work was funded in part by the Natural Sciences and Engineering Research Counci l of Canada. Words cannot express my gratitude to my family for the constant support and en-couragement they have shown me throughout my studies and, indeed, throughout my life. I would especially like to thank my wonderful wife Lisa and my fantastic son Nathan for keeping me sane during this demanding process. This thesis is a testament to their love and support. v i Chapter 1 In t roduct ion It is generally accepted that Quantum Chromodynamics (QCD) is the theory that de-scribes the strong nuclear force. In this theory al l hadrons are composite structures of particles called quarks, and quarks interact with other quarks v ia 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, al 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 in the high energy (weak coupling) l imit , such as in deep inelastic lepton-hadron scattering, and the theory agrees well wi th 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. As well, Q C D is asymptotically free, which means that the perturbative weak coupling l imit applies in the high energy regime while the realm of ordinary nuclear matter is the nonperturbative strong coupling l imit . 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 in the low energy l imit . 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 xAs well as by the exchange of the electroweak force carriers 1 Chapter 1. Introduction 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 Rev iew of Quan tum Chromodynamics and the ^-parameter QCD is an SU(NC) gauge theory2 defined by the QCD Lagrangian: 1 N f C = --G%G*r + (i-fDf - m 1 ^ ) % + G^G?, (1.1) where G^u is the field strength tensor of QCD, WQ are the quark fields, D^f is the covariant derivative and as is the QCD strong coupling parameter. The dual field strength tensor is defined as G^" — Sap^G^ where eappbV 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, Nc. The Greek indices p and v are Lorentz indices and run over the spacetime 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 Nc = 3 and this value is used in this thesis even though we leave the symbol Nc in the formulas so that Nc dependence is explicit. We will be concerned with low energy physics so we 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 - dvAl - QsUcA^Al, (1.2) 2SU(NC) is the group of Nc x Nc unitary matrices with unit determinant Chapter 1. Introduction 3 where g2s = A-Kas and fabc are the structure constants of the gauge group SU(NC). The quark covariant derivative is: where A"^ are the generators of SU(NC). QCD is a SU(NC) gauge theory which means that the Lagrangian is invariant under coupled SU(NC) 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: 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)yector x U ( iV/ ) A x i a l sym-metry. The U(Nf) matrices can be written as the product of the matrices corresponding to subgroups of U(Nf): (1.3) Va -+ e x p ^ A M j (1.4) % -> exp(zAA^ 7 5) \FQ U(Nf) = SU(Nf) ® U(l) (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 then3: Sl^JV/) Vector ® tf(l)vector ® SU(Nf)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 M * ) = 0. (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 3There is also the SU(Nc)c0[or 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: spontaneously breaks the ( / (A r / ) A x i a i subgroup of the symmetry group (1.6). This cor-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 broken4. 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: 4For 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. (1.9) ,3 (1.10) Chapter 1. Introduction 6 where Tr denotes the matrix trace over the generating matrices, A a of SU(NC). In particu-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 d4xC , (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 el{m-n)e{m\0\n). (1.16) m,n 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) = J VA.V^Vty O exp [i j d4x {£ + ^ G ^ G T } ] , (1.18) where we have used g2s = 47ra!s. 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: Q5 = Jd?x J 5 0 . (1.23) Chapter 1. Introduction 8 It can be shown that this charge generates changes of the 9 parameter: eiAQs\9) = \e-2NfA), (1.24) for arbitrary real constant A. This does not contradict our earlier comments because Q5 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) trans-formation, 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 mea-surements of the neutron electric dipole moment[2] which predict a value of 6 < 1 0 ~ 1 0 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 ] ) : 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 [3]. 9 = 9 QCD + arg(det M). (1.25) 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 M o t i v a t i o n 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 Nc Di Vecchia-Veneziano-Witten effective chiral Lagrangian ( W W ECL) for arbitrary 7VC. The W W 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 Nc limit, we reproduce the W W ECL. Furthermore, this approach to the derivation of 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 Nc. 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 (IT0 and rj). We obtain values for the r\ — rj 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 La-grangian, 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 theo-ries in the large Nc limit using the AdS/CFT correspondence[ll]. The same phenomenon 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 evolu-tion 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 for-mation of large scale magnetic fields[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 meaning-ful 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: which was argued for by Veneziano in order for the 17(1) problem to be resolved [19]. 1.3 Overview d2kEvac(9) d92k 2k (1.26) The vacuum properties of QCD are studied starting from an anomalous effective La-grangian. 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 ir0 — r\ — rf mixing angles for non-zero 9. Finally we can show that for 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 The Anomalous Effective Lagrangian and Implications The purpose of this chapter is to introduce the anomalous effective Lagrangian for QCD and discuss some of its implications1. In the first section we discuss the form of the 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 The Anomalous Effective Lagrangian The construction of the effective Lagrangian used in these studies was already accom-plished 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 deriva-tion 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), G2 = G^G^ (gluon scalar current) and tyty (quark condensate). The correlation functions of both gluon operators are fixed in terms of the gluon condensate: 1 M a n y of the results of this chapter were published i n [9]. 15 Chapter 2. The Anomalous Effective Lagrangian and Implications 16 (as/irG2) [20, 21]. Actually, only the potential part of this Lagrangian can be constructed 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 light2 matter fields are described by the unitary matrix, Uij, corresponding to the phases of the chiral condensate: (2.27) with 3. U = exp r ira\a 2 rf tV2 — h z- UU* = 1 , (2.28) fn yJNf frf J where A a are the Gell-Mann matrices of SU(Nf), ira is the octet of pseudoscalar fields (pions, kaons and the eta meson) and rf is the SU(Nf) singlet pseudoscalar field (meson). We use the values fn = 132 MeV[22] and fa = 86 MeV[20] for the meson decay constants4. The indices, % and j, of U^ run over the quark flavour index, i,j = 1 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: Weff(U, U+) = - Jim i log I £ exp v^°° v \l=0 VEcos - - ( 0 - « l o g d e t U) + —l V V + ^VTr(MU + MW) j , (2.29) 2Note 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. 3Note that mixing of the flavour eigenstates is ignored at this level. See Section 2.4 for more details. 4It should be noted that what we refer to as fn/V% is sometimes denoted fw in the literature. As well, what we refer to as fv>, is in the notation of [20]. 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 = bas/(32ix)(G2) 5 , b = y i V c — |7V/ is the first coefficient of the Gell-Mann - Low /^-function of QCD (see 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 branches6 of the potential (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(NC) and q = 0(NC °) = 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 — Nc, • q = 8 and p = 3bQCD = UNC - 27V), 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 non-supersymmetric theories7. The second set of values has been derived from a number of 5We use the values (as/-n G2) = 0.012 GeV 4 , (#^$ u) = (# d* d) = -(240 MeV) 3 and (* s * s ) = -0.8(240MeV) 3, except in Chapter 3 where we use ( $ s * s ) = -(240MeV) 3 . 6 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. 7See Appendix A for brief details and references. 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 / - 1 ) , will be unaffected. This fact was used in [10] to introduce the ^-parameter 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(el(t>i) 8 , becomes a piecewise smooth potential for the phases of the chiral 8 T h i s is not a restr ic t ion since the quark mass ma t r i x can always be diagonalised. Chapter 2. The Anomalous Effective Lagrangian and Implications 19 w, eff -71 / \ / \ / \ / \ / \K 4> Figure 2.1: Branches of the effective potential before the infinite volume limit is taken for 9 = 0, Nf = 1, q — 5 and p = 3bQcD-condensate, with cusp singularities: V P P P J cos</>i , / = 0,1,. . . ,p - 1, (2.31) if ( 2 / - i ) ^ <e-Yi<f>i<(2i + 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 20 we f f 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 The ta 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/Nc in order that the U(l) problem be resolved. With p = 0(NC) the effective potential appears to have the correct form and in fact has 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/Nc, but the fact that in the infinite volume limit only the lowest energy state survives 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: s i n ( ^ - ^ . ) = ^ f s i „ ( « . (2-33) They coincide with those of [10] to lowest order in 1/NC and can be used to find the solutions for the vacuum values of the chiral fields as a function of 9. For general values of Si = pMi/qE an analytic solution is not possible. For three light quarks (Nf = 3) in the realistic case, eu,Ed <C l,ss ~ 1, approximate solutions can be found as solutions of the equations: 0 f 0 ) = 0, 4 = 0 ) + 4 = 0 ) = (2.34) e u s i n ^ / = 0 ) = e d s in0 j = o ) , where 0(eu,Ed) terms have been neglected in (2.33). The solution to these equations is: • j,(i=o) _ mdsm9 S m 0 M = [ml + ml + 2mumd cos 9]V* + ^ £ d ) ' • JI('=O) m usin0 Chapter 2. The Anomalous Effective Lagrangian and Implications 22 sinrf=°> = 0(su,ed). 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 AQCD the lowest energy state is: which, when substituted back into the effective potential for I = 0, produces a formula for the vacuum energy as a function of 9: 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 mq and goes away in the chiral limit.9 It is also important for the results in Chapter 4 that the 0-vacuum energy 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 mq. 9 T h i s was shown i n [10]: physical ^-dependence vanishes i f any quark mass goes to zero. (2.36) (2.37) Chapter 2. The Anomalous Effective Lagrangian and Implications 23 For further discussion of the vacuum structure we will distinguish between two situ-ations: q = 1 and q 1. 2.3.1 q=l In the case where q = 1 we find that the solution of the corresponding equation of min-imisation 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]: mumd > ms\md - mu\, (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 25 2.3.2 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 Phenomenology of the Pseudoscalar Mesons for Non-Zero 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 The chiral condensate is parameterised in the following way: r 9 _ U = UQ exp iV2 — h I- (2.40) fix ' yjNf fv' where U0 is the solution to the minimisation equations, A a are the Ge l l -Mann matrices of SU(Nf), 7r a is the pseudoscalar octet, 7r 9 is the singlet field that is identified with the rf in the absence of mixing. The 7r a and 7r 9 fields have vanishing expectation values. 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 f l 6 7 r a 7 r 6 . (2.41) ^ a,6=3,8,9 Note that the indices do not run over the quark flavour group but over the U(3) indices associated with the Ge l l -Mann matrices 1 0 . The 7r 3 , TT8 and 7r 9 fields are associated wi th the Gel l -Mann matrices: ( 1 0 0 ^ A, = 0 - 1 0 0 0 0 1 0 0 1 0 0 \ ^ 1 0 0 ^ An = 0 1 0 0 0 1 (2.42) which are matrices wi th respect to the quark flavour basis. These fields behave like: ( U M — dd)/y/2, (uu + dd — 2ss)/y/6 and (uu + dd + ss)/\/3, respectively. These are the flavour eigenstates usually associated with the neutral mesons 7 r ° , n and n', respectively, in the absence of mixing. The components of the mass matr ix are: 2 m 3 3 = 2 (Mu cos (j>u + Md cos <f>d), f 2 J 7T 2 3/? (Mu cos 4>u + Md cos <f>d + 4MS cos 4>s), 1 0 We include in this term the identity matrix, denoted by A 9 , which is not properly a Gell-Mann matrix. Chapter 2. The Anomalous Effective Lagrangian and Implications 28 m; 99 2 m3S = m 8 3 = ( M u COS (j)u - Md COS (j)d) , (2.43) v3# ™ 3 9 = m 9 3 = 2\ITT ~T7~ (MU C0S ~ MD C0S ' V Nf Ujrf m89 = m 9 8 = 2\ r^r~F~r (Mu cos-0« + Md C ° S fa - 2 M S COS 0S) , V 37V/ J^Jjji where the fa are the vacuum values of the particular vacuum about which we are expand-ing. Notice that the 77199 component depends on 9 explicit ly and al l of the components depend on 9 impl ic i t ly through the vacuum values of fa. Neglecting mixing effects, mg 9 coincides with the physical mass of the rj. In fact for 9 = 0 and particular values, p = 36QCD and q = 8, we reproduce exactly the relation derived in [8]: fi m * = kNf<^ ° 2 ) ~ w ^ m { q > i % ) + ° { m ^ • ( 2 - 4 4 ) / u , d , s This relation shows that the 77' is heavier than the other pseudoscalar mesons because of the contribution from the gluon condensate. It also predicts that the strange quark contributes 20-40% of the mass of the rf. This suggests that perturbative corrections of 0(m2) (see below) could be quite significant. Taking the formal l imit Nc -> 00 and mq —y 0 leads to the mass relation in [21] where m2, scales like Nf/Nc, in agreement wi th [26]. If we instead take the l imit Nc —> 00 wi th fixed quark mass mq ^ 0, then m2, = 0(mq), as for ordinary pseudo-Goldstone bosons. The mass relation (2.44) gives the numerical result for m g 9 ( ~ m^/)=900 M e V . This result is encouraging as it is of the right order of magnitude, (see below) and is the first connection between the unrelated quantities, mvi and ( G 2 ) 1 1 . The fact that the mass matrix Mab defined by (2.41) has off diagonal terms, as seen in (2.43), means that the physical fields are not the TT0, fields, but linear combinations of 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 c, q = 1, and the standard values mu = 4 MeV, md = 7 MeV and ms = 150 MeV we obtain the values: m,o = 131 MeV, m, = 412 MeV, mv, = 1015 MeV. (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| 8 and m^. The 7 n | 8 term is dominated by the Ms term, because the strange 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>s) term as can be seen by recalling the vacuum solution (2.35). The variation of the ml9 term with respect to 9 is similarly determined by the Ms term because the term involving the larger parameter, E, is constant since Chapter 2. The Anomalous Effective Lagrangian and Implications 30 the vacuum solutions have (Yi 4>i — 9) = 0. The 7T° mass, however, is predominantly determined by the m | 3 term which is determined solely by the up and down quark terms which are strongly dependent on the ^-parameter as illustrated by the vacuum solutions (2.35) for (j)u and (j)d. This 0-dependence of the masses is one of the possible signals of the formation of a nontrivial 0-vacuum discussed in Chapter 4. In the l imit of pure SU(2) flavour symmetry (mu = md = m, which also means 4>u = <t>d) we can consider the simplest case of 77 — rj mixing. For 9 = 0 the mixing matrix is: \ m ( —^(2rns(ss)+rn(uu)) g j y ; (ms(ss) — m(uu)) 77-77-^(ms(ss) - m(uu)) - ^ ^ ( ^ + 7 ' which can easily be shown to coincide with an accuracy 0(m2) wi th the matrix given by Veneziano[26]: V )G2 , (2.46) J { mr)-r)' — \(\m\-m\) V 2V2, (2.47) where we replace the topological susceptibility, x, evaluated in pure gluodynamics, with the value in 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 in Q C D represents a significant improvement in 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: 36 ,OL = -f* K, + m2 - 2 7 7 4 ) + 0(m\ (2.48) This equation is a generalisation of the rj mass relation (2.44) to include the effects of 77 — rj mixing. The 77 — 77' mixing angle as a function of 9 is: Chapter 2. The Anomalous Effective Lagrangian and Implications 31 TC° mass (MeV) 0 -Tt -n/2 0 n/2 TC Figure 2.6: 7r° mass as a function of 9 r\ mass (MeV) 4JJ. 410 \ 40 8 \ 406 / 404 \ / 402 \ e -TC -n/2 n/2 Figure 2.7: r\ mass as a function of 9 r\' mass (MeV) < y ' l014 1013 \ 1012 \ 1011 \ 1010 \ -TC -TC/2 Tt/2 Figure 2.8: 77' mass as a function of 9 Chapter 2. The Anomalous Effective Lagrangian and Implications 32 where m | 9 is given above in (2.43) and, ra2/ and m 2 , are eigenvalues of the mass matrix for a given value of 9. For 9 = 0 and fa given by the vacuum solutions in this case, we obtain an r\ — 77' mixing angle of about —22°. This 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, TC0 — i) — fj' is defined by the eigenvectors of the mass matrix. The orientation of the mass eigenvector coordinate system wi th respect to the SU(3) flavor eigenstate system, #3, 7r 8 and 7r 9 , can be specified by Euler angles. For the specific case where 9 = 0 and mu = 4 M e V , md = 7 M e V and ms = 150 M e V we obtain the Euler angles[28]: 9Euler = l9°, (f>Euler = -0.83°, ^ Euler = 2.5°. (2.50) These values show that at least for 9 = 0 the mixing of the TT° meson wi th 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 in the case of the masses, the mixing 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 Pari ty Operation (P) for non-zero values of 9. For 9 = 0 they are eigenfunctions of Parity with eigenvalue —1, which means they are pure pseudoscalars. For 9^0 they cease to be eigenstates of Pari ty and become mixtures of scalar and pseudoscalar states 1 2 . Recalling the definition (2.27) of the chiral condensate matrix, U, and the diagonal 1 2 T h i s fact was previously noted i n [29]. Chapter 2. The Anomalous Effective Lagrangian and Implications 33 Euler Angle 9 (°) Figure 2.9: The Euler 9 angle for the mass eigenvector coordinate system with respect to the flavor eigenstate coordinate system as a function of 9. Euler Angle cp (°) - l •* -1 .5 \ '• -n/2 -2 .5 n/2 e Figure 2.10: The Euler (f) angle for the mass eigenvector coordinate system with respect to the flavor eigenstate coordinate system as a function of 9. Euler Angle \|/ (°) 9 g • 7 • • 6 5 \ 4 y -71 ^71/2 7C Figure 2.11: The 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: (*Li*'Rj)e = (*Li*Rj)o=o e*', (2.51) can be represented as a chiral rotation of the usual vacuum: $' = UR$Ul , = (tyLltyRj)e=o , (2.52) with: {UR)lk = foe** (Ul)k3 = 8 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: | 7 T ° ) ~ \u'i^5u' — d'i~f5d') = cos^ u {ui^u) — s in0 u \uu) — (u 4-> d). (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 fol-lowing 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 limits1. The existence of this 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 -/Vc limit using the AdS/CFT 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 Nc Yang Mills theory was estimated in [31]. The calculation of the decay rate of the false vacuum presented here differs from [31] because it is valid for finite iVc and because the heavy glueball degrees of freedom2 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 Nc is finite. Our calculation differs from [25] in the inclusion of the singlet rf field and 1 The results of this chapter are contained in [30] 2 We 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 c — 2Nj. As well, both of these estimates only 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 iVc limit was assumed. Even without consideration of these interesting applications, our method of calcula-tion of the determinantal prefactor is useful as an alternative to previous methods. Pre-vious calculations [33] [34] [35] [36] of bubble nucleation rates use a particular method for calculating the determinant ratio of operators of the form: M = -V2+ UJ2+ 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 = 3bQCD = 11ATC — 2Nf and q = 8 are 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 calcula-tion 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 E 1 — cos 1 — cos pf* r p TT Mf(<p) if <T5 > 0 M/(</>) if 0 < 0 / 27T \ cos JnJNf QNf — cos (3.57) The effective potential (3.57) has a global minimum at 4>+ = (TTfn)/'(2q^Nf) and a local minimum at = ~(Trfv)/(2q,jN~f), with a cusp singularity between them (see Figure 3.12). The minima are interpreted as two vacua separated by a high potential barrier (~ G2) which is fairly wide, while the energy splitting, AE, between the states is fairly small in comparison: 2TT ' AE = mqNf (tyty) cos qNf/ + 0 K ) (3.58) 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: Pfir <f>d.w.(x) = 2qJNf •— + 4arctan < tan f — ] exp\a(x — xo)l P I \4PJ if X < X0, Chapter 3. False Vacuum Decay 39 w eff Figure 3.12: Effective potential at 6 = 0 for equal chiral phases 4> = <j)u = (f>d = ^^Z— — — 4 arctan itan ( — | exp [— u(x — x0)] 2q^Nf [P I \4pJ where x0 is the position of the centre of the domain wall and: lcPWeff 2q^N~fVE if x > x0,(3.59) p = PU (3.60) is the inverse width of the wall. The solution (3.59) is shown as a function of x — x0 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(mqfl). (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 i J R <!>- • X - X , 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-B/h[l + 0{h)}. (3.62) The semiclassical approximation at zero temperature tells us that B is given by the Euclidean action of <fi: B = S4 I drd?x (3.63) where U(<f>) is the potential, Weff, for the chiral phases described in the previous section neglecting the energy difference, AE, between the two vacua: E i - c o s t e ^ - z ) if ^ > 0 • E l - c o s i f 0 < O . (3.64) Chapter 3. False Vacuum Decay 41 In order for this to be finite we must have l im^oo <p(r) = = —7rf7T/2qJNf where r — y/x2 + y2 + z2 + t'2. The solution of this problem is the four dimensional equivalent 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>d.w.(R - r) for r « R , (3.65) <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 3 • 2? • T r y ( l ~ cos ^ _ 27 A 2 * ; j ^ G 2 ) 2 4 ^ J " ^ M» ( l _ c o s ^ _ ) 3 - 256 m3 | (M>| 3 ' ( 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 f^fo] , 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 3 M / T ) -.exp ( - ^ ^ ) • (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 Quan tum Correct ions at Zero Temperature The quantum corrections at zero temperature correspond to the coefficient A in Equation(3.62)[43]: det[-d^ + U"(<f>b)} det[-dlld» + 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 LU2 = U"((/>-). These two parts of the spectrum must be analysed separately and it can be shown that this factor separates into three parts: Mwi/<MI)2) The first term comes from the zero eigenvalues. The second term comes from the nega-tive 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 (3.70) det ' [ -<V + U"((j>b) (w- 2) 5 d e t [ - c V + to2 (3.71) 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 contin-uous spectrum we divide the third term by a factor of co2 for each omitted eigenvalue. 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 Posi t ive Eigenvalues The contribution of the positive eigenvalues requires evaluation of the determinant ratio: det't-fyd" + U" {<!>„)] (3.72) (UJ-2)5 det[-d„d» + to2} However, since the (co~2)5 in the denominator corresponds to a set of measure zero in 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>b)\ = exp [Tr log {-d^ + UJ2 + Vpert(r)] = det| — (3.73) 2 I -UJ X exp Tr Vyert(r) -d^d* + CO2 1 2 _ p e r t ( r ) _ N d^ + co2, ' 1 + 77 + / Vpert(r) \ 3 [-8^+CO2) 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 po-tential (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) - c V + UJ2 Tr Tr -d^ + LO2 Vpert(r) j =J J 4 T , >4 Vpert(k) Vpert(-k) P (k2 + co2) ((k + p)2 + UJ2) Vpert(r) V I _ f d 4 p d 4 q Vpertjk) Vpert(p) Vpert(-k - p) -d^d" + LO2 (3-74) (3.75) 43.76) (k2 + UJ2) ( p 2 + UJ2) ((k+p + q)2 + uj2f where Vpert(k) is the Fourier transform of the Vpert(r). We only get a finite answer if we 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(r, <f>, ip) in terms of the 4D 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: ^ -> 1372 dr2 An2 + 8n-l-3' 4r2 r3/2 = J_T> r3/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[—Vn + UJ2]. In order to calculate this determinant we need to solve the eigenvalue equation: [Dn -co2 + 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) = {ujxrf'2jn{ujxr) = ^J^-Jn+1{toxr), (3.81) where UJX = y/X — LO2. jn are the 4D analogs of the spherical Bessel functions and Jn+i(x) 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 2 and indeed there are no solutions for values of A < to2 which are well behaved at r = 0 and r = oo. Therefore the continuous spectrum of eigenvalues can be written as A = UJ2 + co 2 and: det[-Vn + co2] = Y[{UJ2+ uo2x). (3.82) The numerator of (3.72) involves the "potential" U"(fa(r)). For the symmetrised ef-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 pertur-bations 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, cdxr = ujxr + 5(ujx). In this situation, for each value of n, we can obtain the ratio of determinants for a discrete spectrum by [44]: n 4 ± 4 = exp (E In ^ ± 4 ) * exp fc (f * " *') • (3.83) ^^UJZ + UJ{ y—1 uolJrUJlx) y—1 uo l + uj lx J which becomes in the continuum: exp 1 /•«> 2LOx5{UJX) UJ1 + uoi -I TX JO (3.84) Chapter 3. False Vacuum Decay 46 204000 203000 0 R Figure 3.15: The "potential" U"(<f>b(r)). 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: d2 An2 + 8n + 3 dr2 d2 An2 + 8n + 3 AT2 + (A - UJ2) dr2 Ar2 (A - UJ2 + Vpert{r)) u{r) = 0, v(r) = 0. (3.85) (3.86) Multiply (3.85) by v(r) and (3.86) by u(r), subtract, and integrate over r from 0 to oo to obtain: roo poo / (u"(r)v(r) - v"(r)u(r)) dr = (u'(r)v(r) - v'(r)u(r)) |g° = / Vpert(r)u(r)v(r)dr. Jo Jo (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, S(u)\): 1 r°° sin5(ujx) - — / Vpert(r)u(r)v(r)dr. LUX Jo (3.88) Since the perturbed differential equation is extremely difficult to solve, we use perturba-tion theory to obtain: 1 f 0 0 7T / * ° ° 8{ux) ~ 6{ux)0 = — / Vpert{r)u(r)2dr = - / Vpert(r)rJn+l(uuxr)2dr, (3.89) ld\ JO Z JO for small phase shifts. Using this result for the phase shift in Equation (3.84), we obtain the ratio of determinants: det[-2>n+ £/"(&)] det[-X\ + LU2' exp roo re Jo JO dto\dr u)Xr UU2+LU2 Vpert{r)Jn+l{LOXrf (3.90) for each value of n. The complete determinant ratio then becomes: d e t f - d ^ + LU2' exp £ ( 2 n ( n + 2) + l) n=0 °°rCOdLUXdr U x r " <--^2 0 JO uc ~r x n c o , 2VPert{r)Jn+l{LUXry LUZ + LUX (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 2 i i 1 0 > r > n , 0 otherwise (3.92) Chapter 3. False Vacuum Decay 48 Figure 3.16: Comparison of the perturbation Vpert(r) and the approximate perturbation. The perturbation Vpert(r) is the dotted line and the approximate perturbation Vpert(r) is the solid line. Chapter 3. False Vacuum Decay 49 0.03 0.025 0.02 0.015 0.01 0.005 0.01 0.008 >0 0 .006 0 . 004 0.002 400 800 1200 n=0 400 800 n=5 1200 CO, 0 . 006 50 0.004 0 . 002 0 . 006 0 .005 0.004 >0 0.003 0 . 002 0.001 400 800 1200 n=10 c o , 400 800 n=15 1200 CO, Figure 3.17: The phase shift 8Q(ux) shown as a function of tox for various values of n. Notice the changing vertical scale. where B = (£/"(</>(0)) — U"(4>(R)))(see Figure 3.16). With this approximate potential we find the phase shift via (3.89): 80(UJX) = -200-^ \^r2R2 {Jn+i(rujxf - Jn(rujx)Jn+2(rux)) -4-2n^4 r4(rujx) |2+2n , / 3 \ (n + 4)!(„ + T ) T | 2 ( n + 2)(n + 4) (n + n + 4, 2n + 3; -rWx) -r2Lj\ XF2 (n + |, n + 5, 2n + 4; - r 2 ^ ) J 2+2n 11 2-2-2nr3R(rux) 10 (n + f)[(n + l)!]2 ( 3 5 7 \ n + -, n + -, n + 2, n + -, 2n + 3; -r2u2XJ -i 11/lOfl . (3.93) Using this result we can verify that 5(tox) <C 1 which means that the approximations of (3.89) are justified. The phase shift as a function of u>x is shown for a few values of n in Figure (3.17). However, (3.93) is very difficult to work with so we make another approximation. For cox large enough, we can use the asymptotic approximation for the Chapter 3. False Vacuum Decay 50 Bessel function in (3.89), Jn+i{x) = \ — \pn+i(x) V 7TX I COS X -{n+D £]-Gn+i(*)sin [ x - ( n + | ) | ] } , (3.94) where: 1 _ (p-l)(p-9) + ( M - 1 ) ( M - 9)(/i - 25) (ja - 49) 2! (8xy 4! (8x)4 n ^ - k ^ 1 ) (M-1)0"-9) ( / / -25) , Vrc+lW — „ , / „ x „ . / „ x , 1-1! (8a;) / i = 4(n + l ) 2 . 3! (8xY In this way we find: <*oM 307ru; A TT ui\ x > , (3.95) (3.96) (3.97) (3.98) This estimate for 80(LJ\) has an infrared cutoff in u\, below which it is not a good 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 co\ integral at the position of the first peak in 50(id\). Using this approximation for 50(id\) we obtain a first estimate of the ratio of determinants for each value of the quantum number n given by: det[-2?w + U"{<j>h)] det[-Vn + £/"(</>_)] exp IRQ 30™ arctan id (3.99) where: u)(n) «33.16n + 72 = an + b, (3.100) is the infrared cutoff given by the location of the first peak in S0(idX). 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: 4RP n2 30TTU; \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 approx-imation for the phase shift: 1 3 6(LOX) = arcsin(<50(u;A)) = S0(LOX) + - 5 0 ( W A ) 3 + — S0(LOX)5 + (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: \ 3 _ exp TT2 f°° , LOX50(LOX)3 24 /0 d U X -(3.103) co2 + LO2X 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: det[-Vn +U"((f>b) exp TT2 {2RP\3 f LO —— —— <, , \ — arctan 24 V 30vro; LO LO(n) det[-Vn + [/"(</>_)] For large values of n the terms of the sum in (3.91) approach: Mn). (3.104) (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: ,3 7T2 / 2R8 \ 6 X P 24 3 ^ J 3 r £ ( 2 n ( n + 2) + l) LO Ln=0 Lo(n) arctan LO ~ £ 2 n 2 LO 71=1 3a3n 3 3.106) TT2 (2RR X 100 6 X P 24 [SOTTUJ J 100 - j N ( 2 n ( n + 2) + l) 71=1 I where £ is the Riemann Zeta function: LO y"(2n(n + 2) + 1) < - arctan 3 4 w 3 LO LO" 4 W , . . LO" ^ oo -i 71=1 Evaluating gives: exp -5300 TT2 [2RPY exp[-2 x 10~5 . (3.107) (3.108) (3.109) 24 \30TVLOJ Because of the small value of the exponent we can see that the exponential is extremely well approximated by: exp -53007T2 (2Rp\3' 24 \ 307TLOJ _ 53007T2 (2RpV 24 [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 resum-mation of the perturbation expansion (3.73) of the determinant. While we make an ap-proximation through Eq.(3.102), this expansion is nonperturbative since the terms of the expansions do not coincide. Our approximate calculation should therefore have contri-butions 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 and Negative 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 dv d x^ d Xo>(x, V, z, t) = ^Mr) = TrMr) = T ~^Mr), (3.111) where r = y/x2 + y2 + z2 + t2 and: t = r cos(V>) ~ Y100, (3.112) z= r cos(0) sin(V>) ~ Ino, (3.113) y = r cos(0) sin(0) sin(^) ~ Y1U + Yn„u (3.114) x= r sin(0) sin(0) sin(t/>) ~ Yin - r n _ i . (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 this3 and determine its value. We use the method of 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): An2 + 8n + 3 ^pn ^p ~r~ AR2 (3.116) where Xp is a number independent of n. We know that for n = 1 the lowest eigenvalue is zero: A 01 An + 15 AR2 0. Therefore we can obtain the lowest eigenvalue for n = 0: 3 A_ = A, oo A 0 + 15 3 + _3_ R2' (3.117) (3.118) AR2 AR? AR2 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 Rate for Zero Temperature 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 = exp (-S 4) ( u;4 ( ^ LU det'i-d^ + U"(<f>b)] 1 / 277rV\ (27ir2oALu2\2 RLU 1.55 x 10" 1 8MeV 4 * 3 x l O ^ f n T V 1 . {LU-2)5 det[-0„d" + LU2} 5300TT2 (3.119) 24 2RpY 307m; j -1/2 (3.120) (3.121) 3This is most likely a good assumption but not proven due to the cusp in our potential. For more details see [43],[46] and [47]. 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 Quantum Correct ions for H i g h Temperature The quantum corrections to the decay rate at high temperature are[41]: which again factors into three parts corresponding to zero, negative and positive eigen-values respectively. 3.4.1 Posi t ive Eigenvalues The calculation in three dimensions is extremely similar. The expansion of the determi-nant 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: det'[-$# + £/"(&)] 1/2 (3.122) (3.123) Tr (k2 + UJ2) ((k + p)2 + UJ2)' (3.124) Chapter 3. False Vacuum Decay 56 The calculation is done using spherical coordinates and the eigenfunctions are ex-panded in terms of the spherical harmonics: 0 0 7 1 1 1 11 (r) x(r,M) = £ £ E G m ^ Y U M ) . n=0 l = - n m=-\l\ T (3.125) The Laplacian, when operating on each term of (3.125) with quantum number I, becomes: dldl -> -r d2 1(1 + 1) dr2 r2 1 r = -T>ir. r (3.126) Therefore: deti-did1 + U"} = TJ(det[-2?i + U"})2l+l. (3.127) 1=0 The solutions to this differential equation that are well behaved at r = 0 and r = 0 0 are: u(r) = uxrji(ujxr) = ^^Jl+l/2{ooxr), (3.128) where jt are the usual 3D spherical Bessel functions and LOX = y/\ — LO2. We obtain: S(LOX) W — / Vpert{r)u{r)2dr = - / Vpert(r)rJl+1/2(coxr)2dr, (3.129) L0X JO Z JO thus giving4 5 : \det[-didi + U"(<l>b)]\ deti-did1 + a* = exp^(2/ + 1) 1=0 roo /-oo D\r duoxdr A ^ * •<•> ^ 2 ROO RE Jo Jo LO2 +L0l ;Vpert(r)Jl+1/2(ojxrY (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(LOX) W S0(UX) 2Rf3 30TTLOX + 0(l/u2). (3.131) 4Note that we have again dropped the factor of (to 2 ) 3 and the "det"' notation, as the omitted eigenvalues correspond to a set of measure zero. 5The 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]. Chapter 3. False Vacuum Decay 5 7 A first estimate of the ratio of determinants for each value of the quantum number, 1, is given by: det[-Vt + U"(fa)} exp 2R(3 307TCJ arctan LO ( 3 . 1 3 2 ) det[-2?, + U"{4>-)} where 0(1) 4 9 . 5 / + 8 2 = cl + d is the infrared cutoff given by the location of the first peak in 50(LO\) as a function of LO\. For large values of I the terms of the sum in ( 3 . 1 3 0 ) approach: M , ( % ) . (3.133) 307ro; \cL 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 . 1 3 0 ) . Again we adjust our earlier approximation for the phase shift: 1 3 S(ux) = arcsin(<50(a;A)) = ^o(wA) + -S0(LOX)3 + —S0(LOX)5 + • ( 3 . 1 3 4 ) and the correction is implemented as in ( 3 . 1 0 3 ) . Evaluating the contribution to the determinant ratio leads to: 3 det[-Vi + U"((f>b)] det[-Vt + f/"(0_)] exp TT2 (2RP" 2A I 30TTO; LO LO (0 — arctan ui MO, ( 3 . 1 3 5 ) For large values of I the terms of the sum in ( 3 . 1 3 0 ) approach: TT2 ( 2Rp X LO ( 3 . 1 3 6 ) The sum is finite so we do not need to further normalise and the result for the complete determinant ratio ( 3 . 1 3 0 ) is: Tr 2 (2RP\Z rl0° exp 2 4 \30TTU LO W). \ r 0 0 r LO ( 1 0 0 f , ,3 1 O, ,3 , ,3 ( 3 . 1 3 7 ) Chapter 3. False Vacuum Decay 58 Evaluating gives: exp TT2 (2Rfi\ 82— 24 J exp[3 x 10~7;. (3.138) Because of the small value of the exponent, we can see that the exponential is extremely well approximated by: exp 82TT2 (2Rj3\: ~2A V307ra;) 827T2 (2R(3Y 24 \307ruj) (3.139) This result constitutes the contribution of the positive eigenvalues to the quantum cor-rections to the decay rate in the high temperature theory. In the next section we consider the zero and negative eigenvalues. 3.4.2 Zero and 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 Mr)> (3.140) where r = y/x 2 + y2 + z2 and: z = r cos(0) ~ Yin, y = r cos(</>) sin(0) ~ Yn + x = r sin(^) sin(0) ~ Yu — l i - i -(3.141) (3.142) (3.143) 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: 1(1 + 1) AP„ = A p + ( 3 . 1 4 4 ) where A p is a number independent of /. We know that for / = 1 the lowest eigenvalue is zero: Aoi = A 0 + = 0. (3.145) Therefore we can obtain the lowest eigenvalue for I = 0: A_ = A 0 0 = A 0 = - - ^ . (3.146) 3 . 4 . 3 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: 827T 2 (2RB\ 2 4 U O T T U J 31 "1/2 Rw 3 / l 6 7 r g 1 ( T ) 3 y / ^ / 1 6 7 ^ ( T ) 3 \ R / V = ^[U [HAEW) )exprw7] (3.147) We cannot obtain a numerical estimate for the decay rate at this time since the temper-ature 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 realized1. This nontrivial vacuum is very different from the metastable vacuum discussed 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~10[3] so the 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 publ ished 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 = Nc and q = 1 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 mq and negligible at high temperature. 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 (~ mq) in the a 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 ~ mq between the DCC state and the vacuum state. The high energy products of the collision expand outwards at relativistic speeds. In-side 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 temper-ature 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\ (¥LVR) = - | ( * A > | U^. (4.148) The energy density of the DCC is determined by the mass term in an effective Lagrangian: = -\Tr{MU + M^U ]) = -2m|( t f t f ) | cos(^), (4.149) where we put mu = = m for simplicity. Eq.(4.149) implies that any (f> ^ 0(mod 2ir) is not a stable vacuum state because -g^l^o 0; a n d therefore the vacuum is misaligned. On the other hand, the energy difference between the misaligned state and true vacuum (<f> = 0) is small and proportional to mq. Therefore, the probability to create a state with an arbitrary <j> at high temperature T ~ Tc 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 fundamen-tal 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 = e i0, (4.150) where — fa is 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)- . (4.152) q 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 ^> mq\(tyty)\ 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(mg) 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/Nf). (4.153) This shows that the energy cost of creating a nontrivial 9-vacuum is proportional to the much smaller parameter, Mi oc mg, as in the case of the DCC. This is a very important 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 fields2, 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 mq from (4.153). Now we study the evolution of the chiral fields with time. For this study we assume 2 As 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. Chapter 4. Production of Nontrivial Theta Vacua 6 7 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 chiral fields relax 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 Numer i ca l Evo lu t i on 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- V2fa + yfa + ^-V(fa,e) =0, i = l,2 (4.154) where V 2 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 ~ AQCD ~ 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 163 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~bMeV~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 non-vanishing 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> fields is 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> fields were 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 . 1 5 0 . 1 t2=2t 0 . 0 5 t1 =1.6x10-3 MeV"1 o 5 10 15 20 Magnitude of k (arbitrary units) 25 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)3, (16fm)3, and (32fm)3. For each volume, we plotted 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)3 and 9 = 2TT/16 the time for relaxation from the initial non-equilibrium state following the quench to the nontrivial ^ -vacuum is approximately 10 fm/c~ 4 x 10~23s. This is of the same order of magnitude as the time which we might 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 0 , 77 and 77' 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 ?/-production3 could easily explain the excess of dileptons seen at CERN. In summary, 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 - 2 3 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 anoma-lous 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/Nc 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 mq, which is similar to the disoriented chiral condensate. This is despite the fact that 9 appears in the effective Lagrangian associated with a large parameter independent of mq. 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 9-vacua 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 2) and 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 7 T 0 — 77 — 77' mixing angles for non-zero 9. Finally we showed 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. Al 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 9-vacuum will likely be formed in heavy ion collisions. We also calculated a formation time of about 10~23 seconds, which is short enough that this nontrivial ^-vacuum state could 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 Future 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 pseu-doscalar 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. Al l 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' production1 could easily easily explain the excess of dileptons seen at CERN. 1Note 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 Com-pact Halo Objects" (2000) astro-ph/0002007. N. Borghini, W. Cottingham and R. Vinh Mau, J . Phys. G26 (2000) 771, hep-ph/0001284. A. Bhattacharyya, et al, Nucl. Phys. A661 (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. D30 (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), hep-ph/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. D37 (1988) 1380. D. Enstrom, "Gamma-Ray Bursts and Dark Matter - A Joint Origin?" (1998) hep-ph/9810335. B. Kampfer, Annalen Phys. 9 (2000) 1. A. Coley and T. Trappenberg, Phys. Rev. D50 (1994) 4881. M. Hindmarsh, Phys. Rev. D45 (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. B191 (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 C15 (2000) 1. [23] A. Kovner and M. Shifman, Phys. Rev. D56 (1997) 23961, hep-th/9702174. [24] R. Dashen, Phys. Rev. D 3 (1971) 1879. [25] A. Smilga, Phys. Rev. D59 (1999) 114021, hep-ph/9805214. [26] E. Witten, Nucl. Phys. B156 (1979) 269. G. Veneziano, Nucl. Phys. B159 (1979) 213. Bibliography 81 [27] F-G. Cao and A. I. Signal, Phys. Rev. D60 (1999) 114012, hep-ph/9908481. L. Burakovsky and T. Goldman, Phys. Lett. B427 (1998) 361, hep-ph/9802404. T. Feldmann and P. Kroll, Phys. Rev. D58 (1998) 114006, hep-ph/9802409. E.P. Venugopal and B. R. Holstein, Phys. Rev. D57 (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 Ap-proach", hep-ph/9906430. [31] M. Shifman, Phys. Rev. D59 (1999) 021501. [32] D. Kharzeev, R.D. Pisarski, and M.Tytgat, Phys. Rev. Lett. 81 (1998) 512. D. Kharzeev, R.D. Pisarski, Phys. Rev. D61 (2000) 111901, hep-ph/9906401. [33] W. Cottingham, D. Kalafatis and R. Vinh Mau, Phys. Rev B48 (1993) 6788. [34] J. Baacke and V.G. Kiselev, Phys. Rev. D48 (1993) 5648. [35] J. Baacke, Phys. Rev. D52 (1995) 6760. [36] A. Strumia and N. Tetradis, Nucl. Phys. B542 (1999) 719. [37] J. Baacke and A. Surig, Z. Phys. C73 (1997) 369. [38] L. Carson, X. Li , L. McLerran and R.T. Wang, Phys. Rev. D42 (1990) 2127. [39] M.B. Voloshin, I. Yu. Kobzarev and L.B. Okun, Sov. J. Nucl. Phys. 20 (1975) 644. S. Coleman, Phys. Rev. D15 (1977) 2929. [40] I. Kogan, A. Kovner and M. Shifman, Phys. Rev. D57 (1998) 5195. [41] A. Linde, Phys. Lett. B100 (1981) 37. A. Linde, Nuc. Phys. B216 (1983) 421. [42] I. Affleck, Phys. Rev. Lett. 46 (1981) 388. [43] C. Callan and S. Coleman, Phys. Rev. D16 (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. D57 (1998) 5970, hep-ph/9807412. Bibliography 82 [46] S. Coleman, The Uses of Instantons in Aspects of Symmetry, Cambridge Uni-versity 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, hep-ph/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. A7 (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 Equi-librium: 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, hep-ph/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. A p p e n d i x A Der iva t ion of the Effective Poten t ia 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 Anomalous Effective Lagrangian 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 G2 = G ^ G ^ (gluon scalar current). The correlation functions of these operators are fixed in terms of the gluon condensate:(O.S/TT G2) [20, 21]. Actually only the potential part of this 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 = U i M G * + l i i>GA | f i = i / ^ ) G » _ j JL ^ G A , (A.155) where (3(as) = -bYMo?J{2Tx)+0{al) is the Gell-Mann - Low /^-function for gluodynamics1 and the first coefficient of the /3-function is given by bYM = (11/3)ATC. £YM is a gener-ally unknown parameter which, parameterises the correlation function of the topological density: lim i Jdxe^(0\T^GG(x) ^GG(O)} \0)YM = iYM ( ^ G 2 ) 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 dxeiqx{0\T{H{x) if(0)}|0) = -1(H) , lim i J dxeiqx(0\T{H{x) iJ(0)}|0) = - 4 ( F ) , (A.157) lim % J dxeiqx(0\T{H(x) #(0)}|0) = 0 , where the dimensional transmutation formula: /—^±Q2\ — const N 8TT 1 M R exp ( - , 2 8 ^ bgl (A.158) which comes from renormalisation group arguments, has been used. 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 super-potential in the supersymmetric case. These details are mentioned here for completeness 1 I n [9] the one-loop ^- funct ion was used but most of the discussion is formal ly correct keeping the full /3-function. 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)n-1(.r7). Multi-point correlation functions of the operator H are analogously 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: .9 + 27vn' exp [iW(JJ)] = £ [VA exp / dxG2 + i * ™ f dxGG n'ezJ J 327T J +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 " - ™ - hh=ir- (A-ieo) which satisfy the equation: jdxh = (jdxH) , 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) 2Note 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]. 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 fol-lowing set of equations for W(J, J): Qn+l 8Jn+1 Qn+l W\J=j=0 = in jdxdx1...dxn(T{H(x1)...H(xn)H(0)}) = (-A)n J dx(H), W\J=J=0 = inJdxdx1...dxn(T{H(x1)...H(xn)H(0)}) = (-4)nJdx(H), dJn+1 Qk+l W\J=J-=0 = 0. (A.163) dJkdJl The solution to these equations is: W(J, J) = j dx(H)e~iJ - ^ J dx(H)e-iJ, (A.164) which substituted into (A. 160) and solving for J, J gives: J=-HM) • J=-lMw>- (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): TT ,8U TdU 1,, T x u - h T h - h m = - l ( h + h)- (A-166» One solution to this equation is: U^K'h) = \h\og^ + I f c i o g i + 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: U2(h,h) = ~h\og(^\ + - 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: U3{h, h) = \h Log^ + ^h L o g | + | (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]: e-ivwM) = g g 1* L L o g h + FI L O G J L | . T /, q 9 + 2nn\ h-h\ ,. where the constants CYM,CYM 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 Y M = -2eE^M\Q) = -2e(-bYMas/(32Tr)G2), 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 moti-vated 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(NC) and q = 0{NC °) = 0(1). 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 with Nf light flavours and Nc colours. In the effective Lagrangian approach the l ight 3 matter fields are described by the unitary matrix, U^, corresponding to the phases of the chiral condensate: {¥LVR) = -\{tyLVR)\Ulj, (A.171) w i t h 4 : U = exp r- TTa\a 2 rf i\/2 — M - , UU+ = 1 , (A.172) / T T ^JNf fv' where A a are the Gel l -Mann matrices of SU(Nf), i\a is the octet of pseudoscalar fields (pions, kaons and the eta meson), n' is the SU(Nf) singlet pseudoscalar field (meson) and Nf is the number of quark flavours. We use the values /„. = 132 MeV[22] and fv> = 86 MeV[20] for the meson decay constants 5. The indices, i and j , of U^ run over 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 Tr refers to the trace of the matrix. In the sense of anomalous conformal Ward Identities [20], Q C D reduces to pure Yang-Mi l l s theory with the substitutions (G2)QCD —> {G2)YM and b = bQCD —> bYM- Analogously, an effective Lagrangian for Q C D should transform to that of pure Yang-Mil ls 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 o g h + £ L o g h \ . __ (, q 9 - ilogdet U+ 2irn\ h - h i TrrT, /lt rrr , J + ITTV f k + 1 2 _ j _ _ + - y T r ( M U + h.c.) \ , (A.174) 3Note that the T?' is not really very light, but it enters the theory in this way. 4Note that mixing of the flavor eigenstates is ignored at this level. 5It should be noted that what we refer to as fw/V2 is sometimes denoted fn in the literature. As well, what we refer to as fni, is fn'/VH in the notation of [20]. Appendix A. Derivation of the Effective Potential for QCD 90 where M = diac/(m i | ( t f i t f i ) | ) , E = (bas/(32ir)G2), V is the 4-volume and the inte-gers p and q are relatively prime and related to the parameter £ from (A.155) (A.156) by q/p = 2£. The physical inputs to this equation are the values of the vacuum con-densates {as/n G2) = 0.012 G e V 4 , ( tf"tf u ) = ( t f d t f d ) = - (240 M e V ) 3 and ( t f s t f 5 ) = —0.8 (240 M e V ) 3 , and the quark masses (mu, md and ms). As mentioned above, the sum over k prevents an ambiguity due to the multi-valued nature of the log function by summing over all branches of the effective potential[23][8]. The sum over n is the summation over all integer topological charges. The values of the parameters p and q are not definitely known as different proposals for their determining lead to different values which we wi l l discuss in the next section. Wherever possible we wi l l leave the values of p and q arbitrary but we should note that specific results in Chapter 3 use the the values p = l l i V c — 2Nf and q = 8 and the whole basis of the false vacuum decay in that chapter requires q / 1. The numerical results in Chapter 4 use the values p = Nc and q = 1 but the theory of that chapter is valid for arbitrary p and q. Integrating out the heavy "glueball" fields, h and h, by minimising the effective potential with respect to these variables one obtains the effective potential for the light chiral fields U and U 6: which is Eq.(2.29) and the starting point of al l the research presented in 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 in the effective potential. 6Note that the summation variable k in (A. 174) has been renamed / in this formula. exp VE cos — P P -{8-i\ogDet U) + —l Appendix A. Derivation of the Effective Potential for QCD 91 However, the requirement of Veneziano [19] for the resolution of the U(l) problem: d2kEvac(9) I 2fc / -! x 2fc / 9=0 " i=l X1*c' j dxi(GG(x{)... GG{x2k)) ~ ( — ) (A.176) d92k supports the cos(9/Nc) behaviour of the potential. The next section wi l l discuss the values of the integer parameters p and q. A . 2 Determining the values of p and q In the previous section we saw the anomalous effective potential depends on the integer parameters p = O(N0) and q = 0(NC °) = 0(1). There are different approaches to the fixing of these parameters that lead to very different vacuum structure. I wi l l discuss the physical implications of these parameters and mention some of the different proposals for fixing their values. The parameter lp' is the number of branches of the effective potential. It also cor-responds to the number of distinct vacua before the infinite volume l imit is taken. The parameter lq' determines the number of distinct vacua with different energies that exist before the infinite volume l imit 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 in energy. The parameter 'g' is also the number of physically distinct vacua that exist at 9 = 0 in the infinite volume l imit . It can be shown that the vacuum energy of the 0-vacuum has the following form: hrv hcv E v { 9 ) = { 6 l ~ 3 ^ G 2 l 9 ) = ( Q | " 327 ° 2 l 0 ) C O S { m = ~ m C ° s { m > ( A - 1 7 7 ) where the £ is a parameter from the following low energy theorem: l im z Jdxe^(0\T{^GG(x) ^GG(0)} |0> = £ 2 < ^ G 2 > • (A.178) Appendix A. Derivation of the Effective Potential for QCD 92 Comparing equation (A.177) wi th (2.31) for fa = 0 and 1 = 0 we identify 2£ = q/p. Fur-ther, Equation (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. The first method of fixing the values of p and q comes from a self duality hypothesis in [20] where it was suggested that gluodynamics should be a holomorphic in the composite operators: G2 ± iGG. Comparing with (A.155) and using fi(as) = -bYMOs2/(27r), the Ge l l -Mann Low /^-function for gluodynamics with bYM = ( l l / 3 ) i V c , we.find that this corresponds to q/p = 2£ = 12/11TVC. The second approach is based on the analysis of softly broken supersymmetric theories by counting the number of degenerate vacua. A detailed discussion in the case of pure gluodynamics is given in [12]. The gluino condensate is given by: (AA) = exp(i9/Nc + 2-Kk/Nc), k = 0 , . . . , Nc — 1; this corresponds to Nc degenerate vacua. This sets the value of 'p' at Nc. The vacua are al l degenerate which sets the value of V at 1. The problem with this approach is that it is only valid when the gluon mass, mg <C ASYM where A.SYM is the dynamical mass scale; pure gluodynamics corresponds to the opposite l imit , mg 3> A.SYM- We should also mention that a non-standard non-soft supersymmetric breaking approach [62] leads to the value p = 11 Nc. Another approach was motivated by a suggestion in [63] that in Q C D with massless quarks, nonperturbative matrix elements should be holomorphic in the Paul i Vi l lars fermion mass, MR. Assuming this holomorphy, they relate the proton matrix element of the topological density (p\GG\p) wi th the corresponding matr ix element for the gluon scalar current (p\G2\p). A similar approach was used in [21] for pure gluodynamics to relate the zero momentum two-point functions of GG and G2. Pure gluodynamics was considered as the low energy l imit of a theory including a heavy quark and holomorphy in the physical quark mass was argued to hold in the mq —> oo l imi t . Holomorphy requires the values p = 3bYM = H i V c and q = 8. Appendix A. Derivation of the Effective Potential for QCD 93 A n inverse route was used in [64] to determine 'p' and lq\ This 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. The difference is that instead of starting from the coupled theory, in this approach the starting point is the holomorphic effective poten-t ia 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]. The consistency condition used in [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 = 3bYM = H-ZVc, and q = 8. The agreement of this method with the previous approach is encouraging. We do not pretend that the above discussion fully describes the methods or even lists al l the possible ways of setting the values of the parameters lp' and lq'. Our goal was simply to illustrate the uncertainty about the values and give the two main 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 wi th and without a heavy fermion lead to the values p = 3bYM = H i V c and q = 8 or 12. For Q C D the corresponding values are then p = 3bQCD = HNC — 2Nj and q = 8 or 12. The soft S U S Y breaking approach leads to the values p = Nc and q = 1. As 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. A p p e n d i x B Hyperspher ica l Harmonics i n Four Dimensions Hyperspherical coordinates in 4D dimensions are related to the Cartesian coordinates by: x = r sinO/0 sin(0) sin(<£), (B.179) y = r s in(^) sin(0) cos(<£), (B.180) z = r sin(^)cos(0), (B.181) w = r cos(V>). (B.182) The Laplacian in 4D in hyperspherical coordinates is: 1 1 1 1 —dr (r3dr) + . 2 (s in 2 ^ ) + . 2 dg (sin0<9<,) + . 2 , 2 r 3 v 7 r 2 s i n ^ v ' r 2 sin z ?/> sin 0 r 2 s i n ^ s m 0 9 (B.183) Assuming separable solutions and treating 0 and <f> coordinates in exactly the same way as in three dimensions we obtain the differential equation: V"{ib) + 2 cotO/O '^O/O - 1(1 + 1) csc20/0$(VO = A * ( B . 1 8 4 ) W i t h the substitution u = cos(ip) this becomes: (1 - u2)U"(u) - 3uU'(u) - l(; l + l}u(u) = BU(u). (B.185) 1 — ul 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 Geigen-bauer (Ultraspherical) equation: (1 - x2)-^Cia\x) - (2a + l ) ^ ) ( i ) - 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: Ylm{9A) Unl(cos(xp)) 0<l<n , (B.187) [ Y\i\m{6,<f>) Unl-l{cos{ibr)) - n < l < - l where YiM(6,(p) are the usual 3D spherical harmonics and Unl are associated Chebyshev type I I functions defined by: Ynim(9, <j>, ip) = A(n, I, m) I Un l(x) = (1 - x*y/^Un(x) = (1 - x 2 ) ' / 2 ;-l)"(n + l )V i dl < 1-*V / ,-^{<i-* ! ,>" + , / s} dx 0.188) 2"+1(n +1/2)! dx1 for I > 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 ^ ) = ^ ^^'P^-^), (B.190) n=0 l=—n m=—l sin tp sin 6 ROC R-K RN r2-K JO Jo JO JO Y^rnir^i(t)^)Yn'Vm'{r^A,^) = Snn>SW5W. (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
pdf
Page Metadata
Item Metadata
Title | Vacuum structure of QCD in an effective Langrangian approach |
Creator |
Fugleberg, Todd Darwin |
Date Issued | 2000 |
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 |
FileFormat | application/pdf |
Language | eng |
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 |
GraduationDate | 2000-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | 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:
https://iiif.library.ubc.ca/presentation/dsp.831.1-0085708/manifest