Search

The interaction between nucleon and charmonia ($J/\psi$ and $\eta_c$) is expected to deepen our understanding of various aspects in nonperturbative QCD ranging from the origin of nucleon mass to $J/\psi$ mass modification in nuclear medium and properties of hidden-charm pentaquark states. Here, we present the low-energy $NJ/\psi$ and $N\eta_c$ interactions based on ($2+1$) flavor lattice QCD simulations with nearly physical pion mass $m_\pi=146$ MeV. The interactions, extracted from the spacetime correlations of the nucleon and charmonium system by using the HAL QCD method, are found to be attractive in all distances and manifest a characteristic long-range tail consistent with the two-pion exchange interaction. The resulting scattering lengths are around $0.3$ fm, $0.4$ fm and $0.2$ fm for $NJ/\psi$ with spin $3/2$, with spin $1/2$, and $N\eta_c$, respectively. Our results are orders of magnitude larger than those from the photoproduction experiments assuming the vector meson dominance., Comment: 9 pages, 3 figures, contribution to the 41st International Symposium on Lattice Field Theory (Lattice 2024), 28 July - 3 August 2024, Liverpool, UK

The integration of deep learning techniques and physics-driven designs is reforming the way we address inverse problems, in which accurate physical properties are extracted from complex data sets. This is particularly relevant for quantum chromodynamics (QCD), the theory of strong interactions, with its inherent limitations in observational data and demanding computational approaches. This perspective highlights advances and potential of physics-driven learning methods, focusing on predictions of physical quantities towards QCD physics, and drawing connections to machine learning(ML). It is shown that the fusion of ML and physics can lead to more efficient and reliable problem-solving strategies. Key ideas of ML, methodology of embedding physics priors, and generative models as inverse modelling of physical probability distributions are introduced. Specific applications cover first-principle lattice calculations, and QCD physics of hadrons, neutron stars, and heavy-ion collisions. These examples provide a structured and concise overview of how incorporating prior knowledge such as symmetry, continuity and equations into deep learning designs can address diverse inverse problems across different physical sciences., Comment: 14 pages, 5 figures, submitted version to Nat Rev Phys

A new framework is introduced to connect between a charge symmetry breaking (CSB) energy density functional (EDF) and the low-energy constants derived from quantum chromodynamics (QCD). By constructing a QCD-based CSB EDF, this method provides new insights into the Okamoto-Nolen-Schiffer anomaly, a long-standing puzzle in the energy differences of mirror nuclei that lacks a robust microscopic explanation. Using examples such as $ {}^{17} \mathrm{F} $-$ {}^{17} \mathrm{O} $, $ {}^{15} \mathrm{O} $-$ {}^{15} \mathrm{N} $, $ {}^{41} \mathrm{Sc} $-$ {}^{41} \mathrm{Ca} $, and $ {}^{39} \mathrm{Ca} $-$ {}^{39} \mathrm{K} $, we demonstrate that the proposed interaction effectively resolves the anomaly within the range of theoretical uncertainties., Comment: 8 pages, 2 figures, 2 tables; 10th International Conference on Quarks and Nuclear Physics (QNP2024)

We present a realistic lattice QCD study on low-energy $N$-$J/\psi$ and $N$-$\eta_c$ interactions based on (2+1) flavor configurations with nearly physical pion mass $m_\pi=146$ MeV. The interactions, extracted from the spacetime correlations of nucleon and charmonium system by using the HAL QCD method, are found to be attractive in all distances and possess a characteristic long-range tail consistent with the two-pion exchange potential. The resulting $S$-wave scattering lengths are $0.30(2)\left(^{+0}_{-2}\right)$ fm, $0.38(4)\left(^{+0}_{-3}\right)$ fm, and $0.21(2)\left(^{+0}_{-1}\right)$ fm for spin-$3/2$ $N$-$J/\psi$, spin-$1/2$ $N$-$J/\psi$, and spin-$1/2$ $N$-$\eta_c$, respectively. Our results are orders of magnitude larger than those from the photoproduction experiments assuming the vector meson dominance. Our findings may provide deeper understanding of the nonperturbative QCD phenomena ranging from the origin of nucleon mass to the in-medium $J/\psi$ mass modification as well as the properties of hidden-charm pentaquark states., Comment: v2, minor modifications are made to match with the published version

In this study, we develop a deep learning method to learn hadronic interactions unsupervisedly from the correlation functions calculated in lattice QCD simulations. We present our approach of using deep neural networks to model the inter-hadron potentials that are learned from Nambu-Bethe-Salpeter (NBS) wave functions. This enables the incorporation of most general forms of potentials into the Schr\"odinger-type equation for detailed analysis of hadronic interactions. Our results include validations with separable potentials, as well as the local and non-local potentials for the $\Omega_{ccc}-\Omega_{ccc}$ system. The neural networks accurately capture the essential features of these interactions, providing a reliable tool for predicting and analyzing hadron scattering properties, potentially bridging the experimental observables and lattice QCD data., Comment: 8 pages, 4 figures, contribution to the 41st International Symposium on Lattice Field Theory (Lattice 2024), July 28th - August 3rd, 2024, University of Liverpool

Numerous biological processes accelerate as temperatures increase, but the period of circadian rhythms remains constant, known as temperature compensation, while synchronizing with the 24h light-dark cycle. We theoretically explores the possible relevance of waveform distortions in circadian gene-protein dynamics to the temperature compensation and synchronization. Our analysis of the Goodwin model provides a coherent explanation of most of temperature compensation hypotheses. Using the renormalization group method, we analytically demonstrate that the decreasing phase of circadian protein oscillations should lengthen with increasing temperature, leading to waveform distortions to maintain a stable period. This waveform-period correlation also occurs in other oscillators like Lotka-Volterra and van der Pol models. A reanalysis of known data nicely confirms our findings on waveform distortion and its impact on synchronization range. Thus we conclude that circadian rhythm waveforms are fundamental to both temperature compensation and synchronization., Comment: 49 pages, 9 figures, 2 tables

