Topics in Theory and Experiment in Relativistic Heavy-Ion Physics
A dissertation submitted to Kent State University in partial fulfillment of the requirements for the degree of Doctor of Philosophy
by Jeremy Alford May 2015
c Copyright All rights reserved Except for previously published materials Dissertation written by Jeremy Alford BS, Rochester Institute of Technology, 2004 MA, Kent State University, 2009 Ph.D., Kent State University, 2015
Approved by
Michael Strickland, Associate Professor, Ph.D., Department of Physics, Doctoral Advisor Declan Keane, Professor, Ph.D., Department of Physics, Doctoral Advisor Veronica Dexheimer, Assistant Professor, Ph.D., Department of Physics Antal Jakli, Professor, Ph.D., Liquid Crystal Institute John L. West, Professor, Ph.D., Department of Chemical Physics
Accepted by
James Gleeson, Chair, Ph.D., Department of Physics James L. Blank, Dean, Ph.D., College of Arts and Sciences Contents
TABLE OF CONTENTS ...... iii
LIST OF FIGURES ...... vi
LIST OF TABLES...... x
ACKNOWLEDGEMENTS ...... xi
1 Introduction 1 1.1 Heavy-Ion Collisions ...... 1 1.2 Structure of Dissertation ...... 2 1.3 Natural Units ...... 2
2 Quantum Chromodynamics 4 2.1 Confinement ...... 6 2.2 Running Coupling ...... 7 2.3 Lagrangian ...... 7 2.4 Finite Temperature Effects ...... 9
3 Strong Magnetic Field Effects on Quarkonia 13
4 Pseudomomentum 15 4.1 Two Coupled Particles in a Constant Magnetic Field ...... 17 4.2 Two Particles With Equal and Opposite Charge ...... 19 4.3 Heavy-Light System ...... 20 4.4 Particle-Antiparticle Pair ...... 21 4.5 Pseudomomentum and Kinetic COM Momentum ...... 21 4.6 The Motional Stark Effect ...... 22
iii 5 Heavy Quark Model 24 5.1 Heavy Quark Effective Theory ...... 24 5.2 Heavy Quark Lagrangian ...... 25 5.3 Heavy Quark Hamiltonian ...... 32
6 Harmonic Oscillator 38 6.1 Analytic Solution ...... 38 6.2 Sudden Approximation ...... 41
7 Numerical Solution 49 7.1 Solving the Schr¨odinger Equation ...... 51 7.2 Finding the Ground State ...... 52 7.3 Finding the Excited States ...... 52 7.4 Center-of-Mass Kinetic Energy Subtraction ...... 53
8 Vacuum Results 55 8.1 Bottomonia ...... 55 8.2 Charmonia ...... 55
9 Spin-Mixing 58
10 Final Results 64 10.1 Bottomonia ...... 66 10.2 Charmonia ...... 71 10.3 Discussion ...... 76 10.4 What is a Particle? ...... 77
11 Search for Three-body Hypertriton Decay 80 11.1 History of the Hypertriton ...... 80 11.2 Decay Channels ...... 81
iv 11.3 Experimental Search ...... 81 11.4 Previous Measurements ...... 82
12 The STAR Experiment 84 12.1 Time Projection Chamber ...... 85 12.2 Time of Flight ...... 86
13 Identification of Hypertriton Candidates 87 13.1 Background Subtraction ...... 88
14 Locating the Decay Vertex 90
15 Particle Identification 93 15.1 Deuteron Selection ...... 93 15.2 Time-of-Flight Correction ...... 95
16 Background Suppression 99 16.1 TMVA ...... 100 16.2 Simulation ...... 101
17 Yields 105
18 Other Measurements Considered 109 18.1 Embedding ...... 110
19 Summary 112
References 114
v LIST OF FIGURES
Figure 1. Schematic of confinement. The number of lines between the quark (q) and antiquark (¯q) represent the strength of the attractive force. As theq ¯ are forced to separate the force becomes stronger until there is enough energy to create a new qq¯ pair...... 7
Figure 2. Running coupling of the strong force...... 8
Figure 3. Debye screening in classical electrodynamics. Positively charged particles (red) will attract negatively charged particles (blue) from the medium reducing the force between the positively charged particles. The dashed circles represent the Debye screening length ...... 10
Figure 4. The Cornell potential in vacuum and finite temperature. In QCD with 6 flavors,
7g3 β is related to the coupling constant gs by β = − 16π2 . Figure taken from [1] . . 11
Figure 5. a) Rest masses of 1s bottomonium states. b) Rest masses of 1p bottomonium states...... 65
Figure 6. a) Rest masses of 1s charmonia states. b) Rest masses of 1p charmonia states.66
Figure 7. a) Rest masses of the spin-mixed 1s bottomonium states. b) Rest masses of the unmixed degenerate spin-triplet 1s bottomonium states ...... 68
Figure 8. Rest masses of the non-degenerate spin-mixed 1p bottomonium states...... 68
Figure 9. Rest masses of the degenerate spin-mixed 1p bottomonium states ...... 69
Figure 10. a) Rest masses of the unmixed degenerate 1s bottomonium states. b) Rest masses of the unmixed non-degenerate 1p bottomonium states...... 70
Figure 11. a) Rest masses of the spin-mixed 1s charmonium states. b) Rest masses of the unmixed degenerate spin-triplet 1s charmonium states ...... 73
vi Figure 12. Rest masses of the non-degenerate spin-mixed 1p charmonium states...... 74
Figure 13. Rest masses of the degenerate spin-mixed 1p charmonium states ...... 74
Figure 14. a) Rest masses of the unmixed degenerate 1s charmonium states. b) Rest masses of the unmixed non-degenerate 1p charmonium states...... 75
Figure 15. Left: Decay length measurements of the lambda hyperon and hypertriton. The insert shows the χ2 analysis used to obtain the lifetime measurement. Right: Measured lifetimes and theory calculation of the hypertriton. Figure taken from [2]...... 83
Figure 16. Energy loss per centimeter in the TPC vs. particle momenta. The upper band are deuterons, the middle band are protons, and the bottom band are pions. . .85
Figure 17. Schematic of rotation technique used to estimate the combinatorial background.
The transverse momentum (pt) vector and decay vertex position (rt) vectors for a daughter track are rotated about the primary vertex (PV) in the plane transverse to the beam direction (arbitrary scale). Different versions of this technique may use different daughter tracks and different rotation angles...... 89
Figure 18. Left: Hypertriton invariant mass distribution (red) with combinatorial background (blue) estimated by rotating the pion tracks. Right: Hypertriton invariant mass distribution (red) with combinatorial background (blue) estimated by rotating the deuteron tracks...... 89
Figure 19. Schematic of Λ → p+π− decay. The decay is assumed to take place at the midpoint of the line segment representing the distance of closest approach between the daughter particles. A two-body hypertriton decay may be obtained by replacing the proton with a helium-3 nucleus ...... 91
vii Figure 20. Schematic of 3-body hypertriton decay. The decay is assumed to take place at the center of a triangle formed by the midpoints of the line segments representing the distances of closest approach between each pair of daughter particles...... 91
Figure 21. 2D histograms showing deuteron selection for the three matter data sets used. The red box indicates the locations of the track cuts. The antimatter data from the 2011 run at 200 GeV is similar to the corresponding matter data set but with tighter cuts on the horizontal axis. See text for exact cuts ...... 96
Figure 22. Masses of daughter particles with (blue) and without (red) the TOF correction97
Figure 23. Reconstructed hypertriton mass spectrum (red) and estimated combinatorial back- ground (blue) using the TOF correction. This plot was made using the 2011 run, 200 GeV (matter) data set ...... 98
Figure 24. Left: Definition of the pointing angle θ. Right: Distribution of θ. Hypertriton candidates with a large pointing angle cannot be real ...... 100
Figure 25. Transverse momentum vs. the distance of closest approach to the primary vertex of the daughter particles for real and simulated data...... 102
Figure 26. Distance of closest approach between pairs of daughters vs. the reconstructed hypertriton transverse momentum for real and simulated data. Due to limited data storage, no candidate with a DCA larger than 1 cm was saved ...... 102
Figure 27. Distances of closest approach between different daughter particles plotted against each other for real and simulated data. Here we see a weak correlation in the simulated data between the inter-particle DCAs involving the pion but it is not strong enough to define a new cut without significantly reducing the number of accepted hypertriton candidates ...... 103
viii Figure 28. Distance between the decay vertex and the midpoint of the inter-particle DCA vs. the reconstructed hypertriton transverse momentum for real and simulated data...... 103
Figure 29. Top: Distribution of real and simulated deuteron momentum transverse to hyper- triton momentum. Middle: Distribution of real and simulated proton momentum transverse to hypertriton momentum. Bottom: Distribution of real and simulated pion momentum transverse to hypertriton momentum. All data sets have similar distributions for these three quantities. The red vertical lines indicate the locations of the cuts ...... 104
Figure 30. Background-subtracted reconstructed hypertriton invariant mass spectrum with a bin-width of 4 MeV/c2. The blue arrow shows the location of the hypertriton mass (2.991 GeV/c2). There is a narrow peak in the data at this point for the 2011 run data (bottom row). For the 2010 run data (top row), there is a narrow peak located 4 MeV/c2 higher than the known hypertriton mass ...... 106
Figure 31. Background-subtracted reconstructed hypertriton invariant mass spectrum with a bin-width of 4/3 MeV/c2. The blue arrow shows the location of the hypertriton mass (2.991 GeV/c2). There is a narrow peak in the data located 4/3 MeV/c2 higher than this point for the 2011 run data (bottom row) and 8/3 MeV/c2 higher than this point for the 2010 run data (top row) ...... 107
Figure 32. Reconstructed lambda invariant mass spectrum (with combinatorial background) for the decay Λ → p + π−. The large peak is centered on the lambda hyperon mass of 1.116 GeV/c2 ...... 108
ix LIST OF TABLES
Table 1. Unit Conversions ...... 3
Table 2. Expectation Values...... 50
Table 3. Vacuum Bottomonia Results ...... 56
Table 4. Vacuum Charmonia Results ...... 57
Table 5. Vacuum Eigenstates ...... 59
Table 6. Modified Eigenstates ...... 60
Table 7. Bottomonia Eigenstates with Moderate Magnetic Field ...... 67
Table 8. Bottomonia Eigenstates with Strong Magnetic Field ...... 67
Table 9. Charmonia Eigenstates with Moderate Magnetic Field ...... 72
Table 10. Charmonia Eigenstates with Strong Magnetic Field...... 72
Table 11. Eigenstates with Infinite Magnetic Field ...... 78
Table 12. Data Statistics ...... 84
Table 13. Decay Vertex Cuts ...... 92
x Acknowledgements
First, I would like to thank my advisors, Michael Strickland and Declan Keane for their guidance and support. Their generosity and understanding made it possible for me to pursue my education without sacrificing my family life. Although he was never officially my research advisor, Spiros Margetis also provided me with much useful advice. Thanks are also due to the many people who helped me with this research. Specifically, I would like to mention Amilkar Quintero, who did most of the TMVA work; Jinhui Chen, who studied the two-body hypertriton decay and provided me with many valuable suggestions; Yuhui Zhu, who also studied the three-body hypertriton decay and provided useful discussions; Jonathan Bouchet for his help with the hypertriton simulation; and Patrick Huck for his help with the time- of-flight correction. Finally, I would like to thank my family for their enduring faith in me during this long journey, especially my wife, Alana Criswell, and son, Bowen Criswell.
xi 1 Introduction
1.1 Heavy-Ion Collisions
A physics professor once said, “The most primitive way to learn how a locomotive works is to force two of them to crash at high speed and examine the pieces.” This is the basic idea behind relativistic heavy-ion collisions. We learn about nuclear matter by accelerating nuclei to near the speed of light, forcing them to collide, then studying the particles produced. The fundamental theory describing the interaction between nucleons is quantum chromodynamics (QCD). It is a very complicated theory making it very difficult to do most calculations from first principles. Instead, we must rely on lattice calculations and phenomenological models to understand heavy-ion collisions. The hot dense matter formed in these collisions is believed to be very similar to the matter that filled the universe immediately after the big bang. The energy involved in the collision is enough to create new particles that were not present in the colliding nuclei. The newly created particles may be matter or antimatter. Many of the particles produced are never observed directly because they are unstable and quickly decay into lighter particles. In addition to producing many exotic particles, heavy-ion collisions also provide a medium to test predictions of QCD. Because the nuclei are moving at relativistic speeds, Lorentz contraction turns approxi- mately spherical nuclei into an oblate spheroid in the lab frame. Most of the collisions are not head-on, so some of the nucleons, called spectators, are not involved in the collision. At high energy the participant nucleons break into their constituent partons and form a quark- gluon plasma. As the plasma expands and cools, particles are released in all directions.
1 Large particle detectors are placed near the point of the collision. The information from the released particles is used to reconstruct the collision event. In addition to producing particles from the energy of the collision, a short-lived but very intense magnetic field is expected to be generated in the quark-gluon plasma. To fully understand particle production, one must account for these magnetic effects.
1.2 Structure of Dissertation
Over the course of my research I worked on several projects, some experimental and some theoretical. The first part of this dissertation describes a theoretical exploration of the effects of very strong magnetic fields on quarkonia, specifically, bottomonia and charmonia. After the concept of pseudomomentum is introduced and the model Hamiltonian is derived, the results of a numerical calculation with no background magnetic field are presented and com- pared to experimental results. The same numerical algorithm is used to calculate quarkonia energy levels for various values of magnetic field strength neglecting the interaction between the particle spin and the background magnetic field. The spin-field interaction is accounted for and the final results are presented. The second part describes an experimental search for a particular three-body decay of hypertriton particles, the nuclei of an unstable heavy isotope of hydrogen. Hypertritons are produced in heavy-ion collisions but decay into smaller particles so they cannot be observed directly. Instead we must look for the decay products and reconstruct the original parent particle. If enough hypertritons are found, the expected lifetime can be calculated and compared to previous calculations.
1.3 Natural Units
Before making any measurement, one must first choose an appropriate set of units. There is a set of units commonly used by physicists which is called “natural” because it simplifies many calculations by setting fundamental constants, such as the speed of light and Planck’s
2 constant, equal to one. This may seem impossible since, in SI units, the speed of light is a very large number and Planck’s constant is a very small number but the conversion is no more difficult than converting to any other set of units. The conversion factors are listed in Table 1. In natural units, everything has units of energy raised to an integer power. The factors of the speed of light and Planck’s constant are not always written but can easily be replaced using dimensional analysis.
Unit Conversions SI units Natural units Energy 1.6022 × 10−10 J 1 GeV Mass 1.7827 × 10−27 kg 1 GeV/c2 Time 6.5822 × 10−25 s 1 GeV−1~ Distance 1.9733 × 10−16 m 1 GeV−1~c Magnetic Field 1.6904 × 1020 Gauss 1 GeV2(e~c2)−1
Table 1: Conversions between SI and natural units. Formally, the natural units should include factors of Planck’s constant and the speed of light but, in practice, they are usually not written.
3 2 Quantum Chromodynamics
There are many similarities between QCD and the theory of quantum electrodynamics (QED). Instead of two degrees of freedom (positive and negative charges) like QED, QCD has three degrees of freedom, called color charges, which are arbitrarily labeled red, blue, and green. Like QED, the force between different charges is attractive, the force between like charges is repulsive and when all degrees of freedom are combined the overall charge is neutral. Antimatter particles have the opposite electric charge relative to their matter partners. The color charges for antimatter are antired, antiblue and antigreen. When two different color charges are combined, they form the antimatter version of the third color charge. When a color charge is combined with an anticolor charge of the same color, the result is color neutral; this is the case for quarkonia. Under normal circumstances (low temperature, low pressure, ...), quarks always combine in ways that are color neutral. There are only two types of color neutral particles composed of quarks: mesons and baryons. Mesons are bound states of a quark and an antiquark and baryons are bound states of three quarks of different colors. The strong force between quarks is transmitted by gauge bosons called gluons which are analogous to photons from QED. Another important difference between QCD and QED, known as confinement, is the fact that, unlike electric charges and photons, free quarks and gluons are not found in nature. QCD was first proposed in 1963 by Gell-Mann and Zweig [4]. By this time, QED had already been established and there was much experimental data to be explained. The light mesons had the right quantum numbers to be bound states of quarks and antiquarks and baryons were interpreted as bound states of three quarks. The correct electric charges
4 and other quantum numbers were obtained by assuming three flavors of quarks (up, down, strange). The others (charm, bottom, top) were discovered later. Theoretically, there could be many more but they would have to be more massive than the top quark or they would have already been discovered. Nature has already pushed the limit of weirdness with the existence of the top quark, which is about 175 times more massive than the proton, so it seems unlikely that heavier quarks will ever be discovered. To get the correct electric charges of baryons the up, charm and top quarks must have a charge of +2/3e and the down, strange and bottom quarks have a charge of −1/3e, where e is the elementary charge. Now the nucleon charges can be easily explained if the proton is made up of two up quarks and a down quark and the neutron is two down quarks and an up quark. Quarks have spin 1/2, so they must obey Fermi-Dirac statistics but to satisfy the observed baryon spectrum the quark wavefunction must be symmetric under the interchange of spin and flavor quantum numbers. For example, the ∆++ is a spin-3/2 nucleon with charge +2. This particle can be easily understood as a bound state of three up quarks with their spins parallel and no orbital angular momentum but without another quantum number it would not obey the proper spin statistics. The additional quantum number is provided by the color charge. Since there are three quarks in a baryon and three color charges, the baryon wavefunction can be totally symmetric in spin and flavor but totally antisymmetric in color making it totally antisymmetric overall and satisfying the Fermi-Dirac statistics. This is how the color charge was originally identified as the proper degree of freedom to describe the strong nuclear force. Mathematically, we say QED belongs to the group U(1) and there is one type of photon that forms the generator for the group. QCD belongs to the group SU(3) with eight types
2 of gluons as the generators. The number of generators is given by Nc − 1, where Nc is the number of different charges in the theory. QED is an Abelian gauge theory and QCD is a non-Abelian gauge theory. They both have a vector potential and different components of
5 the potential commute in QED but not in QCD.
2.1 Confinement
Confinement is the physical manifestation of the non-Abelian nature of QCD. It is a result of gluons having a color charge causing them to couple to each other and quarks. Confinement does not happen in QED because photons do not have an electric charge. Confinement can be understood qualitatively by imagining the quarks are connected by dense field lines that form a string. The string is described as a color flux tube with constant cross-sectional area and constant energy density. As the quarks are separated, the force holding them together becomes stronger and it becomes energetically favorable to break the string and create a new pair of quarks, see Figure 1. Mathematically, confinement can be understood in terms of the Cornell potential,
4 α V = − s + σr , (1) Cornell 3 r
where αs is the strong coupling constant, σ is the string tension and r is the distance between the quarks. The first term is like the Coulomb potential from classical electrodynamics with a different coupling strength 1. The second term comes from the non-Abelian part of QCD and is responsible for confinement. It can be understood as the potential energy stored in the string which, due to the constant energy density, increases linearly with length. From Eq. (1) it is easy to see that, when the quark separation is small, the Coulomb part dominates and the theory becomes almost identical to QED but, for large quark separation, the linear part dominates and the quarks are never able to become free. Despite much effort, a rigorous mathematical proof of confinement does not exist [5]. However, simple potential models like the one above describe the observed particle spectra very well.
1The factor of 4/3 is called the color factor and is a result of the existence of six different color charges (including anti-colors). The color factor is 4/3 for color-singlet states and −1/6 for color-octet states.
6 Figure 1: Schematic of confinement. The number of lines between the quark (q) and anti- quark (¯q) represent the strength of the attractive force. As the qq¯ are forced to separate the force becomes stronger until there is enough energy to create a new qq¯ pair.
2.2 Running Coupling
The coupling constant in QCD, αs, is not constant. As the distance between interacting quarks decreases, so does the strength of the coupling constant. Since probing the interaction at shorter distances requires higher energy it is also correct to say that the strength of the coupling constant decreases with increasing energy. This is known as asymptotic freedom and is a prediction of perturbative QCD. Figure 2 shows the remarkable agreement between theoretical prediction and experimental results. Asymptotic freedom is also a result of the non-Abelian nature of QCD. In fact, any renomalizable quantum field theory must have non-Abelian gauge fields to be asymptotically free [6].
2.3 Lagrangian
The QCD Lagrangian satisfies the charge and parity symmetries and quark flavor conserva- tion. Although QCD is much more complicated than QED, its Lagrangian can be written in a way that looks like the QED Lagrangian,
1 L = Ψ(¯ iγ Dµ − m)Ψ − F F µν , (2) QCD µ 4 µν
7 Sept. 2013 decays (N3LO) s(Q) Lattice QCD (NNLO) DIS jets (NLO) 0.3 Heavy Quarkonia (NLO) – e+e jets & shapes (res. NNLO) Z pole fit ( N3LO) (–) pp –> jets (NLO) 0.2
0.1 QCD s(Mz) = 0.1185 ± 0.0006 1 10 100 1000 Q [GeV] Figure 2: Running coupling of the strong force. with the covariant derivative defined as
a Dµ = ∂µ − igτaAµ , (3) and the gluon field strength tensor defined as