Classical model for diffusion and thermalization of heavy quarks in a hot medium: memory and out-of-equilibrium effects

We consider a simple model for the diffusion of heavy quarks in a hot bath,modeling the latter by an ensemble of oscillators distributed according to either a thermal distribution or to an out-of-equilibrium distribution with a saturation scale.In this model it is easy to introduce memory effects by changing the distribution of oscillators:we model them by introducing a Gaussian distribution,dN/dω.which can be deformed continuously from a δ-function,giving a Markov dissipation, to a broad kernel with memory.Deriving the equation of motion of the heavy quark in the bath,we remark how dissipation comes out naturally as an effect of the back-reaction of the oscillators on the bath.Moreover,the exact solution of this equation allows to define the thermalization time as the time necessary to remove any memory of the initial conditions.We find that the broadening of the dissipative kernel,while keeping the coupling fixed,lowers the thermalization time.We also derive the fluctuation-dissipation theorem for the bath,and use it to estimate the kinematic regime in which momentum diffusion of the heavy quark dominates over drift.We find that diffusion is more important as long as K_0/ε is small,where K_0 and ε denote the initial energy of the heavy quark and the average energy of the bath,respectively.     

New behaviors of α-particle preformation factors near doubly magic 100Sn

The \begin{document}$ \alpha $\end{document}-particle preformation factors of nuclei above doubly magic nuclei \begin{document}$ ^{100} $\end{document}Sn and \begin{document}$ ^{208} $\end{document}Pb are investigated within the generalized liquid drop model. The results show that the \begin{document}$ \alpha $\end{document}-particle preformation factors of nuclei near self-conjugate doubly magic \begin{document}$ ^{100} $\end{document}Sn are significantly larger than those of analogous nuclei just above \begin{document}$ ^{208} $\end{document}Pb, and they will be enhanced as the nuclei move towards the \begin{document}$ N = Z $\end{document} line. The proton–neutron correlation energy \begin{document}$ E_{p-n} $\end{document} and two protons–two neutrons correlation energy \begin{document}$ E_{2p-2n} $\end{document} of nuclei near \begin{document}$ ^{100} $\end{document}Sn also exhibit a similar situation, indicating that the interactions between protons and neutrons occupying similar single-particle orbitals could enhance the \begin{document}$ \alpha $\end{document}-particle preformation factors and result in superallowed \begin{document}$ \alpha $\end{document} decay. This also provides evidence of the significant role of the proton–neutron interaction on \begin{document}$ \alpha $\end{document}-particle preformation. Also, the linear relationship between \begin{document}$ \alpha $\end{document}-particle preformation factors and the product of valence protons and valence neutrons for nuclei around \begin{document}$ ^{208} $\end{document}Pb is broken in the \begin{document}$ ^{100} $\end{document}Sn region because the \begin{document}$ \alpha $\end{document}-particle preformation factor is enhanced when a nucleus near \begin{document}$ ^{100} $\end{document}Sn moves towards the \begin{document}$ N = Z $\end{document} line. Furthermore, the calculated \begin{document}$ \alpha $\end{document} decay half-lives fit well with the experimental data, including the recent observed self-conjugate nuclei \begin{document}$ ^{104} $\end{document}Te and \begin{document}$ ^{108} $\end{document}Xe [Phys. Rev. Lett. 121, 182501 (2018)].     

Radial excited heavy mesons

In this study, the first radial excited heavy pseudoscalar and vector mesons (\begin{document}$\eta_c(2S)$\end{document}, \begin{document}$\psi(2S)$\end{document}, \begin{document}$B_c(2S)$\end{document}, \begin{document}$B^*_c(2S)$\end{document}, \begin{document}$\eta_b(2S)$\end{document}, and \begin{document}$\varUpsilon(2S)$\end{document}) are investigated using the Dyson-Schwinger equation and Bethe-Salpeter equation approach. It is shown that the effective interactions of the radial excited states are harder than those of the ground states. With the interaction well determined by fitting the masses and leptonic decay constants of \begin{document}$\psi(2S)$\end{document} and \begin{document}$\varUpsilon(2S)$\end{document}, the first radial excited heavy mesons could be quantitatively described in the rainbow ladder approximation. The masses and leptonic decay constants of \begin{document}$\eta_c(2S)$\end{document}, \begin{document}$B_c(2S)$\end{document}, \begin{document}$B^*_c(2S)$\end{document}, and \begin{document}$\eta_b(2S)$\end{document} are predicted.     