We study the doubly charmed tetraquark state $T_{cc}^+$ by the HAL QCD method applied to the $D^*D$ system in $(2+1)$ flavor lattice QCD at nearly physical pion mass, $m_\pi= 146$ MeV. We obtain the attractive potential at all distances in the $S$-wave of the isoscalar $D^* D$ system, whose long distance behavior is well described by the two-pion exchange (TPE), and it generates a virtual pole near $D^* D$ threshold with a pole position $E_{\rm pole} = -59 (^{+53}_{-99}) (^{+2}_{-67})$ keV and an inverse scattering length $1/a_0=0.05(5)(^{+2}_{-2})$ fm$^{-1}$. The virtual pole turns into a loosely bound state pole if the pion mass in the TPE potential is extrapolated to the physical value, $m_\pi =135$ MeV. The potential at the physical pion mass is shown to give a semi-quantitative description of the $D^0 D^0\pi^+$ mass spectrum at the LHCb., Comment: 9 pages, 7 figures, Talk presented at the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermi National Accelerator Laboratory, Batavia, Illinois, USA

We summarize the recent progress on the determination of the charge symmetry breaking term of nuclear energy density functionals. We point out that the strength of the term determined theoretically is remarkably smaller than that determined phenomenologically, which is still an open question., Comment: The 7th International Conference on Collective Motion in Nuclei under Extreme Conditions (COMEX7)

The emergence of cryptoassets has sparked a paradigm shift in the world of finance and investment, ushering in a new era of digital assets with profound implications for the future of currency and asset management. A recent study showed that during the bubble period around the year, 2018, the price of cryptoasset, XRP has a strong anti correlation with the largest singular values of the correlation tensors obtained from the weekly XRP transaction networks. In this study, we provide a detailed analysis of the method of correlation tensor spectra for XRP transaction networks. We calculate and compare the distribution of the largest singular values of the correlation tensor using the random matrix theory with the largest singular values of the empirical correlation tensor. We investigate the correlation between the XRP price and the largest singular values for a period spanning two years. We also uncover the distinct dependence between XRP price and the singular values for bubble and non-bubble periods. The significance of time evolution of singular values is shown by comparison with the evolution of singular values of the reshuffled correlation tensor. Furthermore, we identify a set of driver nodes in the transaction networks that drives the market during the bubble period using the singular vectors., Comment: Main text: 11 pages, 10 figures. Supplementary information: 3 pages, 3 figures

An approach is proposed to link the charge symmetry breaking (CSB) nuclear interaction and the low-energy constants in quantum chromodynamics (QCD) by matching the CSB effect in nuclear matter. The resulting CSB interaction is applied to study the Okamoto-Nolen-Schiffer anomaly, still lacking a satisfactory microscopic understanding, on the energy differences of mirror nuclei by taking $ {}^{17} \mathrm{F} $-$ {}^{17} \mathrm{O} $, $ {}^{15} \mathrm{O} $-$ {}^{15} \mathrm{N} $, $ {}^{41} \mathrm{Sc} $-$ {}^{41} \mathrm{Ca} $, and $ {}^{39} \mathrm{Ca} $-$ {}^{39} \mathrm{K} $ as typical examples. The magnitude and sign of the QCD-based CSB interactions are found to resolve the anomaly successfully within theoretical uncertainties., Comment: 6 pages, 2 figures, 6 tables in the main text; 6 pages, 2 figures, 3 tables in the supplemental material

Cryptoassets are growing rapidly worldwide. One of the large cap cryptoassets is XRP. In this article, we focus on analyzing transaction data for the 2017-2018 period that consist one of the significant XRP market price bursts. We construct weekly weighted directed networks of XRP transactions. These weekly networks are embedded on continuous vector space using a network embedding technique that encodes structural regularities present in the network structure in terms of node vectors. Using a suitable time window we calculate a correlation tensor. A double singular value decomposition of the correlation tensor provides key insights about the system. The significance of the correlation tensor is captured using a randomized correlation tensor. We present a detailed dependence of correlation tensor on model parameters., Comment: 8 pages, 7 figures

We theoretically examine a continuity between atomic and molecular Fermi superfluids in a Bose-Fermi mixture near the Feshbach resonance. Considering a two-channel model describing the Feshbach resonance between Fermi and Bose atoms, we have constructed the mean-field framework based on the perturbative expansion of the Feshbach atom-dimer coupling. The resulting effective Hamiltonian exhibits not only the continuity between atom-atom to molecule-molecule Cooper pairings but also becomes equivalent to the two-band-superconductor model with Suhl-Matthias-Walker type pair-exchange coupling. We demonstrate how these atomic and molecular Fermi superfluids coexist within the two-band-like superfluid theory. The pair-exchange coupling and resulting superfluid gaps are found to be strongly enhanced near the Feshbach resonance due to the interplay between the infrared singularity of Bogoliubov phonons and their Landau damping arising from the coupling with fermions. The pair-exchange coupling can be probed via the observation of the intrinsic Josephson effect between atomic and molecular superfluids., Comment: 15 pages, 5 figures

