• Nuclear mass parabola and its applications(1709.05451)

Sept. 24, 2019 nucl-th
We propose a method to extract the properties of the isobaric mass parabola based on the total double $\beta$ decay energies of isobaric nuclei. Two important parameters of the mass parabola, the location of the most $\beta$-stable nuclei $Z_{A}$ and the curvature parameter $b_{A}$, are obtained for 251 A values based on the total double $\beta$ decay energies of nuclei compiled in AME2016 database.The advantage of this approach is that we can remove the pairing energy term $P_{A}$ caused by odd-even variation, and the mass excess $M(A,Z_{A})$ of the most stable nuclide for mass number $A$ in the performance process, which are used in the mass parabolic fitting method. The Coulomb energy coefficient $a_{c}=0.6910$ MeV is determined by the mass difference relation of mirror nuclei, and the symmetry energy coefficient is also studied by the relation $a_{\rm sym}(A)=0.25b_{A}Z_{A}$.
• QCD Evolution of Superfast Quarks(1511.06044)

June 29, 2019 hep-ph, nucl-ex, nucl-th
Recent high-precision measurements of nuclear deep inelastic scattering at high x and moderate 6 < Q$^2$ < 9GeV$^2$ give a rare opportunity to reach the quark distributions in the {\it superfast} region, in which the momentum fraction of the nucleon carried by its constituent quark is larger than the total fraction of the nucleon at rest, x>1. We derive the leading-order QCD evolution equation for such quarks with the goal of relating the moderate-Q$^2$ data to the two earlier measurements of superfast quark distributions at large 60 < Q$^2$ < 200~GeV$^2$. Since the high-Q$^2$ measurements gave strongly contradictory estimates of the nuclear effects that generate superfast quarks, relating them to the high-precision, moderate-Q$^2$ data through QCD evolution allows us to clarify this longstanding issue. Our calculations indicate that the moderate-Q$^2$ data at $x\lesssim 1.05$ are in better agreement with the high-Q$^2$ data measured in (anti)neutrino-nuclear reactions which require substantial high-momentum nuclear effects in the generation of superfast quarks. Our prediction for the high-Q$^2$ and x>1.1 region is somewhat in the middle of the neutrino-nuclear and muon-nuclear scattering data.
• The angular momentum dependence of nuclear optical potentials(1302.1097)

March 14, 2019 nucl-ex, nucl-th
The nuclear optical model potential (OMP) is generally assumed to be independent of the orbital angular momentum, $l$, of the interacting nuclei. Nucleon-nucleus and nucleus-nucleus interactions are customarily $l$ independent in calculations of nuclear elastic scattering and in standard reaction codes. The evidence for various forms of $l$ dependence of OMPs is reviewed and the importance of implementing these forms is evaluated. Existing arguments and evidence for $l$ dependence are reviewed and new arguments and calculations are introduced. The relationship is examined between (i) $l$ dependence, and, (ii) the undularity (waviness) of $l$-independent potentials that are $S$-matrix equivalent to $l$-dependent potentials. Such undularity is a property of the dynamic polarisation potential (DPP) generated by the coupling to reaction channels, or by coupling to excited states of the target or projectile nuclei. Various examples, particularly involving weakly bound projectile nuclei, are presented. Such undularity also occurs in $l$-independent potentials that have been found in model independent fits to precise wide angular range elastic scattering angular distributions. Cases of such phenomenological undularity, for both light and heavy ions, are referenced and shown to be related to undulatory properties of the dynamic polarisation potentials (DPPs) arising from channel coupling. Other forms of $l$ dependence, that could be standard options in direct reaction codes, are noted. Of particular importance are $l$ dependencies arising from full antisymmetrization. The case is made that reaction-induced $l$ dependence is a general property of nucleon-nucleus and nucleus-nucleus interactions and represents a valid extension of the nuclear optical model.
• Shear viscosity from nuclear stopping(1612.04874)

March 12, 2019 nucl-ex, nucl-th
Within a Boltzmann transport model, we demonstrate correlation between stopping observables and shear viscosity in central nuclear collisions at intermediate energies (on the order of 10 to 1000 MeV/nucleon). The correlation allows us to assess the viscosity of nuclear matter, by tuning the in-medium nucleon-nucleon cross section in our transport model to agree with nuclear stopping data. We also calculate the ratio of shear viscosity to entropy density to determine how close the system is to the universal quantum lower limit proposed in the context of ultrarelativistic heavy ion collisions.
• Standardized Cumulants of Flow Harmonic Fluctuations(1704.06295)

