Correlation of density fluctuation in a magnetized QCD matter near the critical end point

The dynamical correlation of density fluctuation in quark gluon plasma with a critical end point has been investigated within the scope of the M\"uller-Israel-Stewart theory in the presence of static ultra-high external magnetic field. The dynamic structure factor of the density fluctuation exhibits three Lorentzian peaks in absence of external magnetic field- a central Rayleigh peak and two Brillouin peaks situated symmetrically on the opposite sides of the Rayleigh peak. The spectral structure displays five peaks in presence of the magnetic field due to the coupling of the magnetic field with the hydrodynamic fields in second-order hydrodynamics. The emergence of the extra peaks is due to the asymmetry in the pressure gradient caused by the external magnetic field in the system. Interestingly, it is observed that near the critical end point, all the Brillouin peaks disappear irrespective of the presence or absence of the external magnetic field.


I. INTRODUCTION
One of the primary objectives of heavy-ion collision (HIC) experiments at Relativistic Heavy Ion Collider (RHIC) and at Large Hadron Collider (LHC) is to create and characterize a deconfined state of thermal quarks and gluons, called quark gluon plasma (QGP).It is expected that after a proper time τ 0 ∼ 1f m/c of the collision, the system attains a state of local thermal equilibrium.The evolution of the QGP in spacetime can be studied by using relativistic hydrodynamics [1][2][3][4], which is a low frequency or long wavelength effective theory of many body interacting systems.In the non-central collisions of nuclei at RHIC and LHC energies large electric and magnetic (EM) fields will be produced due to the electric current generated by the accelerated motion of charged spectators, i.e. protons.In non-central collisions of heavy nuclei (Au+Au or Pb+Pb) at RHIC and LHC energies, the transient magnetic field (B) can be as high as ∼ 10 17 -10 18 Gauss) [5][6][7][8][9].Therefore, in such a situation, it is imperative to consider the effect of time varying magnetic and electric fields on the characterization of the QGP.The survival time of the fields crucially depends on the electrical conductivity of the QGP [10][11][12].We theoretically evaluate the dynamic structure factor [13] of the system in the present work in an idealistic scenario of non-expanding QGP with critical point under a static magnetic field.The inclusion of magnetic field in relativistic hydrodynamics is studied within the ambit of relativistic magnetohydrodynamics (MHD) which is a self-consistent macroscopic framework that deals with the evolution of mutually interacting charged fluid along with the EM fields.Recently, several authors have studied the effect of the EM fields on QGP fluid in the context of special relativistic systems [14][15][16][17][18][19][20].
In the relativistic viscous hydrodynamics and relativistic magnetohydrodynamics, the transport coefficients, such as the shear viscosity, bulk viscosity, thermal conductivity, etc. are taken as inputs which can be estimated from an underlying microscopic theory [21][22][23][24].A straightforward extension of the non-relativistic viscous fluid dynamics (Navier-Stokes equation) to relativistic regime (without magnetic field) [25,26] leads to acausal and unstable solutions [27][28][29].These issues were addressed and resolved by Müller, and Israel and Stewart (MIS) [30,31] who developed a causal and stable second-order relativistic hydrodynamics.Here the order of the theory is dictated by different orders of gradients in the expansion hydrodynamic quantities i.e., energy-momentum tensor (EMT).
The dynamic density correlator or the structure factor, S nn (k, ω), has been studied ex-tensively in condensed matter physics by using linear response theory [13] within the ambit of fluid dynamics.The correlation of density fluctuations can be examined through the dynamic spectral structure S nn (k, ω) in wave vector (k), frequency (ω) space.Ordinarily, for a fluid modelled by relativistic hydrodynamics, the S nn (k, ω) consists of three distinct peaks, two of them are Brillouin (B) peaks, created due to pressure fluctuation at constant entropy and the other is the Rayleigh (R) peak, originates due to thermal or entropy fluctuation at constant pressure [32,33].
The spectral structure has been extensively investigated experimentally in condensed matter physics to determine the speed of sound by using the scattering of photons and neutrons.It is observed that the position of the B-Peaks in S nn (k, ω) depends on the speed of sound and the width of both B and R peaks enable us to evaluate various transport coefficients (shear and bulk viscosities, thermal conductivity), and thermodynamic response functions (specific heats).For QCD (quantum chromodynamics) matter, no such external probes are avaiable making the direct detection of the peaks extremely difficult, if not impossible.However, such investigation may shed light on the speed of sound waves, and consequently on the equation of state (EoS) of the QCD matter.In the previous studies [32,33], the relativistic hydrodynamics is used to find out the correlation in density fluctuation in QCD matter without including magnetic field.In the present work, we have considered the relativistic MHD to study the dynamical density correlator in a baryon-rich fluid.We derive S nn (k, ω), which allows us to determine the speed of the perturbation propagating as sound waves through the QCD medium immersed in a constant magnetic field.
The manuscript is organized as follows: In Sec.II, we discuss briefly the form of energymomentum tensor of fluid in the presence of EM fields.In the next section i.e. in Sec.III, we find the linearized equations of MIS hydrodynamics and derive the density-density correlation function.In Sec IV, we discuss results that stem from our analysis, and the summary and conclusion is presented in Sec.V.