The doubly charmed tetraquark $T^+_{cc}$ recently discovered by the LHCb Collaboration is studied on the basis of $(2+1)$-flavor lattice QCD simulations of the $D^*D$ system with nearly physical pion mass $m_\pi=146$ MeV. The interaction of $D^*D$ in the isoscalar and $S$-wave channel, derived from the hadronic spacetime correlation by the HAL QCD method, is attractive for all distances and leads to a near-threshold virtual state with a pole position $E_\text{pole}=-59\left(^{+53}_{-99}\right)\left(^{+2}_{-67}\right)$ keV and a large scattering length $1/a_0=0.05(5)\left(^{+2}_{-2}\right)~\text{fm}^{-1}$. The virtual state is shown to evolve into a loosely bound state as $m_\pi$ decreases to its physical value by using a potential modified to $m_\pi=135$ MeV based on the pion-exchange interaction. Such a potential is found to give a semiquantitative description of the LHCb data on the $D^0D^0\pi^+$ mass spectrum. Future study is necessary to perform physical-point simulations with the isospin-breaking and open three-body-channel effects taken into account., Comment: 7 pages, 4 figures + supplemental material; v3, match the published version in PRL

The existence of a nucleon-$\phi$ (N-$\phi$) bound state has been subject of theoretical and experimental investigations for decades. In this letter, indication of a \pphi bound state is found, using for the first time two-particle correlation functions as alternative to invariant mass spectra. Newly available lattice calculations for the spin 3/2 \Nphi interaction by the HAL QCD collaboration are used to constrain the spin 1/2 counterpart from the fit of the experimental \pphi correlation function measured by ALICE. The corresponding scattering length and effective range are $f_0^{(1/2)}=\left(-1.54^{+0.53}_{-0.53}(\mathrm{stat.})^{+0.16}_{-0.09}(\mathrm{syst.})+i\cdot0.00^{+0.35}_{-0.00}(\mathrm{stat.})^{+0.16}_{-0.00}(\mathrm{syst.})\right)$~fm and $d_0^{(1/2)}=\left(0.39^{+0.09}_{-0.09}(\mathrm{stat.})^{+0.02}_{-0.03}(\mathrm{syst.})+i\cdot0.00^{+0.00}_{-0.04}(\mathrm{stat.})^{+0.00}_{-0.02}(\mathrm{syst.})\right)$~fm, respectively. The results imply the appearance of a \pphi bound state with an estimated binding energy in the range of $12.8-56.1$ MeV.

Cryptoassets are becoming essential in the digital economy era. XRP is one of the large market cap cryptoassets. Here, we develop a novel method of correlation tensor spectra for the dynamical XRP networks, which can provide an early indication for XRP price. A weighed directed weekly transaction network among XRP wallets is constructed by aggregating all transactions for a week. A vector for each node is then obtained by embedding the weekly network in continuous vector space. From a set of weekly snapshots of node vectors, we construct a correlation tensor. A double singular value decomposition of the correlation tensors gives its singular values. The significance of the singular values is shown by comparing with its randomize counterpart. The evolution of singular values shows a distinctive behavior. The largest singular value shows a significant negative correlation with XRP/USD price. We observe the minimum of the largest singular values at the XRP/USD price peak during the first week of January 2018. The minimum of the largest singular value during January 2018 is explained by decomposing the correlation tensor in the signal and noise components and also by evolution of community structure., Comment: 12 pages 8 figures, Supplementary Information: 4 pages, 6 figures

First results on the interaction between the $\phi$-meson and the nucleon ($N$) are presented based on the ($2+1$)-flavor lattice QCD simulations with nearly physical quark masses. Using the HAL QCD method, the spacetime correlation of the $N$-$\phi$ system in the spin 3/2 channel is converted into the $N$-$\phi$ scattering phase shift through the interaction potential. The $N$-$\phi$ potential appears to be a combination of a short-range attractive core and a long-range attractive tail. The latter is found to be consistent with the two-pion exchange (TPE) obtained from the interaction between a color-dipole and the nucleon. The resultant scattering length and effective range for $m_{\pi}=$ 146.4 MeV are $ a^{(3/2)}_0=-1.43(23)_{\rm stat.}\left(^{+36}_{-06}\right)_{\rm syst.} {\rm fm}$ and $ r^{(3/2)}_{\rm eff}=2.36(10)_{\rm stat.}\left(^{+02}_{-48}\right)_{\rm syst.} {\rm fm}$, respectively. The magnitude of the scattering length is shown to have nontrivial dependence of $m_{\pi}$ and is sensitive to the existence of the long-range tail from TPE., Comment: 8 pages, 6 figures, matches version published in PRD

Fully general-relativistic binary-neutron-star (BNS) merger simulations with quark-hadron crossover (QHC) equations of state (EOSs) are studied for the first time. In contrast to EOSs with purely hadronic matter or with a first-order quark-hadron phase transition (1PT), in the transition region QHC EOSs show a peak in sound speed, and thus a stiffening. We study the effects of such stiffening in the merger and post-merger gravitational (GW) signals. Through simulations in the binary-mass range $2.5 < M/M_{\odot} < 2.75$, characteristic differences due to different EOSs appear in the frequency of the main peak of the post-merger GW spectrum ($f_2$), extracted through Bayesian inference. In particular, we found that (i) for lower-mass binaries, since the maximum baryon number density ($n_{\rm max}$) after the merger stays below $3\text{--}4$ times the nuclear-matter density ($n_0$), the characteristic stiffening of the QHC models in that density range results in a lower $f_2$ than that computed for the underlying hadronic EOS and thus also than that for EOSs with a 1PT, (ii) for higher-mass binaries, where $n_{\rm max}$ may exceed $4\text{--}5 n_0$ depending on the EOS model, whether $f_2$ in QHC models is higher or lower than that in the underlying hadronic model depends on the height of the sound-speed peak. Comparing the values of $f_2$ for different EOSs and BNS masses gives important clues on how to discriminate different types of quark dynamics in the high-density end of EOSs and is relevant to future kilohertz GW observations with third-generation GW detectors., Comment: 6+8 pages, 5+6 figures