March 7, 2019 hep-ph, nucl-th
The distribution of flow harmonics in heavy ion experiment can be characterized by standardized cumulants. We first model the ellipticity and power parameters of the elliptic-power distribution by employing MC-Glauber model. Then we use the elliptic-power distribution together with the hydrodynamic linear response approximation to study the two dimensional standardized cumulants of elliptic and triangular flow ($v_2$ and $v_3$) distribution. For the second harmonic, it turns out that finding two dimensional cumulants in terms of $2q$-particle correlation functions $c_2\{2q\}$ is limited to the skewness. We also show that $c_3\{2\}$, $c_3\{4\}$, and $c_3\{6\}$, are related to the second, fourth, and sixth standardized cumulants of the $v_3$ distribution, respectively. The cumulant $c_{n}\{2q\}$ can be also written in terms of $v_n\{2q\}$. Specifically, $-(v_3\{4\}/v_3\{2\})^4$ turns out to be the kurtosis of the $v_3$ event-by-event fluctuation distribution. We introduce a new parametrization for the distribution $p(v_3)$ with $v_3\{2\}$, kurtosis and sixth-order standardized cumulant being its free parameters. Compared to the Gaussian distribution, it indicates a more accurate fit with experimental results. Finally, we compare the kurtosis obtained from simulation with that of extracted from experimental data for the $v_3$ distribution.
• Isolating the chiral magnetic effect from backgrounds by pair invariant mass(1705.05410)

March 3, 2019 nucl-ex, nucl-th
Topological gluon configurations in quantum chromodynamics induce quark chirality imbalance in local domains, which can result in the chiral magnetic effect (CME)--an electric charge separation along a strong magnetic field. Experimental searches for the CME in relativistic heavy ion collisions via the charge-dependent azimuthal correlator ($\Delta\gamma$) suffer from large backgrounds arising from particle correlations (e.g. due to resonance decays) coupled with the elliptic anisotropy. We propose differential measurements of the $\Delta\gamma$ as a function of the pair invariant mass ($m_{\rm inv}$), by restricting to high $m_{\rm inv}$ thus relatively background free, and by studying the $m_{\rm inv}$ dependence to separate the possible CME signal from backgrounds. We demonstrate by model studies the feasibility and effectiveness of such measurements for the CME search.
• Integral transforms of the quantum mechanical path integral: hit function and path averaged potential(1709.04984)

We introduce two new integral transforms of the quantum mechanical transition kernel that represent physical information about the path integral. These transforms can be interpreted as probability distributions on particle trajectories measuring respectively the relative contribution to the path integral from paths crossing a given spatial point (the hit function) and the likelihood of values of the line integral of the potential along a path in the ensemble (the path averaged potential).
• Beyond Wigner's isobaric multiplet mass equation: Effect of charge-symmetry-breaking interaction and Coulomb polarization(1601.05167)

Feb. 27, 2019 nucl-th
The quadratic form of the isobaric multiplet mass equation (IMME), which was originally suggested by Wigner and has been generally regarded as valid, is seriously questioned by recent high-precision nuclear mass measurements. The usual resolution to this problem is to add empirically the cubic and quartic $T_z$-terms to characterize the deviations from the IMME, but finding the origin of these terms remains an unsolved difficulty. Based on a strategy beyond the Wigner's first-order perturbation, we derive explicitly the cubic and quartic $T_z$-terms. These terms are shown to be generated by the effective charge-symmetry breaking and charge-independent breaking interactions in nuclear medium combined with the Coulomb polarization effect. Calculations for the $sd$- and lower $fp$-shells explore a systematical emergence of the cubic $T_z$-term, suggesting a general deviation from the original IMME. Intriguingly, the magnitude of the deviation exhibits an oscillation-like behavior with mass number, modulated by the shell effect.
• Probing the hot and dense nuclear matter with $K^*,\bar{K}^*$ vector mesons(1707.00060)

