Arxiv:1910.08133V1 [Nucl-Th] 17 Oct 2019 and the Observed Modiﬁcations of Bound Proton Emffs at Ties Have Been Made Mostly in Vacuum, but Not in Medium

APCTP Pre2019-012, LFTC-19-7/45

Electroweak properties of kaons in a nuclear medium

Parada T. P. Hutauruk1, ∗ and K. Tsushima2, 1, † 1Asia Pacific Center for Theoretical Physics, Pohang, Gyeongbuk 37673, Korea 2Laboratório de Física Teórica e Computacional, Universidade Cruzeiro do Sul / Universidade Cidade de Sao Paulo, 01506-000, São Paolo, SP, Brazil (Dated: October 21, 2019) Kaon electroweak properties in symmetric nuclear matter are studied in the Nambu–Jona-Lasinio model using the proper-time regularization. The valence quark properties in symmetric nuclear matter are calculated in the quark-meson coupling model, and they are used as inputs for studying the in-medium kaon properties in the NJL model. We evaluate the kaon decay constant, kaon-quark coupling constant, and K+ electromagnetic form factor by two different approaches. Namely, by two different ways of calculating the in-medium constituent quark masses of the light quarks. We predict that, in both approaches, the kaon decay constant and kaon-quark coupling constant decrease as nuclear density increases, while the K+ charge radius increases by 20-25% at normal nuclear density.

I. INTRODUCTION following reasons; (i) the constituent quark masses in vacuum are the input parameters in Ref. [32], but it is preferable to Kaons play special roles in strong and electroweak inter- calculate them dynamically, (ii) the vacuum in the light-front action phenomena, and in low-energy hadronic reactions as- approach is generally believed to be "trivial", and there is no sociated with quantum chromodynamics (QCD) at low ener- clearly defined quark chiral condensates in the light-front ap- gies [1–14]. They appear as Nambu-Goldstone bosons like pi- proach. Thus, the model does not have the dynamical chiral ons, but show the flavor SU(3) symmetry breaking by the heav- symmetry breaking mechanism nor does the model have direct ier strange quark mass. Thus, kaons can provide us with useful connection with the emergence of (pseudo-)Goldstone bosons information on chiral symmetry and its restoration, SU(3) fla- such as kaons, and (iii) the in-medium kaon decay constant as vor breaking effect as well as the nonperturbative bound state well as the kaon-quark coupling constant are assumed to be nature of QCD. The internal structure of kaons can be explored density-independent. The purpose of the present work is to by the electromagnetic form factors (EMFFs), which reflect improve further the work of Ref. [32]. We extend the approach the underlying quark-gluon dynamics [15, 16]. Many attempts used for studying the pion properties in symmetric nuclear mat- have been made to understand the internal structure of kaons ter in Ref. [33], and study the in-medium kaon properties using using various approaches, e.g., the Dyson-Schwinger equation the Nambu–Jona-Lasinio (NJL) model [34, 35], supplemented (DSE), Nambu–Jona-Lasinio (NJL) model, chiral perturbation by the QMC model inputs. The NJL model is a powerful chi- theory (ChPT), and instanton vacuum-based models [15–29]. ral effective quark theory of low-energy QCD. Importantly, Experimentally, K-meson EMFFs are poorly known, except the model describes the dynamical chiral symmetry breaking, for the K− in the low Q2 region (Q2 < 0.2 GeV2) [30], where the origin of the pseudoscalar Goldstone bosons such as pions Q2 = −q2 and q is the four-momentum transfer. Recently, K+ and kaons. Furthermore, the model satisfies the chiral limit EMFF was extracted from the kaon electroproduction data at as QCD dictates. Thus, the NJL model, which has several JLab up to Q2 = 2.07 GeV2 [31], but the data uncertainties are improved aspects as addressed above, is suitable to study the still large. We await the new coming data with better precision in-medium kaon properties. up to Q2 ' 5.5 GeV2. The experimental evidences, such as the EMC effect [36, 37] On the other hand, theoretical studies for K-meson proper- arXiv:1910.08133v1 [nucl-th] 17 Oct 2019 and the observed modifications of bound proton EMFFs at ties have been made mostly in vacuum, but not in medium. JLab [38], suggest that the internal structure of hadrons would For example, kaon properties and K+ EMFF in symmetric nu- be modified in a nuclear medium. The phenomena of in- clear matter were studied only recently in Ref. [32] based on medium modifications of hadron properties [39–50] are tightly a light-front constituent quark model. There, the in-medium connected with partial restoration of chiral symmetry [51–54]. valence light-quark properties were calculated in the quark- The order parameters of chiral symmetry in QCD are the light- meson coupling (QMC) model, and they were used as inputs quark chiral condensates, and their changes are expected to be to study the in-medium kaon properties. Although this study one of the most important driving forces for the change of provides us with some insights on the kaon properties in a hadron properties in a nuclear medium. nuclear medium, more complete studies are necessary for the Spontaneous breaking of chiral symmetry generates the nonet of massless pseudoscalar Goldstone bosons. But the explicit breaking of U(1) axial symmetry selectively shifts up ∗ [email protected] the η0-singlet mass, leaving the SU(3) flavor octet of pions, † [email protected] K-mesons, and η8 to be massless. Then the explicit chiral 2 symmetry breaking by non-vanishing current quark masses erties in the QMC model and then calculate the in-medium leads to the experimentally observed low-lying pseudoscalar dynamical (constituent) quark masses in the NJL model using meson mass spectrum [55, 56]. Since the chiral symmetry the QMC model inputs. This approach is called the NJL-based has such a big impact on the low-lying pseudoscalar meson approach. A similar approach was adopted in the study of the mass spectrum with its explicit breaking, partial restoration of in-medium pion properties in Ref. [33]. The second approach chiral symmetry in a strongly interacting medium is inevitable is called the QMC-based approach, uses the constituent quark to study the change of hadron properties in a nuclear medium. mass values in vacuum determined in the NJL model, and the Studying of kaon property changes in medium allows us to values are plugged in the QMC model to calculate the density explore the nature of the bound state with the light and strange dependence of the constituent quark masses. quarks in different environment from vacuum. In the present work, we study the spacelike K+-meson EMFF, kaon leptonic (weak) decay constant, and kaon-quark A. In-medium quark properties in the QMC model coupling constant in symmetric nuclear matter using the NJL model. The NJL model has been very successful in studying We briefly review the in-medium properties of the valence pseudoscalar meson properties [55–57]. Recently, the model light and strange quarks in the QMC model [61]. The QMC was applied to study the K+ EMFF in vacuum [25, 28, 58], model has been successfully applied for many topics of nuclear kaon valence quark distributions in medium [59], and in- and hadronic physics. (See, for example, Refs. [43, 62], and medium π+ EMFF [33]. In these works the in-medium valence references therein.) It has also been used to study the medium light-quark properties were calculated in the QMC model, and modifications of the nucleon weak and electromagnetic form they were used as inputs for studying the in-medium pion and factors for the neutrino scattering in dense matter [63, 64]. kaon properties in the NJL model. Extending the works of Recently, the model has been further applied for studying the Refs. [28, 33], we study here the electroweak properties of in-medium properties of low-lying strange, charm, and bottom kaons in symmetric nuclear matter using the NJL model. baryons [65]. For studying the in-medium kaon properties microscop- In the QMC model medium effects arise from self-consistent ically, we need the valence quark properties in a nuclear exchange of the Lorentz-scalar (σ) and Lorentz-vector (ω, ρ) medium. These are calculated by the QMC model, and used meson fields that couple directly only to the light quarks con- as inputs for studying the in-medium kaons properties. For fined in hadrons. The physics behind this is that the light-quark this purpose, we adopt two approaches in this work. The first chiral condensates in medium are more sensitively change than approach is that, the in-medium constituent quark masses are those of the strange and heavier quarks. Here, we work with calculated by the QMC model [43], which will be denoted symmetric nuclear matter in its rest frame in the Hartree mean by “QMC-based approach”. The second approach is that, the field approximation. (For the Fock term effects in the QMC in-medium constituent quark masses are calculated by solving model, see Ref. [66].) the NJL model gap equations [60] using the in-medium “current quark properties” calculated by the QMC model as inputs. Effective Lagrangian for symmetric nuclear matter in the This will be denoted by “NJL-based approach”. Although the QMC model is given by [67, 68] two approaches yield different in-medium constituent quark = ψ¯ iγ · ∂ − M∗ (σ) − g ωµγ ψ + , masses of the light quarks, they lead to very similar predic- LQMC N N ω µ N Lmeson tions for the K+ EMFF and its charge radius in symmetric (1) nuclear matter. where ψN , σ, and ω are respectively the nucleon, scalar meson This paper is organized as follows. In Sec. II, we present σ, and vector meson ω fields. Note that the isospin-dependent the in-medium constituent valence-quark and kaon properties ρ-meson field is absent, since we use the Hartree approxima- based on the two approaches, the QMC-based and NJL-based + tion for symmetric nuclear matter (isospin saturated), and the approaches. In section III we present how the in-medium K ∗ ρ mean field vanishes. The effective nucleon mass MN (σ) is EMFF in the NJL model is calculated, while the numerical defined by results are given in Sec. IV. We summarize and conclude in section V. ∗ MN (σ) = MN − gσ(σ)σ. (2)