A set of optimized interpolating operators which are dominantly coupled to each eigenstate of two baryons on the lattice is constructed by the HAL QCD method. To test its validity, we consider heavy dibaryons $\Omega_{3Q}\Omega_{3Q}$ ($Q=s,c$) calculated by (2+1)-flavor lattice QCD simulations with nearly physical pion mass. The optimized two-baryon operators are shown to provide effective energies of the ground and excited states separately stable as a function of the Euclidean time. Also they agree to the eigenenergies in a finite lattice box obtained from the leading-order HAL QCD potential $V({\boldsymbol{r}})$ within statistical errors. The overlapping factors between the optimized sink operators and the state created by the wall-type source operator indicate that $V( \boldsymbol{r})$ can be reliably extracted, no matter whether the spacetime correlation of two baryons is dominated by the ground state or the excited state. It is suggested that the optimized set of operators is useful for variational studies of hadron-hadron interactions., Comment: Accepted for publication in Physical Review D (8 pages, 7 figures)

We study the convergence of the derivative expansion in HAL QCD method from the finite volume analysis. Employing the (2+1)-flavor lattice QCD data obtained at nearly physical light quark masses $(m_\pi, m_K) \simeq (146, 525)$ MeV and the physical charm quark mass, we study two representative systems, $\Omega\Omega$ and $\Omega_{ccc}\Omega_{ccc}$ in the $^1S_0$ channel, where both systems were found to have a shallow bound state in our previous studies. The HAL QCD potentials are determined at the leading-order in the derivative expansion, from which finite-volume eigenmodes are obtained. Utilizing the eigenmode projection, we find that the correlation functions are dominated by the ground state (first excited state) in the case of $\Omega\Omega$ ($\Omega_{ccc}\Omega_{ccc}$). In both $\Omega\Omega$ and $\Omega_{ccc}\Omega_{ccc}$, the spectra obtained from eigenmode-projected temporal correlators are found to be consistent with those from the HAL QCD potential for both the ground and first excited state. These results show that the derivative expansion is well converged in these systems, and also provide a first explicit evidence that the HAL QCD method enables us to reliably extract the binding energy of the ground state even from the correlator dominated by excited scattering states., Comment: 8 pages, 4 figures, Talk given at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 Jul 2021, Zoom/Gather@MIT, USA

We present a first study on a pair of triply charmed baryons, $\Omega_{ccc}\Omega_{ccc}$ in the $^1S_0$ channel, on the basis of the HAL QCD method. The measurements are perfomed on the $(2+1)$-flavor lattice QCD configurations with nearly physical light-quark masses and physical charm-quark mass. We show that the system with the Coulomb repulsion taking into account the charge form factor of $\Omega_{ccc}$ leads to the scattering length $a^\mathrm{C}_0\simeq-19$ fm and the effective range $r^\mathrm{C}_\mathrm{eff}\simeq0.45$ fm, which indicates $\Omega_{ccc}\Omega_{ccc}$ is located in the unitary regime., Comment: 8 pages, 4 figures, Talk presented at the 38th International Symposium on Lattice Field Theory, LATTICE2021, July, 26th-30th 2021, Zoom/Gather@Massachusetts Institute of Technology

We investigate the distribution of energy-momentum tensor (EMT) around a static quark in the deconfined phase of SU(3) Yang-Mills theory. The EMT defined through the gradient-flow formalism is used for the numerical analysis of the EMT distribution around the Polyakov loop with the continuum extrapolation. Using EMT, one can study the mechanical distortion of the color gauge field induced by the static charge. We find substantial separation in the absolute values of the EMT eigenvalues which is not observed in Maxwell theory. The separation grows as temperature is lowered toward the critical temperature. The lattice data also indicate the thermal screening at long distance and the perturbative behavior at short distance., Comment: 8 pages, 11 figures, Proceedings of the 38th International Symposium on Lattice Field Theory, 26-30 July 2021

The recent NICER measurement of the radius of the neutron star PSR J0740+6620, and the inferred small variation of radii from 1.4$M_\odot$ to 2.1$M_\odot$, reveal key features of the equation of state of neutron star matter. The pressure rises rapidly in the regime of baryon density $n \sim$ 2-4 times nuclear saturation density, $n_0$ -- the region where we expect hadronic matter to be undergoing transformation into quark matter -- and the pressure in the nuclear regime is greater than predicted by microscopic many-body variational calculations of nuclear matter. To incorporate these insights into the microscopic physics from the nuclear to the quark matter regimes, we construct an equation of state, QHC21, within the framework of quark-hadron crossover (QHC). We include nuclear matter results primarily based on the state-of-the-art chiral effective field theory, but also note results of using nuclear matter variational calculations based on empirical nuclear forces. We employ explicit nuclear degrees of freedom only up to $n \sim 1.5n_0$, in order to explore the possibility of further physical degrees of freedom than nucleonic here. The resulting QHC21, which has a peak in sound velocity in $\sim 2$-$4 n_0$, is stiffer than the earlier QHC19 below 2$n_0$, predicting larger radii in substantial agreement with the NICER data., Comment: 10 pages, 11 figures; v2) 12 pages, 11 figures