Feb. 25, 2019 hep-ph, nucl-th
We investigate probing the hot and dense nuclear matter with strange vector mesons ($K^*, \bar{K}^*$). Our analysis is based on PHSD which incorporates partonic and hadronic dof and describes the full dynamics of HICs. This allows to study the $K^*$ and $\bar{K}^*$ meson formation from the QGP and the in-medium effects related to the modification of their properties during the propagation in dense and hot matter. We employ relativistic Breit-Wigner spectral functions for the $K^*,\bar{K}^*$ mesons with self-energies obtained from a G-matrix approach to study the role of in-medium effects on the $K^*$ and $\bar{K}^*$ meson dynamics in HIC from FAIR/NICA to LHC energies. According to our analysis most of the final $K^*/\bar{K}^*$s, that can be observed experimentally, are produced during the late hadronic phase and stem dominantly from the $K (\bar{K}) + \pi \to K^*(\bar{K}^*)$ formation channel. The amount of $K^*/\bar{K}^*$s originating from the QGP channel is comparatively small even at LHC energies and such $K^*/\bar{K}^*$s can hardly be reconstructed experimentally due to the rescattering of final pions and (anti-)kaons. This mirrors the results from our previous study on the strange vector-meson production in HICs at RHIC energies. The influence of the in-medium effects on the dynamics of the $K^*/\bar{K}^*$ is rather small since they are mostly produced at low baryon densities. Additional cuts on the shape of the observed signal and the range of the invariant mass region of the $K^*/\bar{K}^*$ also affect the final spectra. We demonstrate that the $K^*/\bar{K}^*$ in-medium effects are more visible at lower beam energy, e.g. FAIR/NICA and BES RHIC energies, where the production of $K^*/\bar{K}^*$s occurs at larger baryon densities. Finally, we present the experimental procedures to extract information on the in-medium masses and widths by fitting final mass spectra at LHC energies.
• Tritium $\beta$-decay in pionless effective field theory(1611.10004)

Nov. 15, 2019 nucl-th
We calculate the $\beta$-decay of tritium at next-to-leading order in pionless effective field theory. At this order, a low-energy parameter $L_{1, A}$ enters the calculation that is also relevant for a high-accuracy prediction of the solar proton-proton fusion rate. We use the tritium half-life to determine this parameter and provide uncertainty estimates. We show proper renormalization of our calculation analyzing the residual cutoff dependence of observables. We find that next-to-leading order corrections contribute about $4\%$ to the triton decay Gamow-Teller strength. We show that these conclusions are insensitive to different arrangements of the effective range expansion.
• The relationship between undularity and L dependence of the nuclear optical model potential(1705.07003)

Feb. 14, 2019 nucl-th
The contribution of collective or reaction channels to a local optical model potential, OMP, can be readily calculated as a dynamical polarization potential, DPP. The resulting local DPPs commonly have undulatory (wavy') features, often including local regions of emissivity in the imaginary component. We show here that this undularity arises from $l$-dependence of the underlying formal non-local and $l$-dependent DPP. The $l$-independent proton OMPs, that have the same $S$-matrix $S_{lj}$ as phenomenological $l$-dependent potentials, exhibit undulations that are qualitatively similar to undulations of local DPPs generated by channel coupling. The $l$-dependent phenomenological potentials studied are the potentials that give the best existing fits to the relevant elastic scattering data and the undulatory potentials presented here, being $S$-matrix equivalent (i.e.\ having the same $S$-matrix, $S_{lj}$) give exactly the same scattering. In addition, we present calculations strongly suggesting that undularity (waviness') is a generic property of $l$-independent potentials that are $S$-matrix equivalent to $l$-dependent potentials. Implications for the validity of folding models based on a local density model are noted.
• Jet axes and universal transverse-momentum-dependent fragmentation(1612.04817)

Feb. 8, 2019 hep-ph, nucl-th
We study the transverse momentum spectrum of hadrons in jets. By measuring the transverse momentum with respect to a judiciously chosen axis, we find that this observable is insensitive to (the recoil of) soft radiation. Furthermore, for small transverse momenta we show that the effects of the jet boundary factorize, leading to a new transverse-momentum-dependent (TMD) fragmentation function. In contrast to the usual TMD fragmentation functions, it does not involve rapidity divergences and is universal in the sense that it is independent of the type of process and number of jets. These results directly apply to sub-jets instead of hadrons. We discuss potential applications, which include studying nuclear modification effects in heavy-ion collisions and identifying boosted heavy resonances.
• Iso-vector axial form factors of the nucleon in two-flavour lattice QCD(1705.06186)

