{"id":"9ca99dd7-a062-4b21-836b-0d5440180781","arxiv_id":"2506.00604","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Potassium-decorated graphenyldiene is predicted by DFT to store 8.82 wt% hydrogen with reversible adsorption energies of -0.11 to -0.14 eV per H2, desorbing at 300 K in AIMD.","lead":"This study predicts that adding potassium to a new carbon sheet called graphenyldiene lets it store hydrogen, reaching 8.82% hydrogen by weight, above the U.S. Department of Energy target. The authors also show through simulations that the hydrogen can detach at room temperature, which could make the material reusable for clean energy storage.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central reversibility claim rests on uncorrected PBE-D2 adsorption energies of only -0.11 to -0.14 eV; a standard zero-point energy correction could push them out of the reversible window, so the paper needs that correction before the DOE-target and AIMD conclusions are quantitative.","rationale":"I agree with the reader's identification of ZPE as the weakest load-bearing assumption. The paper's own internal evidence (weak adsorption, low desorption temperatures, AIMD desorption) is consistent but does not remove the need for zero-point corrections; without them, the central quantitative claim lacks support. The other issues (Eq. (5) units, DOE target level, Kubas interpretation) are important but would be corrected by rewriting and do not threaten the core physics as directly. The proposed check is computationally inexpensive and would settle whether the adsorption energies remain in the reversible window after the standard correction.","tokens_in":14339,"tokens_out":4147,"duration_ms":42013,"concrete_test":"Compute harmonic vibrational frequencies for isolated H2 and for the saturated K@GPD–18H2 structure (or at least K@GPD–2H2 and K@GPD–18H2) using the same PBE-D2 setup, and form E_ad^ZPE = E_ad + (ZPE_adsorbed - ZPE_gas) per H2. If any coverage's ZPE-corrected value is above -0.1 eV or positive, the reversibility claim fails; if all remain below -0.1 eV, the omission is not fatal. Running the same correction with a DFT-D3(BJ) or vdW-DF functional would additionally test whether the -0.11 to -0.14 eV window is an artifact of D2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"K@GPD's central claim—reversible near-ambient H2 storage exceeding the DOE target—is anchored by Table 2's E_ads = -0.11 to -0.14 eV per H2, which the paper places in the -0.5 to -0.1 eV 'reversible window'. These values are static PBE-D2 electronic energies with no zero-point energy correction. For molecular physisorption, the difference in zero-point energy between adsorbed and gas-phase H2 is typically tens of meV; at these shallow binding energies a +50 meV correction would push most entries above -0.1 eV, and a +100 meV correction would make them positive, removing the thermodynamic basis for both 'reversible' and the desorption temperatures in Table 2. The omission is not documented or justified in Section 2. The desorption temperatures themselves are also not trustworthy as printed: Eq. (5) is dimensionally inconsistent (R and kB both appear while ΔS is quoted as the liquid-vapor entropy change, 75.44 J mol^-1 K^-1), so the 163–185 K column is not an independent validation. A secondary but real issue is that 8.82 wt% is a material-level ratio, while the 5.5 wt% DOE target is system-level; claiming the material 'exceeds' the target conflates the two. The AIMD desorption at 300 K under effectively zero H2 pressure is consistent with weak binding but does not by itself establish reversible uptake at ambient pressure. The ZPE omission is the load-bearing concern because it directly undermines the quantitative 'optimal adsorption energy' premise from which reversibility and capacity are inferred.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles DFT calculations of potassium-decorated graphenyldiene (K@GPD) as a candidate for reversible hydrogen storage. The authors find that K binds to the GPD monolayer with energies between -2.58 and -2.62 eV, that H2 adsorption energies are -0.11 to -0.14 eV per molecule, that the material-level gravimetric capacity reaches 8.82 wt% at 18 H2 molecules, and that AIMD simulations at 300 K show structural stability of K@GPD and H2 desorption from the saturated system. A grand canonical analysis yields a usable capacity of 5.88 wt% between 25 °C/30 atm and 100 °C/3 atm. The conclusions claim that K@GPD is a reversible, near-ambient hydrogen storage material exceeding the U.S. DOE gravimetric target.","tokens_in":14643,"tokens_out":4492,"duration_ms":41860,"significance":"If the central claims survive the thermochemical corrections discussed below, K@GPD would be a reasonable addition to the alkali-metal-decorated carbon family, with a material-level capacity above the commonly cited DOE gravimetric target. The paper has several concrete strengths: the adsorption energies are direct DFT total-energy differences (Eq. 3) rather than fitted parameters; K binding is benchmarked against the bulk K cohesive energy; GPD is checked for dynamical stability via phonons; the H2 coverage is varied systematically from 2 to 18 molecules; and the grand canonical analysis in Fig. 10 provides a thermodynamic estimate of usable capacity. The main weaknesses are the omission of zero-point energy corrections for H2, the dimensional inconsistency in the desorption-temperature formula, and the conflation of material-level capacity with the system-level DOE target.","major_comments":[{"comment":"The central reversibility claim rests on static PBE-D2 adsorption energies of -0.11 to -0.14 eV per H2, which are placed in the '-0.5 to -0.1 eV reversible window' without any zero-point energy (ZPE) correction. For molecular physisorption, the ZPE difference between adsorbed and gas-phase H2 is typically on the order of tens of meV; a +50 meV correction would shift most values in Table 2 to -0.06 to -0.09 eV, outside the cited reversible window, and a +100 meV correction would make them positive. The manuscript neither computes ZPE corrections nor justifies neglecting them. Please provide ZPE-corrected adsorption enthalpies (or an explicit quantitative justification) and re-evaluate the reversibility claim and the desorption temperatures derived from them.","section":"§2, Table 2, Eq. (3)"},{"comment":"Equation (5), as printed, is dimensionally inconsistent: it places both R and k_B in the denominator together with ΔS, so the expression |E_ads|/(R k_B ΔS) does not have units of temperature. The numerical values in Table 2 (e.g., 0.13 eV giving roughly 163 K) suggest the intended formula is T_des = |E_ads|/(k_B (ΔS/R)) or, equivalently, T_des = |E_ads|/(R ΔS) with ΔS expressed per molecule. Please correct the equation, define all symbols precisely, and justify the use of the liquid-vapor entropy change (75.44 J mol^-1 K^-1) for a surface desorption process rather than the gas-phase translational entropy at the relevant pressure.","section":"§2, Eq. (5)"},{"comment":"The abstract and Section 3 state that the 8.82 wt% hydrogen storage capacity 'exceeds the U.S. DOE target' of 5.5 wt%. However, the capacity in Eq. (4) is a material-level ratio (denominator includes only the masses of C, K, and H in the monolayer), whereas the cited DOE target is for an onboard hydrogen storage system that includes the tank, balance of plant, and other system components. The grand canonical usable capacity of 5.88 wt% reported in Section 3 is likewise a material-level value. Please rephrase the claims to distinguish explicitly between material-level capacity and the system-level DOE target.","section":"Abstract and §3, Eq. (4)"},{"comment":"The AIMD simulations at 300 K show H2 desorption from K@GPD-18H2 under effectively zero applied H2 pressure (the simulation cell contains only the monolayer and the adsorbed molecules). Desorption in vacuum is consistent with weak binding, but it does not by itself establish reversible uptake and release under technologically relevant pressures. The grand canonical analysis in Fig. 10 is the appropriate route to quantify reversibility, and the authors should connect the AIMD desorption observation to that analysis rather than presenting the AIMD result alone as 'demonstrating excellent reversibility'.","section":"§3, Fig. 9, AIMD discussion"}],"minor_comments":[{"comment":"Equation (1) contains a typo: the last term should be the charge density of the pristine GPD monolayer, ρ(GPD), not ρ(K@GPD) as written, otherwise the left-hand side is identically zero.","section":"§2, Eq. (1)"},{"comment":"The text states that desorption temperatures range from 163 to 285 K, but Table 2 lists a maximum of 184.62 K; the upper bound appears to be a typo and should read approximately 185 K.","section":"§3, text near Table 2"},{"comment":"The sentence 'For instance, K@PHE-graphene, despite storing 7.47 wt%, K@PHE-graphene requires desorption temperatures as high as 423 K' repeats the subject 'K@PHE-graphene'; this should be rephrased for readability.","section":"§3, comparison paragraph"},{"comment":"The caption 'H2 saturation on K@GPD pathway' would be clearer as 'H2 coverage progression on K@GPD' or similar, since the figure shows sequential addition of H2 molecules rather than a reaction pathway.","section":"§3, Fig. 6 caption"},{"comment":"The data access statement says that data can be accessed by contacting the corresponding author; for a computational study, depositing input files (e.g., VASP structures and INCAR/POSCAR/KPOINTS) in a public repository would improve reproducibility and is recommended.","section":"Data access statement"}],"recommendation":"major_revision","confidential_remarks":"This is a competent DFT screening study whose main quantitative claims depend on the zero-point energy correction that is currently missing. The manuscript is likely to be acceptable after the authors add ZPE-corrected adsorption energies (or justify their omission), correct Eq. (5), and reframe the DOE-target comparison as material-level rather than system-level. The AIMD desorption claim should also be softened or tied explicitly to the grand canonical isotherm data. I would not recommend rejection because the underlying framework is standard and the core calculations appear internally consistent; the requested revisions are within the scope of a normal revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things before reading this one. First, the substrate is genuinely new — graphenyldiene (GPD) was proposed by the same group last year, and this is the first K-decoration study on it. Second, the central claim of “reversible storage exceeding the DOE target” is not backed by the numbers as presented. The raw DFT is internally consistent, but the paper is built on static PBE-D2 adsorption energies of only -0.11 to -0.14 eV per H2, and the zero-point energy correction is missing. For molecular physisorption that correction is typically tens of meV, so it could push the binding out of the -0.5 to -0.1 eV reversible window the authors invoke. That is a load-bearing issue, not a nitpick.\n\nWhat the paper does well: the K binding energies are compared against the K cohesive energy, which is the right clustering check; the phonon stability of GPD is shown; the AIMD at 300 K confirms K stays on the surface and that H2 desorbs; and the H2 loading study is systematic, with Bader charge analysis and CDDs. There is no fabrication or circular fitting here — the adsorption energies come straight from total-energy differences. The authors also honestly cite the prior K-decorated carbon literature, and the GPD self-citation is appropriate since the substrate is theirs.\n\nThe soft spots, in order of severity. (1) The ZPE omission is the main one; it directly affects whether -0.11 to -0.14 eV really sits in the reversible window. (2) Eq. (5) is dimensionally nonsense as printed — R and kB both appear while ΔS is a molar entropy — and the desorption temperatures in Table 2, while approximately consistent with |Eads|/ΔS, are not a reliable validation. (3) The 8.82 wt% is a material-level number compared to a system-level DOE target; the thermodynamic analysis in Fig. 10 gives a usable capacity of 5.88 wt%, which is a more honest comparison but still material-level. (4) Calling the interaction “Kubas-type” is a stretch for potassium, an alkali metal with no d states; the H–H elongation and charge transfer suggest polarization, not Kubas. None of these is fatal — they are addressable with corrections — but together they mean the paper currently overstates its quantitative conclusions.\n\nWho this is for: computational materials scientists working on metal-decorated 2D carbon for hydrogen storage. It is a competent incremental data point, not a breakthrough. A serious referee should see it because the substrate is new and the workflow is standard enough that the issues can be fixed. My advice: send it to review, but condition acceptance on adding ZPE corrections, fixing Eq. (5), and softening the DOE target language.","headline":"A competent DFT screening of a new carbon monolayer whose reversibility claims rest on uncorrected adsorption energies and a flawed desorption formula; worth a conditional review, not acceptance as-is.","tokens_in":15265,"tokens_out":2398,"would_cite":false,"duration_ms":25722,"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":"Potassium-decorated graphenyldiene binds hydrogen reversibly and reaches 8.82 wt% storage capacity, exceeding the U.S. DOE target.","keywords":["hydrogen storage","2D materials","density functional theory","graphenyldiene","potassium decoration","physisorption","Kubas interaction","reversible storage"],"falsifier":"Measure the isosteric heat of adsorption of H2 on K@GPD (or compute it with a more accurate dispersion-corrected functional that includes zero-point motion). If the measured binding energy per H2 falls outside roughly 0.1–0.5 eV, or if a longer AIMD run at 300 K does not show H2 desorption while the K adatoms stay bound, the central reversibility claim is falsified.","tokens_in":14101,"feed_emoji":"💧","tokens_out":5338,"duration_ms":47172,"temperature":0.7,"pith_summary":"This paper predicts that decorating the porous carbon monolayer graphenyldiene with potassium atoms turns it into a practical hydrogen storage material. At full coverage the K@GPD system stores 8.82 wt% hydrogen, above the U.S. DOE 5.5 wt% target for onboard storage. Calculated H2 adsorption energies of -0.11 to -0.14 eV fall in the window for reversible adsorption and desorption, and molecular dynamics at 300 K shows hydrogen desorbing while the potassium stays put. If the predictions hold, K@GPD could offer a lightweight, reusable 2D platform for hydrogen storage under near-ambient conditions.","feed_headline":"8.82 wt% hydrogen storage predicted in K-decorated carbon","feed_subtitle":"DFT plus molecular dynamics show the material binds and releases H2 near ambient conditions, beating the DOE target.","key_machinery":"The central device is the K@GPD composite: four potassium adatoms per unit cell at the pore sites of the porous graphenyldiene monolayer. The potassium transfers charge into the carbon pi system, creating active sites where H2 binds through weak Kubas-type interactions (donation from K into the H2 sigma* orbital) without dissociating. The adsorption energy per H2 ($-0.11$ to $-0.14$ eV) sits inside the canonical reversible-storage window, and the van't Hoff desorption temperatures ($163$–$185$ K) plus grand-canonical thermodynamic weighting translate that binding strength into predicted gravimetric capacities.","core_discovery":"The authors show that potassium atoms bind strongly to the GPD monolayer (-2.62 eV at the favored pore site, with no tendency to cluster) and donate about 0.5 e each to the carbon network, making the system metallic. Hydrogen molecules then adsorb with energies between -0.11 and -0.14 eV per H2, in the range typically associated with reversible storage, with H–H bonds lengthening from 0.75 Å to as much as 0.86 Å, which they attribute to Kubas-type interactions. At 18 H2 per cell the gravimetric capacity reaches 8.82 wt% (exceeding the DOE target), and thermodynamic analysis under practical uptake and release conditions (25 °C/30 atm versus 100 °C/3 atm) yields a usable capacity of 5.88 wt%. AIMD at 300 K shows H2 desorption while the K@GPD framework remains intact, indicating reusability.","pith_inferences":["The gap between 8.82 wt% saturation capacity and 5.88 wt% usable capacity implies that roughly a third of stored hydrogen remains bound at the release conditions, so a pressure-swing rather than purely thermal-swing cycle may be needed in practice.","Because the H2 binding is attributed to Kubas-type interactions with only about 0.08 e transferred per H2, the same decoration is likely tunable by substituting other alkali metals or applying strain to adjust the reversible window.","The low calculated desorption temperatures (163–185 K) suggest that true ambient-temperature operation would require pressurization for uptake; a testable extension is to measure isotherms to see if the material works at room temperature under moderate pressure.","The thermodynamic analysis assumes ideal-gas H2 and a fixed entropy change; real porous materials often show pore-confinement effects that shift van't Hoff estimates, so experimental adsorption isotherms on synthesized K@GPD would be the decisive check."],"forward_implications":["If the DFT predictions translate to experiment, K@GPD offers material-level capacity above 8 wt%, exceeding the DOE 5.5 wt% target in a lightweight carbon host with abundant, inexpensive potassium.","Reversible desorption near ambient temperature means spent material could be regenerated without high-temperature energy input, a key practical advantage for mobile storage.","The moderate adsorption energies and low desorption temperatures (163–185 K) indicate a physisorption-type material that releases H2 under mild conditions.","The retention of potassium during adsorption/desorption cycles, shown in AIMD, predicts a reusable substrate without metal clustering or degradation.","The computational workflow (PBE-D2, AIMD, grand-canonical analysis) provides a template for screening other alkali-metal-decorated porous 2D carbons for hydrogen storage."],"supporting_citations":[{"why":"Supplies the GPD monolayer structure, stability, and cohesive properties that K@GPD builds on.","marker":"[47]"},{"why":"The DFT-D2 dispersion correction used for all adsorption energies and H2 binding.","marker":"[51]"},{"why":"The van't Hoff equation used to convert adsorption energies into desorption temperatures.","marker":"[53]"},{"why":"The grand canonical partition function approach used for the uptake/release thermodynamic analysis.","marker":"[55]"},{"why":"K-decorated DHP-graphene: prior AIMD evidence of K retention and comparable capacity supporting the K@GPD result.","marker":"[59]"},{"why":"K-decorated C9N4: comparative reversible storage system with similar non-cooperative H2 adsorption behavior.","marker":"[60]"},{"why":"K@PHE-graphene: defines the reversible adsorption-energy window and provides desorption-temperature comparison.","marker":"[61]"},{"why":"Na@B7N5: supplies the reversible-range criterion and a comparative capacity and desorption benchmark.","marker":"[62]"},{"why":"Kubas-type interaction mechanism used to interpret H–H bond elongation and binding.","marker":"[66]"}],"fun_headline_variants":["K-decorated carbon hits 8.82 wt% hydrogen storage","Reversible H2 storage in K@GPD beats DOE target","DFT design: K-decorated monolayer stores 8.82 wt% H2","Potassium on carbon achieves 8.82 wt% reversible H2","Graphenyldiene + K yields 8.82 wt% hydrogen capacity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions rely on PBE-D2 adsorption energies that omit zero-point energy corrections; if the true H2 binding enthalpy is much weaker (above about -0.1 eV) or different in sign, the reversibility and desorption temperatures lose their quantitative basis.","fun_headline_variants_meta":{"raw":{"variants":["K-decorated carbon hits 8.82 wt% hydrogen storage","Reversible H2 storage in K@GPD beats DOE target","DFT design: K-decorated monolayer stores 8.82 wt% H2","Potassium on carbon achieves 8.82 wt% reversible H2","Graphenyldiene + K yields 8.82 wt% hydrogen capacity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000149,"raw_usage":{"total_tokens":1152,"prompt_tokens":865,"completion_tokens":287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":186}},"tokens_in":481,"tokens_out":287,"duration_ms":3036,"temperature":1.0,"reasoning_tokens":186,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:02:23.217381+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the isosteric heat of adsorption of H2 on K@GPD (or compute it with a more accurate dispersion-corrected functional that includes zero-point motion). If the measured binding energy per H2 falls outside roughly 0.1–0.5 eV, or if a longer AIMD run at 300 K does not show H2 desorption while the K adatoms stay bound, the central reversibility claim is falsified.","supporting_citations":[{"cited_title":"Enhanced hydrogen storage in k and na decorated dhp-graphene monolayer: Dft and gcmc study.Chemical Physics Letters, 865:141932, 2025","cited_arxiv_id":null,"evidence_quote":"K-decorated DHP-graphene: prior AIMD evidence of K retention and comparable capacity supporting the K@GPD result."},{"cited_title":"Sambrano","cited_arxiv_id":null,"evidence_quote":"Supplies the GPD monolayer structure, stability, and cohesive properties that K@GPD builds on."},{"cited_title":"Semiempirical gga-type density functional con- structed with a long-range dispersion correction.Journal of compu- tational chemistry, 27(15):1787–1799, 2006","cited_arxiv_id":null,"evidence_quote":"The DFT-D2 dispersion correction used for all adsorption energies and H2 binding."},{"cited_title":"Reviewofhydrogen storage techniques for on board vehicle applications.International journal of hydrogen energy, 38(34):14595–14617, 2013","cited_arxiv_id":null,"evidence_quote":"The van't Hoff equation used to convert adsorption energies into desorption temperatures."},{"cited_title":"Ultra-high capacity hydrogen storage in a li decorated two- dimensionalc2nlayer","cited_arxiv_id":null,"evidence_quote":"The grand canonical partition function approach used for the uptake/release thermodynamic analysis."},{"cited_title":"Reversible hydrogen storage tendency of light- metal (li/na/k) decorated carbon nitride (c9n4) monolayer.Interna- tional Journal of Hydrogen Energy, 48(67):26301–26313, 2023","cited_arxiv_id":null,"evidence_quote":"K-decorated C9N4: comparative reversible storage system with similar non-cooperative H2 adsorption behavior."},{"cited_title":"Laranjeira, Nicolas F","cited_arxiv_id":null,"evidence_quote":"K@PHE-graphene: defines the reversible adsorption-energy window and provides desorption-temperature comparison."},{"cited_title":"Kaewmaraya, N","cited_arxiv_id":null,"evidence_quote":"Kubas-type interaction mechanism used to interpret H–H bond elongation and binding."}],"review_version":1}