Source author record

Igor Mishustin

Igor Mishustin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2020arXiv

Skyrme-Hartree-Fock calculations of nuclear properties in the drip-point region of neutron star crust

In the present paper we explore the neutron-drip region of cold non-rotating isolated neutron stars. We have performed extended nuclear-structure calculations for nuclei embedded in the electron gas. For modeling the outer crust we use a set of Wigner-Seitz cells, where every cell contains one nucleus surrounded by a cloud of relativistic electrons. Above the drip point a non-relativistic neutron gas occurs in the cell. These calculations are carried out within the Hartree-Fock approach in combination with Skyrme effective interactions. For every baryon density we have determined the configuration with a minimal total energy. The drip elements and corresponding drip densities have been determined for about 240 different parametrizations of Skyrme forces used in the literature. We demonstrate that the calculated drip-point densities depend essentially on the Skyrme parametrization used. Even the drip elements and the occupied shells in the drip region differ for different parametrizations. We have found that the number of neutrons building the neutron gas at the drip point also depends essentially on the Skyrme force chosen. Nevertheless, the number density of the neutron gas in the drip-region is more or less the same ($\sim 10^{-5} \frac{neutrons}{{fm}^{3}}$). The drip densities obtained within our approach are generally lower than predicted earlier.

preprint2015arXiv

An effective model for the QCD phase transitions at finite baryon density

We introduce an effective quark-meson-nucleon model for the QCD phase transitions at finite baryon density. The nucleon and the quark degrees of freedom are described within a unified framework of a chiral linear sigma model. The deconfinement transition is modeled through a simple modification of the distribution functions of nucleons and quarks, where an additional auxiliary field, the bag field, is introduced. The bag field plays a key role in converting between the nucleon and the quark degrees of freedom. The model predicts that the chiral and the deconfinement phase transitions are always separated. Depending on the model parameters, the chiral transition occurs in the baryon density range of $(1.5-15.5)n_0$, while the deconfinement transition occurs above $5 n_0$, where $n_0$ is the saturation density.

preprint2015arXiv

Glueballs amass at RHIC and LHC Colliders! - The early quarkless 1st order phase transition at $T=270$ MeV - from pure Yang-Mills glue plasma to GlueBall-Hagedorn states

The early stage of high multiplicity pp, pA and AA collider is represented by a nearly quarkless, hot, deconfined pure gluon plasma. According to pure Yang-Mills Lattice Gauge Theory, this hot pure glue matter undergoes, at a high temperature, $T_c = 270$ MeV, a first order phase transition into a confined Hagedorn-GlueBall fluid. These new scenario should be characterized by a suppression of high $p_T$ photons and dileptons, baryon suppression and enhanced strange meson production. We propose to observe this newly predicted class of events at LHC and RHIC.

preprint2015arXiv

Nonlinear oscillations of compact stars in the vicinity of the maximum mass configuration

We solve the dynamical GR equations for the spherically symmetric evolution of compact stars in the vicinity of the maximum mass, for which instability sets in according to linear perturbation theory. The calculations are done with the analytical Zeldovich-like EOS P=a(rho-rho_0) and with the TM1 parametrisation of the RMF model. The initial configurations for the dynamical calculations are represented by spherical stars with equilibrium density profile, which are perturbed by either (i) an artificially added inward velocity field proportional to the radial coordinate, or (ii) a rarefaction corresponding to a static and expanded star. These configurations are evolved using a one-dimensional GR hydro code for ideal and barotropic fluids. Depending on the initial conditions we obtain either stable oscillations or the collapse to a black hole. The minimal amplitude of the perturbation, needed to trigger gravitational collapse is evaluated. The approximate independence of this energy on the type of perturbation is pointed out. At the threshold we find type I critical behaviour for all stellar models considered and discuss the dependence of the time scaling exponent on the baryon mass and EOS.

preprint2014arXiv