Jan. 30, 2019 hep-ph, nucl-th, hep-lat
We present a lattice calculation of the nucleon iso-vector axial and induced pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical flavours of non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions and an $\mathcal{O}(a)$-improved axial current together with the pseudoscalar density. Excited-state effects in the extraction of the form factors are treated using a variety of methods, with a detailed discussion of their respective merits. The chiral and continuum extrapolation of the results is performed both using formulae inspired by Heavy Baryon Chiral Perturbation Theory (HBChPT) and a global approach to the form factors based on a chiral effective theory (EFT) including axial vector mesons. Our results indicate that careful treatment of excited-state effects is important in order to obtain reliable results for the axial form factors of the nucleon, and that the main remaining error stems from the systematic uncertainties of the chiral extrapolation. As final results, we quote $g_{\rm A} = 1.278 \pm 0.068\genfrac{}{}{0pt}{1}{+0.000}{-0.087}$, $\langle r_{\rm A}^2\rangle = 0.360 \pm 0.036\genfrac{}{}{0pt}{1}{+0.080}{-0.088}~\mathrm{fm}^2$, and $g_{\rm P} = 7.7 \pm 1.8 \genfrac{}{}{0pt}{1}{+0.8}{-2.0}$ for the axial charge, axial charge radius and induced pseudoscalar charge, respectively, where the first error is statistical and the second is systematic.
• Pinning down the linearly-polarised gluons inside unpolarised protons using quarkonium-pair production at the LHC(1710.01684)

Jan. 29, 2019 hep-ph, hep-ex, nucl-ex, nucl-th
We show that the production of J/psi or Upsilon pairs in unpolarised pp collisions is currently the best process to measure the momentum-distribution of linearly-polarised gluons inside unpolarised protons through the study of azimuthal asymmetries. Not only the short-distance coefficients for such reactions induce the largest possible cos 4 phi modulations, but analysed data are already available. Among the various final states previously studied in unpolarised pp collisions within the TMD factorisation approach, di-J/psi production exhibits by far the largest asymmetries, up to 50% in the region studied by the ATLAS and CMS experiments. In addition, we use the very recent LHCb data at 13 TeV to perform the first fit of the unpolarised transverse-momentum-dependent gluon distribution.
• Thermodynamic and transport properties in Au + Au collisions at RHIC energies from the clustering of color strings(1708.06689)

Jan. 29, 2019 hep-ph, hep-ex, nucl-ex, nucl-th, hep-lat
In this work, we have extracted the initial tempearture from the transverse momentum spectra of charged particles in $Au+Au$ collisions using STAR data at RHIC energies from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. The initial energy density ($\varepsilon$), shear viscosity to entropy density ratio ($\eta/s$), trace anomaly ($\Delta$), the squared speed of sound ($C_s^2$), entropy density, and bulk viscosity to entropy density ratio ($\zeta/s$) are obtained and compared with the lattice QCD calculations for (2+1) flavor. The initial temperatures obtained are compared with various hadronization and chemical freeze-out temperatures. The analysis of the data shows that the deconfinement to confinement transition possibly takes place between $\sqrt{s_{NN}}$ = 11.5 and 19.6 GeV.
• Addendum: Muon-neutrino-induced charged current pion production(1708.04528)

March 27, 2019 hep-ex, nucl-ex, nucl-th
This short paper is an addendum to a recent publication on charged current neutrino-induced pion production (Phys.\ Rev.\ C96 (2017) no.1, 015503). It presents comparisons of pion production cross sections measured at the T2K near detector for a CH target.
• Features of $\omega$ photoproduction off nucleon target at forward angles : Dominance of $\pi$ exchange with Regge cuts and scaling of differential cross sections(1710.04511)

