arxiv:hep-ph/ v3 10 Sep 2002

Transcription

1 Formation and Collapse of False Vacuum Bubbles in Relativistic Heavy-Ion Collisions Rajarshi Ray, Soma Sanyal, and Ajit M. Srivastava Institute of Physics, Sachivalaya Marg, Bhubaneswar , India arxiv:hep-ph/ v3 10 Sep 2002 Abstract It is possible that under certain situations, in a relativistic heavy-ion collision, partons may expand out forming a shell like structure. We analyze the process of hadronization in such a picture for the case when the quark-hadron transition is of first order, and argue that the inside region of such a shell must correspond to a supercooled (to T = 0) deconfined vacuum. Hadrons from that region escape out, leaving a bubble of pure deconfined vacuum with large vacuum energy. This bubble undergoes relativistic collapse, with highly Lorentz contracted bubble walls, and may concentrate the entire energy into extremely small regions. Eventually different portions of bubble wall collide, with the energy being released in the form of particle production. Thermalization of this system can lead to very high temperatures. With a reasonably conservative set of parameters, at LHC, the temperature of the hot spot can reach as high as 3 GeV, and well above it with more optimistic parameters. Such a hot spot can leave signals like large P T partons, dileptons, and enhanced production of heavy quarks. We also briefly discuss a speculative possibility where the electroweak symmetry may get restored in the highly dense region resulting from the decay of the bubble wall via the phenomenon of non-thermal symmetry restoration (which is usually employed in models of pre-heating after inflation). If that could happen then the possibility may arise of observing sphaleron induced baryon number violation in relativistic heavy-ion collisions. PACS numbers: q, Mh, Cq Typeset using REVTEX rajarshi@iopb.res.in sanyal@iopb.res.in ajit@iopb.res.in 1

2 I. INTRODUCTION There are strong reasons to believe that in ultra-relativistic collisions of heavy nuclei, a hot dense region of quark-gluon plasma (QGP) may get created. There is a wealth of data which strongly suggests that already at CERN SPS this transient QGP state may have been achieved [1]. With new data already coming out from RHIC at BNL [2], it may be only a matter of time that conclusive evidence of QGP detection would be obtained. Certainly, one expects that at LHC (in next few years) QGP will be produced routinely. There is no question that detection of this new phase of matter will be of utmost importance, not just for testing predictions of QCD, but also in providing us glimpses of how our universe may have looked like at few microseconds of age. It seems to be an appropriate stage that one starts looking beyond the detection of the QGP phase. Many interesting possibilities have been discussed in the literature about physical processes which may become observable once QGP is produced in laboratory. For example, one will have the opportunity of studying phase transitions in a relativistic field theory system under controlled laboratory situations. The richness of QCD phase diagram may become available for probing with quark-hadron transitions occurring in these experiments [3]. In this paper we propose a novel possibility of forming a pure false vacuum bubble in a relativistic heavy-ion collision, whose collapse and decay can lead to extremely hot, tiny regions, with temperatures well above the initial temperature of the QGP system. We assume that the confinement-deconfinement phase transition is of first order, and show that such a false vacuum bubble will form if the parton system (formed in the collision) expands out in a shell like structure, leaving behind the supercooled deconfined phase with large vacuum energy density. Eventually, as the partons hadronize and escape out, the interface, separating this false vacuum bubble from the confined phase outside, undergoes relativistic collapse. Due to absence of any plasma inside, the motion of this interface becomes ultra-relativistic, with highly Lorentz contracted interface thickness. Entire energy of the false vacuum bubble gets converted to the kinetic energy of the wall. Eventually different portions of bubble wall collide, converting the entire energy into particles. These particles may thermalize and lead to an extremely hot, very tiny region (apart from possible effects of quantum fluctuations of the bubble wall, as we discuss later). We find that even at RHIC energies, the temperature at these hot spots can be more than 1 GeV. For LHC the temperature can reach several GeVs. This will have important signals such as increased production of heavy quarks, anomalously large values of P T of some hadrons (those coming from the hot spot) etc. Also, the net energy of the false vacuum bubble may be very large even at RHIC ( 1 TeV), hence the possibility of Higgs and top quark production may also arise (depending on the decay products of the bubble wall). At LHC energies, the net energy of this false vacuum bubble is very large which should lead to very dense parton system after the bubble wall decay. Due to initial very large density of partons at such a spot (forming a non-equilibrium system) we speculate on the possibility that electroweak symmetry may get restored there via the phenomenon of nonthermal symmetry restoration, as in the models of pre-heating after inflation[4]. If that could happen then it raises the possibility that sphaleron transitions may occur in these regions (depending on the size of such a region), which will lead to baryon number non-conservation. It is needless to say that any possibility of detection of electroweak baryon number violation 2