A pair of triply charmed baryons, $\Omega_{ccc}\Omega_{ccc}$, is studied as an ideal dibaryon system by (2+1)-flavor lattice QCD with nearly physical light-quark masses and the relativistic heavy quark action with the physical charm quark mass. The spatial baryon-baryon correlation is related to their scattering parameters on the basis of the HAL QCD method. The $\Omega_{ccc}\Omega_{ccc}$ in the ${^1S_0}$ channel taking into account the Coulomb repulsion with the charge form factor of $\Omega_{ccc}$ leads to the scattering length $a^{\rm C}_0\simeq -19~\text{fm}$ and the effective range $r^{\rm C}_{\mathrm{eff}}\simeq 0.45~\text{fm}$. The ratio $r^{\rm C}_{\mathrm{eff}}/a^{\rm C}_0 \simeq -0.024$, whose magnitude is considerably smaller than that of the dineutron ($-0.149$), indicates that $\Omega_{ccc}\Omega_{ccc}$ is located in the unitary regime., Comment: 6 pages, 4 figures; Supplemental Material(1 page, 1 figure); Accepted for publication in Physical Review Letters

Energy momentum tensor (EMT) characterizes the response of the vacuum as well as the thermal medium under the color electromagnetic fields. We define the EMT by means of the gradient flow formalism and study its spatial distribution around a static quark in the deconfined phase of SU(3) Yang-Mills theory on the lattice. Although no significant difference can be seen between the EMT distributions in the radial and transverse directions except for the sign, the temporal component is substantially different from the spatial ones near the critical temperature $T_c$. This is in contrast to the prediction of the leading-order thermal perturbation theory. The lattice data of the EMT distribution also indicate the thermal screening at long distance and the perturbative behavior at short distance., Comment: 12 pages, 12 figures

The $\Delta\Delta$ dibaryon resonance $d^\ast (2380)$ with $(J^P, I)=(3^+, 0)$ is studied theoretically on the basis of the 3-flavor lattice QCD simulation with heavy pion masses ($m_\pi =679, 841$ and $1018$ MeV). By using the HAL QCD method, the central $\Delta$-$\Delta$ potential in the ${}^7S_3$ channel is obtained from the lattice data with the lattice spacing $a\simeq 0.121$ fm and the lattice size $L\simeq 3.87$ fm. The resultant potential shows a strong short-range attraction, so that a quasi-bound state corresponding to $d^\ast (2380)$ is formed with the binding energy $25$-$40$ MeV below the $\Delta\Delta$ threshold for the heavy pion masses. The tensor part of the transition potential from $\Delta\Delta$ to $NN$ is also extracted to investigate the coupling strength between the $S$-wave $\Delta\Delta$ system with $J^P=3^+$ and the $D$-wave $NN$ system. Although the transition potential is strong at short distances, the decay width of $d^\ast (2380)$ to $NN$ in the $D$-wave is kinematically suppressed, which justifies our single-channel analysis at the range of the pion mass explored in this study., Comment: 8 pages, 6 figures

A novel quantum-classical hybrid scheme is proposed to efficiently solve large-scale combinatorial optimization problems. The key concept is to introduce a Hamiltonian dynamics of the classical flux variables associated with the quantum spins of the transverse-field Ising model. Molecular dynamics of the classical fluxes can be used as a powerful preconditioner to sort out the frozen and ambivalent spins for quantum annealers. The performance and accuracy of our smooth hybridization in comparison to the standard classical algorithms (the tabu search and the simulated annealing) are demonstrated by employing the MAX-CUT and Ising spin-glass problems., Comment: 14 pages, 10 figures; v2: 17 pages, 14 figures, supplementary note added, to be published in Scientific Reports

The $S$-wave $\Lambda\Lambda$ and $N \Xi$ interactions are studied on the basis of the (2+1)-flavor lattice QCD simulations close to the physical point ($m_\pi \simeq 146{\rm{MeV}}$ and $m_K \simeq 525{\rm{MeV}}$). Lattice QCD potentials in four different spin-isospin channels are extracted by using the coupled-channel HAL QCD method and are parametrized by analytic functions to calculate the scattering phase shifts. The $\Lambda \Lambda$ interaction at low energies shows only a weak attraction, which does not provide a bound or resonant dihyperon. The $N\Xi$ interaction in the spin-singlet and isospin-singlet channel is most attractive and lead the $N\Xi$ system near unitarity. Relevance to the strangeness=$-2$ hypernuclei as well as to two-baryon correlations in proton-proton, proton-nucleus and nucleus-nucleus collisions is also discussed., Comment: 16 pages, 8 figures

We explore the distribution of the energy momentum tensor (EMT) around quark--anti-quark and single quark at nonzero temperature in SU(3) Yang-Mills gauge theory by extending our previous study on the EMT distribution in static quark--anti-quark systems in vacuum. We discuss the disappearance of the flux tube structure observed in the vacuum simulation. We investigate the total stress acting on the mid-plane between a quark and an anti-quark and show that it agrees with the force obtained from the derivative of the free energy. The color Debye screening effect in the deconfined phase is also discussed in terms of the EMT distribution., Comment: 7 pages, 3 figures, Proceedings for the 37th Annual International Symposium on Lattice Field Theory (Lattice 2019)

The momentum correlation functions of baryon pairs, which reflects the baryon-baryon interaction at low energies, are investigated for multi-strangeness pairs ($\Omega\Omega$ and $N\Omega$) produced in relativistic heavy-ion collisions. We calculate the correlation functions based on an expanding source model constrained by single-particle distributions. The interaction potentials are taken from those obtained from recent lattice QCD calculations at nearly physical quark masses. Experimental measurements of these correlation functions for different system sizes will help to disentangle the strong interaction between baryons and to unravel the possible existence of strange dibaryons., Comment: 11 pages, 8 figures (reference format is modified in v2)

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., Comment: 8 pages, 4 figures

