Next Article in Journal
DA-LSTM-VAE: Dual-Stage Attention-Based LSTM-VAE for KPI Anomaly Detection
Next Article in Special Issue
Field Theory Approaches to Relativistic Hydrodynamics
Previous Article in Journal
Chaos-Based Lightweight Cryptographic Algorithm Design and FPGA Implementation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isotropization of a Rotating and Longitudinally Expanding ϕ4 Scalar System

by
Margaret E. Carrington
1,2,*,
Gabor Kunstatter
2,3,4,
Christopher D. Phillips
1,† and
Marcelo E. Rubio
1,2,‡
1
Department of Physics, Brandon University, Brandon, MB R7A 6A9, Canada
2
Winnipeg Institute for Theoretical Physics, Winnipeg, MB R3T 2N2, Canada
3
Department of Physics, University of Winnipeg, Winnipeg, MB R3M 2E9, Canada
4
Department of Physics, Simon Fraser University, Burnaby, BC V5A 1S6, Canada
*
Author to whom correspondence should be addressed.
Current address: Department of Electrical and Computer Engineering, University of Waterloo, Waterloo, ON N2L3G1, Canada.
Current address: SISSA, 34136 Trieste, Italy.
Entropy 2022, 24(11), 1612; https://doi.org/10.3390/e24111612
Submission received: 10 October 2022 / Revised: 1 November 2022 / Accepted: 2 November 2022 / Published: 5 November 2022
(This article belongs to the Special Issue Causal Relativistic Hydrodynamics for Viscous Fluids)

Abstract

:
We study numerically the evolution of an expanding system of scalar fields. The initial configuration is non-isotropic and rotating. We calculate the energy–momentum tensor and angular momentum vector of the system. We compare the time scales associated with the isotropization of the transverse and longitudinal pressures, and the decay of the initial angular momentum. We show that even a fairly large initial angular momentum decays significantly faster than the pressure anisotropy.

1. Introduction

In this paper, we study the time evolution of an expanding system of rotating massless real scalar fields with quartic coupling. Our calculation is based on the method developed in [1,2]. Observables calculated in a loop expansion exhibit divergences at next-to-leading order, which originate from instabilities in the classical solutions. The effect is seen in a calculation of the energy–momentum tensor at the next-to-leading order, where the energy density and pressures of the system diverge rapidly with increasing time. Gelis and his collaborators have shown that this problem can be cured using a resummation scheme that collects the leading secular terms at each order of an expansion in the coupling constant. This resummation can be performed by allowing the initial condition for the classical field to fluctuate, and averaging over these fluctuations. They have shown that a system of scalar fields isotropizes when this resummation is performed [2].
The motivation behind the development of this approach is to study the thermalization of the glasma phase of the matter created in a relativistic heavy ion collision. It is known that a hydrodynamic description, which is valid when the system is fairly close to thermal equilibrium, works well at very early times (∼1 fm/c). Approaches that are based on kinetic theory descriptions of the scattering of quasi-particles cannot explain this rapid thermalization. Another possibility that has been studied extensively is that the system is strongly coupled, even at very high energies. The proposal of Gelis et al. is that rapid thermalization could be achieved by a resummation of quantum fluctuations. The color glass condensate (CGC) effective theory provides a natural framework for this formulation [3,4,5]. At very early times, the system is best described as a system of strong classical fields that can be obtained from solutions of the Yang–Mills equation using a CGC approach. The spectrum of quantum fluctuations was derived in [6]. The success of the resummation method was demonstrated in [7], where the authors showed that pressure isotropiztion occurs in an SU(2) analog of QCD.
Our ultimate goal is to use the Gelis et al. approach to study the creation and evolution of angular momentum in a glasma. This is interesting in the context of recent proposals that the glasma is produced in a rapidly rotating state, which could be detected by looking for the polarization of produced hyperons. There have been calculations that predict very large values for the initial angular momentum of the system [8,9,10], but significant final state polarization effects have not been observed [11,12]. In this paper, we develop a formulation to calculate the angular momentum of a system of real scalar fields. We present preliminary results that indicate the angular momentum relaxes to a small value on a time scale significantly smaller than the time scale for pressure isotropization. If a similar result is obtained in a QCD glasma, it would be consistent with the observations in [11,12]. We also comment that a calculation of angular momentum in glasma was conducted in [13], using a CGC approach with a proper time expansion, and it was found also that large amounts of angular momentum were not produced.
Since computations in a gauge theory are considerably more complicated, we will work with a scalar theory. While it is true that QCD and scalar ϕ 4 theory are different in many ways, they have important similarities in the context of this calculation because they both have unstable modes and are scale invariant at the classical level. In addition, we will mimic the kinematics of a relativistic nuclear collision by working in Milne coordinates with a rapidity independent background field. Milne coordinates are suitable because in a nuclear collision, there is a preferred spatial direction provided by the collision axis, and in the high energy limit, one expects invariance under Lorentz boosts in the z-direction.
This paper is organized as follows. In Section 2, we describe the method, and in Section 3, we formulate the calculation of the energy–momentum tensor and angular momentum. Some details of our numerical procedure are discussed in Section 4. In Section 5, we present our results, and in Section 6 we make some concluding remarks.
Throughout this paper, the spacetime is always taken to be Minkowski, with the signature ( + , , , ) . In addition to standard inertial coordinates ( t , x , y , z ) , we also use Milne coordinates ( τ , x , y , η ) , where τ is the proper time and η is the spacetime rapidity. Finally, we choose units such that c = k B = = 1 , where c is the speed of light in a vacuum, k B is the Boltzmann constant, and is the Planck constant divided by 2 π .