3 in laboratory experiments deserves attention, especially with its implications for the theories of baryogenesis in the early universe. We mention here that the possibility of baryon number violation in collider experiments has been discussed earlier, see, ref.[5] and references therein. However, the discussion in these works was about possibility of baryon number violation at high energies, and not due to electroweak symmetry restoration. It is fair to say that at this stage it is not clear whether it is possible to get baryon number violating interactions at high energies. However, baryon number violation in the electroweak symmetric phase is on rather strong foundations [6,7,8], and our reference is to this type of baryon number violation. The paper is organized as follows. In section II, the basic physical picture of the model is discussed where it is argued how a shell enclosing deconfined vacuum may form in ultrarelativistic heavy-ion collisions. Section III discusses the properties of the shell, the surface energy, volume energy etc. Section IV discusses the evolution of this shell and its final ultra-relativistic collapse leading to high concentration of energy in a tiny region. In section V we discuss the issue of thermalization of the decay products of bubble wall. Section VI presents results for various ranges of parameters for LHC and RHIC energies, as well as for the possibility of non-thermal restoration of the electroweak symmetry. Discussion of results and the conclusion is presented in section VII. II. PHYSICAL PICTURE OF THE MODEL The main aspect of our model is based on the observation[9] that under certain situations, the expanding parton system may form a shell like structure. This type of picture emerges under a variety of conditions. For example, in a hydrodynamical expansion, the rarefaction wave reaching center gives rise to a shell like structure [9]. (We mention, that shell like structure with a maximum of density at center has also been discussed in the literature [10].) It has been shown that shell structure will arise generically in an expanding parton system when partons are ultra-relativistic and particle collisions are not dominant [9]. This happens for the simple reason that with all partons having velocity c, the partons pile up in a shell of radius c t, with the thickness of the shell being of the order of the size of the initial region. Indeed, the original so called baked Alaska model for the disoriented chiral condensates, proposed by Bjorken et al. [11] utilizes a shell like picture for the expanding parton systems (see also, ref. [12] where a similar picture has been discussed). Our model is based on this generic shell structure of the parton system. However, there is one important difference between our model and other discussions in the literature where a shell like structure of partons has been discussed. In refs. [11,12], the discussion was in the context of a second order chiral phase transition where the larger potential energy of the inside region slowly decreases as the chiral field relaxes towards the true vacuum. In contrast, our model is based on a first order confinement-deconfinement phase transition, where the vacuum energy of the inside region can only decrease by collapse of the interface, or via nucleation of true vacuum bubbles in this region, which we will argue to be very suppressed. Though we have presented the discussion in terms of a first order deconfinement-confinement transition, our entire discussion will also be valid for a first order chiral phase transition. Basic energy scales of these two transitions being roughly of same 3

4 order, even quantitative aspects of the discussion will not change much. In discussions of [9] it is argued that the hadronization in the shell will proceed from the inner boundary of the shell, as well as the outer boundary of the shell. This assumes that the empty region inside the shell will be in the confined vacuum. We argue below that there is no reason to expect that. Rather, one expects that the shell of partons will enclose a region of deconfined vacuum. Let us start with the initial stage of central collision of two nuclei. For very early stages after the collision, the region between the two receding nuclei will be populated by a dense system of partons, with a temperature which is expected to be well above the critical temperature T c of deconfinement-confinement phase transition (at RHIC and LHC). (We will consistently refer to the phase transition as the deconfinement-confinement (D-C) phase transition rather than using the conventional terminology of quark-hadron transition. This is because we will utilize the difference between the vacua of the two phases of QCD, irrespective of the fact whether these vacua are populated by quark-gluon, or hadron degrees of freedom.) Subsequently, plasma will expand longitudinally for early times (for τ < R A, where τ is the proper time and R A is the radius of the nucleus), and will undergo threedimensional expansion for larger times. For heavy nuclei with large A, and for center of mass energies at RHIC and LHC, it is expected that the plasma will undergo deconfinementconfinement phase transition during this three dimensional expansion stage, and later the expanding hadronic system will freezeout. In various models of freezeout, it is argued that the freezeout happens soon after the D-C phase transition and that this stage is achieved at proper time hypersurface which is close to the hypersurface when three-dimensional expansion commences [13] (see, also ref. [14]). The expected value of central energy density ǫ i at the initial stage (hence the initial temperature T i ), increases with A and with s. An estimate of the dependence of ǫ i on A and s can be obtained from the scaling relations proposed in ref. [15] (see also, ref. [16]), ǫ i = 0.103A ( s) GeVfm 3, (1) where s is in GeV. These scaling exponents reproduce the expected values of initial energy density for SPS, RHIC and LHC. Thus, for LHC, with A = 208 and s = 5.5 TeV (we will use s to refer to s per nucleon), we find ǫ i 1.3 TeV/fm 3 which is within the range of other estimates [17]. We note that dependence on A is somewhat weaker than the dependence on s. Thus, if one decreases value of A and increases s suitably, one can still get large ǫ i (suitable for having a QGP state initially at temperatures well above T c ). However, with lower value of A, R A will be smaller (R A 1.1A 1/3 ), implying that the three-dimensional expansion will commence at earlier values of proper time. This will have two effects. Temperature will decrease faster (T τ 1 for three-dimensional expansion while T τ 1/3 for longitudinal expansion [18]). More importantly, with rapidly expanding system, thermodynamic equilibrium will be maintained for a shorter duration of proper time, implying that freezeout will happen earlier. In principle, there is no reason why freezeout cannot precede the D-C phase transition. In such a situation, initial thermalized QGP system expands longitudinally for a short time, then undergoes three dimensional expansion. Due to rapid 3-dimensional expansion, the system falls out of equilibrium and keeps expanding as a parton system in a 4