We investigate superfluid phase transitions of asymmetric nuclear matter at finite temperature ($T$) and density ($\rho$) with a low proton fraction ($Y_{\rm p} \le 0.2$) which is relevant to the inner crust and outer core of neutron stars. A strong-coupling theory developed for two-component atomic Fermi gases is generalized to the four-component case and is applied to the system of spin-$1/2$ neutrons and protons. The empirical phase shifts of neutron-neutron (nn), proton-proton (pp) and neutron-proton (np) interactions up to $k = 2$ ${\rm fm}^{-1}$ are described by multi-rank separable potentials. We show that (i) the critical temperature of the neutron superfluidity $T_{\rm c}^{\rm nn}$ at $Y_{\rm p}=0$ agrees well with Monte Carlo data at low densities and takes a maximum value $T_{\rm c}^{\rm nn}=1.68$ MeV at $\rho/\rho_0 = 0.14$ with $\rho_0=0.17$ fm$^{-3}$, (ii) the critical temperature of the proton superconductivity $T_{\rm c}^{\rm pp}$ for $Y_{\rm p} \le 0.2$ is substantially suppressed at low densities due to np-pairing fluctuations and starts to dominate over $T_{\rm c}^{\rm nn}$ only above $\rho/\rho_0 = 0.70$ $(0.77)$ for $Y_p =0.1$ $(0.2)$, and (iii) the deuteron condensation temperature $T_{\rm c}^{\rm d}$ is suppressed at $Y_{\rm p}\le 0.2$ due to the large mismatch of the two Fermi surfaces., Comment: 23 pages, 12 figures

A moderately strong vector repulsion between quarks in dense quark matter is needed to explain how a quark core can support neutron stars heavier than two solar masses. We study this repulsion, parametrized by a four-fermion interaction with coupling g_V, in terms of non-perturbative gluon exchange in QCD in the Landau gauge. Matching the energy of quark matter, g_V n_q^2 (where n_q is the number density of quarks) with the quark exchange energy calculated in QCD with a gluon propagator parametrized by a finite gluon mass m_g and a frozen coupling alpha_s, at moderate quark densities, we find that gluon masses m_g in the range 200 - 600 MeV and alpha_s = 2 - 4 lead to a g_V consistent with neutron star phenomenology. Estimating the effects of quark masses and a color-flavor-locked (CFL) pairing gap, we find that g_V can be well approximated by a flavor-symmetric, decreasing function of density. We briefly discuss similar matchings for the isovector repulsion and for the pairing attraction., Comment: 9 pages, 6 figures

We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations., Comment: 8 pages, 11 figures; v3) published version, appendix was added

There exist two methods to study two-baryon systems in lattice QCD: the direct method which extracts eigenenergies from the plateaux of the temporal correlator and the HAL QCD method which extracts observables from the non-local potential associated with the tempo-spatial correlator. Although the two methods should give the same results theoretically, qualitatively different results have been reported. Recently, we pointed out that the separation of the ground state from the excited states is crucial to obtain sensible results in the former, while both states provide useful signals in the latter. In this paper, we identify the contribution of each state in the direct method by decomposing the two-baryon correlators into the finite-volume eigenmodes obtained from the HAL QCD method. We consider the $\Xi\Xi$ system in the $^1$S$_0$ channel at $m_\pi = 0.51$ GeV in 2+1 flavor lattice QCD using the wall and smeared quark sources. We demonstrate that the "pseudo-plateau" at early time slices (t = 1~2 fm) from the smeared source in the direct method indeed originates from the contamination of the excited states, and the true plateau with the ground state saturation is realized only at t > 5~15 fm corresponding to the inverse of the lowest excitation energy. We also demonstrate that the two-baryon operator can be optimized by utilizing the finite-volume eigenmodes, so that (i) the finite-volume energy spectra from the HAL QCD method agree with those from the optimized temporal correlator and (ii) the correct spectra would be accessed in the direct method only if highly optimized operators are employed. Thus we conclude that the long-standing issue on the consistency between the L\"uscher's finite volume method and the HAL QCD method for two baryons is now resolved: They are consistent with each other quantitatively only if the excited contamination is properly removed in the former., Comment: 34 pages, 17 figurs, published version

We study the spatial distribution of the stress tensor around static quark-anti-quark pair in SU(3) lattice gauge theory. In particular, we reveal the transverse structure of the stress tensor distribution in detail by taking the continuum limit. The Yang-Mills gradient flow plays a crucial role to make the stress tensor well-defined and derivable from the numerical simulations on the lattice., Comment: 7 pages, 4 figures, Fig.3 is corrected, Proceedings for the 36th Annual International Symposium on Lattice Field Theory (Lattice 2018)

The nucleon($N$)-Omega($\Omega$) system in the S-wave and spin-2 channel ($^5$S$_2$) is studied from the (2+1)-flavor lattice QCD with nearly physical quark masses ($m_\pi \simeq 146$~MeV and $m_K \simeq 525$~MeV). The time-dependent HAL QCD method is employed to convert the lattice QCD data of the two-baryon correlation function to the baryon-baryon potential and eventually to the scattering observables. The $N\Omega$($^5$S$_2$) potential, obtained under the assumption that its couplings to the D-wave octet-baryon pairs are small, is found to be attractive in all distances and to produce a quasi-bound state near unitarity: In this channel, the scattering length, the effective range and the binding energy from QCD alone read $a_0= 5.30(0.44)(^{+0.16}_{-0.01})$~fm, $r_{\rm eff} = 1.26(0.01)(^{+0.02}_{-0.01})$~fm, $B = 1.54(0.30)(^{+0.04}_{-0.10})$~MeV, respectively. Including the extra Coulomb attraction, the binding energy of $p\Omega^-$($^5$S$_2$) becomes $B_{p\Omega^-} = 2.46(0.34)(^{+0.04}_{-0.11})$~MeV. Such a spin-2 $p\Omega^-$ state could be searched through two-particle correlations in $p$-$p$, $p$-nucleus and nucleus-nucleus collisions., Comment: 16 pages, 6 figures, a reference added