Jan. 16, 2019 hep-ph, nucl-th
We investigate photoproduction of $\omega$ off a nucleon target $\gamma N\to\omega N$ by using a Reggeized model where $\pi(135)+\sigma(500)+f_1(1285)+f_2(1270)$+Pomeron exchanges are included in the $t$-channel for the reaction at forward angles. The reaction mechanism at low energy is featured by the dominance of the $\pi$ exchange with the absorptive cuts introduced to modulate the pion contribution to both $\gamma p\to \omega p$ and $\gamma n\to \omega n$ reactions. Necessity of $\sigma$ exchange is illustrated in the analysis of the cross section from natural parity exchanges. Cross sections for differential, total, and spin density matrix with beam polarization are reproduced and compared with existing data on $\gamma p\to \omega p$. Scaled differential cross sections of Jefferson Lab data on $\gamma p\to \omega p$ are investigated at the production angle $\theta=90^\circ$ with the nonlinear trajectories for a saturation. Differential and total cross sections for the $\gamma n\to\omega n$ reaction are analyzed to compare with recent experimental data at the CBELSA/TAPS Collaboration. An application to the $\gamma p\to\omega\Delta^+$ reaction is presented to demonstrate the significance of the $\pi$ exchange in the total and differential cross sections.
• Off-forward gluonic structure of vector mesons(1703.08220)

Jan. 16, 2019 hep-ph, nucl-th, hep-lat
The spin-independent and transversity generalised form factors (GFFs) of the $\phi$ meson are studied using lattice QCD calculations with light quark masses corresponding to a pion mass $m_\pi\sim450(5)$ MeV. One transversity and three spin-independent GFFs related to the lowest moments of leading-twist spin-independent and transversity gluon distributions are obtained at six non-zero values of the momentum transfer up to 1.2 GeV$^2$. These quantities are compared with the analogous spin-independent quark GFFs and the electromagnetic form factors determined on the same lattice ensemble. The results show quantitative distinction between the spatial distribution of transversely polarised gluons, unpolarised gluons, and quarks, and point the way towards further investigations of the gluon structure of nucleons and nuclei.
• Improved version of the PHOBOS Glauber Monte Carlo(1408.2549)

Jan. 14, 2019 nucl-ex, nucl-th
Glauber models are used to calculate geometric quantities in the initial state of heavy ion collisions, such as impact parameter, number of participating nucleons and initial eccentricity. Experimental heavy-ion collaboration, in particular at RHIC and LHC, use Glauber Model calculations for various geometric observables. In this document, we describe the assumptions inherent to the approach, and provide an updated implementation (v2) of the Monte Carlo based Glauber Model calculation, which originally was used by the PHOBOS collaboration. The main improvement w.r.t. the earlier version (arXiv:0805.4411) are the inclusion of Tritium, Helium-3, and Uranium, as well as the treatment of deformed nuclei and Glauber-Gribov fluctuations of the proton in p+A collisions. A users' guide (updated to reflect changes in v2) is provided for running various calculations.
• Neutrinoless double-$\beta$ decay of $^{124}$Sn, $^{130}$Te, and $^{136}$Xe in the Hamiltonian-based generator-coordinate method(1709.05313)

Jan. 7, 2019 nucl-th
We present a generator-coordinate method for realistic shell-model Hamiltonians that closely approximates the full shell model calculations of the matrix elements for the neutrinoless double-$\beta$ decay of $^{124}$Sn, $^{130}$Te, and $^{136}$Xe. We treat not only quadrupole deformations but also the proton-neutron pairing amplitudes as generator coordinates. We validate this method by calculating and comparing spectroscopic quantities with the exact shell model results and experimental data. Our Hamiltonian-based generator-coordinate method produces $0\nu\beta\beta$ matrix elements much closer to the shell model ones, compared to the existing energy-density-functional-based generator-coordinate approaches. The remaining overestimation of $0\nu\beta\beta$ nuclear matrix element suggests that additional correlations may be needed to be taken into account for $^{124}$Sn, $^{130}$Te, and $^{136}$Xe when calculating with the Hamiltonian-based generator-coordinate method. The validation of this method may open the possibility of calculating $0\nu\beta\beta$ matrix element of $^{150}$Nd in a large shell-model space.
• Kondo phase diagram of quark matter(1604.07208)