5 non-equilibrium state (still with high density of partons so that no hadronization takes place yet). Eventually the parton system will be dilute enough to enter into the non-perturbative regime, and will undergo hadronization. This process of hadronization will be in a nonequilibrium state, where a parton system out of equilibrium converts to a hadron system (again, out of equilibrium). Possibility of freezeout preceding hadronization has also been discussed previously in ref. [19]. A quantitative discussion of the limiting values of s, and that of A below which freezeout may precede hadronization, can only be given using elaborate numerical computations (such as those in ref. [13]). In such a picture of expansion of parton system, it is possible to have the central region (in the shell) depleted of partons while the partons get accumulated in the shell. Thus, central region does not cool to become zero temperature, zero density QCD matter, as would happen if the system was always in thermal equilibrium (and that would then require that the interior of the shell be in the confined vacuum, as assumed in [9]). Rather, in the above picture, the interior of the shell is depleted of partons because of rapid expansion of the parton system. As thermal equilibrium is not maintained, parton collision rate lags behind the expansion rate. As mentioned above, in such a situation the analysis of refs. [9] suggests that partons, due to their relativistic velocities, pile up in a shell of radius R c t where t is the time in the center of mass frame of the collision. The thickness of the shell r will be expected to be of the order of the size of the region when thermal equilibrium is first broken such that parton system undergoes almost free expansion after that. We mention here that it is conceivable that the initial parton system never thermalizes (for example for very small A), but nevertheless, deconfinement phase is achieved. This possibility can originate from the discussions of non-thermal symmetry restoration in the context of inflation in the early universe [4]. In these works it is shown that rapid particle production due to parametric resonance can lead to modification of the effective potential, even without any notion of thermal equilibrium. Under certain situations the symmetry restoration can be achieved even when thermal symmetry restoration is not expected. The early stages in nucleus-nucleus collisions resemble the stage of parametric resonance at least in the sense of rapid particle production. It is then possible that in this case also deconfinement phase is achieved with parton system always remaining out of equilibrium. This system will then undergo expansion forming a shell structure [9] before hadronizing eventually. (Later, in Sec.VI, we will discuss whether a similar picture can also emerge for the decay products of bubble wall, leading to non-thermal restoration of the electroweak symmetry even when the system remains in a non-equilibrium state.) Thus, we assume that the expanding parton system forms a shell, with the interior of the shell being in the deconfined vacuum. At the outer boundary of the shell there must be an interface separating the deconfining vacuum in the interior from the confining vacuum in the outer region [20]. (We are assuming the D-C transition to be of first order, with a barrier separating the metastable deconfined vacuum from the true confining vacuum, even at temperatures approaching zero. We will discuss this issue in detail below.) Partons will be piled up inside the shell, within a thickness r of the interface. The center of the shell interior thus represents a supercooled (to T = 0) region trapped in the metastable deconfining vacuum. It is somewhat uncommon to talk about zero temperature deconfining phase of QCD. Normally one associates deconfinement phase with the quark-gluon plasma at temperatures 5

6 above T c. However, it is important to distinguish the vacuum of a theory from the particle excitations about that vacuum. Confining phase and the deconfining phase of QCD should be taken to correspond to the two different vacua of some effective theory of QCD, with the deconfining vacuum being metastable at zero temperature. The two vacua may be characterized by an order parameter, such as the expectation value of the Polyakov line [21]. That means that certain type of gauge field background leads to one vacuum, while a different gauge field background gives the different vacuum. These two vacua are thus defined irrespective of the presence or absence of quarks and gluons, or hadrons, as particle degrees of freedom (except as test particles in order to define confinement). It then makes sense to talk about deconfining vacuum without any quark or gluon being present as particle degrees of freedom, though this vacuum will be metastable. (For example, if the chiral phase transition is of first order, with a chirally symmetric vacuum even at T = 0, then it will make perfect sense to talk about a supercooled chirally symmetric phase with T = 0.) The type of picture we are taking in our model, with the shell interior depleted of quarks and gluons but still trapped in the metastable deconfining vacuum, is very similar to what happens in the inflationary theories of the universe [22]. For example, in the old inflation, the scalar field φ gets trapped in the false vacuum and supercools as the universe undergoes exponential expansion due to non-zero potential energy of φ. Initially the universe is filled with scalar particles (as well as other particles) in thermal equilibrium corresponding to a temperature above the transition temperature. However, due to exponential expansion, initially the system cools maintaining equilibrium, but then falls out of equilibrium with all particles undergoing depletion due to exponential expansion of the universe. Eventually the particles are completely diluted away, while the universe still remains in the metastable vacuum. This is precisely the same picture as we are using in our model. The only thing required for such a picture is that QCD admit a metastable vacuum even at zero temperature. The exponential expansion of the universe is replaced by rapid outward expansion of partons in our case. Possibility of supercooling and a brief inflationary phase at quark-hadron transition in the universe has also been discussed in the literature [23]. III. PHYSICAL PROPERTIES OF THE SHELL It is perfectly sensible to believe that QCD admits a metastable vacuum even at zero temperature. The entire physics of Bag model is within this type of framework. Inside of the bag is in the deconfining vacuum which is metastable with higher potential energy (the bag constant). The region outside the bag is in the confining vacuum with zero potential energy, with a wall separating the two regions. This is entirely consistent with the picture of a first order deconfinement-confinement phase transition at finite temperature where one usually writes down a finite temperature effective potential as a function of an appropriate order parameter. The expectation value of the Polyakov line < L > has been used to describe the physics of this phase transition at finite temperature [21,24], where < L > is given by, < L >= 1 N c < trpe i β 0 A 0dτ >, (2) where N c is the number of colors, β = 1/T, P denotes path ordering, and A 0 is the time 6

