# Nuclear Theory (nucl-th)

• We present a symmetry-based approach for shape coexistence in nuclei, founded on the concept of partial dynamical symmetry (PDS). The latter corresponds to a situation when only selected states (or bands of states) of the coexisting configurations preserve the symmetry while other states are mixed. We construct explicitly critical-point Hamiltonians with two or three PDSs of the type U(5), SU(3), ${\overline{\rm SU(3)}}$ and SO(6), appropriate to double or triple coexistence of spherical, prolate, oblate and $\gamma$-soft deformed shapes, respectively. In each case, we analyze the topology of the energy surface with multiple minima and corresponding normal modes. Characteristic features and symmetry attributes of the quantum spectra and wave functions are discussed. Analytic expressions for quadrupole moments and $E2$ rates involving the remaining solvable states are derived and isomeric states are identified by means of selection rules.
• An optical potential constructed using the ${\bar p}N$ scattering amplitudes derived from the 2009 version of the Paris ${\bar N}N$ potential is applied in calculations of ${\bar p}$ quasi-bound states in selected nuclei across the periodic table. A proper self-consistent procedure for treating energy dependence of the amplitudes in a nucleus appears crucial for evaluating ${\bar p}$ binding energies and widths. Particular attention is paid to the role of $P$-wave amplitudes. While the $P$-wave potential nearly does not affect calculated ${\bar p}$ binding energies, it reduces considerably the corresponding widths. The Paris $S$-wave potential supplemented by a phenomenological $P$-wave term yields in dynamical calculations ${\bar p}$ binding energies $B_{\bar p}\approx 200$~MeV and widths $\Gamma_{\bar p}\sim~200 - 230$~MeV, which is very close to the values obtained within the RMF model consistent with ${\bar p}$-atom data.
• The existence of neutron stars with $2M_\odot$ requires the strong stiffness of the equation of state (EoS) of neutron-star matter. We introduce a multi-pomeron exchange potential (MPP) working universally among 3- and 4-baryons to stiffen the EoS. Its strength is restricted by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoSs are derived using the Brueckner-Hartree-Fock (BHF) and the cluster variational method (CVM) with the nuclear interactions ESC and AV18. The mass-radius relations are derived by solving the Tolmann-Oppenheimer-Volkoff (TOV) equation, where the maximum masses over $2M_\odot$ are obtained on the basis of the terrestrial data. Neutron-star radii $R$ at a typical mass $1.5M_\odot$ are predicted to be $12.3\!\sim\!13.0$ km. The uncertainty of calculated radii is mainly from the ratio of 3- and 4-pomeron coupling constants, which cannot be fixed by any terrestrial experiment. Though values of $R(1.5M_\odot)$ are not influenced by hyperon-mixing effects, finely-observed values for them indicate degrees of EoS softening by hyperon mixing in the region of $M\!\sim\!2M_\odot$. If $R(1.5M_\odot)$ is less than about 12.4 km, the softening of EoS by hyperon mixing has to be weak. Useful information can be expected by the space mission NICER offering precise measurements for neutron-star radii within $\pm 5\%$.
• We study $\beta$ decays within an effective theory that treats nuclei as a spherical collective core with an even number of neutrons and protons that can couple to an additional neutron and/or proton. First we explore Gamow-Teller $\beta$ decays of parent odd-odd nuclei into low-lying ground-, one-, and two-phonon states of the daughter even-even system. The low-energy constants of the effective theory are adjusted to data on $\beta$ decays to ground states or Gamow-Teller strengths. The corresponding theoretical uncertainty is estimated based on the power counting of the effective theory. For a variety of medium-mass and heavy isotopes the theoretical matrix elements are in good agreement with experiment within the theoretical uncertainties. We then study the two-neutrino double-$\beta$ decay into ground and excited states. The results are remarkably consistent with experiment within theoretical uncertainties, without the necessity to adjust any low-energy constants.
• A finite-range separable pairing force is implemented in the tilted axis cranking covariant density functional theory for describing rotational bands in open shell nuclei. The developed method is used to investigate the yrast sequence of 109Ag for an example. The experimental energy spectrum, angular momenta, and electromagnetic transition probabilities are well reproduced by taking into account pairing correlations with the separable paring force. An abrupt transition of the rotational axis from the long-intermediate plane to the long-short one is obtained and discussed in detail.
• A Fermi liquid model for hadrons has been suggested for the hadrons in medium. The hadrons are supposed to behave like quasi particle as Fermi excitation while in the medium and the effective mass of the hadrons have been estimated using Fermi liquid model. Considering a momentum dependent potential inside the medium to describe the interaction, the effective masses of the hadrons are estimated. The temperature dependence of effective masses has also been studied. The possibility of describing masses of the quarks as Fermi excitation has been investigated. Compressibility, specific heats, density of states in medium has been studied. The potential depth for light and singly heavy baryons in medium has been extracted. The results are found to be very interesting and compared with the other studies available in literature.
• Oscillation modes of isolated compact stars can, in principle, be a fingerprint of the equation of state (EoS) of dense matter. We study the non-radial high-frequency l=2 spheroidal modes of neutron stars and strange quark stars, adopting a two-component model (core and crust) for these two types of stars. Using perturbed fluid equations in the relativistic Cowling approximation, we explore the effect of a strangelet or hadronic crust on the oscillation modes of strange stars. The results differ from the case of neutron stars with a crust. In comparison to fluid-only configurations, we find that a solid crust on top of a neutron star increases the p-mode frequency slightly with little effect on the f-mode frequency, whereas for strange stars, a strangelet crust on top of a quark core significantly increases the f-mode frequency with little effect on the p-mode frequency.
• Background: The interacting boson model (IBM) has been used extensively to calculate the matrix elements governing neutrinoless double-beta decay. Studies within other models indicate that a good description of neutron-proton pairing is essential for accurate calculations of those matrix elements. The usual interacting boson model is based only on like-particle pairs, however, and the extent to which it captures neutron-proton pairing is not clear. Purpose: To determine whether neutron-proton pairing should be explicitly included as neutron-proton bosons in IBM calculations of neutrinoless double-beta decay matrix elements. Method: An isospin-invariant version of the nucleon-pair shell model is applied to carry out shell-model calculations in a large space and in a collective subspace, and to define effective operators in the latter. A democratic mapping is then used to define corresponding boson operators for the IBM, with and without an isoscalar neutron-proton pair boson. Results: IBM calculations with and without the isoscalar boson are carried out for nuclei near the beginning of the $pf$ shell, with a realistic shell-model Hamiltonian and neutrinoless double-beta-decay operator as the starting point. Energy spectra and double-beta matrix elements are compared to those obtained in the underlying shell model. Conclusions: The isoscalar boson does not improve energy spectra but does improve double-beta matrix elements. To be useful at the level of precision we need, the mapping procedure must be further developed to better determine the dependence of the boson Hamiltonian and decay operator on particle number and isospin. But the benefits provided by the isoscalar boson suggest that through an appropriate combination of mappings and fitting, it would make IBM matrix elements more accurate for the heavier nuclei used in experiments.
• We perform quantum calculations of fluctuations of the electromagnetic fields in $AA$ collisions at RHIC and LHC energies. We find that in the quantum picture the field fluctuations are much smaller than predictions of the classical Monte-Carlo simulation with the Woods-Saxon nuclear density.
• High precision data of lepton angular distributions for $\gamma^*/Z$ production in $pp$ collisions at the LHC, covering broad ranges of dilepton transverse momenta ($q_T$) and rapidity ($y$), were recently reported. Strong $q_T$ dependencies were observed for several angular distribution coefficients, $A_i$, including $A_0 - A_4$. Significant $y$ dependencies were also found for the coefficients $A_1$, $A_3$ and $A_4$, while $A_0$ and $A_2$ exhibit very weak rapidity dependence. Using an intuitive geometric picture we show that the $q_T$ and $y$ dependencies of the angular distributions coefficients can be well described.
• We have computed the isospin and spin averaged cross sections of the processes $\pi K^*\to \rho K$ and $\rho K^*\to \pi K$, which are crucial in the determination of the abundances of $K^*$ and $K$ in heavy ion collisions. Improving previous calculations, we have considered several mechanisms which were missing, such as the exchange of axial and vector resonances ($K_1(1270)$, $K^*_2(1430)$, $h_1(1170)$, etc...) and also other processes such as $\pi K^*\to \omega K, \phi K$ and $\omega K^*,\,\phi K^*\to \pi K$. We find that some of these mechanisms give important contributions to the cross section. Our results also suggest that, in a hadron gas, $K^*$ production might be more important than its absorption.
• We discuss how photoproduction and electroproduction processes will provide a precise tool to study nuclear structure in the small-$x$ region in the future Electron-Ion Collider. In particular, we emphasize how exclusive vector meson production at the EIC can be used to accurately map the gluonic density profile of nuclei, including its event-by-event fluctuations. In addition, we show how measurements of the total diffractive cross section will be sensitive to saturation effects.
• A new approach is described to help improve the foundations of relativistic viscous fluid dynamics and its coupling to general relativity. Focusing on neutral conformal fluids constructed solely in terms of hydrodynamic variables, we derive the most general viscous energy-momentum tensor yielding equations of motion of second order in the derivatives, which is shown to provide a novel type of resummation of the relativistic Navier-Stokes equations for which causality holds. We rigorously prove existence, uniqueness, and causality of solutions of this theory (in the full nonlinear regime) both in a Minkowski background and also when the fluid is dynamically coupled to Einstein's equations. Linearized disturbances around equilibrium in Minkowski spacetime are stable in this causal theory and the second law of thermodynamics is satisfied. A numerical study reveals the presence of a hydrodynamic attractor for a rapidly expanding fluid. Further properties are also studied and a brief discussion of how this approach can be generalized to non-conformal fluids is presented.