The properties of open strange meson $K_1^\pm$ in nuclear matter are estimated in the QCD sum rule approach. We obtain a relation between the in-medium mass and width of $K_1^-$ ($K_1^+$) in nuclear matter, and show that the upper limit of the mass shift is as large as -249 (-35) MeV. The spectral modification of the $K_1$ meson is possible to be probed by using kaon beams at J-PARC. Such measurement together with that of $K^*$ will shed light on how chiral symmetry is partially restored in nuclear matter., Comment: 12 pages, 3 figures, accepted for publication in Phys. Lett. B

The $\Xi\Xi$ interaction in the $^1$S$_0$ channel is studied to examine the convergence of the derivative expansion of the non-local HAL QCD potential at the next-to-next-to-leading order (N$^2$LO). We find that (i) the leading order potential from the N$^2$LO analysis gives the scattering phase shifts accurately at low energies, (ii) the full N$^2$LO potential gives only small correction to the phase shifts even at higher energies below the inelastic threshold, and (iii) the potential determined from the wall quark source at the leading order analysis agrees with the one at the N$^2$LO analysis except at short distances, and thus, it gives correct phase shifts at low energies. We also study the possible systematic uncertainties in the HAL QCD potential such as the inelastic state contaminations and the finite volume artifact for the potential and find that they are well under control for this particular system., Comment: 23 pages, 13 figures

The spatial distribution of the stress tensor around the quark--anti-quark ($Q\bar{Q}$) pair in SU(3) lattice gauge theory is studied. The Yang-Mills gradient flow plays a crucial role to make the stress tensor well-defined and derivable from the numerical simulations on the lattice. The resultant stress tensor with a decomposition into local principal axes shows, for the first time, the detailed structure of the flux tube along the longitudinal and transverse directions in a gauge invariant manner. The linear confining behavior of the $Q\bar{Q}$ potential at long distances is derived directly from the integral of the local stress tensor., Comment: 7 pages, 4 figures, Fig.3 is corrected, version published in Phys. Lett. B

We study how vortices in dense superfluid hadronic matter can connect to vortices in superfluid quark matter, as in rotating neutron stars, focusing on the extent to which quark-hadron continuity can be maintained. As we show, a singly quantized vortex in three-flavor symmetric hadronic matter can connect smoothly to a singly quantized non-Abelian vortex in three-flavor symmetric quark matter in the color-flavor locked (CFL) phase, without the necessity for boojums appearing at the transition., Comment: 10 pages, 2 figures; discussions on flavor non-singlet vortices elaborated; comments on color flux and center vortices added

We invalidate the arguments given in [T.Yamazaki and Y.Kuramashi, Phys. Rev. D96, 114511 (2017)] over the HAL QCD method for hadron-hadron interactions on the lattice. We also pose questions on the practical usefulness of the method proposed in this reference., Comment: 3 pages. Version accepted for publication in Physical Review D

We present our recent results of baryon interactions with strangeness $S=-1$ based on Nambu-Bethe-Salpeter (NBS) correlation functions calculated from lattice QCD with almost physical quark masses corresponding to $(m_\pi,m_K)\approx(146,525)$ MeV and large volume $(La)^4=(96a)^4\approx$ (8.1 fm)$^4$. In order to perform a comprehensive study of baryon interactions, a large number of NBS correlation functions from NN to $\Xi\Xi$ are calculated simultaneously by using large scale computer resources. In this contribution, we focus on the strangeness $S=-1$ channels of the hyperon interactions by means of HAL QCD method. Four sets of three potentials (the $^3S_1-^3D_1$ central, $^3S_1-^3D_1$ tensor, and the $^1S_0$ central potentials) are presented for the $\Sigma N - \Sigma N$ (the isospin $I=3/2$) diagonal, the $\Lambda N - \Lambda N$ diagonal, the $\Lambda N \rightarrow \Sigma N$ transition, and the $\Sigma N - \Sigma N$ ($I=1/2$) diagonal interactions. Scattering phase shifts for $\Sigma N$ $(I=3/2)$ system are presented., Comment: Talk presented at the 35th International Symposium on Lattice Field Theory (Lattice 2017), 18-24 June 2017, Granada, Spain, 8 pages, 17 figures, references added, typos corrected, arXiv admin note: text overlap with arXiv:1702.00734

We present the latest lattice QCD results for baryon interactions obtained at nearly physical quark masses. $N_f = 2+1$ nonperturbatively ${\cal O}(a)$-improved Wilson quark action with stout smearing and Iwasaki gauge action are employed on the lattice of $(96a)^4 \simeq (8.1\mbox{fm})^4$ with $a^{-1} \simeq 2.3$ GeV, where $m_\pi \simeq 146$ MeV and $m_K \simeq 525$ MeV. In this report, we study the two-nucleon systems and two-$\Xi$ systems in $^1S_0$ channel and $^3S_1$-$^3D_1$ coupled channel, and extract central and tensor interactions by the HAL QCD method. We also present the results for the $N\Omega$ interaction in $^5S_2$ channel which is relevant to the $N\Omega$ pair-momentum correlation in heavy-ion collision experiments., Comment: Talk given at 35th International Symposium on Lattice Field Theory (Lattice 2017), Granada, Spain, 18-24 Jun 2017, 8 pages, 9 figures. arXiv admin note: text overlap with arXiv:1702.01600