Here, gσ(σ) and gω are the σ-dependent nucleon-σ coupling II. IN-MEDIUM KAON PROPERTIES strength and the nucleon−ω coupling constant, respectively. N We deﬁne gσ ≡ gσ(σ = 0) for later convenience. The In this section we discuss the details of the quark and kaon mesonic Lagrangian density Lmeson in Eq. (1) is given by properties in symmetric nuclear matter with both the QMC- 1 1 + µ 2 2 µ ν ν µ based and NJL-based approaches. To study the K EMFF in Lmeson = (∂µσ∂ σ − mσ σ ) − ∂µων(∂ ω − ∂ ω ) 2 2 the NJL model, we need the in-medium valence quark prop- 1 erties as inputs calculated in the QMC model. There are two + m2 ωµω . (3) 2 ω µ ways of adapting the inputs. It is concerned about how the in-medium constituent quark masses are calculated and used. In the Hartree approximation for symmetric nuclear matter, Namely, we ﬁrst evaluate the in-medium “current quark” prop- the nucleon Fermi momentum kF is related with the baryon 3

(ρB) and scalar (ρs) densities respectively, will be shown to be Lorentz scalar, is calculated as ∗ ¹ 3 nj Ω − zh 4 4kF ∗ Õ j 4 ∗3 ρ = dk θ (kF − |k|) = , m = + πR B, (10) B (2π)3 3π2 h R∗ 3 h j=l,l¯,s,s¯ h ¹ ∗ 4 MN (σ) ρ = dk θ (k − |k|) . (4) ∗ s 3 F q and the in-medium bag radius R is determined by the stability (2π) ∗2( ) 2 h MN σ + k condition for a given baryon density self-consistently,

∗ dmh ∗ = 0. (11) In the QMC model [43, 67, 68], nuclear matter is treated as dRh the collection of nonoverlapping nucleon-MIT bags [69]. The Dirac equations for the quarks and antiquarks in a hadron bag In Eq. (10) zh is assumed to be density independent. It is re- are given by lated with the sum of the center-of-mass and gluon fluctuation corrections [67], which is determined by the hadron mass in l 0 l 1 l ψu(x) vacuum. The bag constant B is fixed by the inputs for the nu- iγ · ∂x − ml − Vσ ∓ γ Vω + Vρ = 0, 2 ψu¯ (x) cleon in vacuum, namely, RN = 0.8 fm with MN = 939 MeV, for a chosen value of ml. l 0 l 1 l ψd(x) iγ · ∂x − ml − Vσ ∓ γ Vω − Vρ = 0, The ground state wave function of a quark in hadron h 2 ψd¯(x) satisfies the boundary condition at the bag surface, ψ (x) [iγ · ∂ − m ] s = 0, (5) x s ψ (x) s¯ j0(x) = βq j1(x) (12) where l stands for “light”, i.e., l = u or d, and the effective with q = l (= u, d) or s, and j and j are the spherical Bessel ∗ 0 1 light-quark mass ml is defined as functions, and

m∗ ≡ m − Vl . (6) sΩ∗ − ∗ ∗ l l σ q mq Rh βq = ∗ ∗ ∗ . (13) Ωq + mq R Here we assume SU(2) symmetry, ml = mu = md for the h valence quark masses throughout this study (thus m∗ = m∗ = l u The scalar σ and vector ω meson fields at the nucleon level m∗ ), and m is the strange valence quark mass in vacuum, and d s are related as −Vl is the Lorentz scalar potential, which couples only to the σ g ρ light quarks. The strange quark is decoupled from the scalar ω = ω B , 2 and vector potentials so that its effective mass does not change mω ∗ ∗ in a nuclear medium, ms = ms. 4gN C (σ) ¹ M (σ) σ = σ N dk θ(k − |k|) N , (14) 3 2 F q (2π) mσ ∗2( ) 2 The scalar and vector mean fields in symmetric nuclear mat- MN σ + k ter are defined, respectively by the mean expectation values, where CN (σ) is defined by Vl ≡ gl σ = gl hσi , Vl ≡ gl ω = gl δµ,0 hωµi , (7) σ σ σ ω ω ω ∗ −1 ∂MN (σ) where the light-quark-meson coupling constants, gl and gl , CN (σ) = , (15) σ ω gN ∂σ are defined later through Eq. (17). The eigenenergies in units σ of 1/R∗ with the bag radius of the hadron h, are given by h with CN (σ) = 1 for the pointlike nucleon case such as in Quantum Hadrodynamics (QHD) [70, 71]. The σ-dependent 1 u ∗ ∗ l l coupling gσ(σ) or CN (σ) is the origin of the novel satura- = Ωl ± Rh Vω + Vρ , u 2 tion properties in the QMC model, and the quark dynamics is included in the effective nucleon mass M∗ (σ) via a self- d ∗ ∗ l 1 l N = Ω ± R Vω − Vρ , ∗ → ∗ l h 2 consistent manner as in Eqs. (10) [with mh MN ] and (14). d By solving the self-consistent equation for the σ mean field Ω∗. s = s = s (8) Eq. (14), the total energy per nucleon is calculated by where 4 ¹ q Etot/A = dk θ(k − |k|) M∗2(σ) + k2 3 F N h i1/2 (2π) ρB Ω∗ = Ω∗ = x2 + R∗ m∗2 , l l l h l m2 σ2 g2 ρ + σ + ω B . (16) m∗ = m − Vl = m − gl σ, 2 l l σ l σ 2ρB 2mω Ω∗ = Ω∗ = [x2 + (R∗ m )2]1/2. (9) s s s h s N N The coupling constants gσ (= gσ(σ = 0)) and gω = gω ∗ N The effective mass of hadron h in nuclear medium mh, which in Eq. (16) (gσ is implicit) are determined to reproduce the 4 binding energy of symmetric nuclear matter 15.7 MeV at the given by saturation density ρ = 0.15 fm−3, and they are, respectively, 0 ∗ ∗ N l ( ) l k/ + M gσ = 3gσ SN σ = 0 and gω = 3gω. In this study, we will use S∗(k∗) = l , l ( ∗)2 − ( ∗)2 ml = 16.4 and 400 MeV, and thus we give the explicit values of k Ml + i the coupling constants for the both quark mass values. Namely, k/ + Ms gl ' 5.6251 [4.7009] for m = 16.4 [400] MeV, where S (σ) S∗(k∗) = S (k) = , (19) σ l N s s 2 2 is defined through [44, 61] k − Ms + i dM∗ (σ) ¹ where the quantities with the asterisk above denote those in N = −3gl d3 y ψ¯ (y) ψ (y) medium, and medium effects enter in the light-quark mass M∗ dσ σ l l l bag and momentum kµ by k∗µ = kµ + V µ due to the vector mean µ 0 l d fieldV = (V , 0). The modification of the space component of ≡ −3g SN (σ) = − g (σ)σ , (17) σ dσ σ the light-quark momentum k∗µ is neglected, since it is known very small [66]. which gives SN (σ = 0) = 0.4899 [0.6950] for ml = . [ ] R . ψ 16 4 400 MeV and N = 0 8 fm with l being the lowest In Eq. (18) the vector potential entering in the light-quark mode bag wave function in medium. Note that, the right hand propagator was eliminated by the shift of the integral vari- side of Eq. (17) is the quark scalar charge, which is Lorentz able [74]. Note that, the in-medium strange quark propagator scalar, and thus the left-hand-side of Eq. (17) is Lorentz scalar, is the same as that in vacuum, since the strange quark is decou- and thus M∗ (σ) as well. These relations show that the in- N pled from the scalar and vector mean fields in the QMC model, medium quark dynamics is explicitly included in the QMC ∗ and the same assumption is adopted. Namely, Ms = Ms is as- model. Since the light quarks in any hadrons should generally sumed also in the present NJL model calculation. feel the same scalar and vector mean fields as those in the nucleon (after one chooses a fixed light-quark mass value in The description of kaon as the dressed quark-antiquark vacuum), we can systematically study the hadron properties in bound state is obtained by solving the Bethe-Salpeter equa- medium without introducing any new coupling constants for tion (BSE) in the random phase approximation (RPA). The the σ and ω mean fields for different hadrons. solution to the BSE in each meson channel is given by the two-body scattering t-matrix, that depends on the interaction channel. We also introduce below the reduced t-matrix for vector mesons, since it will be used later in Eqs. (33) and (34). The reduced t-matrices in the K-meson and vector meson V take the forms [28], B. Kaon properties in the NJL model

