Source author record

Su Houng Lee

Su Houng Lee 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

56works
6topics
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

56 published item(s)

preprint2026arXiv

Exploring bound states and interactions of the nucleon-antinucleon system in a constituent quark model

In this work, we study the nucleon-antinucleon system in a constituent quark model. We first construct the nucleon-antinucleon wave function such that the multiquark and multiantiquark components each satisfy the Pauli exclusion principle, and then investigate the possibility of a bound state using a Hamiltonian that includes color-color, spin-dependent interactions and three-quark potentials. Our results indicate that, in specific channels, the nucleon-antinucleon interaction exhibits significant attraction, suggesting a strong possibility of a bound state.

preprint2026arXiv

Role of $Σ^*(1385)$ on $Λ$ hyperon polarization in relativistic heavy ion collisions

The effect of $Σ^*(1385)$ baryon resonance on the time evolution of the $Λ$ hyperon polarization in hadronic matter is studied using a kinetic approach. This approach explicitly includes the production of the $Σ^*$ resonance from the $Λ-π$ and $Σ(1192)-π$ scatterings as well as its decay into $Λ+π$ or $Σ+π$. The resulting coupled kinetic equations governing the time evolution of $Λ$, $Σ$ and $Σ^*$ numbers and polarizations are solved for Au-Au collisions at $\sqrt{s_{NN}}=7.7$ GeV and 20-50\% centrality, using initial values determined by thermal yields and the thermal vorticity at chemical freeze-out temperature. As the hadronic matter expands and cools, the $Λ$ polarization is found to increase slightly during early times and then decreases very slowly afterwards, while the $Σ$ polarization remains nearly constant and the $Σ^*$ polarization continuously decreases. Including feed-down contributions to the $Λ$ polarization from the decays of partially polarized $Σ^0$, $Σ^*$, and $Ξ(1322)$ hyperons, where the $Ξ$ polarization is obtained by solving coupled kinetic equations for the $Ξ$ and $Ξ^*(1532)$ system, the resulting $Λ$ polarization becomes smaller and decreases over time. In both cases, however, the time variation of the $Λ$ polarization is sufficiently small to support the assumption of an early freeze-out of $Λ$ spin degree of freedom in relativistic heavy ion collisions.

preprint2022arXiv

$ϕ$ meson properties in nuclear matter from QCD sum rules with chirally separated four-quark condensates

The modification of the $ϕ$ meson spectrum in nuclear matter is studied in an updated QCD sum rule analysis, taking into account recent improvements in properly treating the chiral invariant and breaking components of four-quark condensates. Allowing both mass and decay width to change at finite density, the QCD sum rule analysis determines certain combinations of changes for these parameters that satisfy the sum rules equally well. A comprehensive error analysis, including uncertainties related to the behavior of various condensates at linear order in density, the employed renormalization scale and perturbative corrections of the Wilson coefficients, is used to compute the allowed ranges of these parameter combinations. We find that the $ϕ$ meson mass shift in nuclear matter is especially sensitive to the strange sigma term $σ_{sN}$, which determines the decrease of the strange quark condensate in nuclear matter. Specifically, we obtain a linear relation between the width $Γ_ϕ$ and mass shift $Δm_ϕ$ given as $ Γ_ϕ = aΔm_ϕ + bσ_{sN}+c$ with $a = (3.947^{+0.139}_{-0.130})$, $b = (0.936^{+0.180}_{-0.177} )$ and $c = -(7.707^{+4.791}_{-5.679}) \mathrm{MeV}$.

preprint2022arXiv

Algebraic approach to quarkyoniclike configuration and stable diquarks in dense matter

We study the color-spin interaction energy of a quark, a diquark and a baryon with their surrounding baryons and/or quark matter. This is accomplished by classifying all possible flavor and spin states of the resulting multiquark configuration in both the flavor SU(2) and SU(3) symmetric cases. We find that while the baryon has the lowest interaction energy when there is only a single surrounding baryon, the quark has the lowest interaction energy when the surrounding has more than three baryons or becomes a quark gas. As the short range nucleon-nucleon interactions are dominated by the color-spin interactions, our finding suggests that the baryon modes near other baryons are suppressed due to larger repulsive energy compared to that of a quark and thus provides a quark model basis for the quarkyoniclike phase in dense matter. At the same time, when the internal interactions are taken into account, and the matter density is high so that the color-spin interaction becomes the dominant interaction, the diquark becomes the lowest energy configuration and will thus appear in both the dense baryonic and/or quark matter.

preprint2022arXiv

Production of $P_{c}$(4312) state in electron-proton collisions

We study the cross sections for the electro-production of $P_c(4312)$ particle, a recently discovered pentaquark state, in electron-proton collisions assuming possible quantum numbers to be $J^{P}=\frac{1}{2}^\pm, \frac{3}{2}^\pm$. $\sqrt{s}$ is set to the energy of the future Electron Ion Collider at Brookhaven National Laboratory, in order to asses the possibility of the measurement in this facility. One can discriminate the spin of $P_c(4312)$ by comparing the pseudorapidity distribution in two different polarization configurations for proton and electron beams. Furthermore, the parity of $P_c(4312)$ can be discerned by analyzing the decay angle in the $P_c \rightarrow p +J/ψ$ channel. As the multiplicity of $P_c$ production in our calculation is large, the EIC can be considered as a future facility for precision measurement of heavy pentaquarks.

preprint2022arXiv

Transport coefficients of heavy quarkonia comparing with heavy quark coefficients