Production of bottomonium-like Z_b states in e-h and ultraperipheral h-h collisions

The photoproduction of the bottomonium-like states Z_b(10610) and Z_b(10650) via γ p scattering is studied within an effective Lagrangian approach and the vector-meson-dominance model. The Regge model is employed to calculate the photoproduction of Z_b states via the t-channel with π exchange. The numerical results show that the values of the total cross-sections of Z_b(10610) and Z_b(10650) can reach 0.09 nb and 0.02 nb, respectively,near the center-of-mass energy of 22 GeV. Experimental measurements and studies of the photoproduction of Z_b states near the energy region around W ■ 22 GeV are suggested. Moreover, with the help of eSTARlight and STARlight programs, we have obtained the cross-sections and numbers of events for Z_b(10610) production in electron-ion collisions(EIC) and ultraperipheral collisions(UPCs). The results show that a considerable number of Z_b(10610)events can be produced in the relevant experiments on EICs and UPCs. We have also calculated the rates and kinematic distributions for γp → Z_bn in ep and pA collisions via EICs and UPCs. The results will provide an important reference for the RHIC, LHC, EIC-US, LHeC, and FCC experiments in searching for bottomonium-like Z_b states.     

Systematic calculations of cluster radioactivity half-lives in trans-lead nuclei

In the present work, based on the Wentzel-Kramers-Brillouin (WKB) theory, considering the cluster preformation probability (\begin{document}$ P_{c} $\end{document}), we systematically investigate the cluster radioactivity half-lives of 22 trans-lead nuclei ranging from ²²¹Fr to ²⁴²Cm. When the mass number of the emitted cluster \begin{document}$ A_{c} $\end{document} \begin{document}$ < $\end{document} 28, \begin{document}$P_{c} $\end{document} is obtained by the exponential relationship of \begin{document}$ P_{c} $\end{document} to the α decay preformation probability (\begin{document}$ P_{\alpha} $\end{document}) proposed by R. Blendowskeis \begin{document}$ et $\end{document} \begin{document}$ al. $\end{document} [Phys. Rev. Lett. 61, 1930 (1988)], while \begin{document}$ P_{\alpha} $\end{document} is calculated through the cluster-formation model (CFM). When \begin{document}$ A_{c} $\end{document} \begin{document}$ \ge $\end{document} 28, \begin{document}$ P_{c} $\end{document} is calculated through the charge-number dependence of \begin{document}$ P_{c} $\end{document} on the decay products proposed by Ren \begin{document}$ et $\end{document} \begin{document}$ al. $\end{document} [Phys. Rev. C 70, 034304 (2004)]. The half-lives of cluster radioactivity have been calculated by the density-dependent cluster model [Phys. Rev. C 70, 034304 (2004)] and by the unified formula of half-lives for alpha decay and cluster radioactivity [Phys. Rev. C 78, 044310 (2008)]. For comparison, a universal decay law (UDL) proposed by Qi \begin{document}$ et $\end{document} \begin{document}$ al. $\end{document} [Phys. Rev. C 80, 044326 (2009)], a semi-empirical model for both α decay and cluster radioactivity proposed by Santhosh [J. Phys. G: Nucl. Part. Phys. 35, 085102 (2008)], and a unified formula of half-lives for alpha decay and cluster radioactivity [Phys. Rev. C 78, 044310 (2008)] are also used. The calculated results of our work, Ni's formula , and the UDL can well reproduce the experimental data and are better than those of Santhosh's model. In addition, we extend this model to predict the half-lives for 51 nuclei, whose cluster radioactivity is energetically allowed or observed but not yet quantified in NUBASE2020.     

Investigating S-wave bound states composed of two pseudoscalar mesons