Comparative study of RBE and cell survival fractions for $^{1}$H, $^{4}$He, $^{12}$C and $^{16}$O beams using Geant4 and Microdosimetric Kinetic model

Beams of $^{4}$He and $^{16}$O nuclei are considered for ion-beam cancer therapy as alternative options to protons and $^{12}$C nuclei. Spread-out Bragg peak (SOBP) distributions of physical dose and relative biological effectiveness for 10% survival are calculated by means of our Geant4-based Monte Carlo model for Heavy Ion Therapy (MCHIT) and the modified microdosimetric kinetic model. The depth distributions of cell survival fractions are calculated for $^{1}$H, $^{4}$He, $^{12}$C and $^{16}$O for tissues with normal (HSG cells), low and high radiosensitivity. In each case the cell survival fractions were compared separately for the target volume, behind and in front of it. In the case of normal radiosensitivity $^{4}$He and $^{12}$C better spare tissues in the entrance channel compared to protons and $^{16}$O. The cell survival fractions calculated, respectively, for the entrance channel and target volume are similar for $^{4}$He and $^{12}$C. When it is important to spare healthy tissues located after the distal edge of the SOBP plateau, $^{4}$He can be recommended due to reduced nuclear fragmentation of these projectiles. No definite advantages of $^{16}$O with respect to $^{12}$C were found, with the except of an enhanced impact of these heavier projectiles on radioresistant tumors.

preprint2014arXiv

EMMI Rapid Reaction Task Force Meeting on 'Quark Matter in Compact Star'

The recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

preprint2014arXiv

Formation of droplets with high baryon density at the QCD phase transition in expanding matter

We consider the (3+1) dimensional expansion and cooling of the chirally-restored and deconfined matter at finite net-baryon densities as expected in heavy-ion collisions at moderate energies. In our approach, we consider chiral fields and the Polyakov loop as dynamical variables coupled to a medium represented by a quark-antiquark fluid. The interaction between the fields and the fluid leads to dissipation and noise, which in turn affect the field fluctuations. We demonstrate how inhomogeneities in the net-baryon density may form during an evolution through the spinodal region of the first-order phase transition. For comparison, the dynamics of transition through the crossover and critical end point is also considered.

preprint2014arXiv

Monte Carlo modeling of spallation targets containing uranium and americium

Neutron production and transport in spallation targets made of uranium and americium are studied with a Geant4-based code MCADS (Monte Carlo model for Accelerator Driven Systems). A good agreement of MCADS results with experimental data on neutron- and proton-induced reactions on $^{241}$Am and $^{243}$Am nuclei allows to use this model for simulations with extended Am targets. It was demonstrated that MCADS model can be used for calculating the values of critical mass for $^{233,235}$U, $^{237}$Np, $^{239}$Pu and $^{241}$Am. Several geometry options and material compositions (U, U+Am, Am, Am$_2$O$_3$) are considered for spallation targets to be used in Accelerator Driven Systems. All considered options operate as deep subcritical targets having neutron multiplication factor of $k \sim 0.5$. It is found that more than 4 kg of Am can be burned in one spallation target during the first year of operation.

preprint2014arXiv

Phase transition in compact stars due to a violent shock

In this letter we study the dynamics of a first order phase transition from nucleonic to quark matter in neutron stars. Using standard equations of state for these two phases we find the density range where such a transition is possible. Then we study the transformation of the star assuming that the quark core is formed via a spherical shock wave. The thermodynamical conditions in the quark core are found from the conservation laws across the transition region. Their dependence on the density and velocity of the incoming nuclear matter are studied. It is found that the shock is especially violent in the beginning of the conversion process when the velocity of the infalling matter is especially high. As the shock propagates further from the center the front velocity first increases and reaches a maximum value when the incoming velocity is around $0.2$. Finally, the front velocity quickly goes to zero when incoming matter velocity approaches zero. We have shown that the density and pressure jumps are especially large in the begining of the transition process.

preprint2014arXiv

Synthesis of neutron-rich transuranic nuclei in fissile spallation targets

