Search

Rapid recent development in laser technology and methods learned from relativistic heavy ion physics led to new possibilities for fusion. Using a Hydrogen rich UDMA-TEGDMA polymer fusion target, laser irradiation ionizes the target. If we implant nanoantennas into the target resonating to the laser light frequency massive number of electrons of the ionized plasma resonate within the nanoantenna forming a so called nanoplasmonic wave. Our kinetic model simulation with a Hydrogen target indicates that the field of these resonating electrons attracts and accelerates the surrounding protons of the plasma to multi-MeV energy. These protons are then energetic enough to achieve nuclear transmutation and fusion reactions. Without resonating nanoantenna there is no such collective proton acceleration, no energetic protons, and nuclear reactions at 30 mJ laser pulse energy., Comment: arXiv admin note: text overlap with arXiv:2212.03638

The effects of a finite system volume on thermodynamic quantities, such as the pressure, energy density, specific heat, speed of sound, conserved charge susceptibilities and correlations, in hot and dense strongly interacting matter are studied within the parity-doublet Chiral Mean Field (CMF) model. Such an investigation is motivated by relativistic heavy-ion collisions, which create a blob of hot QCD matter of a finite volume, consisting of strongly interacting hadrons and potentially deconfined quarks and gluons. The effect of the finite volume of the system is incorporated by introducing a lower momentum cut-offs in the momentum integrals appearing in the model, the numerical value of the momentum cut-off being related to the de Broglie wavelength of the given particle species. It is found that some of these quantities show a significant volume dependence, in particular those sensitive to pion degrees of freedom, and the crossover transition is generally observed to become smoother in finite volume. These findings are relevant for the effective equation of state used in fluid dynamical simulations of heavy-ion collisions and efforts to extract the freeze out properties of heavy-ion collisions with susceptibilities involving electric charge and strangeness., Comment: 10 pages, 4 figures

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS., Comment: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Using relativistic supernova simulations of massive progenitor stars with a quark-hadron equation of state (EoS) and a purely hadronic EoS, we identify a distinctive feature in the gravitational-wave signal that originates from a buoyancy-driven mode (g-mode) below the proto-neutron star convection zone. The mode frequency lies in the range $200\lesssim f\lesssim 800\,\text{Hz}$ and decreases with time. As the mode lives in the core of the proto-neutron star, its frequency and power are highly sensitive to the EoS, in particular the sound speed around twice saturation density., Comment: Accepted Phys. Rev. Lett. (29 September 2023)

Recently laser induced fusion with simultaneous volume ignition, a spin-off from relativistic heavy ion collisions, was proposed, where implanted nanoantennas regulated and amplified the light absorption in the fusion target. Studies of resilience of the nanoantennas was published recently in vacuum. These studies are extended to nanoantennas embedded into a polymer, which modifies the nanoantenna's lifetime and absorption properties., Comment: 4 pages, 4 figures

A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $m_T$. We also show that an increase in $v_2$ is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $n_{\mathrm{CE}}<4 n_0$, show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $4n_0$., Comment: 12 pages, 5 figures, version accepted for publication

The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\texttt{-}100\,\text{M}_\odot$. We find explosions only for two low-compactness models ($14 \text{M}_\odot$ and $16\,\text{M}_\odot$) with the DD2F\_SF EoS, both with low explosion energies of $\mathord{\sim}10^{50}\,\mathrm{erg}$. These weak explosions are characterised by a neutrino signal with several mini-bursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as $^{90}\mathrm{Zr}$ for the $14\,\text{M}_\odot$ zero-metallicity model and $^{94}\mathrm{Zr}$ for the $16\,\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F\_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the proto-compact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to $4k_\text{B}/\text{baryon}$ and may have a distinctive gravitational wave signature, also for a second-order phase transition.

Binary neutron-star mergers and heavy-ion collisions are related through the properties of the hot and dense nuclear matter formed during these extreme events. In particular, low-energy heavy-ion collisions offer exciting prospects to recreate such {extreme} conditions in the laboratory. However, it remains unexplored to what degree those collisions can actually reproduce hot and dense matter formed in binary neutron star mergers. As a way to understand similarities and differences between these systems, we {discuss their geometry and }perform a direct numerical comparison of the thermodynamic conditions probed in both collisions. To enable a direct comparison, we employ a finite-temperature equation of state able to describe the entire high-energy phase diagram of Quantum Chromodynamics. Putting side by side the evolution of both systems, we find that laboratory heavy-ion collisions at the energy range of $E_{\mathrm{lab}}=0.4 - 0.6\ A$ MeV probe (thermodynamic) states of matter that are very similar to those created in binary neutron-star mergers. These results can inform future low-energy heavy-ion collisions probing this regime., Comment: 12 pages, 4 figures, accepted version

