**1. Introduction**

A complete understanding of multi-particle production dynamics continues to be a challenge for theory in high energy collisions. The full description of real-time dynamical evolution in a strongly coupled non-Abelian gauge theory can be notoriously difficult. The availability of large, diverse, and high quality accumulated proton-proton (*pp*) and heavy ion (HI) collision data from the energy frontier at the Large Hadron Collider (LHC) and at Brookhaven National Laboratory's Relativistic Heavy Ion Collider (RHIC) are providing new insights into puzzling behavior observed in these strong interaction processes. Additionally, studies of similar intensity frontier anomalies in (anti)neutrino-nucleus scattering gathered recently also benefit from newly accumulated event statistics at Fermilab. The transverse momentum distribution of differential cross sections in *pp* collisions shows process-dependent behavior that requires a more subtle explanation compared to previous ideas.

For example, the differential distribution of charged hadrons resulting from *pp* collisions are observed to exhibit both a well-understood quark and gluon "hard scattering" component at high transverse momenta that can best be described by a power-law fit to the distribution, and a component at low transverse momenta that exhibits a less well-understood "thermal" behavior that is best represented by an exponential fit to the data [1, 2]. See [3, 4] for a review. The sum of these two contributions is necessary to properly characterize the transverse momentum differential distribution. And yet, for diffractive events where there is a large rapidity gap in the transverse momentum distribution of the differential cross section for photo-produced muon pairs, there is only a hard-scattering contribution; the thermal component disappears in this class of processes even though these inelastic events produce large number of hadrons. Charged-current weak interactions also exhibit differential cross section momentum distributions that have process-dependent results. As in the case of strong interactions, for anti-neutrino scattering from a nucleon inside a nucleus, there is also a "hard scattering" component at high particle momenta that is manifested by a power-law fit to the momentum distribution of the differential cross section. And at low momenta, the differential cross section behavior is best described by an exponential fit to the momentum distribution in this region. Furthermore, in coherent anti-neutrino scattering from a nucleus, where the nucleus remains intact, there is no thermal component to the differential cross section momentum distribution, only the hardscattering component [5].

It is interesting to consider the possibility that these two interactions of vastly differing collision or scattering energies are different manifestations of a single underlying fundamental process. There is growing interest in the link between quantum entanglement, entanglement entropy (EE), and high energy physics presently. Here, we describe a relationship between quantum entanglement in the nucleon wave functions associated with the hadron collisions at the LHC experiments and with the electroweak scattering in Fermilab experiments. There are several examples of these relationships in theoretical physics: Research on the dynamics of quantum entanglement and entanglement entropy in the regime of small Bjorken-*x* in deep inelastic scattering, in electromagnetic interactions [6]. Additionally in [1] the case is made that quantum entanglement between partons inside a nucleon can be probed by deep inelastic lepton scattering. Deeper insights into these dynamics involving entanglement entropy in the regions of black holes is now provided by the AdS/CFT correspondence [7]. Well-understood quarkantiquark correlations are now improved even more due to considerations of entanglement entropy in Lattice Gauge Theory [8] that are or soon will be augmented by new and planned parton distribution functions data in particle and nuclear physics. Newer intuition into this relationship that involves quantum entanglement, entanglement entropy, thermal behavior, is gained from heavy ion and proton-proton collisions [9, 10], nuclear shadowing effects [11], and chiral symmetry breaking [12].

In [1, 2, 13], it is proposed that the thermal component is a result of entanglement between causally disconnected parts of the nucleon in the interaction. For this reason inelastic *pp* collisions exhibit a thermal component, while diffractive collisions, where the nucleon as a whole is probed, give rise to only the hard-scattering component, and no thermal behavior. That the thermal component is a consequence of quantum entanglement between different regions of the colliding protons wave functions in *pp* colisions is proposed in [14]. A thermal behavior can, in the collisions where there are large number of particle produced in the final state, be the result of rescattering among the produced final state particles. However, the thermal component was shown to not only be present in the transverse momentum distribution of charged hadrons, but also in the transverse momentum spectrum of Higgs bosons production and decay differential cross section, resulting from the collisions. For Higgs boson production in *pp* collisions, there are very few final state particles for rescattering compared to heavy ion collisions for example.

