{"id":"1cd623f5-b1b7-47b4-bab3-2c500227468c","arxiv_id":"2506.08707","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A hybrid light-front plus quark mean-field calculation predicts density-, temperature-, and isospin-dependent masses, decay constants, and distribution amplitudes for D and D* mesons.","lead":"This paper computes how D and D* mesons change inside dense, neutron-rich nuclear matter, using a two-model pipeline: a quark mean-field model supplies in-medium quark masses and a light-front quark model turns those masses into meson masses, decay constants, and distribution amplitudes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"In-medium results do not re-minimize the LFQM variational parameter beta; since beta enters the mass, decay constant, and DA formulas, the reported shifts may be sensitive to this omission and require a consistency check.","rationale":"The reader's weakest_assumption is that in-medium quark masses can be directly fed into the LFQM while keeping potential parameters fixed at vacuum values. My concern is a concrete version of that assumption: beta, the one variational parameter in the LFQM part, is not re-fitted in medium even though the paper states it is determined by a variational principle. This is not a claim of internal mathematical error in the equations, but an incomplete application of the model's own logic. A concrete test can settle it: re-minimize beta in medium and check whether the results shift. If they do, the central quantitative claims are weaker than presented; if they do not, the concern is resolved and the paper's numbers stand. This supports the reader's CONDITIONAL verdict: the model is plausible and qualitatively interesting, but a key parameter-dependence check is missing.","tokens_in":25123,"tokens_out":1947,"duration_ms":25245,"concrete_test":"For rho_B = rho_0, T = 0, eta = 0 and 0.5, recompute beta for D0 and D+ by minimizing <H_qbarq> with the CQMF in-medium m*_u and m*_d. Then recompute the mass shifts and f*/f ratios with the new beta. If the shifts change by more than ~10% relative to the vacuum-beta values, the reported numbers should be revised; if they change negligibly, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The hybrid scheme computes in-medium quark masses m*_q from CQMF (Eq. 13) and inserts them into the LFQM Hamiltonian (Eqs. 17-19, 28), but the Gaussian width beta in the trial wavefunction (Eq. 26) is kept at its vacuum value from Table II. The paper explicitly states beta is fixed by variational minimization; in medium, m*_q changes, so the variational minimum should shift. Since M*_qbarq (Eq. 28), f*_P (Eq. 31), f*_V (Eq. 32), and the DAs (Eqs. 35-36) all depend on beta, a medium-dependent beta could modify the headline mass shifts (-0.097 GeV for D0) and decay constant ratios (0.936) by an amount that may be comparable to the isospin splitting the paper emphasizes. No sensitivity estimate is given. This is an internal consistency gap rather than an external model disagreement: the calculation does not apply its own variational principle consistently.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the in-medium properties of pseudoscalar D (D0, D+, Ds) and vector D* (D0*, D+*, Ds*) mesons in isospin asymmetric nuclear matter using a hybrid approach that combines the light-front quark model (LFQM) with the chiral SU(3) quark mean field (CQMF) model. In-medium quark masses obtained from the CQMF model are inserted into the LFQM Hamiltonian, and the effective masses, weak decay constants, and distribution amplitudes are computed as functions of baryon density, isospin asymmetry, and temperature. The paper reports significant density-dependent mass shifts and decay-constant suppression for mesons containing u/d quarks, smaller effects for strange mesons, and weak temperature sensitivity. Vacuum masses and decay constants are compared with experimental data and with results from other models.","tokens_in":25341,"tokens_out":5369,"duration_ms":62854,"significance":"The results, if taken at face value, provide phenomenological predictions for open-charm meson modifications in dense matter that are relevant for CBM and PANDA at FAIR, and they extend earlier hybrid LFQM+CQMF studies from light and bottom mesons to D and D* mesons. The paper is built on explicit formulas and gives comparisons with several other approaches. However, the quantitative central claim rests on a variational consistency step that is not carried out in the medium, and on an unstated transferability assumption for the CQMF quark masses. These points need to be addressed before the quoted mass shifts and decay-constant ratios can be regarded as robust predictions of the hybrid framework.","major_comments":[{"comment":"The variational parameter beta is not re-minimized in the medium. The paper states in Sec. II B that beta is fixed by minimizing the expectation value of the QCD-motivated Hamiltonian, and Table II lists only vacuum values. When the in-medium quark masses m*_q from Eq. (13) are inserted into Eqs. (17)-(19) and (25), the Hamiltonian changes, so the variational minimum for beta should shift with density, isospin asymmetry, and temperature. Since beta enters not only the mass formula in Eq. (28) but also the decay constants f*_P and f*_V in Eqs. (31)-(32) and the distribution amplitudes in Eqs. (35)-(36), the reported medium modifications, including the headline shift of -0.097 GeV for the D0 mass at rho_B = rho_0, eta = 0, T = 0, are not the outcome of the model's own variational principle. The authors should recompute beta for each in-medium mass configuration and quantify how the mass shifts, decay-constant ratios, and DAs change. A sensitivity estimate is needed even if the effect turns out to be small.","section":"Sec. II B (Eqs. (26)-(28)), Sec. III A"},{"comment":"The transferability of the CQMF in-medium quark masses to the D-meson bound state is an explicit model assumption that is not stated or tested. The effective quark masses m*_q in Sec. II A are computed for quarks confined inside nucleons via the Dirac equation in Eq. (12), but they are then used as the constituent quark masses of the quark-antiquark pair inside a D meson in the LFQM Hamiltonian, while the confinement parameters (a, b, alpha_s) and the smearing parameter kappa are kept at their vacuum values. This assumes that a quark inside a D meson experiences the same scalar and vector mean fields as a quark inside a nucleon, and that the confining interaction is medium-independent. I ask the authors to state this assumption explicitly and to provide a quantitative check of its robustness, for example by comparing with an alternative in-medium mass prescription or by varying the input m*_q within a reasonable range and reporting the resulting spread in the predicted mass shifts and decay-constant ratios.","section":"Sec. II A (Eqs. (12)-(13)), Sec. II B (Eqs. (17)-(19))"}],"minor_comments":[{"comment":"The caption of Table III reads \"Predicted ground-state mass spectra,\" but the D0 and D0* masses are fitted inputs used to fix m_u/d, m_c, alpha_s, a, b, and beta; only the D_s and D_s* rows are genuine predictions and should be labeled accordingly.","section":"Table III"},{"comment":"The text quotes vacuum masses of 2.008 GeV for D_s and 2.111 GeV for D_s*, which are inconsistent with the values 2.010 and 2.112 GeV listed in Table III.","section":"Sec. III A"},{"comment":"There are several typographical and labeling errors, including \"Additonally,\" \"isosppin asymmetry,\" and the figure panel labels \"B = 0\" and \"B = 3 0\" in Figs. 8 and 9, which should read rho_B = rho_0 and rho_B = 3 rho_0.","section":"Throughout"},{"comment":"The paper does not provide numerical tables of the in-medium quark masses m*_q or of the beta values used for the in-medium calculations; providing these numbers, at least for representative densities, would improve reproducibility and allow readers to check the LFQM inputs.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of a hadronic physics journal and addresses a topic of current experimental interest. The main technical concern is internal consistency: the LFQM variational parameter beta is not re-minimized in the medium. If the authors can demonstrate, by recomputation or a controlled sensitivity test, that a medium-dependent beta leaves the central mass shifts and decay-constant ratios essentially unchanged, the paper would be publishable after a revision. The CQMF-to-LFQM transferability assumption is a legitimate modeling choice, but it should be acknowledged and tested; I would not recommend rejection solely on that basis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent extension of the group's hybrid LFQM+CQMF program to D and D* mesons in isospin asymmetric matter. The new physics is in the vector D* channel and the isospin splitting from the delta field. The machinery is standard for this group, and the paper is honest about what is fitted.\n\nWhat is actually new: previous LFQM studies (Ref. [49]) covered only symmetric matter; here the CQMF in-medium quark masses acquire u-d splitting through the delta mean field, and that feeds through to D0 vs D+ masses, decay constants, and DAs. The predictions are cleanly presented: at rho0 and T=0 the D0 mass shifts by -97 MeV at eta=0 and -85 MeV at eta=0.5, with a smaller effect for Ds. Those qualitative trends are credible within the model.\n\nWhat is good: the formalism is self-consistent, the vacuum mass spectra and decay constants are compared to a broad set of lattice/QCD sum rule/BS results, and the authors do not overclaim. The D* mass trend (decrease at low density, then rise due to the hyperfine term) is physically explained.\n\nWhere I would push: the main internal-consistency gap, which the stress-test note identifies, is real. The Gaussian width beta is fixed at its vacuum variational minimum in Table II, but the same variational principle applied to the in-medium Hamiltonian with m*_q would move beta. Since beta enters the mass, decay constants, and DAs in Eqs. (28), (31)-(36), a medium-dependent beta could shift the headline numbers by a noticeable fraction of the 12 MeV isospin splitting they emphasize. They should either re-minimize beta in medium, or give a sensitivity estimate showing the effect is small. Also, the vacuum D0 and D0* masses are fitted by construction, so the vacuum agreement proves less than it might; the in-medium direction is built in, though the magnitude is a genuine prediction. No error bars, no code, no tabulated in-medium quark masses, so independent checks are harder. The comparison with other models in Section III A is useful but brief; a table of competing predictions would clarify how model-dependent the 100 MeV shift is.\n\nThis is not a paper with a load-bearing flaw. It is a solid model prediction that needs a consistency check before publication. I would send it to a referee; the beta issue is addressable and the results are of interest to the FAIR/CBM community.","headline":"A solid, incremental extension of the group's LFQM+CQMF program to D/D* mesons in isospin asymmetric matter; the vector D* results and delta-field isospin splitting are new, but the fixed variational width beta is an addressable internal-consistency gap.","tokens_in":142,"tokens_out":1439,"would_cite":true,"duration_ms":27149,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper predicts that in isospin-asymmetric nuclear matter the masses, weak decay constants, and distribution amplitudes of $D$ and $D^*$ mesons are significantly modified, with the $D^0$ mass shifting down by about $0.097$ GeV at…","keywords":["D mesons","D* mesons","isospin asymmetric nuclear matter","light-front quark model","chiral quark mean field model","weak decay constants","distribution amplitudes","in-medium hadron properties"],"falsifier":"Measure the $D^0$ mass shift in cold nuclear matter at $\\rho_B = \\rho_0$ and $T \\approx 0$, for example from $D$-mesic nuclei or $D$-meson production in heavy-ion collisions; a shift much smaller than $-0.097$ GeV would contradict the central prediction. A lattice calculation of the in-medium $D$-meson mass and weak decay constant at $\\rho_0$ would provide an independent check.","tokens_in":24887,"feed_emoji":"⚛️","tokens_out":10232,"duration_ms":96005,"temperature":0.7,"pith_summary":"The paper predicts how open-charm $D$ and $D^*$ mesons change inside dense, isospin-asymmetric nuclear matter by feeding in-medium quark masses from a chiral SU(3) quark mean field model into a light-front quark model. It claims that the medium lowers the $D^0$ mass by about $0.097$ GeV at nuclear saturation density and zero temperature, with smaller shifts for mesons containing a strange quark. It also predicts that weak decay constants are suppressed by about 6% for the $u/d$-containing $D$ mesons at saturation density, and that distribution amplitudes shift and split between $D^0$ and $D^+$ when isospin asymmetry is introduced. Baryon density, rather than temperature or isospin asymmetry, is identified as the dominant driver of these medium modifications.","feed_headline":"Charm meson mass drops 97 MeV in dense nuclear matter","feed_subtitle":"Density, not temperature, drives shifts in charm meson masses, decays, and distribution amplitudes.","key_machinery":"The load-bearing object is the in-medium constituent quark mass $m_q^* = -g_q^\\sigma \\sigma - g_q^\\zeta \\zeta - g_q^\\delta I_{3q}\\delta + \\Delta m$, obtained by minimizing the CQMF thermodynamic potential and then used directly as the quark mass in the light-front Hamiltonian for the meson. The meson mass eigenstate is computed variationally with a Gaussian trial wave function of width $\\beta$ through Eq. (28), and the same wave function determines the weak decay constants (Eqs. (31)-(32)) and the distribution amplitudes (Eqs. (35)-(36)). The hyperfine term in Eq. (19), proportional to $\\langle S_q \\cdot S_{\\bar q}\\rangle$, distinguishes pseudoscalar from vector mesons and is the reason their mass-versus-density trends differ. The confinement parameters $a$, $b$, $\\alpha_s$, and $\\kappa$ are kept at their vacuum values; only the quark masses are modified by the medium.","core_discovery":"The paper claims that a hybrid treatment, in which in-medium constituent quark masses from the chiral SU(3) quark mean field model are inserted into the light-front quark model, captures the dominant medium effects on $D$ and $D^*$ mesons. At $\\rho_B = \\rho_0$, $T = 0$ and $\\eta = 0$, the $D^0$ mass shifts downward by $-0.097$ GeV, and at $\\eta = 0.5$ the shift is $-0.085$ GeV. $D^+$ behaves similarly, while $D_s$ shifts by roughly $-0.029$ GeV at $\\eta=0.5$ and $D_s^*$ by $-0.024$ GeV. Vector $D^{0*}$ and $D^{+*}$ masses first drop and then rise with baryon density because the hyperfine term acts with opposite signs for vector and pseudoscalar mesons. Weak decay constant ratios $f_M^*/f_M$ fall with density, reaching $0.936$ for $D^0$ at $\\rho_0$, $T=0$, $\\eta=0$, and the distribution amplitudes develop a visible $D^0$/$D^+$ splitting under isospin asymmetry.","pith_inferences":["If the same mean-field assumption holds outside the nucleon sector, the predicted $-0.097$ GeV shift at $\\rho_0$ could be used to estimate $D$-meson binding energies and the formation threshold of $D$-mesic nuclei.","The predicted splitting of the $D^0$ and $D^+$ distribution amplitudes under isospin asymmetry could be probed in processes sensitive to the light-quark momentum fraction, such as $D$-meson production in asymmetric heavy-ion collisions.","A lattice QCD computation of the $D$-meson mass and decay constant in uniform nuclear matter at $\\rho_0$ would test whether the quark inside a meson feels the same mean field as the quark inside a nucleon; a null result for the decay-constant suppression would point to that assumption as the weak link."],"forward_implications":["At nuclear saturation density the $D^0$ mass shift is predicted to be about $-0.097$ GeV in symmetric matter, a signal that could be searched for in $D$-meson spectral measurements.","Mesons containing a strange quark ($D_s$, $D_s^*$) are nearly unaffected, so isospin asymmetry acts mainly on the doublet members with $u/d$ quarks.","The weak decay constants of the $u/d$-containing $D$ mesons are suppressed by roughly 6% at $\\rho_0$ and by about 15% at $3\\rho_0$, which would alter predicted $D$-meson production and decay rates in dense matter.","Vector $D^*$ mesons show a non-monotonic mass shift with density, initially decreasing and then increasing, a signature of the hyperfine interaction that separates them from their pseudoscalar partners.","Raising temperature from $0$ to $0.15$ GeV partially cancels the medium-induced shifts, so cold dense matter gives the cleanest signal."],"supporting_citations":[{"why":"Supplies the light-front quark model formulas for in-medium decay constants and distribution amplitudes that the calculation is built on.","marker":"[49]"},{"why":"Defines the chiral SU(3) quark mean field Lagrangian and quark-meson couplings that produce the in-medium quark masses.","marker":"[52]"},{"why":"Provides the parameter set and the in-medium quark mass definition used in the CQMF calculation.","marker":"[59]"},{"why":"Provides the spin-orbit wave functions that distinguish pseudoscalar $D$ from vector $D^*$ mesons.","marker":"[48]"},{"why":"Supplies the Gaussian smearing and light-front quark model vacuum fits used as the starting point.","marker":"[67]"},{"why":"Provides the experimental vacuum masses and decay constants used as the baseline in Tables III and IV.","marker":"[68]"},{"why":"Gives QCD sum rule values for masses and decay constants that the vacuum results are compared against.","marker":"[70]"},{"why":"Provides a prior quark-meson coupling model prediction of a negative $D$-meson mass shift used as a comparison benchmark.","marker":"[36]"}],"fun_headline_variants":["D0 mass falls 97 MeV in dense nuclear matter","Baryon density, not temperature, shifts D meson properties","Isospin asymmetry splits D0 and D+ meson behaviors","Strange quark reduces nuclear medium impact on charm mesons","Hybrid model predicts density-driven D and D* meson shifts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the quark inside a $D$ meson experiences the same scalar and vector mean fields as a quark confined inside a nucleon, so the in-medium quark masses from the chiral quark mean field model can be used directly in the meson wave function while the confining interaction parameters stay at their vacuum values.","fun_headline_variants_meta":{"raw":{"variants":["D0 mass falls 97 MeV in dense nuclear matter","Baryon density, not temperature, shifts D meson properties","Isospin asymmetry splits D0 and D+ meson behaviors","Strange quark reduces nuclear medium impact on charm mesons","Hybrid model predicts density-driven D and D* meson shifts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000731,"raw_usage":{"total_tokens":3323,"prompt_tokens":1046,"completion_tokens":2277,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":2190}},"tokens_in":662,"tokens_out":2277,"duration_ms":19243,"temperature":1.0,"reasoning_tokens":2190,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:04:09.872217+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $D^0$ mass shift in cold nuclear matter at $\\rho_B = \\rho_0$ and $T \\approx 0$, for example from $D$-mesic nuclei or $D$-meson production in heavy-ion collisions; a shift much smaller than $-0.097$ GeV would contradict the central prediction. A lattice calculation of the in-medium $D$-meson mass and weak decay constant at $\\rho_0$ would provide an independent check.","supporting_citations":[{"cited_title":"Chhabra and A","cited_arxiv_id":null,"evidence_quote":"Supplies the light-front quark model formulas for in-medium decay constants and distribution amplitudes that the calculation is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the chiral SU(3) quark mean field Lagrangian and quark-meson couplings that produce the in-medium quark masses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the parameter set and the in-medium quark mass definition used in the CQMF calculation."},{"cited_title":"Chhabra and A","cited_arxiv_id":null,"evidence_quote":"Provides the spin-orbit wave functions that distinguish pseudoscalar $D$ from vector $D^*$ mesons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Gaussian smearing and light-front quark model vacuum fits used as the starting point."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental vacuum masses and decay constants used as the baseline in Tables III and IV."},{"cited_title":"Cassing, E","cited_arxiv_id":null,"evidence_quote":"Provides a prior quark-meson coupling model prediction of a negative $D$-meson mass shift used as a comparison benchmark."}],"review_version":1}