A possibility of synthesizing neutron-reach super-heavy elements in spallation targets of Accelerator Driven Systems (ADS) is considered. A dedicated software called Nuclide Composition Dynamics (NuCoD) was developed to model the evolution of isotope composition in the targets during a long-time irradiation by intense proton and deuteron beams. Simulation results show that transuranic elements up to Bk-249 can be produced in multiple neutron capture reactions in macroscopic quantities. However, the neutron flux achievable in a spallation target is still insufficient to overcome the so-called fermium gap. Further optimization of the target design, in particular, by including moderating material and covering it by a reflector will turn ADS into an alternative source of transuranic elements in addition to nuclear fission reactors.

preprint2013arXiv

Chiral fluid dynamics with explicit propagation of the Polyakov loop

We present a fully dynamical model to study nonequilibrium effects in both the chiral and the deconfinement phase transition. The sigma field and the Polyakov loop as the corresponding order parameters are propagated by Langevin equations of motion. The locally thermalized background is provided by a fluid of quarks and antiquarks. Allowing for an exchange of energy and momentum through dissipative and stochastic processes we ensure that the total energy of the coupled system remains conserved. We study its relaxational dynamics in different quench scenarios and are able to observe critical slowing down as well as the enhancement of long wavelength modes at the critical point. During the fluid dynamical expansion of a hot plasma fireball typical nonequilibrium effects like supercooling and domain formation occur when the system evolves through the first order phase transition.

preprint2013arXiv

Implementation of chromomagnetic gluons in Yang-Mills thermodynamics

Motivated by the recent high-precision lattice data on Yang-Mills equations of state, we propose an effective theory of SU(3) gluonic matter. The theory is constructed based on the center and scale symmetries and their dynamical breaking, so that the interplay between color-electric and color-magnetic gluons is included coherently. We suggest, that the magnetic gluon condensate changes its thermal behavior qualitatively above the critical temperature, as a consequence of the matching to dimensionally-reduced magnetic theories. We consider thermodynamics in the mean field approximation and discuss the properties and interpretation of the trace anomaly.

preprint2013arXiv

Microdosimetry of radiation field from therapeutic C-12 beam in water: a study with Geant4 toolkit

We model the responses of Tissue-Equivalent Proportional Counters (TEPC) to radiation fields of therapeutic C-12 beams in a water phantom and to quasi-monoenergetic neutrons in a PMMA phantom. Simulations are performed with the Monte Carlo model for Heavy Ion Therapy (MCHIT) based on the Geant4 toolkit. The shapes of the calculated lineal energy spectra agree well with measurements in both cases. The influence of fragmentation reactions on the TEPC response to a narrow pencil-like beam with its width smaller than the TEPC diameter is investigated by Monte Carlo modeling. It is found that total lineal energy spectra are not very sensitive to the choice of the nuclear fragmentation model used. The calculated frequency-mean lineal energy y_f differs from the data on the axis of a therapeutic beam by less than 10% and by 10-20% at other TEPC positions. The validation of MCHIT with neutron beams gives us confidence in estimating the contributions to lineal energy spectra due to secondary neutrons produced in water by C-12 nuclei. As found, the neutron contribution at 10 cm distance from the beam axis amounts to ~ 50% close the entrance to the phantom and decreases to ~ 25% at the depth of the Bragg peak and beyond it. The presented results can help in evaluating biological out-of-field doses in carbon-ion therapy.

preprint2013arXiv

Microdosimetry spectra and RBE of H-1, He-4, Li-7 and C-12 nuclei in water studied with Geant4