*Quantum Information Science in High Energy Physics DOI: http://dx.doi.org/10.5772/intechopen.98577*

If quantum entanglement and EE are responsible for the thermal behavior in charged hadron as well as Higgs boson production in *pp* collisions at the LHC, then it should also be observed in neutrino-nucleus scattering, where only a fraction of the struck nucleon in the interaction is probed by the exchanged charged-current probe. As explained in [5] the observation of a thermal component exists in the momentum distribution of neutral pions from antineutrinonucleon scattering, while no such component is present in coherent antineutrinonucleon scattering. This absence of a thermal component in the latter case is due to the fact that antineutrinos probe the nucleus as a whole in coherent scattering; there is no un-probed region of the nucleus which can be entangled with the probed region. So with no quantum entanglement in the interaction, the thermal component in the momentum distribution of the differential cross section is absent, as expected.

These topics are presented and discussed in this chapter. A brief description of the theory motivating this proposed link between quantum entanglement, entanglement entropy, and thermal behavior in *pp* collisions is given in Section 2 and subSection 2.1. This is followed by experimental results of transverse momentum distributions in *pp* collisions at 13 TeV collisions energy. Charged hadron production where (and why) both hard scattering and thermal components are present in the differential distribution are described in subSection 2.2. The absence of the thermal component in diffractive production of muon pairs in the reaction *pp* ! *μ*þ*μ*�*X* is explained in subSection 2.3. The interesting need for superposition of both the hard scattering and thermal components to describe the transverse momentum distribution of the Higgs boson is presented in subSection 2.4. Section 3 includes the presentation of this phenomena in charged current weak interactions. Concluding remarks are given in Section 4.

## **2. Entanglement entropy and thermal behavior in the strong interaction**

We begin by considering the possibility that the observed thermalization in *pp* collisions is the result of a sudden perturbation or rapid "quench" due to the high degree of entanglement inside the protons involved in the collision [15]. The link between quantum thermalization and quantum entanglement is shown to exist in an experimental quench in Bose-Einstein condensates of Rb atoms in atomic and condensed matter physics [16]; the rapid eigenstate thermalization was found to be the result of a quantum entanglement. In *pp* process described here, low momenta correspond to late times after the collision. The thermal behavior begins to dominate over the hard scattering component in the transverse momentum distribution at late times. This is consistent with theoretical studies in 1ð Þ þ 1 dimensional conformal field theories of quenches in entangled quantum systems [17–19] where a system can be described by a generalized thermal Gibbs ensemble at late times.

Since a high-energy collision can be viewed as a rapid quench of the entangled partonic state [15], it is thus possible that the effective temperature inferred from the transverse momentum distributions of the secondaries in a collision can depend upon the momentum transfer, that is an ultraviolet cutoff on the quantum modes resolved by the collision. In analyzing the high-energy collisions with different characteristic momentum transfer *Q* we thus expect to find different effective temperatures *T* � *Q*. We can also look at the inelastic events characterized by a rapidity gap, where the proton is probed as a whole, and no entanglement entropy arises [15]. In this case, if the quantum entanglement is responsible for the thermalization, we expect no thermal radiation.

The presence of both a thermal and a hard scattering component in inclusive deep-inelastic scattering at HERA has been observed [20]. And the absence of this thermal component in processes characterized by a rapidity gap is also manifested in these studes. In diffractive events where there is a rapidity gap, the entire proton wave function is involved in the scattering process. In diffractive scattering, the proton remains fully intact in the central part of the collider detector where scattering takes place. There is no entanglement entropy due to different regions of the proton wave function being involved in the scattering in different ways. This observation points to a connection between this thermalization and quantum entanglement between different parts of the proton wave function. This link is described in the next section.

#### **2.1 Entanglement and thermalization in high energy collisions - theory**

A brief summary of these proposition is as follows. The hard process in *pp* collisions probes only the part of the proton wave function that is localized in a region of space denoted here as **A**. For a hard process such as the one shown in **Figure 1** this region has a transverse size that is less that the full proton diameter and, in the proton's rest frame, longitudinal size in terms of Bjorken-x is � ð Þ *mx* �<sup>1</sup> , where *m* is the proton mass.