7 component of the vector potential. < L > vanishes in the confining phase, while it is non-zero in the deconfining phase, breaking the Z(3) symmetry spontaneously [21,24]. On the other hand, at T = 0, the physics of bag model can be captured by postulating a color dielectric field χ with the following form of the effective potential for χ [25,26], V(χ) = m2 gb 2 χ2 [1 2(1 2 α ) χ +(1 3 σ v α )χ2 ]. (3) with corresponding Lagrangian density being L = 1 2 ( µχ) 2 V(χ). Here, σ v = 2αB/m 2 gb, m gb is the glueball mass and B is the bag constant. We take α = 24, B 1/4 = MeV, and m gb = MeV. These are the parameter values used in ref. [25], and we use these as a sample. The relevance of these parameters for us is in determining the values of surface tension σ and false vacuum energy density ρ of the deconfining bubble. We will present results for a wide range of σ and ρ. The plot of V(χ) is shown in Fig.1. χ = 0 in the confining region outside the bag, while it has a non-zero value in the deconfining region inside the bag. χ thus has the same physical behavior as < L >, the expectation value of the Polyakov line L. It would seem rather superfluous to have two entirely different characterizations of same physical phenomenon, i.e. two different phases of the same system, separated by an interface. One would like to think of the color dielectric field χ as capturing the physics of < L >. With this identification, one has a single, unifying picture where the effective potential for a strongly interacting system is always of the form given in Eq.(3), with χ being interpreted as either the color dielectric field of the bag model, or the expectation value of the Polyakov line L (neglecting the symmetry aspects of the order parameter). At finite temperature, the coefficients in V(χ) will become temperature dependent which will change the values of surface tension of the interface, latent heat etc. compared to the values obtained from the parameter values given above. The important thing is that the barrier between the deconfining vacuum and the confining vacuum must survive even at T = 0 if the physics of the bag model has to be preserved. We mention that the forms of effective potential such as the one used in ref. [27] do not capture this physics of deconfinement-confinement phase transition. For example, in ref. [27], an effective potential has been used with a Tφ 3 term (with φ being the order parameter for D-C phase transition) to model the first order nature of the transition. Thus in this model, barrier between the deconfining vacuum and the confining vacuum disappears below certain temperature, which is not in the spirit of the physics of bag model. For the same reason we do not use the parameterization used in ref. [24] as that model also is appropriate for describing physics near and above T c. The final picture of our model can now be presented. We start with the initial time τ i, with dense partons filling the collision region. This parton region will correspond to the region between the two nuclei, receding after overlap, with transverse radius of the region being equal tor A, thenuclear radius. Much ofthe collisionenergy will becarriedby partons in the fragmentation region. We are interested in the total energy available to the expanding parton system in the central region, which will eventually lead to the formation of the false vacuum bubble. For this, we take the total energy E tot to be, σ 2 v E tot = ǫ i πr 2 A z i, (4) 7

8 where ǫ i is the initial energy density expected in the collision (as given in Eq.(1)). R A 1.1A 1/3 fm, and z i is the initial thickness of the central region. For central region near z = 0, z i corresponding to a rapidity interval Y at the initial time τ i can be taken as z i = τ i Y (see, ref. [16]). Estimates forτ i, initial QGPformationtime, varyfromτ 0.15 fm (for LHC) to τ 0.22 fm for RHIC. In view of various uncertainties [15,16,17] in the estimates of ǫ i, and τ i, (as well as in determining the relation between the plasma formation time τ i and the initial longitudinal extent of the plasma), we will present results for several values of z i. We will take z i = 0.15 fm for LHC, and equal to 0.22 fm for RHIC, and an optimistic value z i = 1 fm for both LHC and RHIC which will correspond to a large value of E tot. It is important to mention that our results only depend on the value of E tot. Thus, large value of E tot may arise from uncertainties in ǫ i, or z i. For Pb-Pb collision at LHC with s = 5.5 TeV, expected initial energy density in the central region is about 1.3 TeV/fm 3. This gives E tot 25, and 175 TeV corresponding to z i = 0.15, and 1 fm respectively. For Au-Au collision at RHIC with s = 200 GeV, the energy is much smaller and we find E tot 3, and 10 TeV for z i = 0.22, and 1 fm respectively. This parton system then expands and cools (for non-equilibrium expansion, we say that the energy density of parton system decreases as temperature may not be defined). Important thing is that at the initial time τ i, the parton system must be in the deconfining phase (either due to temperature > T c, or due to non-thermal transition with extremely high energy density of partons). The region outside the initial parton system being in the confining phase, there must be an interface separating the two regions. Fig.2a shows such a situation. Initial expansion will be longitudinal for proper times τ < R A, and will become threedimensional expansion for larger times. It is expected that energy density (or the temperature) still remains sufficiently high so that by the time expansion becomes three-dimensional, the parton system still remains in the deconfining phase. Fig. 2b denotes such a stage of the central parton system at τ R A. (We mention here that it is not crucial for our model whether the fragmentation region is inside this region shown in Fig.2b, or falls outside it.) For subsequent expansion, we assume (as explained in detail above) that a shell structure starts developing. For times much larger than R A, one expects an almost complete depletion of partons in the central region [9], with all partons piled up near the shell boundary. Fig. 2c denotes this intermediate stage. (We can take the proper time τ as measured in a frame at rest at some average point in the middle of the shell. Velocity there will not be very close to c, so τ will not be much different from the time t in the center of mass frame. In any case, the 3-dimensional expansion is not expected to be ultra-relativistic, so at this stage, one can use τ and lab time t interchangeably. We also neglect any Lorentz contraction in the thickness of the parton shell.) The stage shown in Fig.2c is the non-trivial part of our model. Here, the entire region inside the sphere of radius ct is taken to be in the deconfining vacuum. There is then the interface at the boundary. The partons accumulate near the boundary of the shell from inside. The thickness of this region containing the partons will be expected to be of order of about 2 R A as that is the total extent of the region at the time at which three dimensional expansion begins. The region with r < ct 2R A has no partons, but is still in the deconfining vacuum. This somewhat unconventional picture can be further justified as follows. Starting from a stage as shown in Fig.2b, where the entire interior region was in the deconfining vacuum, 8