It is shown that the initial compression in central heavy ion collisions at beam energies of $E_\mathrm{lab}=1-10A$~GeV depends dominantly on the underlying equation of state and only marginally on the model used for the dynamical description. To do so, a procedure to incorporate any equation of state in the UrQMD transport model is introduced. In particular we compare the baryon density, temperature and pressure evolution as well as produced entropy in a relativistic ideal hydrodynamics approach and the UrQMD transport model, where the same equation of state is used in both approaches. Not only is the compression similar if the same equation of state is used in either dynamical model, but it also strongly depends on the actual equation of state. These results indicate that the equation of state can be studied with observables which are sensitive to the initial compression phase and maximum compression achieved in heavy ion collisions at these beam energies., Comment: 11 pages, 7 figures, version accepted by: The European Physical Journal C

The recent discovery of binary neutron star mergers has opened a new and exciting venue of research into hot and dense strongly interacting matter. For the first time this elusive state of matter, described by the theory of quantum chromo dynamics, can be studied in two very different environments. On the macroscopic scale in the collisions of neutron stars and on the microscopic scale in collisions of heavy ions at particle collider facilities. We will discuss the conditions that are created in these mergers and the corresponding high energy nuclear collisions. This includes the properties of QCD matter, i.e. the expected equation of state as well as expected chemical and thermodynamic properties of this exotic matter. To explore this matter in the laboratory - a new research prospect is available at the Facility for Antiproton and Ion Research, FAIR. The new facility is being constructed adjacent to the existing accelerator complex of the GSI Helmholtz Center for Heavy Ion Research at Darmstadt/Germany, expanding the research goals and technical possibilities substantially. The worldwide unique accelerator and experimental facilities of FAIR will open the way for a broad spectrum of unprecedented research supplying a variety of experiments in hadron, nuclear, atomic and plasma physics as well as biomedical and material science which will be briefly described., Comment: 10 pages, 8 figures, accepted contribution to a 'special issue' of the Astronomische Nachrichten

The thermal fit to preliminary HADES data of Au+Au collisions at $\sqrt{s_{_{NN}}}=2.4$ GeV shows two degenerate solutions at $T\approx50$ MeV and $T\approx70$ MeV. The analysis of the same particle yields in a transport simulation of the UrQMD model yields the same features, i.e. two distinct temperatures for the chemical freeze-out. While both solutions yield the same number of hadrons after resonance decays, the feeddown contribution is very different for both cases. This highlights that two systems with different chemical composition can yield the same multiplicities after resonance decays. The nature of these two minima is further investigated by studying the time-dependent particle yields and extracted thermodynamic properties of the UrQMD model. It is confirmed, that the evolution of the high temperature solution resembles cooling and expansion of a hot and dense fireball. The low temperature solution displays an unphysical evolution: heating and compression of matter with a decrease of entropy. These results imply that the thermal model analysis of systems produced in low energy nuclear collisions is ambiguous but can be interpreted by taking also the time evolution and resonance contributions into account., Comment: 8 pages, 6 figures

Second-order susceptibilities $\chi^{11}_{ij}$ of baryon, electric, and strangeness, $B$, $Q$, and $S$, charges, are calculated in the Chiral Mean Field (CMF) model and compared to available lattice QCD data. The susceptibilities are sensitive to the short range repulsive interactions between different hadron species, especially to the hardcore repulsion of hyperons. Decreasing the hyperons size, as compared to the size of the non-strange baryons, does improve significantly the agreement of the CMF model results with the Lattice QCD data. The electric charge-dependent susceptibilities are sensitive to the short range repulsive volume of mesons. The comparison with lattice QCD data suggests that strange baryons, non-strange mesons and strange mesons have significantly smaller excluded volumes than non-strange baryons. The CMF model with these modified hadron volumes allows for a mainly hadronic description of the QCD susceptibilities significantly above the chiral pseudo-critical temperature. This improved CMF model which is based on the lattice QCD data, has been used to study the properties of both cold QCD matter and neutron star matter. The phase structure in both cases is essentially unchanged, i.e. a chiral first-order phase transition occurs at low temperatures ($T_{\rm CP}\approx 17$ MeV), and hyperons survive deconfinement to higher densities than non-strange hadrons. The neutron star maximal mass remains close to 2.1$M_\odot$ and the mass-radius diagram is only modified slightly due to the appearance of hyperons and is in agreement with astrophysical observations., Comment: 11 pages, 9 figures