We revisit the transport coefficients of heavy quarkonia moving in high-temperature QCD plasmas. The thermal width and mass shift for heavy quarkonia are closely related to the momentum diffusion coefficient and its dispersive counterpart for heavy quarks, respectively. For quarkonium at rest in plasmas the longitudinal gluon part of the color-singlet self-energy diagram is sufficient to determine the leading-order thermal width, whereas the momentum dependence is obtained from the transverse gluon channel. Using the quarkonium-gluon effective vertex based on the dipole interaction of color charges, we discuss the damping rate, the effective rest and kinetic mass shifts of slowly moving quarkonia and compare with the corresponding coefficients of heavy quarks.

preprint2021arXiv

The masses of hadrons in the chiral symmetry restored vacuum

We calculate the masses of the vector and axial-vector mesons as well as the nucleon and the delta resonance in the chiral symmetry restored vacuum. This is accomplished by separating the quark operators appearing in the QCD sum rules for these hadrons into the chiral symmetric and symmetry breaking parts depending on the contributions of the fermion zero modes. We then extract the vacuum expectation values of all the separated parts of the quark operators using the QCD sum rule relations for these hadrons with their vacuum masses and widths. By taking the chiral symmetry breaking parts to be zero while keeping the symmetric operators to their vacuum values, we obtain the chiral symmetric part of the hadron masses. We find that the masses of chiral partners, such as the $(ρ,a_1)$ and $(K^*,K_1)$, become degenerate to values between 500 and 600 MeV in the chiral symmetry restored vacuum, while parity partners $(ω,f_1)$ that are chiral partners only in the limit where the disconnected diagrams are neglected remain non-degenerate with masses $(655,1060)$ MeV, respectively. The masses of the nucleon and the Delta are also found to reduce to about 500 and 600 MeV, respectively, in the chiral symmetric vacuum. This shows that while chiral symmetry breaking is responsible for the mass difference between chiral partner, both the meson and baryon retain non-trivial fraction of their masses in the chiral symmetry restored vacuum.

preprint2020arXiv

$J/ψ$ near $T_c$

We calculate the mass shift and thermal decay width of the $J/ψ$ near the QCD transition temperature $T_c$ by imposing two independent constraints on these variables that can be obtained first by solving the Schrödinger equation and second from the QCD sum rule approach. While the real part of the potential is determined by comparing the QCD sum rule result for charmonium and the D meson to that from the potential model result, the imaginary potential is taken to be proportional to the perturbative form multiplied by a constant factor, which in turn can be determined by applying the two independent constraints. The result shows that the binding energy and the thermal width becomes similar in magnitude at around $T=1.09T_c$, above which the sum rule analysis also becomes unstable, strongly suggesting that the $J/ψ$ will melt slightly above $T_c$.

preprint2020arXiv

Charmed hadron production in an improved quark coalescence model

We study the production of charmed hadrons $D^{0}$ and $Λ_c^+$ in relativistic heavy-ion collisions using an improved quark coalescence model. In particular, we extend the usual coalescence model by letting a produced hadron to have the same velocity as the center-of-mass velocity of coalesced constituent quarks during hadronization to take into account the effect of collective flow in produced quark-gluon plasma. This results in a shift of charmed resonances of higher masses to larger transverse momenta ($p_T^{}$). Requiring all charm quarks of very low $p_T^{}$ to be converted to hadrons via coalescence and letting charm quarks not undergoing coalescence to hadronize by independent fragmentation, we obtain a good description of the measured yield ratio $Λ_c^+/D^0$ as a function of $p_T^{}$ in $\text{Au} + \text{Au}$ collisions at $\sqrt{s_{NN}}^{}=200$~GeV by the STAR Collaboration at the Relativistic Heavy Ion Collider.

preprint2020arXiv

D meson mass and heavy quark potential at finite temperature

Based on the observation that the heavy quark-antiquark potential value at infinity corresponds to twice the D meson mass, we constrain the asymptotic value of the heavy quark potential in a hot medium through a QCD sum rule calculation of the D meson at finite temperature. We find that to correctly reproduce the QCD sum rule results as well as a recent model calculation for the D meson mass near the critical temperature, the heavy quark potential should be composed mostly of the free energy with an addition of a small but nontrivial fraction of the internal energy. Combined with a previous study comparing potential model results for the $J/ψ$ to a QCD sum rule calculation, we conclude that the composition of the effective heavy quark potential should depend on the inter-quark distance. Namely, the potential is dominated by the free energy at short distance, while at larger separation, it has a fraction of about 20% of internal energy.

preprint2020arXiv

Signatures of the vortical quark-gluon plasma in hadron yields

We investigate the hadron production from the vortical quark-gluon plasma created in heavy-ion collisions. Based on the quark-coalescence and statistical hadronization models, we show that total hadron yields summed over the spin components are enhanced by the local vorticity with quadratic dependence. The enhancement factor amounts to be a few percent and may be detectable within current experimental sensitivities. We also show that the effect is stronger for hadrons with larger spin, and thus propose a new signature of the local vorticity, which may be detected by the yield ratio of distinct hadron species having different spins such as $ϕ$ and $η'$. The vorticity dependence of hadron yields seems robust, with consistent predictions in both of the hadron production mechanisms for reasonable values of the vorticity strength estimated for heavy-ion collisions.

preprint2019arXiv

Baryon-baryon interactions at short distances -- constituent quark model meets lattice QCD

