{"id":"44334892-3654-4437-af7c-41c3961bdedf","arxiv_id":"2606.12272","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Plane-wave orbital-optimized DFT gives reliable dipole moments for diffuse Rydberg states where atom-centered basis sets fail, and PBE0 is the best functional tested.","lead":"This computational chemistry paper tests how well orbital-optimized density functional theory predicts the electric dipole moments of Rydberg excited states. It finds these properties are far more sensitive to the choice of basis set than excitation energies, with plane waves outperforming atom-centered basis sets.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No box-size convergence test supports the plane-wave reference for the most diffuse Rydberg states; the claimed LCAO sign errors may be a PW artifact.","rationale":"The paper's central claim is that dipole moments are a stricter test than excitation energies and that plane waves reveal failures of atom-centered basis sets. The most convincing evidence would be an independent, converged reference. Instead, the PW calculations are presented without any convergence study of the simulation cell size or kinetic-energy cutoff. The reader's weakest assumption identified exactly this missing validation. I agree that this is the load-bearing concern: if PW results are biased by confinement, then the sign and magnitude differences attributed to LCAO inflexibility could be artifacts. The paper's own discussion concedes that high-level references are unavailable for the most diffuse states, which makes the PW reference even more critical. I do not see a fatal flaw in the paper's overall methodology; the qualitative pattern (aug underestimates variance, d-aug improves but not always) is internally consistent and the composite-basis experiment provides supporting evidence for genuine LCAO flexibility limitations. However, those internal comparisons still use the unconverged PW values as the target. A box-size convergence test is cheap and would either validate or invalidate the main qualitative conclusion. Since the reader already conditioned the verdict on this issue, the verdict should remain unchanged pending that test.","tokens_in":30439,"tokens_out":5532,"duration_ms":65328,"concrete_test":"Repeat the OO-PBE plane-wave calculations for water S4/S5, ammonia S3, and methanol S3/S4 using the same 1200 eV cutoff and 0.16 Å grid but larger simulation cells with 13 Å and 16 Å of vacuum. Record dipole components and σ(r). If any reported value shifts by more than 0.15 D (or changes sign), the published PW reference is not converged with respect to box size; the LCAO-vs-PW comparisons in Section 3.1 would need to be repeated with a converged cell. Also test one state with a 1500 eV cutoff to rule out cutoff effects.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 2.4 the PW setup is specified only as 'a minimum of 10.5 Å of vacuum' and a 0.16 Å grid, with no convergence test with respect to cell size or cutoff. The central qualitative result—that atom-centered aug/d-aug bases give large dipole errors, even sign errors, for diffuse Rydberg states—is established by comparing against this PW reference, not against external high-level data: the states driving the conclusion (water S4/S5, ammonia S3, methanol S3/S4) have no CC reference in Table 2, and the paper itself argues that CC/d-aug references can be unreliable for the most diffuse states (Section 4). Thus the PW result is the only anchor for the most dramatic claims. A 10.5 Å vacuum corresponds to roughly 20 bohr from atom to cell boundary; the most diffuse states have σ(r) ≈ 80 bohr² (RMS radius ≈ 9 bohr), so the Rydberg density tail can reach the cell boundary. Without a box-size test, the reported PW dipole components (e.g., water S4 µz = −1.13 D vs d-aug +0.03 D) could be biased by density truncation or periodic-image effects. That would undermine the claim that LCAO is deficient rather than PW being unconverged. The paper's validation of PW against CC values is also selective: only states where d-aug and PW already agree are used, which cannot validate the pathological cases.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports orbital-optimized (OO) density functional calculations of excited-state dipole moments for Rydberg states of water, formaldehyde, ammonia, and methanol. The authors compare plane-wave/PAW and atom-centered (aug-cc-pVDZ+sz, d-aug-cc-pVDZ+sz) basis sets with PBE, PBE0, and globally scaled Perdew–Zunger SIC functionals, and benchmark against high-level coupled-cluster reference values where available. The central claim is that dipole moments are a much stricter test than excitation energies for OO calculations of diffuse Rydberg states: a single-augmented basis can give large errors or even sign errors in the dipole moment while the excitation energy remains close to the plane-wave result, and residual errors persist with double augmentation for the most diffuse states. The authors further report that PBE/PW agrees well with reference values, PBE0 improves the statistics, and global SIC overestimates dipole magnitudes despite restoring the correct asymptotic potential.","tokens_in":30748,"tokens_out":8721,"duration_ms":98889,"significance":"If the central claim holds, the paper establishes an important new benchmark dimension for OO-DFT Rydberg-state calculations, showing that energy-only validation is insufficient for property predictions. The study is original, uses multiple molecules and functionals, provides state-specific analysis via the variance of the electron density and difference densities, and is implemented in the open-source GPAW code. The plane-wave/PAW approach to diffuse Rydberg dipoles is a useful methodological contribution. However, the strongest conclusions rely on the plane-wave results being converged with respect to simulation cell size, and that convergence is not demonstrated; the functional ranking also rests on a selected subset of states. These issues need to be resolved before the quantitative conclusions can be fully trusted.","major_comments":[{"comment":"The plane-wave results are used as the reference for the most dramatic claims (e.g., sign errors in water S4/S5, ammonia S3, methanol S3/S4), but no convergence test is reported for the simulation cell size or the grid/cutoff. The most diffuse states have σ(r) ≈ 80–90 bohr², i.e., an RMS radius of ≈9 bohr, while the stated 'minimum of 10.5 Å of vacuum' puts the cell boundary only ≈20 bohr from the atoms. The Rydberg tail can therefore be sensitive to periodic confinement, and the dipole moment is especially sensitive to the tail density. Since most of the states driving the conclusion have no coupled-cluster reference in Table 2, the PW values are the sole anchor. Please add explicit box-size convergence tests (e.g., vacuum of 10.5, 12, and 14 Å) and cutoff tests for the most diffuse states, and report whether the qualitative LCAO-vs-PW differences survive. Also state precisely how the '","section":"§2.4 / Tables 1–2"},{"comment":"The statistical comparison of functionals in Figure 5 is restricted to states for which the d-aug and PW PBE results differ by less than 5%. This excludes precisely the most diffuse states where basis-set flexibility is the issue and where functional dependence is largest (water S4/S5, ammonia S3, methanol S3/S4). The reported median absolute errors (PBE ≈20%, PBE0 ≈6%, SIC ≈24–26%) therefore characterize a favorable subset, not the Rydberg states that are the paper's main focus. Please report the exact set and number of states used in the statistics, and provide a reference-free analysis (e.g., PW-vs-LCAO differences as a function of functional) to support the claim that PBE0 improves and SIC worsens dipole moments for diffuse Rydberg states.","section":"§3.2 / Fig. 5"}],"minor_comments":[{"comment":"The spin-purification formula for the dipole moment is stated without derivation or citation. Since it is not immediately implied by the energy formula Eq. (4), please provide a justification or a reference showing that the mixed-spin density is the average of the singlet and triplet densities in the OO framework.","section":"§2.3, Eq. (8)"},{"comment":"The text says the ammonia S3 dipole is overestimated by 'more than 4.5 D' with aug, but the table gives aug = 4.79 D and PW = 1.07 D, i.e., a difference of 3.72 D. Please correct this numerical value.","section":"§3.1 / Table 1"},{"comment":"The statement about two PBE-SIC solutions for the formaldehyde 2p_y→3p_z state gives spin-purified dipole moments of 0.03 and 0.72 D, but Table 2 reports PBE-SIC S3 as 1.21 D. The numbers 0.03 and 0.72 appear to belong to a different state (S4). Please clarify the state labels and the values reported in Table 2.","section":"§4"},{"comment":"There are several typographical errors: 'Unfortunatley' in §3.2, 'The xcomponent' in §3.1, and inconsistent use of 'SIC/2' vs 'PBE-SIC/2' in a few places. A careful proofreading pass is needed.","section":"§3.2 / §4"},{"comment":"The ammonia reference values are attributed to 'personal correspondence' with the authors of Ref. [12]. For reproducibility, please include these reference values and the underlying calculations in the Supporting Information or make them publicly available in a persistent form.","section":"Table 2 / Ref. [85]"},{"comment":"The conclusion states that 'PWs are not affected by confinement effects.' This is only true if the cell is sufficiently large; the manuscript itself acknowledges this in the Introduction. Please rephrase to 'PWs are not affected by confinement effects for the cell sizes used here, as verified by ...' once the convergence test is added.","section":"§5"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable after the requested convergence tests are added. The central qualitative claim is plausible and the dataset is valuable, but the missing box-size convergence test for the plane-wave reference is load-bearing because the most diffuse states have no independent high-level reference. The numerical inconsistencies in §3.1 and §4 should be checked carefully; if they are typographical, they can be fixed in revision, but if they reflect a data handling issue, the tables need re-verification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the punchline: this is a useful benchmark paper, and its central claim holds up. Dipole moments are a stricter test than excitation energies for OO-DFT Rydberg states, and a well-chosen atom-centered basis can give large dipole errors—even the wrong sign—while the excitation energy is nearly converged. That's worth knowing.\n\nWhat's new: prior OO-DFT Rydberg benchmarks mostly reported energies, with dipoles only for the lowest state. Here they cover multiple higher Rydberg states, use a plane-wave representation in GPAW, and systematically compare against aug- and d-aug basis sets. The analysis is careful: they use the variance σ(r) to track the spatial extent, and the orbital-density plots make the failure mode concrete. The conclusion that plane waves are a better choice for diffuse Rydberg states is practical and probably right.\n\nThe functional comparison is also sensible as a first pass: PBE0 has the smallest median error, PBE is okay, and globally scaled SIC overestimates dipole moments. The SIC side is honest about convergence problems.\n\nNow the soft spots, in proportion. The big one: the plane-wave reference is not tested for box-size convergence. The text states a minimum of 10.5 Å of vacuum and a 0.16 Å grid, but no cutoff or box-size study. For states with σ around 80 bohr², the Rydberg tail reaches out to ~9 bohr RMS, and a box edge at ~20 bohr could bias the dipole. The most dramatic cases—water S4 and S5, methanol S3/S4, ammonia S3—have no CC reference, so the PW value is the only anchor. If a larger box shifts those dipoles significantly, the reported LCAO sign errors could turn out to be artifacts. This is fixable: a box-size convergence test for the most diffuse state of each molecule is a few extra calculations.\n\nThe functional-error medians (PBE0 ~6%) are computed on a filtered subset, and the formaldehyde S3 SIC result was selected from two solutions. That's not disqualifying, but it means the quantitative rankings are less robust than the qualitative message.\n\nThe TDDFT comparison is suggestive, but the sets aren't identical; worth a careful note.\n\nOverall: solid, honestly written, aware of its limitations. The authors explicitly say that CC references with atom-centered bases may be unreliable for the most diffuse states—they know the issue. They just need to apply the same scrutiny to their own PW reference. With the box-size test added, this is a good benchmark that deserves publication. I would send it to peer review; the referee should ask for that test and for a clear statement about which PW results are robust enough to anchor the sign-error claim.\n\nMy own verdict: this paper is for people working on OO-DFT or Rydberg properties. I'd bring it to reading group, and I'd cite it once I'm comfortable with the box-size question.","headline":"A worthwhile OO-DFT benchmark: dipole moments of Rydberg states are much harder to converge than energies, and plane waves expose LCAO limits—but the plane-wave reference itself lacks a box-size test, so the dramatic sign-error claims are not yet fully anchored.","tokens_in":31245,"tokens_out":5024,"would_cite":true,"duration_ms":60292,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dipole moments are a stricter test than excitation energies for orbital-optimized DFT calculations of Rydberg states, and plane-wave bases reveal biases hidden in atom-centered basis sets.","keywords":["Rydberg states","excited-state dipole moment","orbital-optimized DFT","plane-wave basis","basis-set effects","self-interaction correction","PBE0","spin purification"],"falsifier":"Recompute the most diffuse Rydberg states (e.g., water S4 and S5) with a much larger simulation cell (e.g., 15–20 Å vacuum) and a tighter grid; if the plane-wave dipole moments shift by more than ~0.1 D relative to the 10.5 Å results, the claimed PW reference would be compromised.","tokens_in":30348,"feed_emoji":"⚛️","tokens_out":2661,"duration_ms":31305,"temperature":0.7,"pith_summary":"This paper shows that a common atom-centered basis set can badly misestimate the dipole moments of Rydberg excited states—sometimes even flipping their sign—while the excitation energy looks fine. Plane waves, being delocalized, avoid the overconfinement of diffuse Rydberg orbitals and provide a more reliable reference. The authors test several density functionals and find PBE0 performs best, whereas explicit self-interaction corrections, despite fixing the long-range potential, systematically overestimate dipole moments. If correct, this means benchmarking Rydberg-state calculations on excitation energies alone is insufficient for property prediction.","feed_headline":"Dipole moments expose basis-set bias hidden in Rydberg excitation energies","feed_subtitle":"A common atom-centered basis can flip the sign of a Rydberg state's dipole even when its energy looks right; plane waves fix it.","key_machinery":"The central object is the excited-state dipole moment evaluated from state-specific orbital-optimized densities, with open-shell singlet values obtained by spin purification (µ_S = 2µ_M − µ_T). The key tool is the plane-wave/projector-augmented-wave (PAW) representation, which provides a flexible, non-atom-centered description of diffuse Rydberg orbitals. The basis-set comparison is aided by the variance of the electron position operator, σ(r), which quantifies the spatial extent of the excited-state density and shows that matching σ(r) does not guarantee an accurate dipole moment.","core_discovery":"Using orbital-optimized DFT with a plane-wave basis, the authors compute dipole moments for a set of Rydberg excited states of water, formaldehyde, ammonia, and methanol. They find that the dipole moment is far more sensitive to basis-set choice than the excitation energy. A singly-augmented atom-centered basis (aug-cc-pVDZ) overconfines the Rydberg density, causing large magnitude errors and sometimes wrong orientation of the dipole, even when the excitation energy is nearly basis-set independent. Adding a second diffuse function set (d-aug) improves radial extent but leaves persistent errors for the most diffuse states, because the atom-centered form cannot fully capture the anisotropic de","pith_inferences":["The same atom-centered basis bias likely affects other properties (e.g., oscillator strengths, polarizabilities) of diffuse excited states, not just dipole moments.","Locally scaled self-interaction corrections, which reduce the correction in regions of overlapping orbital densities, could plausibly resolve the overestimation seen here, but this is a testable conjecture beyond the paper.","The finding that matching σ(r) does not ensure accurate dipole moments suggests that anisotropic basis flexibility—not just diffuseness—should be a design criterion for future basis sets for Rydberg states.","The spin-purification procedure assumes the mixed-spin and triplet states have identical spatial orbitals; the paper notes SIC-induced symmetry breaking that prevents purification in some cases, implying this assumption can fail and may affect other states where it was applied."],"forward_implications":["Benchmarking OO-DFT Rydberg states solely by excitation energies is insufficient; dipole moments (and likely other one-electron properties) expose basis-set deficiencies hidden by variational energy stationary points.","Plane-wave or otherwise delocalized basis representations should be preferred for computing properties of diffuse Rydberg states, especially when accurate reference values are sought.","Adding more diffuse Gaussian functions does not systematically eliminate dipole errors; the atom-centered anchoring itself is a limitation, not just the radial extent.","PBE0 emerges as a reliable low-cost functional for Rydberg-state dipole moments, while a globally scaled self-interaction correction, despite its correct asymptotic potential, systematically overestimates dipole magnitudes.","High-level coupled-cluster references that use atom-centered bases may themselves be biased for the most diffuse Rydberg states, so new benchmark data with flexible basis representations are needed."],"fun_headline_variants":["Excitation energy hides basis-set errors that flip dipole moments","Plane waves fix dipole moments that atom-centered bases get wrong","Rydberg dipole moments more sensitive to basis than excitation energy","Orbital-optimized DFT: plane waves beat atom-centered for Rydberg dipoles","Basis-set overconfinement distorts Rydberg dipole moments"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The plane-wave calculations are treated as the converged reference without a demonstrated test that the chosen vacuum size and grid spacing are fully converged for dipole moments, and the coupled-cluster references themselves use atom-centered basis sets that the paper argues can be unreliable for the most diffuse states.","fun_headline_variants_meta":{"raw":{"variants":["Excitation energy hides basis-set errors that flip dipole moments","Plane waves fix dipole moments that atom-centered bases get wrong","Rydberg dipole moments more sensitive to basis than excitation energy","Orbital-optimized DFT: plane waves beat atom-centered for Rydberg dipoles","Basis-set overconfinement distorts Rydberg dipole moments"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00057,"raw_usage":{"total_tokens":2588,"prompt_tokens":855,"completion_tokens":1733,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":1639}},"tokens_in":599,"tokens_out":1733,"duration_ms":12386,"temperature":1.0,"reasoning_tokens":1639,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T11:43:38.862854+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the most diffuse Rydberg states (e.g., water S4 and S5) with a much larger simulation cell (e.g., 15–20 Å vacuum) and a tighter grid; if the plane-wave dipole moments shift by more than ~0.1 D relative to the 10.5 Å results, the claimed PW reference would be compromised.","supporting_citations":[],"review_version":2}