2. Formalism

2.1. Preliminaries

We consider a massless self-interacting real scalar field ϕ with quartic coupling. The Lagrangian density is given by
L = 1 2 μ ϕ μ ϕ g 2 4 ! ϕ 4
where g is the coupling constant. To mimic the kinematics of a high energy nuclear collision, we work in Milne coordinates ( τ , η , x ) with
τ = t 2 z 2 η = 1 2 ln t + z t z .
Under a Lorentz boost in the z-direction, the proper time is unchanged and η is shifted by a constant. The metric in Milne coordinates is
g μ ν d x μ d x ν = d τ 2 τ 2 d η 2 d x 2 d y 2 .
Figure 1 shows the curves of constant τ and η .

2.2. The Resummation Procedure

As explained in [1,2], observables calculated in a loop expansion exhibit secular divergences at next-to-leading order that originate from instabilities of the classical solutions. Gelis et al. proposed to cure this problem using a resummation scheme that collects the leading secular terms at each order of an expansion in the coupling constant, by averaging over an ensemble of initial conditions. The energy–momentum tensor is ultraviolet divergent, but the divergence corresponds to a vacuum contribution and can be removed by repeating the calculation with the background field set to zero, and subtracting the results. This vacuum subtraction was performed for all the calculations presented in this paper.
The equation of motion for the scalar field obtained from the Lagrangian (1) is
ϕ ¨ ( τ , η , x ) 1 τ ϕ ˙ 1 τ 2 η 2 ϕ Δ ϕ + g 2 6 ϕ 3 = 0
where the “dot” indicates a derivative with respect to τ , and Δ is the transverse Laplacian operator. The initial field is written as the sum of a background field contribution, φ , which is assumed to be boost invariant and therefore independent of η , and an η -dependent fluctuation, which we call α
ϕ χ ( τ 0 , η , x ) = φ ( τ 0 , x ) + α χ ( τ 0 , η , x ) .
The initial time τ 0 is chosen to be small but nonzero (see Section 4.3 for further discussion). The initial background field φ ( τ 0 , x ) is discussed in Section 4.5. The index χ in Equation (4) indicates that we have a Gaussian ensemble of initial conditions defined as
α χ ( τ 0 , η , x ) = d K c K χ a K + c K χ * a K * .
The index K labels the momentum variables ( ν , k ) that are conjugate to the coordinate-space variables ( η , x ) , respectively. The notation c K χ indicates an element in a Gaussian-distributed ensemble of N χ random numbers, with variance
c K * c L = 1 2 δ K L .
We use the momentum space integration measure
d K = d ν 2 π d k ( 2 π ) 2
and the delta function in Equation (6) is defined so that d K δ K L = 1 . The mode functions a K a ν k ( τ 0 , η , x ) are obtained from the linearized equations of motion
a ¨ K + 1 τ a ˙ K 1 τ 2 η 2 a K Δ a K + g 2 2 φ 2 ( τ 0 , x ) a K = 0
and normalized so that d K ( a K , a L ) = 1 with
( a K , a L ) = i τ d η d 2 x a K * τ a L ( τ a K * ) a L .
Separating variables and performing the normalization, one finds
a K a ν k ( τ 0 , η , x ) = 1 2 π e π ν / 2 e i ν η χ k ( x ) H i ν ( 2 ) ( λ k τ 0 )
where the χ k is the solution of the eigenvalue equation
Δ + g 2 2 φ 2 ( τ 0 , x ) χ k ( x ) = λ k 2 χ k ( x ) .
The field ϕ χ ( τ , η , x ) at finite proper time is obtained by solving Equation (3) with the initial condition ϕ χ ( τ 0 , η , x ) obtained from Equations (4), (5)–(7), (10) and (11). From this point on, we drop the subscript χ .

