{"id":"b3cb9419-305e-433a-ab7a-855310e48ef4","arxiv_id":"1908.09265","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"DFT predicts tetragonal C8 and B4N4 with diamond-like incompressibility and a layered C8 with a ferromagnetic ground state and magnetization collapse around 12 GPa.","lead":"This paper uses density functional theory to predict new carbon structures: a tetragonal 3D form of octacarbon C8, its boron-nitride analog B4N4, and a layered 2D C8. The 3D forms are claimed to be nearly as stiff as diamond, while the 2D form is predicted to be a soft ferromagnet.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proposed C8 and B4N4 phases are not shown to be dynamically stable; the elastic-constant and cohesive-energy checks in Sections 3.2 and 3.5 do not rule out imaginary phonon modes or competing lower-energy phases, so the central claims of near-diamond hardness and a ferromagnetic ground state…","rationale":"At root, the paper's claims are about new carbon polymorphs that are exceptionally hard (3D) or magnetically ordered (2D). For any such claim to hold, the structure itself must be a stable or at least a metastable local minimum. The authors attempt to establish this through cohesive energies and elastic constants, but both checks are insufficient: cohesive energy relative to isolated atoms does not compare against known carbon polymorphs, and elastic stability does not imply dynamic stability. The paper contains no phonon calculation, the standard test for dynamic stability, and no formation enthalpy relative to graphite or diamond, the standard test for thermodynamic accessibility. The concern is made concrete by the fact that the reported 2D-C8 elastic constants violate the hexagonal identity C66=(C11-C12)/2, so the elastic data are not self-consistent. A phonon calculation will settle the issue directly: if imaginary modes appear, the phases are not real local minima; if all modes are real, the structural claims are substantially strengthened, though hardness and magnetism would still need separate validation. Because the reader's conditional verdict is already keyed to exactly this missing evidence, the stress-test does not move the verdict but reinforces it.","tokens_in":9198,"tokens_out":10660,"duration_ms":103507,"concrete_test":"Compute the phonon dispersion of the three proposed structures—tetragonal C8 (P4/mmm), B4N4, and 2D-C8 (P6/mmm, both nonmagnetic and ferromagnetic)—using density functional perturbation theory or the finite-displacement method as implemented in VASP/Phonopy. If any imaginary-frequency mode exists at any k-point, the structure is dynamically unstable and the central claims collapse. As a secondary numerical check, recompute the elastic constants of 2D-C8 to determine whether C66 equals (C11-C12)/2; if the identity is not satisfied, the reported elastic data are erroneous.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sections 3.2 and 3.5 infer the existence and stability of tetragonal C8, B4N4, and 2D-C8 from cohesive energies and from positive elastic constants that satisfy three mechanical inequalities. These are necessary but not sufficient conditions. No phonon dispersion is calculated, and no formation enthalpy relative to graphite or diamond is reported, leaving open the possibility that the proposed structures are not local minima on the Born-Oppenheimer surface or are thermodynamically inaccessible. This gap is load-bearing because every headline property—near-diamond incompressibility and hardness for 3D-C8, the B4N4 band gap, and the ferromagnetic ground state of 2D-C8—presupposes that the structures are realizable. The elastic-stability argument is further undermined by an internal inconsistency: for a hexagonal lattice, C66 must equal (C11-C12)/2, but the reported 2D-C8 values (C11=477, C12=117, C66=13 GPa) violate this identity, indicating that the elastic data cannot be trusted as a stability proof.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports first-principles DFT (VASP/GGA-PAW) calculations for new octacarbon allotropes: a tetragonal P4/mmm 3D C8, an isoelectronic distorted tetragonal B4N4, and a hexagonal P6/mmm 2D C8. The authors claim that tetragonal C8 is cohesive and has incompressibility near diamond, that B4N4 is also ultra-hard and shows a band gap comparable to diamond, and that 2D-C8 has a ferromagnetic ground state with a magnetization-collapse pressure of 12 GPa, suggesting 'soft' ferromagnetic behavior. Stability is asserted from cohesive energies, Birch-Murnaghan EOS fits, and positive elastic constants satisfying mechanical inequalities. Electronic band structures and ELF maps are used to discuss bonding and spin polarization.","tokens_in":9452,"tokens_out":6733,"duration_ms":66064,"significance":"If the claims were fully substantiated, the paper would contribute new potential carbon allotropes with high incompressibility and a rare example of a 2D carbon ferromagnet, with possible relevance for spintronics and hard materials. The work uses standard DFT settings, and the EOS fits appear internally consistent; the cohesive-energy table provides a useful comparison among the phases. However, the current evidence for structural stability is incomplete, the 'hardness' language overstates what bulk-modulus data can support, and an internal inconsistency in the reported elastic constants of the hexagonal phase undermines the stability argument as presented. Because the headline claims depend on these points, the paper is not yet ready for publication.","major_comments":[{"comment":"For a hexagonal lattice in Voigt notation the elastic stiffness constants must satisfy C66 = (C11 - C12)/2. The reported 2D-C8 values (C11=477, C12=117, C66=13 GPa) violate this identity, since (477 - 117)/2 = 180 GPa, not 13 GPa. This inconsistency directly affects the section's conclusion that the positive elastic constants establish mechanical stability, and it must be corrected or explained.","section":"§3.5"},{"comment":"The paper repeatedly equates hardness with the bulk modulus B0 obtained from EOS fits. Hardness is a distinct mechanical property, typically correlated with shear modulus and plastic deformation resistance, not bulk modulus alone. No shear modulus, Vickers hardness, or other hardness indicator is computed, so the statement that tetragonal C8 has 'hardness close to diamond' is not supported by the presented data and should be rephrased in terms of incompressibility or supplemented with actual hardness estimates.","section":"§3.4 and Abstract"},{"comment":"Structural stability is inferred from cohesive energies and from positive elastic constants satisfying three inequalities. These conditions are necessary but not sufficient: a structure can pass them while still possessing imaginary phonon modes or being thermodynamically inaccessible. The manuscript presents no phonon dispersion calculation and no formation enthalpy relative to graphite or diamond. The cohesive energies in Table 2 are referenced to isolated atoms and do not establish relative phase stability. Without these checks, the proposed structures cannot be considered established local minima on the Born-Oppenheimer surface.","section":"§3.2 and §3.5"},{"comment":"The ferromagnetic ground state of 2D-C8 rests on a total-energy difference of about 0.15–0.2 eV per cell between spin-polarized and non-spin-polarized states. This energy difference is within the typical accuracy of GGA for magnetic ordering energies, so the identification is fragile without a more robust treatment (e.g., a Hubbard U, hybrid functional, or systematic antiferromagnetic ordering search). The AFM calculation, described only as assigning opposite spins to various atoms, is too briefly documented to serve as a decisive check.","section":"§3.1.2 and Fig. 4"}],"minor_comments":[{"comment":"The text calls cubic C8 (Ia-3) 'experimentally identified,' but reference [12] is a theoretical prediction; please correct the wording and the citation in Table 1, which currently points to reference [4].","section":"Abstract and Table 1"},{"comment":"The critical-pressure formula is typeset incorrectly; it should read P_C = (B0/B')[(V0/V1)^B' - 1].","section":"§3.4"},{"comment":"The text cites B0=416 GPa for C8 Ia-3, but the EOS fit in Fig. 3 gives B0=413 GPa; please unify the values.","section":"§3.4 and Fig. 3"},{"comment":"The convergence criterion for forces is given as '0.02 eV/Å^3'; force convergence should be in units of eV/Å.","section":"§2"},{"comment":"The caption is garbled ('...with a magnetization of 1.8states diamond...') and should be rewritten.","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's English and notation require substantial editorial work before resubmission. The central claims are interesting but currently depend on insufficient stability evidence and on an elastic-constant inconsistency. If the authors can supply phonon calculations, formation enthalpies relative to known carbon phases, and corrected elastic data, the paper may become publishable in a specialized condensed-matter journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper proposes three new carbon/BN structures — tetragonal C8, distorted-tetragonal B4N4, and 2D P6/mmm C8 — and characterizes them with standard DFT. The structures are genuinely new as far as I can tell; cubic Ia-3 C8 and Si8 were known, but these variants are not. The EOS fits are internally consistent, the cohesive-energy trends are sensible, and the ELF maps and band structures are presented in enough detail to follow the argument. Credit where due: this is a straightforward, reproducible DFT pipeline applied to plausible compositions, and the author does not appear to have tuned any target property.\n\nThe soft spots are load-bearing. First, the headline claims of near-diamond hardness rest on bulk modulus, not hardness; a high B0 does not imply high hardness, and the paper even conflates the two in the abstract and conclusions. Second, structural stability is inferred only from positive elastic constants and cohesive energies. No phonon calculation is shown, and no formation enthalpy relative to graphite or diamond is reported. That leaves open the real possibility that these structures are dynamically unstable or thermodynamically inaccessible — which matters because every headline property presupposes the structures can actually exist. Third, and more specifically, the hexagonal elastic constants for 2D-C8 violate the identity C66 = (C11−C12)/2: the paper reports C11=477, C12=117, C66=13 GPa, while (C11−C12)/2 = 180 GPa. That is a clear internal inconsistency, so the elastic-stability check for that phase cannot be trusted as reported. Finally, the ferromagnetic ground state of 2D-C8 rests on an energy difference of ~0.15 eV per cell between SP and NSP, which is right at the GGA error bar, and the estimated PC = 12 GPa comes from a formula that looks dimensionally suspect as written.\n\nThe citation pattern is clean; the self-citations are contextual, and the Johnston/Hoffmann C8 precedent is properly cited. The writing has typos and some awkward phrasing, but the content is understandable.\n\nWho is this for? People working on carbon allotropes and 2D magnets might find it a useful pointer, but only after the stability questions are answered. I would not cite it in its current form.\n\nRecommendation: this deserves a serious referee rather than a desk reject — the structures are new and the calculations are, in principle, checkable. But a referee should insist on phonon dispersions, formation enthalpies relative to graphite/diamond, a proper hardness measure, and a corrected elastic-constant table for 2D-C8.","headline":"Plausible new C8/B4N4 structures with overclaimed stability and hardness; the hexagonal elastic constants are internally inconsistent and no phonon or formation-enthalpy check is provided.","tokens_in":9998,"tokens_out":2047,"would_cite":false,"duration_ms":21467,"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":"The paper proposes new carbon phases: a tetragonal 3D C8 nearly as incompressible as diamond, a B4N4 analog, and a 2D C8 calculated to be ferromagnetic.","keywords":["octacarbon","carbon allotrope","ultra-hard material","boron nitride","two-dimensional ferromagnetism","density functional theory","magneto-volume effect","spin chemistry"],"falsifier":"A phonon calculation of tetragonal C8 at its relaxed equilibrium volume would settle the structural claim: imaginary vibrational frequencies would show the phase is dynamically unstable, directly contradicting the stability conclusion.","tokens_in":9002,"feed_emoji":"💎","tokens_out":14776,"duration_ms":122094,"temperature":0.7,"pith_summary":"This paper proposes that carbon can take two new three-dimensional forms and one two-dimensional form with properties beyond those of known allotropes: a tetragonal C8 (space group P4/mmm) calculated to be cohesive and nearly as incompressible as diamond, a boron-nitride counterpart B4N4 with the same electron count, slightly lower stiffness, and a diamond-scale band gap, and a layered hexagonal C8 whose calculated ground state is ferromagnetic with a magnetization of 1.8 $\\mu_\\mathrm{B}$ per cell. A careful reader would care because the predictions place a carbon-only three-dimensional lattice near diamond's bulk modulus and suggest a carbon-based two-dimensional magnet that works without transition metals. The evidence comes from density functional total energies, energy-volume equations of state, elastic constants, electron-localization maps, and spin-polarized band structures; stability is asserted from positive elastic constants and favorable cohesive energies.","feed_headline":"Two new carbon phases: one near-diamond hard, one ferromagnetic","feed_subtitle":"Calculations place 3D C8 and B4N4 near diamond's stiffness and predict a layered C8 whose magnetism collapses at 12 GPa.","key_machinery":"The central machinery is plane-wave density functional theory in the generalized gradient approximation, used to relax each proposed structure, followed by three cross-checks. Birch–Murnaghan energy-volume equations of state (a standard fit used to extract elastic response) provide equilibrium volumes, bulk moduli, and the pressure derivative; finite-distortion elastic constants produce the full stiffness matrix, the Voigt bulk modulus, and the mechanical stability inequalities ($C_{11}>C_{12}$, $C_{11}C_{33}>C_{13}^2$, $(C_{11}+C_{12})C_{33}>2C_{13}^2$); and electron-localization maps visualize the covalent C–C versus ionocovalent B–N character. For 2D-C8 the decisive comparison is between non-spin-polarized and spin-polarized total energies, whose roughly 0.15 eV difference favors the ferromagnetic state; the critical pressure $P_\\mathrm{C} = (B_0/B')[(V_0/V_1)^{B'}-1]$ converts the spin-polarized to non-spin-polarized volume difference into a magnetization-collapse pressure.","core_discovery":"On the paper's own terms, the discovery is that a body-centered tetragonal arrangement of eight carbon atoms—C8 in space group P4/mmm—is a cohesive, ultra-incompressible allotrope: its calculated bulk modulus is about 395 GPa, close to cubic Ia-3 C8 (413 GPa) and diamond (about 429 GPa), and its cohesive energy per atom sits between the two. Substituting boron and nitrogen pairwise for carbon gives a distorted tetragonal B4N4 with ionocovalent B–N bonding, a bulk modulus near 350 GPa, and a band gap of roughly the same magnitude as diamond's. In two dimensions, a P6/mmm C8 sheet built from two interpenetrating carbon substructures—magnetic C1 pairs and a semiconducting honeycomb C2 layer—is calculated to have a ferromagnetic ground state with 1.8 $\\mu_\\mathrm{B}$ per cell (0.9 $\\mu_\\mathrm{B}$ per C1 atom), a large volume increase upon magnetization, and a critical pressure of 12 GPa for the collapse of magnetization, behavior the paper classifies as soft ferromagnetic.","pith_inferences":["Beyond the paper, the small energy difference between the spin-polarized and non-spin-polarized states (about 0.15 eV per cell) suggests the ferromagnetic order of 2D-C8 may be tunable by strain, doping, or gating, not only by applied pressure.","The two-substructure motif—a magnetic metal threaded through a semiconducting honeycomb layer—could serve as a design template for carbon-based spintronic monolayers, though transport properties and exchange coupling are not computed here.","If the P4/mmm motif is stable for carbon, the same isoelectronic substitution logic used to construct B4N4 could be carried to other group-IV pairs, such as silicon and phosphorus analogues, widening the search for hard or magnetic low-dimensional allotropes.","Because the magnetic transition is tied to a large c/a change, growing 2D-C8 on a lattice-mismatched substrate might strain the layer enough to move the ferromagnetic-to-nonmagnetic transition to pressures far below 12 GPa."],"forward_implications":["If tetragonal C8 is real and synthesizable, carbon gains a new three-dimensional allotrope with a bulk modulus near 395 GPa and a cohesive energy closer to diamond than that of cubic C8, making it a candidate hard material.","B4N4 would offer a wide-gap (about 5 eV), hard, ionocovalent counterpart to c-BN, potentially useful where diamond-like hardness and a large electronic gap are desired together.","2D-C8 would be a carbon-only two-dimensional ferromagnet whose magnetization can be switched off by 12 GPa of pressure, a rare combination of magnetism and pressure sensitivity in a light-element layer.","The internal separation into a metallic magnetic C1 substructure and a semiconducting honeycomb C2 substructure implies the same monolayer could carry both magnetic and semiconducting functionality.","The agreement between equation-of-state and elastic-constant bulk moduli supports the reported mechanical trends for all three phases."],"supporting_citations":[{"why":"Supplies the cubic C8 allotrope against which the new tetragonal phase's hardness is calibrated.","marker":"[12]"},{"why":"Supplies the body-centered structural motif that is adapted to build P4/mmm C8 and B4N4.","marker":"[13]"},{"why":"Gives the generalized-gradient exchange-correlation functional used in all total-energy calculations.","marker":"[17]"},{"why":"Provides the prior layered carbon-nitride system whose magnetic instability motivates the spin-polarized treatment of 2D-C8.","marker":"[21]"},{"why":"Supplies the Birch equation of state used to extract bulk moduli and the critical-pressure relation.","marker":"[23]"},{"why":"Provides reference bulk modulus of c-BN and diamond band-gap values used for comparison.","marker":"[24]"},{"why":"Provides the CrO2 and Fe-Ni critical-pressure benchmarks used to classify 2D-C8 as a soft ferromagnet.","marker":"[25]"},{"why":"Supplies the Voigt elastic-constant averaging and the stability inequalities used to assess mechanical stability.","marker":"[26]"}],"fun_headline_variants":["Near-diamond hard C8 and ferromagnetic 2D C8 predicted","Tetragonal C8 rivals diamond; 2D C8 is ferromagnetic","Ferromagnetic 2D C8 collapses at 12 GPa; 3D C8 near-diamond","C8: diamond-hard 3D, ferromagnetic 2D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the proposed structures are stable because their calculated elastic constants are positive and their cohesive energies are favorable; the paper does not compute a vibrational spectrum or a formation enthalpy relative to graphite and diamond, so a dynamically unstable or thermodynamically inaccessible phase would pass its checks.","fun_headline_variants_meta":{"raw":{"variants":["Near-diamond hard C8 and ferromagnetic 2D C8 predicted","Tetragonal C8 rivals diamond; 2D C8 is ferromagnetic","Ferromagnetic 2D C8 collapses at 12 GPa; 3D C8 near-diamond","C8: diamond-hard 3D, ferromagnetic 2D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000959,"raw_usage":{"total_tokens":4128,"prompt_tokens":1032,"completion_tokens":3096,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":648,"completion_tokens_details":{"reasoning_tokens":3006}},"tokens_in":648,"tokens_out":3096,"duration_ms":22059,"temperature":1.0,"reasoning_tokens":3006,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:17:06.866006+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A phonon calculation of tetragonal C8 at its relaxed equilibrium volume would settle the structural claim: imaginary vibrational frequencies would show the phase is dynamically unstable, directly contradicting the stability conclusion.","supporting_citations":[{"cited_title":"Solid State Comm","cited_arxiv_id":null,"evidence_quote":"Supplies the body-centered structural motif that is adapted to build P4/mmm C8 and B4N4."},{"cited_title":"Perdew, K","cited_arxiv_id":null,"evidence_quote":"Gives the generalized-gradient exchange-correlation functional used in all total-energy calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the prior layered carbon-nitride system whose magnetic instability motivates the spin-polarized treatment of 2D-C8."},{"cited_title":"Birch, J","cited_arxiv_id":null,"evidence_quote":"Supplies the Birch equation of state used to extract bulk moduli and the critical-pressure relation."},{"cited_title":"Mattesini, S.F","cited_arxiv_id":null,"evidence_quote":"Provides reference bulk modulus of c-BN and diamond band-gap values used for comparison."},{"cited_title":"Matar, G","cited_arxiv_id":null,"evidence_quote":"Provides the CrO2 and Fe-Ni critical-pressure benchmarks used to classify 2D-C8 as a soft ferromagnet."}],"review_version":1}