{"id":"337ef9d1-b4bf-4776-b71e-2cb39504066b","arxiv_id":"2411.13253","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The vDZP basis set, originally built for the ωB97X-3c composite method, also yields accurate and fast DFT results with B97-D3BJ, r2SCAN-D4, B3LYP-D4, M06-2X, and ωB97X-D4 without reparameterization.","lead":"A small, recently introduced basis set called vDZP is shown to work well with several common density functionals, not just the one it was originally designed for. The paper benchmarks these combinations and finds accuracy close to expensive composite methods at much lower cost.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Omitted GMTKN55 subsets are exactly the heavy-element/large-system cases where vDZP could fail; full-benchmark verification is required.","rationale":"The reader's conditional verdict is exactly right, and the omitted-subset concern is the load-bearing one. The paper is honest and internally consistent, and the per-functional comparisons are valid as far as they go; the revMOBH35 and ROT34 results provide independent support for the basis in transition-metal and geometry contexts. What is missing is a direct test of vDZP on the five excluded GMTKN55 subsets. Since vDZP is ECP-heavy and deeply contracted, the excluded subsets are a targeted stress test rather than an arbitrary gap. A concrete check in an independent code would settle the issue in either direction. I would not move the verdict; the conditional status should remain until the full benchmark is provided. The custom fluorine file is a smaller, secondary implementation risk that should be checked alongside.","tokens_in":7129,"tokens_out":5261,"duration_ms":55297,"concrete_test":"Recompute NBPRC, FH51, DC13, C60ISO, and HEAVY28 with the published vDZP basis in a code with an independent ECP implementation (ORCA or PySCF), using the same functionals and grids as Table 1. Recompute full-55 WTMAD2 values for at least B97-D3BJ and r2SCAN-D4. Also verify the custom fluorine basis file against the published vDZP for F by comparing one- and two-electron energy components on a small fluorine-containing molecule. If the full WTMAD2 values change by less than ~0.5 kcal/mol relative to the truncated Table 1 values, the omission is benign; otherwise the claim degrades.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that vDZP is a generally accurate low-cost basis across many functionals—rests almost entirely on the WTMAD2 values in Table 1. Those values omit five GMTKN55 subsets (NBPRC, FH51, DC13, C60ISO, HEAVY28) because of Psi4 ECP errors. These are not randomly omitted: they cover heavy main-group elements, halogen-containing systems, and large conjugated molecules, exactly the systems most likely to expose a flaw in a basis built from large-core ECPs and deep contraction. The manuscript gives no estimate, bound, or independent-code check of vDZP on these subsets. If vDZP errors are concentrated there, the reported WTMAD2 values would be flattering, and the conclusion that vDZP 'can produce highly accurate methods' across main-group chemistry would be overstated. Until a full-55 WTMAD2 is produced with a correct ECP implementation, the breadth of the general-applicability claim is conditional, not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports that the vDZP basis set, originally developed for the composite method ωB97X-3c, can be combined with several unmodified density functionals (B97-D3BJ, r2SCAN-D4, B3LYP-D4, M06-2X, and ωB97X-D4) to yield accuracies that approach those of the much larger (aug)-def2-QZVP basis set and of purpose-built composite methods such as B97-3c and r2SCAN-3c. The evidence consists of GMTKN55 weighted total mean absolute deviations (WTMAD2, with five subsets omitted), revMOBH35 transition-metal barrier heights, ROT34 rotational constants, TorsionNet206 torsional profiles, and n-alkane/perbromo-n-alkane timing studies. The central claim is that vDZP is a generally applicable low-cost basis set that does not require functional- or correction-specific reparameterization.","tokens_in":7306,"tokens_out":3977,"duration_ms":42379,"significance":"If the claim holds, this is a practically important result: it would overturn the usual assumption that double-zeta basis sets are inadequate except inside tightly coupled composite schemes, and it would give computational chemists a cheap, general-purpose basis for main-group thermochemistry, geometries, and conformational energies. The paper's strengths are its use of standard external benchmark sets, its reliance on previously optimized external parameters (vDZP, D4, and the functionals themselves) rather than any fit performed in this work, and the clear head-to-head comparisons against conventional double-zeta basis sets. The main quantitative finding—that vDZP substantially outperforms 6-31G(d), def2-SVP, and pcseg-1 while approaching composite-method accuracy—is useful even if the general-applicability claim is later bounded by the omitted benchmark subsets.","major_comments":[{"comment":"The general-applicability conclusion rests on GMTKN55 WTMAD2 values in Table 1, but five subsets (NBPRC, FH51, DC13, C60ISO, HEAVY28) are omitted because of \"documented errors in Psi4's effective-core-potential implementation.\" These subsets are not random: they contain heavy main-group elements, halogen-containing systems, and large conjugated molecules, exactly the regimes in which a basis built from large-core ECPs and deep contraction could be most fragile. As written, the WTMAD2 numbers are not the full GMTKN55 WTMAD2, and the statement that vDZP is \"generally applicable\" is therefore conditional. The authors should either compute these subsets with a code having a correct ECP implementation (e.g., ORCA, Q-Chem, or a newer Psi4 version) or report per-subset errors/bounds for the omitted sets to demonstrate that no systematic failure is hidden.","section":"Methodology and Table 1"},{"comment":"The paper states that \"a custom basis-set file was used which adds the missing basis functions for fluorine\" because of a documented absence of fluorine in Psi4's internal vDZP implementation. No content or provenance of this custom file is given. Fluorine appears in many GMTKN55 subsets and in the other benchmarks; if the added functions are not exactly the published vDZP fluorine functions, the reported numbers are not vDZP results for those species. The custom basis file, its source, and a verification against an independent implementation should be provided, and the affected subset results should be identified.","section":"Methodology"},{"comment":"The abstract and conclusion claim that vDZP works \"without any method- or correction-specific reparameterization,\" but the ωB97X-D4 row in Table 1 carries the footnote \"utilizes a refit D4 correction.\" If this is the D4 parametrization optimized for ωB97X-3c/vDZP, then that row is not a test of an unmodified functional and should not be cited as evidence of transferability. The claim should be restricted to the four truly unmodified functionals, or the refit should be justified as a standard parameter set rather than a bespoke correction.","section":"Table 1, footnote a"},{"comment":"The paper provides only summary spreadsheets in the Supporting Information and no input files, scripts, or machine-readable protocol. Given that the central claim is a numerical benchmark comparison, the absence of input files for the GMTKN55, revMOBH35, ROT34, and TorsionNet206 calculations—especially the custom fluorine basis file and the exact D4 version/parameters—prevents independent verification and limits the utility of the results as a community resource. The authors should archive the complete input set.","section":"Supporting Information and Methodology"}],"minor_comments":[{"comment":"The text says reference values were obtained with \"(aug)-def2-QZVP,\" but the table column header reads \"def2-QZVP\"; this inconsistency should be resolved.","section":"Results and discussion, Table 1"},{"comment":"There is a typo in the list of omitted subsets: \"HEA VY28\" should be \"HEAVY28,\" and a footnote explaining that Table 1's WTMAD2 excludes these subsets should be added directly to the table.","section":"Methodology"},{"comment":"The footnote \"ωB97X-3x utilizes a refit D4 correction\" appears to contain a typo; it should presumably refer to ωB97X-D4 or ωB97X-3c.","section":"Table 1, footnote a"},{"comment":"The phrase \"remove core elections\" should be \"remove core electrons.\"","section":"Introduction"},{"comment":"The timing section reports that vDZP-based methods were on average 40% slower than the corresponding composite methods, yet Figure 1 and the discussion conclude \"comparable efficiency.\" The 40% slowdown and the hardware-specific nature of the result should be stated more prominently so that the Pareto-efficiency claim is not overstated.","section":"Results and discussion, timing"},{"comment":"For the TorsionNet206 comparison, the reference method used to score the benchmark (CCSD(T)/def2-TZVP) is mentioned, but the protocol for obtaining the MAE values (e.g., single-point energies at a fixed geometry versus relaxed scans) is not fully specified; a sentence clarifying this would improve reproducibility.","section":"Results and discussion, TorsionNet206"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central claim is plausible and potentially valuable, but the two technical gaps—the omitted GMTKN55 subsets and the undocumented custom fluorine basis file—sit exactly where the claim could break. Both are fixable within the scope of the paper, so I recommend major revision rather than rejection. I would also suggest that the editor ask for the complete input archive before a second round of review, because the absence of machine-readable inputs makes the numerical results difficult to audit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper to know about: Wagen and Vandezande test the vDZP basis from ωB97X-3c across four other functionals and find it reproduces most of the accuracy of the large QZVP basis at double-zeta cost, without refitting anything. The evidence supports that for the systems they actually ran. But the headline GMTKN55 numbers are averaged over 50 of the 55 subsets: NBPRC, FH51, DC13, C60ISO, and HEAVY28 were dropped because of Psi4 ECP errors. Those are heavy-element, halogen, and large-system cases, exactly where a deeply contracted, ECP-heavy basis is most at risk. The paper is honest about the omission, but it does not quantify what those subsets would contribute. That makes the general-applicability claim conditional rather than proven.\n\nWhat is genuinely new: this is the first systematic test of vDZP outside the ωB97X-3c framework it was built for. The computational protocol is straightforward and sensible: GMTKN55, revMOBH35, ROT34, TorsionNet206, plus timing on n-alkanes. They compare against def2-QZVP and common double-zeta sets, and vDZP wins clearly. The revMOBH35 and ROT34 results show the basis works for transition-metal barriers and geometries, not just main-group energetics. The TorsionNet206 numbers are a useful practical addition. The timing discussion is also fair: vDZP methods are about 40% slower than the composite methods on alkanes, but faster for heavy-atom systems because of the ECPs.\n\nSoft spots, in order. First, the five omitted subsets are a real hole in the central evidence, and no independent code check or error bound is given. Second, the custom fluorine basis file needed because Psi4 lacks vDZP for F is not documented enough to reproduce, and a bad F function could quietly shift results. Third, the underlying data is promised as a spreadsheet, but no input files or code are provided. These are all fixable, and none of them undermines the main trend in the subsets that were run.\n\nWho this is for: anyone choosing a low-cost basis for DFT on large molecules, and anyone working on composite methods. It deserves a serious referee, but the referee should ask for the missing subsets to be addressed before the general claim is accepted.\n\nRecommendation: send to peer review, with the explicit request that the authors either run the five omitted subsets in another program or show they do not change the conclusion.","headline":"Solid transferability benchmark for vDZP, but the headline GMTKN55 numbers rest on 50/55 subsets — the five omitted are the heavy-element/large-system cases most likely to expose trouble.","tokens_in":7892,"tokens_out":3549,"would_cite":true,"duration_ms":33998,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.Ew"],"model":"deepseek-v4-flash","headline":"The vDZP basis set, originally developed for a single composite method, works across many density functionals without reparameterization, giving near-composite accuracy in thermochemistry, barrier heights, geometries, and torsional…","keywords":["vDZP basis set","density functional theory","composite methods","basis set superposition error","GMTKN55","effective core potentials","transition metal barrier heights","torsional energy profiles"],"falsifier":"Run the five omitted GMTKN55 subsets with vDZP and all five functionals in a quantum chemistry program that implements effective core potentials correctly, then compare the weighted mean absolute deviations back to the (aug)-def2-QZVP references: if any subset's error is dramatically larger than the trends in Tables 1–3, the paper's general-applicability claim fails.","tokens_in":6869,"feed_emoji":"⚛️","tokens_out":8254,"duration_ms":75642,"temperature":0.7,"pith_summary":"Computational chemists have long assumed that accurate density functional theory needs triple-zeta basis sets, and that small-basis methods only work if the functional, basis, and empirical corrections are reparameterized together into composite schemes. This paper claims that vDZP, a polarized double-zeta basis set with effective core potentials and deeply contracted functions originally built for the composite method ωB97X-3c, is not overfit to that one method. When paired with four additional functionals (B97-D3BJ, r2SCAN-D4, B3LYP-D4, and M06-2X), vDZP yields main-group thermochemistry errors close to those of the huge (aug)-def2-QZVP basis and roughly half the errors of conventional double-zeta sets, with no reparameterization beyond the standard D4 dispersion correction. On transition-metal barriers, rotational constants, and drug-like torsional profiles, vDZP methods match or beat fine-tuned composite methods. If right, vDZP gives a general low-cost basis that breaks the speed-accuracy tradeoff that motivated bespoke composite schemes.","feed_headline":"A double-zeta basis set rivals fine-tuned composite DFT methods","feed_subtitle":"No reparameterization needed: vDZP matches composite accuracy across four benchmark families.","key_machinery":"vDZP is a polarized valence double-zeta basis set whose defining features are large-core effective core potentials that remove core electrons, deeply contracted valence basis functions, and parameters optimized on molecular systems rather than free atoms. These features suppress the two classic failures of small basis sets, basis-set incompleteness error (the density is too rigid) and basis-set superposition error (fragments borrow each other's basis functions), down to near triple-zeta levels. The basis thus acts as a drop-in replacement: any functional can be combined with it, and the only extra ingredient needed is the now-standard D4 empirical dispersion correction for functionals that lack dispersion.","core_discovery":"At the paper's center is the demonstration that vDZP is a general-purpose basis rather than a bespoke component. The authors combine vDZP with five functionals spanning GGA, meta-GGA, hybrid, and range-separated hybrid classes, and evaluate them on GMTKN55 (minus five subsets that the software's effective-core-potential implementation cannot handle), revMOBH35, ROT34, and TorsionNet206. The result: weighted mean absolute deviations for vDZP are only moderately worse than those of (aug)-def2-QZVP, clearly better than 6-31G(d), def2-SVP, and pcseg-1, and comparable to the purpose-built composite methods B97-3c, r2SCAN-3c, and ωB97X-3c. For geometry predictions, vDZP-based methods actually produce the lowest mean deviations among the low-cost methods tested, and for torsions they are within 0.02–0.06 kcal/mol of triple-zeta hybrids. The paper argues this breaks the assumption that the tradeoff between speed and accuracy can only be resolved by tight coupling of methods, basis sets, and empirical corrections.","pith_inferences":["The paper's omission of five GMTKN55 subsets means the claim of general applicability is directly tested only for the main-group organic subset; a user targeting heavy-element or large-system chemistry should look for a dedicated test before relying on vDZP.","A natural next experiment is to run vDZP with a functional not in the test set, such as a double-hybrid, to see if the 'no reparameterization' pattern holds when the functional itself has a different error profile.","The reported timing data suggest vDZP's real speed advantage will emerge only after integral codes are optimized for deeply contracted, low-angular-momentum basis functions; until then, practical speed parity with composite methods is partly hardware- and software-dependent.","If the omitted heavy-element subsets do degrade, the paper's stronger claim of high generality would fall back to 'highly effective for organic main-group chemistry,' a useful but narrower result."],"forward_implications":["vDZP can replace conventional double-zeta bases such as 6-31G(d), def2-SVP, and pcseg-1 for routine DFT, roughly halving weighted mean absolute errors on GMTKN55 without increasing runtime.","For drug-design workflows, torsional scans with B97-D3BJ/vDZP or r2SCAN-D4/vDZP approach triple-zeta accuracy at a fraction of the cost, making broad conformer screening feasible.","The five functionals tested all retain most of their large-basis accuracy with vDZP, so users are not locked into a single bespoke composite method.","For heavy-element-rich systems, vDZP's extensive use of effective core potentials can make it substantially faster than composite methods, with a 3.4-fold speedup observed for perbromo-n-pentane.","Basis-set design that optimizes molecular performance directly, rather than atomic energies, is a promising route to further Pareto improvements in quantum chemistry."],"supporting_citations":[{"why":"Defines the ωB97X-3c composite method and its refit D4 correction, the reference point that vDZP was designed to serve.","marker":"[13]"},{"why":"Reports the vDZP basis set itself, with its ECP-heavy, deeply contracted, molecule-optimized design that the authors test.","marker":"[14]"},{"why":"Supplies the GMTKN55 main-group thermochemistry benchmark set used for the core accuracy comparisons.","marker":"[18]"},{"why":"Provides the B97-3c composite method baseline that the new B97-D3BJ/vDZP combination is compared against.","marker":"[11]"},{"why":"Provides the r2SCAN-3c composite method baseline for the r2SCAN-D4/vDZP comparison.","marker":"[12]"},{"why":"Supplies the ROT34 rotational-constant benchmark used to test geometries of optimized structures.","marker":"[21]"},{"why":"Provides the TorsionNet206 drug-like torsional energy benchmark used to test conformational accuracy.","marker":"[23]"},{"why":"Supplies the revMOBH35 transition-metal barrier-height benchmark used to test organometallic accuracy.","marker":"[19]"}],"fun_headline_variants":["vDZP: universal basis, composite accuracy","No tuning needed: vDZP matches composite accuracy","Off-the-shelf vDZP rivals bespoke composite DFT methods","vDZP competes with composite methods without per-method tuning","vDZP: one basis for many functionals, composite accuracy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that vDZP is generally applicable depends on the five omitted GMTKN55 subsets (NBPRC, FH51, DC13, C60ISO, and HEAVY28) not hiding a systematic failure, particularly for heavy elements and large systems where an ECP-heavy, deeply contracted basis might be most fragile.","fun_headline_variants_meta":{"raw":{"variants":["vDZP: universal basis, composite accuracy","No tuning needed: vDZP matches composite accuracy","Off-the-shelf vDZP rivals bespoke composite DFT methods","vDZP competes with composite methods without per-method tuning","vDZP: one basis for many functionals, composite accuracy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001406,"raw_usage":{"total_tokens":5656,"prompt_tokens":894,"completion_tokens":4762,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":4673}},"tokens_in":510,"tokens_out":4762,"duration_ms":28826,"temperature":1.0,"reasoning_tokens":4673,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:38:24.926448+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the five omitted GMTKN55 subsets with vDZP and all five functionals in a quantum chemistry program that implements effective core potentials correctly, then compare the weighted mean absolute deviations back to the (aug)-def2-QZVP references: if any subset's error is dramatically larger than the trends in Tables 1–3, the paper's general-applicability claim fails.","supporting_citations":[{"cited_title":"B97X-3c: A composite range-separated hybrid DFT method with a molecule-optimized polarized valence double- basis set","cited_arxiv_id":null,"evidence_quote":"Defines the ωB97X-3c composite method and its refit D4 correction, the reference point that vDZP was designed to serve."},{"cited_title":"Optimal Small Basis Set and Geometric Counterpoise Correction for DFT Computations","cited_arxiv_id":null,"evidence_quote":"Reports the vDZP basis set itself, with its ECP-heavy, deeply contracted, molecule-optimized design that the authors test."},{"cited_title":"A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry , kinetics and noncovalent interactions","cited_arxiv_id":null,"evidence_quote":"Supplies the GMTKN55 main-group thermochemistry benchmark set used for the core accuracy comparisons."},{"cited_title":"G.; Bannwarth, C.; Hansen, A.; Grimme, S","cited_arxiv_id":null,"evidence_quote":"Provides the B97-3c composite method baseline that the new B97-D3BJ/vDZP combination is compared against."},{"cited_title":"Swiss army knife","cited_arxiv_id":null,"evidence_quote":"Provides the r2SCAN-3c composite method baseline for the r2SCAN-D4/vDZP comparison."},{"cited_title":"Implementation of nuclear gradients of range-separated hybrid density functionals and benchmarking on rotational constants for organic molecules","cited_arxiv_id":null,"evidence_quote":"Supplies the ROT34 rotational-constant benchmark used to test geometries of optimized structures."},{"cited_title":"A machine learning-based high-precision density functional method for drug-like molecules","cited_arxiv_id":null,"evidence_quote":"Provides the TorsionNet206 drug-like torsional energy benchmark used to test conformational accuracy."},{"cited_title":"A.; Janes, T","cited_arxiv_id":null,"evidence_quote":"Supplies the revMOBH35 transition-metal barrier-height benchmark used to test organometallic accuracy."}],"review_version":1}