{"id":"92b10f39-0de2-4441-aa14-17e44e209020","arxiv_id":"2606.11064","paper_version":1,"verdict":"UNVERDICTED","confidence":"UNKNOWN","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"The UZH protocol decomposes CP2K simulation errors into basis-set and pseudopotential components via molecular calibration, crystal validation, and external benchmarks, then produces revised MOLOPT basis sets and GTH pseudopotentials.","lead":"The paper presents the UZH protocol, a closed-loop workflow that uses three-way comparisons between CP2K, plane-wave SIRIUS, and all-electron references to separate Gaussian basis set errors from pseudopotential errors in density-functional calculations. A smart generalist might read it to see how computational chemistry tools can be systematically improved for more reliable simulations of molecules and materials.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"SIRIUS PW+GTH and AE FP-LAPW benchmarks assumed converged without shown cutoff/k-point tests","rationale":"The reader's weakest assumption directly identifies the same reference-accuracy premise that must hold for the decomposition to be valid. No other internal inconsistency is visible from the abstract-level description; the concern is therefore the one that would need to be settled by the concrete convergence check above.","tokens_in":1847,"tokens_out":359,"duration_ms":14406,"concrete_test":"From the manuscript or SI, extract the SIRIUS PW cutoff and k-point settings used for one unary-crystal EOS (e.g., the transition-metal case flagged as pseudopotential-limited); recompute that EOS with 1.5× cutoff and denser k-mesh; if total energy or equilibrium volume shifts by >2 meV/atom or 0.01 Å, the benchmark convergence is insufficient to support the claimed separation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The protocol's error decomposition rests on the three-way comparison (CP2K-Gaussian+GTH vs SIRIUS-PW+GTH vs AE-FP-LAPW). This correctly isolates Gaussian-basis error only if the SIRIUS plane-wave results have reached the complete-basis limit for the fixed GTH pseudopotential (i.e., PW cutoff and k-mesh errors << chemical accuracy) and the AE reference is itself free of significant numerical error. The abstract states \"systematic plane-wave representation\" but supplies no explicit convergence thresholds, residual energies, or cross-checks against higher cutoffs; if those residuals are comparable to the reported CP2K discrepancies, the diagnosis of \"basis-limited\" vs \"pseudopotential-limited\" cases becomes unreliable.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents the UZH protocol, a closed-loop workflow for CP2K that performs molecular calibration of MOLOPT Gaussian basis sets, validates them via unary-crystal equation-of-state benchmarks, and uses three-way comparisons (CP2K Gaussian+GTH vs. SIRIUS PW+GTH vs. all-electron FP-LAPW) to decompose practical errors into Gaussian-basis and pseudopotential components. The diagnosis then guides targeted revisions of basis sets and GTH pseudopotentials, producing improved parameter files as explicit outputs.","tokens_in":1987,"tokens_out":442,"duration_ms":15414,"significance":"If the SIRIUS benchmarks are demonstrably converged, the protocol supplies a constructive, non-circular route to separate and reduce the two dominant numerical approximations in CP2K, directly yielding revised, publicly usable MOLOPT and GTH files that improve both molecular and condensed-phase accuracy. The explicit use of an external all-electron reference lowers circularity relative to purely internal fitting.","major_comments":[{"comment":"The central error-decomposition claim rests on the assumption that SIRIUS plane-wave+GTH calculations have reached the complete-basis limit for the fixed GTH pseudopotential (i.e., PW cutoff and k-mesh residuals ≪ chemical accuracy). The abstract states only that a “systematic plane-wave representation” is used; no cutoff values, k-mesh densities, residual energies, or convergence plots are supplied. Without these data the assignment of cases as “basis-limited” versus “pseudopotential-limited” cannot be verified and the subsequent revisions lack a firm diagnostic foundation.","section":"Abstract (protocol description)"},{"comment":"The molecular calibration and crystal validation steps are described as independent, yet the manuscript does not report whether the revised parameters were tested on an external hold-out set of molecules or crystals outside the calibration loop. If the same systems or similar fitting targets are reused, the reported improvements may partly reflect re-optimization rather than genuine transferability.","section":"Validation section (implied by abstract)"}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive comments, which help clarify the presentation of the UZH protocol. We address each major point below and indicate the revisions we will make.","responses":[{"response":"We agree that explicit convergence data for the SIRIUS PW+GTH calculations are required to substantiate the error decomposition. Although the manuscript describes the use of systematic plane-wave representations and three-way comparisons with all-electron references, the abstract and main text do not tabulate the specific cutoffs, k-meshes, or residual energies. In the revised manuscript we will add these parameters together with convergence plots demonstrating that PW residuals lie well below chemical accuracy for the systems considered. This will allow readers to verify the assignment of basis-limited versus pseudopotential-limited cases.","revision_made":"yes","referee_comment":"[Abstract (protocol description)] The central error-decomposition claim rests on the assumption that SIRIUS plane-wave+GTH calculations have reached the complete-basis limit for the fixed GTH pseudopotential (i.e., PW cutoff and k-mesh residuals ≪ chemical accuracy). The abstract states only that a “systematic plane-wave representation” is used; no cutoff values, k-mesh densities, residual energies, or convergence plots are supplied. Without these data the assignment of cases as “basis-limited” versus “pseudopotential-limited” cannot be verified and the subsequent revisions lack a firm diagnostic foundation."},{"response":"The protocol deliberately separates the calibration (small molecules) from the validation (unary-crystal equations of state) to reduce circularity, and the three-way comparison with all-electron references provides an external anchor. Nevertheless, we acknowledge that the manuscript does not explicitly document performance on a fully disjoint hold-out set. In the revision we will add a dedicated subsection reporting results on an external test set of molecules and crystals not used in the calibration or validation loops, thereby quantifying transferability.","revision_made":"yes","referee_comment":"[Validation section (implied by abstract)] The molecular calibration and crystal validation steps are described as independent, yet the manuscript does not report whether the revised parameters were tested on an external hold-out set of molecules or crystals outside the calibration loop. If the same systems or similar fitting targets are reused, the reported improvements may partly reflect re-optimization rather than genuine transferability."}],"tokens_in":1511,"tokens_out":502,"duration_ms":19041,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that this paper supplies a closed-loop workflow that diagnoses whether a CP2K discrepancy is coming from the atom-centered Gaussian basis or from the norm-conserving pseudopotential, then uses the diagnosis to generate revised MOLOPT sets and GTH pseudopotentials.\n\nWhat is new is the explicit three-way split: production CP2K-GTH, SIRIUS plane-wave with the identical GTH, and all-electron FP-LAPW. That split lets them label systems as basis-limited (noble gases, heavy elements) or pseudopotential-limited (transition metals) and feed the label into the existing optimizers. The outputs are concrete parameter files rather than just error tables.\n\nThe paper does the calibration on small molecules and the validation on unary-crystal equations of state, which is a reasonable scope. Anchoring the references outside CP2K reduces the risk of circular fitting.\n\nThe soft spot is the unshown convergence of the SIRIUS benchmarks. The error separation only holds if the plane-wave cutoff and k-mesh errors are well below the CP2K discrepancies being diagnosed. The abstract calls the plane-wave representation “systematic” but gives no residual energies or cutoff tests. If those residuals are comparable to the reported differences, the basis-versus-pseudopotential assignment becomes unreliable. The all-electron references also need to be demonstrably tight. This is a fixable gap rather than a fatal one.\n\nThe work is aimed at CP2K users who want tighter numerical control without changing the underlying theory. Anyone who maintains basis or pseudopotential libraries would get direct value from the revised files and the diagnostic method.\n\nIt deserves peer review. The idea is practical and the outputs are usable; referees can ask for the missing convergence data and check the revised parameters on a few extra test cases.","headline":"The UZH protocol gives a workable way to split Gaussian-basis error from GTH pseudopotential error in CP2K and then revise the files, but the SIRIUS reference convergence is not shown.","tokens_in":2531,"tokens_out":458,"would_cite":false,"duration_ms":24346,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"The UZH protocol decomposes CP2K errors into separate Gaussian-basis and pseudopotential contributions using three-way comparisons.","keywords":["Gaussian basis sets","pseudopotentials","error decomposition","density functional theory","basis set optimization","equation of state benchmarks","molecular calibration"],"falsifier":"Independent tests on additional molecules or crystals in which the revised basis sets and pseudopotentials produce larger errors than the originals, or in which the three-way comparison misattributes the dominant error source.","tokens_in":2732,"feed_emoji":"","tokens_out":702,"duration_ms":19941,"temperature":0.7,"pith_summary":"The paper describes a workflow that runs production calculations with Gaussian basis sets and norm-conserving pseudopotentials, then compares them to plane-wave calculations that reuse the identical pseudopotential and to all-electron references. The differences isolate whether the dominant remaining error arises from the basis-set representation or from the pseudopotential. Once the limiting source is identified, the workflow feeds the diagnosis into optimizers that revise the parameter files. A reader would care because the combined use of both approximations is common in simulations, and knowing which one sets the accuracy floor allows targeted rather than blanket improvements.","feed_headline":"Protocol splits CP2K errors into basis-set and pseudopotential parts","feed_subtitle":"Three-way benchmarks identify the limiting approximation and produce revised parameter files for molecules and crystals.","key_machinery":"The three-way comparison between production Gaussian-basis calculations, plane-wave calculations that reuse the identical pseudopotential, and all-electron references, which isolates the separate error contributions.","core_discovery":"The UZH protocol is a closed-loop workflow that first calibrates molecularly optimized Gaussian basis sets on small molecules, validates the settings through unary-crystal equation-of-state benchmarks, and then performs a three-way comparison among production CP2K calculations, systematic plane-wave calculations that employ the same Goedecker-Teter-Hutter pseudopotential, and all-electron full-potential linearized augmented-plane-wave references. This comparison decomposes the practical error into a Gaussian-basis component and a pseudopotential component, distinguishes basis-limited cases from pseudopotential-limited cases, and supplies the diagnosis to basis and pseudopotential optimizers","pith_inferences":["Similar closed-loop comparisons could be constructed for other codes that combine atom-centered bases with pseudopotentials.","Repeated application across the periodic table might reveal systematic patterns in which elements require the next round of pseudopotential refinement.","The improved parameters could serve as a more consistent starting point when different simulation packages are compared on the same physical system."],"forward_implications":["The protocol distinguishes basis-limited noble-gas and heavy-element cases from pseudopotential-limited transition-metal cases.","It directs targeted revisions through the existing basis-set and pseudopotential optimizers.","The resulting parameter files are validated for use across both molecular and condensed-phase simulations.","Verification outliers can be converted directly into updated, usable parameter files rather than remaining as diagnostic warnings."],"fun_headline_variants":["UZH protocol isolates basis and pseudopotential errors in CP2K","Three-way benchmarks separate CP2K Gaussian basis errors from GTH pseudopotentials","Closed-loop UZH method refines MOLOPT bases and GTH pseudopotentials","Protocol distinguishes basis-limited from pseudopotential-limited CP2K cases","UZH workflow improves CP2K parameters through error decomposition"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The plane-wave calculations that reuse the same pseudopotential and the all-electron references are accurate enough to correctly attribute the observed discrepancy to either the Gaussian basis or the pseudopotential.","fun_headline_variants_meta":{"raw":{"variants":["UZH protocol isolates basis and pseudopotential errors in CP2K","Three-way benchmarks separate CP2K Gaussian basis errors from GTH pseudopotentials","Closed-loop UZH method refines MOLOPT bases and GTH pseudopotentials","Protocol distinguishes basis-limited from pseudopotential-limited CP2K cases","UZH workflow improves CP2K parameters through error decomposition"]},"model":"grok-4.3","cost_usd":0.004682,"raw_usage":{"total_tokens":2289,"prompt_tokens":779,"num_sources_used":0,"completion_tokens":94,"cost_in_usd_ticks":46815500,"prompt_tokens_details":{"text_tokens":779,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1416,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":779,"tokens_out":94,"duration_ms":10098,"temperature":1.0,"reasoning_tokens":1416,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T11:24:52.588310+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Independent tests on additional molecules or crystals in which the revised basis sets and pseudopotentials produce larger errors than the originals, or in which the three-way comparison misattributes the dominant error source.","supporting_citations":[],"review_version":1}