3. Observables

3.1. Energy Momentum Tensor

The energy–momentum tensor of theory (1) is
T μ ν = μ ϕ ν ϕ g μ ν 1 2 α ϕ α ϕ g 2 4 ! ϕ 4 .
The invariance of the Lagrangian under the conformal transformation
g μ ν Ω 2 g μ ν ; ϕ Ω 1 ϕ
implies that T μ ν is traceless on shell.
The expressions for the energy and pressure are
ϵ = T 00 = 1 2 ( τ ϕ ) 2 + ( η ϕ ) 2 τ 2 + ( x ϕ ) 2 + ( y ϕ ) 2 + V ( ϕ ) p L = τ 2 T 11 = 1 2 ( τ ϕ ) 2 + ( η ϕ ) 2 τ 2 ( x ϕ ) 2 ( y ϕ ) 2 V ( ϕ ) p T = 1 2 T 22 + T 33 = 1 2 ( τ ϕ ) 2 ( η ϕ ) 2 τ 2 V ( ϕ )
where V ( ϕ ) = g 2 ϕ 4 / 4 ! . In terms of the energy and pressure, the trace condition is
ϵ = 2 p T + p L .

3.2. Angular Momentum

We use the standard Pauli–Lubanski formalism [14,15] to obtain an expression for the angular momentum in terms of the energy–momentum tensor. We define the tensor field
M μ ν λ = T μ ν R λ T μ λ R ν
where R μ is the coordinate vector. Using Stokes’ theorem, one obtains a set of six conserved quantities
J ν λ = Σ d 3 y | γ | n μ M μ ν λ ,
where n μ is a unit vector perpendicular to the hypersurface Σ , γ i j is the induced metric on this hypersurface, and d 3 y is the corresponding volume element. The angular momentum is obtained from the Pauli–Lubanski vector
L μ = 1 2 ϵ μ α β ρ J α β u ρ
where u ρ is the vector that denotes the rest frame of the system. Equations (16)–(18) will give
L μ = 1 2 ϵ μ α β ρ d 3 y γ n σ u ρ ( T σ α R β T σ β R α )
where the energy–momentum tensor is given in Equation (12).
To find the angular momentum on a surface of constant τ , we define
n μ = τ x μ .
In Minkowski coordinates, this gives n μ = ( cos h ( η ) , 0 , 0 , sin h ( η ) ) , and it is easy to verify that n μ Milne = ( 1 , 0 , 0 , 0 ) , as expected. The fluid velocity is the local rest frame in comoving coordinates, which is written u Milne ρ = ( 1 , 0 , 0 , 0 ) . In Minkowski coordinates, this becomes u ρ = ( cos h ( η ) , 0 , 0 , sin h ( η ) ) . We could calculate the angular momentum directly in Minkowski coordinates, or alternatively, we could perform the calculation in Milne coordinates and perform a coordinate transformation to obtain the Minkowski space result. We checked our computations by verifying that both calculations give the same result. The components of the angular momenta about each of the Minkowski coordinate axes are
L t = τ d 2 x d η sin h ( η ) ϕ ˙ x y ϕ y x ϕ L x = d 2 x d η ϕ ˙ y η ϕ L y = d 2 x d η ϕ ˙ x η ϕ L z = τ d 2 x d η cos h ( η ) ϕ ˙ y x ϕ x y ϕ .
We note that all components of the angular momentum are dimensionless (in natural units, with = 1 ).

