{"id":"12f978f3-076b-4605-bc86-e16bf0dd51a8","arxiv_id":"2607.13846","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In all six Cu2NiXY4 kesterites modeled, replacing half the Ni with Mn narrows the band gap while preserving the tetragonal structure.","lead":"Using computer simulations, the authors replaced half the nickel atoms in six copper-nickel sulfur/selenium compounds with manganese. They found that this always shrinks the material's electronic band gap, which could help engineers tune how these cheap semiconductors absorb sunlight.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Non-spin-polarized DFT for Mn2+ leaves the bandgap trend and hybridization claim unsubstantiated; spin-polarized recalculation required.","rationale":"The reader's weakest_assumption correctly identifies the missing magnetic treatment. The paper's Conclusion explicitly states magnetic properties were not considered, and the Methods section provides no ISPIN or spin-polarization settings. For Mn2+ (d5), this is not a minor omission: without spin splitting, the Mn-3d states are positioned incorrectly, directly compromising the stated mechanism (Mn-3d/Cu-3d/S-p hybridization near the band edges) and the resulting bandgap values. The computed narrowing is only 0.02–0.10 eV across the six compounds, while the method's error versus experiment is up to ~0.24 eV (Cu2NiGeS4), so the trend could be within numerical uncertainty. A spin-polarized recalculation with magnetic ordering is the minimal check that would settle whether the trend persists. We agree with the reader's CONDITIONAL verdict; the paper is internally consistent and the pure-compound gaps agree reasonably with experiment, but the central claim is not yet robust. Therefore UNCHANGED verdict is appropriate.","tokens_in":8497,"tokens_out":7428,"duration_ms":71816,"concrete_test":"Re-run the full SCAN geometry optimization and mBJ+U electronic-structure calculation for Cu2Ni0.5Mn0.5SnS4 with spin polarization (ISPIN=2) and with ferromagnetic and antiferromagnetic orderings, using U=5 eV for Cu, Ni, Mn. Compare the resulting lattice parameters, band gap, Mn magnetic moment, and DOS near the Fermi level to the reported unpolarized results. If the band gap differs by more than 0.05 eV, or if Mn-3d states appear in the gap, the paper's central narrowing trend and hybridization mechanism are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's method section never mentions spin polarization, and the Conclusion explicitly states magnetic properties were not considered. Mn2+ (d5) is an open-shell ion; a nonmagnetic treatment artificially pins up- and down-spin occupations, which shifts the Mn-3d states relative to the band edges and alters the d-p hybridization that the authors claim drives the bandgap narrowing. The reported gap changes are all small (0.02–0.10 eV), comparable to the method's deviations from experiment (e.g., Cu2NiGeS4 calc. 2.039 eV vs exp. 1.8 eV). Without spin-polarized calculations and specified magnetic ordering (FM/AFM), the DFT solution may not even be the ground state; the Mn magnetic moment and the positions of the Mn 3d states could change substantially, possibly placing deep in-gap levels and violating the paper's conclusion that semiconducting behavior is preserved.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a DFT study (SCAN for relaxation, mBJ+U for electronic structure) of 50% Ni→Mn substitution in six kesterite compounds Cu2NiXY4 (X=Sn, Ge, Si; Y=S, Se). The authors find that the tetragonal kesterite structure is preserved upon substitution, lattice parameters follow expected ionic-radius trends, and the band gap narrows in all six systems (e.g., Cu2NiSnS4: 1.592→1.487 eV). They attribute the narrowing to hybridization between Mn-3d, Cu-3d, and S/Se-p states near the band edges. The paper is framed as the first systematic comparative study of Mn substitution across this family, and the authors present it as showing a route to band-gap engineering for photovoltaic and optoelectronic applications.","tokens_in":8813,"tokens_out":3555,"duration_ms":34967,"significance":"If correct, the paper provides a systematic dataset for a previously unexplored compositional space in kesterite-type chalcogenides, with a clear and internally consistent trend. The use of SCAN relaxation and mBJ+U is a reasonable state-of-the-art protocol for this class of materials, and the validation of the pure-compound band gaps against experiment (Table 2) gives some confidence in the computational setup. The explicit acknowledgement of limitations (50% substitution only, magnetic properties deferred) is honest. However, the central physical claim—that the band-gap narrowing arises from Mn-3d hybridization—rests on a non-spin-polarized treatment of an open-shell Mn2+ (d5) ion, and the magnitude of the effect (0.02–0.10 eV) is comparable to known DFT band-gap errors. The significance for applications would therefore be considerably strengthened by spin-polarized calculations and a robustness test of the Hubbard U.","major_comments":[{"comment":"The calculations do not specify or include spin polarization for Mn2+, which has a half-filled 3d shell. The Conclusion explicitly defers “magnetic properties … were not considered” to future work. This omission is load-bearing: the reported band-gap narrowing is 0.02–0.10 eV (Table 2), comparable to both the method's deviations from experiment (e.g., Cu2NiGeS4: 2.039 eV calc. vs 1.8 eV exp.) and to the shifts expected from different spin states. A nonmagnetic treatment of Mn2+ artificially pins up/down occupations and shifts the Mn-3d states relative to the band edges, directly affecting the d–p hybridization that the authors claim drives the narrowing. The authors should either repeat the calculations with spin polarization for at least a representative subset (e.g., Cu2Ni0.5Mn0.5SnS4 and Cu2Ni0.5Mn0.5SnSe4) for FM and AFM orderings, or provide a clear justification for why the nonmagn","section":"§2 (Method) and §4 (Conclusion)"},{"comment":"A single Hubbard U=5 eV is adopted for Cu, Ni, and Mn, taken from literature values that were not derived for Mn in a chalcogenide environment. Because the band-gap changes are small and U directly controls the position of Mn-3d states, the quantitative values in Table 2 are not robust against reasonable variations in U. I request a U-dependence test for at least one Mn-substituted compound (e.g., U = 3, 5, 7 eV for Mn) to show that the narrowing trend and the hybridization picture are not artifacts of the chosen parameter. This is particularly important because the paper validates the method only against pure compounds, where Mn is absent.","section":"§2 (Method), Table 2"},{"comment":"The central claim that “Mn substitution leads to a narrowing of the band gap in all studied compounds” is based solely on mBJ+U calculations with no experimental or independent theoretical data for any Mn-substituted composition. The small gap reductions (0.02–0.10 eV) are comparable to the known accuracy limits of the method. While the absence of prior data is acknowledged, the conclusion is stated without an explicit error estimate or a discussion of how sensitive the trend is to the magnetic state and U. I recommend at least a conservative statement of uncertainty and, ideally, a test of the trend with a different functional or with spin-polarized calculations.","section":"§3 (Results and Discussion), Table 2"}],"minor_comments":[{"comment":"The abstract says “the effect of partial substitution of Mn by Ni”; this should read “substitution of Ni by Mn.”","section":"Abstract"},{"comment":"Typo: “VASР” contains a Cyrillic “Р”; should be “VASP.” Also, “Monhost-Pak” should be “Monkhorst-Pack.”","section":"§2 (Method)"},{"comment":"“three-dimensional localized electrons” is unclear; the authors mean “d electrons,” and “self-interference errors” should be “self-interaction errors.”","section":"§2 (Method)"},{"comment":"“primary structure” should be “primitive cell.” There is also a stray “Å” before “compounds” in the sentence “general trend remains the same for all Å compounds.”","section":"§3 (Results and Discussion)"},{"comment":"The column header “This Work Literature mBJ+U Exp Calc” is confusing; separate the columns (e.g., “This work,” “Experimental,” “Other calc.”). Also, the entry for Cu2NiGeSe4 has no literature values; this should be explicitly stated as “No data found.”","section":"Table 2"},{"comment":"The caption lists six panels (a)–(f), but only one DOS plot is shown in the extracted text. Please ensure all six panels are actually included in the figure, with axis labels (energy in eV and DOS in states/eV) clearly readable.","section":"Figure 2"},{"comment":"Reference [14] (ZrO2/YSZ surface adsorption) appears unrelated to the kesterite context; please verify the citation or replace it with an appropriate reference on Cu2NiXY4.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's main finding is a systematic but small band-gap narrowing trend. The lack of any spin-polarized treatment for Mn2+ is a genuine correctness risk, not just a missing nicety, because the proposed mechanism is d–p hybridization and the effect size is comparable to methodological error. The single U value transferred from Cu/Ni literature to Mn also needs a sensitivity test. I would encourage the editor to require these additional calculations before publication, rather than accepting the present version, because the central claim is currently not well-grounded. The paper is otherwise clearly written and the comparative systematic approach is a useful contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a straightforward, first systematic scan of 50% Mn-for-Ni substitution in the Cu2NiXY4 kesterite family, and the dataset is genuinely new. The authors compute SCAN geometries and mBJ+U gaps for six doped compounds and report consistent structural trends — smaller IV cations shrink the lattice, Se expands it — which look right. The bandgap narrowing is uniform, though small (0.02–0.10 eV). If the underlying electronic structure is trustworthy, this is a useful reference set for a niche material family.\n\nThe soft spot is real and load-bearing: nothing in the methods says spin polarization was enabled, and the conclusion explicitly punts on magnetic properties. Mn2+ is d5. A nonmagnetic calculation forces equal up/down occupations, which misplaces the Mn-3d states and can put artificial states near the band edges. The paper’s central mechanistic claim — the gap shrinks because Mn-3d hybridizes with Cu-3d and S/Se-p near the edges — is exactly the kind of claim that depends on that magnetic treatment. The observed shifts are also small enough (e.g., 0.021 eV in Cu2NiSnSe4) to be within the method’s deviation from experiment, so the trend needs spin-polarized confirmation, preferably with a couple of magnetic orderings and a magnetic moment check.\n\nLesser issues: a single U=5 V is used for three different 3d ions, which is convenient but not well justified; validation against experiment rests on three compounds, one of which is off by ~0.24 eV (Cu2NiGeS4, calc 2.039 vs exp 1.8), and is still called good agreement. The DOS figures are only shown for two systems, so the hybridization story is asserted more than demonstrated.\n\nWhat is good: the authors are explicit about the 50% substitution being a modeling choice, about the lack of lower-concentration data, and about magnetic properties being left for future work. That honesty is not common, and it makes the paper a decent starting point. A spin-polarized recalculation could well preserve the qualitative trend, since Mn substitution usually does reduce gaps in these chalcogenides, but the current evidence is not enough to publish the trend as established.\n\nWho should read it: computational people working on kesterite doping or on Mn in tetrahedral chalcogenides. It is not a broad-impact result. My recommendation: send it to peer review with the expectation of heavy revision, specifically asking for spin-polarized calculations and a statement of the magnetic ground state. If the authors can confirm the trend with proper Mn spin treatment, the dataset becomes credible; if not, the conclusion will need to be retracted to a weaker 'geometry-only' claim.","headline":"A genuinely new dataset of Mn-substituted Cu2NiXY4 gaps, but the central trend rests on a non-spin-polarized treatment of Mn2+ and needs a targeted recalculation before it can be believed.","tokens_in":9194,"tokens_out":4147,"would_cite":false,"duration_ms":56049,"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 shows that 50% Mn-for-Ni substitution narrows the band gap in every Cu2NiXY4 (X=Sn, Ge, Si; Y=S, Se) kesterite, from 1.028–3.397 eV to 1.007–3.333 eV, by hybridizing Mn-3d states with Cu-3d and S/Se-p bands near the edges.","keywords":["density functional theory","kesterite","band gap engineering","Mn substitution","electronic structure","Cu2NiXY4","mBJ+U","chalcogenide semiconductors"],"falsifier":"Recompute Cu2Ni0.5Mn0.5SnS4 and Cu2Ni0.5Mn0.5GeSe4 with spin-polarized DFT+U, testing ferromagnetic and antiferromagnetic Mn orders, and compare the band gaps and density of states with the present results. If the band gap changes by more than about 0.1 eV or the gap-narrowing reverses in any compound, the central trend is not robust. Alternatively, measure the optical absorption edge of a well-characterized 50%-Mn-substituted thin film; a redshift matching the predicted ~0.1 eV lowering would confirm the claim.","tokens_in":8432,"feed_emoji":"☀️","tokens_out":5983,"duration_ms":46582,"temperature":0.7,"pith_summary":"Replacing half of the nickel in the kesterite semiconductors Cu2NiXY4 (X = Sn, Ge, Si; Y = S, Se) with manganese is predicted by density functional theory to shrink the band gap in every member of the family, from 1.028–3.397 eV down to 1.007–3.333 eV, while retaining the tetragonal kesterite structure. The mechanism, according to the paper, is hybridization between Mn-3d, Cu-3d, and S/Se-p orbitals near the band edges, which redistributes the density of states around the Fermi level without creating deep mid-gap levels. The central example is Cu2NiSnS4, whose calculated gap falls from 1.592 eV to 1.487 eV, in line with the experimentally known undoped value of about 1.6 eV. The authors frame this as the first systematic comparison of Mn substitution across the whole Cu2NiXY4 family, and as a practical route to tunable band gaps in earth-abundant absorber layers for photovoltaics and optoelectronics.","feed_headline":"Mn doping narrows kesterite band gaps across six compounds","feed_subtitle":"DFT predicts a small but systematic gap reduction in Cu2NiXY4, keeping the absorber structure intact.","key_machinery":"The central mechanism is d-p hybridization around the band edges: impurity Mn-3d states sit close to the Fermi level and mix with Cu-3d and S/Se-p states, shrinking the gap without forming deep mid-gap states. The computational machinery is density functional theory with the SCAN meta-GGA for geometry relaxation and the mBJ+U method with Hubbard U = 5 eV on Cu, Ni, and Mn for the electronic structure. This combination yields band gaps in good agreement with available experimental values for the undoped compounds, e.g., 1.592 eV calculated vs 1.6 eV measured for Cu2NiSnS4, giving confidence that the Mn-induced trend is not an artifact of the functional.","core_discovery":"Mn substitution at the 50% level (one of the two Ni sites per primitive cell replaced by Mn) preserves the tetragonal kesterite structure in every compound studied, while consistently reducing the band gap from 1.028–3.397 eV to 1.007–3.333 eV. The mechanism is electronic rather than purely structural: the localized Mn-3d states appear near the band edges, hybridize with Cu-3d and S/Se-p orbitals, and redistribute the density of states around the Fermi level, lowering the energy difference between valence-band maximum and conduction-band minimum. The Mn–Y bond lengths form a distorted tetrahedron and vary systematically with cation and anion, but the gap-narrowing trend is consistent across","pith_inferences":["The paper does not include spin polarization for Mn2+, which has a half-filled 3d shell; because the gap closing is attributed to Mn-3d states, a magnetic ground state could shift the gap sizes or even the sign of the change.","Since only the 50% substitution level is accessible in the primitive cell, the claimed trend is not yet established for dilute concentrations; a supercell study below 50% would test whether the gap keeps shrinking monotonically or saturates.","The same d-p hybridization argument suggests a testable extension: substituting Mn at the Cu or group-IV site in the same kesterite family should produce a different gap response, which would clarify whether the Ni-site substitution is special.","An experimental falsifier is straightforward: optical absorption measurements on Mn-substituted Cu2NiSnS4 thin films should show a redshifted edge relative to the undoped compound; a blueshift or a change larger than about 0.1 eV would challenge the predicted mechanism."],"forward_implications":["The band gaps of Cu2NiSnS4, Cu2NiGeS4, and Cu2NiSnSe4 after Mn substitution remain in the 1.0–2.0 eV window considered optimal for thin-film solar absorbers.","Because the gap narrowing is systematic across the entire family, Mn doping can serve as a general band-gap tuning strategy for these earth-abundant semiconductors, not just a special case.","The Mn-derived states stay near the band edges rather than forming deep gap levels, so the semiconducting character is retained after substitution.","The preserved tetragonal kesterite structure means the substituted compounds remain compatible with existing device fabrication approaches.","The increased density of states near the band edges could improve visible-light absorption and photoexcited carrier generation, though the paper notes transport properties require separate calculations."],"fun_headline_variants":["Mn substitution narrows bandgaps across six kesterite compounds","DFT predicts Mn doping tunes electronic gaps in Cu2NiXY4 family","Manganese doping preserves structure while cutting bandgap in kesterites","Half-Mn substitution lowers band edges in Cu2Ni chalcogenide kesterites","Cu2NiXY4 bandgaps shrink with Mn, DFT shows"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the calculated Mn-doped band gaps are reliable without specifying the magnetic ordering of the half-filled Mn2+ 3d shell; the paper explicitly defers magnetic properties to future work, so the reported narrowing trend could change under spin-polarized treatment.","fun_headline_variants_meta":{"raw":{"variants":["Mn substitution narrows bandgaps across six kesterite compounds","DFT predicts Mn doping tunes electronic gaps in Cu2NiXY4 family","Manganese doping preserves structure while cutting bandgap in kesterites","Half-Mn substitution lowers band edges in Cu2Ni chalcogenide kesterites","Cu2NiXY4 bandgaps shrink with Mn, DFT shows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000686,"raw_usage":{"total_tokens":3032,"prompt_tokens":913,"completion_tokens":2119,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":2020}},"tokens_in":657,"tokens_out":2119,"duration_ms":15228,"temperature":1.0,"reasoning_tokens":2020,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T03:31:34.923315+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute Cu2Ni0.5Mn0.5SnS4 and Cu2Ni0.5Mn0.5GeSe4 with spin-polarized DFT+U, testing ferromagnetic and antiferromagnetic Mn orders, and compare the band gaps and density of states with the present results. If the band gap changes by more than about 0.1 eV or the gap-narrowing reverses in any compound, the central trend is not robust. Alternatively, measure the optical absorption edge of a well-characterized 50%-Mn-substituted thin film; a redshift matching the predicted ~0.1 eV lowering would confirm the claim.","supporting_citations":[],"review_version":1}