2025-05-13 | | Total: 19
Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. We further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.
The astrophysical direct nuclear capture reaction ^{12}{\rm C}(p, \gamma)^{13}{\rm N} is studied within the framework of a potential model. Parameters of the nuclear p-^{12}C interaction potentials of the Woods-Saxon form are adjusted to reproduce experimental p-^{12}C scattering phase shifts, as well as the binding energies and empirical values of the asymptotic normalization coefficient (ANC) for the ^{13}N(1/2^-) ground state from the literature. The reaction rates are found to be very sensitive to the description of the value of the ANC of the ^{13}N(1/2^{-}) ground state and width of the ^{13}N(1/2^+) resonance at the E_x=2.365 MeV excitation energy. The potential model, which yields the ANC value of 1.63 fm^{-1/2} for the ^{13}N(1/2^{-}) ground state and a value \Gamma=39 keV for the ^{13}N(1/2^+) resonance width, is able to reproduce the astrophysical S factor in the energy interval up to 2 MeV, the empirical values of the reaction rates in the temperature region up to T=10^{10} K of the LUNA Collaboration and the results of the R-matrix fit. The astrophysical factor S(0)=1.35 keV b is found using the asymptotic expansion method of D. Baye. The obtained value is in a good agreement with the Solar Fusion II result. At the same time, the calculated value of 1.44 keV b of the astrophysical S factor at the Solar Gamow energy is consistent with the result of the R-matrix fit of S(25~\rm{keV})=1.48 \pm 0.09 keV b by Kettner {\it et al.}, but slightly less than the result of S(25~\rm{keV})=1.53 \pm 0.06 keV b the LUNA Collaboration.
The earliest phase of an ultrarelativistic heavy ion collision can be described as a highly populated system of gluons called glasma. We study some glasma characteristics related to the system's angular momentum. The first one is the global angular momentum perpendicular to the reaction plane, which is spanned by the beam axis and impact parameter vector. We show that only a small fraction of the enormous initial angular momentum of the colliding ultrarelativistic nuclei is transferred to the glasma. Our main focus is the glasma angular momentum directed along the beam axis. This quantity has a local character and results from the inhomogeneous velocity field generated in the glasma. We calculate the vorticity and local angular momentum which show noticeably different behaviour. The results are analyzed in detail and discussed in the context of existing experimental data. We argue that neither vorticity nor thermal vorticity but instead the local angular momentum controls the polarization of final-state hadrons.
Using the Beyond-Skyrme-Hartree-Fock approach with various Skyrme-type NN effective interactions, we explore the large spectroscopic quadrupole moment Q_s of the first excited 2_{1}^{+} state in ^{20}Ne. Our calculated Q_s(2_{1}^{+}) \approx -0.20~e\textrm{b} aligns closely with the experimental value of -0.22(2)~e\textrm{b}. Consistent with other theoretical methods, our results confirm the presence of an \alpha cluster structure in this state. We also investigate the robustness of the \alpha cluster against tensor-force, demonstrating that incorporating tensor components significantly enhances Q_s(2_{1}^{+}) and strengthens the cluster structure along the symmetrical axis.
A microscopic investigation of the rotational properties in ^{129}\mathrm{Cs} and ^{131}\mathrm{La} was carried out using the three-dimensional tilted axis cranking covariant density functional theory (3DTAC-CDFT). The calculations reveal the coexistence of magnetic and chiral rotation built on identical qusiparticle configurations \pi {h}_{11/2}^{1}\otimes \nu {h}_{11/2}^{-2} in ^{129}\mathrm{Cs} and ^{131}\mathrm{La}, establishing a new type of shape coexistence. The calculations predict the deformation parameters for magnetic rotation in ^{129}\mathrm{Cs} (\beta\approx0.23, \gamma \approx41^\circ) and ^{131}\mathrm{La} (\beta \approx0.25, \gamma \approx42^\circ), along with those for possible chiral rotation in ^{129}\mathrm{Cs} (\beta \approx0.20, \gamma \approx29^\circ) and in ^{131}\mathrm{La} (\beta \approx0.20, \gamma \approx27^\circ). The energy spectra, the relation between the spin and the rotational frequency, and the reduced M1 and E2 transition probabilities are obtained with the various configurations. The experimental characteristics of band B8 in ^{129}\mathrm{Cs} and band 13 in ^{131}\mathrm{La} are well reproduced. Moreover, a distinctive rotational mode transition is uncovered in this work, progressing from the principal-axis rotation to the planar rotation and finally to the chiral rotation, as the rotational frequency evolves through different regimes.
We investigate the number-of-constituent-quark (NCQ) scaling of elliptic flow for various hadrons in non-central Au+Au collisions at \sqrt{s_{NN}} = 3.0\text{--}7.7\,\mathrm{GeV} using the AMPT model with string melting (SM) and pure hadron cascade (HC) modes. For the SM case, NCQ scaling is absent at \sqrt{s_{NN}} =3.0\,\mathrm{GeV} but is largely restored by \sqrt{s_{NN}} =3.9\,\mathrm{GeV}. Although quark coalescence occurs at \sqrt{s_{NN}} =3.0\,\mathrm{GeV}, the lack of NCQ scaling is attributed to the insufficient development of quark elliptic flow and the limited production of strange quarks and antiquarks. This finding suggests that NCQ scaling could not be considered a definitive signature of quark-gluon plasma (QGP) formation in the RHIC fixed-target energy region. For the HC case, as expected, no NCQ scaling is observed. However, a mass ordering in the elliptic flow emerges at \sqrt{s_{NN}} =4.5\,\mathrm{GeV}, indicating that full thermalization may not be a prerequisite for mass ordering.
We suggest a new look onto the thermodynamics of cold isospin QCD matter described within the quarkyonic framework, when confining forces are essential only for the momentum states close to the Fermi sphere and lead to formation of the quark-antiquark correlations with quantum numbers of pions. Within this picture cold isospin QCD matter is constituted by the Bose-Einstein condensate of pions and free (anti)quarks. We develop a simple confining model of pions existing as the lowest energy P-state of in-medium three-dimensional harmonic oscillator and use it to show that under the conditions of isospin QCD quarkyonic matter can onset at densities significantly smaller than the overlap density of pions. Using the experimentally established pion charge radius we demonstrate that the onset density can be even lower than two nuclear saturation densities. To analyze the lattice QCD data we develop a field-theory inspired model of quark-antiquark-pion matter with vector-isovector repulsion among pions and (anti)quarks. The quark substructure of pions, which is one of the key elements of the model, is accounted for within the one-loop no-sea approximation and leads to an effective medium dependence of the pion mass. The latter is shown to be crucial for providing asymptotic vanishing of the pion density and pion condensate, which reflects pion dissociation due to asymptotic freedom of QCD. The proposed quarkyonic picture of isospin QCD agrees with the results of perturbative calculations on reaching the conformal limit of QCD and with the data of lattice simulations, which supports its viability and physical relevance.
We calculate Gamow-Teller strength distributions for \beta \beta-decay nuclei ^{76}Ge and ^{82}Se using the deformed pn-QRPA model. We use a deformed Nilsson basis and consider pairing correlations within the deformed BCS theory. Ground state correlations and two-particle and two-hole mixing states were included in our pn-QRPA model. Our calculated strength distributions were compared with experimental data and previous calculation. The total Gamow-Teller strength and centroid placement calculated in our model compares well with the measured value. We calculate \beta-decay and positron capture rates on ^{76}Ge and ^{82}Se in supernovae environments and compare them to those obtained from experimental data and previous calculation. Our study shows that positron capture rates command the total weak rates at high stellar temperatures. We also calculate energy rates of \beta-delayed neutrons and their emission probabilities.
We theoretically investigate Lambda production via a (K-, pi-) reaction on a 6Li target, using the distorted-wave impulse approximation (DWIA) with a Fermi-averaged K-n --> pi-Lambda amplitude. We calculate Lambda production spectra using the Green's function method for a 5Li nuclear core + Lambda system, employing a Lambda folding-model potential based on the 5Li nuclear density, which is constructed from alpha-p and 3He-d cluster wave functions. The results show that the calculated spectrum, which includes Lambda bound, resonance, and continuum states, agrees well with the experimental data from the (K-, pi-) reaction at p_K- = 790 MeV/c (0 deg.), where substitutional (0p_Lambda, 0p^{-1}_n) and (0s_Lambda, 0s^{-1}_n) configurations dominate in the near-recoilless reactions. A narrow peak corresponds to a high-lying excited state with spin-parity J^P= 1+ at E_Lambda= 13.8 MeV near the 3He + d + Lambda threshold, arising from interference effects between 5Li(3/2+)x(0s1/2)_Lambda and 5Li(1/2+)x(0s1/2)_Lambda components. This study offers a valuable framework for extracting essential information on the structure and production mechanisms of hypernuclear states from experimental data.
We probe the superfluid-superconductor dynamics of the rotating outer core of neutron stars through half-vortex states. By means of a generalized hydrodynamic model, where proton and neutron fluids are coupled by both dynamic entrainment and Skyrme SLy4 nucleon-nucleon interactions, we analyze single flux tubes in the proton-superconductor component of the system that thread proton vortices located faraway from neutron vortices. It is shown how they give rise to hydrodynamic perturbations in the coexisting neutron superfluid, and its structure remains unaltered for varying rotation rates and magnetic fields within ranges of observational values.
We calculate temperature-dependent nuclear partition functions (TDNPFs) and nuclear abundances for 728 nuclei assuming nuclear statistical equilibrium (NSE). The theories of stellar evolution support NSE. Discrete nuclear energy levels have been calculated \textit{microscopically}, using the pn-QRPA theory, up to an excitation energy of 10 MeV in the calculation of TDNPFs. This feature of our paper distinguishes it from previous calculations. Experimental data is also incorporated wherever available to ensure reliability of our results. Beyond 10 MeV we employ simple Fermi gas model and perform integration over the nuclear level densities to approximate the TDNPFs. We calculate nuclidic abundances, using the Saha equation, as a function of three parameters: stellar density, stellar temperature and lepton-to-baryon content of stellar matter. All these physical parameters are considered to be extremely important in stellar interior. Results obtained in this paper show that the equilibrium configuration of nuclei remains unaltered by increasing stellar density (only calculated nuclear abundances increases by roughly same order of magnitude). Increasing the stellar temperature smooths the equilibrium configuration showing peaks at neutron-number magic nuclei.
The allowed Gamow-Teller (GT) transitions are the most common weak nuclear processes of spin-isospin (\sigma\tau) type. These transitions play a key role in numerous processes in the domain of nuclear physics. Equally important is their contribution in astrophysics, particularly in nuclear synthesis and supernova-explosions. In situations where allowed GT transitions are not favored, first-forbidden transitions become significant, specifically in medium heavy and heavy nuclei. For neutron-rich nuclei, first-forbidden transitions are favored mainly due to the phase-space amplification for these transitions. In this work we calculate the allowed GT as well as unique first-forbidden (U1F) |\DeltaJ| = 2 transitions strength in odd-odd and even-even nuclei in mass range 70\leq A \leq214. Two different pn-QRPA models were used with a schematic separable interaction to calculate GT and U1F transitions. The inclusion of U1F strength improved the overall comparison of calculated terrestrial \beta-decay half-lives in both models. The \textit{ft} values and reduced transition probabilities for the 2^-\longleftrightarrow 0^+ transitions were also calculated. We compared our calculations with the previously reported correlated RPA calculation and experimental results. Our calculations are in better agreement with measured data. For stellar applications we further calculated the allowed GT and U1F weak rates. These include \beta^{\pm}-decay rates and electron/positron capture rates of heavy nuclei in stellar matter. Our study shows that positron and electron capture rates command the total weak rates of these heavy nuclei at high stellar temperatures.
This thesis aims to elucidate the role of initial baryon stopping and its diffusion in heavy-ion collisions (HIC) using hydrodynamic model. In this regard, we have studied the observable-directed flow (v_1) of identified hadrons, particularly the v_1 of baryons and antibaryons, as well as the splitting observed between them in detail. We propose a new ansatz for the initial baryon distribution. By employing this initial baryon deposition model alongside a tilted energy distribution as inputs to a hybrid framework, we successfully describe the rapidity-odd v_1 of identified hadrons, including the elusive baryon-antibaryon splitting of v_1 across a wide range of \sqrt{s_{NN}}. Our model, incorporating baryon stopping and it's subsequent diffusion within a relativistic hydrodynamic framework and employing a crossover equation of state derived from lattice QCD calculations, establishes a non-critical baryonic baseline. Moreover, we demonstrate that recent STAR measurements of the centrality and system-size dependence of v_1 splitting between oppositely charged hadrons-attributed to electromagnetic field effects-are significantly influenced by background contributions from baryon stopping and its diffusion. Furthermore, we show that the rapidity dependence of the splitting of the rapidity-even component of v_1 between p and \bar{p} is highly sensitive to the initial baryon deposition scheme. If measured experimentally, this could constraint the rapidity dependence of the initial baryon deposition profile. Moreover, it could offer valuable phenomenological insights into the baryon junction picture and help refine constraints on the baryon diffusion coefficient of the medium. Notably, utilizing this phenomenologically successful baryon deposition model, we present the first estimation of the baryon diffusion coefficient for the strongly interacting QCD matter created in HIC.
Neural control variates (NCVs) have emerged as a powerful tool for variance reduction in Monte Carlo (MC) simulations, particularly in high-dimensional problems where traditional control variates are difficult to construct analytically. By training neural networks to learn auxiliary functions correlated with the target observable, NCVs can significantly reduce estimator variance while preserving unbiasedness. However, a critical but often overlooked aspect of NCV training is the role of autocorrelated samples generated by Markov Chain Monte Carlo (MCMC). While such samples are typically discarded for error estimation due to their statistical redundancy, they may contain useful information about the structure of the underlying probability distribution that can benefit the training process. In this work, we systematically examine the effect of using correlated configurations in training neural control variates. We demonstrate, both conceptually and numerically, that training on correlated data can improve control variate performance, especially in settings with limited computational resources. Our analysis includes empirical results from U(1) gauge theory and scalar field theory, illustrating when and how autocorrelated samples enhance NCV construction. These findings provide practical guidance for the efficient use of MCMC data in training neural networks.
We compute photon emission rate enhancement near the QCD critical point using an effective theory of dynamic critical phenomena and derive a universal photon spectrum. The emission rate scales similarly to conductivity, increasing with the correlation length (\xi), diverging at the critical point. The spectrum exhibits \omega dN_{\gamma}/d^3k \propto \omega^{-1/2} in the scaling regime, with the transition occurring at a frequency comparable to shear damping rate \omega \sim \gamma_{\eta}/\xi^2, reflecting the nonequilibrium properties of the near-critical liquid.
The effective mass and decay width of the \phi meson in the isospin asymmetric hot and dense resonance matter are studied using the effective Lagrangian framework considering the \phi K \bar K interactions at one-loop level. In addition to spin-1/2 octet baryons, we consider the effect of resonances \Delta^{++,+,0,-}, \Sigma^{*\pm,0},\Xi^{*0,-}, \Omega^{-}, on the properties of \phi meson. The in-medium effects on the \phi meson properties are simulated through the effective masses of kaons and antikaons computed using the chiral SU(3) hadronic mean field model in the presence of resonance baryons. The loop integral appearing in the computation of \phi meson self energies is regularized using the dipole form factor with a cutoff parameter. The presence of resonance baryons within the medium at finite temperature is observed to significantly modify the effective mass and decay width of \phi mesons. Examining the \phi meson masses and decay width within a dense medium is anticipated to be essential for understanding experimental results from heavy-ion collision experiments.
In this paper, we investigate the gluon distributions for the kaon and pion, as well as the improvement of the valence-quark distributions, in the framework of the gauge-invariant nonlocal chiral quark model (NL\chiQM), where the momentum dependence is taken into account. We then compute the gluon distributions for the kaon and pion that are dynamically generated from the splitting functions in the DGLAP QCD evolution. In a comparison with the recent lattice QCD and JAM global analysis results, it is found that the results for the pion gluon distributions at Q = 2 GeV, which is set based on the lattice QCD, have a good agreement with the recent lattice QCD data, which is followed by the up valence-quark distribution of the pion results at Q = 5.2 GeV in comparison with the reanalysis experimental data. Our prediction on the kaon gluon distributions at Q = 2 GeV is consistent with the recent lattice QCD calculation.
In our previous paper [JHEP 10 (2020) 144], we found that the complex Langevin (CL) method works for QCD at finite density on the 16^3 \times 32 lattice in the low-temperature high-density regime within the range \mu / T = 1.6 - 9.6 with \mu and T being the quark chemical potential and the temperature, which enabled us to see a clear trend towards the formation of the Fermi sphere. Here we investigate the validity of the CL method on the 24^3 \times 12 lattice in the deconfined phase near the deconfinement phase transition. As before, we use four-flavor staggered fermions and judge the validity using the criterion based on the probability distribution of the drift term. The spatial extent is L = (1.3 - 2.7 {\rm ~fm} )> \Lambda_{\rm LQCD}^{-1} \sim 1{\rm ~fm}, in contrast to our previous study with L < \Lambda_{\rm LQCD}^{-1}. We find that the CL method works in a broad region up to \mu / T = 4.8, while it starts to fail as we approach the phase boundary due to the singular drift problem, which can be understood qualitatively by extending the Banks-Casher relation to the case at finite density.
Nuclear modifications to collinear parton distribution functions are conventionally quantified by the ratios r^{\textrm{A}}_i(x,Q^2) = f^\textrm{A,proton}_i(x,Q^2) / f^\textrm{proton}_i(x,Q^2) . For a given nucleus A, these ratios generally depend on the parton momentum fraction x , the probing scale Q^{2} , and the parton species i . Determining these dependencies relies on a global analysis of diverse experimental data. However, in realistic observables, these dependencies are intricately intertwined, making their extraction challenging. In this paper, we propose a novel approach to effectively image the nuclear modification factors r^{\textrm{A}}_i(x,Q^2) at the observable level in proton-nucleus collisions at the Large Hadron Collider. Specifically, through a combined study of Z -boson production, Z +jet production, and Z+c -jet production, we separately enhance signals arising from light-quark, gluon, and heavy-flavor (charm) distributions in nuclei. This enables us to effectively image the r^{\textrm{A}}_i(x,Q^2) for specific parton species. The feasibility of this method is validated through perturbative calculations at next-to-leading order in the strong coupling constant, employing three sets of nuclear PDF parametrizations: EPPS21, nCTEQ15, and TUJU19. Future measurements of these observables are expected to provide more efficient constraints on nuclear PDFs and yield new insights into the detailed partonic structures of nuclei.