−2i G τ∗ (p∗) = π , Based on the SU(3) NJL model in vacuum described in K ∗ ∗2 1 + 2 Gπ ΠK (p ) Ref. [28], we explore the in-medium properties of dynamical ∗µ ∗ν −2i Gρ p p quark and kaon properties in this section. The medium effect is τ∗µν(p∗) = gµν + 2 G Π∗ (p∗2) , V Π∗ ( ∗2) ρ V ∗2 implemented through the in-medium potentials for the valence 1 + 2 Gρ V p p quarks generated by the QMC model. The dynamical quark (20) ∗ mass in medium, Mq (q = u, d, s), in the NJL model with the proper-time regularization scheme is given by [33] where Gρ is the four-fermion coupling constant for ρ channel, and the bubble diagrams in nuclear medium give ∗ ∞ 3Gπ Mq ¹ dτ ∗2 M∗ = m∗ + e−τMq , (18) ¹ d4k q q 2 2 ∗ ∗2 ∗ ∗ ∗ ∗ ∗ π 1/Λ2 τ Π (p ) = 6i TrD γ5S (k )γ5S s(k + p ) , UV K (2π)4 l ∗ ∗ ¹ 4 where Mq and mq are the in-medium dynamical and current ∗qq ∗2 ∗µν d k µ ∗ ∗ ν ∗ ∗ ∗ ΠV (p ) PT = 6i TrD γ S q(k )γ S q(k + p ) , quark masses, respectively. Gπ is the four-fermion coupling (2π)4 constant, which is taken the same as in vacuum, and ΛUV is the (21) ultraviolet cutoff. For the infrared cutoff ΛIR, we take τIR = /Λ2 ∞ Λ ∗ ∗ ∗ll ∗ ∗ss 1 IR = + in nuclear medium. (This change in IR from with Πρ = Πω = Πv , and Πφ = Πv . The trace is only for that in Ref. [28] gives negligible difference in results.) Note ∗µν µν ∗µ ∗ν ∗2 ∗ the Dirac indices, and PT = g − p p /p . that the expression of Mq is different from that in Refs. [72, 73], since the information of “nucleon” Fermi momentum and The in-medium meson mass is defined by the pole in the ∗ nuclear matter saturation properties are included in mq (and corresponding t-matrix as in the vacuum case, Vl ) calculated in the QMC model, and we should drop the ω ∗ ∗2 ∗2 density dependent term. 1 + 2 Gπ ΠK (p = mK ) = 0, (22) The in-medium dressed quark propagators in Eq. (18) are where the similar conditions determine other meson masses in 5 medium. This leads to the in-medium kaon mass, TABLE I. Parameters in the QMC model and some quantities obtained −3 m∗ m∗ 1 at saturation density ρ0 = 0.15 fm for the two approaches, the NJL- m∗2 = s + l + (M∗ − M∗)2, (23) based and QMC-based. They respectively correspond to the quark K M∗ M∗ I ( ∗2) s l s l Gπ ls mK masses in vacuum, ml = 16.4 (NJL-based), and 400 MeV (QMC- ∗ based). The light-quark mass ml, effective nucleon mass MN , nuclear ∗ ∗ ∗ ∗2 1/4 where ml = mu = md and Ils(p ) is defined by, incompressibility K, and B are given in units of MeV. The other parameters are obtained by using the inputs for the nucleon bag radius 3 ¹ 1 ¹ dτ R = 0.8 fm in vacuum, nucleon mass in vacuum M = 939 MeV, I (p∗2) = dz N N ab 2 and the nuclear matter binding energy of 15.7 MeV at the saturation π 0 τ density ρ . −τ[−z(1−z) p∗2+2pV z(1−z)−z(1−z)V 2+x M∗2+(1−z) M∗2] 0 × e 0 0 b a . 2 2 1/4 ∗ (24) ml gσ/4π gω/4π B zN MN K 16.4 5.438 5.412 169.2 3.334 751.95 281.50 This demonstrates the Goldstone boson nature of the kaon in 400 7.645 10.136 98.5 5.321 625.98 339.10 the chiral limit. The residue at the pole in the qq¯ t-matrices ∗ defines the kaon-quark coupling constant (gKqq) in medium: 400

∂ Π∗ (p∗2) (g∗ )−2 = − K . (25) Kqq ∂p∗2 ∗2 ∗2 p =mK 300 (MeV)

Following Refs. [28, 58, 63], we have chosen ΛIR = 240 ∗ q MeV, which is of the order of ΛQCD, and Ml = 400 MeV in m or m vacuum. These values were ﬁxed to give = 140 MeV and ∗ q π 200 mK = 495 MeV, together with the pion decay constant fπ = M 93 MeV. Furthermore, this gives ΛUV = 645 MeV, Gπ = 19.0 Ml∗, ml = 16.4 MeV −2 GeV . Fit to the physical masses of the vector mesons, mρ = ml∗, ml = 400 MeV −2 770 MeV and mω = 782 MeV, gives Gρ = 11.0 GeV , and 100 −2 0 0.2 0.4 0.6 0.8 1 1.2 Gω = 10.4 GeV , respectively. These parameters are used in 2 ρB/ρ0 the present study, except for (1/ΛIR) → ∞ in medium. FIG. 1. In-medium constituent quark mass M∗ in the NJL-based l∗ approach (solid line) with ml = 16.4 MeV, and m in the QMC-based C. QMC-based and NJL-based approaches l approach (dashed line) with ml = 400 MeV.