The interaction energies between two baryons at short distance in different flavor channels are calculated from the constituent quark model (CQM) and are compared with the recent lattice QCD (LQCD) results for baryon-baryon potentials at short distance. We consider the six-quark system with two strange quarks and focus on the quantum numbers, (Flavor,Spin)=(1,0),(8,1),(10,1),($\overline{10}$,1) and (27,0). The interaction energy is defined by subtracting out isolated baryon masses and relative kinetic energy of two baryons from the total energy of a compact six-quark state. We introduce interaction energy ratio between different flavors as a useful measure to test the prediction of CQM. We find that the ratios in CQM show good agreement with those in LQCD, which indicates that the short range part of the baryon-baryon interaction can be understood qualitatively in terms of the Pauli principle and spin-dependent color interaction among constituent quarks.

preprint2019arXiv

Production of multi-charmed hadrons by recombination in heavy ion collisions

We study the production of multi-charmed hadrons by recombination in heavy ion collisions by focusing on the production of $Ξ_{cc}$, $Ξ_{cc}^*$, $Ω_{scc}$, $Ω_{scc}^*$, $Ω_{ccc}$ baryons and X(3872) mesons. Starting from the estimation of yields for those hadrons at chemical freeze-out in both the statistical and coalescence model, we evaluate their transverse momentum distributions at mid-rapidity in the coalescence model. We show that yields of multi-charmed hadrons in heavy ion collisions at RHIC and LHC are large enough, and thereby not only multi-charmed hadrons observed so far, e.g., the $Ξ_{cc}$ but also those which have not been observed yet, can be discovered sufficiently in heavy ion collisions. We also find that the transverse momentum distribution ratio between various multi-charmed hadrons sensitively reflects the interplay between quark contents of corresponding hadrons as well as the transverse momentum distribution of charm quarks at the hadronization point, and therefore we insist that studying both the transverse momentum distributions of multi-charmed hadrons themselves and transverse momentum distribution ratios between various multi-charmed hadrons provide us with useful information on hadron production mechanism involving charm quarks in heavy ion collisions.

preprint2019arXiv

Theory on Hadrons in Nuclear Medium

After decades-long attempts to measure the mass shift and understand the origin of hadron mass, it became clear that one has to analyze hadrons with small vacuum width. Also, to identify the effect of chiral symmetry breaking, one has to start by looking at chiral partners. In this talk, I will review why such consideration inevitably led us to consider $K^*$ and $K_1$ in nuclear matter [T. Song, T. Hatsuda, S H Lee, PLB792 (2019) 160-169]. With the kaon beam at JPARC, one could observe the mass shift of both particles in a nuclear target experiment. Once the masses and mass difference of $K^*$ and $K_1$ mesons are measured, we will be closer to understanding the origin of the hadron masses and the effects of chiral symmetry breaking in them.

preprint2016arXiv

D mesons in a magnetic field

We investigate the mass spectra of open heavy flavor mesons in an external constant magnetic field within QCD sum rules. Spectral ansätze on the phenomenological side are proposed in order to properly take into account mixing effects between the pseudoscalar and vector channels, and the Landau levels of charged mesons. The operator product expansion is implemented up to dimension-5 operators. As a result, we find for neutral D mesons a significant positive mass shift that goes beyond simple mixing effects. In contrast, charged D mesons are further subject to Landau level effects, which together with the mixing effects almost completely saturate the mass shifts obtained in our sum rule analysis.

preprint2016arXiv

Dibaryons with two strange quarks and one heavy flavor in a constituent quark model

We investigate the symmetry property and the stability of dibaryons containing two strange quarks and one heavy flavor with $I=\frac{1}{2}$. We construct the wave function of the dibaryon in two ways. First, we directly construct the color and spin state of the dibaryon starting from the four possible SU(3) flavor state. Second, we consider the states composed of five light quarks, and then construct the wave function of the dibaryon by adding one heavy quark. The stability of the dibaryon against the strong decay into two baryons is discussed by using variational method in a constituent quark model with confining and hyperfine potential. We find that for all configurations with S=0,1,2, the ground states of the dibaryons are the sum of two baryons, and there are no compact bound state that is stable against the strong decay.

preprint2016arXiv

Dibaryons with two strange quarks and total spin zero in a constituent quark model

We investigate the symmetry property and construct the wave function of the dibaryon states containing two strange quarks with S=0 in both the flavor SU(3) symmetric and breaking cases. We discuss how the color $\otimes$ isospin $\otimes$ spin states of dibaryon in the symmetry broking case of flavor SU(3) can be extracted from the fully antisymmetric states in flavor SU(3). The stability of the dibaryon against the strong decay into two baryons are then discussed, by using the variational method within a constituent quark model with a confining and color-spin interactions. To compare our results with that from lattice QCD in flavor SU(3) limit, we search for the stable H-dibaryon in a wide range of $π$ meson mass. We find that with the given potential, there is no compact six quark dibaryon state in the SU(3) flavor symmetry broken case with realistic quark masses as well as in flavor SU(3) symmetric case in a wide range of quark masses.

preprint2016arXiv

QCD sum rules for the neutron, $Σ$, and $Λ$ in neutron matter