In this same figure, the spatial region **B** is complementary to **A**, that is, the entire space is **A**∪**B**. Hard processes have their origin in the physical states inside the region **A**. These are states in a Hilbert space H*<sup>A</sup>* of dimension *nA*. Unobserved states (not part of hard scattering) in the region **B** belong to the Hilbert space H*<sup>B</sup>* of dimension *nB*. With this picture, the protons prior to the collision, both composite systems in **A**∪**B** (the entire proton in each case) are then separately described by the vector represented as, for example, ∣*ψAB*i in a tensor product of the two spaces H*<sup>A</sup>* ⊗ H*B*:

$$|\psi\_{AB}\rangle = \sum\_{i,j} c\_{ij} \ |\rho\_i^A\rangle \otimes |\rho\_j^B\rangle,\tag{1}$$

where *cij* are the elements of the matrix *C* that has a dimension *nA* � *nB*. In the case where there are states <sup>∣</sup>*φ<sup>A</sup>*<sup>i</sup> and <sup>∣</sup>*φ<sup>B</sup>*<sup>i</sup> that <sup>∣</sup>*ψAB*i ¼ <sup>∣</sup>*φ<sup>A</sup>*<sup>i</sup> ⊗ ∣*φ<sup>B</sup>*i, where that the sum (Eq. (1)) contains only one term, then the state ∣*ψAB*i is product state that is separable. In the case where it is not separable, ∣*ψAB*i is entangled.

∣*ψAB*i is called a bi-partite system that, making use of the Schmidt decomposition theorem, can be expanded as a *single* sum in *n* instead of a double sum over *ij*.

#### **Figure 1.**

*Characterization of the entanglement entropy in pp collisions. In the leftmost depiction, the collider protons, before collision, are both pure states. In the rightmost depiction, during the pp scattering, there exists the proton overlap collision region (A) and the overlap spectator region (B).*

*Quantum Information Science in High Energy Physics DOI: http://dx.doi.org/10.5772/intechopen.98577*

$$|\psi\_{AB}\rangle = \sum\_{n} a\_n |\psi\_n^A\rangle |\psi\_n^B\rangle \tag{2}$$

Here ∣*ψ<sup>A</sup> <sup>n</sup>* <sup>i</sup> and <sup>∣</sup>*ψ<sup>B</sup> <sup>n</sup>*i are orthonormal sets of states (properly chosen) localized in the domains **A** and **B**, respectively. And *α<sup>n</sup>* are real, positive numbers that are the square roots of the eigenvalues of matrix *CC*† .

The density matrix formalism is now a better tool to use in the discussion. For a mixed state that is probed in region *A*, the density matrix can be expressed as

$$\rho\_A = \text{tr}\_B \,\, \rho\_{AB} = \sum\_n \alpha\_n^2 |\psi\_n^A\rangle\langle\psi\_n^A|,\tag{3}$$

where the symbol *α*<sup>2</sup> *<sup>n</sup>* � *pn* denotes the probability of an *n*-parton state. The basis ∣*ψ<sup>A</sup> <sup>n</sup>* i in (Eq. (2)) with states having a fixed number of *n* partons does not have interference between states with different number of partons due to the fact that this sort of interference is absent in the parton model. In this Schmidt decomposition, there can be an infinite number of terms (the Schmidt rank) in the sum shown in (Eq. (2)). A Schmidt rank one state is then a pure product state that does not include entanglement.

In the case of a mixed state, the probabilities corresponding to the different states described above can be used to define the von Neumann entropy of the mixed state given by

$$S = -\sum\_{n} p\_n \ln p\_n. \tag{4}$$

It is the entanglement between regions **A** and **B** defined above that gives rise to what is called entanglement entropy here, and is related to Shannon entropy in information theory shown in (Eq. (4)). Hence the entanglement entropy can be determined form the QCD evolution equations that are used to evaluate the probabilities *pn*. After the hard scattering takes place in the collision, the mixed quantum state characterized by the entanglement entropy (Eq. (4)) undergoes the evolution towards the final asymptotic state of hadrons that are measured by the detectors. This final state is characterized by the Boltzmann entropy. Further discussions of the relationship between Boltzmann entropy and entanglement entropy can be found in [14].

