Finite Temperature Asymmetric Nuclear Matter Using Pion Dressing

We study the asymmetric nuclear matter in a nonperturvative manner at finite temperatures using thermofield dynamics method. The nucleon-meson interaction is taken to examine the binding energy (EB), pressure (P) for various proton fractions.


Introduction
The properties of hot dense nuclear matter are very important in the context of neutron stars [1] as well as in high-energy heavy-ion collisions experiments.In those experiments, nuclei undergo collisions [2] and produce states of nuclear matter at different temperatures at conditions far from those normally encountered in low energy collisions.At low temperatures, this can be attained in medium energy heavy-ion reactions; no such exotic state can be produced but there is the possibility of interesting liquid-gas phase transitions leading to the breakup of heated nuclei into small clusters or droplets of nucleons [3].
Such phase transitions have been identified in multifragmentation experiments and in the crusts of neutron stars [4] [5].The possibilities of such phase transitions were previously considered by several authors using different approaches [6]- [8].Walecka [9] in a attempt, developed a method consisting of interactions of nucleons with a neutral scalar field σ and ρ mesons as well as neutral vector meson ω.Variations of the same model have also been considered including cubic and quartic terms in the σ fields to reproduce correct bulk modulus of nuclear matter [10].These calculations however use meson fields as classical, and use a σ-field which is not observed.
A model for infinite nuclear matter consisting of interacting nucleons and pions was considered in [11].The Pion Nucleon interaction was taken with off mass shell of pion.The scalar isoscalar pion condensates simulated the effects of σ mesons [12] with the short distance repulsion arising from composite structure of nucleons and/or through vector meson exchanges.This appeals aesthetically as classical σ fields arise from quantum mechanical structures and also has a stronger phenomenological appeal as σ mesons have not been found in nature.With a similar approach we shall reconsider nuclear matter at finite temperatures with pion dressing [13] with off mass shell.The methods of thermofield dynamics [14] fit naturally for this purpose because here statistical average is done through an expectation value over a "thermal vacuum" [11] in an extended Hilbert space.
The article is organised as follows.In Section 2, we review the thermofield dynamics to consider hot and dense nuclear matter and obtain expressions for temperature-dependant pressure P, binding energy B and nuclear density ρ.In Section 3, we evaluate numerically the above applying variational technique, study their characteristics and discuss the results so obtained.