The nuclear density dependencies of the neutron and $Σ$ and $Λ$ hyperons are important inputs in the determination of the neutron star mass as the appearance of hyperons coming from strong attractions significantly changes the stiffness of the equation of state (EOS) at iso-spin asymmetric dense nuclear matter. In-medium spectral sum rules have been analyzed for the nucleon, $Σ$, and $Λ$ hyperon to investigate their properties up to slightly above the saturation nuclear matter density by using the linear density approximation for the condensates. The construction scheme of the interpolating fields without derivatives has been reviewed and used to construct a general interpolating field for each baryon with parameters specifying the strength of independent interpolating fields. Optimal choices for the interpolating fields were obtained by requiring the sum rules to be stable against variations of the parameters and the result to be consistent with known phenomenology. The optimized result shows that Ioffe's choice is not suitable for the $Λ$ hyperon sum rules. It is found that, for the $Λ$ hyperon interpolating field, the up and down quark combined into the scalar diquark structure $u^T C γ_5 d$ should be emphasized to ensure stable sum rules. The quasi-$Σ$ and -$Λ$ hyperon energies are always found to be higher than the quasineutron energy in the region $0.5 <ρ/ρ_0<1.5 $ where the linear density approximation in the sum-rule analysis is expected to be reliable.

preprint2016arXiv

The mass of heavy-light mesons in a constituent quark picture with partially restored chiral symmetry

We probe effects of the partial chiral symmetry restoration to the mass of heavy-light mesons in a constituent quark model by changing the constituent quark mass of the light quark. Due to the competing effect between the quark mass and the linearly rising potential, whose contribution to the energy increases as the quark mass decreases, the heavy-light meson mass has a minimum value near the constituent quark mass typically used in the vacuum. Hence, the meson mass increases as one decreases the constituent quark mass consistent with recent QCD sum rule analyses, which show an increasing $D$ meson mass as the chiral order parameter decreases.

preprint2015arXiv

Charmonium Spectroscopy in Strong Magnetic Fields by QCD Sum Rules: S-Wave Ground States

We investigate quarkonium mass spectra in external constant magnetic fields by using QCD sum rules. We first discuss a general framework of QCD sum rules necessary for properly extracting meson spectra from current correlators computed in the presence of strong magnetic fields, that is, a consistent treatment of mixing effects caused in the mesonic degrees of freedom. We then implement operator product expansions for pseudoscalar and vector heavy-quark current correlators by taking into account external constant magnetic fields as operators, and obtain mass shifts of the lowest-lying bound states $η_c$ and $J/ψ$ in the static limit with their vanishing spatial momenta. Comparing results from QCD sum rules with those from hadronic effective theories, we find that the dominant origin of mass shifts comes from a mixing between $η_c$ and $J/ψ$ with a longitudinal spin polarization, accompanied by other subdominant effects such as mixing with higher excited states and continua.

preprint2015arXiv

Dibaryons in a constituent quark model

We investigate the properties of dibaryons containing u and d quarks in the constituent quark model. In constructing the ground state wave function, we choose the spatial part to be fully symmetric and the remaining color, isospin and spin part to be antisymmetric so as to satisfy the the Pauli principle. By adapting the IS coupling scheme that combine the isospin basis function with the spin basis function, and subsequently coupling this to the color singlet basis function, we construct the color $\otimes$ isospin $\otimes$ spin states compatible with the physical states of the dibaryon. By using the variational method, we then calculate the mass of the dibaryon in a nonrelativistic potential model, involving Coulomb, color confinement and color-spin hyperfine interaction. In particular, to asses the stability for different types of the confinement potential, we introduce one that is linearly proportional to the interquark distance and another to its square root. For all cases considered, we find that there are no compact bound states against the strong decay.

preprint2015arXiv

New determination of $\mathcal{S} \mathcal{T} \langle N| \overline{q} D_μ D_ν q |N \rangle$ based on recent experimental constraints

The symmetric and traceless part of the matrix element $\mathcal{S} \mathcal{T} \langle N| \overline{q} D_μ D_ν q |N \rangle$ can be determined from the second moment of the twist-3 parton distribution function $e(x)$. Recently, novel experimental data on $e(x)$ have become available, which enables us to evaluate the magnitude of the above matrix element with considerably reduced systematic uncertainties. Based on the new experimental data, we show that $\mathcal{S} \mathcal{T} \langle N| \overline{q} D_μ D_ν q |N \rangle$ is likely to be at least an order of magnitude smaller than what previous model-based estimates have so far suggested. We discuss the consequences of this observation for the analysis of deep inelastic scattering and QCD sum rules studies at finite density for the vector meson and the nucleon, in which this matrix element is being used as an input parameter.

preprint2015arXiv

Quarkonium formation time in relativistic heavy-ion collisions

We calculate the quarkonium formation time in relativistic heavy-ion collisions from the space-time correlator of heavy quark vector currents in a hydrodynamics background with the initial nonequilibrium stage expanding only in the longitudinal direction. Using in-medium quarkonia properties determined with the heavy quark potential taken to be the free energy from lattice calculations and the fact that quarkonia can only be formed below their dissociation temperatures due to color screening, we find that $Υ$(1S), $Υ$(2S), $Υ$(3S), $J/ψ$ and $ψ^\prime$ are formed, respectively, at 1.2, 6.6, 8.8, 5.8, and 11.0 fm/c after the quark pair are produced in central Au+Au collisions at the top energy of Relativistic Heavy Ion Collider (RHIC), and these times become shorter in semi-central collisions. We further show, as an example, that including the effect of formation time enhances appreciably the survivability of $Υ$(1S) in the produced hot dense matter.

preprint2015arXiv

Symmetry energy in cold dense matter