4. Numerical Implementation

4.1. Lattice Discretization

We discretize in both directions in the transverse plane with L grid points and lattice spacing set to 1, which effectively means that we define all dimensionful quantities in terms of the transverse lattice grid spacing. The rapidity variable η is discretized with N grid points and lattice spacing h. We consider a unit slice of rapidity, and therefore take h = 1 / N .
The discretization of the transverse variables is straightforward. The discretized version of Equation (11) is
D i j ; k l χ k l = λ 2 χ i j
with
D i j ; k l = ( 4 + V i j ) δ i k δ j l ( δ i + 1 k + δ i 1 k ) δ j l δ i k ( δ j + 1 l + δ j 1 l ) .
Since D is a rank 4 tensor with L 4 components, we obtain L 2 eigenfunctions χ i j e , and L eigenvalues ( λ 2 ) e , with e ( 1 , L 2 ) . The normalized eigenfunctions are
i j χ i j * e χ i j e ¯ = L 2 δ e e ¯
and the momentum integration is discretized as
d 2 k ( 2 π ) 2 1 L 2 e = 1 L 2 .
Since the spatial lattice spacing is set to 1, an integral over transverse coordinates is discretized as
d 2 x i = 0 L 1 j = 0 L 1 .
The discretization of the longitudinal variables is a little more subtle. The constraint
η 2 e i ν η = ν 2 e i ν η
gives
ε v 2 : = ν 2 = 2 h sin π v N 2
and we replace ν ε v in every factor e π ν / 2 and in the Hankel functions. For the complex exponential, we use e i ν η e 2 π i v n N . The integral over ν becomes a sum over v using
d ν 2 π 1 N h v = 0 N 1 .
Combining these expressions, we find the discretized versions of Equations (4), (6) and (10):
α n i j ( τ ) = 1 N L 2 h v = 0 N 1 p = 1 L 2 c v p a n i j v p ( τ ) + c . c . a n i j v p ( τ ) = 1 2 π e 2 π i v n N χ i j p e π ν / 2 H i ν ( 2 ) ( λ k τ ) c v e c u e ˜ * = 1 2 N L 2 h δ v u δ e e ˜ .
To verify that discretization is performed correctly, we checked the discretized version of the normalization condition (9).

4.2. Boundary Conditions

We use periodic boundary conditions, which means that the indices ( i , j ) that correspond to the transverse spatial coordinates are defined as modulo L, and the index n for the rapidity is modulo N. The boundary conditions satisfy the self-adjointness condition
F ϕ ( x ) = ϕ ( i + 1 ) ϕ ( i ) B ϕ ( x ) = ϕ ( i ) ϕ ( i 1 ) i f ( i ) F g ( i ) = i B f ( i ) g ( i ) .

4.3. Hankel Functions