9 one reaches the stage shown in Fig.2c by parton expansion. The inner region (devoid of partons) can turn into the confining vacuum only if there is nucleation of bubbles of confining vacuum in that region, which upon coalescing will convert the inner region into the confining vacuum. This could happen if the inner region cooled to T = 0 (with zero parton density, assuming zero baryon number in the central region) maintaining thermal equilibrium. However, thermal nucleation of bubbles could not happen since we are assuming that freezeout occurred before the hadronization transition (or, due to large time scales of bubble nucleation [28]). Even if the system did not go through finite temperature bubble nucleation when temperature decreased below T c, the inside region can still be in the confining vacuum if quantum nucleation of confining vacuum bubbles (at T = 0) could happen. As we will discuss below, the action S 0 of such a bubble is about 100 in natural units. The probability of nucleation of such a bubble is proportional to e S 0, and hence is completely negligible for the relevant time scales, allowing for reasonable dimensional estimates for the pre-factor for the nucleation probability. (As in our model a significant barrier separating the two vacua remains even at T = 0, there is no possibility of the phase transition occurring via spinodal decomposition.) Thus, we conclude that the inner region of the shell must correspond to T = 0 supercooled deconfining vacuum. Eventually, as parton shell keeps expanding, its energy density falls below a critical value so that hadronization must take place. Partons will then form hadrons as they leave the spherical region bounded by the interface. This is the stage when the interface will stop moving outward, and will start shrinking. As the interface passes through the parton system, it converts partons into hadrons, as shown in Fig.2d. (Alternatively, hadronization may happen at various places in the parton shell. This will also amount to an effective shrinking of the interface.) The motion of the interface will be dissipative at this stage due to its interaction with partons. Typical velocity of the interface, therefore, will be less than the speed of sound [29,27,28]. However, as all (or most) of the partons hadronize, the interface will continue to move inward due to negative pressure of the metastable vacuum. Now there are no partons to impede the motion of the interface. What one is left with is a spherical bubble of pure false vacuum. To obtain the profile of the wall of this false vacuum bubble, we first obtain the profile of the true vacuum bubble as follows. Given V(χ) in Eq.(3), one can obtain the instanton solution for the tunneling through the barrier from χ 0 metastable vacuum to the χ = 0 true vacuum. This is given by the solution of the following equation [30], d 2 χ dr + κ dχ 2 r dr V (χ) = 0, (5) where V(χ) is the effective potential in Eq.(3) and r is the radial coordinate in the Euclidean space. In the Minkowski space initial profile for this true vacuum bubble is obtained by putting t = 0 in the solution of theabove equation. Parameter κ = 3 forquantum nucleation of bubble (at T = 0), while κ = 2 for thermal nucleation of bubble at finite temperature. Initially the parton region had large temperature, so the profile of the interface at that stage would be given by the finite temperature bubble (i.e. κ = 2). This will be true up to the stage shown in Fig.2b. After freezeout, or once the partons leave the shell in the form of hadrons as in Fig.2e, the profile of the interface will be given by T = 0 bubble, that is κ = 3 case. We have solved Eq.(5) using a fourth order Runge Kutta algorithm for the 9