Jan. 7, 2019 hep-ph, nucl-th, cond-mat.str-el
We discuss the ground state of a quark matter containing heavy quarks as impurities in a simple model which exhibits the QCD Kondo effect. The model includes a current-current interaction with the color exchange between a light quark ($\psi$) and a heavy quark ($\Psi$). We introduce a gap function $\Delta \sim \langle \bar \psi \Psi \rangle$ which represents the correlation between $\psi$ and $\Psi$, and perform the mean-field approximation assuming that heavy quarks are uniformly distributed. Values of the gap $\Delta$ measure the strength of mixing between $\psi$ and $\Psi$. The gap equation obtained from the minimum of the thermodynamical potential together with the condition for the heavy-quark number conservation turns out to allow for nonzero values of the gap as the most stable state. We draw a phase diagram in $\mu$ (the light-quark chemical potential) and $\lambda$ (an analog of the heavy-quark chemical potential) plane, and identify the region where the QCD Kondo effect occurs.
• $D D^{*}$ potentials in chiral perturbation theory and possible molecular states(1708.06918)

Jan. 6, 2019 hep-ph, hep-ex, nucl-th, hep-lat
The $DD^{*}$ potentials are studied within the framework of heavy meson chiral effective field theory. We have obtained the effective potentials of the $DD^{*}$ system up to $O(\epsilon^2)$ at one loop level. In addition to the one-pion exchange contribution, the contact and two-pion exchange interactions are also investigated in detail. Furthermore, we have searched for the possible molecular states by solving Schr\"odinger equation with the potentials. We notice that the contact and two-pion exchange potentials are non-negligible numerically and important for the existence of a bound state. In our results, no bound state is founded in the $I=0$ channel within a wide range of cutoff parameter, while there exists a bound state in the $I=1$ channel as cutoff is near $m_\rho$ in our approach.
• Transverse expansion of hot magnetized Bjorken flow in heavy ion collisions(1710.01037)

June 7, 2019 hep-ph, nucl-th
We argue that the existence of an inhomogeneous external magnetic field can lead to radial flow in transverse plane. Our aim is to show how the introduction of a magnetic field generalizes the Bjorken flow. We investigate the effect of an inhomogeneous weak external magnetic field on the transverse expansion of in-viscid fluid created in high energy nuclear collisions. In order to simplify our calculation and compare with Gubser model, we consider the fluid under investigation to be produced in central collisions, at small impact parameter; azimuthal symmetry has been considered. In our model, we assume an inhomogeneous external magnetic field following the power-law decay in proper time and having radial inhomogeneity perpendicular to the radial velocity of the in-viscid fluid in the transverse plane; then the space time evolution of the transverse expansion of the fluid is obtained. We also show how the existence of an inhomogeneous external magnetic field modifies the energy density. Finally we use the solutions for the transverse velocity and energy density in the presence of a weak magnetic field, to estimate the transverse momentum spectrum of protons and pions emerging from the Magneto-hydrodynamic solutions.
• $\Psi(2s)$ production in $p+A$ collisions(1703.04580)

Dec. 27, 2018 hep-ph, nucl-th
We have analyzed the data on $\psi'$ production in proton-nucleus ($p+A$) collisions, available from the NA50 Collaboration in the SPS energy domain. The investigated data sets include the absolute production cross sections as well as $\psi'$-to-Drell Yan (DY) cross section ratios. An adapted version of two component QVZ model has been employed to calculate $\psi'$ production cross sections. For both $\psi'$ and DY production, nuclear modifications to the free nucleon structure functions are taken into account. For $\psi'$, final state interaction of the produced $c\bar{c}$ pairs with the nuclear medium is also taken into account, in accordance with the previously analyzed $J/\psi$ data. A satisfactory description of the data in $p+A$ collisions is obtained. Model calculations are extrapolated to low energy regime to predict the $\psi'$ suppression in proton induced collisions at near threshold energies.
• Neutron skin and halo in medium and heavy nuclei within the extended Thomas-Fermi theory(1505.04625)

Dec. 17, 2018 nucl-th
The neutron skin and halo distributions in medium and heavy nuclei are calculated within the extended Thomas-Fermi approximation. Calculations are carried out for the effective Skyrme-like forces using the direct variational method. The analytical expression for the isovector shift of the rms radii $\Delta r_{np}$ as a sum of skin- and halo-like terms is obtained. The contribution of halo and skin terms to $\Delta r_{np}$ are found to be approximately equal.