The interaction between $\Lambda_c$ and a nucleon ($N$) is investigated by employing the HAL QCD method in the (2+1)-flavor lattice QCD on a $(2.9~\mathrm{fm})^3$ volume at $m_\pi \simeq 410,~570,~700$ MeV. We study the central potential in $^1S_0$ channel as well as central and tensor potentials in $^3S_1 - $$^3D_1$ channel, and find that the tensor potential for $\Lambda_c N$ is negligibly weak and central potentials in both $^1S_0$ and $^3S_1 - $$^3D_1$ channels are almost identical with each other except at short distances. Phase shifts and scattering lengths calculated with these potentials show that the interaction of $\Lambda_c N$ system is attractive and has a similar strength in $^1S_0$ and $^3S_1$ channels at low energies (i.e. the kinetic energy less than about $40$ MeV). While the attractions are not strong enough to form two-body bound states, our results lead to a possibility to form $\Lambda_c$ hypernuclei for sufficiently large atomic numbers ($A$). To demonstrate this, we derive a single-folding potential for $\Lambda_c$ hypernuclei from the $\Lambda_c$-nucleon potential obtained in lattice QCD, and find that $\Lambda_c$ hypernuclei can exist for $A \ge 12$ with the binding energies of a few MeV. We also estimate the Coulomb effect for the $\Lambda_c$ hypernuclei., Comment: 27 pages, 10 figures, 4 tables

Deriving accurate energy density functional is one of the central problems in condensed matter physics, nuclear physics, and quantum chemistry. We propose a novel method to deduce the energy density functional by combining the idea of the functional renormalization group and the Kohn-Sham scheme in density functional theory. The key idea is to solve the renormalization group flow for the effective action decomposed into the mean-field part and the correlation part. Also, we propose a simple practical method to quantify the uncertainty associated with the truncation of the correlation part. By taking the $\varphi^4$ theory in zero dimension as a benchmark, we demonstrate that our method shows extremely fast convergence to the exact result even for the highly strong coupling regime., Comment: 6 pages, 3 figures, and 1 table

The $\Omega\Omega$ system in the $^1S_0$ channel (the most strange dibaryon) is studied on the basis of the (2+1)-flavor lattice QCD simulations with a large volume (8.1 fm)$^3$ and nearly physical pion mass $m_{\pi}\simeq 146$ MeV at a lattice spacing $a\simeq 0.0846$ fm. We show that lattice QCD data analysis by the HAL QCD method leads to the scattering length $a_0 = 4.6 (6)(^{+1.2}_{-0.5}) {\rm fm}$, the effective range $r_{\rm eff} = 1.27 (3)(^{+0.06}_{-0.03}) {\rm fm}$ and the binding energy $B_{\Omega \Omega} = 1.6 (6) (^{+0.7}_{-0.6}) {\rm MeV}$. These results indicate that the $\Omega\Omega$ system has an overall attraction and is located near the unitary regime. Such a system can be best searched experimentally by the pair-momentum correlation in relativistic heavy-ion collisions., Comment: 6 pages and 4 figures

Euclidean two-point correlators of the energy-momentum tensor (EMT) in SU(3) gauge theory on the lattice are studied on the basis of the Yang-Mills gradient flow. The entropy density and the specific heat obtained from the two-point correlators are shown to be in good agreement with those from the one-point functions of EMT. These results constitute a first step toward the first principle simulations of the transport coefficients with the gradient flow., Comment: 5 pages, 12 figures

We review the equation of state of matter in neutron stars from the solid crust through the liquid nuclear matter interior to the quark regime at higher densities. We focus in detail on the question of how quark matter appears in neutron stars, and how it affects the equation of state. After discussing the crust and liquid nuclear matter in the core we briefly review aspects of microscopic quark physics relevant to neutron stars, and quark models of dense matter based on the Nambu--Jona-Lasinio framework, in which gluonic processes are replaced by effective quark interactions. We turn then to describing equations of state useful for interpretation of both electromagnetic and gravitational observations, reviewing the emerging picture of hadron-quark continuity in which hadronic matter turns relatively smoothly, with at most only a weak first order transition, into quark matter with increasing density. We review construction of unified equations of state that interpolate between the reasonably well understood nuclear matter regime at low densities and the quark matter regime at higher densities. The utility of such interpolations is driven by the present inability to calculate the dense matter equation of state in QCD from first principles. As we review, the parameters of effective quark models -- which have direct relevance to the more general structure of the QCD phase diagram of dense and hot matter -- are constrained by neutron star mass and radii measurements, in particular favoring large repulsive density-density and attractive diquark pairing interactions. We describe the structure of neutron stars constructed from the unified equations of states with crossover. Lastly we present the current equations of state -- called "QHC18" for quark-hadron crossover -- in a parametrized form practical for neutron star modeling., Comment: v2, 42 pages, 36 figures, 3 tables; to be published in Reports on Progress in Physics; new sections for cooling, X-ray analyses, and gravitational waves are added; the results for tidal deformability are included; equations of state and the numerical tables are updated; v3, typos corrected in eq.8

On the basis of the $p\Xi^-$ interaction extracted from (2+1)-flavor lattice QCD simulations at the physical point, the momentum correlation of $p$ and $\Xi^-$ produced in relativistic heavy ion collisions is evaluated. $C_{\rm SL}(Q)$ defined by a ratio of the momentum correlations between the systems with different source sizes is shown to be largely enhanced at low momentum due to the strong attraction between $p$ and $\Xi^-$ in the $I=J=0$ channel. Thus, measuring this ratio at RHIC and LHC and its comparison to the theoretical analysis will give a useful constraint on the $p\Xi^-$ interaction., Comment: 4 pages, 2 figures; proceedings contribution for Quark Matter 2017