pith. machine review for the scientific record. sign in

arxiv: 2603.15788 · v2 · submitted 2026-03-16 · ⚛️ physics.chem-ph

Recognition: 1 theorem link

· Lean Theorem

On the performance of QTP functionals applied to second-order response properties II: Dynamic polarizability and long-range C₆ coefficients

Rodrigo A. Mendes , Peter R. Franke , Ajith Perera , Rodney J. Bartlett
Authors on Pith no claims yet

Pith reviewed 2026-05-15 09:50 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords dynamic polarizabilityC6 dispersion coefficientsdensity functional theoryQTP functionalsEOM-CCSDCasimir-Polder equationresponse propertieshybrid functionals
0
0 comments X

The pith

TPSS0 and QTP01 rank highest among Kohn-Sham functionals for dynamic polarizabilities while O3LYP leads for C6 coefficients.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper evaluates twenty-five density functionals for computing dynamic polarizabilities at five wavelengths and C6 dispersion coefficients via the Casimir-Polder equation. It benchmarks the results against EOM-CCSD reference values and identifies the strongest performers. A reader would care because these frequency-dependent properties control optical responses and long-range forces that become prohibitively expensive to compute with wavefunction methods for larger molecules. The study shows that selected hybrid functionals deliver accuracy close to the reference at practical cost.

Core claim

For dynamic polarizabilities at wavelengths 632.99 nm, 594.10 nm, 543.52 nm, 514.50 nm, and 325.13 nm, TPSS0 and QTP01 give the best overall agreement with EOM-CCSD among Kohn-Sham calculations, and the QTP functionals reproduce the pole structure of CO adequately. For C6 coefficients obtained from the Casimir-Polder equation, O3LYP ranks first overall, yet the first eleven functionals show very similar accuracy, and within the QTP family QTP01 and LC-QTP perform best.

What carries the argument

Linear-response time-dependent Kohn-Sham calculations with 25 exchange-correlation functionals benchmarked to EOM-CCSD, followed by Casimir-Polder integration for C6 coefficients.

If this is right

  • TPSS0 and QTP01 can be used in place of EOM-CCSD for routine dynamic polarizability calculations at visible and near-UV wavelengths.
  • O3LYP and the other top-ranked functionals deliver practical accuracy for C6 dispersion coefficients.
  • The narrow spread among the leading eleven functionals for C6 indicates that functional choice is less critical once a good hybrid is selected.
  • QTP01 emerges as a consistent performer across both dynamic polarizabilities and C6 coefficients within its family.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The close match to CC3 at lower frequencies supports applying these same functionals to other frequency-dependent response properties such as hyperpolarizabilities.
  • Good reproduction of CO poles suggests the QTP functionals may handle certain excited-state features more reliably than their ground-state design implies.
  • Repeating the ranking on a broader set of molecules would test whether the current ordering generalizes beyond the systems examined.

Load-bearing premise

EOM-CCSD supplies a sufficiently accurate reference to rank the functionals, including near excitation poles where CC3 and EOM-CCSD begin to differ.

What would settle it

Gas-phase experimental dynamic polarizabilities measured for CO at 325 nm would show whether the functional rankings and pole positions match real data.

read the original abstract

This work is the second in the series "On the performance of QTP functionals applied to second-order response properties." In the first paper (J. Chem. Phys. 162, 054105, 2025), we demonstrated the good performance of Quantum Theory Project functionals in predicting static perturbed second-order properties, such as static polarizabilities, nuclear magnetic resonance (NMR) spin-spin coupling constants, and NMR chemical shifts. In the present study, we focus on frequency-dependent properties, namely dynamic polarizabilities and C$_6$ dispersion coefficients. For completeness, a total of 25 exchange-correlation (XC) functionals were investigated. Dynamic polarizabilities were evaluated at five different perturbation wavelengths: 632.99 nm, 594.10 nm, 543.52 nm, 514.50 nm, and 325.13 nm. This property was also computed using HF and EOM-CCSD. In general, EOM-CCSD results are very close to those obtained with linear-response CC3, except at the highest frequency. Among Kohn-Sham calculations, TPSS0 and QTP01 showed the best overall performance for dynamic polarizabilities. We also assessed how well QTP functionals reproduce the pole structure of the CO molecule. For the C$_6$ dispersion coefficients, calculations were performed using the Casimir-Polder equation. The best overall performance was obtained with O3LYP; however, the first eleven ranked functionals show very similar accuracy. Within the QTP family, QTP01 and LC-QTP provide the best results for C$_6$ coefficients.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. This manuscript is the second in a series on QTP functionals for second-order response properties. It benchmarks 25 XC functionals for dynamic polarizabilities at five wavelengths (632.99 nm to 325.13 nm) against HF and EOM-CCSD, notes EOM-CCSD closeness to linear-response CC3 except at the highest frequency, assesses QTP reproduction of the CO pole structure, and computes C6 coefficients via the Casimir-Polder equation. Among KS methods, TPSS0 and QTP01 show best overall performance for polarizabilities; O3LYP ranks highest overall but the top eleven functionals have very similar accuracy. Within the QTP family, QTP01 and LC-QTP perform best for C6 coefficients.