As mentioned already, we adopt two approaches to estimate the in-medium constituent quark masses. In the NJL-based In Fig. 1 we compare the in-medium constituent quark approach we use the in-medium “current quark” properties masses calculated in the two approaches. The solid line is ∗ evaluated by the QMC model, and calculate the in-medium dy- the Ml (ml = 16.4 MeV) in the NJL-based approach, while ∗ namical (constituent) quark masses in the NJL model. While the dashed line is the ml (ml = 400 MeV) in the QMC-based in the QMC-based approach, we use the vacuum constituent approach. This shows the two approaches give similar density quark mass values fixed by the fit in the NJL model, and the dependence for the constituent quark mass of the light-quark. QMC model calculates the density dependence of the con- The largest difference between the two approaches is about stituent quark masses. Note that, the in-medium kaon proper- 10%. Relative to the vacuum value, the light-quark constituent ∗ ties including effective kaon mass mK and quark EMFFs, are mass is reduced by 30–35% at nuclear matter saturation density all calculated in the NJL model for both approaches. ρ0. The negative of binding energy per nucleon (Etot/A− MN ) We first discuss the parameters of the NJL-based approach. for symmetric nuclear matter is also shown in Fig. 2 for the two approaches. Following Ref. [33], we adopt ml = mu = md = 16.4 MeV and ms = 356 MeV, for the vacuum current quark mass val- Results for the density dependence of the light-quark chiral ues. These values are optimally fixed by the NJL model in condensate are given in Table II in the two approaches. We find Refs. [16, 28]. The other parameters in the QMC model are that the in-medium light-quark chiral condensate decreases fixed by the nucleon mass MN = 939 MeV, nucleon bag ra- about 13–16% relative to that in vacuum at normal nuclear dius RN = 0.8 fm, and the nuclear matter binding energy density ρ0, and the difference in the two approaches is small. per nucleon of 15.7 MeV at the saturation density ρ0 = 0.15 In Tables III and IV we list the density dependence of the fm−3. For the QMC-based approach, the NJL-model vacuum kaon properties including the kaon-quark coupling constant constituent quark mass values are used as inputs to calculate calculated in the two approaches. We find that the in-medium the density dependence of the constituent quark masses in the light-quark dynamical quark mass, effective kaon mass, and ∗ QMC model. In this case we use Ml = Mu = Md = 400 MeV kaon-quark coupling constant gKqq decrease in medium as and Ms = 611 MeV. They are summarized in Table I. nuclear density increases. The in-medium kaon leptonic de- 6

60 ml = 16.4 MeV TABLE IV. Same as Table III, but in the QMC-based approach (ml = 400 MeV). ml = 400 MeV ρ /ρ m∗ m∗ f ∗ g∗ (MeV) 40 B 0 l K K Kqq N 0.00 0.400 0.495 0.091 4.570 M

− 20 0.25 0.357 0.452 0.091 4.520 A)