In this proceeding, we review our recent work using deep convolutional neural network (CNN) to identify the nature of the QCD transition in a hybrid modeling of heavy-ion collisions. Within this hybrid model, a viscous hydrodynamic model is coupled with a hadronic cascade "after-burner". As a binary classification setup, we employ two different types of equations of state (EoS) of the hot medium in the hydrodynamic evolution. The resulting final-state pion spectra in the transverse momentum and azimuthal angle plane are fed to the neural network as the input data in order to distinguish different EoS. To probe the effects of the fluctuations in the event-by-event spectra, we explore different scenarios for the input data and make a comparison in a systematic way. We observe a clear hierarchy in the predictive power when the network is fed with the event-by-event, cascade-coarse-grained and event-fine-averaged spectra. The carefully-trained neural network can extract high-level features from pion spectra to identify the nature of the QCD transition in a realistic simulation scenario., Comment: 4 pages, 4 figures, Quark Matter 2019, will be published in Nucl. Phys. A

Recently NAano-Plasmonic, Laser Inertial Fusion Experiments (NAPLIFE) were proposed, as an improved way to achieve laser driven fusion. The improvement is the combination of two basic research discoveries: (i) The possibility of detonations on space-time hyper-surfaces with time-like normal (i.e. simultaneous detonation in a whole volume) and (ii) to increase this volume to the whole target, by regulating the laser light absorption using nano-shells or nano-rods as antennas. These principles can be realized in an in-line, one dimensional configuration, in the simplest way with two opposing laser beams as in particle colliders. Such, opposing laser beam experiments were also performed recently. Here we study the consequences of the Laser Wake Field Acceleration (LWFA) if we experience it in a colliding laser beam set up. These studies can be applied to laser driven fusion, but also to other rapid phase transition, combustion, or ignition studies in other materials., Comment: 7 pages, 5 figures. arXiv admin note: text overlap with arXiv:2008.09847

The thermodynamic properties of high temperature and high density QCD-matter are studied using the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The CMF model provides a proper description of lattice QCD data, heavy-ions physics, and static neutron stars. The behavior of lines of constant pressure with increase of baryon density is discussed. The rapid change of pressure behavior at $\mu_B/T\approx3$ suggests a strong contribution of baryons to thermodynamic properties at this region. The position of this region is very close to the radius of convergence for a Taylor expansion of the QCD pressure. The role of mesons and unstable hadrons in the hydrodynamic expansion of strongly interacting matter is also discussed., Comment: 4 pages, 2 figures, Proceedings of the XXVIIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2019), Wuhan, China, November 3-9 2019

Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade "after-burner". Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments., Comment: 16 pages, 21 figures. Updated version for the publication in Eur. Phys. J. C

A method to determine the kinetic freeze-out temperature in heavy-ion collisions from measured yields of short-lived resonances is presented. The resonance production is treated in the framework of thermal model with an evolution between chemical and kinetic freeze-outs. The yields of many short-lived resonances are suppressed at $T = T_{\rm kin} < T_{\rm ch}$. We determine the values of $T_{\rm kin}$ and $T_{\rm ch}$ for various centralities in Pb--Pb collisions at $\sqrt{s_{_{NN}}} = 2.76$ TeV by fitting the abundances of both the stable hadrons and the short-lived resonances such as $\rho^0$ and $ \text{K}^{*0}$, that were measured by the ALICE collaboration. This allows to extract the kinetic freeze-out temperature from the measured hadron and resonance yields alone, independent of assumptions about the flow velocity profile and the freeze-out hypersurface. The extracted $T_{\rm ch}$ values exhibit a moderate multiplicity dependence whereas $T_{\rm kin}$ drops, from $T_{\rm kin} \simeq T_{\rm ch} \simeq 155$ MeV in peripheral collisions to $T_{\rm kin} \simeq 110$ MeV in 0-20% central collisions. Predictions for other short-lived resonances are presented. A potential (non-)observation of a suppressed $f_0(980)$ meson yield will allow to constrain the lifetime of that meson., Comment: 6 pages, 3 figures. Version published in Physical Review C