We calculate the symmetry energy in cold dense matter both in the normal quark phase and in the 2-color superconductor (2SC) phase. For the normal phase, the thermodynamic potential is calculated by using hard dense loop (HDL) resummation to leading order, where the dominant contribution comes from the longitudinal gluon rest mass. The effect of gluonic interaction to the symmetry energy, obtained from the thermodynamic potential, was found to be small. In the 2SC phase, the non-perturbative BCS paring gives enhanced symmetry energy as the gapped states are forced to be in the common Fermi sea reducing the number of available quarks that can contribute to the asymmetry. We used high density effective field theory to estimate the contribution of gluon interaction to the symmetry energy. Among the gluon rest masses in 2SC phase, only the Meissner mass has iso-spin dependence although the magnitude is much smaller than the Debye mass. As the iso-spin dependence of gluon rest masses is even smaller than the case in the normal phase, we expect that the contribution of gluonic interaction to the symmetry energy in the 2SC phase will be minimal. The different value of symmetry energy in each phase will lead to different prediction for the particle yields in heavy ion collision experiment.

preprint2015arXiv

Temperature dependence of dimension 6 gluon operators and their effects on Charmonium

Starting from an earlier representation of the independent dimension 6 gluon operators in terms of color electric and magnetic fields, we estimate their changes near the critical temperature $T_c$ using the temperature dependence of the dimension 4 electric and magnetic condensates extracted from pure gauge theory on the lattice. We then improve the previous QCD sum rules for the $J/ψ$ mass near $T_c$ based on dimension 4 operators, by including the contribution of the dimension 6 operators to the OPE. We find an enhanced stability in the sum rule and confirm that the $J/ψ$ will undergo an abrupt change in the property across $T_c$.

preprint2014arXiv

QCD sum rules for magnetically induced mixing between $η_c$ and $J/ψ$

We investigate the properties of charmonia in strong magnetic fields by using QCD sum rules. We show how to implement the mixing effects between $η_c$ and $J/ψ$ on the basis of field-theoretical approaches, and then show that the sum rules are saturated by the mixing effects with phenomenologically determined parameters. Consequently, we find that the mixing effects are the dominant contribution to the mass shifts of the charmonia in strong magnetic fields.

preprint2013arXiv

Color spin wave functions of heavy tetraquark states

Using the variational method, we calculate the mass of the J^P=1^+ (ud)bar(bb) tetraquark containing two identical heavy antiquarks in a nonrelativistic potential model with color confinement and spin hyperfine interaction. In particular, we extend a previous investigation of the model by Brink and Stancu by investigating the effect of including the color anti-sextet component of the diquark configuration as well as using several more Gaussian parametrization for the L=0 part of the spatial wave function. We find that for the heavy tetraquark, the 6bar{6} component among the color singlet bases is negligible and that the previously used specific Gaussian spatial configuration is good enough in obtaining the ground state energy.

preprint2013arXiv

Free energy versus internal energy potential for heavy quark systems at finite temperature

Using the QCD sum rule with its operator product expansion reliably estimated from lattice calculations for the pressure and energy density of the quark-gluon plasma, we calculate the strength of the wave function of J/psi at the origin and find that it decreases with temperature above the transition temperature. This result is shown to follow exactly that obtained from the solution of the Schrödingier equation for a charm and anticharm quark pair using the free energy from lattice calculations as the potential and is in sharp contrast to that using the deeper potential associated with the internal energy, which shows an enhanced strength of the J/psi wave function at the origin. Our result thus has resolved the long-standing question of whether the free energy potential or the internal energy potential should be used in analyzing the spectrum of heavy quark systems at finite temperature.

preprint2013arXiv

Hadronic effects on the X(3872) meson abundance in heavy ion collisions

We study the hadronic effects on the X(3872) meson abundance in heavy ion collisions. We evaluate the absorption cross sections of the X(3872) meson by pions and rho mesons in the hadronic stage of heavy ion collisions, and further investigate the variation in the X(3872) meson abundance during the expansion of the hadronic matter for its two possible quantum number states; $J^P=1^+$ and $2^-$. We show that the absorption cross sections and the time evolution of the X(3872) meson abundance are strongly dependent on the structure and quantum number of the X(3872) meson. We thus suggest that studying the abundance of the X(3872) meson in relativistic heavy ion collisions provides a chance to infer its quantum number as well as its structure.

preprint2013arXiv

Nuclear Symmetry Energy from QCD sum rules

We calculate the nucleon self-energies in isospin-asymmetric nuclear matter using QCD sum rules. Taking the difference of these for the neutron and proton enables us to express the potential part of the nuclear symmetry energy in terms of local operators. We find that the scalar (vector) self-energy part gives a negative (positive) contribution to the nuclear symmetry energy which is consistent with the results from relativistic mean-field theories. Moreover, we find that an important contribution to the negative contribution of the scalar self-energy comes from the twist-4 matrix elements, whose leading density dependence can be extracted from deep inelastic scattering experiments. This suggests that the twist-4 contribution partly mimics the exchange of the $δ$ meson and that it constitutes an essential part in the origin of the nuclear symmetry energy from QCD. Our result also extends an early success of QCD sum rule method in understanding the symmetric nuclear matter in terms of QCD variables to the asymmetric nuclear matter case.

preprint2012arXiv

Charmonium mass in hot and dense hadronic matter

We investigate mass shifts of charmonia driven by change of the gluon condensate below but near transition temperatures at finite baryonic chemical potential. Extending previous prescription on the relation between gluon condensates and thermodynamic quantities, we model the gluon condensates of hadronic matter at finite temperature and baryonic chemical potential such that the scalar gluon condensate fits with the latest lattice QCD data. By making use of the QCD sum rule and the second order Stark effect, we find that the smoother transition in the full QCD can lead to moderate mass shifts of charmonia even below the transition temperature. We also find larger mass shift at fixed temperature as chemical potential increases. Existing data on charmonium-charmonium ratio is found to be consistent with the statistical hadronization scenario including the obtained mass shift.