The differential equation for the mode function was solved by separating variables, which gives the solution in (10). The time-dependent part of the equation is of the second order, and has two independent solutions, which are the Hankel functions H i ν ( 1 ) ( λ τ ) and H i ν ( 2 ) ( λ τ ) . We use only the second because it has positive frequency behavior at large times
lim τ H i ν ( 2 ) ( τ ) = 2 π τ e i ( τ i π ν / 2 π / 4 ) .
From now on, we suppress the superscript ( 2 ) on the Hankel function. When τ 0 , the Hankel function oscillates like e ± i τ ν and the derivative diverges. Numerically, we must start the evolution at a small positive time, which we choose as τ 0 = 10 2 . One can check that the value chosen for this small initial time does not change the results at finite times.
We describe below our method to calculate the Hankel functions. First, we define the scaled function
h i ν ( λ τ ) = e π ν / 2 H i ν ( λ τ )
which is easier to calculate numerically. At large times, one can obtain the scaled Hankel function for given values of ν and λ from the asymptotic series
h i ν ( λ τ ) = 2 π λ τ e i ( λ τ π / 4 ) k = 0 n t k + O ( τ ( n + 1 ) ) t k = ( 1 ) k k ! ( 2 i λ τ ) k s = 1 k ν 2 + ( 2 s 2 ) 2 4 .
This expression must be used carefully because the series does not converge for arbitrarily large values n. We proceed as follows. For a given value of ν and λ , choose some value of τ and look for a value of k max so that t k max + 1 < 10 9 and Max ( t k k max ) < 10 6 . If this k max can be found, use Equation (33) with n = k max . If k max does not exist, then increase the chosen value of τ and try again. Using this procedure, we can find h i ν ( λ τ ) and its first derivative for each value of ν and λ , for some (possibly very large) time. We then use adaptive fifth-order Runge–Kutta to find each Hankel function at the initial time τ 0 .

4.4. Discretized Derivatives

The conservation equation
ϵ τ + ϵ + p L τ = 0
is an exact equation that should be satisfied whether or not the system is in equilibrium. Additionally, we should have that the trace of the energy–momentum tensor is zero, so that Equation (15) is satisfied. It is easy to show analytically that these conditions are satisfied for background fields if we use forward derivatives: x f ( x ) f ( i + 1 ) f ( i ) . The point is that while centered derivatives are not wrong, much larger lattices must be used to achieve the same numerical accuracy.
For angular momentum, the situation is different. All contributions to the angular momentum have an integral of the form d x φ ˙ x ϕ . If the initial value of φ ˙ is constant, the integrand is a total derivative and therefore the integral will give zero. However, this is not well satisfied numerically with forward derivatives. In the calculation of angular momentum, it is therefore better to use centered derivatives: x f ( x ) ( f ( i + 1 ) f ( i 1 ) ) / 2 .

4.5. Initial Conditions

The initial conditions that we use for the background field and its derivative are
φ ( τ 0 , i , j ) = φ 0 cos ( k x i + k y j ) φ ˙ ( τ 0 , i , j ) = φ ˙ 0 sin i L + 1 2 π L 1 .
The argument of the sine function is π / 2 at i = 1 and i = L , and zero at i = ( L + 1 ) / 2 , so the field has negative φ ˙ 0 on the left side of the lattice and positive φ ˙ 0 on the right side.
The astute reader will note that our initial classical field is not periodic, and therefore does not respect our boundary conditions. The reason is that we wish to avoid problems that may arise when resonant modes are considered, which in the present model would correspond to the normal modes of the finite spatial lattice. For a large enough lattice, all modes are effectively periodic, and it is therefore expected that the precise form of the initialization is not important.

5. Results and Discussion

