{"id":"680aaa77-5e6e-442b-bd09-f20e9885eff9","arxiv_id":"2506.05021","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"DFT calculations trace proton pathways and energy barriers across g-C3N4/TiO2-B(001) interfaces, reporting a 1.103 eV (pristine) and 0.999 eV (Li-F doped) surface-to-interface barrier that is called rate-determining yet not rate-inhibiting.","lead":"This computational study uses density functional theory to simulate water splitting and proton migration at the interface between g-C3N4 and TiO2-B(001), with and without lithium-fluorine doping. It reports migration barriers around 1 eV and claims the materials are efficient hydrogen evolution photocatalysts, a conclusion the paper's own numbers undermine.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Paper's own ~1 eV proton migration barrier implies negligible HER rate at 300 K, contradicting the 'efficient photocatalyst' claim and the statement that the barrier is 'not sufficiently high to inhibit'.","rationale":"The reader's verdict rejects the paper on multiple grounds, listing van der Waals corrections as the weakest assumption but also noting the internal contradiction between the ~1 eV barrier and the efficiency claim, and the non-optimal ΔG(*H). I focus on the internal numerical contradictions because they are more directly load-bearing than the functional choice: they are the paper's own computed values contradicting its headline performance claim. A 1.103 eV rate-determining barrier gives an essentially zero thermal rate at room temperature, and the paper's own ΔG(*H) values are far from the Sabatier optimum, with the text even conceding that stronger adsorption 'could impede H2 desorption'. These inconsistencies do not rely on a specific functional and would survive most dispersion corrections. The van der Waals issue is real and would affect quantitative barriers and geometries, but it is secondary; even a 0.2–0.3 eV shift from dispersion corrections would not transform a ~1 eV barrier into a readily accessible room-temperature step, nor would it move ΔG from −0.6 eV to near zero. Therefore the verdict should remain REJECT, and the reader's REJECT verdict is unchanged. Agreement is 'partial' because the reader's explicitly identified weakest assumption (PBE without vdW) is not the concern I consider most load-bearing, though the reader's rationale does include the internal inconsistencies.","tokens_in":15067,"tokens_out":3820,"duration_ms":48324,"concrete_test":"Estimate the proton-transfer rate constant at 300 K using the Eyring equation k = (kBT/h)exp(−ΔG‡/kBT) with the reported 1.103 eV barrier, including zero-point and tunneling corrections if vibrational frequencies are available. Compare the resulting turnover frequency per active site to the range typically required for practical photocatalytic HER (e.g., >10⁻³ s⁻¹ per site). If the computed rate is below ~10⁻¹⁰ s⁻¹, then the statement that a ~1 eV barrier is 'not sufficiently high to inhibit' HER is quantitatively false, and the central claim of an efficient HER photocatalyst fails regardless of functional or dispersion corrections.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper identifies proton migration from TiO2-B(001) to the interface as the rate-determining step, with barriers of 1.103 eV (pristine) and 0.999 eV (Li-F doped), yet concludes the material has 'remarkable potential as efficient HER photocatalyst'. This is internally inconsistent. At 300 K, Eyring transition-state theory gives k ≈ (kBT/h)·exp(−Ea/kBT) ≈ 1.3×10⁻¹⁹ s⁻¹ for Ea = 1.1 eV; even with a generous prefactor of 10¹³ s⁻¹, the rate is ~10⁻⁶ s⁻¹ per site — one proton transfer per site every ~12 days. The authors' assertion that 'the diffusion energy barrier of approximately 1 eV is not sufficiently high to inhibit the HER process' is quantitatively unsupported and would require extreme heating, electrochemical overpotential, or tunneling enhancements not discussed. Similarly, the computed ΔG(*H) values of −0.586 eV and −0.646 eV are far from the Sabatier optimum near 0 eV, and the paper itself admits the more negative value for Li-F doping 'could impede H2 desorption', directly contradicting the abstract's 'optimal proton adsorption and desorption characteristics'. These are internal numerical contradictions that do not depend on the choice of exchange-correlation functional; they undermine the central efficiency claim even if the van der Waals issue were resolved.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports DFT-PBE calculations of water adsorption/dissociation, proton adsorption and migration barriers, and HER free energies for g-C3N4/TiO2-B(001) heterostructures with and without Li-F co-doping. It concludes that both heterostructures are efficient HER photocatalysts, identifies proton migration from the TiO2-B(001) surface to the interface as the rate-determining step (barriers 1.103 eV and 0.999 eV), and reports hydrogen adsorption free energies of -0.586 eV and -0.646 eV.","tokens_in":15312,"tokens_out":3873,"duration_ms":46090,"significance":"If the quantitative results were reliable, the systematic CI-NEB study of proton pathways, the Bader-charge/PDOS analysis, and the comparison with literature water-splitting energies would provide a useful mechanistic reference for g-C3N4/TiO2-B interfaces. The manuscript also includes machine-checkable reproducible data availability statements and verifies transition states by vibrational analysis. However, the central efficiency claim is contradicted by the paper's own computed barriers and adsorption free energies, and the absence of dispersion corrections undermines the quantitative interface energetics. As it stands, the contribution is not reliable enough to support the stated conclusions.","major_comments":[{"comment":"The claim that an ~1 eV barrier (1.103 eV pristine, 0.999 eV doped) \"is not sufficiently high to inhibit the HER process\" is quantitatively unsupported. A transition-state-theory estimate with the kBT/h prefactor gives ~10^-19 s^-1 at 300 K; even with a generous prefactor of 10^13 s^-1, the rate is ~10^-6 s^-1 per site, meaning one proton transfer every ~12 days. The manuscript provides no rate model, no overpotential or electrochemical driving force, and no tunneling analysis. This is an internal inconsistency independent of the choice of exchange-correlation functional and directly undermines the paper's efficiency conclusion.","section":"Proton migration from TiO2-B(001) to g-C3N4 surface (Fig. 4)"},{"comment":"The abstract claims the heterojunction surface exhibits \"optimal proton adsorption and desorption characteristics,\" but the computed ΔG(*H) values are -0.586 eV and -0.646 eV, far from the Sabatier optimum near 0 eV. The text itself admits that Li-F doping \"could impede H2 desorption\" and that \"excessively strong adsorption can inhibit H2 production.\" These statements directly contradict the abstract's central claim; the authors should quantify the effect of these free energies on HER rate through a microkinetic or exchange-current-density analysis rather than relying on qualitative 'optimal' language.","section":"HER process on different g-C3N4 surfaces; Abstract"},{"comment":"The methods section specifies the PBE functional with no van der Waals correction for a layered g-C3N4/TiO2-B(001) heterostructure. Dispersion interactions dominate the interlayer binding and affect the interface geometry, adsorption energies, and NEB barriers along the proton migration path. Since every calculated barrier and adsorption energy feeds directly into the mechanistic claims, the authors should provide at least a benchmark of the key IS/TS/FS energies with a dispersion-corrected method (e.g., DFT-D3, dDsC, or optB88-vdW). Without this, the quantitative barriers, especially the rate-determining 1 eV step, are not reliable.","section":"Calculation methods"}],"minor_comments":[{"comment":"The text states that \"the adsorption energy of the adsorbed water system is 0.442 eV lower than that of the dissociated system\" and then later refers to a \"dissociation adsorption energy of 0.177 eV, with an increase of 0.441 eV after dissociation.\" These statements are ambiguous and appear inconsistent; please report the reaction energy explicitly as endothermic by a single number and define all quantities used.","section":"Water-splitting on different TiO2-B(001) surfaces"},{"comment":"The Ti13 site is introduced as the most stable adsorption site, but the manuscript does not define which titanium atom is labeled Ti13 in the structural figures; please indicate it clearly in the figure or caption.","section":"Fig. 1(b)"},{"comment":"The manuscript says all transition states were verified by imaginary frequency analysis, but no imaginary frequencies or vibrational analysis details are reported anywhere; please provide the values or cite the verification explicitly.","section":"Calculation methods"},{"comment":"The manuscript contains numerous typographical errors, including \"emloyed,\" \"ehxibits,\" \"adsoprtion,\" \"favorbale,\" \"hte,\" \"remarkbale,\" \"stesps,\" \"chanllenge,\" \"efficiecy,\" \"stratey,\" and \"the the.\" A thorough language edit is needed.","section":"Throughout"},{"comment":"The author list includes Hanyang Jia, but the contributions section lists \"Hanyang Ji\"; please correct this inconsistency.","section":"Author contributions"}],"recommendation":"reject","confidential_remarks":"The paper's own computed numbers contradict its central claim of \"efficient HER photocatalyst\": a 1 eV rate-determining barrier and ΔG(*H) values near -0.6 eV are both far from the values expected for high HER activity. The missing van der Waals correction is an additional load-bearing issue for a layered heterostructure. Even if the authors were to repeat the calculations with a dispersion-corrected functional, the internal numerical contradictions would still require a fundamental change in the conclusions. The manuscript is therefore not suitable for publication in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper adds real mechanistic numbers to a system the authors already published, but the headline claim dies on the authors' own numbers. The ~1 eV proton-migration barrier is the rate-determining step they identify, and at 300 K that barrier gives a per-site rate of roughly 10^-6 s^-1 or slower. Calling it 'not sufficiently high' is not supported by any kinetic estimate in the paper. The paper would need tunneling, hot-carrier effects, or an overpotential argument to rescue it, and it provides none.\n\nWhat is new: the water dissociation energies for the heterostructures (0.441 and 0.338 eV), the CI-NEB pathways that separate surface-to-interface (1.103/0.999 eV) from interface-to-g-C3N4 (0.168/0.226 eV) transport, and the ΔG(*H) values (-0.586 and -0.646 eV). These are concrete numbers absent from the group's prior work, and the charge-density/PDOS analysis of the transition states is competently done.\n\nSoft spots, in order: first, the internal contradiction. The abstract claims 'optimal proton adsorption and desorption characteristics,' but in the HER section the authors admit that the more negative ΔG(*H) 'could impede H2 desorption.' Both cannot be true, and the computed values are far from the Sabatier optimum near zero. Second, the 1 eV barrier is treated as an acceptable hurdle without a single rate estimate. Eyring with a generous prefactor gives k ~ 10^-6 s^-1 per site at 300 K, i.e., one proton hop per site every several days; the paper's 'heating and stirring' remark does not address the photochemical context. Third, the lack of any van der Waals correction for a stacked layered interface is a real methodological weakness. I consider it secondary, because the rate problem would survive a re-relaxation, but it should be fixed if the numbers are to be trusted.\n\nMinor items: zero-point energies are not included for the barriers, which matters for proton transfer; the adsorption-energy definition for proton includes the H2 binding energy, which is a bit unusual; and the polarization direction is taken from the authors' ref. 13, which is fine as background.\n\nWho gets value: computational researchers working on g-C3N4/TiO2 interfaces could use the pathway analysis as a starting template, but not the conclusion that Li-F doping yields an efficient HER photocatalyst.\n\nI would send this to review rather than desk-reject, because the mechanistic data deserve scrutiny and the contradiction is exactly what a referee should catch. But I would not accept it in current form; the efficiency claim needs either quantitative kinetic support or removal.","headline":"Routine follow-up DFT study whose central HER efficiency claim is contradicted by its own ~1 eV proton-migration barrier and non-optimal adsorption free energies.","tokens_in":15880,"tokens_out":3293,"would_cite":false,"duration_ms":33889,"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 claims that proton migration from the TiO2-B(001) surface to the interface is the rate-determining step for hydrogen evolution, with a 1.103 eV barrier that Li–F co-doping lowers to 0.999 eV.","keywords":["photocatalytic water splitting","hydrogen evolution reaction","g-C3N4/TiO2-B heterostructure","Li-F co-doping","proton migration barrier","CI-NEB","density functional theory","Gibbs free energy of hydrogen adsorption"],"falsifier":"Recompute the same CI-NEB pathways with a dispersion-corrected functional, for example DFT-D3 or optB88-vdW, and check whether the 1.103 eV and 0.999 eV surface-to-interface barriers, and the minimum-energy path itself, survive; a shift of more than about 0.2 eV would overturn the rate-determining-step assignment.","tokens_in":14854,"feed_emoji":"⚫","tokens_out":8328,"duration_ms":85740,"temperature":0.7,"pith_summary":"Using density functional theory, this paper traces the full hydrogen evolution path on g-C3N4/TiO2-B(001) heterostructures: water adsorption and splitting on the titania surface, proton migration across the interface, and H2 formation on the g-C3N4 surface. The central claim is that the step controlling the overall rate is proton migration from the TiO2-B(001) surface to the interface, which must break and reform a polar O–H covalent bond and costs 1.103 eV in the pristine heterostructure. Li–F co-doping weakens that bond and strengthens interfacial polarization, lowering the barrier to 0.999 eV, while interfacial O–H···N hydrogen bonds make the later hop from the interface to g-C3N4 much easier (0.168 eV pristine, 0.226 eV doped). The paper concludes that both heterostructures are promising HER photocatalysts, with the doped system also requiring an interfacial proton-diffusion step of 0.701 eV that it argues does not determine the overall rate.","feed_headline":"Li-F doping cuts the 1.1 eV proton barrier limiting H2 production","feed_subtitle":"DFT maps the full water-splitting path on g-C3N4/TiO2-B; the bottleneck drops from 1.103 to 0.999 eV.","key_machinery":"The load-bearing machinery is the proton-transfer coordinate O29→O12→N16 (pristine) and O29→O9→O12→N16 (doped), resolved with climbing-image nudged elastic band calculations. The barrier heights are controlled by two competing interactions: a polar O–H covalent bond whose cleavage and recombination costs about 1 eV, and an O–H···N hydrogen bond to pyridine-type nitrogen that stabilizes the interface and lowers the later hop to 0.168 or 0.226 eV. Bader charge analysis, charge-density differences, and projected density of states are used to show that the proton loses about 0.6|e| to oxygen and that hybridization between the proton and its neighbors weakens at the transition state.","core_discovery":"The authors establish that the rate-determining elementary step in the g-C3N4/TiO2-B(001) heterostructure HER cycle is not water splitting (0.442 eV, endothermic with no transition state) nor the final proton hop to g-C3N4 (0.168 eV), but the migration of a proton from the TiO2-B(001) surface to the interface, where oxygen acts as both proton donor and acceptor through a polar O–H covalent bond; this barrier is 1.103 eV. Li–F co-doping weakens the O–H bond and adds interfacial polarization, reducing that barrier to 0.999 eV, and introduces an interfacial diffusion step (0.701 eV) that the authors argue does not control overall HER. The thermodynamic driving force is a monotonic decrease in proton adsorption energy from the TiO2-B surface (0.041 eV) to the interface (−0.380 eV) to g-C3N4 (−0.586 eV; −0.646 eV after doping), and the computed hydrogen adsorption free energies $ΔG(*H)$ are −0.586 eV and −0.646 eV, which the paper reads as a favorable balance of adsorption and desorption.","pith_inferences":["The paper stops at electronic-structure barriers; a microkinetic model translating the 1.103 and 0.999 eV barriers into turnover rates would clarify whether stirring and mild heating are truly enough to drive HER at practical rates.","Because the calculations use the PBE functional without dispersion corrections, the absolute barriers are more uncertain than the relative trend; the comparison between pristine and doped systems is likely more robust than any single barrier value.","The monotonically decreasing adsorption energies suggest that the interface itself could act as a proton reservoir; one testable design extension is to introduce dopants that stabilize the interface state without adding a 0.701 eV diffusion step.","If the rate-limiting step is O–H bond cleavage during the surface-to-interface hop, a normal kinetic isotope effect for H/D substitution should appear in the HER rate; measuring it would provide an experiment-level check of the mechanism."],"forward_implications":["Water splitting on these surfaces is endothermic and has no transition state, so its energy cost (0.442, 0.441, and 0.338 eV) is a thermodynamic hurdle rather than a kinetic one, and Li–F co-doping cuts that cost by about 0.1 eV.","The proton's adsorption energy falls monotonically from the TiO2-B surface to the interface to the g-C3N4 surface, so after water splitting the proton is thermodynamically pulled across the junction.","The rate-limiting step is the surface-to-interface hop (1.103 eV), not the interface-to-g-C3N4 hop (0.168 eV), because the first hop requires breaking a polar O–H bond while the second is assisted by an O–H···N hydrogen bond.","Li–F co-doping lowers the rate-limiting barrier to 0.999 eV by weakening the O–H bond and enhancing interfacial polarization, and it adds a separate interfacial diffusion barrier (0.701 eV) that does not determine overall HER.","The calculated $ΔG(*H)$ values (−0.586 eV pristine, −0.646 eV doped) imply stable proton adsorption; the paper concludes this is a favorable balance while noting that stronger adsorption can hinder H2 desorption."],"supporting_citations":[{"why":"Supplies the parent g-C3N4/TiO2-B heterostructure design and its predicted HER activity, which this paper extends.","marker":"12"},{"why":"Established that Li–F co-doping enhances photocatalytic hydrogen evolution in this heterostructure, the baseline for the doped system studied here.","marker":"13"},{"why":"The VASP plane-wave DFT code used for all total-energy and electronic-structure calculations.","marker":"14"},{"why":"The PBE exchange-correlation functional defining the level of theory for every adsorption energy and barrier.","marker":"16"},{"why":"The climbing-image nudged elastic band method used to locate transition states and compute all proton migration barriers.","marker":"20"},{"why":"Provides the literature value (0.42 eV) used to benchmark the computed water-splitting energy on TiO2-B(001).","marker":"22"},{"why":"Underpins the Sabatier-principle interpretation of $ΔG(*H)$ as the HER activity descriptor.","marker":"30"}],"fun_headline_variants":["Li-F doping lowers the 1.1 eV proton barrier that controls HER","Proton migration to interface, not water split, is the HER rate limit","DFT: Li-F co-doping reduces the proton transfer barrier from 1.10 to 0.999 eV","Polar O–H bond sets 1.1 eV proton hurdle; Li-F doping weakens it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole barrier picture rests on the assumption that the computer model describes the weak forces between the stacked layers accurately enough, even though the calculation uses a standard approximation that often underestimates such forces.","fun_headline_variants_meta":{"raw":{"variants":["Li-F doping lowers the 1.1 eV proton barrier that controls HER","Proton migration to interface, not water split, is the HER rate limit","DFT: Li-F co-doping reduces the proton transfer barrier from 1.10 to 0.999 eV","Polar O–H bond sets 1.1 eV proton hurdle; Li-F doping weakens it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001244,"raw_usage":{"total_tokens":5196,"prompt_tokens":1130,"completion_tokens":4066,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":746,"completion_tokens_details":{"reasoning_tokens":3968}},"tokens_in":746,"tokens_out":4066,"duration_ms":31152,"temperature":1.0,"reasoning_tokens":3968,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:27:55.868881+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same CI-NEB pathways with a dispersion-corrected functional, for example DFT-D3 or optB88-vdW, and check whether the 1.103 eV and 0.999 eV surface-to-interface barriers, and the minimum-energy path itself, survive; a shift of more than about 0.2 eV would overturn the rate-determining-step assignment.","supporting_citations":[],"review_version":1}