A Geant4-based Monte Carlo model for Heavy-Ion Therapy (MCHIT) is used to study radiation fields of H-1, He-4, Li-7 and C-12 beams with similar ranges (~160-180 mm) in water. Microdosimetry spectra are simulated for wall-less and walled Tissue Equivalent Proportional Counters (TEPCs) placed outside or inside a phantom, as in experiments performed, respectively, at NIRS, Japan and GSI, Germany. The impact of fragmentation reactions on microdosimetry spectra is investigated for He-4, Li-7 and C-12, and contributions from nuclear fragments of different charge are evaluated for various TEPC positions in the phantom. The microdosimetry spectra measured on the beam axis are well described by MCHIT, in particular, in the vicinity of the Bragg peak. However, the simulated spectra for the walled TEPC far from the beam axis are underestimated. Relative Biological Effectiveness (RBE) of the considered beams is estimated using a modified microdosimetric-kinetic model. Calculations show a similar rise of the RBE up to 2.2-2.9 close to the Bragg peak for helium, lithium and carbon beams compared to the modest values of 1-1.2 at the plateau region. Our results suggest that helium and lithium beams are also promising options for cancer therapy.

preprint2011arXiv

The impact of dissipation and noise on fluctuations in chiral fluid dynamics

We investigate the nonequilibrium evolution of the sigma field coupled to a fluid dynamic expansion of a hot fireball to model the chiral phase transition in heavy-ion collisions. The dissipative processes and fluctuations are allowed under the assumption that the total energy of the coupled system is conserved. We use the linear sigma model with constituent quarks to investigate the effects of the chiral phase transition on the equilibration and excitation of the sigma modes. The quark fluid acts as a heat bath in local thermal equilibrium and the sigma field evolves according to a semiclassical stochastic Langevin equation of motion. The effects of supercooling and reheating in a first order phase transition are observed when the sigma field relaxes to equilibrium with the quark fluid. Nonequilibrium fluctuations at the first order phase transition lead to an increase in the intensity of sigma fluctuations in comparison to a scenario with a critical point.

preprint2011arXiv

The phase structure of a chiral model with dilatons in hot and dense matter

We explore the phase structure of a chiral model of constituent quarks and gluons implementing scale symmetry breaking at finite temperature and chemical potential. In this model the chiral dynamics is intimately linked to the trace anomaly saturated by a dilaton field. The thermodynamics is governed by two condensates, thermal expectation values of sigma and dilaton fields, which are the order parameters responsible for the phase transitions associated with the chiral and scale symmetries. Within the mean field approximation, we find that increasing temperature a system experiences a chiral phase transition and then a first-order phase transition of partial scale symmetry restoration characterized by a melting gluon-condensate takes place at a higher temperature. There exists a region at finite chemical potential where the scale symmetry remains dynamically broken while the chiral symmetry is restored. We also give a brief discussion on the sigma-meson mass constrained from Lattice QCD.

preprint2010arXiv

How large is "large $N_c$" for Nuclear matter?

We argue that a so far neglected dimensionless scale, the number of neighbors in a closely packed system, is relevant for the convergence of the large $N_c$ expansion at high chemical potential. It is only when the number of colors is large w.r.t. this new scale ($\sim \order{10}$) that a convergent large $N_c$ limit is reached. This provides an explanation as to why the large $N_c$ expansion, qualitatively successful in in vacuum QCD, fails to describe high baryo-chemical potential systems, such as nuclear matter. It also means that phenomenological claims about high density matter based on large $N_c$ extrapolations should be treated with caution.

preprint2010arXiv

On possibility of synthesizing superheavy elements in nuclear explosions

The possibility to produce superheavy elements in the course of low-yield nuclear explosions is analyzed within a simple kinetic model which includes neutron capture, gamma-emission, fission and particle evaporation from excited nuclei. We have calculated average numbers of absorbed neutrons as well as mass distributions of U and Cm nuclei exposed to an impulsive neutron flux, as functions of its duration. It is demonstrated that detectable amounts of heavy nuclei absorbing from 20 to 60 neutrons may be produced in this process. According to an optimistic scenario, after multiple beta-decay such nuclei may reach the long-living elements of the predicted ``island of stability''.

preprint2010arXiv

The nuclear liquid-gas phase transition at large $N_c$ in the Van der Waals approximation

