{"id":"e855b259-a6cd-4be9-b86d-b7a92732d37c","arxiv_id":"2505.04804","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"HOP-graphene is predicted by DFT to be a high-capacity anode for Li/Na-ion batteries, with capacities of 1338 and 1227 mAh/g and diffusion barriers of 0.70 and 0.39 eV.","lead":"A computational study predicts that HOP-graphene, a 2D carbon allotrope with 5-6-8-membered rings, can store lithium and sodium ions with capacities of 1338 and 1227 mAh/g, far above graphite. The results suggest a candidate anode material for next-generation batteries, though experimental synthesis and testing are still needed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The -0.5 THz imaginary mode near Gamma in Fig. 1(b) is dismissed with a screening threshold rather than a convergence test; if physical, the stability premise for all subsequent battery claims fails.","rationale":"The reader identified the imaginary phonon mode as the weakest assumption, and I agree. This is the only concern whose truth would falsify the central claim: an unstable host invalidates the capacity, OCV, and diffusion results computed on the idealized symmetric lattice. I also inspected the Arrhenius diffusion coefficients in Sec. 3.3 and Fig. 8; applying Eq. (3) with the stated barriers, a typical jump length, and nu0 = 1e13 Hz gives values several orders of magnitude lower than those reported, suggesting a unit or prefactor error in the diffusion-coefficient numbers. This is a genuine reproducibility issue, but it is secondary because corrected coefficients would still leave Li and Na diffusion faster than the graphite reference values cited in the paper. The paper's acknowledged limitation in Section 3.1 is therefore the load-bearing point, and the existing CONDITIONAL verdict should stand, conditioned on the phonon convergence check and, as a lesser recommendation, release of input files and data.","tokens_in":14877,"tokens_out":15774,"duration_ms":152069,"concrete_test":"Recompute the HOP-graphene phonon dispersion with PHONOPY (or equivalent) on a 2x2 supercell using the same PBE+D2/VASP setup, enforce the acoustic sum rule, and repeat with a second PAW pseudopotential; if cost permits, also test a 3x3 supercell. If the -0.5 THz imaginary mode at Gamma vanishes, the stability premise is cleared. If it persists, displace the atoms along that mode, relax all degrees of freedom, and compare total energies; a lower-energy distorted structure would invalidate all subsequent battery calculations based on the symmetric lattice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim treats HOP-graphene as a viable anode host, which requires the monolayer to be dynamically stable. Section 3.1 and Fig. 1(b) report a -0.5 THz imaginary mode near Gamma, and the authors dismiss it by citing Wang et al. [55]'s -2 THz screening cutoff. That cutoff is a high-throughput operational threshold, not a proof that a Gamma-point acoustic mode is a numerical artifact; no supercell-size, pseudopotential, or acoustic-sum-rule convergence tests are reported. A real mode of this kind would mean the symmetric structure is a saddle point rather than a local minimum, so the adsorption energies, NEB barriers, capacities, and OCVs computed on that structure do not describe a stable ground-state host. The 10 ps AIMD at 300 K is only weak supporting evidence, because a shallow soft-mode distortion can take longer than 10 ps to appear, and the AIMD was run on the same idealized cell. The concern is not that the material is definitely unstable; the original Mandal et al. prediction and the small mode magnitude point the other way. The issue is that the paper's own evidence for the load-bearing stability premise is incomplete, and the cited -2 THz threshold is being used as a substitute for a dedicated convergence test.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a first-principles study of HOP-graphene, a predicted 2D carbon allotrope with 5-6-8-membered rings, as an anode material for Li- and Na-ion batteries. Using PBE+DFT-D2, the authors compute the phonon dispersion, mechanical properties, adsorption energies, Bader charges, diffusion barriers via NEB, and open-circuit voltages. They report theoretical capacities of 1338 mAh/g for Li and 1227 mAh/g for Na, diffusion barriers of 0.70 eV (Li) and 0.39 eV (Na), and average OCVs of 0.42 V and 0.33 V, concluding that HOP-graphene is a promising high-capacity, safe anode.","tokens_in":15143,"tokens_out":7368,"duration_ms":67522,"significance":"If the reported properties are correct, HOP-graphene would compare favorably with previously proposed 2D carbon anodes, offering capacities more than three times that of graphite for Li and much higher for Na, while maintaining low diffusion barriers and safe OCV. The DFT calculations are performed with standard protocols, and the central quantities are derived directly from total energies with no fitting to target properties; the only empirical input is the attempt frequency in the Arrhenius diffusivity, which does not affect the barriers. The main weakness is the incomplete evidence for dynamic stability of the host structure, which underpins the adsorption, diffusion, capacity, and OCV calculations.","major_comments":[{"comment":"The -0.5 THz imaginary frequency near the Gamma point is dismissed by citing the -2 THz screening threshold of Wang et al. [55]. That threshold is an operational high-throughput criterion, not a demonstration that a Gamma-point acoustic mode is a numerical artifact. No convergence tests with respect to supercell size, k-point sampling, pseudopotential, or acoustic sum rule are reported. Because the dynamic stability of the pristine monolayer is the foundation for all subsequent adsorption, diffusion, capacity, and OCV calculations, the authors should provide a targeted phonon investigation (e.g., larger supercell, finite-displacement method with acoustic sum correction) to establish that the mode is an artifact. As written, the manuscript does not rule out the possibility that HOP-graphene is a saddle-point structure.","section":"Section 3.1, Fig. 1(b)"},{"comment":"The AIMD simulations used as evidence of thermal stability are inconsistently described: Section 2 states a total simulation time of 5 ps, Section 3.1 reports 10 ps for the pristine monolayer, and Fig. 12 indicates 5 ps for the fully loaded systems. Moreover, a 5-10 ps trajectory is short relative to the period of a soft mode at -0.5 THz (about 2 ps) and cannot exclude a slow structural distortion. The authors should unify the reported times and, if the phonon convergence tests are inconclusive, extend the AIMD or perform replica exchange to justify the stability conclusion.","section":"Sections 2 and 3.1/3.4"},{"comment":"The capacity calculation is not fully specified: the text reports a maximum loading of 24 Li and 22 Na atoms but never states the supercell size or the number of carbon atoms used for the host weight. The notation Li24C40 implies a 2x2 supercell of the 10-atom unit cell, but this is not stated explicitly, and the text also writes Na24C40 while reporting 22 Na atoms. Please specify the supercell, correct the Na formula, and confirm that the capacity values correspond to the stated loading.","section":"Section 3.4, Eq. (4)"}],"minor_comments":[{"comment":"The extrema for Young's modulus and Poisson's ratio appear to be interchanged (313.45 vs. 347.67 N/m; 0.15 vs. 0.19), and the sentence describing negligible variation followed by exceptionally high anisotropy is contradictory. Please correct the numbers and wording.","section":"Section 3.1"},{"comment":"The OCV formula is written with an unclear subscript E_substrate+xE_Li(Na) and appears to be missing the proper reference-state terms. Please rewrite the formula with consistent notation, for example OCV = -(E_{substrate+xA} - E_{substrate} - xE_A)/(xze).","section":"Eq. (5)"},{"comment":"The figure callouts Figures 4(b) and 4(b) should be Figures 5(b) and 5(c).","section":"Section 3.2"},{"comment":"The ring composition is described as 4-, 5-, and 8-membered in the Introduction but as 5-, 6-, and 8-membered in the abstract and conclusions. Please harmonize these descriptions.","section":"Introduction and Sections 3.1/4"},{"comment":"The statement that data can be accessed by contacting the corresponding author is not a reproducible data-sharing plan. Please provide a repository link or a detailed data availability statement.","section":"Data access statement"}],"recommendation":"major_revision","confidential_remarks":"The phonon stability issue is the key risk to the central claim. I would recommend that the editor require a dedicated convergence test for the -0.5 THz mode before publication. The manuscript also contains several presentation errors (swapped extrema, garbled OCV formula, inconsistent AIMD times, and inconsistent ring-member descriptions) that indicate a careful proofreading pass is needed. The paper is otherwise methodologically straightforward and the results are plausible, so a major revision appears appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a standard but competent DFT screening study, the first to evaluate the known 2013 HOP-graphene allotrope as a Li/Na-ion battery anode. The headline numbers – capacities of 1338 mAh/g (Li) and 1227 mAh/g (Na), barriers of 0.70 eV (Li) and 0.39 eV (Na), average OCVs of 0.42 V and 0.33 V – are derived consistently from the computed energies. The identification of pentagonal bridges as fast diffusion channels is a nice, concrete insight. The mechanical stability check via Born–Huang criteria is standard, and the comparisons against other 2D carbon anodes are useful. I do not see a major red flag in the batch of calculations.\n\nThe soft spots are real but not fatal. The most important is the -0.5 THz imaginary phonon mode near Gamma, reported in Fig. 1(b) and dismissed with a -2 THz screening threshold from Wang et al. That threshold is a high-throughput operational rule, not proof that this particular mode is a numerical artifact. The AIMD run (10 ps) is weak supporting evidence because a soft mode can take longer to manifest. The paper would be stronger with a supercell-size or pseudopotential convergence test, or at least a more explicit caveat. The data availability statement is also thin – only \"contact the corresponding author\" – which is becoming less acceptable for computational screening papers. There are minor typos (e.g., the text refers to Na24C40 despite 22 Na atoms; figure callouts reference nonexistent panels), but they do not affect the conclusions. The attempt frequency in the Arrhenius diffusivity is a free parameter, though it only scales D, not the barriers or capacities.\n\nWho is this for? Researchers working on 2D carbon anodes will read it for another data point, and it is a reasonable candidate for a journal like Applied Surface Science. The phonon situation means I would not want it published exactly as is, but the central argument holds up well enough that a serious referee should see it. I would send it to peer review with a request for a convergence test on the imaginary mode and fuller data deposition.","headline":"Competent, internally consistent DFT screening of HOP-graphene as a Li/Na anode, but the dismissal of a -0.5 THz imaginary phonon mode needs a real convergence test before the stability premise is fully trusted.","tokens_in":15685,"tokens_out":1591,"would_cite":false,"duration_ms":17017,"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":"HOP-graphene is a high-capacity anode for Li- and Na-ion batteries, with DFT-predicted capacities of 1338 and 1227 mAh/g.","keywords":["HOP-graphene","2D carbon allotrope","lithium-ion anode","sodium-ion anode","density functional theory","diffusion barrier","theoretical capacity","open-circuit voltage"],"falsifier":"A phonon calculation in larger supercells with different pseudopotentials that still shows the -0.5 THz mode would mean the free-standing monolayer is dynamically unstable; on the experimental side, a synthesized HOP-graphene with measured Li or Na capacity far below the predicted values would contradict the central claim.","tokens_in":1757,"feed_emoji":"🔋","tokens_out":2034,"duration_ms":122085,"temperature":0.7,"pith_summary":"With a porous two-dimensional carbon allotrope that has not yet been synthesized, the paper aims to show that HOP-graphene can outperform graphite as an anode in both lithium- and sodium-ion batteries. Density functional theory predicts theoretical capacities of 1338 mAh/g for Li and 1227 mAh/g for Na, roughly 3.6 and 35 times the graphite values, by adsorbing up to 24 Li or 22 Na atoms per 40-carbon cell. The same calculations report fast ion mobility, with lowest energy barriers of 0.70 eV for Li and 0.39 eV for Na, and average open-circuit voltages of 0.42 V and 0.33 V, within the safe anodic range. The sheet is described as metallic and mechanically stable both before and after full loading, which would let the electrode keep its conductivity during charge and discharge. If these first-principles predictions hold, HOP-graphene would be a concrete candidate for next-generation lithium- and sodium-ion anodes.","feed_headline":"3.6x lithium, 35x sodium vs graphite in one carbon sheet","feed_subtitle":"First-principles calculations predict a 5-6-8 ring carbon sheet can hold 24 Li or 22 Na atoms per cell.","key_machinery":"The carrying object is the HOP-graphene lattice: a rectangular 10-carbon unit cell that tiles into octagons, hexagons, and pentagons. The octagonal hollow is the deepest adsorption site for both Li and Na, the pentagonal bridges connecting adjacent octagons form the lowest-energy migration paths, and a large vacuum gap isolates a true monolayer in the calculations. The argument itself is carried by a density-functional-theory workflow, including a generalized-gradient exchange-correlation functional, projector-augmented-wave potentials, a dispersion correction, nudged-elastic-band transition searches, and Bader charge decomposition, which converts the ring geometry into adsorption energies, diffusion barriers, capacities, open-circuit voltages, and charge transfer. The key structural statement is that the pentagonal motif, not the hexagon, is what gives fast diffusion, and the octagonal pore is what gives high storage.","core_discovery":"The paper claims that HOP-graphene, a rectangular two-dimensional carbon allotrope made of five-, six-, and eight-membered rings, is a high-performance anode material for both lithium- and sodium-ion batteries. On the basis of density functional theory, it reports that the monolayer binds single Li and Na atoms most strongly at the octagonal hollow sites, with adsorption energies of -2.92 eV and -2.32 eV, and that full double-sided loading reaches 24 Li or 22 Na atoms per 40-carbon cell. From these loadings it derives theoretical capacities of 1338 mAh/g (Li) and 1227 mAh/g (Na), against 372 and 35 mAh/g for graphite. Diffusion calculations give minimum barriers of 0.70 eV for Li and 0.39 eV for Na, both along a path through pentagonal bridges, and Bader analysis shows +0.89 e and +0.87 e transferred from the metal atoms to the carbon sheet. Average open-circuit voltages of 0.42 V and 0.33 V place the electrode inside the safe operating window, and the metallic character of the sheet is reported to survive full loading.","pith_inferences":["The same 5-6-8 ring geometry could be screened for potassium, magnesium, calcium, or aluminum ions; the pore size and charge-transfer pattern reported here suggest those chemistries are a natural next test, but the paper does not run it.","The -0.5 THz mode is the cheapest place to check the prediction: a phonon calculation in larger supercells with different pseudopotentials would either confirm the artifact or invalidate the stability premise.","The present numbers describe a pristine, defect-free monolayer; a real electrode will contain grain boundaries, vacancies, and contact layers, so practical capacities are likely lower until those effects are included."],"forward_implications":["If the central claim is correct, the Li capacity of 1338 mAh/g is 3.6 times the 372 mAh/g of graphite anodes, so a HOP-graphene anode could cut anode mass roughly by that factor at full lithiation.","The Na capacity of 1227 mAh/g is 35 times the 35 mAh/g of graphite, making sodium-ion chemistry far more competitive than it is on carbon.","A room-temperature Na diffusion coefficient of 2.78e-6 cm2/s, with a 0.39 eV barrier, is orders of magnitude above Na in graphite, so fast charging would not be bottlenecked by sodium mobility.","Average open-circuit voltages near 0.42 V for Li and 0.33 V for Na sit below the 1 V threshold for practical anodes, balancing cell voltage against the risk of metal plating.","Because the monolayer is reported to remain metallic at full loading, the electrode would not lose its electrical conductivity as the battery charges."],"supporting_citations":[{"why":"Defines the HOP-graphene structure used throughout as the substrate.","marker":"[41]"},{"why":"Supplies the generalized-gradient exchange-correlation functional used for all total-energy calculations.","marker":"[47]"},{"why":"Provides the projector-augmented-wave treatment of electron-ion interactions in the DFT setup.","marker":"[48]"},{"why":"Adds the dispersion correction that controls the Li/Na adsorption energetics.","marker":"[51]"},{"why":"Justifies treating the -0.5 THz imaginary frequency as within a stable range; it is the cited authority for the stability premise.","marker":"[55]"},{"why":"Sets the lithium reference capacity of graphite (372 mAh/g) that the paper's capacity must beat.","marker":"[32]"},{"why":"Sets the sodium reference capacity of graphite (35 mAh/g) that motivates the sodium-ion part.","marker":"[33]"},{"why":"Provides the nudged-elastic-band method used to compute the diffusion barriers.","marker":"[61]"},{"why":"Gives the graphite Li/Na diffusion coefficients used as baselines for the mobility claims.","marker":"[71]"}],"fun_headline_variants":["3.6x Li, 35x Na capacity vs graphite in a carbon monolayer","DFT predicts carbon allotrope anode with 1338 mAh/g Li capacity","Low barriers: 0.70 eV (Li), 0.39 eV (Na) on HOP-graphene anode","HOP-graphene: 1338 mAh/g Li, 1227 Na from DFT","New 2D carbon holds 24 Li or 22 Na atoms per cell"],"cache_read_input_tokens":17792,"weakest_assumption_plain":"The load-bearing premise is that the small negative phonon frequency near -0.5 THz is a harmless numerical artifact; this is accepted on the authority of a cited -2 THz threshold instead of a convergence test, and if that mode is real the HOP-graphene sheet is unstable and none of the battery findings stand.","fun_headline_variants_meta":{"raw":{"variants":["3.6x Li, 35x Na capacity vs graphite in a carbon monolayer","DFT predicts carbon allotrope anode with 1338 mAh/g Li capacity","Low barriers: 0.70 eV (Li), 0.39 eV (Na) on HOP-graphene anode","HOP-graphene: 1338 mAh/g Li, 1227 Na from DFT","New 2D carbon holds 24 Li or 22 Na atoms per cell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001553,"raw_usage":{"total_tokens":6244,"prompt_tokens":1018,"completion_tokens":5226,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":5105}},"tokens_in":634,"tokens_out":5226,"duration_ms":32324,"temperature":1.0,"reasoning_tokens":5105,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:20:23.507833+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A phonon calculation in larger supercells with different pseudopotentials that still shows the -0.5 THz mode would mean the free-standing monolayer is dynamically unstable; on the experimental side, a synthesized HOP-graphene with measured Li or Na capacity far below the predicted values would contradict the central claim.","supporting_citations":[{"cited_title":"Perdew, Kieron Burke, and Matthias Ernzerhof","cited_arxiv_id":null,"evidence_quote":"Supplies the generalized-gradient exchange-correlation functional used for all total-energy calculations."},{"cited_title":"Theoretical prediction of a new two-dimensional carbon al- lotropeandndrbehaviourofitsone-dimensionalderivatives","cited_arxiv_id":null,"evidence_quote":"Defines the HOP-graphene structure used throughout as the substrate."},{"cited_title":"High- throughputcomputationalscreeningoftwo-dimensionalsemiconduc- tors","cited_arxiv_id":null,"evidence_quote":"Justifies treating the -0.5 THz imaginary frequency as within a stable range; it is the cited authority for the stability premise."},{"cited_title":"Building better batteries.nature, 451(7179):652–657, 2008","cited_arxiv_id":null,"evidence_quote":"Sets the lithium reference capacity of graphite (372 mAh/g) that the paper's capacity must beat."},{"cited_title":"Emerging non- lithium ion batteries.Energy Storage Materials, 4:103–129, 2016","cited_arxiv_id":null,"evidence_quote":"Sets the sodium reference capacity of graphite (35 mAh/g) that motivates the sodium-ion part."},{"cited_title":"Nudgedelastic band method for finding minimum energy paths of transitions","cited_arxiv_id":null,"evidence_quote":"Provides the nudged-elastic-band method used to compute the diffusion barriers."},{"cited_title":"Diffusionofalkalimetalsinthefirststagegraphiteintercalation compounds by vdw-dft calculations","cited_arxiv_id":null,"evidence_quote":"Gives the graphite Li/Na diffusion coefficients used as baselines for the mobility claims."}],"review_version":1}