II. RELATIVISTIC MAGNETOHYDRODYNAMICS
In this section, first we discuss the equation of motion of electromagnetic (EM) field and magnetohydrodynamics.

A. Equation of motion of the EM field
We start by discussing the relativistically covariant formulation of classical electrodynamics.The second rank anti-symmetric EM field tensor F µν , can be defined as, in terms of electric field four-vector E µ (= F µν u ν ), magnetic field four-vectors B µ = ( 1 2 ϵ µναβ u ν F αβ = F µν u ν ), and the fluid four-velocity u µ [34][35][36], and its dual counter part is represented by, where, ϵ µναβ is the Levi-Civita tensor.It is quite straight forward to see by using the antisymmetric property of F µν that both E µ and B µ are orthogonal to u µ i.e., E µ u µ = B µ u µ = 0.
It may be noted that in the rest frame u µ = (1, 0), E µ = (0, E), and B µ = (0, B), where E and B correspond to the electric and the magnetic three vectors fields with E i = F i0 and B i = − 1 2 ϵ ijk F jk .The indices i, j, k run over 1, 2, 3.The E µ and B µ can be interpreted as the electric and magnetic fields respectively measured in a frame in which the fluid moves with velocity u µ .
Using the EM field tensor and its dual, we can express the Maxwell's equations in a covariant way as, where, J ν represents the electric charge four-current and it acts as the source of EM field.
The electric charge four-current (J µ ) can be decomposed as follows: where, j µ is the conduction current and d µ = ∆ µ ν J ν , is the charge diffusion current, n q = u µ J µ the proper net charge density.If we consider a linear relation between j µ and E µ (Ohm's law) then we can write j µ = σ µν E ν , where σ µν is the second rank conductivity tensor.The construction of j µ follows u µ j µ = 0.It implies that the conduction current exists even in the absence of any net charge.The solutions of Eqs. ( 3) and (4) for given electric charge four current J µ in Eq. ( 5) completely determine the evolution of electromagnetic field.It also acts as a coupling between the fluid and EM fields as it contains the fluid information e.g.fluid conductivity σ µν , net charge density n q etc.We consider a system consists of quarks, anti-quarks and gluons.All the relevant thermodynamic quantities and the EoS used in this work have been calculated by taking these degrees of freedom into account (we do not repeat the discussions on EoS here but refer to Ref. [37] for details).Therefore, the net (quark-antiquark) charge density is equivalent to net (quark-antiquark) number density and the relation n q = qn f holds, where n f corresponds to net number density.
At first, we assume that the fluid does not possess any polarization or magnetization.In that case the EM field stress-energy tensor (EMT) can be written as, If we take the partial derivative of the field stress-energy tensor, we obtain the equation of motions, In the above equation the current density due to external source is ignored.In presence of an external source (J µ ext ), the total current is given by, In such case, the external current acts as a source term in energy-momentum conservation equation.
In this article, we consider an ideal magneto-hydrodynamic limit which resembles very large magnetic Reynolds number R m = Lvσµ >> 1, where L is the characteristic macroscopic length or time scale, v is the characteristic velocity of the flow, σ is the isotropic electrical conductivity (σ µν = σg µν ), and µ is the magnetic permeability of the QGP.It is evident that R m increases with σ.From the induced current density equation, J µ ind = j µ = σE µ , it is clear that when σ attains a very large value (→ ∞) then the electric field will approach zero (E µ → 0) to keep J µ ind but finite.This simplifies the EM tensor F µν to the following form, Using Eqs. ( 8) and ( 9) in the Maxwell's equations Eq. ( 3) we obtain, The EMT in absence of electric field can be obtained from Eqs. ( 6) and ( 9) as, where Using Eq. ( 9), one can show that B µν B µν = 2B 2 .We can define another anti-symmetric tensor as:

B. Equation of motion of magnetohydrodynamics
In this section, we derive the equation of motion for relativistic fluid under the influence of external EM field.

Conservation of energy and momentum of fluid and electromagnetic field
In the absence of any magnetic field, the EMT and the particle currents are conserved separately according to the following conservation laws, In presence of external EM field the total EMT tensor (fluid+field), T µν is given by, where T µν f and T µν EM are the contributions from the fluid and EM field respectively.The total EMT in Ref. [35] contains additional terms which can not be unambiguously attributed to the fluid or to the field.But for constant susceptibility and vanishing E µ such contributions vanish and then Eq. ( 14) becomes a good approximation.As the electric charge is conserved, the charge current of the fluid is individually conserved too, If we have an external charge current, then it will act as a source term of the EMT: The conservation equation for the electromagnetic field Eq. ( 7) with external source can be written as, Using Eqs ( 14), ( 16) and ( 17) we get, Usually, the total EMT of an isolated system remains conserved but in the presence of an external charge current an appropriate source term should be taken into account.In this case, the fluid evolution depends on the fluid charge current density through Eq. (18).
It is also useful to express the conservation equations by taking projection along and perpendicular to fluid four velocity in an alternative way.If we take the parallel projection of Eq. ( 17) and Eq. ( 18), then we obtain, In presence of current density given by Eq. ( 8), the perpendicular projection of Eqs.(17) and ( 18) can be written as: It implies that the momentum density of fluid depends on diffusion current/magnetic field, momentum density of the field, external current and fluid diffusion current.

Ideal and dissipative non-resistive magnetohydrodynamics
The EMT for non resistive magnetohydrodynamics is given by, In case the ideal fluid, the total EMT takes the form- where, ϵ and P denote energy density and the thermodynamic pressure respectively.For dissipative fluid with non-zero shear and bulk viscosities, and non-zero thermal conductivity, the EMT becomes- where, Π, q µ and π µν denote bulk pressure, heat flux and shear stress respectively, and are expressed in Eckart's frame of reference [25,31].The system of equations can be closed with the help of constitutive relation of charged-current and with an Equation of State (EoS), which relates the thermodynamic pressure to energy and number density P = P (ϵ, n f ).
The inclusion of the polarization modifies the energy-momentum tensor as follows [14,38,39]: where, M µν is the polarization tensor.In the non-dissipative limit, the entropy is also conserved along with the charge.The charge and entropy currents can be expressed in the non-dissipative hydrodynamics as, where, n and s are the electric charge density and entropy density measured in the local rest frame.It is convenient to decompose the tensor F µν into components parallel and perpendicular to u µ as: The anti-symmetric polarization tensor M µν represents the response of matter to F µν .It is given by M µν ≡ −∂Φ/∂F µν , where Φ is the thermodynamic potential function.It is also convenient to define the in-medium field strength tensor H µν ≡ F µν − M µν .Similar to F µν , we decompose M µν and H µν as follows: with There are several physical system where the electric field is much weaker than the magnetic field.The interior of a neutron star is one such example.In the following discussions, we will omit the contribution from the electric field.We introduced the four-vector In the absence of electric field, we have, where M ≡ |M| and H ≡ |H|.In absence of electric fields, the matter and field contributions to the EMT (26) can now be written in terms of b µ and b µν as (see, e.g., Refs.[14,40]) where Ξ µν ≡ ∆ µν + b µ b ν is a new projection tensor with Ξ µν u µ = Ξ µν b µ = 0.The transverse and longitudinal pressures relative to b µ can be defined as P ⊥ = P − M B and P ∥ = P relative to the vector b µ .In the absence of magnetic field, the fluid is isotropic and P ⊥ = P ∥ = P , where P is the thermodynamic pressure defined in Eq. ( 28).In the local rest frame of the fluid, the direction of the magnetic field is chosen as the z-axis without loss of generality, so we have b µ = (0, 0, 0, 1).Then the EMT takes the form, T µν F0 = diag(ϵ, P ⊥ , P ⊥ , P ∥ ).In the presence of polarization (or non-zero magnetization), the full energy momentum tensor (EMT) for relativistic MHD can be expressed as, The form of q µ , Π, π µν used in MIS hydrodynamics contains additional coupling and relaxation coefficients arising due to inclusion of second order gradients [27,41,42]: where, ζ, η, and κ are the coefficient of bulk viscosity, shear viscosity, and thermal conductivity respectively, and Here, β 0 , β 1 , β 2 are relaxation coefficients, α 0 and α 1 are coupling coefficients.The relaxation times for the bulk pressure (τ Π ), the heat flux (τ q ) and the shear tensor (τ π ) are defined as [43] The relaxation lengths which couple to heat flux and bulk pressure (l Πq , l qΠ ), the heat flux and shear tensor (l qπ , l πq ) are defined as follows: In the ultra-relativistic limit, β(= m/T ) → 0, where m is the mass of the particle.We also have [31], To check the consistency of the terms in T µν F0 involving electromagnetic fields, we use the thermodynamic relation, where, µ is the chemical potential introduced to constrain the conservation of net (baryon) number.By using the conservation equations for n µ 0 and s µ 0 in the ideal hydrodynamics, it is straight forward to show that the hydrodynamic equation u ν ∂ µ T µν 0 = 0 along with the Maxwell equation Eq. ( 3) implies, which is in accordance with the standard thermodynamic relation, From Eqs. ( 44) and ( 46), we obtain the Gibbs-Duhem relation, The complete set of non-dissipative hydrodynamic equations in presence of an external magnetic field is formed by Eqs.( 39) and ( 40) along with the following two equations: Contracting Eq. ( 48) with b µ , one obtains the induction equation, where, θ ≡ ∂ µ u µ is the velocity divergence and D ≡ u µ ∂ µ is the co-moving derivative.All of these facts confirm that the form of EMT, T µν F0 is consistent with thermodynamic formulas for matter in the presence of electromagnetic fields [14].

III. HYDRODYNAMIC EQUATIONS IN LINEARISED FORM TO ESTIMATE THE DYNAMIC STRUCTURE FACTOR
Presently, we aim to evaluate the dynamic structure factor [S nn (k, ω)] by taking the correlation of dynamical density fluctuations in ω − k space.A small deviation from the equilibrium state of a thermodynamic variable can be accompanied by linearized form of the hydrodynamic equations.The density fluctuation can be obtained from the linearised equation to estimate the structure factor.Let the equilibrium state and the away from the equilibrium state of a thermodynamic quantity (n, ϵ, u α , q α , s, Π, π αβ etc.) are generically denoted by Q 0 and Q respectively.Therefore, any state (away from equilibrium) can be expressed as Q = Q 0 + δQ, where δQ is some tiny perturbation to the equilibrium value, Q 0 .
We can thus express the hydrodynamic equations around the equilibrium in the linearized form as: We can further simplify the set of Eqs.(51a)-(51f) by decomposing the fluid four velocity along the directions parallel and perpendicular to the direction of wave vector, k, and are termed as longitudinal component (δv || ) and transverse component (δv ⊥ ) respectively.The longitudinal and the transverse components of the fluid four velocity completely decouple the linearized hydrodynamic equations and we can get two sets linearly independent solutions for the two components.But the density perturbations appears only in the longitudinal component, which is our main focus here.Therefore, we consider the longitudinal component only for the present analysis.The hydrodynamic equations can be solved for a given set of initial condition, n(0), v || (0), T (0), q(0), Π(0) and π(0), by using the Fourier-Laplace transformation as: The δP and δs can be written in terms of the independent variables n and T as follows by using the thermodynamic relations: We use Eqs.( 52) and ( 53) to write down the longitudinal linearized hydrodynamic equation as: where, and here, A = (B 2 − M B).
The solution for density fluctuation is obtained by solving the set of above algebraic equations which gives, The appearance of A in the density fluctuation indicates the presence of magnetic field in the system.The expression of density fluctuation is a function of the other fluctuating hydrodynamic variables such as δT, δv || , δq || , δΠ, and δπ || .The independence among hydrodynamic variables does not necessarily imply the absence of correlations [44].However, in this work, we have specifically assumed that there are no correlations between density fluctuations and other fluctuating hydrodynamic variables [32,37].Therefore, the correlation between two independent thermodynamic variables, say, Q i and Q j vanishes i.e., The required correlator, S ′ nn (k, ω) is obtained as: Finally, the S nn (k, ω) is defined as: The S nn (k, ω) contains the transport coefficients such as η, ζ, κ and other thermodynamic response functions, which will be used to study the behaviour of the structure factor.It is well known that the relativistic Navier-Stokes (NS) hydrodynamic equation can be obtained by setting the various coupling (α 0 , α 1 ) and relaxation (β 0 , β 1 , β 2 ) coefficients to zero.
Therefore, the S nn (k, ω) for NS hydrodynamics can be obtained under the similar limits as: The appearance of higher order derivatives in MIS hydrodynamics makes the dispersion relation a quintic equation in ω and the corresponding dispersion equation in NS hydrodynamics is a cubic equation in ω.
The possibility of the existence of the critical end point (CEP) in the QCD phase diagram [45,46] is considered as one of the most interesting development in the field of relativistic heavy ion research.The location of the CEP in T − µ plane is not known from first principle.
The model dependent predictions vary widely as indicated in Ref. [47].In the present work we chose (T c , µ c ) = (154, 367)MeV to study the effects of the CEP on the spectral function in the presence of external magnetic field, B. The effects of CEP are included in the calculation of spectral function through the equation of state [48,49] and scaling behaviour of the transport coefficients [50][51][52].The details of this part of the calculations are given in Refs. [37,53].Therefore, we refer to these references for details to avoid repetition.

IV. RESULTS
Before we present the results it is important to mention the following points.In the context nuclear collisions at relativistic energies one needs to solve the relativistic viscous second order hydrodynamical equations in the presence of the CEP and ultra-high magnetic field.Development of such a numerical code is highly time consuming which will be considered in future publications.In the present work we focus on the study of the dynamical density correlation in a non-expanding QGP in the presence of an ultra-high magnetic field with and without the inclusion of the CEP.This study is important to check whether the magnetic field alters the nature of the correlation near the CEP.We find that the nature of the correlation near the CEP marked by the absence of the Brillouin peaks remains unaltered in presence of the magnetic field.This is one of the important observation of the present work as shown below.
It may also be mentioned here that the matter produced after nuclear collisions at LHC and top RHIC energies will be gluon dominated initially, however, quarks-antiquarks pairs will be produced dynamically due to the interactions among the gluons.Therefore, if the magnetic field (B) does not decay (which depends on the electrical conductivity of the QGP) completely before the generation of quark-antiquark pairs then the effects of B will be realized.Moreover, at lower RHIC energies (Beam Energy Scan Program) the quarks and antiquarks will be present from the very early stage and will get affected by the magnetic field.
We aim to study the nature of the dynamic structure factor [S nn (k, ω)] in the presence of a static magnetic field, B here.The S nn (k, ω) in the absence of B has been studied in earlier works [37,53], where it has been shown that the S nn (k, ω) admits three identifiable and the integrated intensities of those peaks are associated with various thermodynamic quantities such as the isothermal compressibilities and specific heats of the system [13].
Fig. 1 shows the S nn (k, ω) for eB = m 2 π at r = (T − T c )/T c = 0.2, i.e. when the system is away from the CEP.The transport coefficients are taken as η/s = ζ/s = χ T /s = 1/4π.We see three different peaks which are recognized as the R-peak (central), and the other two peaks as the B-peaks.The B-peaks are positioned symmetrically about ω = 0, but their heights are not identical.The unequal height of the B-peaks may occur due to the local inhomogeneity present in the system.The B-peaks arise from propagating sound modes associated with pressure fluctuations at constant entropy.In condensed matter physics, the asymmetry of the B-peaks is identified from the fact that two sound modes with different ω values, −c s k and +c s k originate from different temperature zones [54,55].Furthermore, we see that the widths of the peaks are very narrow, implying the slow relaxation of the thermal as well as the pressure fluctuations, which keeps the system to linger at out-of-equilibrium state for a longer time.
The introduction of B results in splitting the pressure into transverse, P ⊥ = P − M B and longitudinal, P L = P components.The expression of the magnetization is adapted from Ref. [14].The value of A determines the effects of the magnetic field on the spectral peak is found, and they start appearing after a threshold value of B th .Fig. 5 shows the S nn (k, ω) in the presence of the magnetic field (eB = 3 m 2 π ) in NS theory, which is recovered from the MIS model by putting the relaxation and coupling coefficients to zero.We see only three peaks, where the extra peaks due to the presence of the magnetic field are not appearing, but the presence of the magnetic field is evident from the magnitude as well as positional asymmetry of the B-peaks.Therefore, we argue that the extra near-side B-peaks are generated due to the coupling of hydrodynamic fields with the magnetic field when we consider second-order hydrodynamics.This is understood because the dispersion relations of MIS and NS hydrodynamics are quintic and cubic equations respectively.
The behaviour of the structure factor near the CEP (r = 0.01) is shown in Fig. 6.The transport coefficients and the thermodynamics response functions that appear in Eq.( 60) are taken from the scaling laws [37,[50][51][52][53].As the near side B-peaks appear only after a threshold intensity of the magnetic field, B th , we have plotted the S nn (k, ω) for two different values of magnetic field: a) B = m 2 π , and b) B = 3 m 2 π .In both cases, we do not observe any B-peaks but only the R-peak with a larger magnitude.This happens due to the absorption of sound near the CEP and hence B-peaks vanish.Comparing both the cases, we also observe that the magnitude of the R-peak in Fig. 6(b) is larger compared to the R-peak in Fig. 6(a).The most important observation here is that the B-peaks disappear near the CEP irrespective of the value of the external fields (zero or non-zero).