We examine the nuclear liquid-gas phase transition at large number of colors ($N_c$) within the framework of the Van Der Waals (VdW) model. We argue that the VdW equation is appropriate at describing inter-nucleon forces, and discuss how each parameter scales with $N_c$. We demonstrate that $N_c=3$ (our world) is not large with respect to the other dimensionless scale relevant to baryonic matter, the number of neighbors in a dense system $N_N$. Consequently, we show that the liquid-gas phase transition looks dramatically different at $N_c \to \infty$ with respect of our world: The critical point temperature becomes of the order of $\lqcd$ rather than below it. The critical point density becomes of the order of the baryonic density, rather than an order of magnitude below it. These are precisely the characteristics usually associated with the "Quarkyonic phase". We therefore conjecture that quarkyonic matter is simply the large $N_c$ limit of the nuclear liquid, and the interplay between $N_c$ and $N_N$ is the reason why the nuclear liquid in our world is so different from quarkyonic matter. We conclude by suggesting ways our conjecture can be tested in future lattice measurements.

preprint2010arXiv

Thermodynamics of dense hadronic matter in a parity doublet model

We study thermodynamics of nuclear matter in a two-flavored parity doublet model within the mean field approximation. Parameters of the model are chosen to reproduce correctly the properties of the nuclear ground state. The model predicts two phase transitions in nuclear matter, a liquid-gas phase transition at normal nuclear density and a chiral transition at higher density. At finite temperature the pion decay constant exhibits a considerable reduction at intermediate values of chemical potential, which is traced back to the presence of the liquid-gas transition, and approaches zero at higher chemical potential associated with the chiral symmetry restoration. A "transition" from meson-rich to baryon-rich matter is also discussed.

preprint2009arXiv

Bulk viscosity-driven freeze-out in heavy ion collisions

We give an review the HBT puzzle, and argue that its resolution requires the introduction of new physics close to the phase transition scale. We argue that a candidate for this new physics is bulk viscosity, recently postulated to peak, and even diverge, close to the phase transition temperature. We show that such a viscosity peak can force the system created in heavy ion collisions to become unstable, and filament into fragments whose size is weakly dependent on the global size of the system, thereby triggering freeze-out.

preprint2009arXiv

Nuclear fragmentation reactions in extended media studied with Geant4 toolkit

It is well known from numerous experiments that nuclear multifragmentation is a dominating mechanism for production of intermediate-mass fragments in nucleus-nucleus collisions at energies above 100 A MeV. In this paper we investigate the validity and performance of the Fermi break-up model and the statistical multifragmentation model implemented as parts of the Geant4 toolkit. We study the impact of violent nuclear disintegration reactions on the depth-dose profiles and yields of secondary fragments for beams of light and medium-weight nuclei propagating in extended media. Implications for ion-beam cancer therapy and shielding from cosmic radiation are discussed.

preprint2009arXiv

Universality of the Diffusion Wake from Stopped and Punch-Through Jets in Heavy-Ion Collisions

We solve (3+1)-dimensional ideal hydrodynamical equations with source terms that describe punch-through and fully stopped jets in order to compare their final away-side angular correlations in a static medium. For fully stopped jets, the backreaction of the medium is described by a simple Bethe-Bloch-like model which leads to an explosive burst of energy and momentum (Bragg peak) close to the end of the jet's evolution through the medium. Surprisingly enough, we find that the medium's response and the corresponding away-side angular correlations are largely insensitive to whether the jet punches through or stops inside the medium. This result is also independent of whether momentum deposition is longitudinal (as generally occurs in pQCD energy loss models) or transverse (as the Bethe-Bloch formula implies). The existence of the diffusion wake is therefore shown to be universal to all scenarios where momentum as well as energy is deposited into the medium, which can readily be understood in ideal hydrodynamics through vorticity conservation. The particle yield coming from the strong forward moving diffusion wake that is formed in the wake of both punch-through and stopped jets largely overwhelms their weak Mach cone signal after freeze-out.