10 effective potential as given in Eq.(3). The numerical technique is the same as used earlier for standard false vacuum decay [31] with the obvious difference that now the false vacuum occurs at non-zero value of χ while the true vacuum occurs at χ = 0. Thus the search for the bounce solution has to be done differently in the present case. In Fig.3a we have shown the solution of the T = 0 quantum bubble of true vacuum. The radius of this bubble is about 5 fm. (The finite temperature bubble has a smaller radius 3.8 fm.) We have calculated the action for the true vacuum bubble shown in Fig.3a and we find it to be about 100 in natural units. As mentioned above, this large value of the action suppresses the nucleation of true vacuum bubbles in the the central deconfined (T = 0) region, so that this region converts to the confining vacuum only via collapse of the spherical interface shown in Fig.2e. Fig.3a shows the bubble of confining vacuum embedded in the deconfining vacuum (as obtained by the instanton solution of Eq.(5)). However, the experimental situation we have discussed above is exactly the reverse. As shown in Fig.2e, we have a large bubble of deconfining vacuum which is embedded in the confining vacuum. Such a bubble cannot be obtained as a classical solution of Eq.(5). The situation shown in Fig.2 arises because of changing temperature, and is similar to the false vacuum bubble formation due to coalescence of true vacuum bubbles as discussed in the context of quark-hadron transition in the early universe [32]. Main difference between our case and the case discussed in ref. [32] (see, also ref. [28]) is that there the motion of bubble wall remains dissipative due to presence of QGP in the interior of the bubble. In contrast, in our case, the special geometry of collision gives the shell like structure which leads to a pure false vacuum bubble with no QGP inside. The motion of the interface, therefore, will not be dissipative. Fig.3b shows this false vacuum bubble which is simply obtained by inverting the profile of the bubble as given in Fig.3a. This is a somewhat approximate way of generating the appropriate bubble profile. However, our interest is only in the values of surface tension σ of the interface and the false vacuum energy density ρ. (Note that ρ is the same as the bag constant B. We use a different notation for it due to the method by which we calculate it by generating the confining vacuum bubble and then inverting it.) We will present results for a range of values of these parameters, with the profile given in Fig.3b corresponding to one set of these values. Note that the radius of this deconfining vacuum bubble (Fig.3b) is not obtained from equation of motion. In contrast, radius of the confining vacuum bubble of Fig.3a was fixed by the solution of Eq.(5). Inverted profile in Fig.3b is generated only to estimate the values of σ and ρ corresponding to parameter values in Eq.(3). Its radius will be determined by the physics of the problem, such as total energy of collision, as we discuss below. In the same way, the radius of the false vacuum bubble in ref. [32] is determined by separation between the nucleation sites of the hadronic bubbles, and the manner in which they coalesce. To obtain values of surface tension σ and the false vacuum energy density ρ of the deconfining vacuum bubble, we obtain total energy E of the false vacuum bubble by numerically 1 integrating the energy density, 2 ( χ)2 + V(χ), for the bubble profile in Fig.3b. Bubble energies are obtained in this manner for two different values of bubble radius. σ and ρ are then determined by using the following equation for the two values of bubble energy and radius, E = 4πr 2 σ + 4π 3 r3 ρ. (6) 10

11 The values of σ and ρ obtained by using this equation for two values of bubble radius are tested for other bubble radii and it is found that bubble energies are accurately reproduced. For the parameters used for Eq.(3) in ref. [25], we find σ = 64.8MeV/fm 3 and ρ = 27.8MeV/fm 3 ( (122MeV) 4 ). Note that this value of ρ is the same as the bag constant B in Eq.(3) as it should be. This gives us confidence that this procedure of inverting profile of confining vacuum bubble to get the profile of the deconfining vacuum bubble is reasonably accurate. For comparison, we repeated this procedure for the finite temperature bubble, i.e. solution of Eq.(5) with κ = 2, (pretending that the effective potential in Eq.(3) corresponds to the finite temperature case). As we mentioned, resulting bubble radius of the confining vacuum bubble (similar to the T = 0 bubble in Fig.3a) is about 3.8 fm. We find that the inverted bubble (deconfining vacuum bubble, as in Fig.3b) has σ = 64.4 MeV/fm 2 and ρ = 27.8 MeV/fm 3. These values are essentially the same as in the case of T = 0 bubble. Since the only thing we need from the above calculations is the values of these parameters, it is immaterial whether we think of the interface as corresponding to the finite temperature bubble or to T = 0 bubble. IV. EXPANSION AND SUBSEQUENT ULTRA-RELATIVISTIC COLLAPSE OF THE SHELL We need to make one further specification about the expansion of the parton system. If the expansion is adiabatic, then total entropy will be conserved. In that situation, total energy will decrease due to the work done by the pressure of the plasma. This description is appropriate when expansion happens maintaining equilibrium, as in the hydrodynamical models of parton expansion. In the situation of free expansion, or if there are dissipative effects present, then expansion could be energy conserving [33]. We have argued above that shell picture is more reasonable in the case when parton system freezes out early. Thus energy conserving expansion may be more appropriate for our model. We consider this case first. We will also briefly discuss results for the adiabatic expansion, and we will see that the temperatures of the hot spot are much smaller for that case. (Note that the only thing relevant to us about the parton expansion is the development of a shell structure with central region devoid of partons, but still in the deconfined vacuum. We have argued that this seems natural when partons freezeout early. However, all of our arguments are valid if this picture can be justified even by taking partons to be in equilibrium, through the stages in Fig.2a-2d.) Assuming the energy to be conserved, the energy density of the partons in the shell, ǫ(r), at any stage shown in Fig.2c, is determined by, 4πr 2 σ + 4π 3 (r3 (r r) 3 )ǫ(r)+ 4π 3 r3 ρ = E tot, (7) where E tot is given in Eq.(4). Here ǫ(r) denotes the parton energy density in the shell of thickness r 2R A. The shell will expand to a largest size r max at which stage parton system will hadronize. We determine r max by taking ǫ(r) = ǫ c, the critical value of energy density ofpartonsbelowwhich hadronizationisexpected totakeplace. Wetakeǫ c π T4 c, where T c is the critical temperature for the quark-hadron transition. ǫ c is not to be taken necessarily as the energy density of plasma in equilibrium (with temperature = T c ), as we 11