The production of backward nucleons, $N(180^\circ)$, at $180^\circ$ in the nuclear target rest frame in proton-nucleus ($\mathrm{p}+A$) collisions is studied. The backward nucleons appearing outside of the kinematically allowed range of proton-nucleon ($\mathrm{p}+N$) reactions are shown to be due to secondary reactions of heavy baryonic resonances produced inside the nucleus. Baryonic resonances $R$ created in primary $\mathrm{p}+N$ reactions can change their masses and momenta due to successive collisions $R+N\rightarrow R +N $ with other nuclear nucleons. Two distinct mechanisms and kinematic restrictions are studied: the reaction $R+N\rightarrow N(180^\circ)+N$ and the resonance decay $R\rightarrow N(180^\circ)+\pi$. Simulations of $\mathrm{p}+A$ collisions using the Ultra-relativistic Quantum Molecular Dynamics model support these mechanisms and are consistent with available data on proton backward production., Comment: 16 pages, 15 figures

The gravitational waves emitted from a binary neutron star merger, as predicted from general relativistic magneto-hydrodynamics calculations, are sensitive to the appearance of quark matter and the stiffness of the equation of state of QCD matter present in the inner cores of the stars. This is a new messenger observable from outer space, which does provide direct signals for the phase structure of strongly interacting QCD matter at high baryon density and high temperature. These astrophysically created extremes of thermodynamics do match, to within 20\%, the values of densities and temperatures which we find in relativistic hydrodynamics and transport theory of heavy ion collisions at the existing laboratories, if though at quite different rapidity windows, impact parameters and bombarding energies of the heavy nuclear systems. We demonstrate how one unified equation of state can be constructed and used for both neutron star physics and hot QCD matter excited at laboratory facilities. The similarity in underlying QCD physics allows the gravitational wave signals from future advanced LIGO and Virgo events to be combined with the analysis of high multiplicity fluctuations and flow measurements in heavy ion detectors in the lab to pin down the EoS and the phase structure of dense matter., Comment: 16 pages, 10 figures; contribution to the proceedings of the CPOD2018 conference, Corfu Island, Greece, September 24-28, 2018

The thermodynamic properties of high temperature and high density QCD-matter are explored within the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The quark sector of the CMF model is tuned to describe the $\mu_B=0$ thermodynamics data of lattice QCD. The resulting lines of constant physical variables as well as the baryon number susceptibilities are studied in some detail in the temperature/chemical potential plane. The CMF model predicts three consecutive transitions, the nuclear first-order liquid-vapor phase transition, chiral symmetry restoration, and the cross-over transition to a quark-dominated phase. All three phenomena are cross-over, for most of the $T-\mu_B$-plane. The deviations from the free ideal hadron gas baseline at $\mu_B=0$ and $T\approx 100-200$ MeV can be attributed to remnants of the liquid-vapor first order phase transition in nuclear matter. The chiral crossing transition determines the baryon fluctuations at much higher $\mu_B\approx1.5$ GeV, and at even higher baryon densities $\mu_B\approx2.4$ GeV, the behavior of fluctuations is controlled by the deconfinement cross-over. The CMF model also describe well the static properties of high $\mu_B$ neutron stars as well as the new neutron star merger observations. The effective EoS presented here describes simultaneously lattice QCD results at $\mu_B=0$, as well as observed physical phenomena (nuclear matter and neutron star matter) at $T\cong0$ and high densities, $\mu_B>1$ GeV., Comment: 13 pages, 11 figures

Thermodynamic properties of systems with repulsive interactions, are considered in the grand canonical ensemble. The analytic structure of the excluded-volume model in the complex plane of the system chemical potential (fugacity) is elaborated, based on the fact that the pressure function can be given in terms of the Lambert W-function. Even though the excluded volume model has no phase transitions at real values of the chemical potential, it does exhibit a branch cut singularity in the complex plane, thus limiting the convergence range of the Taylor expansion in the chemical potential. Close similarities to analytic properties of the other models with repulsive interactions, such as a cluster expansion model, the mean-field model, and the ideal Fermi gas model, are pointed out. As an example, repulsive baryonic interactions in a hadron gas, with a focus on the fugacity/virial and Taylor expansion methods used in lattice QCD, are presented. The asymptotic behavior of the Fourier expansion coefficients in these various models suggests that the singular part of net baryonic density can to leading order be universally expressed in terms of polylogarithms., Comment: 11 pages, 4 figures

