-
Statistical approach of nuclear multifragmentation with realistic nuclear equation of state
Authors:
S. Mallik
Abstract:
In this work, Canonical Thermodynamical model for nuclear multifragmentation has been updated with realistic nuclear equation of state. Mass distribution, intermediate mass fragment multiplicity as well as isospin sensitive observables have been investigated with semi-microscopic approach of determining nuclear binding and excitation energies. Production of neutron rich isotopes as well as isoscal…
▽ More
In this work, Canonical Thermodynamical model for nuclear multifragmentation has been updated with realistic nuclear equation of state. Mass distribution, intermediate mass fragment multiplicity as well as isospin sensitive observables have been investigated with semi-microscopic approach of determining nuclear binding and excitation energies. Production of neutron rich isotopes as well as isoscaling and isobaric yield ratio parameters have been significantly modified due to inclusion of this realistic nuclear equation of state.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions
Authors:
Hermann Wolter,
Maria Colonna,
Dan Cozma,
Pawel Danielewicz,
Che Ming Ko,
Rohit Kumar,
Akira Ono,
ManYee Betty Tsang,
Jun Xu,
Ying-Xun Zhang,
Elena Bratkovskaya,
Zhao-Qing Feng,
Theodoros Gaitanos,
Arnaud Le Fèvre,
Natsumi Ikeno,
Youngman Kim,
Swagata Mallik,
Paolo Napolitani,
Dmytro Oliinychenko,
Tatsuhiko Ogawa,
Massimo Papa,
Jun Su,
Rui Wang,
Yong-Jia Wang,
Janus Weil
, et al. (27 additional authors not shown)
Abstract:
Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with…
▽ More
Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.
△ Less
Submitted 4 May, 2022; v1 submitted 14 February, 2022;
originally announced February 2022.
-
Binding energy shifts from heavy-ion experiments in a nuclear statistical equilibrium model
Authors:
S. Mallik,
H. Pais,
F. Gulminelli
Abstract:
Chemical constants extracted from $^{124}$Xe+ $^{124}$Sn collisions at 32 AMeV are compared to the predictions of an extended Nuclear Statistical Equilibrium model including mean-field interactions and in-medium binding energy shifts for the light ($Z\leq 2$) clusters. The ion species and density dependence of the in-medium modification is directly extracted from the experimental data. We show tha…
▽ More
Chemical constants extracted from $^{124}$Xe+ $^{124}$Sn collisions at 32 AMeV are compared to the predictions of an extended Nuclear Statistical Equilibrium model including mean-field interactions and in-medium binding energy shifts for the light ($Z\leq 2$) clusters. The ion species and density dependence of the in-medium modification is directly extracted from the experimental data. We show that the shift increases with the mass of the cluster and the density of the medium, and we provide a simple linear fit for future use in astrophysical simulations in the framework of the CompOSE data base. The resulting mass fractions are computed in representative thermodynamic conditions relevant for supernova and neutron star mergers. A comparison to the results of a similar analysis of the same data performed in the framework of a relativistic mean-field model shows a good agreement at low density, but significant discrepancies close to the Mott dissolution of clusters in the dense medium.
△ Less
Submitted 7 September, 2021;
originally announced September 2021.
-
Comparison of Heavy-Ion Transport Simulations: Mean-field Dynamics in a Box
Authors:
Maria Colonna,
Ying-Xun Zhang,
Yong-Jia Wang,
Dan Cozma,
Pawel Danielewicz,
Che Ming Ko,
Akira Ono,
Manyee Betty Tsang,
Rui Wang,
Hermann Wolter,
Jun Xu,
Zhen Zhang,
Lie-Wen Chen,
Hui-Gan Cheng,
Hannah Elfner,
Zhao-Qing Feng,
Myungkuk Kim,
Youngman Kim,
Sangyong Jeon,
Chang-Hwan Lee,
Bao-An Li,
Qing-Feng Li,
Zhu-Xia Li,
Swagata Mallik,
Dmytro Oliinychenko
, et al. (4 additional authors not shown)
Abstract:
Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared…
▽ More
Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. The significance of these results for the description of real heavy-ion collisions is discussed.
△ Less
Submitted 23 June, 2021;
originally announced June 2021.
-
New signatures of phase transition from Statistical Models of Nuclear multifragmentation
Authors:
G. Chaudhuri,
S. Mallik,
P. Das,
S. Das Gupta
Abstract:
The study of liquid-gas phase transition in heavy ion collisions has generated a lot of interest amongst the nuclear physicists in the recent years. In heavy ion collisions, there is no direct way of measuring the state variables like entropy, pressure, energy and hence unambiguous characterization of phase transition becomes difficult. This work proposes new signatures of phase transition that ca…
▽ More
The study of liquid-gas phase transition in heavy ion collisions has generated a lot of interest amongst the nuclear physicists in the recent years. In heavy ion collisions, there is no direct way of measuring the state variables like entropy, pressure, energy and hence unambiguous characterization of phase transition becomes difficult. This work proposes new signatures of phase transition that can be extracted from the observables which are easily accessible in experiments. It is observed that the temperature dependence of the first order derivative of the order parameters in nuclear liquid gas phase transition exhibit similar behavior as that of the variation of specific heat at constant volume Cv which is an established signature of first order phase transition. This motivates us to propose these derivatives as confirmatory signals of liquid-gas phase transition. The measurement of these signals in easily feasible in most experiments as compared to the other signatures like specific heat, caloric curve or bimodality. Total multiplicity, size of largest cluster are some of the order parameters which have been studied. Statistical Models based on canonical ensemble and lattice gas model has been used for the study. This temperature where the peak appears is designated to be the transition temperature and the effect of certain parameters on this has also been examined. The multiplicity derivative signature proposed in this work has been further confirmed by other theoretical models as well as in experimental study.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
Statistical and dynamical bimodality in multifragmentation reactions
Authors:
S. Mallik,
G. Chaudhuri,
F. Gulminelli,
S. Das Gupta
Abstract:
The bimodal behavior of the order parameter is studied in the framework of Boltzmann-Uehling-Uhlenbeck (BUU) transport model. In order to do that, simplified yet accurate method of BUU model is used which allow calculation of fluctuations in systems much larger than what was considered feasible in a well-known and already existing model. It is observed that depending on the projectile energy and c…
▽ More
The bimodal behavior of the order parameter is studied in the framework of Boltzmann-Uehling-Uhlenbeck (BUU) transport model. In order to do that, simplified yet accurate method of BUU model is used which allow calculation of fluctuations in systems much larger than what was considered feasible in a well-known and already existing model. It is observed that depending on the projectile energy and centrality of the reaction, both entrance channel and exit channel effects can be at the origin of the experimentally observed bimodal behavior. Both dynamical and statistical bimodality mechanisms are associated in the theoretical model to different time scales of the reaction, and to different energy regimes.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
Constraining the density dependence of the symmetry energy: the isospin transport ratio revisited
Authors:
S. Mallik,
F. Gulminelli
Abstract:
The isospin diffusion of the quasi-projectile formed in the $^{64,58}Ni$ on $^{64,58}Ni$ reactions in the Fermi energy domain is investigated in the framework of the Boltzmann-Uehling-Uhlenbeck transport model. The well known isospin transport ratio observable is revisited, with the aim of insuring an optimal comparison between experimental data and theoretical calculations and reducing the presen…
▽ More
The isospin diffusion of the quasi-projectile formed in the $^{64,58}Ni$ on $^{64,58}Ni$ reactions in the Fermi energy domain is investigated in the framework of the Boltzmann-Uehling-Uhlenbeck transport model. The well known isospin transport ratio observable is revisited, with the aim of insuring an optimal comparison between experimental data and theoretical calculations and reducing the present uncertainties in the extraction of empirical equation of state parameters. We show that isospin transport ratios are sensitive to all the low order isovector parameters ($E_{sym}$, $L_{sym}$ and $K_{sym}$). We demonstrate that realistic models of the equation of state, covering the uncertainty that presently affects the theoretical description of neutron stars static observables, can be effectively discriminated by isospin diffusion experiments, provided the neutron to proton ratio of the projectile remnant is precisely measured as a function of centrality.
△ Less
Submitted 27 October, 2020;
originally announced October 2020.
-
Nuclear Multifragmentation: Basic Concepts
Authors:
G. Chaudhuri,
S. Mallik,
S. Das Gupta
Abstract:
We present a brief overview of nuclear multifragmentation reaction. Basic formalism of canonical thermodynamical model based on equilibrium statistical mechanics is described. This model is used to calculate basic observables of nuclear multifragmentation like mass distribution, fragment multiplicity, isotopic distribution and isoscaling. Extension of canonical thermodynamical model to a projectil…
▽ More
We present a brief overview of nuclear multifragmentation reaction. Basic formalism of canonical thermodynamical model based on equilibrium statistical mechanics is described. This model is used to calculate basic observables of nuclear multifragmentation like mass distribution, fragment multiplicity, isotopic distribution and isoscaling. Extension of canonical thermodynamical model to a projectile fragmentation model is outlined. Application of the projectile fragmentation model for calculating average number of intermediate mass fragments and the average size of largest cluster at different $Z_{bound}$, differential charge distribution and cross-section of neutron rich nuclei of different projectile fragmentation reactions at different energies are described. Application of nuclear multifragmentation reaction in basic research as well as in other domains is outlined.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Isospin dependent hybrid model for studying isoscaling in heavy ion collisions around the Fermi energy domain
Authors:
S. Mallik,
G. Chaudhuri
Abstract:
Investigation of observables from nuclear multifragmentation reactions depending on isospin led to the development of a hybrid model. The mass and charge distribution as well as isotopic distribution was studied using this model for $^{112}$Sn+$^{112}$Sn reaction as well as $^{124}$Sn+$^{124}$Sn reactions at different energies. The agreement of the results obtained from the model with those from e…
▽ More
Investigation of observables from nuclear multifragmentation reactions depending on isospin led to the development of a hybrid model. The mass and charge distribution as well as isotopic distribution was studied using this model for $^{112}$Sn+$^{112}$Sn reaction as well as $^{124}$Sn+$^{124}$Sn reactions at different energies. The agreement of the results obtained from the model with those from experimental data confirms the accuracy of the model. Isoscaling coefficients were extracted from these observables which can throw light on the symmetry energy coefficient. Another important facet of this model is that temperature of the studied reaction can be directly extracted using this model.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Effect of liquid drop model parameters on nuclear liquid gas phase transition
Authors:
G. Chaudhuri,
S. Mallik
Abstract:
The phenomenon of liquid-gas phase transition occurring in heavy ion collisions at intermediate energies is a subject of contemporary interest. In statistical models of fragmentation, the liquid drop model is generally used to calculate the ground state binding energies of the fragments. It is well known that the surface and symmetry energy of the hot fragments at the low density freeze out can be…
▽ More
The phenomenon of liquid-gas phase transition occurring in heavy ion collisions at intermediate energies is a subject of contemporary interest. In statistical models of fragmentation, the liquid drop model is generally used to calculate the ground state binding energies of the fragments. It is well known that the surface and symmetry energy of the hot fragments at the low density freeze out can be considerably modified. In addition to this, the level density parameter also has a wide variation. The effect of variation of these parameters is studied on fragmentation observables which are related to the nuclear liquid gas phase transition. The canonical thermodynamical model which has been very successful in describing the phenomenon of fragmentation is used for the study. The shift in transition temperature owing to the variation in liquid drop model parameters has been examined.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Statistical and Dynamical Model Studies of Nuclear Multifragmentation Reactions at Intermediate Energies
Authors:
S. Mallik
Abstract:
Nuclear multifragmentation is an important phenomenon, the study of which can throw light on reaction mechanism in heavy ion collisions at intermediate and high energies. Based on statistical and dynamical model studies, this thesis is concentrated mainly on, the following three aspects of nuclear multifragmentation reactions namely (i) production of exotic nuclei which are normally not available…
▽ More
Nuclear multifragmentation is an important phenomenon, the study of which can throw light on reaction mechanism in heavy ion collisions at intermediate and high energies. Based on statistical and dynamical model studies, this thesis is concentrated mainly on, the following three aspects of nuclear multifragmentation reactions namely (i) production of exotic nuclei which are normally not available in the laboratory (ii) nuclear symmetry energy from heavy ion collisions at intermediate energies and (iii) Nuclear liquid-gas phase transition. In addition to these equivalence of statistical ensembles under different conditions is also studied in the framework of multifragmentation.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Statistical treatment of nuclear clusters in the continuum
Authors:
S. Mallik,
F. Gulminelli
Abstract:
The evaluation of the sub-saturation nuclear equation of state at finite temperature requires a proper state counting of the internal partition sum of nuclei which are immersed in the background of their continuum states. This classical statistical problem is addressed within the self-consistent mean-field approximation, which naturally accounts for isospin and effective mass effects in the nuclea…
▽ More
The evaluation of the sub-saturation nuclear equation of state at finite temperature requires a proper state counting of the internal partition sum of nuclei which are immersed in the background of their continuum states. This classical statistical problem is addressed within the self-consistent mean-field approximation, which naturally accounts for isospin and effective mass effects in the nuclear density of states. The nuclear free energy is decomposed into bulk and surface terms, allowing a simple analytical prescription for the subtraction of gas states from the nuclear partition sum, that avoids double counting of unbound single particle states. We show that this correction leads to a sizeable effect in the composition of matter at high temperature and low proton fractions, such as it is formed in supernova collapse, early proto-neutron star evolution, as well as laboratory experiments. Specifically, the energy stored in the internal nuclear degrees of freedom is reduced, as well as the mass fraction of heavy clusters in the statistical equilibrium. The gas subtraction prescription is compared to different phenomenological methods proposed in the literature, based on a high energy truncation of the partition sum. We show that none of these methods satisfactorily reproduces the gas subtracted level density, if the temperature overcomes ~4 MeV.
△ Less
Submitted 14 January, 2021; v1 submitted 26 August, 2020;
originally announced August 2020.
-
Sensitivity of the evaporation residue observables to the symmetry energy
Authors:
S. Mallik,
G. Chaudhuri,
F. Gulminelli
Abstract:
The static properties of the heaviest residue and unbound particles produced in central $^{64,58}Ni$ on $^{64,58}Ni$ collisions at 50 MeV/nucleon are predicted within the BUU transport model, in order to explore the sensitivity of those observables to the density dependence of the symmetry energy. We include fluctuations in the collision integral and use a meta-modelling for the mean-field which a…
▽ More
The static properties of the heaviest residue and unbound particles produced in central $^{64,58}Ni$ on $^{64,58}Ni$ collisions at 50 MeV/nucleon are predicted within the BUU transport model, in order to explore the sensitivity of those observables to the density dependence of the symmetry energy. We include fluctuations in the collision integral and use a meta-modelling for the mean-field which allows an independent variation of the different empirical parameters of the equation of state. We find that the isospin ratio of pre-equilibrium particles is a good estimator of the stiffness of the symmetry energy, in agreement with previous works. In addition to that, whatever be the functional form of the equation of state, we show that a higher symmetry energy at subsaturation densities leads to an increased size and isotopic ratio for the heaviest residue. This is understood in terms of energy sharing between the pre-equilibrium particles and the (quasi)fused system. The combination of the two observables might be an interesting tool to constrain the different density dependence below and above saturation, which is linked to the relatively poorly known parameter $K_{sym}$.
△ Less
Submitted 1 September, 2020; v1 submitted 15 April, 2019;
originally announced April 2019.
-
Comparison of heavy-ion transport simulations: Collision integral with pions and $Δ$ resonances in a box
Authors:
Akira Ono,
Jun Xu,
Maria Colonna,
Pawel Danielewicz,
Che Ming Ko,
Manyee Betty Tsang,
Yong-Jia Wang,
Hermann Wolter,
Ying-Xun Zhang,
Lie-Wen Chen,
Dan Cozma,
Hannah Elfner,
Zhao-Qing Feng,
Natsumi Ikeno,
Bao-An Li,
Swagata Mallik,
Yasushi Nara,
Tatsuhiko Ogawa,
Akira Ohnishi,
Dmytro Oliinychenko,
Jun Su,
Taesoo Song,
Feng-Shou Zhang,
Zhen Zhang
Abstract:
We compare ten transport codes for a system confined in a box, aiming at improved handling of the production of $Δ$ resonances and pions, which is indispensable for constraining high-density symmetry energy from observables such as the $π^-/π^+$ yield ratio in heavy-ion collisions. The system in a box is initialized with nucleons at saturation density and at 60 MeV temperature. The reactions…
▽ More
We compare ten transport codes for a system confined in a box, aiming at improved handling of the production of $Δ$ resonances and pions, which is indispensable for constraining high-density symmetry energy from observables such as the $π^-/π^+$ yield ratio in heavy-ion collisions. The system in a box is initialized with nucleons at saturation density and at 60 MeV temperature. The reactions $NN\leftrightarrow NΔ$ and $Δ\leftrightarrow Nπ$ are implemented, but the Pauli blocking and the mean-field potential are deactivated in the present comparison. Results are compared to those from the two reference cases of a chemically equilibrated ideal gas mixture and of the rate equation. In the results of the numbers of $Δ$ and $π$, deviations from the reference values are observed in many codes, and they depend significantly on the size of the time step. These deviations are tied to different ways in ordering the sequence of collisions and decays, that take place in the same time step. Better agreements are seen in the reaction rates and the number ratios among the isospin species of $Δ$ and $π$. These are, however, affected by the correlations, which are absent in the Boltzmann equation, but are induced by the way particle scatterings are treated in transport calculations. The uncertainty in the transport-code predictions of the $π^-/π^+$ ratio for the system initialized at n/p = 1.5, after letting the existing $Δ$ resonances decay, is found to be within a few percent, which is sufficiently small so that it does not strongly impact constraining the high-density symmetry energy from heavy-ion collisions. Most of the sources of uncertainties have been understood, and individual codes may be further improved. This investigation will be extended in the future to heavy-ion collisions to ensure the problems identified here remain under control.
△ Less
Submitted 27 September, 2019; v1 submitted 5 April, 2019;
originally announced April 2019.
-
Further studies of multiplicity derivative in models of heavy ion collision at intermediate energies as a probe for phase transition
Authors:
S Das Gupta,
S. Mallik,
G. Chaudhuri
Abstract:
In conjunction with models, the experimental observable total multiplicity can be used to check if the data contain the signature of phase transition and if it is first order. Two of the models reach similar conclusions. The third one is quite different.
In conjunction with models, the experimental observable total multiplicity can be used to check if the data contain the signature of phase transition and if it is first order. Two of the models reach similar conclusions. The third one is quite different.
△ Less
Submitted 16 April, 2018;
originally announced April 2018.
-
First-order derivative of cluster size as a new signature of phase transition in heavy ion collisions at intermediate energies
Authors:
P. Das,
S. Mallik,
G. Chaudhuri
Abstract:
The phenomenon of liquid-gas phase transition occurring in heavy ion collisions at intermediate energies is a subject of contemporary interest. Phase transition is usually characterized by the specific behaviour of state variables like pressure, density, energy etc. In heavy ion collisions there is no direct way of accessing these state variables and hence unambiguous detection of phase transition…
▽ More
The phenomenon of liquid-gas phase transition occurring in heavy ion collisions at intermediate energies is a subject of contemporary interest. Phase transition is usually characterized by the specific behaviour of state variables like pressure, density, energy etc. In heavy ion collisions there is no direct way of accessing these state variables and hence unambiguous detection of phase transition becomes difficult. This work establishes that signatures of phase transition can be extracted from the observables which are easily accessible in experiments and these have similar behaviour as the state variables. The temperature dependence of the first order derivative of the order parameters related to the largest and second largest cluster size (produced in heavy ion collisions) exhibit similar behavior as that of the variation of specific heat at constant volume Cv which is an established signature of first order phase transition. This motivates us to propose these derivatives as confirmatory signals of liquid-gas phase transition. The measurement of these signals in easily feasible in most experiments as compared to the other signatures like specific heat, caloric curve or bimodality. This temperature where the peak appears is designated to be the transition temperature and the effect of certain parameters on this has also been examined.
△ Less
Submitted 24 July, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Dynamical and statistical bimodality in nuclear fragmentation
Authors:
S. Mallik,
G. Chaudhuri,
F. Gulminelli
Abstract:
The origin of bimodal behavior in the residue distribution experimentally measured in heavy ion reactions is reexamined using Boltzmann-Uehling-Uhlenbeck simulations. We suggest that, depending on the incident energy and impact parameter of the reaction, both entrance channel and exit channel effects can be at the origin of the observed behavior. Specifically, fluctuations in the reaction mechanis…
▽ More
The origin of bimodal behavior in the residue distribution experimentally measured in heavy ion reactions is reexamined using Boltzmann-Uehling-Uhlenbeck simulations. We suggest that, depending on the incident energy and impact parameter of the reaction, both entrance channel and exit channel effects can be at the origin of the observed behavior. Specifically, fluctuations in the reaction mechanism induced by fluctuations in the collision rate, as well as thermal bimodality directly linked to the nuclear liquid-gas phase transition are observed in our simulations. Both phenomenologies were previously proposed in the literature, but presented as incompatible and contradictory interpretations of the experimental measurements. These results indicate that heavy ion collisions at intermediate energies can be viewed as a powerful tool to study both bifurcations induced by out-of-equilibrium critical phenomena, as well as finite size precursors of thermal phase transitions.
△ Less
Submitted 14 February, 2018;
originally announced February 2018.
-
Comparison of heavy-ion transport simulations: Collision integral in a box
Authors:
Ying-Xun Zhang,
Yong-Jia Wang,
Maria Colonna,
Pawel Danielewicz,
Akira Ono,
Betty Tsang,
Hermann Wolter,
Jun Xu,
Lie-Wen Chen,
Dan Cozma,
Zhao-Qing Feng,
Subal Das Gupta,
Natsumi Ikeno,
Che-Ming Ko,
Bao-An Li,
Qing-Feng Li,
Zhu-Xia Li,
Swagata Mallik,
Yasushi Nara,
Tatsuhiko Ogawa,
Akira Ohnishi,
Dmytro Oliinychenko,
Massimo Papa,
Hannah Petersen,
Jun Su
, et al. (5 additional authors not shown)
Abstract:
Simulations by transport codes are indispensable to extract valuable physics information from heavy ion collisions. In order to understand the origins of discrepancies between different widely used transport codes, we compare 15 such codes under controlled conditions of a system confined to a box with periodic boundary, initialized with Fermi-Dirac distributions at saturation density and temperatu…
▽ More
Simulations by transport codes are indispensable to extract valuable physics information from heavy ion collisions. In order to understand the origins of discrepancies between different widely used transport codes, we compare 15 such codes under controlled conditions of a system confined to a box with periodic boundary, initialized with Fermi-Dirac distributions at saturation density and temperatures of either 0 or 5 MeV. In such calculations, one is able to check separately the different ingredients of a transport code. In this second publication of the code evaluation project, we only consider the two-body collision term, i.e. we perform cascade calculations. When the Pauli blocking is artificially suppressed, the collision rates are found to be consistent for most codes (to within $1\%$ or better) with analytical results, or completely controlled results of a basic cascade code after eliminating the correlations within the same pair of colliding particles. In calculations with active Pauli blocking, the blocking probability was found to deviate from the expected reference values. The reason is found in substantial phase-space fluctuations and smearing tied to numerical algorithms and model assumptions in the representation of phase space. This results in the reduction of the blocking probability in most transport codes, so that the simulated system gradually evolves away from the Fermi-Dirac towards a Boltzmann distribution. As a result of this investigation, we are able to make judgements about the most effective strategies in transport simulations for determining the collision probabilities and the Pauli blocking. Investigation in a similar vein of other ingredients in transport calculations, like the mean field propagation or the production of nucleon resonances and mesons, will be discussed in the future publications.
△ Less
Submitted 17 April, 2018; v1 submitted 16 November, 2017;
originally announced November 2017.
-
Multiplicity Derivative: A new signature of first order phase transition in intermediate energy heavy ion collision
Authors:
S. Mallik,
G. Chaudhuri,
P. Das,
S Das Gupta
Abstract:
Measurement of M, the total multiplicity, for central collision between comparable mass heavy ions can provide a signature for first order phase transition. The derivative of M with respect to E*/A where E* is the excitation energy in the centre of mass and A the total mass of the dissociating system is expected to go through maximum as a function of E*. Theoretical modelling shows that this is th…
▽ More
Measurement of M, the total multiplicity, for central collision between comparable mass heavy ions can provide a signature for first order phase transition. The derivative of M with respect to E*/A where E* is the excitation energy in the centre of mass and A the total mass of the dissociating system is expected to go through maximum as a function of E*. Theoretical modelling shows that this is the energy where the specific heat Cv maximizes which typically happens at first order phase transition. The measurement of total M is probably feasible in more than one laboratory.
△ Less
Submitted 23 June, 2017; v1 submitted 23 May, 2017;
originally announced May 2017.
-
Effect of hyperons on phase coexistence in strange matter
Authors:
P. Das,
S. Mallik,
G. Chaudhuri
Abstract:
The study of liquid gas phase transition in fragmentation of nuclei in heavy ion collisions has been extended to the strangeness sector using the statistical model for multifragmentation. Helmholtz's free energy, specific heat and few other thermodynamic observables have been analyzed in order to examine the occurence of phase transition in the strange matter. The bimodal behaviour of the largest…
▽ More
The study of liquid gas phase transition in fragmentation of nuclei in heavy ion collisions has been extended to the strangeness sector using the statistical model for multifragmentation. Helmholtz's free energy, specific heat and few other thermodynamic observables have been analyzed in order to examine the occurence of phase transition in the strange matter. The bimodal behaviour of the largest cluster formed in fragmentation also strongly indicates coexistence of both the phases. The presence of hyperons strengthens the signals and also shifts the transition temperature to lower values.
△ Less
Submitted 23 May, 2017; v1 submitted 4 April, 2016;
originally announced April 2016.
-
Liquid gas phase transition in hypernuclei
Authors:
S. Mallik,
G. Chaudhuri
Abstract:
The fragmentation of excited hypernuclear system formed in heavy ion collisions has been described by the canonical thermodynamical model extended to three component systems. The multiplicity distribution of the fragments has been analyzed in detail and it has been observed that the hyperons have the tendency to get attached to the heavier fragments. Another important observation is the phase coex…
▽ More
The fragmentation of excited hypernuclear system formed in heavy ion collisions has been described by the canonical thermodynamical model extended to three component systems. The multiplicity distribution of the fragments has been analyzed in detail and it has been observed that the hyperons have the tendency to get attached to the heavier fragments. Another important observation is the phase coexistence of the hyperons, a phenomenon which is linked to liquid gas phase transition in strange matter.
△ Less
Submitted 4 April, 2016;
originally announced April 2016.
-
Effect of particle fluctuation on isoscaling and isobaric yield ratio of nuclear multifragmentation
Authors:
S. Mallik,
G. Chaudhuri
Abstract:
Isoscaling and isobaric yield ratio parameters are compared from canonical and grand canonical ensembles when applied to multifragmentation of finite nuclei. Source dependence of isoscaling parameters and source and isospin dependence of isobaric yield ratio parameters are examined in the framework of the canonical and the grand canonical models. It is found that as the nucleus fragments more, res…
▽ More
Isoscaling and isobaric yield ratio parameters are compared from canonical and grand canonical ensembles when applied to multifragmentation of finite nuclei. Source dependence of isoscaling parameters and source and isospin dependence of isobaric yield ratio parameters are examined in the framework of the canonical and the grand canonical models. It is found that as the nucleus fragments more, results from both the ensembles converge and observables calculated from the canonical ensemble coincide more with those obtained from the formulae derived using the grand canonical ensemble.
△ Less
Submitted 1 April, 2016;
originally announced April 2016.
-
Finite size effects on the phase diagram of the thermodynamical cluster model
Authors:
S. Mallik,
F. Gulminelli,
G. Chaudhuri
Abstract:
The thermodynamical cluster model is known to present a first-order liquid-gas phase transition in the idealized case of an uncharged, infinitely extended medium. However, in most practical applications of this model, the system is finite and charged. In this paper we study how the phase diagram is modified by finite size and Coulomb effects.
We show that the thermodynamic anomalies which are as…
▽ More
The thermodynamical cluster model is known to present a first-order liquid-gas phase transition in the idealized case of an uncharged, infinitely extended medium. However, in most practical applications of this model, the system is finite and charged. In this paper we study how the phase diagram is modified by finite size and Coulomb effects.
We show that the thermodynamic anomalies which are associated to the finite system counterpart of first order phase transitions, are correctly reproduced by this effective model. However, approximations in the calculation of the grandcanonical partition sum prevent obtaining the exact mapping between statistical ensembles which should be associated to finite systems. The ensemble inequivalence associated to the transition persists in the presence of Coulomb, but the phase diagram is deeply modified with respect to the simple liquid-gas phase transition characteristic of the neutral system.
△ Less
Submitted 1 April, 2016;
originally announced April 2016.
-
Bimodality emerges from transport model calculations of heavy ion collisions at intermediate energy
Authors:
S. Mallik,
S. Das Gupta,
G. Chaudhuri
Abstract:
This work is a continuation of our effort [Phys. Rev. C 91, 034616 (2015)] to examine if signatures of a phase transition can be extracted from transport model calculations of heavy ion collisions at intermediate energy. A signature of first order phase transition is the appearance of a bimodal distribution in $P_m(k)$ in finite systems. Here P_m(k) is the probability that the maximum of the multi…
▽ More
This work is a continuation of our effort [Phys. Rev. C 91, 034616 (2015)] to examine if signatures of a phase transition can be extracted from transport model calculations of heavy ion collisions at intermediate energy. A signature of first order phase transition is the appearance of a bimodal distribution in $P_m(k)$ in finite systems. Here P_m(k) is the probability that the maximum of the multiplicity distribution occurs at mass number k. Using a well-known model for event generation (BUU plus fluctuation), we do two cases of central collision: mass 40 on mass 40 and mass 120 on mass 120. Bimodality is seen in both the cases. The results are quite similar to those obtained in statistical model calculations.
An intriguing feature is seen. We observe that at the energy where bimodality occurs. other phase transition like signatures appear. There are breaks in certain first order derivatives. We then examine if such breaks appear in standard Botzmann-Uehling-Uhlenbeck (BUU) calculations without fluctuations. They do. The implication is interesting. If first order phase transition occurs, it may be possible to recognise that from ordinary BUU calculation. Probably the reason this was not seen already is because this aspect was not investigated before.
△ Less
Submitted 2 February, 2016;
originally announced February 2016.
-
Event simulations in a transport model for intermediate energy heavy ion collisions: Applications to multiplicity distributions
Authors:
S. Mallik,
S. Das Gupta,
G. Chaudhuri
Abstract:
We perform transport model calculations for central collisions of mass 120 on mass 120 at laboratory beam energy in the range 20 MeV/nucleon to 200 MeV/nucleon. A simplified yet accurate method allows calculation of fluctuations in systems much larger than what was considered feasible in a well-known and already existing model. The calculations produce clusters. The distribution of clusters is rem…
▽ More
We perform transport model calculations for central collisions of mass 120 on mass 120 at laboratory beam energy in the range 20 MeV/nucleon to 200 MeV/nucleon. A simplified yet accurate method allows calculation of fluctuations in systems much larger than what was considered feasible in a well-known and already existing model. The calculations produce clusters. The distribution of clusters is remarkably similar to that obtained in equilibrium statistical model.
△ Less
Submitted 19 March, 2015;
originally announced March 2015.
-
A hybrid model for studying nuclear multifragmentation around Fermi energy domain: Case for central collision of Xe on Sn
Authors:
S. Mallik,
G. Chaudhuri,
S. Das Gupta
Abstract:
Experimental data for central collisions of $^{129}$Xe on $^{119}$Sn at beam energies of (a) 32 MeV/nucleon, (b) 39 MeV/nucleon, (c) 45 MeV/nucleon and (d) 50 MeV/nucleon are compared with results calculated using a hybrid model. We use a transport model (BUU) to obtain the excitation energy per nucleon in the center of mass of the multifragmenting system. The canonical thermodynamic model is then…
▽ More
Experimental data for central collisions of $^{129}$Xe on $^{119}$Sn at beam energies of (a) 32 MeV/nucleon, (b) 39 MeV/nucleon, (c) 45 MeV/nucleon and (d) 50 MeV/nucleon are compared with results calculated using a hybrid model. We use a transport model (BUU) to obtain the excitation energy per nucleon in the center of mass of the multifragmenting system. The canonical thermodynamic model is then used to determine the temperature which would lead to this excitation energy. With this temperature we use the canonical thermodynamic model to calculate various experimental data such as multiplicities of different composites, probability distribution of the largest cluster etc. Agreement with data establishes the validity of the model.
△ Less
Submitted 17 March, 2015;
originally announced March 2015.
-
Estimates for temperature in projectile like fragment in geometric and transport models
Authors:
S. Mallik,
S. Das Gupta,
G. Chaudhuri
Abstract:
Projectile like fragments emerging from heavy ion collision have an excitation energy which is often labeled by a temperature. This temperature was recently calculated using a geometric model. We expand the geometric model to include also dynamic effects using a transport model. The temperatures so deduced agree quite well with values of temperature needed to fit experimental data.
Projectile like fragments emerging from heavy ion collision have an excitation energy which is often labeled by a temperature. This temperature was recently calculated using a geometric model. We expand the geometric model to include also dynamic effects using a transport model. The temperatures so deduced agree quite well with values of temperature needed to fit experimental data.
△ Less
Submitted 16 September, 2013;
originally announced September 2013.
-
Transformation between statistical ensembles in the modelling of nuclear fragmentation
Authors:
G. Chaudhuri,
F. Gulminelli,
S. Mallik
Abstract:
We explore the conditions under which the particle number conservation constraint deforms the predictions of fragmentation observables as calculated in the grand canonical ensemble. We derive an analytical formula allowing to extract canonical results from a grand canonical calculation and vice versa. This formula shows that exact canonical results can be recovered for observables varying linearly…
▽ More
We explore the conditions under which the particle number conservation constraint deforms the predictions of fragmentation observables as calculated in the grand canonical ensemble. We derive an analytical formula allowing to extract canonical results from a grand canonical calculation and vice versa. This formula shows that exact canonical results can be recovered for observables varying linearly or quadratically with the number of particles, independent of the grand canonical particle number fluctuations. We explore the validity of such grand canonical extrapolation for different fragmentation observables in the framework of the analytical Grand Canonical or Canonical Thermodynamical Model [(G)CTM] of nuclear multifragmentation. It is found that corrections to the grand canonical expectations can be evaluated with high precision, provided the system does not experience a first order phase transition. In particular, because of the Coulomb quenching of the liquid-gas phase transition of nuclear matter, we find that mass conservation corrections to the grand canonical ensemble can be safely computed for typical observables of interest in experimental measurements of nuclear fragmentation, even if deviations exist for highly exclusive observables.
△ Less
Submitted 24 May, 2013;
originally announced May 2013.
-
Temperature of projectile like fragments in heavy ion collisions
Authors:
S. Das Gupta,
S. Mallik,
G. Chaudhuri
Abstract:
A model in which a projectile like fragment can be simply regarded as a remnant after removal of some part of the projectile leads to an excited fragment. This excitation energy can be calculated with a Hamiltonian that gives correct nuclear matter binding, compressibility and density distribution in finite nuclei. In heavy ion collisions the model produces a dependence of excitation energy on imp…
▽ More
A model in which a projectile like fragment can be simply regarded as a remnant after removal of some part of the projectile leads to an excited fragment. This excitation energy can be calculated with a Hamiltonian that gives correct nuclear matter binding, compressibility and density distribution in finite nuclei. In heavy ion collisions the model produces a dependence of excitation energy on impact parameter which appears to be correct but the magnitude of the excitation energy falls short. It is argued that dynamic effects left out in the model will increase this magnitude. The model can be directly extended to include dynamics but at the expense of increased computation. For many calculations for observables, a temperature is an easier tool to use rather than an excitation energy. Hence temperature dependences on impact parameter in heavy ion collisions are displayed.
△ Less
Submitted 27 September, 2013; v1 submitted 6 May, 2013;
originally announced May 2013.
-
Symmetry energy from nuclear multifragmentation
Authors:
Swagata Mallik,
Gargi Chaudhuri
Abstract:
The ratio of symmetry energy coefficient to temperature $C_{sym}/T$ is extracted from different prescriptions using the isotopic as well as the isobaric yield distributions obtained in different projectile fragmentation reactions. It is found that the values extracted from our theoretical calculation agree with those extracted from the experimental data but they differ very much from the input val…
▽ More
The ratio of symmetry energy coefficient to temperature $C_{sym}/T$ is extracted from different prescriptions using the isotopic as well as the isobaric yield distributions obtained in different projectile fragmentation reactions. It is found that the values extracted from our theoretical calculation agree with those extracted from the experimental data but they differ very much from the input value of the symmetry energy used. The best possible way to deduce the value of the symmetry energy coefficient is to use the fragment yield at the breakup stage of the reaction and it is better to use the grand canonical model for the fragmentation analysis. This is because the formulas that are used for the deduction of the symmetry energy coefficient are all derived in the framework of the grand canonical ensemble which is valid only at the break-up (equilibrium) condition. The yield of "cold" fragments either from the theoretical models or from experiments when used for extraction of the symmetry energy coefficient using these prescriptions might lead to the wrong conclusion.
△ Less
Submitted 23 January, 2013;
originally announced January 2013.
-
Evaluating transport coefficients in real time thermal field theory
Authors:
S. Mallik,
Sourav Sarkar
Abstract:
Transport coefficients in a hadronic gas have been calculated earlier in the imaginary time formulation of thermal field theory. The steps involved are to relate the defining retarded correlation function to the corresponding time-ordered one and to evaluate the latter in the conventional perturbation expansion. Here we carry out both the steps in the real time formulation.
Transport coefficients in a hadronic gas have been calculated earlier in the imaginary time formulation of thermal field theory. The steps involved are to relate the defining retarded correlation function to the corresponding time-ordered one and to evaluate the latter in the conventional perturbation expansion. Here we carry out both the steps in the real time formulation.
△ Less
Submitted 12 November, 2012;
originally announced November 2012.
-
Conditions for equivalence of Statistical Ensembles in Nuclear Multifragmentation
Authors:
Swagata Mallik,
Gargi Chaudhuri
Abstract:
Statistical models based on canonical and grand canonical ensembles are extensively used to study intermediate energy heavy ion collisions. The underlying physical assumption behind canonical and grand canonical models is fundamentally different, and in principle agree only in the thermodynamical limit when the number of particles become infinite. Nevertheless, we show that these models are equiva…
▽ More
Statistical models based on canonical and grand canonical ensembles are extensively used to study intermediate energy heavy ion collisions. The underlying physical assumption behind canonical and grand canonical models is fundamentally different, and in principle agree only in the thermodynamical limit when the number of particles become infinite. Nevertheless, we show that these models are equivalent in the sense that they predict similar results if certain conditions are met even for finite nuclei. In particular, the results converge when nuclear multifragmentation leads to the formation of predominantly nucleons and low mass clusters. The conditions under which the equivalence holds are amenable to present day experiments.
△ Less
Submitted 30 October, 2012; v1 submitted 11 September, 2012;
originally announced September 2012.
-
A model for projectile fragmentation
Authors:
G. Chaudhuri,
S. Mallik,
S. Das Gupta
Abstract:
A model for projectile fragmentation is developed whose origin can be traced back to the Bevalac era. The model positions itself between the phenomenological EPAX parametrization and transport models like "Heavy Ion Phase Space Exploration" (HIPSE) model and antisymmetrised molecular dynamics (AMD) model. A very simple impact parameter dependence of input temperature is incorporated in the model w…
▽ More
A model for projectile fragmentation is developed whose origin can be traced back to the Bevalac era. The model positions itself between the phenomenological EPAX parametrization and transport models like "Heavy Ion Phase Space Exploration" (HIPSE) model and antisymmetrised molecular dynamics (AMD) model. A very simple impact parameter dependence of input temperature is incorporated in the model which helps to analyze the more peripheral collisions. The model is applied to calculate the charge, isotopic distributions, average number of intermediate mass fragments and the average size of largest cluster at different Z_{bound} of different projectile fragmentation reactions at different energies.
△ Less
Submitted 28 March, 2013; v1 submitted 3 September, 2012;
originally announced September 2012.
-
Fluctuations in the number of intermediate mass fragments in small projectile like fragments
Authors:
S. Mallik,
G. Chaudhuri,
S. Das Gupta
Abstract:
The origin of fluctuations in the average number of intermediate mass fragments seen in experiments in small projectile like fragments is discussed. We argue that these can be explained on the basis of a recently proposed model of projectile fragmentation.
The origin of fluctuations in the average number of intermediate mass fragments seen in experiments in small projectile like fragments is discussed. We argue that these can be explained on the basis of a recently proposed model of projectile fragmentation.
△ Less
Submitted 2 November, 2011;
originally announced November 2011.
-
Improvements to model of projectile fragmentation
Authors:
S. Mallik,
G. Chaudhuri,
S. Das Gupta
Abstract:
In a recent paper [Phys. Rev. C 044612 (2011)] we proposed a model for calculating cross-sections of various reaction products which arise from disintegration of projectile like fragment resulting from heavy ion collisions at intermediate or higher energy. The model has three parts: (1) abrasion, (2) disintegration of the hot abraded projectile like fragment (PLF) into nucleons and primary composi…
▽ More
In a recent paper [Phys. Rev. C 044612 (2011)] we proposed a model for calculating cross-sections of various reaction products which arise from disintegration of projectile like fragment resulting from heavy ion collisions at intermediate or higher energy. The model has three parts: (1) abrasion, (2) disintegration of the hot abraded projectile like fragment (PLF) into nucleons and primary composites using a model of equilibrium statistical mechanics and (3) possible evaporation of hot primary composites. It was assumed that the PLF resulting from abrasion has one temperature T. Data suggested that while just one value of T seemed adequate for most cross-sections calculations, it failed when dealing with very peripheral collisions. We have now introduced a variable T=T(b) where b is the impact parameter of the collision. We argue there are data which not only show that T must be a function of b but, in addition, also point to an approximate value of T for a given b. We propose a very simple formula: T(b)=D_0+D_1(A_s(b)/A_0) where A_s(b) is the mass of the abraded PLF and A_0 is the mass of the projectile; D_0 and D_1 are constants. Using this model we compute cross-sections for several collisions and compare with data.
△ Less
Submitted 22 August, 2011;
originally announced August 2011.
-
Model for projectile fragmentation: case study for Ni on Ta, Be and Xe on Al
Authors:
S. Mallik,
G. Chaudhuri,
S. Das Gupta
Abstract:
For projectile fragmentation we work out details of a model whose origin can be traced back to the Bevalac era. The model positions itself between the phenomenological EPAX parametrization and microscopic transport models like "Heavy Ion Phase Space Exploration Model" (HIPSE) and antisymmetrised molecular dynamics(AMD). We apply the model to some recent data of projectile fragmentation of Ni on Ta…
▽ More
For projectile fragmentation we work out details of a model whose origin can be traced back to the Bevalac era. The model positions itself between the phenomenological EPAX parametrization and microscopic transport models like "Heavy Ion Phase Space Exploration Model" (HIPSE) and antisymmetrised molecular dynamics(AMD). We apply the model to some recent data of projectile fragmentation of Ni on Ta and Be at beam energy 140 MeV/nucleon and some older data of Xe on Al at beam energy 790 MeV/nucleon. Reasonable values of cross-sections for various composites populated in the reactions are obtained.
△ Less
Submitted 30 December, 2010;
originally announced January 2011.
-
Effect of secondary decay on isoscaling: Results from the canonical thermodynamical model
Authors:
Gargi Chaudhuri,
Swagata Mallik
Abstract:
The projectile fragmentation reactions using $^{58}Ni$ $\&$ $^{64}Ni$ beams at
140 MeV/n on targets$^{9}Be$ $\&$ $^{181}Ta$ are studied using the canonical thermodynamical model coupled with an evaporation code. The isoscaling property of the fragments produced is studied using both the primary and the secondary fragments and it is observed that the secondary fragments also respect isoscaling thou…
▽ More
The projectile fragmentation reactions using $^{58}Ni$ $\&$ $^{64}Ni$ beams at
140 MeV/n on targets$^{9}Be$ $\&$ $^{181}Ta$ are studied using the canonical thermodynamical model coupled with an evaporation code. The isoscaling property of the fragments produced is studied using both the primary and the secondary fragments and it is observed that the secondary fragments also respect isoscaling though the isoscaling parameters $α$ and $β$ changes. The temperature needed to reproduce experimental data with the secondary fragments is less than that needed with the primary ones. The canonical model coupled with the evaporation code successfully explains the experimental data for isoscaling for the projectile fragmentation reactions.
△ Less
Submitted 13 July, 2010;
originally announced July 2010.
-
Relativistic spectral function of nucleon in hot nuclear matter
Authors:
Sabyasachi Ghosh,
S. Mallik,
Sourav Sarkar
Abstract:
We present a simple calculation of the nucleon self-energy in nuclear matter at finite temperature in a relativistic framework, using the real time thermal field theory. The imaginary parts of one-loop graphs are identified with discontinuities across the unitary and the Landau cuts. We find that in general both the cuts contribute significantly to the spectral function in the region of (virtual)…
▽ More
We present a simple calculation of the nucleon self-energy in nuclear matter at finite temperature in a relativistic framework, using the real time thermal field theory. The imaginary parts of one-loop graphs are identified with discontinuities across the unitary and the Landau cuts. We find that in general both the cuts contribute significantly to the spectral function in the region of (virtual) nucleon mass usually considered, even though the unitary cut is ignored in the literature. Also our relativistic spectral function differs from the one in non-relativistic approximation, used in some earlier calculations.
△ Less
Submitted 28 September, 2010; v1 submitted 13 April, 2010;
originally announced April 2010.
-
Analytic structure of rho meson propagator at finite temperature
Authors:
Sabyasachi Ghosh,
S. Mallik,
Sourav Sarkar
Abstract:
We analyse the structure of one-loop self-energy graphs for the rho meson in real time formulation of finite temperature field theory. We find the discontinuities of these graphs across the unitary and the Landau cuts. These contributions are identified with different sources of medium modification discussed in the literature. We also calculate numerically the imaginary and the real parts of the s…
▽ More
We analyse the structure of one-loop self-energy graphs for the rho meson in real time formulation of finite temperature field theory. We find the discontinuities of these graphs across the unitary and the Landau cuts. These contributions are identified with different sources of medium modification discussed in the literature. We also calculate numerically the imaginary and the real parts of the self-energies and construct the spectral function of the rho meson, which are compared with an earlier determination. A significant contribution arises from the unitary cut of the pi-omega loop, that was ignored so far in the literature.
△ Less
Submitted 28 September, 2010; v1 submitted 18 November, 2009;
originally announced November 2009.
-
S-matrix approach to equation of state of nuclear matter
Authors:
S. Mallik,
J. N. De,
S. K. Samaddar,
Sourav Sarkar
Abstract:
We calculate the equation of state of nuclear matter based on the general analysis of the grand canonical partition function in the $S$-matrix framework. In addition to the low mass stable particles and their two-body scattering channels considered earlier, the calculation includes systematically all the higher mass particles and their exited states as well as the scattering channels formed by a…
▽ More
We calculate the equation of state of nuclear matter based on the general analysis of the grand canonical partition function in the $S$-matrix framework. In addition to the low mass stable particles and their two-body scattering channels considered earlier, the calculation includes systematically all the higher mass particles and their exited states as well as the scattering channels formed by any number of these species. We estimate the latter contribution by resonances in all the channels. The resulting model-independent virial series for pressure gets substantial contribution from the heavy particles and the channels containing them. The series converges for larger values of baryon density than found earlier.
△ Less
Submitted 18 March, 2008; v1 submitted 3 January, 2008;
originally announced January 2008.
-
QCD sum rules for nucleon two-point function in nuclear medium
Authors:
S. Mallik,
Hiranmaya Mishra
Abstract:
We derive QCD sum rules from the nucleon two-point function in nuclear medium, calculating its specral function in chiral perturbation theory to one loop. Our calculation shows the inadequacy of the commonly used ansatz to represent the two-point function as the sum of a nucleon pole term and a spectral integral over the high energy continuum. We also point out that it is the energy variable and…
▽ More
We derive QCD sum rules from the nucleon two-point function in nuclear medium, calculating its specral function in chiral perturbation theory to one loop. Our calculation shows the inadequacy of the commonly used ansatz to represent the two-point function as the sum of a nucleon pole term and a spectral integral over the high energy continuum. We also point out that it is the energy variable and not its square that must be used in writing the dispersion integrals. The sum rules are obtained by equating terms linear in nucleon number density of the spectral and operator representations of the correlation function. We use them to predict the curent-nucleon coupling, in good agreement with the value known from the vacuum sum rules. They also predict a definite deviation of the four-quark condensate from its ground state saturation value.
△ Less
Submitted 28 October, 2004;
originally announced October 2004.
-
On the nucleon self-energy in nuclear matter
Authors:
S. Mallik,
Andreas Nyffeler,
M. C. M. Rentmeester,
Sourav Sarkar
Abstract:
We consider the nucleon self-energy in nuclear matter in the absence of Pauli blocking. It is evaluated using the partial-wave analysis of $NN$ scattering data. Our results are compared with that of a realistic calculation to estimate the effect of this blocking. It is also possible to use our results as a check on the realistic calculations.
We consider the nucleon self-energy in nuclear matter in the absence of Pauli blocking. It is evaluated using the partial-wave analysis of $NN$ scattering data. Our results are compared with that of a realistic calculation to estimate the effect of this blocking. It is also possible to use our results as a check on the realistic calculations.
△ Less
Submitted 13 July, 2004;
originally announced July 2004.
-
Pion parameters in nuclear medium from chiral perturbation theory and virial expansion
Authors:
S. Mallik,
Sourav Sarkar
Abstract:
We consider two methods to find the effective parameters of the pion traversing a nuclear medium. One is the first order chiral perturbation theoretic evaluation of the pion pole contribution to the two-point function of the axial-vector current. The other is the exact, first order virial expansion of the pion self-energy. We find that, although the results of chiral perturbation theory are not…
▽ More
We consider two methods to find the effective parameters of the pion traversing a nuclear medium. One is the first order chiral perturbation theoretic evaluation of the pion pole contribution to the two-point function of the axial-vector current. The other is the exact, first order virial expansion of the pion self-energy. We find that, although the results of chiral perturbation theory are not valid at normal nuclear density, those from the virial expansion may be reliable at such density. The latter predicts both the mass-shift and the in-medium decay width of the pion to be small, of about a few MeV.
△ Less
Submitted 11 September, 2003;
originally announced September 2003.
-
QCD sum rules for rho mesons in nuclear matter
Authors:
S. Mallik,
Andreas Nyffeler
Abstract:
We investigate QCD sum rules for vector currents in the rho meson channel in the nuclear medium. For increased sensitivity, we subtract out the vacuum contributions. With a saturation scheme often considered in the literature, we find these subtracted sum rules to be unstable. It indicates the importance of other contributions neglected in the saturation scheme. These include more Landau singula…
▽ More
We investigate QCD sum rules for vector currents in the rho meson channel in the nuclear medium. For increased sensitivity, we subtract out the vacuum contributions. With a saturation scheme often considered in the literature, we find these subtracted sum rules to be unstable. It indicates the importance of other contributions neglected in the saturation scheme. These include more Landau singularities at low energy, additional operators in the medium, and possibly the in-medium width of the rho meson.
△ Less
Submitted 31 May, 2001; v1 submitted 6 February, 2001;
originally announced February 2001.