preprint2012arXiv

Contributions of hyperon-hyperon scattering to subthreshold cascade production in heavy ion collisions

Using a gauged flavor SU(3)-invariant hadronic Lagrangian, we calculate the cross sections for the strangeness-exchange reactions YY to NΞ(Y=Λ, Σ) in the Born approximation. These cross sections are then used in the Relativistic Vlasov-Uehling-Uhlenbeck (RVUU) transport model to study Ξproduction in Ar+KCl collisions at incident energy of 1.76A GeV and impact parameter b=3.5 fm. We find that including the contributions of hyperon-hyperon scattering channels strongly enhances the yield of Ξ, leading to the abundance ratio Ξ^{-}/(Λ+Σ^{0})=3.38E-3, which is essentially consistent with the recently measured value of $(5.6 \pm 1.2_{-1.7}^{+1.8})\times 10^{-3}$ by the HADES collaboration at GSI.

preprint2011arXiv

A shear spectral sum rule in a non-conformal gravity dual

A sum rule which relates a stress-energy tensor correlator to thermodynamic functions is examined within the context of a simple non-conformal gravity dual. Such a sum rule was previously derived using AdS/CFT for conformal $\mathcal{N} = 4$ Supersymmetric Yang-Mills theory, but we show that it does not generalize to the non-conformal theory under consideration. We provide a generalized sum rule and numerically verify its validity. A useful byproduct of the calculation is the computation of the spectral density in a strongly coupled non-conformal theory. Qualitative features of the spectral densities and implications for lattice measurements of transport coefficients are discussed.

preprint2011arXiv

Diquarks: a QCD sum rule perspective

We propose a phenomenological QCD sum rule with an explicit diquark field to investigate the essential ingredients inside the hadrons. Introducing the mass (m_ϕ) and the condensate (ϕ^2) for the diquark field as parameters in the model, we find that the sum rule works well for Lambda, Lambda_c and Lambda_b. This implies that these Lambda baryons can be represented by a diquark and a quark configuration. We also find that there is a duality relation among the parameters (m_ϕ, ϕ^2), for which the sum rule is equally good. In the limit when ϕ^2 =0 we find m_ϕ=0.4 GeV, which can be thought as the constituent diquark mass.

preprint2011arXiv

Exotics from Heavy Ion Collisions

Discriminating hadronic molecular and multi-quark states is a long standing problem in hadronic physics. We propose here to utilize relativistic heavy ion collisions to resolve this problem, as exotic hadron yields are expected to be strongly affected by their structures. Using the coalescence model, we find that the exotic hadron yield relative to the statistical model result is typically an order of magnitude smaller for a compact multi-quark state, and larger by a factor of two or more for a loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic states could be produced and realistically measured at RHIC and LHC.

preprint2011arXiv

Multi-quark hadrons from Heavy Ion Collisions

Identifying hadronic molecular states and/or hadrons with multi-quark components either with or without exotic quantum numbers is a long standing challenge in hadronic physics. We suggest that studying the production of these hadrons in relativistic heavy ion collisions offer a promising resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their structures. Using the coalescence model for hadron production, we find that compared to the case of a non-exotic hadron with normal quark numbers, the yield of an exotic hadron is typically an order of magnitude smaller when it is a compact multi-quark state and a factor of two or more larger when it is a loosely bound hadronic molecule. We further find that due to the appreciable numbers of charm and bottom quarks produced in heavy ion collisions at RHIC and even larger numbers expected at LHC, some of the newly proposed heavy exotic states could be produced and realistically measured in these experiments.

preprint2011arXiv

Studying Exotic Hadrons in Heavy Ion Collisions

We investigate the possibilities of using measurements in present and future experiments on heavy ion collisions to answer some longstanding problems in hadronic physics, namely identifying hadronic molecular states and exotic hadrons with multiquark components. The yields of a selected set of exotic hadron candidates in relativistic heavy ion collisions are discussed in the coalescence model in comparison with the statistical model. We find that the yield of a hadron is typically an order of magnitude smaller when it is a compact multiquark state, compared to that of an excited hadronic state with normal quark numbers. We also find that some loosely bound hadronic molecules are formed more abundantly than the statistical model prediction by a factor of two or more. Moreover, due to the significant numbers of charm and bottom quarks produced at RHIC and even larger numbers expected at LHC, some of the proposed heavy exotic hadrons could be produced with sufficient abundance for detection, making it possible to study these new exotic hadrons in heavy ion collisions.

preprint2011arXiv

The flow of heavy flavor in hydrodynamics

The flow of charm is calculated in 2+1 ideal hydrodynamics by introducing the charge of $c\bar{c}$ pair assuming that the number of $c\bar{c}$ pairs is conserved in relativistic heavy-ion collisions. It is found that the mean radial flow velocity of charm quarks is smaller than that of bulk matter by 10$\sim$15 \% and the measured $v_2$ of heavy-flavor electrons is reproduced up to $p_T^e=$ 1.5 GeV/c in Au+Au collision at RHIC. The same flow is applied to regenerated $J/ψ$ and its $v_2$ is discussed.

preprint2010arXiv

$J/ψ$ production and elliptic flow in relativistic heavy-ion collisions