Significance. If the rankings hold, the work supplies a practical benchmark extending the prior static-property study to frequency-dependent cases important for optical response and dispersion. Direct multi-wavelength comparisons to wavefunction references and the explicit CO pole-structure check are positive features that aid functional selection in response-property calculations.

major comments (1)
  1. [Dynamic polarizabilities and pole-structure assessment] Dynamic polarizabilities section (EOM-CCSD vs. CC3 comparison): The performance rankings identifying TPSS0 and QTP01 as best among KS methods, and O3LYP overall, rest on EOM-CCSD as primary reference. The manuscript acknowledges EOM-CCSD differs from CC3 at 325.13 nm and evaluates the CO pole structure, yet the small accuracy spread among the top eleven functionals means any non-uniform reference error near resonances could reorder the claimed best performers. Direct inclusion of CC3 values for the leading functionals or a sensitivity analysis would confirm ranking stability.
minor comments (1)
  1. [Abstract] Abstract: the statement that 'the first eleven ranked functionals show very similar accuracy' would be clearer with a quantitative metric (e.g., maximum error range or reference to a supporting table).

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive overall assessment and the constructive comment regarding reference-method robustness. We address the major comment below.

read point-by-point responses

  1. Referee: Dynamic polarizabilities section (EOM-CCSD vs. CC3 comparison): The performance rankings identifying TPSS0 and QTP01 as best among KS methods, and O3LYP overall, rest on EOM-CCSD as primary reference. The manuscript acknowledges EOM-CCSD differs from CC3 at 325.13 nm and evaluates the CO pole structure, yet the small accuracy spread among the top eleven functionals means any non-uniform reference error near resonances could reorder the claimed best performers. Direct inclusion of CC3 values for the leading functionals or a sensitivity analysis would confirm ranking stability.

    Authors: We appreciate the referee pointing out the potential sensitivity of the rankings to the choice of reference near resonance. The manuscript already states that EOM-CCSD is very close to linear-response CC3 except at the highest frequency (325.13 nm) and includes an explicit check of the CO pole structure. Given the small spread among the top eleven functionals, we agree that additional confirmation is useful. We will therefore add CC3 values (where computationally feasible) for the leading functionals TPSS0, QTP01, and O3LYP at 325.13 nm, together with a short sensitivity discussion comparing rankings obtained with EOM-CCSD versus CC3. This will be included in the revised manuscript to demonstrate that the ordering remains stable. revision: partial

Circularity Check

1 steps flagged

Minor self-citation to prior series paper; rankings rest on independent EOM-CCSD/CC3 benchmarks

specific steps
  1. self citation load bearing [Abstract]
    "This work is the second in the series... In the first paper (J. Chem. Phys. 162, 054105, 2025), we demonstrated the good performance of Quantum Theory Project functionals in predicting static perturbed second-order properties"

    The reference supplies only background continuity for the QTP family; the current results on dynamic properties and C6 coefficients are obtained by direct comparison to EOM-CCSD and CC3, which are independent of the cited prior work.

full rationale

The paper performs a standard benchmarking study of 25 XC functionals against external high-level references (EOM-CCSD, linear-response CC3) for dynamic polarizabilities at five wavelengths and C6 coefficients via Casimir-Polder. No derivation chain exists that reduces to fitted inputs, self-definitions, or author-specific uniqueness theorems. The single self-citation to the 2025 JCP predecessor paper supplies series context only and is not invoked to justify any performance claim or to forbid alternatives. All reported rankings derive from direct numerical comparison to the cited ab initio methods, which lie outside the present work's fitted values or normalizations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

No new free parameters, axioms beyond standard DFT, or invented entities; the study relies on existing functionals and established reference methods.

axioms (2)
  • domain assumption Kohn-Sham DFT framework for linear response properties
    Invoked for all XC functional calculations of polarizabilities and dispersion coefficients.
  • domain assumption EOM-CCSD as reliable benchmark reference
    Used to assess functional performance and pole structure.

pith-pipeline@v0.9.0 · 5614 in / 1240 out tokens · 35722 ms · 2026-05-15T09:50:11.987969+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.