All of our results are obtained with L = 41 spatial grid points, N = 120 points for the rapidity coordinate, and N χ = 256 configurations. The initial conditions for the background field are obtained from (35) with φ 0 = 15 , k x = k y = 1 / 2 and φ ˙ 0 = 10 .
To investigate if the system obeys Equation (15), we compare the energy density and the sum of the pressures. This is shown in Figure 2. One sees that after some initial oscillations have damped out, the condition ϵ = 2 p T + p L is well satisfied.
To see if the system approaches an isotropic state, and if it obeys an equation of state, we look at the transverse and longitudinal pressures. The left panel of Figure 3 shows that, after some initial oscillations have disappeared, the transverse and longitudinal pressures approach each other up to a time of about τ 160 . The right panel shows the two pressures normalized by the energy density, both approaching 1/3, again up to τ 160 . For large times, the simulation breaks down, which is not unexpected when one studies the dynamics of an expanding system inside a box of finite size.
In Figure 4, we show the three components of the angular momentum in Equation (21). The z component, which depends weakly on the rapidity, is averaged over the unit slice of rapidity that we consider. In comparison with the energy and pressure, the oscillatory behavior is more severe and does not completely disappear. To obtain a better idea of the overall behavior, we also plot the accumulated average for each component, which is shown in Figure 4 with the thick lines. In each case, the darker color corresponds to the average of the component with the same but lighter color. The figure shows that even a fairly large initial angular momentum decays very quickly.
We want to compare the time scales for the isotropization of the pressures, and the decay of the initial angular momentum. In Figure 5, we show in blue the curve in the left panel of Figure 3 over the range of τ for which the decay is strongest. To produce the light green points, we took the data for | L | versus τ with τ > 12.0 , where the large initial fluctuations are mostly gone, and shifted the first point (which was (12.0, 9.10)) so that it sits on top of the first point of the data that made the blue curve. The dark green line is a fit obtained for these data using the function A + B / τ + C e D τ . The plot shows clearly that the initial angular momentum decays much more quickly than the pressure anisotropy, although the dispersion in the data is large. To quantify this dispersion, we calculated
σ = 1 N i = 1 N L i L fit ( τ i ) 2 L i + L fit ( τ i ) 2
where the sum is over the points shown in green in Figure 5 and L fit ( τ ) is the fitted function that gives the green curve. The result is σ = 0.22 .

6. Conclusions

In this paper, we presented some preliminary results from our study of the angular momentum in an expanding system of rotating massless scalar fields. Our results indicate that even when a large amount of angular momentum is put into the system, it decays very rapidly. Future work will include an investigation of how much these results depend on the exact form of the initialization and the boundary conditions, and possibly the extension of the calculation to physical theories, such as QCD.

Author Contributions

Conceptualization, M.E.C., G.K., C.D.P. and M.E.R.; formal analysis, M.E.C., G.K., C.D.P. and M.E.R.; writing—original draft, M.E.C., G.K., C.D.P. and M.E.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Natural Sciences and Engineering Research Council of Canada Discovery Grant program from grants 2017-00028 and 2018-04090. This research was enabled in part by support provided by WestGrid (www.westgrid.ca, accessed on 1 November 2022) and the Digital Research Alliance of Canada (alliancecan.ca, accessed on 1 November 2022).

Data Availability Statement

Not applicable.

Acknowledgments