V. SUMMARY AND CONCLUSION
The dynamic structure factor, S nn (k, ω), is evaluated theoretically from the correlation of the density fluctuation, when the system with a CEP is subjected to an external static magnetic field.Without any magnetic field, the S nn (k, ω) admits a Rayleigh peak and two Brillouin peaks.But in presence of the magnetic field, we find two extra peaks appearing closer to the R-peak which are identified also as the B-peaks due to the adiabatic transverse pressure fluctuation.The asymmetry in magnitudes of the B-peaks (away side) are realized due to the local inhomogeneity of the system, whereas, the positional asymmetry appears due to the presence of the magnetic field.The extra B-peaks at lower ω exclusively appear beyond a threshold value of the magnetic field.The magnetic field splits the pressure into transverse, P ⊥ = P − M B and longitudinal, P L = P components.The away side and  the near side B-peaks are caused by the pressure fluctuations in longitudinal and transverse direction respectively.To understand the role of the relaxation and the coupling coefficients, we have evaluated S nn (k, ω) for NS theory and find two away side B-peaks, positioned asymmetrically with uneven heights.Therefore, we argue that the near side B-peaks appear due to the coupling of the magnetic field with the hydrodynamic fields in the second-order hydrodynamic theory.It will be interesting to investigate the role of this asymmetry in a realistic scenario with time varying electromagnetic field in an expanding QGP with CEP, and D µ ≡ H µν u ν .In the local rest frame of the fluid, the non-trivial components of these tensors are(F 10 , F 20 , F 30 ) = E, (F 32 , F 13 , F 21 ) = B, (M 10 , M 20 , M 30 ) = −P, (M 32 , M 13 , M 21 ) = M, (H 10 , H 20 , H 30 ) = D,and (H 32 , H 13 , H 21 ) = H, where P and M are the electric polarization and magnetization vector respectively.In linear domain, they are related to the fields E and B by the expressions P = χ e E and M = χ m B, where χ e and χ m represents the electric and magnetic susceptibilities respectively.Here E µ , B µ are space-like, i.e.E µ E µ = −E 2 and B µ B µ = −B 2 and orthogonal to u µ i.e.E µ u µ = 0, B µ u µ = 0 with E ≡ |E| and B ≡ |B|.

FIG. 3 .
FIG. 3. (Color online) The variation of the S nn (k, ω) with ω is shown in the presence of a magnetic field (B = 3 m 2 π ) for k = 0.1 fm −1 at r = 0.2.It shows one R-peak, two away side B-peaks asymmetrically located about ω = 0, and two near side B-peaks, originated from the coupling of the magnetic field with the other thermodynamic fields.

FIG. 5 .
FIG. 5. (Color online) The variation of the S nn (k, ω) with ω is shown in the presence of a magnetic field in Navier-Stokes theory.It shows one R-peak and two away side B-peaks asymmetrically located about ω = 0.

FIG. 6 .
FIG. 6. (Color online) The variation of the S nn (k, ω) with ω near the critical end point for k = 0.1 fm −1 at r = 0.01 is shown in the presence of a magnetic field in MIS theory for a) magnetic field strength B = m 2 π < B th , and b) magnetic filed strength B = 3 m 2 π ≥ B th .