Studies of quantum entanglement and thermalization in atomic and condensed matter physics were shown to depend upon the quench properties, and that there is evidence for quantum propogation and information propagation [16, 21–23]. It is instructive to compare this with a quench induced by a high energy collision. The quench associated with the latter [17, 18] leads to the following interpretation. A quench produces a highly excited state of a Hamiltonian *H* ¼ *H*<sup>0</sup> þ *V t*ð Þ from what was the ground state of an unperturbed Hamiltonian *H*<sup>0</sup> originally. Here *V t*ð Þ is the term induced by the inelastic collision. Gluon exchange in the strong interaction induces the inelastic interaction, so the term *V t*ð Þ is seen to represent an effect of the pulse of the color field. The onset of this pulse in a hard scattering with a hardness scale *Q*, by the uncertainty principle, is *τ* � 1*=Q* where *τ* is the proper time. Since this time is short on the QCD scale, *τ* ≪ 1*=*Λ, the quench creates a highly excited multi-particle state. A short pulse of (chromo)electric field produces particles that have a thermal-like exponential spectra. The thermal spectrum in this case can be attributed to the emergence of an event horizon formed due to the acceleration induced by the electric field. Associated with this system is an effective temperature of *<sup>T</sup>* <sup>≃</sup>ð Þ <sup>2</sup>*πτ* �<sup>1</sup> <sup>≃</sup> *<sup>Q</sup>=*ð Þ <sup>2</sup>*<sup>π</sup>* [24–27].

As shown in [17, 18] for a rapid quench (such as the one that occurs in a highenergy collision) in a 1ð Þ þ 1 dimensional CFT the entanglement entropy of a segment of length *L* first grows linearly in time, until *t* ≃*L=*2, and then saturates at the value

$$S(t) \simeq \frac{c}{3} \ln \tau\_0 + \frac{\pi c L}{12 \tau\_0}. \tag{5}$$

where *c* is the conformal charge of the CFT, and *τ*�<sup>1</sup> <sup>0</sup> is the energy cutoff for the ultraviolet modes [14]. A sketch of the picture of the resulting thermalization from entanglement caused by the quench is shown in **Figure 2**.

The interpretation of the result (5) is the following [17, 18]. The quench leads to the production of entangled (quasi)particle pairs, since what used to be the ground state of the undisturbed Hamiltonian *H*<sup>0</sup> is a highly excited state of the Hamiltonian after the quench, *H* ¼ *H*<sup>0</sup> þ *V t*ð Þ. The entangled pairs produced by the quench propagate along the light cone, and contribute to the entanglement entropy of the segment of length *L* if only one particle of the pair is detected within this segment. Shortly after the quench, only particle pairs produced near the boundary of the segment thus contribute to the entanglement, and the entanglement entropy is not extensive in the length *L*. However, at times *t*> *L=*2, even in the center of the segment one can detect a particle whose entangled partner is outside of the segment – this means that the entanglement entropy receives contributions from the entire segment, and should scale extensively in *L* in accord with the result (5). This scaling is a necessary condition for an effective thermalization.

For a quench induced by a high-energy collision, we sketch the resulting picture of thermalization from entanglement in **Figure 2**. Note that the hardest quasiparticle modes that propagate along the light cone thermalize first. For the softer

#### **Figure 2.**

*An illustration of the onset of quantum thermalization through entanglement in a high energy pp collision. Time runs along the vertical axis, while space runs along the horizontal axis. The outermost lines define the light cone. The variables used are defined in the text. Entangled particle pairs that are produced at a proper time τ* <*τ*<sup>0</sup> *contribute to the entanglement entropy in the interval of length L shown by the hashed segment of the curve. Figure from [14].*

*Quantum Information Science in High Energy Physics DOI: http://dx.doi.org/10.5772/intechopen.98577*

particles that propagate in the interior of the light cone, it takes a longer time to thermalize, that is, to exhibit an extensive scaling of the entropy. The detection of particles is assumed to be performed within the interval of length *L* (see Eq. (5)), corresponding to a limited range in (pseudo)rapidity. While (Eq. (5)) has been obtained in the framework of CFT, the simple physical interpretation of this result makes its broader validity quite likely.