12 have argued that the partons may be in a non-equilibrium state. We take ǫ c as giving the scale below which partons should convert to hadrons. The value of T c is related to the bag constant (for an equilibrium transition) using Gibbs criterion of equal pressure at transition temperature, T c = (B/(g q g h )) 1/4. Here, g q = 37π 2 /90 and g h = 3π 2 /90, corresponding to 2 massless quark flavors and 8 gluons in the QGP phase, and 3 pions in the hadronic phase. Again, the value of B appropriate at the scale given by T c may be different from the value of B at T = 0. For simplicity, we ignore any scale dependence of B and use the same value of B which equals the false vacuum energy density ρ, as also determining T c via above relation. Therefore, we present results for a large range of values of B, not limiting to only those values which correspond to realistic values of T c. r max gives the largest radius of the shell, with interface being at the outer boundary (i.e. at r = r max ) initially. As this is the stage when hadrons start forming, interface starts moving inwards. Initial motion of the interface is highly dissipative [29,27], as long as it traverses the region which is filled with the partons, i.e. a thickness of r 2R A. This stage is shown in Fig.2d. Once the interface has shrunk below this parton filled region, it is free to undergo unimpeded, relativistic collapse. The region bounded by this interface at a later stage (shown in Fig.2e) represents a pure false vacuum bubble with no partons inside. The initial radius of this pure false vacuum bubble, r f will be approximately given by, r f r max R A, (8) where r max is determined by solving Eq.(7) with ǫ(r) = ǫ c. In writing this expression for r f, we have taken that when ǫ(r) = ǫ c is achieved during the shell expansion, the interface starts collapsing with almost speed of light, as the parton shell keeps expanding relativistically. However, due to dissipative motion of the wall through the parton system, velocity of the wall will be typically much smaller than the speed of light [29]. If one takes the wall to be almost static, as compared to the outward velocity of the shell, then one will expect r f to be almost equal to r max. Thus, the value of r f as given above is an underestimate. Resulting total energy of the false vacuum bubble, and hence the final temperature of the hot spot will also be somewhat underestimated in using the above equation. This should hopefully compensate, to some degree, the effects of assuming all partons to be ultra-relativistic, (e.g., if the parton distribution extends beyond the assumed shell thickness of 2r max, it will reduce the value of r f ). Total energy E f of this pure false vacuum bubble is, E f = 4πr 2 fσ + 4π 3 r3 fρ. (9) Further evolution of this false vacuum bubble is well understood. It will undergo relativistic collapse. Due to surface tension, bubble will become more spherical as it collapses. This is important asthere arevariousfactorsbecause ofwhich theinitial shape of thebubble wall may not be spherical. First, the collision geometry itself will not be expected to give rise to completely spherical structure. Secondly, the interface motion through the parton shell itself may not be very isotropic. However, during the free collapse of the interface, the shape should become more spherical due to surface tension. Any remaining asphericity will ultimately affect the final radius to which the bubble can collapse before the bubble walls 12