Using a two-component model for charmonium production, which includes contributions from both the initial hard nucleon-nucleon scattering and from the regeneration in the quark-gluon plasma, we study the nuclear modification factor $R_{AA}$ and elliptic flow $v_2$ of $J/ψ$ in relativistic heavy ion collisions. For the expansion dynamics of produced hot dense matter, we introduce a schematic fireball model with its transverse acceleration determined from the pressure gradient inside the fireball and azimuthally anisotropic expansion parameterized to reproduce measured $v_2$ of light hadrons. We assume that light hadrons freeze out at the temperature of 120 MeV while charmonia at 160 MeV, similar to the kinetic and chemical freeze-out temperatures in the statistical model, respectively. For the properties of charmonia in the quark-gluon plasma, we use the screening mass between their charm and anticharm quarks and their dissociation cross sections given by the perturbative QCD (pQCD) in the leading order and up to the next-to-leading order, respectively. For the relaxation time of charm and anticharm quarks in the quark-gluon plasma, we also use the one calculated in the leading order of pQCD. Modeling the effect of higher-order corrections in pQCD by introducing multiplicative factors to the dissociation cross section of charmonia and the elastic scattering cross sections of charm and anticharm quarks, we find that this effect is small for the $R_{AA}$ of $J/ψ$ as they suppress the number of initially produced $J/ψ$ but enhance the number of regenerated ones. The higher-order corrections increase, however, the $v_2$ of $J/ψ$. Our results suggest that the $v_2$ of $J/ψ$ can play an important role in discriminating between $J/ψ$ production from the initial hard collisions and from the regeneration in the quark-gluon plasma.

preprint2010arXiv

$R_{AA}$ of $J/ψ$ near mid-rapidity in heavy ion collisions at $\sqrt{s_{NN}}=200$ GeV

We build up a model to reproduce the experimentally measured $R_{AA}$ of $J/ψ$ near midrapidty in Au+Au collision at $\sqrt{s_{NN}}=200$ GeV. The model takes into account the $J/ψ$ suppression from the quark-gluon plasma and hadron gas as well as the nuclear absorption of primordial charmonia and the regeneration effects at the hadronization stage, and hence is a generalization of the two component model introduced by Grandchamp and Rapp. The improvements in this work are twofold; the addition of the initial local temperature profile and a consistent use of QCD NLO formula for both the dissociation cross section in the hadron gas and the thermal decay widths in the quark-gluon plasma for the charmonium states. The initial local temperature profile is determined from the assumption that the local entropy density is proportional to a formula involving the number densities of the number of participants and of the binary collisions that reproduces the multiplicities of charged particles at chemical freeze-out. The initial local temperature profile brings about a kink in the $R_{AA}$ curve due to the initial melting of $J/ψ$. The initially formed fireball, composed of weakly interacting quarks and gluons with thermal masses that are extracted from lattice QCD, follows an isentropic expansion with cylindrical symmetry. The fit reproduces well the Au+Au as well as the Cu+Cu data. The same method is applied to predict the $R_{AA}$ expected from the Pb+Pb collision at LHC energy.

preprint2010arXiv

Heavy quarkonium correlators at finite temperature: QCD sum rule approach

We investigate the properties of heavy quarkonia at finite temperature in detail using QCD sum rules. Extending previous analyses, we take into account a temperature dependent effective continuum threshold and derive constraints on the mass, the width, and the varying effective continuum threshold. We find that at least one of these quantities of a charmonium changes abruptly in the vicinity of the phase transition. We also calculate the ratio of the imaginary time correlator to its reconstructed one, $G/G_{\text{rec}}$, by constructing a model spectral function and compare it to the corresponding lattice QCD results. We demonstrate that the almost constant unity of $G/G_{\text{rec}}$ can be obtained from the destructive interplay of the changes in each part of the spectral modification which are extracted from QCD sum rules.

preprint2010arXiv

New Charmonium States in QCD Sum Rules: a Concise Review

In the past years there has been a revival of hadron spectroscopy. Many interesting new hadron states were discovered experimentally, some of which do not fit easily into the quark model. This situation motivated a vigorous theoretical activity. This is a rapidly evolving field with enormous amount of new experimental information. In the present report we include and discuss data which were released very recently. The present review is the first one written from the perspective of QCD sum rules (QCDSR), where we present the main steps of concrete calculations and compare the results with other approaches and with experimental data.

preprint2009arXiv

Sigma meson in QCD sum rules using a two quark current with derivatives

We study the $σ$ meson in QCD sum rules using a two quark interpolating field with derivatives. In the constituent quark model, the $σ$ meson is composed of a quark and an antiquark in the relative p-wave state and is thus expected to have a larger overlap with an interpolating field that measures the derivative of the relative quark wave-function. While the sum rule with a current without derivatives gives a pole mass of around 1 GeV, the present sum rule with a derivative current gives a mass of around 550 MeV and a width of 400 MeV, that could be identified with the $σ$ meson.

preprint2008arXiv

$D_sD^*$ molecule as an axial meson

We use QCD sum rules to study the possible existence of a $D_s\bar{D}^*+ D_s^*\bar{D}$ molecule with the quantum number $J^P=1^+$. We consider the contributions of condensates up to dimension eight and work at leading order in $α_s$. We obtain $m_{D_sD^*}=(3.96\pm0.10) \GeV$ around 100 MeV above the mass of the meson X(3872). The proposed state is a natural generalized state to the strangeness sector of the X(3872), which was also found to be consistent with a multiquark state from a previous QCD sum rule analysis.

preprint2004arXiv

Spin asymmetries in gamma + N --> K*bar + Theta^+

