{"id":"453a0e86-25d4-47f2-a49b-efa9fdde2bf5","arxiv_id":"2412.03360","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"DFT simulations predict low exfoliation energies (0.24 to 0.28 J/m2) and moderate tensile strengths (5.6 to 10.2 GPa) for Nb3X8 (X=Cl, Br, I) monolayers.","lead":"Using computer simulations of atoms, this paper calculates how strong and how easy to peel apart three new two-dimensional crystals made of niobium and halogen atoms. The results suggest these sheets are fairly strong and can be separated from their parent crystals more easily than graphene.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dynamic stability claim rests on MTP phonons without per-system DFT validation; a missed imaginary mode would invert the central claim.","rationale":"The reader's identified thickness-convention issue is real: the absolute GPa values depend on an effective monolayer thickness that is not first-principles, so the abstract's numbers are convention-dependent while the halogen-size trends are not. However, I see the MTP-based phonon calculation as the more load-bearing concern because it underpins the qualitative claim of dynamical stability, which is one of the paper's headline results. A missed imaginary mode would not just rescale numbers; it would change the physical conclusion. The paper does provide some independent support: structural lattice constants agree closely with earlier DFT+U results [15], the exfoliation-energy protocol is plausible, and the mechanical trends are internally consistent. The MTP method is also established in the authors' prior work. But because no per-system validation or reproducible data are supplied, the stability claim cannot currently be independently checked from the manuscript alone. My proposed direct DFT phonon calculation is a concrete, low-cost check: if it confirms the MTP dispersions, the central claim holds; if it reveals an imaginary mode, the paper's conclusion is overturned. I therefore keep the verdict CONDITIONAL, with the condition being direct DFT verification of the phonon results.","tokens_in":11006,"tokens_out":3897,"duration_ms":43452,"concrete_test":"Compute the phonon dispersion of the Nb3Cl8 monolayer directly with DFT using the same PBE+D3 settings (plane-wave cutoff 300 eV, 5x5x1 k-grid, 4x4x1 supercell, finite-displacement method with 0.01 Å displacements). Compare the branches at Gamma, M, and K with Fig. 2. If any imaginary frequency appears or the low-frequency branches differ substantially, the MTP-based stability claim is not supported; if no imaginary modes appear and the branches match within a few wavenumbers, the claim is confirmed post hoc. Repeating for Nb3I8 would strengthen the check because the heavier system is closest to any soft-mode boundary.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The first central claim is that Nb3Cl8, Nb3Br8 and Nb3I8 monolayers are dynamically stable. Section 2 states that phonon dispersions were obtained with moment tensor potentials fitted to DFT data and evaluated with PHONOPY, not by direct DFT phonon calculations. The only validation offered is a citation to prior work [22,24] and a qualitative statement that the dispersions agree with Jiang et al. [15]. No per-system metrics are given: no MTP-vs-DFT force/energy errors for these specific compounds, no comparison of the MTP dynamical matrices with DFT data, and no raw input or output files. Phonon instabilities are precisely the kind of soft, anharmonic features that surrogate potentials can over-stabilize or miss, especially when fitted to a limited configuration set. If any imaginary branch exists for one of the three monolayers, the qualitative conclusion changes from 'dynamically stable' to 'unstable,' which is a central advertised result. The mechanical and exfoliation trends could still be correct, but the stability claim is the least secured part of the paper because its evidence is indirect.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports spin-polarized DFT-D3 calculations of the structural parameters, phonon dispersions, exfoliation energies, and uniaxial stress-strain response of monolayer Nb3Cl8, Nb3Br8, and Nb3I8. The lattice constants are in good agreement with earlier DFT+U results. The authors find no imaginary phonon frequencies for the three monolayers, based on moment tensor potential (MTP) fitted phonon dispersions; exfoliation energies of 0.24, 0.27, and 0.28 J/m2; isotropic in-plane elasticity with elastic moduli of 98, 89, and 77 GPa; and anisotropic ultimate tensile strengths of 10.2/8.1, 8.1/7.0, and 5.9/5.6 GPa along the zigzag/armchair directions. All mechanical and phonon-frequency indicators decrease as the halogen mass increases.","tokens_in":11161,"tokens_out":7337,"duration_ms":71827,"significance":"The exfoliation and mechanical-response calculations are direct DFT total-energy and stress calculations, and the systematic halogen-size trends are a useful, falsifiable prediction for this recently fabricated kagome family. If the MTP-based phonon validation is strengthened, the paper would provide a compact reference dataset for Nb3X8 monolayers. The two main caveats are that the dynamic-stability claim is currently supported only by a surrogate potential with no in-paper quantitative validation, and that the GPa-scale quantities depend on a non-unique monolayer-thickness convention. Both issues are fixable and should be addressed before the quantitative abstract claims are taken as first-principles.","major_comments":[{"comment":"The dynamical-stability conclusion is derived from phonon dispersions computed with moment tensor potentials rather than from DFT phonon calculations. The text reports no per-system validation: there are no MTP force/energy errors, no training-set statistics, no comparison of MTP dynamical matrices with DFT, and the agreement with Jiang et al. [15] is only qualitative. A missed imaginary branch in any of the three compounds would invert the headline stability claim. Please add direct DFT phonon calculations (finite-displacement or DFPT) for the three monolayers, or at a minimum a quantitative MTP-vs-DFT comparison of forces and phonon frequencies, and show an overlay of the dispersions against [15].","section":"§2, Fig. 2"},{"comment":"The absolute elastic moduli and tensile strengths are obtained by converting 2D in-plane stresses to GPa using an effective monolayer thickness defined as the distance between boundary halogen atoms plus their effective van der Waals diameter. This thickness is a modeling convention; alternative reasonable definitions (bulk interlayer spacing, electron-density extent, or the vdW diameter alone) would rescale all reported GPa values by tens of percent. Please report the primary 2D stiffness and strength in N/m, and present GPa values only together with an explicit statement that they are convention-dependent. This affects the numerical values quoted in the abstract and conclusions, although the halogen trends would likely survive.","section":"§3, Fig. 4"}],"minor_comments":[{"comment":"The phrase 'The plane wave and self-consistent loop cutoff energies were set to 300 and 10-5 eV' conflates the plane-wave energy cutoff with the electronic self-consistency criterion; please distinguish the 300 eV cutoff from the 10-5 eV convergence threshold.","section":"§2"},{"comment":"The deformed Nb3Cl8 panels are listed as (d), (e), and (e); the last panel should be labeled (f) to match the in-text references to Fig. 4e and 4f.","section":"Fig. 4 caption"},{"comment":"The compressed notation 'C11(C12) ... 105 (27) GPa' should be written explicitly, for example C11 = 105 GPa and C12 = 27 GPa for Nb3Cl8, so that the elastic constants are unambiguous.","section":"§3"},{"comment":"The armchair and zigzag loading directions are not defined in the text; please add a small schematic or a sentence relating these directions to the hexagonal lattice vectors.","section":"§3"},{"comment":"The MTP training details (training-set size, active-learning iterations, interaction order, and displacement settings used with PHONOPY) are missing; providing them, or a data-availability statement with the trained potentials, would make the phonon results reproducible.","section":"§2"}],"recommendation":"major_revision","confidential_remarks":"The phonon-validation chain leans heavily on self-citations [22,24] for the accuracy of the MTP approach; this is not disqualifying, but it makes the requested per-system validation more important. The manuscript also lacks a data availability statement, which is increasingly expected for MLIP-based results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read for you on Mortazavi et al. The genuinely new stuff is the cleavage energies (0.24–0.28 J/m2) and the full uniaxial stress-strain curves with failure modes for the three Nb3X8 monolayers. Those weren't in Jiang et al.'s earlier DFT+U paper, and they're consistent with the experimental success exfoliating Nb3Cl8. The lattice constants match prior work closely, and the paper gives the relaxed structures in the text, which is a reproducibility plus.\n\nThe main soft spot is exactly the one the stress-test flags: the phonon dispersions used to claim dynamical stability come from MTP surrogate potentials, with no per-system DFT validation. The authors cite their own prior MTP benchmark paper and say the dispersions look like Jiang's, but there are no force/energy errors for these specific compounds and no direct DFT phonon comparison. If any of the three monolayers had a hidden imaginary branch, the headline stability claim flips. I'd want the authors to run DFT phonons for at least one system, or provide MTP-vs-DFT validation, before I'd trust that specific claim.\n\nSecond issue, also fair: the GPa values for stress depend on an effective monolayer thickness (halogen vdW diameter). That's a modeling convention, not first-principles. The trends across Cl→Br→I are robust, but the absolute moduli and strengths could shift by tens of percent under a different convention. The abstract presents them as if unambiguous.\n\nMinor: no convergence tests or input files beyond the structures. That's increasingly a norm for computational papers, and the lack is a small dent.\n\nOverall, I think the reader's conditional verdict is right. The exfoliation energy and tensile trends are likely correct and practically useful—they tell experimentalists that Br and I monolayers are worth trying to peel. The phonon stability claim and the absolute stress numbers need a firmer basis.\n\nFor peer review: yes, send it out, but ask for MTP validation and a statement on thickness sensitivity.","headline":"Useful DFT add-on for a hot 2D family—new exfoliation and tensile numbers are plausible, but the phonon stability claim rests on surrogate potentials and the absolute stress scale is convention-dependent.","tokens_in":11720,"tokens_out":1652,"would_cite":true,"duration_ms":15909,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.15.Mb","62.20.-x","68.35.Gy"],"model":"deepseek-v4-flash","headline":"Spin-polarized density functional theory predicts that Nb3Cl8, Nb3Br8, and Nb3I8 monolayers are dynamically stable, easily exfoliable (0.24–0.28 J/m2), and have isotropic elastic moduli of 77–98 GPa with anisotropic tensile strengths of…","keywords":["Nb3Cl8","kagome lattice","two-dimensional materials","density functional theory","exfoliation energy","elastic modulus","tensile strength","phonon dispersion"],"falsifier":"Indentation of a suspended Nb3Cl8 monolayer: extracting an in-plane modulus that, when converted with the paper's 6.16 Å thickness, deviates substantially from 98 GPa would falsify the quantitative prediction.","tokens_in":10785,"feed_emoji":"🧊","tokens_out":7647,"duration_ms":69713,"temperature":0.7,"pith_summary":"Using spin-polarized density functional theory, this paper predicts the structural, vibrational, exfoliation, and mechanical properties of three two-dimensional kagome crystals: Nb3Cl8, Nb3Br8, and Nb3I8. The central claim is that all three monolayers are dynamically stable, that their cleavage energies (0.24, 0.27, and 0.28 J/m2) are lower than graphene's 0.37 J/m2, and that they behave as isotropic-elastic but anisotropically strong materials, with elastic moduli from 98 down to 77 GPa. These results matter because a monolayer of Nb3Cl8 with topological flat bands has already been made by exfoliation, so the same route should work for the bromide and iodide cousins, and the predicted stiffness and strength offer concrete numbers for nanodevice design. The paper also finds that replacing Cl with Br and then I systematically lowers stiffness, strength, and phonon group velocities.","feed_headline":"Kagome Nb3X8 sheets predicted stable, peelable, 77–98 GPa","feed_subtitle":"DFT maps exfoliation and strength for Cl, Br, and I versions of the Nb3X8 kagome monolayers.","key_machinery":"The workhorse is a chain of first-principles and machine-learning calculations: spin-polarized DFT (generalized-gradient approximation plus dispersion correction) for geometry and stresses; moment tensor potentials, a class of machine-learned interatomic potentials, fitted to DFT data to evaluate phonon dispersions over supercells; a cleavage-energy procedure that gradually separates one layer from a six-layer slab to compute exfoliation energy; and uniaxial tensile simulations in which true stress is obtained by dividing force by the deformed real volume, using an effective monolayer thickness defined as the distance between boundary halogen atoms plus their van der Waals diameter. That thickness convention is what turns the computed in-plane forces into gigapascals, so it carries the absolute calibration of all reported moduli and strengths.","core_discovery":"The paper establishes, by spin-polarized DFT with a van der Waals correction, that isolated Nb3X8 monolayers are stable in their kagome geometry, since phonon dispersions contain no imaginary frequencies. It then quantifies layer separation: exfoliation energies of 0.24, 0.27, and 0.28 J/m2 for Cl, Br, and I, all below the 0.37 J/m2 benchmark of graphene, indicating weak interlayer coupling and practical mechanical exfoliation. In the monolayer plane, the elastic moduli come out at 98, 89, and 77 GPa, independent of loading direction, while ultimate tensile strengths under uniaxial load are 10.2 (zigzag) and 8.1 (armchair) GPa for Nb3Cl8, 8.1 and 7.0 GPa for Nb3Br8, and 5.9 and 5.6 GPa for Nb3I8. The armchair direction fails more abruptly, with one Nb–X bond elongating far more than the others and with greater thickness reduction, explaining the anisotropic strength and more brittle armchair fracture.","pith_inferences":["The absolute gigapascal values inherit the paper's thickness convention; expressing the same results as in-plane force per unit length (N/m) would give convention-independent numbers and preserve the Cl→Br→I ordering, but would change how these sheets compare with graphene and other 2D materials on an absolute scale.","Because heavier halogens lower both stiffness and phonon group velocity, the same trend suggests lower lattice thermal conductivity in Nb3I8 than in Nb3Cl8; computing or measuring thermal transport would test this extension.","If the topological flat bands of Nb3Cl8 persist under strain, the predicted anisotropy suggests that uniaxial strain along the zigzag direction could tune the electronic structure over a wider range before fracture than armchair strain.","The cleavage-energy curves could also be used to estimate interlayer sliding barriers or the work of adhesion for heterostructures, since the same dispersion-corrected DFT setup would give the needed energy landscape."],"forward_implications":["Nb3Br8 and Nb3I8 monolayers are predicted to be isolable by mechanical exfoliation, with cleavage energies within 0.03 J/m2 of the already-exfoliated Nb3Cl8.","In the linear regime, all three monolayers have isotropic elasticity, so measured in-plane stiffness should not depend on the loading direction.","Under uniaxial tension, the zigzag direction is stronger than the armchair direction in all three compounds, and armchair loading shows a sharper post-peak stress drop, indicating more brittle failure.","Heavier halogens systematically reduce elastic modulus, tensile strength, and phonon group velocities across the Nb3X8 family.","Phonon dispersions free of imaginary frequencies confirm the dynamical stability of all three free-standing monolayers."],"supporting_citations":[{"why":"Reports the experimental exfoliation of Nb3Cl8 monolayers and their topological flat bands, motivating the study.","marker":"[8]"},{"why":"Provides DFT+U lattice constants and elastic moduli for the three Nb3X8 monolayers, used as comparison benchmarks.","marker":"[15]"},{"why":"Gives the cleavage energy of graphene, the baseline against which the predicted exfoliation energies are judged.","marker":"[27]"},{"why":"Introduces moment tensor potentials, the class of machine-learned interatomic potentials used for phonon calculations.","marker":"[21]"},{"why":"Demonstrates the accuracy of moment tensor potentials for phonon dispersions in 2D materials, validating the method used here.","marker":"[22]"},{"why":"Supplies the phonon post-processing routine used to obtain dispersion relations from force constants.","marker":"[23]"}],"fun_headline_variants":["Stable kagome Nb3X8 nanosheets: exfoliation beats graphene","DFT: Nb3X8 kagome sheets stable, exfoliation beat graphene","Exfoliate Nb3Cl8, Nb3Br8, Nb3I8: stable kagome 2D crystals","Kagome Nb3X8 films: elastic 77-98 GPa, exfoliable","2D Nb3X8 kagome: stable, peelable, elastic 77-98 GPa"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the effective monolayer thickness used to convert force to gigapascals—the distance between boundary halogen atoms plus their van der Waals diameter—is the correct thickness to use for real volume, since any other convention rescales every reported modulus and strength by the same factor.","fun_headline_variants_meta":{"raw":{"variants":["Stable kagome Nb3X8 nanosheets: exfoliation beats graphene","DFT: Nb3X8 kagome sheets stable, exfoliation beat graphene","Exfoliate Nb3Cl8, Nb3Br8, Nb3I8: stable kagome 2D crystals","Kagome Nb3X8 films: elastic 77-98 GPa, exfoliable","2D Nb3X8 kagome: stable, peelable, elastic 77-98 GPa"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001096,"raw_usage":{"total_tokens":4628,"prompt_tokens":1051,"completion_tokens":3577,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":3449}},"tokens_in":667,"tokens_out":3577,"duration_ms":22401,"temperature":1.0,"reasoning_tokens":3449,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:28:38.729052+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Indentation of a suspended Nb3Cl8 monolayer: extracting an in-plane modulus that, when converted with the paper's 6.16 Å thickness, deviates substantially from 98 GPa would falsify the quantitative prediction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the experimental exfoliation of Nb3Cl8 monolayers and their topological flat bands, motivating the study."},{"cited_title":"Jiang, Q","cited_arxiv_id":null,"evidence_quote":"Provides DFT+U lattice constants and elastic moduli for the three Nb3X8 monolayers, used as comparison benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the cleavage energy of graphene, the baseline against which the predicted exfoliation energies are judged."}],"review_version":1}