/ 0.50 0.316 0.418 0.090 4.467 tot 0 0.75 0.278 0.390 0.088 4.419 (E 1.00 0.243 0.371 0.087 4.376 20 1.25 0.212 0.360 0.085 4.339 − 0 0.5 1 1.5 2 2.5 3

ρB/ρ0 1000

FIG. 2. Negative of binding energy per nucleon (Etot/A − MN ) for symmetric nuclear matter calculated in the QMC model for the light- 800 quark mass ml = 16.4 MeV (solid line) and that for ml = 400 MeV (dotted line). 600 (MeV) ∗ N 1/3 400 TABLE II. In-medium light-quark chiral condensates, − hl li = M − huui1/3 (GeV), calculated in the two approaches. 200 ml = 16.4 MeV m = 400 MeV ρB/ρ0 NJL-based approach QMC-based approach l 0.00 0.171 0.171 0 0 0.5 1 1.5 2 2.5 3 0.25 0.167 0.165 ρB/ρ0 0.50 0.162 0.158 ∗ 0.75 0.156 0.151 FIG. 3. Eﬀective nucleon mass MN for symmetric nuclear matter cal- 1.00 0.149 0.144 culated in the QMC model, corresponding to the NJL-based approach (solid line), and for the QMC-based approach (dotted line). 1.25 0.136 0.137

∗ III. K+ ELECTROMAGNETIC FORM FACTOR IN cay constant fK decreases as nuclear density increases, but the reduction rate is smaller than that of pion [33], which is MEDIUM consistent with the result in Ref. [75]. In this section, we present the calculation of the K+ EMFF in symmetric nuclear matter following Ref. [28]. Here, the TABLE III. Constituent quark mass of the light quark and kaon prop- in-medium light-quark propagators appearing in the Bethe- erties in symmetric nuclear matter calculated in the NJL-based ap- Salpeter equation are modiﬁed (see Eq. (19)). proach (ml = 16.4 MeV). The masses and decay constant are given ∗ in units of GeV while the coupling constant gKqq is dimensionless. The electromagnetic interaction with the quark is ob- µ tained by the minimal substitution: i∂/ → i∂/ − QAˆ µ γ , ρ /ρ M∗ m∗ f ∗ g∗ B 0 l K K Kqq where Aµ is photon ﬁeld, e is the positron charge and e 1 0.00 0.400 0.495 0.091 4.570 Qˆ = diag [eu, ed, es] = (λ3 + √ λ8) is the quark charge 2 3 0.25 0.370 0.465 0.091 4.536 operator with eu,d,s being the respective quark charges. 0.50 0.339 0.437 0.090 4.495 The matrix element of the in-medium electromagnetic cur- 0.75 0.307 0.411 0.089 4.455 rent for the K+-meson (will be denoted simply as kaon here- 1.00 0.270 0.386 0.088 4.408 after, otherwise stated) is written as 1.25 0.207 0.359 0.084 4.332 µ 0∗ ∗ 0∗µ ∗µ ∗ 2 J (p , p ) = (p + p ) FK (Q ), (26)

∗ 0∗ Illustrated in Fig. 3 are the effective nucleon mass M∗ calcu- where p and p are the initial and final four-momenta of the N 2 ( 0∗ − ∗)2 ≡ − 2 ∗ ( 2) lated in the QMC model, corresponding to the two approaches, kaon, respectively, with q = p p Q and FK Q the NJL-based and the QMC-based. The corresponding effec- denotes the in-medium kaon EMFF. ∗ l l tive quark mass ml , scalar −Vσ, and vector Vω potentials are As the case in vacuum, the in-medium kaon EMFF in the shown in Fig. 4. NJL model is calculated by the sum of the two Feynman dia- 7

400 400

0 0

400 400 − −

m∗ ml = 16.4 MeV m∗ ml = 400 MeV and potentials (MeV) l and potentials (MeV) l 800 l 800 l ∗ l − -Vσ ∗ l − -Vσ m l m l Vω Vω 1200 1200 − 0 0.5 1 1.5 2 2.5 3 − 0 0.5 1 1.5 2 2.5 3

ρB/ρ0 ρB/ρ0

∗ − l l FIG. 4. Eﬀective light-quark mass (ml ) and quark potentials ( Vσ and Vω) used as inputs in the NJL-based approach with ml = 16.4 MeV (left panel), and those in the QMC-based approach with ml = 400 MeV (right panel), respectively. These are calculated in the QMC model.

µ p q ′ p′ p′ k p + k p′ + k = + + p p p p p′ p p′ k p k p′ p′ p′ k − q − µ = + p p FIG. 5. Diagrammatic representation of the electromagnetic interaction for the kaon. FIG. 6. Illustration of the inhomogeneous BSE which gives the dressed quark-photon vertex. The large shaded oval represents the solution to the inhomogeneous BSE, the small dot is the inhomogeneous driving term (Qˆ γµ), and the double-dots represents the qq¯ interaction kernel [28].

indicates the spectator, and explicit expression is given by grams depicted in Fig. 5, which give [28] 3g∗2 ¹ 1 ¹ ¹ d4k ∗ab 2 Kqq dτ ∗2 2 µ 0∗ ∗ ∗2 fα (Q ) = dx exp{−τ[Ma + x(1 − x) Q ]} j (p , p ) = igKqq 2 τ 1,K (2π)4 4π 0 3g∗2 ¹ 1 ¹ 1−x ¹ × † ∗( 0∗ ∗) ˆ µ ∗( ∗ ∗) ∗( ∗) Kqq Tr γ5 λα Sl p + k Q γ Sl p + k γ5 λα Ss k , (27) + dx dz dτ 2 ¹ 4 4π 0 0 µ d k 0∗ ∗ ∗2 ∗ ∗ ∗ ∗ ∗ ∗ j (p , p ) = igKqq ×[( ) 2 ( − )2( ) ( − )] 2,K (2π)4 x + z mK + Ma Mb x + z + 2 Mb Ma Mb ∗ ∗ ∗ ˆ µ ∗ ∗ 0∗ † ∗ ∗ × exp{−τ[(x + z)(x + z − 1) m∗2 + (x + z) M∗2 × Tr γ5 λα Ss (k − p ) Q γ Ss (k − p ) γ5 λα Sl (k ) , (28) K a +(1 − x − z) M∗2 + x z Q2]}. (30) where the trace is over the Dirac, color, and flavor indices. The b index α labels the state and λα are the corresponding flavor These results are denoted as “bare”, because the quark-photon matrices. In flavor space the in-medium quark propagator is vertex is elementary, i.e., Λµ(bare) = Qˆ γµ, with Qˆ being the ∗ ( ∗) [ ∗ ( ∗) ∗ ( ∗) ∗( ) ( )] γq defined by Sq k = diag Su k , Sd k , Ss k = Ss k (see quark-charge operator. We note that these expressions satisfy Eq. (19)). charge conservation exactly. The quark sector EMFFs in medium for a hadron α can be defined similarly to those in vacuum [28]:

We begin with the relation between the in-medium quark ∗ 2 ∗u 2 ∗d 2 ∗s 2 bare EMFFs and kaon EMFF: Fα(Q ) = eu Fα (Q ) + ed Fα (Q ) + es Fα (Q ) + ··· . (31) * (bare) 2 ∗ls 2 ∗sl 2 F + (Q ) = eu f (Q ) − es f (Q ). (29) K K K ∗u 2 The “bare” quark form factor in medium Fα (Q ) in a pseu- ∗ab 2 doscalar meson can easily be related with Eq. (29). In the expression fK (Q ) the first superscript a indicates the struck quark by the photon, and the second superscript b Generally, the quark-photon vertex is not elementary (Qˆ γµ) 8 but dressed. With the dressing, the quark-photon vertex is are calculated using the density dependent effective light-quark given by the inhomogeneous Bethe-Salpeter equation as il- masses determined by the QMC model with the NJL-model lustrated in Fig. 6. For denoting the “dressing”, we replace fixed vacuum values, ml = 400 MeV and ms = 611 MeV (See by (q, u, d, s) → (Q˜, U, D, S) in the rest of this section. With Table. IV). the NJL-model interaction kernel, the solution to the dressed The in-medium K+ EMFF along with the K+ vacuum form quark-photon vertex for a flavor Q˜ quark has the form: factor and its u and s quark form factors are shown in Fig. 7(a)- 2 ∗µ ∗ (d). Figure 7(a) shows the Q dependence of the in-medium p /p + Λ∗µ (p0∗, p∗) = e γµ + γµ − F∗ (Q2) total K EMFF as well as the in-medium quark EMFFs for a γ Q˜ q q∗2 Q˜ few nuclear densities, ρB/ρ0 = 0.00 (solid line), 0.50 (dash- µ ∗ 2 dotted line), and 1.00 (dotted line), where ρ and ρ are re- → γ F ˜ (Q ). (32) B 0 1Q spectively the nuclear density and the normal nuclear density, ∗µ ∗ ∗2 −3 To arrive at the final result, we have used that the p /p /p ρ0 = 0.15 fm . term does not contribute to the electromagnetic current For the vacuum case, we find an excellent agreement with by the current conservation. Note that this form clearly the existing data points in Ref. [30] as well as the empirical ∗µ 0∗ ∗ 2 2 2 −1 satisfies the Ward-Takahashi identity q Λ (p , p ) = monopole function of FK (Q ) = [1 + Q /ΛK ] [30] with µ γ Q˜ Λ2 = 0.687 GeV2. We find the total K+ EMFF, as well as h ∗−1 0∗ ∗−1 ∗ i K e ˜ S (p ) − S (p ) . + Q Q˜ Q˜ the u-quark EMFF in the K -meson in medium (shown in For the dressed l = u, d, and s quarks, we find Fig. 7(b)), decrease faster than those corresponding in vacuum with Q2 as nuclear density increases. As in Figs. 7(b) and 7(c), e F∗ (Q2) = l , (33) we find that the in-medium strange quark EMFF begins to 1 l=U,D ∗ll 2 + 2 2 1 + 2Gρ Πv (Q ) dominate the K EMFF for Q ≥ 1.6 GeV compared to the e u-quark EMFF, and this becomes more pronounced at larger F∗ (Q2) = s , (34) 1 S ∗ss 2 2 1 + 2Gρ Πv (Q ) Q and higher density ρB. The medium effect can clearly be seen, when the K+ EMFF is multiplied by the Q2 as shown in where the explicit form of the in-medium bubble diagram is Fig. 7(d). 2 ¹ 1 ¹ Next, we discuss the results obtained by the NJL-based ∗qq 2 3Q dτ + Πv (Q ) = dx x(1 − x) approach. Results for the in-medium K EMFF are calculated π2 τ ∗ 0 with the NJL model dynamical quark masses in medium, Mu ∗2 2 ∗ × exp −τ Mq + x (1 − x) Q . (35) and Ms = Ms, using the inputs from the QMC model obtained ∗ by ml = 16.4 MeV and ms = 356 MeV in vacuum as well as ml As in the vacuum case [28], we neglect the quark flavor in medium. The kaon EMFF results are shown in Fig. 8(a)-(d). mixing for the in-medium dressed quark form factors, as well 1 Figure 8(a)-(d) show similar features with those shown in as the in-medium dressed quark masses. Fig. 7(a)-(d) for both the in-medium K+ EMFF, in-medium u K+ + 2 2 The final expression for the in-medium EMFF with a and s quark EMFFs in the K , and Q FK (Q ), respectively. dressed quark-photon vertex is given by The result of the NJL-based approach shows the similar trend as that of the QMC-based approach. This can be understood, ∗ ( 2) ∗ ( 2) ∗ls( 2) − ∗ ( 2) ∗sl( 2) FK+ Q = F1U Q fK Q F1S Q fK Q , (36) since the values of the effective quark masses, effective kaon mass, and in-medium kaon-quark coupling constant are nearly where the in-medium quark EMFFs are obtained by Eqs. (31)- the same between the two approaches. Based on this fact, (34). (See also Eq. (29)). the similar features can also be expected for the other mesons by the two approaches adopted in this study. For large Q2, the in-medium K+ EMFF for a few nuclear densities in both IV. NUMERICAL RESULTS approaches have a good agreement with the QCD result in the conformal limit [76] as the vacuum case [28]. We present our results for the K+-meson EMFF in the two In Fig. 9 we show the total K+ EMFF calculated by the approaches: the QMC-based and the NJL-based. First, we dis- QMC-based and NJL-based approaches for a few nuclear den- cuss the results obtained by the QMC-based approach. Those sities, compared with the existing data in vacuum. Results for the K+ charge radius2 for several nuclear densities calculated in the QMC-based and NJL-based approaches are listed in Tables V and VI, respectively. We find that the 1 In the limit Q2 → ∞ these form factors reduce to the elementary quark charge radius in vacuum has an excellent agreement with the charges as expected from asymptotic freedom in QCD. For small Q2 these empirical data in vacuum [77]. Table V clearly shows that results are similar to the expectations from vector meson dominance, where the u and d quarks are dressed by ρ and ω mesons and the s quark by the φ meson. Note, the denominators in Eqs. (33) and (34) are the same as the pole condition obtained by solving the Bethe-Salpeter equation in the ρ, ω or φ channels. Therefore, the dressed u and d quark form factors have 2 The charge radius of the K+ meson is calculated using the relation: 2 2 2 poles at Q = −mρ = −mω , while the dressed s quark form factor has a r 2 ∂FK (Q ) 2 2 hrK i = −6 2 | 2 . pole at Q = −mφ . ∂Q Q =0 9

1 1 (a) (b) [ρB = 0.00ρ0] 0.8 [ρB = 0.50ρ0] 0.8 [ρB = 1.00ρ0] ) ) 2 0.6 2 0.6 Q Q ( ( u K 0.4 K 0.4 F F

0.2 0.2

0 0 0 1 2 3 4 0 1 2 3 4 Q2 (GeV2) Q2 (GeV2)

1 1 (c) (d)

0.8 0.8 ) ) 2 2