The photoproduction processes of the exotic $Θ^+(1540)$ baryon and the $\bar{K}^*$ meson from the nucleon targets, i.e., $γn \to K^{*-} Θ^+$ and $γp \to \bar{K}^{*0} Θ^+$, are investigated in a hadronic model. We consider $K$ and $K^*$ exchanges as well as the $s$ and $u$ channel nucleon and $Θ$ terms. Various spin asymmetries together with cross sections are first computed in order to study the production mechanisms and the parity of the $Θ^+(1540)$ baryon. Within the uncertainties arising from the model-dependence of the production mechanisms and several coupling constants, we find that some target-recoil double spin asymmetries, $C^{\rm TR}_{xx'}$ and $C^{\rm TR}_{xz'}$, are sensitive to the parity of $Θ^+$. In addition, the parity asymmetry of this reaction on the neutron target, which can be obtained by analyzing $K^*$ decay distribution, is found to be useful to estimate the $K^* N Θ$ coupling.

preprint2003arXiv

Charm production from proton-proton collisions

We evaluate the cross sections for charmed hadron production from proton-proton reactions $pp\to\bar D^0pΛ_c^+$ and $pp\to\bar D^{*0}pΛ_c^+$ using a hadronic Lagrangian. With empirical coupling constants and cutoff parameters in the form factors, sum of their cross sections at center-of-mass energy of 11.5 GeV is about 1 $μ$b and is comparable to measured inclusive cross section for charmed hadron production from proton-proton reactions. The cross section decreases to about 1 nb at 40 MeV above threshold.

preprint2003arXiv

Theta^+ baryon production in KN and NN reactions

We study $Θ^+(1540)$ productions in kaon-nucleon ($KN$) and nucleon-nucleon ($NN$) interactions by assuming that the $Θ^+$ is an isosinglet with $J^P = \frac12^+$. Possible $t$-channel diagrams with $K^*$ exchange are considered in both reactions as well as $K$ exchange in $NN$ reaction. The cross section for $np \to Λ^0 Θ^+$, which has not been considered in previous calculations, is found to be about a factor of 5 larger than that for $np \to Σ^0 Θ^+$ due to the large coupling of $KNΛ$ interaction. The cross sections are obtained by setting $g_{KNΘ}^{}=1$ and varying the ratio of $g_{K^*NΘ}^{}/g_{KNΘ}^{}$ so that future experimental data can be used to estimate these couplings. We also find that the isospin relations hold for these reactions.

preprint2001arXiv

Couplings between pion and charmed mesons

We compute the couplings $DD^*π$,$D_1D^*π$,$D^*D^*π$, $D_1 D_1π$ using QCD sum rules. These couplings are important inputs in the meson exchange model calculations used to estimate the amount of $J/ψ$ absorption due to pions and rho mesons in heavy ion collisions. Our sum rules are constructed at the first order in the pion momentum $p_μ$, which give the couplings that are not trivially related to the soft-pion theorem. Our calculated couplings, which somewhat depend upon the values of the heavy meson decay constants, are $g_{DD^*π}=8.2 \pm 0.1$, $g_{D_1 D^* π} = 15.8 \pm 2 {\rm GeV}$, $g_{D^* D^* π}= 0.3 \pm 0.03$ and $g_{D_1 D_1 π}= 0.17 \pm 0.04$.

preprint1999arXiv

Meson-baryon couplings and the $F/D$ ratio from QCD sum rules

Motivated by the recent work on the calculation of the $πNN$ coupling constant using QCD sum rule beyond the chiral limit, we construct the corresponding sum rules for the couplings, $ηNN$, $πΞΞ$, $ηΞΞ$, $πΣΣ$ and $ηΣΣ$. In constructing the $η$-baryon sum rules, we use the second moment of the $η$ wave function, which we obtain from the pion wave function after SU(3) rotation. In the SU(3) symmetric limit, we can identify the term responsible for the $F/D$ ratio in the OPE, which after the sum rule analysis gives $F/D\sim 0.2$, a factor of 3 smaller than from other studies. We also present a qualitative analysis including the SU(3) breaking terms.

preprint1998arXiv

J/psi and eta_c in the nuclear medium: QCD sum rule approach

We investigate the masses of the lowest c\bar{c} states, the J/psi and eta_c, in nuclear matter using QCD sum rules. Up to dimension four, the differences between the operator product expansions in vacuum and in medium arise from the density-dependent change in the gluon condensate and from a new contribution proportional to the nucleon expectation value of the twist-2 gluon operator. Both terms together give an attractive shift of about 5-10 MeV to the J/psi and eta_c masses in nuclear matter.

preprint1997arXiv

Vector mesons in medium with finite three momentum

We formulate a QCD sum rule to find the three momentum dependence of the peak position of the vector meson spectral density in nuclear medium. We find less than 2 % (0.1 %) shift of the peak position at nuclear matter density and at ($q =0.5$GeV/c) for the $ρ,ω(ϕ$) mesons. However, at higher density and three momentum, its effect becomes non-negligible. The relevance of our result in the dilepton spectrum in A-A and p-A reaction is discussed.

preprint1995arXiv

Thermal effects on dilepton production from $π-π$ annihilation

We study finite temperature effects on dilepton production from pion-pion annihilation in hot hadronic matter. The softening of the pion dispersion relation in a medium is found to enhance the production rate of dileptons with invariant masses in the region of $2m^*_π(T)<M<m_ρ$. On the other hand, the reduction of the pion electromagnetic form factor at finite temperature leads to a suppression of the dilepton production rate. Including both effects, we have found that the dilepton yield is slightly enhanced in the invariant mass region $M= 270\sim 600$ MeV but is suppressed around the vector meson resonance. We further discuss the relevance of our results to recent experimental data from the S+Au collisions at CERN/SPS energies by the CERES collaboration.