Formalism
We consider the effective Hamiltonian for pion nucleon interaction at zero temperature [15]- [17] as where ( ) ( ) are the Hamiltonians for the free nucleon part, the pion-interaction part and the free meson part respectively.Here ( ) and the effective Hamiltonian ( ) We have taken ( ) , with M as the mass of the nucleon.The free meson part of the Hamiltonian is given by ( ) ( ) ( ) where m π denotes the mass of the meson and The pion field expansion in Equation (4) in terms of annihilation and creation operators is given by with ( ) . We now study asymmetric nuclear matter at finite temperature using thermofield dynamics method.Here the thermal average of a quantum operator Ô given in [13] [14] with The above reduces to ground state expectation, in zero temperature limit value for the operator Ô , shown below as where 0 is the lowest energy state.In thermofield method, one generalises [13] to the case of finite temperature and for some ( ) where ( ) is defined as the "thermal vacuum".Here corresponding to physical operator a, a "tilde" operator a  is introduced to denote thermal modes.At finite temperature the ground state is replaced by ( ) where ( ) † a −k  in the above corresponds to the extra Hilbert space.The function ( ) β is calculated through minimization of thermodynamic potential density given as ( ) where µ , the chemical potential, N, number of nucleons and the entropy density For zero chemical potential and for free fields, extremization of the free energy then yields with Hamiltonian density as If we substitute this value the free energy density becomes Similarly for the fermionic sector the thermal vacuum is where † I ψ corresponds to the creation of the fermionic thermal modes.The entropy is given as where µ is the chemical potential corresponding to baryon number conservation and ( ) , k M ω β = + k for free fermions of mass M. In this methodology the correct distribution function with extra thermal modes can be obtained and this enables us to have the temperature dependent background off-shell pion pair configuration given as ( ) ( ) exp and In the above, In nucleon sector for fermions with where ( ) will be determined later.We then have the nuclear matter density Here the thermodynamic potential [18] is given by , where p µ and n µ are chemical potentials of proton and neutron respectively, n ρ and p ρ are neutron and proton densities and S is the entropy of the asymmetric nuclear matter ( ) θ for the interacting sys- tem will be determined here from the construction of the thermodynamic potential.We thus have for the nucleon kinetic term temperature dependent as ( ) The temperature dependent kinetic energy due to the mesons is given by where ( ) f k is pion dressing function [13] which extremises energy, pressure and entropy etc of nuclear matter.Then the interaction energy density becomes where we shall now assume a phenomenological term corresponding to meson repulsion due to composite structure of mesons given as Finally, the nucleon repulsion term given as where ρ is as given in Equation (20).Thus the energy density is given by ( ) where as before.
The thermodynamic potential density Ω is given by where the last term corresponds to nucleon number conservation with p µ and n µ as the chemical potentials of proton and neutron respectively.The temperature dependent pressure is given as P = −Ω .The entropy density above is with F S being the entropy in fermion sector given as and similarly the meson sector contribution B S is given as Thus the thermodynamic potential density now is a functional of ( ) as well as the pion dressing function ( ) f  k which will of course depend upon temperature.Extremisation of Equation ( 29) with respect to ( ) which is of the same form as [16] for zero temperature.Similarly minimising the thermodynamic potential with respect to ( ) where Once we substitute the optimised dressing as in Equation (30), the above simplifies to ( ) which is different from ω due to interactions.Further, minimising the thermodynamic potential with respect to ( ) we have the solution where J is given in Equation (25).We may note that the change in F  above from 2 2 k M is also due to interaction.
In Figure 1 we study the nature of variation of E B with different baryon densities at different temperatures.The saturation binding energy increases from −16 MeV at zero temperature to higher values with rise of temperature.This clearly shows that the temperature has a significant effect on symmetric nuclear matter.
In Figure 2, the variation of pressure P with nucleon densities have been studied at different temperatures (T = 0, 4, 8, 12, 16, 20 MeV).It shows that at the saturation density, the pressure is zero for zero temperature as usual.But pressure becomes negative if one goes towards lower density and it becomes zero at zero nuclear density.This peculear behaviour is first order phase transition.This trend disappears at higher temperatures.
In  It shows that the slope becomes negative as the relative nucleon densities increases to a certain value and then increases at higher values.The trend is same at different higher temperatures but terminates to a certain value of slope at certain lower densities.

Summary
In this paper, we study the warm equation of state (EOS) of asymmetric nuclear matter taking pion-nucleon interaction with repulsive effect due to ρ and ω mesons.We observe here that our results are comparable with the Non-Linear Walecka (NLWM) and Quark Meson Coupling (QMC) model [20].The binding energy changes from negative to positive value with increase of temperature around saturation density.The profile of pressure variation with density at different temperatures shows that there is a liquid-gas phase transition and the transition  We also observe that the pressure gradient becomes more negative due to increase of proton fraction leading to a possibility of instability.

λ
[16] are determined variationally.The first three of these four para- meters are evaluated by constraining the binding energy per nucleon P and the compressibility K of the symmetric nuclear matter i.e. of our calculation is evaluated by fixing the value of sym E to the standard value of 31 MeV.The pion-nucleon coupling constant 2 4π 14.

Figure 3 ,Figure 1 .
Figure 1.The binding energy E B as a function of nucleon density ρ.

Figure 2 .In Figure 4 ,
Figure 2. The pressure P as a function of nucleon density ρ at different temperature.

Figure 3 .
Figure 3.The pressure P as a function of nucleon density ρ for different proton fraction.

Figure 4 .
Figure 4.The pressure variation ( ) , p T y P ρ ∂ ∂ as a function of relative nucleon density