0.6 Q Q 0.6 ( ( ¯ K s K F

F 0.4

2 0.4 − Q 0.2 0.2

0 0 0 1 2 3 4 0 1 2 3 4 Q2 (GeV2) Q2 (GeV2)

+ 2 2 FIG. 7. K EMFF in symmetric nuclear matter calculated by the QMC-based approach, (a) Total kaon EMFF, FK (Q ) as a function of Q for u ( 2) 2 s¯ ( 2) 2 a few densities, (b) Up quark EMFF, FK Q as a function of Q for a few densities, (c) Strange quark EMFF, -FK Q as a function of Q for 2 ( 2) 2 / a few densities, and (d) Q FK Q as a function of Q for a few densities. The lines are for ρB ρ0 = 0.00 (solid line), 0.50 (dash-dotted line), and 1.00 (dotted line), respectively. the values of the u-quark charge radius as well as the kaon TABLE VI. Same as Table V, but calculated in the NJL-based ap- charge radius increase as nuclear density increases. At normal proach using the in-medium inputs from the QMC model generated nuclear density, the kaon and the u-quark charge radii increase for mu = 16.4 MeV, and ms = 356 MeV. These mu and ms yield respectively about 20% and 25% relative to those correspond- the NJL model in-vacuum dynamical (constituent) quark masses ing in vacuum. For the s-quark, the value of the charge radius Mu = 400 and Ms = 611 MeV, respectively. Using the QMC ∗ model generated mu from mu = 16.4 MeV, the NJL model calcu- ∗ lates the density dependent dynamical quark mass Mu, and uses this ∗ ∗ Mu together with Ms = Ms. TABLE V. K+ charge radius and quark charge radii calculated in expt the QMC-based approach. The results are calculated using the in- ρB/ρ0 rK ru rs r medium inputs from the QMC model generated for mu = 400 MeV ∗ 0.00 0.59 0.65 0.44 0.56 ± 0.03 and ms = 611 MeV, namely the density dependent mu together with ∗ 0.25 0.61 0.69 0.44 ms = ms. All the charge radii are in units of fm. The empirical result in vacuum is from Ref. [77]. 0.50 0.62 0.69 0.44 expt 0.75 0.66 0.74 0.44 ρB/ρ0 rK ru rs r 0.00 0.59 0.65 0.44 0.56 ± 0.03 1.00 0.69 0.79 0.44 0.25 0.62 0.69 0.44 1.25 0.74 0.86 0.44 0.50 0.65 0.73 0.44 0.75 0.68 0.77 0.44 1.00 0.71 0.81 0.44 in medium is identical to that in vacuum, because the s-quark properties are not modiﬁed in medium in the present approach. 1.25 0.74 0.85 0.44 + Table VI shows that at normal nuclear density ρ0, the K and 10

1 1 (a) (b) [ρB = 0.00ρ0] 0.8 [ρB = 0.50ρ0] 0.8 [ρB = 1.00ρ0] ) ) 2 0.6 2 0.6 Q Q ( ( u K 0.4 K 0.4 F F

0.2 0.2

0 0 0 1 2 3 4 0 1 2 3 4 Q2 (GeV2) Q2 (GeV2)

1 1 (c) (d)

0.8 0.8 ) ) 2 2

0.6 Q 0.6 Q ( ( K ¯ s K F F

0.4 2 0.4 − Q 0.2 0.2

0 0 0 1 2 3 4 0 1 2 3 4 Q2 (GeV2) Q2 (GeV2)

/ FIG. 8. Same as Fig. 7 but for the NJL-based approach. The lines are for ρB ρ0 = 0.00 (dashed line), 0.50 (dash-dotted line), and 1.00 (dotted line), respectively u-quark charge radii increase respectively about 17% and 22% tions. Namely, the NJL-based approach that the in-medium ∗ ∗ relative to those corresponding in vacuum. dynamical light-quark mass Ml = Mu calculated by the NJL Overall we find the results calculated by both the QMC- model using the in-medium input for the current-quark mass ∗ based and NJL-based approaches show very similar behavior ml (in vacuum ml = 16.4 MeV) generated by the quark-meson for the in-medium K+ and quark EMFFs up to the nuclear coupling model, yields nearly the same results as those cal- density ρB ≤ 1.25ρ0, where the in-medium constituent quark culated in the QMC-based approach, that uses the in-medium ∗ masses for the both approaches are almost the same in this constituent quark mass ml (in vacuum ml = 400 MeV) cal- region as in Fig. 1. Consequently, they give nearly the same culated in the QMC model. In both approaches, the strange results for the in-medium K+ EMFF and quark charge radii quark constituent quark masses are kept the same as that in ∗ ∗ in this density region, as can be seen in Tables. V and VI. vacuum, Ms = Ms = ms = ms = 611 MeV. The feature of The difference in the in-medium K+ charge radii in the both yielding nearly the similar predictions in the two approaches, approaches is about 2% at normal nuclear density. holds up to around the normal nuclear densities. The difference in the two approaches appears only in the higher density region (as shown in Fig. 1). At normal nuclear ∗ density, the values of the constituent light-quark mass ml cal- V. SUMMARY AND CONCLUSION culated in the QMC-based approach is smaller about 10% than that of calculated in the NJL-based approach. We have studied the kaon properties and the space-like K+- The NJL model has been used as a complementary model meson electromagnetic form factor (EMFF) in symmetric nu- to the MIT bag model in the variants of the QMC model clear matter in the Nambu-Jona-Lasinio (NJL) model. We in Ref. [73]. With similar motivation, we have studied the in- take into account the in-medium effects in the dressed quark- medium kaon properties in the NJL model using the in-medium photon vertex. The in-medium constituent light-quark masses inputs calculated by the QMC model. These are either using are calculated by the two different approaches to use in the the in-medium current quark properties for the NJL-based NJL model. The results of the two different approaches are approach, or the in-medium constituent quark properties for compared, and they turned out to give very similar predic- the QMC-based approach as mentioned above. 11

(a) (b) 1 ρB = 0.00ρ0 1 ρB = 0.50ρ0

0.8 0.8 ) ) 2 2 Q Q