A unified chiral mean field approach is presented for QCD thermodynamics in a wide range of temperatures and densities. The model simultaneously gives a satisfactory description of lattice QCD thermodynamics and fulfills nuclear matter and astrophysical constraints. The resulting equation of state can be incorporated in relativistic fluid-dynamical simulations of heavy-ion collisions and neutron stars mergers. Access to different regions of the QCD phase diagram can be obtained in simulations of heavy-ion data and observations of neutron star mergers., Comment: 4 pages, 2 figures. Contribution to the Quark Matter 2018 conference proceedings

Event-by-event multiplicity fluctuations in nucleus-nucleus collisions are studied within the relativistic transport models: EPOS, PHSD, and UrQMD. As measures of particle number fluctuations we consider the scaled variances $\omega[X]$ for positive, negative, and all charged hadrons, and the strongly intensive quantities $\Delta[K^+,\pi^+],$ $\Sigma[K^+,\pi^+]$ for $K^+$ and $\pi^+$ yields. At the SPS energy range the fluctuation measures are calculated for proton-proton, Ar+Sc, and Pb+Pb collisions. Comparison with recent NA61/SHINE and older NA49 measurements of the multiplicity fluctuations is done. A validity of the model of independent sources, a role of the experimental acceptance, and the centrality selection procedure are studied.Event-by-event multiplicity fluctuations in nucleus-nucleus collisions are studied within the relativistic transport models: EPOS, PHSD, and UrQMD. As measures of particle number fluctuations we consider the scaled variances $\omega[X]$ for positive, negative, and all charged hadrons, and the strongly intensive quantities $\Delta[K^+,\pi^+],$ $\Sigma[K^+,\pi^+]$ for $K^+$ and $\pi^+$ yields. At the SPS energy range the fluctuation measures are calculated for proton-proton, Ar+Sc, and Pb+Pb collisions. Comparison with recent NA61/SHINE and older NA49 measurements of the multiplicity fluctuations is done. A validity of the model of independent sources, a role of the experimental acceptance, and the centrality selection procedure are studied., Comment: 9 pages, 10 figures, Accepted for publication in Journal of Physics G

The quantum mechanical Beth-Uhlenbeck (BU) approach for repulsive hard-core interactions between baryons is applied to the thermodynamics of a hadron gas. The second virial coefficient $a_2$ -- the "excluded volume" parameter -- calculated within the BU approach is found to be temperature dependent, and it differs dramatically from the classical excluded volume (EV) model result. At temperatures $T =100-200$ MeV, the widely used classical EV model underestimates the EV parameter for nucleons at a given value of the nucleon hard-core radius by large factors of 3-4. Previous studies, which employed the hard-core radii of hadrons as an input into the classical EV model, have to be re-evaluated using the appropriately rescaled EV parameters. The BU approach is used to model the repulsive baryonic interactions in the hadron resonance gas (HRG) model. Lattice data for the second and fourth order net baryon susceptibilities are described fairly well when the temperature dependent BU baryonic excluded volume parameter corresponds to nucleon hard-core radii of $r_c = 0.25-0.3$ fm. Role of the attractive baryonic interactions is also considered. It is argued that HRG model with a constant baryon-baryon EV parameter $v_{NN} \simeq 1$ fm$^3$ provides a simple yet efficient description of baryon-baryon interaction in the crossover temperature region., Comment: 22 pages, 4 figures, accepted for publication in Phys. Rev. C

A generalization of the quantum van der Waals equation of state for a multi-component system in the grand canonical ensemble is proposed. The model includes quantum statistical effects and allows to specify the parameters characterizing repulsive and attractive forces for each pair of particle species. The model can be straightforwardly applied to the description of asymmetric nuclear matter and also for mixtures of interacting nucleons and nuclei. Applications of the model to the equation of state of an interacting hadron resonance gas are discussed., Comment: 20 pages, 5 figures, presentation improved, misprints corrected, version accepted for publication in Physical Review C