It is instructive to point out the difference in the mechanisms of thermalization expected at weak and strong coupling. At weak coupling, the "bottom-up" thermalization mechanism [28] also yields an effective temperature *T* � *Qs* in inelastic high energy collisions. However the thermalization in this picture begins from the soft, low-momentum modes that eventually draw the energy from the harder modes; the thermalization of the hard, high-momentum modes is thus expected to take a parametrically long time proportional to the inverse power of the (small) coupling constant [28]. On the other hand, in strongly coupled entangled systems the process of thermalization is fast and determined by the size of the system and the parameters of the quench; moreover, it starts from the *hardest* modes resolved in the process. In the dual holographic description of conformal field theory, this process is described by the formation of trapped surface near the Minkowski boundary that then falls into the AdS bulk, corresponding to the spreading of thermalization from hard to soft modes [29, 30]. A similar picture emerges from the analysis of entanglement entropy in an expanding string [31], where the entropy has been found to have a thermal form with an effective temperature *T* � 1*=τ* at early time *τ*.

#### **2.2 Charged hadron transverse momentum distribution**

The discussion presented in the previous sections provide motivation to compare with experimental results from inelastic collision events at high energies. It also gives the opportunity to explore the possible relation between effective temperature and the hard scale of the collision. Consider proton-proton collisions data recorded by the LHC ATLAS collaboration at ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> 13 TeV center of mass energy yield multiple charged particles in the final state [32]. The data presented here corresponds to 151 *μ*b�<sup>1</sup> of integrated luminosity for charged particles with greater than 100 MeV/c transverse momenta and absolute pseudorapidity of less than 2.5. Events with two or more final state charged particles were selected in the analysis. Final state hadrons that originate in the primary *pp* interaction and that have a lifetime of greater than 30 ps were excluded from the final selected events in order to remove the presence of charged particles that have strangeness or are from heavier flavors.

The normalized charged hadron transverse momentum distribution is shown in **Figure 3**. The thermal component is shown by the exponential, red dashed curve; we parameterize it as

$$\frac{1}{p\_T} \frac{d^2 N\_{ev}}{dp\_T} = \mathbf{A}\_{therm} \exp\left(-\mathbf{m}\_T/\mathbf{T}\_{\rm th}\right),\tag{6}$$

where *mT*, the hadron transverse mass, is defined as by *mT* � ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *<sup>m</sup>*<sup>2</sup> <sup>þ</sup> *<sup>p</sup>*<sup>2</sup> *T* p (*m* is the hadron mass; dominated by pions it is assumed), and *Tth* is an effective temperature. The hard scattering (power law, green solid curve) component is parameterized similar to [13],

$$\frac{1}{p\_T} \frac{d^2 N\_{ev}}{dp\_T} = \frac{\mathbf{A}\_{hard}}{\left(\mathbf{1} + \frac{\mathbf{m}\_T^2}{\mathbf{T}^2 \cdot \mathbf{n}}\right)^{\mathbf{n}}},\tag{7}$$

where *T* and *n* are parameters to be determined from the fit. The sum of the thermal and hard scattering contribution terms is shown by the blue solid curve in **Figure 3**.

The extracted value of the thermal temperature, T*th* ¼ 0*:*17 GeV describes well the experimental transverse momentum distribution, and it agrees with the temperature expected from the extrapolation of the relation [13] deduced at lower energies;

$$\mathbf{T}\_{th} = \mathbf{0}.\mathbf{0}\mathbf{9} \mathbf{8} \cdot \left(\sqrt{\mathbf{s}/\mathbf{s}\_0}\right)^{0.06} \mathbf{G} \mathbf{e} \mathbf{V} \tag{8}$$

to the LHC 13 TeV collision energy; here *<sup>s</sup>*<sup>0</sup> <sup>¼</sup> 1 GeV<sup>2</sup> . Similarly, the hard scale temperature parameter T is [13]

$$\mathbf{T} = \mathbf{0}.409 \cdot \left(\sqrt{\mathbf{s}/\mathbf{s}\_0}\right)^{0.06} \text{ GeV.} \tag{9}$$

It's interesting that the parameterizations (Eqs. (8) and (9)) imply that the effective thermal temperature T*th* is proportional to the hard scale temperature parameter T, which is in agreement with the Section 2.1 discussion.