13 decay via collision. For simplicity, we will assume that the interface collapses maintaining spherical symmetry. During bubble collapse, the potential energy of the false vacuum simply gets converted into the kinetic energy of the collapsing interface [30]. This is where our model crucially differs from the collapsing QGP bubbles previously discussed in the context of quark-hadron transition [32,28]. There, interface always moves through a QGP system filling the interior of the bubble, which impedes the motion of the interface. The false vacuum energy there gets converted to the heat which raises the temperature of the plasma [32,28]. In contrast, there are no partons in the interior of the deconfining vacuum bubble in our model. The entire false vacuum energy thus converts to the kinetic energy of the interface. The collapsing interface quickly becomes ultra-relativistic, with extremely large Lorentz contraction factor, as indicated by thinner interface in Fig.2e. At any stage during the collapse of the shell, the value of the Lorentz γ factor primarily depends on the initial radius r f of the shell, and the value of ρ (i.e. B). To give an idea of how rapidly the bubble wall Lorentz contracts, we give values of γ at the stage when the bubble has collapsed to a radius of 1 fm. For Pb-Pb collision with s = 5.5 TeV, and with z i = 1 fm, we find that r f varies from 20 to 90 fm as B 1/4 is reduced from 240 MeV to 120 MeV (as we will see below). Value of γ, when bubble has collapsed to a radius of 1 fm, ranges from to For same parameters, but with z i = 0.15 fm, γ ranges from 2000 to At RHIC energy, with s = 200 GeV for Au-Au collision (with z i = 1 fm), r f ranges from 5 fm to 35 fm as B 1/4 is reduced from 250 MeV to 100 MeV. Resulting γ factor (again, when bubble has collapsed to 1 fm radius) ranges from about 500 to For this case, by the stage when bubbles radius has decreased by another order of magnitude, i.e. to a value of 0.1 fm, γ ranges from about to The thickness of the interface initially (i.e. at stages shown in Fig.2a-2d), is about 1 fm. We see that by the time bubble collapses to a radius of about 1 fm, the thickness of the interface is much smaller, about 10 3 fm to 10 5 fm. Interface thickness reduces to about 10 5 fm to 10 7 fm by the time the bubble collapses to a radius of 0.1 fm. Again, it is unconventional to talk about such large length contraction factors in the context of heavyion collisions. We know that due to virtual partons at small x (the so called wee partons), nucleus thickness does not Lorentz contract below about 1 fm even at ultra-high energies [34] (see, also, ref. [35]). However, this limiting Lorentz contraction arises due to virtual particle production in the color field of the partons. In our case, by the time the interface contracts below the partonic shell, there will be no such virtual partons inside the bubble. Thus, there is no compelling reason to believe that the Lorentz contraction of the interface of this pure false vacuum bubble should be limited by the typical QCD scale of 1 fm. Basically, the radial profile of the interface represents the kink solution (Eq.(5)) which can be Lorentz boosted to any velocity in the absence of interactions with partons. Even with quantum fluctuations on the background of this ultra-relativistic bubble wall, there is no reason to expect that these fluctuations will put a stop to the highly Lorentz boosted bubble wall. Though some of the wall energy may get dissipated due to these effects, bubble collapse may still continue to sizes much smaller than 1 fm as long as walls remain ultra-relativistic. This extreme Lorentz contraction of the bubble wall has following important consequence. Normally, the process of the collapse of a bubble wall will stop when its radius is of the order of the thickness of the wall. (See, for example, the numerical study of collapsing domain walls in ref.[36]. Though, for ultra-relativistic collapse numerical errors can build up and one 13

14 needs more sophisticated numerical techniques, see ref. [37].) Without significant Lorentz contraction, we would expect the bubble collapse to halt at the stage when bubble radius is of the order of 1 fm, as in the standard treatments of collapse of QGP bubbles in a quarkhadron transition [32]. Then bubble walls (of different bubble portions) will collide and the energy of walls will be released in the form of particle production. However, in our case the bubble wall undergoes large Lorentz contraction even when bubble radius is about 1 fm (in the center of mass frame). As compared to the wall thickness, bubble radius is about 10 3 to 10 5 times larger at this stage. Thus, there is no reason to expect that bubble collapse will halt at this stage. In fact it is easy to see (using the fact that as bubble collapses, false vacuum energy gets converted to the kinetic energy of the bubble wall) that bubble wall thickness decreases much faster (due to Lorentz contraction) compared to the bubble radius, as bubble collapse proceeds. Of course, eventually this type of classical evolution of bubble must stop, at least by the stage when the net size of the system has become smaller than that allowed by the uncertainty relation. The net energy E f we start with (of the false vacuum bubble at the initial stage when free collapse of the interface commences), is much larger than a TeV. Thus, when collapsing bubble size becomes smaller than Ef 1 (say a TeV 1 ), it seems reasonable that bubble collapse may halt. Entire energy E f of the initial pure false vacuum bubble, which was subsequently converted into the kinetic energy of ultra-relativistic walls, may thus get converted to particles in a region of size less than (E f ) 1. (It will be interesting to work out the exact conditions when bubble collapse will halt. Even if the bubble radius shrinks only to a value of say 0.1 fm, one still gets a hot spot, though with a lower temperature. It is also possible that in some situations, bubble collapse may develop strong anisotropies. In that case the bubble may break into smaller bubbles.) It is important to realize that the above picture only requires that most, or a significant fraction, of the bubble wall energy (after bubble has collapsed to size of order 1 fm) is not dissipated away subsequently until different portions of the (Lorentz contracted) bubble wall collide with each other. This is because essentially all of the energy of the initial false vacuum bubble is already stored in the bubble walls by the time its radius shrinks to a value of order 1 fm (for initial radius much larger than 1 fm). Even with quantum fluctuations on the bubble wall background, there is no reason to expect that these fluctuations will tend to put an early stop on (or significantly dissipate the kinetic energy of) the highly relativistic bubble wall. For example, even for the case of nucleus-nucleus collision, interactions of wee partons do not stop the nuclei from overlapping (or going through) at sufficiently large collision energies. Though, we again emphasize that a nucleus, filled with partons, is a very different object than a pure false vacuum bubble. Thus, the arguments based on wee partons made for nucleus [34] do not extend in a natural way to the Lorentz contracted bubble wall. After all, one does not even know what is the correct field describing such a wall (apart from various possible order parameters), let alone the detailed nature of quantum fluctuations above the bubble wall separating deconfined and confined vacua of QCD (at zero temperature) and their Lorentz transformation properties. From all this, we would like to conclude that there seems a genuine possibility that bubble collapse may continue, preserving most of its energy, down to sizes much smaller than 1 fm. As we will see below, this possibility leads to very interesting implications. The final conclusion is that the entire energy E f of the false vacuum bubble will be 14