In this study, we systematically investigated two-pseudoscalar meson systems with the Bethe-Salpeter equation in the ladder and instantaneous approximations. By solving the Bethe-Salpeter equation numerically with the kernel containing the one-particle exchange diagrams, we found that the \begin{document}$ K\bar{K} $\end{document}, \begin{document}$ DK $\end{document}, \begin{document}$ B\bar{K} $\end{document}, \begin{document}$ D\bar{D} $\end{document}, \begin{document}$ B\bar{B} $\end{document}, \begin{document}$ BD $\end{document}, \begin{document}$ D\bar{K} $\end{document}, \begin{document}$ BK $\end{document}, and \begin{document}$ B\bar{D} $\end{document} systems with \begin{document}$ I=0 $\end{document} can exist as bound states. We also studied the contributions from heavy meson (\begin{document}$ J/\psi $\end{document} and \begin{document}$\Upsilon $\end{document}) exchanges and found that the contributions from heavy meson exchanges cannot be ignored.     

Discovery potentials of double-charm tetraquarks

In this study, we investigate the discovery potential of double-charm tetraquarks \begin{document}$ T^{\{cc\}}_{[\bar{q}\bar{q}']} $\end{document}. We find that their production cross sections at the LHCb with \begin{document}$ \sqrt{s} = 13 $\end{document} TeV reach \begin{document}$ \mathcal{O}(10^4) $\end{document} pb, which indicates that the LHCb has collected \begin{document}$ \mathcal{O}(10^8) $\end{document} such particles. Through the decay channels of \begin{document}$ T^{\{cc\}}_{[\bar{u}\bar{d}]}\to D^{+}K^{-}\pi^{+} $\end{document} or \begin{document}$ D^0D^+\gamma $\end{document} (if stable) or \begin{document}$ T^{\{cc\}}_{[\bar{u}\bar{d}]}\to D^0D^{*+}\to D^0D^0\pi^+ $\end{document} (if unstable), it is highly expected that they get discovered at the LHCb in the near future. We also discuss the productions and decays of the double-charm tetraquarks at future Tera-\begin{document}$ Z $\end{document} factories.     

Hadronic coupling constants of the lowest hidden-charm pentaquark state,using QCD sum rules with rigorous quark-hadron duality

In this article,we illustrate how to calculate the hadronic coupling constants of the pentaquark states with QCD sum rules based on rigorous quark-hadron quality.We then study the hadronic coupling constants of the lowest diquark-diquark-antiquark type hidden-charm pentaquark state with spin-parity J~P=1/2~-in detail,and calculate the partial decay widths.The total width Г(P_c)=14.32±3.31 MeV is compatible with the experimental value T_(P_c(4312))=9.8±2.7_(-4.5)~(+3.7) MeV from the LHCb collaboration and favors assigning the P_c(4312) to be the [ud][uc]c pentaquark state with J~P=1/2~-.The hadronic coupling constants have the relation ■,and favor the hadronic dressing mechanism.The P_c(4312) may have a diquark-diquark-antiquark type pentaquark core with the typical size of the qqq-type baryon states.The strong couplings to the meson-baryon pairs DE_c lead to some pentaquark molecule components,and the P_c(4312) may spend a rather large time as the DE_c molecular state.     

Nucleon effective mass splitting and density-dependent symmetry energy effects on elliptic flow in heavy ion collisions at E_(lab)= 0.09 ～ 1.5 GeV/nucleon

By incorporating an iso spin-dependent form of the momentum-dependent potential in the ultra-relativistic quantum molecular dynamics(UrQMD) model,we systematically investigate effects of the neutron-proton effective mass splitting m*_(n-p)=m*_n-m*_p/m and the density-dependent nuclear symmetry energy E_(sym)(ρ) on the elliptic flow v_2 in~(197)Au+~(197) Au collisions at beam energies from 0.09 to 1.5 GeV/nucleon.It is found that at higher beam energies(≥0.25 GeV/nucleon) with the approximately 75 MeV difference in slopes of the two different E_(sym)(ρ),and the variation of m*_(n-p) ranging from-0.03 to 0.03 at saturation density with isospin asymmetry δ=(ρ_n-ρ_p)/ρ-0.2,the E_(sym)(ρ) has a stronger influence on the difference in v_2 between neutrons and protons,i.e.,v_2~n-v_2~p,than m*_(n-p) has.Meanwhile,at lower beam energies(≤0.25 GeV/nucleon),v_2~n-v_2~p is sensitive to both the E_(sym)(ρ) and the m*_(n-p).Moreover,the influence of m*_(n-p) on v_2~n-v_2~p is more evident with the parameters of this study when using the soft,rather than stiff,symmetry energy.     

Simple Woods-Saxon-type form for Ωα and Ξα interactions using folding model

We derive a simple Woods-Saxon-type form for potentials between Y=Ξ,Ωandαusing a single-folding potential method,based on a separable Y-nucleon Potential.The PotentialsΞ+αandΩ+αare accordingly obtained using the ESC08 c Nijmegens potential(in 3 S1 channel)and HAL QCD collaborationΩN interactions(in lattice QCD),respectively.In deriving the potential between Y andα,the same potential between Y and N is employed.The binding energy,scattering length,and effective range of the Y particle on the alpha particle are approximated by the resulting potentials.The depths of the potentials inΩαandΞαsystems are obtained at-61 MeV and-24.4 MeV,respectively.In the case of theΞαpotential,a fairly good agreement is observed between the single-folding potential method and the phenomenological potential of the Dover-Gal model.These potentials can be used in 3-,4-and 5-body cluster structures ofΩandΞhypernuclei.     

Revisiting D-meson twist-2, 3 distribution amplitudes

Owing to the significant difference between the experimental measurements and the theoretical predictions of the standard model (SM) for the value of \begin{document}$ {\cal{R}}(D) $\end{document} of the semileptonic decay \begin{document}$ B\to D\ell\bar{\nu}_{\ell} $\end{document}, researchers speculate that this decay may be evidence of new physics beyond the SM. Usually, the D-meson twist-2, 3 distribution amplitudes (DAs) \begin{document}$ \phi_{2;D}(x,\mu) $\end{document}, \begin{document}$ \phi_{3;D}^p(x,\mu) $\end{document} , and \begin{document}$ \phi_{3;D}^\sigma(x,\mu) $\end{document} are the main error sources when perturbative QCD factorization and light-cone QCD sum rules are used to study \begin{document}$ B\to D\ell\bar{\nu}_{\ell} $\end{document}. Therefore, it is important to obtain more reasonable and accurate behaviors for these DAs. Motivated by our previous work [Phys. Rev. D 104, no.1, 016021 (2021)] on pionic leading-twist DA, we revisit D-meson twist-2, 3 DAs \begin{document}$ \phi_{2;D}(x,\mu) $\end{document}, \begin{document}$ \phi_{3;D}^p(x,\mu) $\end{document}, and \begin{document}$ \phi_{3;D}^\sigma(x,\mu) $\end{document}. New sum rule formulae for the \begin{document}$\xi $\end{document}-moments of these three DAs are suggested for obtaining more accurate values. The light-cone harmonic oscillator models for the DAs are improved, and their parameters are determined by fitting the values of ξ-moments via the least squares method.     

Charmonia and bottomonia in asymmetric magnetized hot nuclear matter

We investigate the mass-shift of P-wave charmonium(χ_(c0),χ_(c1)),and S and P-wave bottomonium(η_b,■,χ_(b0),and χ_(b1)) states in magnetized hot asymmetric nuclear matter using the unification of QCD sum rules(QCDSR)and the chiral S U(3) model.Within QCDSR,we use two approaches,i.e.,the moment sum rule and the Borel sum rule.The magnetic field induced scalar gluon condensate α_s/πG_(μν)~aG~(aμν) and the twist-2 gluon operatorα_s/πG_(μσ)~aG~a ν~σcalculated in the chiral S U(3) model are utilised in QCD sum rules to calculate the in-medium mass-shift of the above mesons.The attractive mass-shift of these mesons is observed,which is more sensitive to magnetic field in the high density regime for charmonium,however less so for bottomonium.These results may be helpful to understand the decay of higher quarkonium states to the lower quarkonium states in asymmetric heavy ion collision experiments.     

Exact recession velocity and cosmic redshift based on cosmological principle and Yang-Mills gravity

Based on the cosmological principle and quantum Yang-Mills gravity in the super-macroscopic limit, we obtain an exact recession velocity and cosmic redshift z, as measured in an inertial frame F ≡ F(t, x, y, z). For a matter-dominated universe, we have the effective cosmic metric tensor G_(μν)(t) =(B~2(t),-A~2(t),-A~2(t),-A~2(t)),A ∝ B ∝ t~(1/2), where t has the operational meaning of time in F frame. We assume a cosmic action S ≡ S cos involving Gμν(t) and derive the ‘Okubo equation' of motion, G μν(t)?μS ?νS-m2= 0, for a distant galaxy with mass m. This cos-√mic equation predicts an exact recession velocity, ■, where H = A˙(t)/A(t) and Co = B/A, as observed in the inertial frame F. For small velocities, we have the usual Hubble's law r≈ rH for recession velocities. Following the formulation of the accelerated Wu-Doppler effect, we investigate cosmic redshifts z as measured in F. It is natural to assume the massless Okubo equation, G μν(t)?μψe?νψe= 0, for light emitted from accelerated distant galaxies. Based on the principle of limiting continuation of physical laws, we obtain a transformation for covariant wave 4-vectors between and inertial and an accelerated frame, and predict a relationship for the exact recession velocity and cosmic redshift, z = [(1 + V_r)/(1-V_r~2)~(1/2)]-1, where Vr= r˙/Co 1, as observed in the inertial frame F. These predictions of the cosmic model are consistent with experiments for small velocities and should be further tested.     

Gluon-pair-creation production model of strong interaction vertices

By studying the \begin{document}$\eta_c$\end{document} exclusive decay to double glueballs, we introduce a model to phenomenologically mimic the gluon-pair-vacuum interaction vertices, namely the \begin{document}$0^{++}$\end{document} model. Based on this model, we study glueball production in pseudoscalar quarkonium decays, explicitly \begin{document}$\eta_c \to f_0(1500)\eta(1405)$\end{document}, \begin{document}$\eta_b\to f_0(1500)\eta(1405)$\end{document} , and \begin{document}$\eta_b\to f_0(1710)\eta(1405)$\end{document} processes. Among them \begin{document}$f_0(1500)$\end{document} and \begin{document}$f_0(1710)$\end{document} are well-known scalars possessing large glue components, while \begin{document}$\eta(1405)$\end{document} is a potential candidate for a pseudoscalar glueball. The preliminary calculation results indicate that these processes are marginally accessible in the presently running experiments BES III, BELLE II, and LHCb.     

Nuclear mass parabola and its applications

We propose a method for extracting the properties of the isobaric mass parabola based on the total doubleβ-decay energies of isobaric nuclei.Two important parameters of the mass parabola,the location of the most β-stable nuclei Z_A and the curvature parameter b_A,are obtained for 251 A values,based on the total double β-decay energies of nuclei compiled in the AME2016 database.The advantage of this approach is that the pairing energy term P_A caused by the odd-even variation can be removed in the process,as well as the mass excess M(A,Z_A) of the most stable nuclide for the mass number A,which are employed in the mass parabolic fitting method.The Coulomb energy coefficient a_c=0.6910 MeV is determined by the mass difference relation for mirror nuclei,and the symmetry energy coefficient is also studied by the relation a_(sym)(A)=025 b_AZ_A.     

Electromagnetic field produced in high-energy small collision systems within charge density models of nucleons

Recent experiments show that \begin{document}$ \Delta\gamma $\end{document}, an observable designed to detect the chiral magnetic effect (CME), in small collision systems (\begin{document}$ p+A $\end{document}) is similar to that in heavy ion collisions (\begin{document}$ A+A $\end{document}). This introduces a challenge to the existence of the CME because it is believed that no azimuthal correlation exists between the orientation of the magnetic field (\begin{document}$ \Phi_B $\end{document}) and participant plane (\begin{document}$ \Phi_2 $\end{document}) in small collision systems. In this work, we introduce three charge density models to describe the inner charge distributions of protons and neutrons and calculate the electric and magnetic fields produced in small \begin{document}$ p+A $\end{document} collisions at both RHIC and LHC energies. Our results show that the contribution of the single projectile proton is the main contributor to the magnetic field after averaging over all participants. The azimuthal correlation between \begin{document}$ \Phi_B $\end{document} and \begin{document}$ \Phi_2 $\end{document} is small but not vanished. Additionally, owing to the large fluctuation in field strength, the magnetic-field contribution to \begin{document}$ \Delta\gamma $\end{document} may be large.     

Neutron skin thickness of 90Zr and symmetry energy constrained by charge exchange spin-dipole excitations

The charge exchange spin-dipole (SD) excitations of \begin{document}$ ^{90} $\end{document}Zr are studied using the Skyrme Hartee-Fock plus proton-neutron random phase approximation with SAMi-J interactions. The experimental value of the model-independent sum rule obtained from the SD strength distributions of \begin{document}$ ^{90} $\end{document}Zr(p, n)\begin{document}$ ^{90} $\end{document}Nb and \begin{document}$ ^{90} $\end{document}Zr(n, p)\begin{document}$ ^{90} $\end{document}Y is used to deduce the neutron skin thickness. The neutron skin thickness \begin{document}$ \Delta r_{np} $\end{document} of \begin{document}$ ^{90} $\end{document}Zr is extracted as \begin{document}$ 0.083\pm0.032 $\end{document} fm, which is similar to the results of other studies. Based on the correlation analysis of the neutron skin thickness \begin{document}$ \Delta r_{np} $\end{document} and the nuclear symmetry energy J as well as its slope parameter L, a constraint from the extracted \begin{document}$ \Delta r_{np} $\end{document} leads to the limitation of J to \begin{document}$ 29.2 \pm 2.6 $\end{document} MeV and L to \begin{document}$ 53.3 \pm 28.2 $\end{document} MeV.     

Phase structures of neutral dense quark matter and applicationto strange stars

In the contact interaction model,the quark propagator has only one solution,namely,the chiral symmetry breaking solution,at vanishing temperature and density in the case of physical quark mass.We generalize the condensate feedback onto the coupling strength from the 2 flavor case to the 2+1 flavor case,and find the Wigner solution appears in some regions,which enables us to tackle chiral phase transition as two-phase coexistences.At finite chemical potential,we analyze the chiral phase transition in the conditions of electric charge neutrality andβequilibrium.The four chemical potentials,μ_u,μ_d,μ_s and He,are constrained by three conditions,so that one inde-pendent variable remains:we choose the average quark chemical potential as the free variable.All quark masses and number densities suffer discontinuities at the phase transition point.The strange quarks appear after the phase trans-ition since the system needg more energy to produce a d.-quark than an s-quark.Taking the EOS as an input,the TOV equations are solved numerically,and we show that the mass--radius relation is sensitive to the EOS.The max-imum mass of strange quark stars is not susceptible to the parameter Aq we introduced.     

Sequential single pion production explaining the dibaryon "d*(2380)" peak

In this study, we investigate the two step sequential one pion production mechanism, that is, \begin{document}$ np(I=0)\to $\end{document}\begin{document}$ \pi^-pp $\end{document} followed by the fusion reaction \begin{document}$ pp\to \pi^+d $\end{document}, to describe the \begin{document}$ np\to \pi^+\pi^-d $\end{document} reaction with \begin{document}$ \pi^+\pi^- $\end{document} in state \begin{document}$ I=0 $\end{document}. In this reaction, a narrow peak identified with a "\begin{document}$ d(2380) $\end{document}" dibaryon has been previously observed. We discover that the second reaction step \begin{document}$ pp\to \pi^+d $\end{document} is driven by a triangle singularity that determines the position of the peak of the reaction and the high strength of the cross section. The combined cross section of these two mechanisms produces a narrow peak with a position, width, and strength, that are compatible with experimental observations within the applied approximations made. This novel interpretation of the peak accomplished without invoking a dibaryon explains why this peak has remained undetected in other reactions.