The fits to the charged hadron transverse momentum distribution in **Figure 3** yields the hard scale temperature parameter *T* ¼ 0*:*72 GeV and *n* ¼ 3*:*1, in agreement with the extrapolation of (Eq. (9)) to 13 TeV *pp* collision energy, but with a smaller value of *n*. This reflecting the slower fall-off of the transverse momentum distribution at the LHC energy.

The integral of the area under the fit curves carries important information about entanglement in these and other in high energy physics processes. Defining the ratio *R* of the integral under the power law (hard scattering component) curve, *Ip* and the sum of the integrals of the exponential (thermal component) curve, *Ie* and power law curve of the fit in **Figure 3**:

#### **Figure 3.**

*Normalized transverse momentum distribution of charged hadrons from* ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> <sup>13</sup> *TeV pp collisions. The curves shown are exponential (red dashed) representing the thermal component of the distribution, and power law (green solid) corresponding to the hard scattering contribution. The superposition of these two contributions are also shown (blue, thin solid). Figure from [14], and data is from [32].*

*Quantum Information Science in High Energy Physics DOI: http://dx.doi.org/10.5772/intechopen.98577*

$$\mathbf{R} = \frac{\text{power}}{\text{power} + \text{exponential}} = \frac{\mathbf{I\_p}}{\mathbf{I\_p} + \mathbf{I\_e}} \tag{10}$$

The calculation yields the value of *R*≃0*:*16 � 0*:*05, in agreement (within the uncertainty interval) with the ratio calculated from the charged hadron spectra in inelastic proton-proton collisions at ISR energies of ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> 23, 31, 45, and 53 GeV [20] even given the large beam energy difference between the LHC and the ISR accelerators.

### **2.3 Diffractive events and di-muon pair transverse momentum distribution in proton-proton collisions**

Diffractive proton-proton (*pp*) collision events at the LHC can proceed through the photon-photon (*γγ*) interactions shown in (Eq. (11)). Both *X*<sup>0</sup> and *X*<sup>00</sup> can be final state protons from the collision, or the products *X*<sup>0</sup> ,*X*<sup>00</sup> of their diffractive dissociation (single diffraction (in which one of the incident protons dissociates into an inelastic state), and double diffraction (in which both of the incident protons dissociate)). Measurements from the ATLAS collaboration [33] of the reaction

$$pp(\eta\gamma) \to \mu^{+}\mu^{-}\mathbf{X}'\mathbf{X}''\tag{11}$$

at ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> 13 TeV center of mass energy in *pp* collisions are studied. Selection of the exclusive *γγ* ! *μ*þ*μ*� process was implemented by only including events that have both muon tracks (*μ*<sup>þ</sup> and *μ*�) while at the same time excluding events that show additional charged particle activity in the central region of the detector. Transverse momenta of greater than 400 MeV were used in the ATLAS analysis, with pseudorapidity range the same as that of the charged hadron analysis described in subSection 2.2. In the most recent ATLAS analysis of the reaction (11) only diffractive events that proceed through the *γγ* scattering were selected [33].

**Figure 4** shows the transverse momentum distribution in the case of *γγ* production of di-muon pairs in *pp* collisions at ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> 13 TeV center of mass energy. As can

**Figure 4.**

*The (normalized) transverse momentum diffractive scattering event distribution* <sup>1</sup> *pT dNμμ dpT in units of GeV*�<sup>2</sup> *for the reaction of (Eq. (11)) showing the absence of a thermal component to the distribution. The curve shown (green, solid) is the power law contribution corresponding to the hard-scattering process. Data from [33], figure from [14].*

be seen there, the hard scattering term alone describes well the distribution, and the thermal (exponential) component is absent. As discussed already in this chapter and in [13–15], diffractive events are expected to have a suppressed thermal (exponential) component due to the fact that in these diffractive processes the photon interacts coherently with the entire proton, and no entanglement entropy is expected. This was discussed in Section 2.1. As the presence of the thermal component in this approach is the consequence of the entanglement, we expect it to be absent in diffractive events, as confirmed in **Figure 4**. Furthermore, the ratio *R* defined in the previous section in this case is *R*≃1, in agreement with the theoretical expectations and the previous data for *γγ* scattering at OPAL at ffiffi *<sup>s</sup>* <sup>p</sup> <sup>¼</sup> 15 and 35 GeV that also show no thermal component. *R* is then equal to one within experimental uncertainty.
