{"id":"d4525e23-0850-4a44-b6bc-951b82a37d78","arxiv_id":"2412.09883","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A hybrid light-front quark model plus quark-meson coupling calculation predicts that in nuclear matter the charge radii of pseudoscalar mesons grow, driven mainly by the light quark sector.","lead":"Calculates how the electromagnetic form factors and charge radii of pions, kaons, D mesons, and B mesons change when the mesons are embedded in dense nuclear matter. Finds that absolute charge radii grow with density, with the medium effects concentrated in the light-quark sector.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The in-medium charge-radius increase rests on holding the LFQM width parameter beta fixed; a density-dependent beta could change the predicted rates and even the flavor ordering, so the central quantitative claim is not yet secured.","rationale":"The reader's weakest-assumption analysis already identifies the constant-beta approximation as the load-bearing point, and my independent reading agrees. The central claim is not a parameter-free consequence of the model; it relies on inserting QMC-modified quark masses into a free-space LFWF whose width is held fixed. Because beta controls the transverse size of the Gaussian wave function, it directly controls the form factor slope and hence the charge radius. The paper's own caveat in Sec. III B and the restriction to rho <= 1.5 rho0 due to the pion decay constant becoming negative both signal that the in-medium wave function is being used outside its calibrated domain. I do not see a separate internal inconsistency in the derivation of Eqs. (56) and (59): the vector-potential shifts are handled consistently for the averaged radii, and the normalization of the quark-sector form factors at Q^2 = 0 is preserved. The comparison with prior NJL and Bethe-Salpeter-based calculations is a point in the paper's favor, since the qualitative trend is shared, but those calculations have their own model assumptions. The concrete test proposed here is a minimal, feasible check: re-fit or scale beta within the same framework and see whether the predicted radius growth and flavor ordering are stable. If they are stable, the constant-beta assumption is a benign simplification; if not, the quantitative predictions in Table IV should not be presented as robust. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":21944,"tokens_out":9778,"duration_ms":124666,"concrete_test":"Recompute Table IV with a density-dependent beta_{q\\bar q} while keeping all other inputs fixed. Use several controlled variants: (i) scale beta^* = beta (m_q^*/m_q)^n with n = +1/2 and n = -1/2; (ii) scale beta^* by the QMC bag-radius ratio beta R_N/R_N^*(rho); and (iii) re-fit beta^* by requiring the LFQM to reproduce the in-medium decay constants computed in the same framework (the authors' Ref. [58]) for rho <= 1.5 rho0. Compare the resulting |<r^2>| at rho/rho0 = 1.0 and 1.5 with Table IV. If the monotonic increase of all meson radii and the stated flavor ordering survive every variant, the beta-constant assumption is not load-bearing; if any meson's sign or the ordering changes, the central claim must be reframed as conditional on beta remaining fixed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The observable content of the paper is the slope of the in-medium form factor at Q^2 = 0 through Eq. (48). In the LFQM this slope is controlled by the Gaussian width beta_{q\\bar q} in Eq. (7), which sets the transverse-momentum scale of the wave function, and by the combination A* = x m_Q^* + (1-x) m_q^* entering Eqs. (56) and (59). The authors modify only A* through the QMC scalar potential, while in Sec. III B they state: \"for simplicity, we assume it [beta] to be constant in the present work.\" This is a mixed approximation: in free space beta is not an independent physical constant but a parameter fitted to meson masses, decay constants, and charge radii. The QMC model itself yields density-dependent bag radii and effective quark masses, so the transverse-size scale of the meson wave function is not obviously frozen. Since the form factor slope near Q^2 = 0 scales roughly as 1/beta^2, a modest density dependence of beta changes the magnitudes in Table IV by an amount comparable to the predicted 50-100% increases, and flavor-dependent changes in beta can alter the stated ordering (pion fastest, charged D slowest, kaon and B nearly equal). In some variants the sign of a neutral-meson radius change could also be affected. The central qualitative claim may survive, but without quantifying the beta dependence the quantitative predictions in Table IV and Fig. 9 are not robust.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper combines the light-front quark model (LFQM) with the quark-meson coupling (QMC) model to compute the spacelike electromagnetic form factors (EMFFs) and charge radii of the pion, kaon, D mesons, and B mesons in symmetric nuclear matter. The LFQM parameters are fixed to free-space meson observables and the QMC parameters are fixed to nuclear saturation properties. The authors find that with increasing nuclear density the charged-meson EMFFs fall faster with Q^2, the neutral-meson EMFFs rise, and the absolute values of the charge radii increase, with the rate depending on the quark flavor content. They also decompose the EMFFs into light- and heavy-quark sector contributions and conclude that the medium affects mainly the light-quark sector. The calculation is restricted to densities up to 1.5 rho0 because the pion decay constant becomes negative at higher densities, and the Gaussian width beta is assumed to be density independent.","tokens_in":22292,"tokens_out":14565,"duration_ms":144485,"significance":"If the quantitative predictions are robust, this is a useful model study of in-medium hadron structure in a regime where lattice QCD is not directly applicable and where JLab 12 GeV and EIC programs may eventually provide constraints. The paper has clear strengths: the free-space EMFFs are benchmarked against experimental and lattice data, the quark-flavor decomposition is systematic, and the main assumptions and limitations are stated explicitly. However, the central quantitative content is the density dependence of the charge-radii slope near Q^2=0, and this content is controlled by an assumption that is not yet tested in the manuscript. The conclusions should therefore be regarded as conditional until the sensitivity to that assumption is quantified.","major_comments":[{"comment":"The Gaussians width beta_{q qbar} in Eq. (7) controls the transverse size of the meson wave function and therefore the slope of the EMFF at Q^2=0, i.e., the charge radius in Eq. (48). In the in-medium calculation the authors modify the effective light-quark mass m*_q and the vector-potential shifts, but keep beta fixed, stating in Sec. III B that \"for simplicity, we assume it to be constant in the present work.\" No sensitivity estimate is given. Since |<r^2>| scales approximately as 1/beta^2, a density-dependent change of beta of only 10% changes the predicted radius by about 20%, which is the same order as the reported 50-100% increases at rho=rho0 in Table IV. A flavor-dependent beta* could also alter the claimed ordering, for example the near equality of the K and B meson radii. I ask the authors to add a sensitivity analysis, for example by varying beta by plausible amounts or by estimating beta* from the QMC bag radius or the in-medium meson mass, and to state explicitly which conclusions survive.","section":"Sec. III B, Eqs. (56), (59), Table IV, Fig. 9"},{"comment":"The definitions of the plus-momentum variables are ambiguous. Equation (51) defines p*0_i as the in-medium energy including the vector potential, while Eq. (54) adds V_qomega to p*+_q, and Eq. (57) does the same in the numerator of the shifted variable. The text after Eq. (59) states that in the meson rest frame P*+ = M*, which appears inconsistent with Eq. (52), where P*0 = E*_M +/- V_qomega for q Qbar and Q qbar systems. If p*+_q and P*+ already contain the vector potential, the shifts in Eqs. (54) and (57) double-count it; if they do not, this should be stated explicitly. Because the individual K+ and K- radii and the vector-potential-induced splitting shown in Fig. 9 depend on these shifts, the derivation should be clarified before the individual-meson results can be assessed.","section":"Sec. III B, Eqs. (54), (57), (59), and text after Eq. (59)"}],"minor_comments":[{"comment":"The text says the pion charge radius increases to approximately 1.45 times its free-space value, with <r*2_pi> = 0.897 fm^2 at rho=rho0. This factor refers to the radius sqrt(<r^2>), not to <r^2> itself, whose ratio is about 2.1; please phrase this as sqrt(<r*2>) to avoid confusion.","section":"Sec. IV B"},{"comment":"The text refers to the quark-sector curves as \"blue line\" and \"red line,\" while the figure legends and captions use orange and green; please make the color references consistent.","section":"Sec. IV C and Figs. 5, 7"},{"comment":"The in-medium meson masses M* from Ref. [58] enter P*+ and the light-front variable in Eqs. (52) and (57), but the present paper does not tabulate their density dependence. A short table or a reminder of the relevant values would improve reproducibility.","section":"Sec. III B"},{"comment":"The authors note that the calculation is limited to rho <= 1.5 rho0 because the pion decay constant may become negative at higher densities. It would be helpful to state how close the 1.5 rho0 results are to that breakdown, so that the reader knows whether the highest-density predictions are already in a regime where the model is strained.","section":"Sec. IV B"},{"comment":"Reference [25] appears incomplete: it lists the journal and article number but no year or volume; please complete the entry.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the unquantified density independence of beta; this is precisely the input that controls the observable content of the paper. The vector-potential plus-momentum notation in Eqs. (54)-(59) should also be cleaned up, because it currently makes the individual K+ vs K- results difficult to verify. If the authors provide a sensitivity analysis for beta and show that the qualitative trends and the flavor ordering survive, the paper could become acceptable. I do not see a reason to reject the manuscript outright; the framework is transparent and the free-space benchmarks are meaningful."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Put it on the pile. The genuinely new content is the in-medium D and B meson electromagnetic form factors and the flavor decomposition showing medium effects sit almost entirely in the light quark sector. The free-space form factors are anchored to data and lattice results, the QMC sector reproduces saturation, and the authors state their main caveats up front: beta held fixed, densities capped at 1.5 rho0, no error bars. That is honest work, not a dressed-up fit.\n\nThe paper extends the earlier Arifi-Hutauruk-Tsushima framework to heavy-light mesons, so the novelty is real but incremental for the pion and kaon, which already show similar trends in NJL and BSA-LFQM papers. The heavy-light results are what is new.\n\nThe soft spots are what the reader flagged. The constant-beta assumption is explicitly stated, and it is load-bearing because the charge radius slope scales roughly as 1/beta^2. The authors justify it by saying beta is a short-distance scale, but the QMC model itself allows bag radii and effective masses to move with density, so freezing the transverse width is a choice, not a derivation. A density-dependent beta could shift the magnitudes in Table IV by the same order as the predicted 50-100% increases, and flavor-dependent changes could reorder the pion/kaon/B pattern. That does not kill the qualitative claim that charge radii grow with density—the effective light-quark mass dropping pushes that direction in every variant—but it does mean the quantitative table is not secured. The pion decay constant turning negative above 1.5 rho0 is an honest limitation, and the authors point to the BSA treatment for higher densities. The difference from the BSA pion EMFF is noted but not resolved; the reader should know the in-medium pion result is treatment-sensitive.\n\nOne more thing: the paper averages over K+ and K- and takes absolute values to cancel the vector potential. That is a modelling choice that hides the vector-potential dependence except in the radius differences in Fig. 9. It is presented clearly, so no foul, but it means the 'unified' in-medium radii are scalar-only quantities.\n\nBottom line: this deserves a serious referee. The calculation is transparent, the limitations are stated, the heavy-light predictions are new, and the qualitative medium trend is robust even if the quantitative slope is not. A referee should push for a sensitivity study on beta and for error estimates, but this is not a desk reject. I would bring it to reading group and probably cite the heavy-light tables.","headline":"New heavy-light in-medium form factor predictions with a honestly stated but load-bearing constant-beta assumption; worth a serious referee.","tokens_in":22861,"tokens_out":1757,"would_cite":true,"duration_ms":19402,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Meson charge radii grow with nuclear density, fastest for pions.","keywords":["in-medium form factors","light-front quark model","quark-meson coupling model","charge radii","pseudoscalar mesons","symmetric nuclear matter","heavy-light mesons","pion and kaon"],"falsifier":"Measure the pion charge radius in nuclear matter at saturation density, for example through the energy shifts of deeply bound pionic atoms or electron scattering on nuclei: the central claim predicts $r_\\pi^2 \\simeq 0.897$ fm$^2$ at $\\rho_0$, and a measured value close to the free-space $0.427(10)$ fm$^2$ would falsify it.","tokens_in":21707,"feed_emoji":"⚛️","tokens_out":8348,"duration_ms":85761,"temperature":0.7,"pith_summary":"This paper predicts how the electromagnetic form factors of pseudoscalar mesons change when the meson sits in symmetric nuclear matter, using a light-front quark model fed by quark-meson coupling model mean fields. It claims that with rising density the form factors of charged mesons fall off faster with momentum transfer, neutral meson form factors rise, and the absolute charge radius of every meson grows. The growth rate is flavor dependent: pions change most, charged D mesons least, kaons and B mesons at nearly the same rate. The authors argue the medium acts mainly on light quarks, leaving heavy-quark contributions almost untouched, which gives a clean quark-level signature of nuclear medium effects.","feed_headline":"Nuclear density inflates meson charge radii, fastest for pions","feed_subtitle":"A quark-model calculation predicts charged meson form factors steepen and neutral ones rise, changing measurable radii.","key_machinery":"The engine is a light-front quark model (a constituent-quark description of mesons as quark-antiquark bound states on a fixed light-front plane) with a Gaussian radial wave function of width $\\beta_{q\\bar q}$, combined with the quark-meson coupling model (a relativistic mean-field description in which nuclear scalar and vector fields act on light quarks). The scalar mean field reduces the light-quark mass to $m_q^*$, and the vector mean field shifts quark energies and the longitudinal momentum variable $x$ to $\\tilde{x}^*$; the vector shift cancels for equal-mass $q\\bar q$ pairs and does not affect heavy quarks. The form factor is a convolution of initial and final wave functions with the shifted variables, and its slope at $Q^2=0$ yields the charge radius. The paper keeps $\\beta_{q\\bar q}$ fixed at its free-space value, so the entire density dependence enters through the quark mass and energy shifts.","core_discovery":"The paper's central claim is that the in-medium modifications of the spacelike electromagnetic form factors are controlled by the light-quark sector. In symmetric nuclear matter the light quark acquires an effective mass $m_q^*$ smaller than its free value, and this alone makes the form factor of a charged meson drop more steeply with $Q^2$ as the density rises, while a neutral meson's form factor rises with density. The corresponding absolute charge radii therefore increase with density: at saturation density the pion radius squared reaches $0.897$ fm$^2$ against $0.427(10)$ fm$^2$ in free space, and the kaon reaches $0.697$ fm$^2$. The flavor decomposition shows the heavy-quark ($s$, $c$, $b$) contribution is nearly density independent, so the ordering of radii and the charged-versus-neutral asymmetry are carried by the light quark.","pith_inferences":["Our inference: the constant-$\\beta_{q\\bar q}$ assumption is the main lever; if the medium changes the transverse size of the wave function, the direction of the radius change could survive but the flavor ordering and the kaon-B coincidence could shift.","Our inference: applying the same machinery to distribution amplitudes or generalized parton distributions would predict that medium effects concentrate in the light-quark partonic content, linking these radius changes to nuclear partonic modifications.","Our inference: because the model's pion decay constant turns negative above $1.5\\rho_0$, the density trend beyond that point is not trustworthy; testing the trend would require a model that stabilizes the decay constant.","Our inference: a finite-density lattice calculation of the neutral kaon form factor, once the sign problem is bypassed, would be a sharp check, since the predicted increase with density is opposite to naive screening intuition."],"forward_implications":["At saturation density the pion's charge radius squared is about 1.45 times its free-space value, so measurements of pionic atoms or electron scattering on nuclei could see a clearly inflated pion in matter.","The charged-versus-neutral asymmetry (faster fall-off for charged, rising form factor for neutral kaons and D mesons) offers a distinctive experimental signature that does not depend on the overall normalization.","Because heavy-quark form factors are nearly unchanged, observables built from heavy-light mesons in nuclei isolate the light-quark medium response.","Medium effects fade at higher $Q^2$, so upcoming electron scattering measurements should target the low-$Q^2$ region to test these predictions.","The nearly identical growth of kaon and B-meson radii, despite very different quark masses, is a parameter-free feature of the model that a measurement could confirm or rule out."],"supporting_citations":[{"why":"It supplies the free-space light-front quark model wave functions and parameter set (quark masses and $\\beta$) used throughout.","marker":"[26]"},{"why":"It establishes the previous in-medium decay constants and distribution amplitudes, and the shifted longitudinal momentum variable used here.","marker":"[58]"},{"why":"It provides the quark-meson coupling model mean-field framework that yields the in-medium light-quark mass and energy shifts.","marker":"[69]"},{"why":"It gives the in-medium pion form factors in a related light-front approach, used as a comparison for the density trend.","marker":"[51]"},{"why":"It gives in-medium kaon form factors and the behavior at densities above $1.5\\rho_0$, used as a comparison.","marker":"[52]"},{"why":"It provides a recent field-theoretic model calculation of in-medium meson structures, used to cross-check the density dependence.","marker":"[21]"},{"why":"It supplies free-space lattice QCD form factors at high $Q^2$, used to benchmark the model in free space.","marker":"[46]"},{"why":"It supplies the lattice charge radius of the kaon quoted in the comparison table.","marker":"[44]"},{"why":"It supplies lattice QCD data for D-meson couplings used to benchmark the heavy-light predictions.","marker":"[41]"}],"fun_headline_variants":["Meson charge radii swell with nuclear density","Pion radius grows fastest in dense nuclear matter","Light quarks inflate meson radii in medium","Nuclear medium modifies meson form factors via light quarks","Dense matter puffs up mesons, pions most"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the Gaussian width $\\beta_{q\\bar q}$ of the meson wave function is unchanged in the medium; if the medium alters that transverse size, the predicted growth of the charge radii and its flavor ordering could change.","fun_headline_variants_meta":{"raw":{"variants":["Meson charge radii swell with nuclear density","Pion radius grows fastest in dense nuclear matter","Light quarks inflate meson radii in medium","Nuclear medium modifies meson form factors via light quarks","Dense matter puffs up mesons, pions most"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000474,"raw_usage":{"total_tokens":2357,"prompt_tokens":953,"completion_tokens":1404,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":1328}},"tokens_in":569,"tokens_out":1404,"duration_ms":9998,"temperature":1.0,"reasoning_tokens":1328,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:37:29.951837+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the pion charge radius in nuclear matter at saturation density, for example through the energy shifts of deeply bound pionic atoms or electron scattering on nuclei: the central claim predicts $r_\\pi^2 \\simeq 0.897$ fm$^2$ at $\\rho_0$, and a measured value close to the free-space $0.427(10)$ fm$^2$ would falsify it.","supporting_citations":[{"cited_title":"Suzuki et al., Precision spectroscopy of pionic 1s states of Sn nuclei and evidence for partial restoration of chiral symmetry in the nuclear medium, Phys","cited_arxiv_id":null,"evidence_quote":"It supplies the free-space light-front quark model wave functions and parameter set (quark masses and $\\beta$) used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes the previous in-medium decay constants and distribution amplitudes, and the shifted longitudinal momentum variable used here."},{"cited_title":"Bozkır, A","cited_arxiv_id":null,"evidence_quote":"It provides the quark-meson coupling model mean-field framework that yields the in-medium light-quark mass and energy shifts."},{"cited_title":"Koponen, F","cited_arxiv_id":null,"evidence_quote":"It gives the in-medium pion form factors in a related light-front approach, used as a comparison for the density trend."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives in-medium kaon form factors and the behavior at densities above $1.5\\rho_0$, used as a comparison."},{"cited_title":"Fornetti, E","cited_arxiv_id":null,"evidence_quote":"It provides a recent field-theoretic model calculation of in-medium meson structures, used to cross-check the density dependence."},{"cited_title":"Arrington et al., Physics with CEBAF at 12 GeV and future opportunities, Prog","cited_arxiv_id":null,"evidence_quote":"It supplies free-space lattice QCD form factors at high $Q^2$, used to benchmark the model in free space."},{"cited_title":"Maris and P","cited_arxiv_id":null,"evidence_quote":"It supplies the lattice charge radius of the kaon quoted in the comparison table."},{"cited_title":"Chang, I","cited_arxiv_id":null,"evidence_quote":"It supplies lattice QCD data for D-meson couplings used to benchmark the heavy-light predictions."}],"review_version":1}