( 0.6 ( 0.6 K K

F 0.4 F 0.4 NJL NJL 0.2 QMC 0.2 QMC

0 0 0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 Q2 (GeV2) Q2 (GeV2)

0.8 ) 2 Q

( 0.6 K

F 0.4 NJL 0.2 QMC

0 0 0.1 0.2 0.3 0.4 Q2 (GeV2)

+ / FIG. 9. In-medium K EMFF in the QMC-based and the NJL-based approaches, compared with the existing data in vacuum (ρB ρ0 = 0.00) at low Q2. Experimental data are in vacuum are from Ref. [30].

By the NJL-based approach we have calculated the effective In contrast, the in-medium u-quark EMFF as well as that of ∗ + kaon mass mK , kaon decay constant, and kaon-quark coupling K , decrease appreciably as nuclear density increases. This constant in symmetric nuclear matter. We predict that, at leads to a larger u-quark charge radius as well as that of K+ normal nuclear density, the constituent quark mass, effective in symmetric nuclear matter. At normal nuclear density, the + kaon mass, kaon-quark coupling constant, and kaon leptonic K charge radius rK increases about 20% relative to that in decay constant decrease respectively about 33%, 22%, 4%, vacuum in the QMC-based approach, while it yields about 17 and 3% relative to those corresponding in vacuum. % increase in the NJL-based approach. These results indicate Alternatively, in the QMC-based approach using the ef- that the valence s-quark gives a significant contribution for the + fective constituent quark mass of the light quark m∗ directly K EMFF in vacuum as well as in symmetric nuclear matter, l 2 calculated in the QMC model, we have calculated the same for larger Q and larger nuclear density ρB. Our predictions quantities using the NJL model as those calculated in the NJL- can be verified by experiments such as the Compressed Bary- based approach. In this case we predict, at normal nuclear onic Matter (CBM) experiment at GSI [78]. density, the constituent light-quark mass, effective kaon mass, kaon-quark coupling constant, and kaon lepton decay constant decrease respectively about 39%, 25%, and 4%, and 4% com- ACKNOWLEDGMENTS pared to those corresponding in vacuum. Based on the in-medium kaon properties calculated in the P.T.P.H. and K.T. thank Yongseok Oh for comments on the both approaches, the QMC-based and the NJL-based, we have manuscript and useful conversations. K.T. thanks the Asia Pa- studied the K+ EMFF in symmetric nuclear matter using the cific Center for Theoretical Physics (APCTP) and Kyungpook NJL model. We have found that the in-medium s-quark EMFF National University for warm hospitality and excellent sup- increases as nuclear density increases, which is unexpected. ports during his visit. P.T.P.H. was supported by the Ministry However, the increasing rate against the increase of nuclear of Science, ICT and Future Planning, Gyeongsangbuk-do, and density is very small. This leads to nearly unmodified in- Pohang City. The work of K.T. was supported by the Con- medium s-quark charge radius relative to that in vacuum, irre- selho Nacional de Desenvolvimento Científico e Tecnológico spective of nuclear densities studied. (CNPq) Process, No. 313063/2018-4, and No. 426150/2018-0, 12 and Fundação de Amparo à Pesquisa do Estado de São Paulo clear Physics and Applications (INCT-FNA), Brazil, Process. (FAPESP) Process, No. 2019/00763-0, and was also part of No. 464898/2014-5, and FAPESP Temático, Brazil, Process, the projects, Instituto Nacional de Ciência e Tecnologia – Nu- No. 2017/05660-0.

[1] J. Aichelin and C. M. Ko, Subthreshold Kaon Production as 371, 163 (1996). a Probe of the Nuclear Equation of State, Phys. Rev. Lett. 55, [21] J. Bijnens and P. Talavera, Pion and kaon electromagnetic form 2661 (1985). factors, JHEP 03, 046 (2002). [2] K. Tsushima, S. W. Huang and A. Faessler, The Role of the Delta [22] S.-I. Nam and H.-C. Kim, Electromagnetic form factors of the (1920) resonance for kaon production in heavy ion collisions, pion and kaon from the instanton vacuum, Phys. Rev. D 77, Phys. Lett. B 337, 245 (1994). 094014 (2008). [3] K. Tsushima, S. W. Huang and A. Faessler, Resonance model [23] U. Raha and A. Aste, Electromagnetic pion and kaon form of pi Delta —> Y K for kaon production in heavy ion collisions, factors in light-cone resummed perturbative QCD, Phys. Rev. J. Phys. G 21, 33 (1995). D 79, 034015 (2009). [4] T. Waas, N. Kaiser, and W. Weise, Effective kaon masses in [24] E. O. da Silva, J. P. B. C. de Melo, B. El-Bennich, and V. S. dense nuclear and neutron matter, Phys. Lett. B 379, 34 (1996). Filho, Pion and kaon elastic form factors in a refined light-front [5] R. Barth et al. [KaoS Collaboration], Subthreshold production model, Phys. Rev. C 86, 038202 (2012). of kaons and anti-kaons in nucleus-nucleus collisions at equiv- [25] Y. Ninomiya, W. Bentz, and I. C. Cloët, Dressed quark mass alent beam energies, Phys. Rev. Lett. 78, 4007 (1997). dependence of pion and kaon form factors, Phys. Rev. C 91, [6] K. Tsushima, A. Sibirtsev and A. W. Thomas, Resonance model 025202 (2015). study of strangeness production in p p collisions, Phys. Lett. B [26] A. F. Krutov, S. V. Troitsky, and V. E. Troitsky, The K-meson 390, 29 (1997). form factor and charge radius: Linking low-energy data to future [7] A. Sibirtsev, K. Tsushima and A. W. Thomas, A Clue to the Jefferson Laboratory measurements, Eur. Phys. J. C 77, 464 mechanism of Lambda K+ production in p p reactions, Phys. (2017). Lett. B 421, 59 (1998). [27] Z. Abidin and P. T. P. Hutauruk, Kaon form factor in holographic [8] K. Tsushima, K. Saito, A. W. Thomas and S. V. Wright, In- QCD, Phys. Rev. D 100, no. 5, 054026 (2019). medium kaon and antikaon properties in the quark meson cou- [28] P. T. P. Hutauruk, I. C. Cloët, and A. W. Thomas, Flavor depen- pling model, Phys. Lett. B 429, 239 (1998), Erratum: [Phys. dence of the pion and kaon form factors and parton distribution Lett. B 436, 453 (1998)]. functions, Phys. Rev. C 94, 035201 (2016). [9] R. Shyam, p p → p K+ Lambda reaction in an effective La- [29] P. T. P. Hutauruk, Kaon structure in the confining Nambu– grangian model, Phys. Rev. C 60, 055213 (1999). Jano-Lasinio model, talk at the 45th International Meeting on [10] K. Tsushima, A. Sibirtsev, A. W. Thomas and G. Q. Li, Res- Fundamental Physics, arXiv:1801.06299. onance model study of kaon production in baryon baryon re- [30] S. R. Amendolia et al., A measurement of the kaon charge actions for heavy ion collisions, Phys. Rev. C 59, 369 (1999), radius, Phys. Lett. B 178, 435 (1986). Erratum: [Phys. Rev. C 61, 029903 (2000)]. [31] M. Carmignotto et al., JLAB FPI-2 and E93-018 Collaboration, [11] F. Laue et al. [KaoS Collaboration], Medium effects in kaon Separated kaon electroproduction cross section and the kaon and anti-kaon production in nuclear collisions at subthreshold form factor from 6 GeV JLab data, Phys. Rev. C 97, 025204 beam energies, Phys. Rev. Lett. 82, 1640 (1999). (2018). [12] A. Sibirtsev, K. Tsushima, W. Cassing and A. W. Thomas, The [32] G. H. S. Yabusaki, J. P. B. C. de Melo, W. de Paula, K. Tsushima, Role of the P(11)(1710) in the N N —> N Sigma K reaction, and T. Frederico, In-medium K+ electromagnetic form factor Nucl. Phys. A 646, 427 (1999). with a symmetric vertex in a light front approach, Few-Body [13] K. Tsushima, A. Sibirtsev and A. W. Thomas, Strangeness Syst. 59, 37 (2018). production from pi N collisions in nuclear matter, Phys. Rev. C [33] P. T. P. Hutauruk, Y. Oh, and K. Tsushima, Electroweak prop- 62, 064904 (2000). erties of pions in a nuclear medium, Phys. Rev. C 99, 015202 [14] K. Tsushima, A. Sibirtsev and A. W. Thomas, Kaon properties (2019). and cross-sections in nuclear medium, J. Phys. G 27, 349 (2001). [34] Y.Nambu and G. Jona-Lasinio, Dynamical model of elementary [15] F. Gao, L. Chang, Y.-X. Liu, C. D. Roberts, and P. C. Tandy, Ex- particles based on an analogy with superconductivity. I, Phys. posing strangeness: Projections for kaon electromagnetic form Rev. 122, 345 (1961). factors, Phys. Rev. D 96, 034024 (2017). [35] Y.Nambu and G. Jona-Lasinio, Dynamical model of elementary [16] P. T. P. Hutauruk, W. Bentz, I. C. Cloët, and A. W. Thomas, particles based on an analogy with superconductivity. II, Phys. Charge symmetry breaking effects in pion and kaon structure, Rev. 124, 246 (1961). Phys. Rev. C 97, 055210 (2018). [36] J. J. Aubert et al., The ratio of the nucleon structure functions N [17] N. Zovko, Pion and kaon form factors and heavy vector mesons, F2 for iron and deuterium, Phys. Lett. B 123B, 275 (1983). Phys. Lett. 51B, 54 (1974). [37] D. F. Geesaman, K. Saito and A. W. Thomas, The nuclear EMC [18] V. Bernard and U.-G. Meissner, Electromagnetic structure of effect, Ann. Rev. Nucl. Part. Sci. 45, 337 (1995). the pion and the kaon, Phys. Rev. Lett. 61, 2296 (1988), 61, [38] S. Strauch et al., Polarization transfer in the 4He(eì, e0pì)3H 2973(E) (1988). reaction up to Q2 = 2.6 (GeV/c)2, Phys. Rev. Lett. 91, 052301 [19] W. W. Buck, R. A. Williams, and H. Ito, Elastic charge form (2003). factors for K mesons, Phys. Lett. B 351, 24 (1995). [39] G. E. Brown and M. Rho, Scaling effective Lagrangians in a [20] C. J. Burden, C. D. Roberts, and M. J. Thomson, Electromag- dense medium, Phys. Rev. Lett. 66, 2720 (1991). netic form factors of charged and neutral kaons, Phys. Lett. B [40] R. S. Hayano and T. Hatsuda, Hadron properties in the nuclear 13

medium, Rev. Mod. Phys. 82, 2949 (2010). medium, arXiv:1908.02406 [hep-ph]. [41] S. Leupold, V. Metag, and U. Mosel, Hadrons in strongly [60] S. P. Klevansky, The Nambu–Jona-Lasinio model of Quantum interacting matter, Int. J. Mod. Phys. E 19, 147 (2010). Chromodynamics, Rev. Mod. Phys. 64, 649 (1992). [42] K. Tsushima, K. Saito, A. W. Thomas, and S. V. Wright, In- [61] P.A. M. Guichon, A possible quark mechanism for the saturation medium kaon and antikaon properties in the quark-meson cou- of nuclear matter, Phys. Lett. B 200, 235 (1988). pling model, Phys. Lett. B 429, 239 (1998), 436, 453(E) (1998). [62] P. A. M. Guichon, J. R. Stone, and A. W. Thomas, Quark- [43] K. Saito, K. Tsushima, and A. W. Thomas, Nucleon and hadron meson-coupling (QMC) model for finite nuclei, nuclear matter structure changes in the nuclear medium and the impact on and beyond, Prog. Part. Nucl. Phys. 100, 262 (2018). observables, Prog. Part. Nucl. Phys. 58, 1 (2007). [63] P. T. P. Hutauruk, Y. Oh, and K. Tsushima, Impact of medium [44] G. Krein, A. W. Thomas, and K. Tsushima, Nuclear-bound modifications of the nucleon weak and electromagnetic form quarkonia and heavy-flavor hadrons, Prog. Part. Nucl. Phys. factors on the neutrino mean free path in dense matter, Phys. 100, 161 (2018). Rev. D 98, 013009 (2018). [45] A. Hosaka, T. Hyodo, K. Sudoh, Y. Yamaguchi, and S. Yasui, [64] P. T. P. Hutauruk, Y. Oh and K. Tsushima, Effects of medium Heavy hadrons in nuclear matter, Prog. Part. Nucl. Phys. 96, 88 modifications of nucleon form factors on neutrino scattering in (2017). dense matter, arXiv:1901.05986 [nucl-th]. [46] V. Metag, M. Nanova, and K.-T. Brinkmann, In-medium prop- [65] K. Tsushima, In-medium properties of the low-lying strange, erties of mesons, EPJ Web Conf. 134, 03003 (2017). charm, and bottom baryons in the quark-meson coupling model, [47] T. Horikawa and W. Bentz, Medium modifications of nucleon Phys. Rev. C 99, 014026 (2019). electromagnetic form factors, Nucl. Phys. A 762, 102 (2005). [66] G. Krein, A. W. Thomas, and K. Tsushima, Fock terms in the [48] D. H. Lu, A. W. Thomas, K. Tsushima, A. G. Williams, and quark-meson coupling model, Nucl. Phys. A 650, 313 (1999). K. Saito, In-medium electron-nucleon scattering, Phys. Lett. B [67] P. A. M. Guichon, K. Saito, E. Rodionov, and A. W. Thomas, 417, 217 (1998). The role of nucleon structure in finite nuclei, Nucl. Phys. A 601, [49] E. Ruiz Arriola, Chr. V. Christov, and K. Goeke, Medium effects 349 (1996). on nucleon properties, Phys. Lett. B 225, 22 (1989). [68] P. A. M. Guichon, A. W. Thomas, and K. Tsushima, Binding [50] U. Vogl and W. Weise, The Nambu and Jona-Lasinio model: Its of hypernuclei in the latest quark-meson coupling model, Nucl. applications for hadrons and nuclei, Prog. Part. Nucl. Phys. 27, Phys. A 814, 66 (2008). 195 (1991). [69] A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, and V. F. [51] G. E. Brown and M. Rho, Chiral restoration in hot and/or dense Weisskopf, New extended model of hadrons, Phys. Rev. D 9, matter, Phys. Rep. 269, 333 (1996). 3471 (1974). [52] J. W. Holt and K. E. Rockcliffe, Partial restoration of chiral [70] B. D. Serot and J. D. Walecka, The relativistic nuclear many symmetry in hot and dense neutron matter, Int. J. Mod. Phys. E body problem, Adv. Nucl. Phys. 16, 1 (1986). 26, 1740010 (2017). [71] B. D. Serot and J. D. Walecka, Recent progress in Quantum [53] D. Jido, T. Hatsuda, and T. Kunihiro, In-medium ππ correlation Hadrodynamics, Int. J. Mod. Phys. E 6, 515 (1997). induced by partial restoration of chiral symmetry, Phys. Rev. D [72] S. Maedan, Influence of current mass on the spatially inhomo- 63, 011901(R) (2000). geneous chiral condensate, Prog. Theor. Phys. 123, 285 (2010). [54] D. Jido, T. Hatsuda, and T. Kunihiro, In-medium pions and [73] D. L. Whittenbury, M. E. Carrillo-Serrano, and A. W. Thomas, partial restoration of chiral symmetry: A model-independent Quark-meson coupling model based upon the Nambu–Jano- analysis, Prog. Theor. Phys. Suppl. 168, 478 (2007). Lasinio model, Phys. Lett. B 762, 467 (2016). [55] V. Bernard and U.-G. Meissner, Properties of vector and axial- [74] I. C. Cloët, G. A. Miller, E. Piasetzky, and G. Ron, Neutron vector mesons from a generalized Nambu–Jona-Lasinio model, properties in the medium, Phys. Rev. Lett. 103, 082301 (2009). Nucl. Phys. A 489, 647 (1988). [75] V. Bernard and U.-G. Meissner, Meson properties at finite [56] S. Klimt, M. Lutz, U. Vogl, and W. Weise, Generalized SU(3) density from SU(3)f quark dynamics, Phys. Rev. D 38, 1551 Nambu–Jona-Lasinio model. (I) Mesonic modes, Nucl. Phys. A (1988). 516, 429 (1990). [76] G. P. Lepage and S. J. Brodsky, Exclusive processes in Quantum [57] V. Bernard, R. L. Jaffe, and U.-G. Meissner, Strangeness mixing Chromodynamics: Evolution equations for hadronic wavefunc- and quenching in the Nambu–Jona-Lasinio model, Nucl. Phys. tions and the form factors of mesons, Phys. Lett. 87B, 359 B 308, 753 (1988). (1979). [58] I. C. Cloët, W. Bentz, and A. W. Thomas, Role of diquark [77] M. Tanabashi et al., Particle Data Group, Review of particle correlations and the pion cloud in nucleon elastic form factors, physics, Phys. Rev. D 98, 030001 (2018). Phys. Rev. C 90, 045202 (2014). [78] W. Cassing and E. L. Bratkovskaya, Hadronic and electromag- [59] P. T. P. Hutauruk, J. J. Cobos-Martinez, Y. Oh and K. Tsushima, netic probes of hot and dense nuclear matter, Phys. Rept. 308 Valence-quark distributions of pions and kaons in a nuclear (1999) 65.