preprint2007arXiv

Comparative study of depth-dose distributions for beams of light and heavy nuclei in tissue-like media

We study the energy deposition by light and heavy nuclei in tissue-like media as used for cancer therapy. The depth-dose distributions for protons, $^{3}$He, $^{12}$C, $^{20}$Ne, and $^{58}$Ni nuclei are calculated within a Monte Carlo model based on the GEANT4 toolkit. These distributions are compared with each other and with available experimental data. It is demonstrated that nuclear fragmentation reactions essentially reduce the peak-to-plateau ratio of the dose profiles for deeply penetrating energetic ions heavier than $^{3}$He. On the other hand, all projectiles up to $^{20}$Ne were found equally suitable for therapeutic use at low penetration depths.

preprint2007arXiv

MCHIT - Monte Carlo model for proton and heavy-ion therapy

We study the propagation of nucleons and nuclei in tissue-like media within a Monte Carlo Model for Heavy-ion Therapy (MCHIT) based on the GEANT4 toolkit (version 8.2). The model takes into account fragmentation of projectile nuclei and secondary interactions of produced nuclear fragments. Model predictions are validated with available experimental data obtained for water and PMMA phantoms irradiated by monoenergetic carbon-ion beams. The MCHIT model describes well (1) the depth-dose distributions in water and PMMA, (2) the doses measured for fragments of certain charge, (3) the distributions of positron emitting nuclear fragments produced by carbon-ion beams, and (4) the energy spectra of secondary neutrons measured at different angles to the beam direction. Radial dose profiles for primary nuclei and for different projectile fragments are calculated and discussed as possible input for evaluation of biological dose distributions. It is shown that at the periphery of the transverse dose profile close to the Bragg peak the dose from secondary nuclear fragments is comparable to the dose from primary nuclei.

preprint2007arXiv

PET monitoring of cancer therapy with He-3 and C-12 beams: a study with the GEANT4 toolkit

We study the spatial distributions of $β^+$-activity produced by therapeutic beams of $^3$He and $^{12}$C ions in various tissue-like materials. The calculations were performed within a Monte Carlo model for Heavy-Ion Therapy (MCHIT) based on the GEANT4 toolkit. The contributions from $^{10,11}$C, $^{13}$N, $^{14,15}$O, $^{17,18}$F and $^{30}$P positron-emitting nuclei were calculated and compared with experimental data obtained during and after irradiation. Positron emitting nuclei are created by $^{12}$C beam in fragmentation reactions of projectile and target nuclei. This leads to a $β^+$-activity profile characterised by a noticeable peak located close to the Bragg peak in the corresponding depth-dose distribution. On the contrary, as the most of positron-emitting nuclei are produced by $^3$He beam in target fragmentation reactions, the calculated total $β^+$-activity during or soon after the irradiation period is evenly distributed within the projectile range. However, we predict also the presence of $^{13}$N, $^{14}$O, $^{17,18}$F created in charge-transfer reactions by low-energy $^3$He ions close to the end of their range in several tissue-like media. The time evolution of $β^+$-activity profiles was investigated for both kinds of beams. Due to the production of $^{18}$F nuclide the $β^+$-activity profile measured 2 or 3 hours after irradiation with $^{3}$He ions will have a distinct peak correlated with the maximum of depth-dose distribution. We found certain advantages of low-energy $^{3}$He beams over low-energy proton beams for reliable PET monitoring during particle therapy of shallow located tumours. In this case the distal edge of $β^+$-activity distribution from $^{17}$F nuclei clearly marks the range of $^{3}$He in tissues.

preprint1999arXiv

Reconstruction of the Proton Source in Relativistic Heavy Ion Collisions

We describe a direct method to reconstruct the transverse proton source formed in a relativistic heavy ion collision, making use of experimentally measured proton and deuteron spectra and assuming that deuterons are formed via two-nucleon coalescence. We show that an ambiguity with respect to the source temperature still persists and we indicate a possible solution to the problem.