{"id":"9c7b564c-dc1c-4484-9d56-b88beb65c11e","arxiv_id":"2608.01213","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the Kek-alpha model, two Dirac cones with different Fermi velocities generate a two-branch particle-hole continuum that shrinks the region of stable plasmons as alpha increases.","lead":"This paper calculates how electrons in a proposed 2D material, the Kek-alpha model, respond to electric fields and what collective oscillations, called plasmons, look like. It finds that two electron speed scales create extra ways for plasmons to decay, so undamped plasmons survive only for weak coupling or very small wave vectors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model validity at α=0.5 and α=1 is the load-bearing assumption: Sec. II says the effective Hamiltonian is accurate only for small α, but the fast-cone dominance and dice-limit plasmon suppression are computed precisely in that regime; if the Hamiltonian is invalid there, the central claim is unsuppo","rationale":"The reader's weakest assumption is exactly the one I find most load-bearing: the paper applies a Hamiltonian it explicitly labels as a small-α approximation to α=0.5 and α=1 and then bases its headline physics on those regimes. My independent reading of Sec. II confirms the quoted caveats, and the central claim cannot survive if the Hamiltonian is quantitatively wrong there. The overlap-factor error in Eq. (18)—the stated result does not follow from the squared modulus of Eq. (19), which contains cos(2δ), not cos(δ)—is a separate correctness problem that further undermines confidence in the numerics, but it is localized and could be fixed without changing the qualitative framework. The model-validity issue is more fundamental because it questions whether the authors are studying the claimed lattice at all for the decisive parameter values. I would keep the reader's CONDITIONAL verdict: the paper needs a direct comparison to the parent tight-binding model at α=0.5 and α=1 before its strong α-dependent conclusions can be accepted. No change to the verdict is therefore needed.","tokens_in":24690,"tokens_out":7250,"duration_ms":76880,"concrete_test":"Reconstruct the tight-binding Hamiltonian of Ref. [65] for the Kek-Y α-T3 lattice with one-sublattice Kekulé distortion and diagonalize it on a finite lattice for α=0.5 and α=1. Compare the two low-energy Dirac velocities and the six low-energy eigenvectors with Eq. (8) and Eqs. (9)–(12). If either velocity deviates by more than about 10% at α=0.5 or α=1, repeat the polarization calculation of Sec. IV using the exact lattice band structure; if the fast-cone dominance or the dice-limit suppression disappears, the abstract's α-dependence conclusion is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central claim—that the fast Dirac cone dominates the particle-hole continuum and that stable plasmons survive only for small α or very small q—rests on the 6×6 Hamiltonian Eq. (1) taken from Ref. [65]. In Sec. II the authors themselves state that the approximation 'only works well for small α≪1' and that 'the precision of this model is much higher for a small α, close to graphene.' Yet the distinctive predictions are presented at α=0.5 and α=1 (Figs. 2, 4, and Appendix C). The existence of a formal α→1 limit of the effective matrix does not establish that the matrix describes the Kek-Y α-T3 lattice at α=1; it only shows the model has a dice-like limit. If the α-dependent fast-cone velocity r_α v_F and the wave functions in Eqs. (9)–(12) are artifacts of the small-α expansion, then the two-branch particle-hole continuum and the enhanced damping at α=0.5 and α=1 are not properties of the lattice model but of an uncontrolled approximation. This is the single most load-bearing concern because it affects both the qualitative band structure and every subsequent response calculation, not just a localized overlap factor.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Kek-Y strained α–T3 (Kek-α) model defined by the 6×6 effective Hamiltonian of Ref. [65]. It derives the band structure (two inequivalent Dirac cones with Fermi velocities v_F and sqrt(1+4α²)v_F, plus two flat bands), lists all 36 wave-function overlap factors, and uses them to compute the dynamical polarization, RPA dielectric function, plasmon dispersions, Landau damping, and static screening. The main claims are that the particle-hole continuum has two branches associated with the two cones, that the fast-cone contribution dominates as α increases, and that undamped plasmons survive only for small α or very small wave vectors. The α→0 limit is benchmarked against graphene.","tokens_in":25036,"tokens_out":4189,"duration_ms":40500,"significance":"If the results are correct, the paper identifies a qualitatively new collective-response feature in a flat-band two-cone Dirac model: a two-branch particle-hole continuum with an α-controlled crossover and a strong suppression of stable plasmons. This is a useful addition to the plasmonics literature on α-T3 and Kek-Y systems. The paper also contains analytically derived projector-operator checks and an explicit graphene-limit verification, which are strengths. However, the central quantitative statements rest on a demonstrably wrong printed overlap factor and on computations at α values for which the adopted effective Hamiltonian is stated to be uncontrolled. These issues must be fixed before the qualitative conclusions can be accepted.","major_comments":[{"comment":"Eq. (18) as printed is algebraically impossible: for α=0.5 and δ_{k,k+q}=π it gives O_{2,2}=(1−1.25)/(1.5)²=−1/9, a negative value for a squared inner product. The error originates in Eq. (19), where the squared modulus of [(1+α²)e^{iδ}+α²e^{-iδ}] is not (1+4α²(1+α²)cosδ). The correct result is [α⁴+(1+α²)²+2α²(1+α²)cos(2δ)]/(1+2α²)². The same problem appears in Appendix C, Eq. (C6), where the dice-limit expression 1/9[1+8cosδ] is also negative at δ=π. Since O_{2,2} enters the polarization sum, the numerical results in Figs. 2–4 may be affected unless the code happened to use the correct expression; the manuscript must be corrected and the calculations re-examined.","section":"Sec. III, Eq. (18)"},{"comment":"The model-validity issue is load-bearing. The paper states in Sec. II that the Hamiltonian from Ref. [65] 'only works well for small α≪1' and that 'the precision of this model is much higher for a small α, close to graphene.' Yet the central claims—fast-cone dominance, enhanced damping, and dice-lattice-limit suppression of plasmons—are presented at α=0.5 and α=1 (Figs. 2 and 4, and Appendix C). A formal α→1 limit of the effective matrix does not establish that the matrix quantitatively describes the Kek-Y α-T3 lattice at large α. The authors should either restrict the quantitative claims to the regime where the Hamiltonian is valid, or provide a validation (e.g., against the parent lattice model or a controlled expansion) showing that the α=0.5 and α=1 response functions are reliable.","section":"Sec. II and Figs. 2, 4"},{"comment":"Because overlap factors are the only nontrivial input to the polarization function beyond the band dispersions, any incorrect overlap changes both the real and imaginary parts of Π⁰ and therefore the plasmon damping boundaries. The manuscript should state which overlap factors were actually implemented in the numerics. If the printed Eq. (18) was not used, the text should be corrected; if it was used, the numerical results and the abstract's claim that plasmons are suppressed for larger α need to be re-derived. This is not a cosmetic typo, since O_{2,2} describes transitions within the α-dependent flat band and contributes to the particle-hole continuum.","section":"Eqs. (20)–(30) and polarization sum"}],"minor_comments":[{"comment":"The caption for Fig. 3 refers to 'upper panels (a), (b) and (c)' and 'three lower plots', but the figure contains only two panels. The caption and panel labels should be reconciled.","section":"Fig. 3 caption"},{"comment":"The notation '6X' in Eq. (13) is not standard; it presumably denotes a summation symbol. Please typeset all sums properly.","section":"Sec. III, text near Eq. (13)"},{"comment":"The phrase 'for a Keck-alpha model' contains a typo ('Keck' should be 'Kek').","section":"Sec. III, last paragraph"},{"comment":"The Hamiltonian in Eq. (5) and Eq. (6) should be explicitly checked against each other: the phase factors e^{±iθ_k} in Eq. (5) should correspond to k±/k in Eq. (6). A brief statement verifying this equivalence would help the reader.","section":"Eqs. (5)–(6)"},{"comment":"The notation 'σ r_ρ^α (ħv_F k)' is confusing because σ is a band index, not a sign factor; please clarify the ordering of the ± signs so that σ=±1 corresponds to conduction/valence bands.","section":"Sec. II, Eq. (8)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a demonstrable error in a central analytic result and a stated model-validity limitation that directly bears on the headline claims. Both are fixable in revision: correct Eq. (18) (and Eq. (C6)), rerun the numerics, and either validate the model at large α or soften the dice-limit conclusions. I would not reject, because the framework and the graphene-limit check are sound, and the qualitative mechanism may survive the correction. However, the authors must clearly report whether the numerical code used the printed or corrected overlap; otherwise the figures cannot be trusted. The manuscript also needs careful copyediting for typos and inconsistent captions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The polarization, plasmon, and static screening calculation for the Kek-α model is genuinely new, and the qualitative picture—two unequal Dirac cones, two-branch particle-hole continuum, enhanced damping from flat-band and inter-cone transitions—is believable. Credit is due for systematically deriving all 36 overlap factors, checking the α=0 graphene limit against the known result, and providing a Green's-function/projector derivation that cross-checks the overlap formalism. The comparison with the unstrained α-T₃ model adds useful context. The Hamiltonian from Ref. [65] is properly credited, and the citation pattern looks honest.\n\nThe soft spots, in order of seriousness. First, Eq. (18) as printed is wrong: for α=0.5 and δ=π it gives a negative value, which a squared inner product cannot do. The line just below shows the origin—cos(δ) appears instead of cos²(δ), and the final expression should be nonnegative. This is a localized algebraic error, but it sits in a central overlap factor and needs to be corrected, with the numerics rechecked where that factor contributes. Second, the stress-test note about model validity is on target. The authors themselves say in Sec. II that the effective Hamiltonian only works well for α≪1. Yet the abstract and the headline results emphasize α=0.5 and α=1. The formal α→1 limit of the matrix does not establish that it correctly describes the dice-lattice limit of the Kek-α lattice; it only shows the model has a dice-like limit. This does not destroy the small-α conclusions, but the strong statements about suppressed plasmons at large α are not yet grounded in a controlled approximation. Third, the Fermi-energy scale E_F⁰ is never properly defined; μ is later set to 1.0E_F⁰ and q is measured in terms of Q⁰=E_F⁰/(ħv_F), but the reader cannot reconstruct what this corresponds to in the band structure. Minor but needs cleanup. Fourth, there is no code or data to reproduce the numerics—annoying for an RPA calculation, not disqualifying.\n\nThe central argument holds up if you read the paper as a small-α calculation with plausible extrapolation. The large-α claims need either a benchmark against the full lattice model or explicit hedging. This paper deserves a serious referee: the errors are localized and the qualitative result is likely right. Send it to peer review, but expect major revision.","headline":"New and mostly plausible RPA response calculation for the Kek-α model, with a real algebraic error in Eq. (18) and a load-bearing validity problem at the large-α values where the physics is claimed to be most distinctive.","tokens_in":704,"tokens_out":762,"would_cite":false,"duration_ms":26923,"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":"In the Kek-α model, the fast Dirac cone takes over the particle-hole continuum as α grows, leaving stable plasmons only at small α and very small wave vectors.","keywords":["Kek-α model","α-T3 lattice","plasmon damping","Landau damping","polarization function","Dirac cones","flat bands","static screening"],"falsifier":"Compute the dynamic polarization from a full tight-binding model of a one-sublattice Kekulé lattice at α=0.5 and compare the imaginary part: if the high-frequency branch (fast cone) does not dominate the response, or if an undamped plasmon appears well above q/k_F ≈ 0.5, the paper's central claim would be contradicted. Alternatively, measure the plasmon dispersion by electron energy loss spectroscopy on a Kek-α sample and look for the predicted rapid onset of Landau damping.","tokens_in":24604,"feed_emoji":"⚡","tokens_out":6325,"duration_ms":61366,"temperature":0.7,"pith_summary":"The paper studies the Kek-α model, a hybrid of Kekulé-distorted graphene and the α-T3 lattice, and works out its collective electronic response. It shows that the low-energy spectrum carries two inequivalent Dirac cones with different Fermi velocities plus two flat bands, and that these give rise to two distinct branches of particle-hole excitations. As the hopping parameter α grows, the fast-cone branch dominates and pushes most plasmon modes into the damping region; stable plasmons survive only for small α or very small wave vectors. The paper also finds that static screening and Friedel oscillations keep the same qualitative form as in the unstrained α-T3 model. If correct, the model provides a tunable platform for controlling plasmon damping by changing α.","feed_headline":"Two Dirac cones squeeze plasmons in Kek-α material","feed_subtitle":"Fast-cone transitions and flat bands damp collective modes, leaving stable plasmons only at small α.","key_machinery":"The central object is the 6×6 low-energy Hamiltonian (Eq. 1) built from two 3×3 blocks and the relative hopping parameter $\\alpha=\\tan\\phi$, whose spectrum gives two degenerate flat bands and two Dirac cones with velocities $v_F$ and $\\sqrt{1+4\\alpha^2}\\,v_F$. The argument runs through the full 36-element set of wave-function overlap factors $O_{\\lambda,\\lambda'}(\\mathbf{k},\\mathbf{q})$, which enter the RPA polarization function; the two distinct velocity scales set the two particle-hole branches, and the flat-band transitions add damping channels. An alternative derivation via band projectors confirms the overlap results.","core_discovery":"Starting from the 6×6 effective Hamiltonian derived in Ref. [65], the paper derives all 36 wave-function overlap factors and evaluates the dynamical polarization function analytically and numerically. The central finding is that the particle-hole continuum splits into two branches set by the two Dirac-cone velocities $v_F$ and $\\sqrt{1+4\\alpha^2}\\,v_F$. Transitions involving the fast cone, including those with the flat bands, create a broad damping region above the main diagonal; as α increases this fast-cone contribution dominates the response, so the zero of the dielectric function $\\epsilon(q,\\omega|\\alpha)=1-v_C(q)\\Pi^{(0)}(q,\\omega|\\alpha)$ that defines the plasmon lies inside the Landa","pith_inferences":["A testable extension is to compute the same polarization function from a full tight-binding model at α=0.5; if the fast-cone branch does not dominate, the model's large-α predictions would need revision.","The two-velocity splitting mechanism is generic; any Dirac material with unequal cone velocities, e.g., strained graphene, should show similar split particle-hole branches and enhanced damping.","The non-monotonic α dependence of the plasmon frequency suggests that a Kek-α sample could be tuned through a maximum collective response by small changes in the Kekulé coupling strength."],"forward_implications":["For α up to roughly 0.5, undamped plasmon modes exist only for wave vectors well below the Fermi wave vector, making the Kek-α model a tunable absorber or damper at moderate α.","The position of the kink in the static polarizability (and hence Friedel-oscillation decay ~1/r⁴) survives strain, so Kek-α screening behaves like α-T3, not like graphene.","The α-dependence of the plasmon frequency is non-monotonic, unlike the unstrained α-T3, so the same material could show a minimal or maximal collective response depending on α.","The two-peak structure of the polarization function offers a direct spectroscopic signature to identify the fast and slow cones in experiments, for instance through electron energy loss or optical absorption.","The graphene limit α=0 checks against known results, validating the overlap-factor machinery for the new model."],"supporting_citations":[{"why":"Supplies the effective 6×6 Hamiltonian that defines the Kek-α model and its band structure.","marker":"[65]"},{"why":"Provides the conventional α-T3 Hamiltonian and pseudospin-1 matrices that the present model generalizes.","marker":"[36]"},{"why":"Defines α=tan φ as the tunable hopping parameter interpolating between graphene and the dice lattice.","marker":"[33]"},{"why":"Supplies the RPA polarization-function formalism for Dirac materials used as the method baseline.","marker":"[68]"},{"why":"Gives the graphene polarization function and particle-hole continuum used to check the α=0 limit.","marker":"[72]"},{"why":"Provides the static-screening and Friedel-oscillation results for unstrained α-T3 that the paper compares against.","marker":"[40]"}],"fun_headline_variants":["Two Dirac cones squeeze plasmons in Kek-α model","Fast-cone damping shrinks stable plasmon region in Kek-α","Kek-α stable plasmons only at small α, flat bands block","Plasmon damping by fast cone narrows Kek-α stability window","Kek-α fast cone and flat bands kill stable plasmons"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"In Sec. II the authors state that the effective Hamiltonian is accurate mainly for small α ≪ 1, close to graphene; the central claim about suppressed plasmons at larger α, including the dice-limit α=1, depends on that Hamiltonian remaining quantitatively valid there.","fun_headline_variants_meta":{"raw":{"variants":["Two Dirac cones squeeze plasmons in Kek-α model","Fast-cone damping shrinks stable plasmon region in Kek-α","Kek-α stable plasmons only at small α, flat bands block","Plasmon damping by fast cone narrows Kek-α stability window","Kek-α fast cone and flat bands kill stable plasmons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000142,"raw_usage":{"total_tokens":1030,"prompt_tokens":796,"completion_tokens":234,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":153}},"tokens_in":540,"tokens_out":234,"duration_ms":3173,"temperature":1.0,"reasoning_tokens":153,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T00:27:21.392546+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the dynamic polarization from a full tight-binding model of a one-sublattice Kekulé lattice at α=0.5 and compare the imaginary part: if the high-frequency branch (fast cone) does not dominate the response, or if an undamped plasmon appears well above q/k_F ≈ 0.5, the paper's central claim would be contradicted. Alternatively, measure the plasmon dispersion by electron energy loss spectroscopy on a Kek-α sample and look for the predicted rapid onset of Landau damping.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the effective 6×6 Hamiltonian that defines the Kek-α model and its band structure."},{"cited_title":"Hwang and author S","cited_arxiv_id":null,"evidence_quote":"Supplies the RPA polarization-function formalism for Dirac materials used as the method baseline."},{"cited_title":"Wunsch , author T","cited_arxiv_id":null,"evidence_quote":"Gives the graphene polarization function and particle-hole continuum used to check the α=0 limit."}],"review_version":1}