Margaret Carrington gratefully acknowledges the helpful discussions with François Gelis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dusling, K.; Epelbaum, T.; Gelis, F.; Venugopalan, R. Role of quantum fluctuations in a system with strong fields: Onset of hydrodynamical flow. Nucl. Phys. A 2011, 850, 69–109. [Google Scholar] [CrossRef] [Green Version]
  2. Dusling, K.; Epelbaum, T.; Gelis, F.; Venugopalan, R. Instability induced pressure isotropization in a longitudinally expanding system. Phys. Rev. D 2012, 86, 085040. [Google Scholar] [CrossRef] [Green Version]
  3. McLerran, L.; Venugopalan, R. Computing quark and gluon distribution functions for very large nuclei. Phys. Rev. D 1994, 49, 2233–2241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. McLerran, L.; Venugopalan, R. Gluon distribution functions for very large nuclei at small transverse momentum. Phys. Rev. D 1994, 49, 3352–3355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. McLerran, L.; Venugopalan, R. Green’s function in the color field of a large nucleus. Phys. Rev. D 1994, 50, 2225–2233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Epelbaum, T.; Gelis, F. Fluctuations of the initial color fields in high energy heavy ion collisions. Phys. Rev. D 2013, 88, 085015. [Google Scholar] [CrossRef] [Green Version]
  7. Epelbaum, T.; Gelis, F. Pressure isotropization in high energy heavy ion collisions. Phys. Rev. Lett. 2013, 111, 232301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gao, J.H.; Chen, S.W.; Deng, W.T.; Liang, Z.T.; Wang, Q.; Wang, X.N. Global quark polarization in noncentral A+A collisions. Phys. Rev. C 2008, 77, 044902. [Google Scholar] [CrossRef] [Green Version]
  9. Becattini, F.; Piccinini, F.; Rizzo, J. Angular momentum conservation in heavy ion collisions at very high energy. Phys. Rev. C 2008, 77, 024906. [Google Scholar] [CrossRef] [Green Version]
  10. Liang, Z.T. Spin Effects In Heavy Ion Collisions at High Energies. In Proceedings of the 24th International Symposium on Spin Physics, Matsue, Japan, 18–22 October 2021. [Google Scholar]
  11. Adam, J.; Adamczyk, L.; Adams, J.R.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Ajitan, N.N.; Alekseev, I.; Anderson, D.M.; et al. Global polarization of Λ hyperons in Au-Au collisions at sNN=200 GeV. Phys. Rev. C 2018, 98, 014910. [Google Scholar] [CrossRef]
  12. Acharya, S.; Adamova, D.; Adhya, S.P.; Adler, A.; Adolfsson, J.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agrawal, N.; Ahammed, Z.; et al. Global polarization of Λ and Λ¯ hyperons in Pb-Pb collisions at the LHC. Phys. Rev. C 2020, 101, 044611. [Google Scholar] [CrossRef] [Green Version]
  13. Carrington, M.E.; Czajka, A.; Mrówczyński, S. Physical characteristics of glasma from the earliest stage of relativistic heavy ion collisions. Phys. Rev. C 2022, 106, 034904. [Google Scholar] [CrossRef]
  14. Lubański, J. Sur la theorie des particules élémentaires de spin quelconque. I. Physica 1942, 9, 310–324. [Google Scholar] [CrossRef]
  15. Lubanski, J. Sur la théorie des particules élémentaires de spin quelconque. II. Physica 1942, 9, 325–338. [Google Scholar] [CrossRef]
Figure 1. Representation of hypersurfaces of constant τ and η . The rapidity determines the location of a particle along a surface of fixed τ .
Figure 1. Representation of hypersurfaces of constant τ and η . The rapidity determines the location of a particle along a surface of fixed τ .
Entropy 24 01612 g001
Figure 2. The energy and sum of the pressures as functions of τ .
Figure 2. The energy and sum of the pressures as functions of τ .
Entropy 24 01612 g002
Figure 3. The transverse and longitudinal pressures, normalized by the energy density.
Figure 3. The transverse and longitudinal pressures, normalized by the energy density.
Entropy 24 01612 g003
Figure 4. The three components of the angular momentum vector and their accumulated averages.
Figure 4. The three components of the angular momentum vector and their accumulated averages.
Entropy 24 01612 g004
Figure 5. A comparison of ( p T p L ) / ϵ (blue) and | L | (green); see text for details.
Figure 5. A comparison of ( p T p L ) / ϵ (blue) and | L | (green); see text for details.
Entropy 24 01612 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Carrington, M.E.; Kunstatter, G.; Phillips, C.D.; Rubio, M.E. Isotropization of a Rotating and Longitudinally Expanding ϕ4 Scalar System. Entropy 2022, 24, 1612. https://doi.org/10.3390/e24111612

AMA Style

Carrington ME, Kunstatter G, Phillips CD, Rubio ME. Isotropization of a Rotating and Longitudinally Expanding ϕ4 Scalar System. Entropy. 2022; 24(11):1612. https://doi.org/10.3390/e24111612

Chicago/Turabian Style

Carrington, Margaret E., Gabor Kunstatter, Christopher D. Phillips, and Marcelo E. Rubio. 2022. "Isotropization of a Rotating and Longitudinally Expanding ϕ4 Scalar System" Entropy 24, no. 11: 1612. https://doi.org/10.3390/e24111612

APA Style

Carrington, M. E., Kunstatter, G., Phillips, C. D., & Rubio, M. E. (2022). Isotropization of a Rotating and Longitudinally Expanding ϕ4 Scalar System. Entropy, 24(11), 1612. https://doi.org/10.3390/e24111612

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop