{"id":"9ef3991e-eb01-499d-ad2b-5c678bc00b03","arxiv_id":"2607.18071","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"DFT+U+V is equivalent, to first order in V/U, to DFT+U with density-dependent Hubbard projectors hybridized with neighboring sites.","lead":"This paper shows that the inter-site +V correction in DFT+U+V is mathematically the same, to first order in V/U, as an on-site DFT+U correction using projectors rotated to mix in neighboring orbitals. The result gives a formal basis for why +V improves covalent bonding, and it questions how V should be computed.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (5)'s rotation generator uses the conjugated inter-site occupancy: for complex Λ the first-order equivalence fails.","rationale":"The reader's accepted verdict identified the orthonormal-projector assumption as the weakest point, but there is a more direct algebraic defect in the central formula. The paper's main theorem depends on Eq. (5) producing the correct rotation S; for complex inter-site occupations — which occur in spin-polarized systems away from Γ, in the presence of spin-orbit coupling, or in general k-space — the formula as written does not satisfy the defining commutator relation (Eq. 20) and does not reproduce the DFT+U+V energy to first order. The numerical NiO test at Γ with real off-diagonal elements cannot detect this. The conceptual claim ('+V is +U on redefined projectors') is likely recoverable by conjugating the numerator in Eq. (5), but as published the explicit generator is incorrect, and the derivation in the SI inherits the inconsistency. A conditional acceptance requiring this correction and a re-verification of the proof is appropriate. This is a specific algebraic flaw, not a disagreement with consensus or a personal criticism.","tokens_in":22168,"tokens_out":39248,"duration_ms":307677,"concrete_test":"Take a 2-site, 1-orbital model with λ_i=0.8, λ_j=0.2, U_I=7 eV, U_J=0, V=1 eV, and off-diagonal occupancy Λ=i0.5 (pure imaginary). Compute to first order the energy change from the paper's Eq. (5): ΔE = Tr[V_U[ρ,S]] = -V Re(Λ²) = +0.25 eV, while E^V = -V|Λ|² = -0.25 eV. The claimed equality fails. Repeating with the corrected numerator Λ^{IJ}_{ij} gives ΔE = E^V. This is a few-line symbolic check.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Eq. (11) defines the inter-site potential as V_V = -∑ V^{IJ} Λ^{JI}_{ji}|φ_i^I><φ_j^J|. But from E^V (Eq. 17) the correct potential is δE^V/δρ = -∑ V^{IJ} P^I ρ P^J = -∑ V^{IJ} Λ^{IJ}_{ij}|φ_i^I><φ_j^J|. Since Λ^{JI}_{ji} = (Λ^{IJ}_{ij})^*, Eq. (11) is the Hermitian conjugate of the correct operator for complex inter-site density-matrix elements. This error propagates to Eq. (5)/(21): the generator becomes S_{ij} = V Λ^*/D instead of V Λ/D, where D = U_I(1-2λ_i)-U_J(1-2λ_j). With the wrong S, [S,V_U]_{ij} = -V Λ^*/2, not -V Λ/2, so Eq. (20) fails unless Λ is real. Consequently the first-order energy satisfies E^U[S]-E^U = -V Re(Λ^2), while E^V = -V|Λ|^2. For imaginary Λ these differ in sign. The SI S6 projector-response cancellation (which relies on Eq. (20)) also breaks down. The central formula is therefore not correct for a general density matrix, although a simple conjugation repair restores the proof.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims that, to first order in V/U, the inter-site Hubbard correction DFT+U+V is exactly equivalent to a standard on-site DFT+U functional evaluated on a density-dependent, redefined set of Hubbard projectors hybridized with neighboring-site orbitals. The central result is an explicit rotation-like generator (Eq. (5)) and an energy/potential identity (Eqs. (22)-(23)). The authors also analyze the frozen-projector case, show that only a partial equivalence holds there, and illustrate the formalism on bulk NiO.","tokens_in":22534,"tokens_out":33323,"duration_ms":306075,"significance":"If the equivalence is correct, it gives the first formal interpretation of the widely used +V correction: +V does not introduce a new interaction beyond the +U occupation-curvature penalty, but rather changes the subspace on which that penalty is imposed, thereby making the on-site correction aware of covalency. This has direct implications for how V is computed and for the broader question of what the Hubbard projectors represent. The paper is clearly written, includes open data and reproducible scripts, and provides a numerical demonstration. However, the central algebraic derivation contains a sign/conjugation error that must be corrected before the claim is reliable.","major_comments":[{"comment":"The definition of the inter-site potential is incorrect. From Eq. (18), δE^V/δρ = -∑ V^{IJ} P^I ρ P^J, so the matrix element in the φ basis is -V^{IJ} Λ^{IJ}_{ij} = -V^{IJ} ⟨φ_i^I|ρ|φ_j^J⟩. Eq. (11) instead writes -V^{IJ} Λ^{JI}_{ji}, which is the complex conjugate. This is not a harmless reordering: for a complex density matrix the two operators differ. This error propagates into the generator: solving Eq. (20) with the correct V_V gives S_{ij}^{IJ} = V^{IJ} Λ^{IJ}_{ij}/D for the (I,J) block, not V^{IJ} Λ^{JI}_{ji}/D as in Eq. (5)/(21).","section":"Appendix, Eq. (11)"},{"comment":"The generator as printed is not anti-Hermitian. For real Λ, Eq. (5) gives S_{ij}^{IJ}=S_{ji}^{JI}, i.e. a symmetric matrix, whereas the derivation assumes an anti-Hermitian T so that e^T preserves orthonormality. With this symmetric S, [S,V_U] is antisymmetric (e.g. in a two-site/one-orbital model, [S,V_U]_{12}=-aD/2, [S,V_U]_{21}=+aD/2), so Eq. (20) cannot be satisfied with the symmetric V_V of Eq. (11). Consequently the first-order energy change Tr[V_U[ρ,S]] vanishes for real Λ, and Eq. (22) does not hold as stated. The proof can be repaired by choosing an antisymmetric S, e.g. S_{ij}^{IJ}=V^{IJ}Λ^{IJ}_{ij}/D for an ordered pair and the negative for the swapped pair, but this is a substantive correction, not a typo.","section":"Appendix, Eqs. (5), (20)"},{"comment":"The numerical validation uses NiO at the Γ point, where the inter-site density-matrix elements are real; it therefore cannot distinguish Λ^{IJ} from Λ^{JI}. Given that the paper's claim is general, the authors should either state explicitly that the equivalence holds only for real density matrices in the current derivation, or provide a test with complex inter-site elements (e.g. spin-orbit coupling or a non-Γ k-point) after correcting the generator. The conjugation issue is otherwise hidden by the test system.","section":"§4, Fig. 2 and SI S4"}],"minor_comments":[{"comment":"The abstract says 'to first order in V/U'; the expansion in fact also involves the occupations through the denominator D, and the authors themselves note the failure as λ_i→1/2. This caveat appears only later in the text; a sentence in the abstract or introduction would help set expectations.","section":"Abstract and §2"},{"comment":"The symbols Λ^{JI}_{ji} are easy to misread. After correcting the generator, define the ordered-pair convention explicitly (e.g. for I<J and I>J separately) so that the reader can verify anti-Hermiticity.","section":"Notation, Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The reported error in Eq. (11) is a transpose/conjugation mistake that invalidates the printed proof, but the underlying idea is sound and the fix is local. The numerical agreement in Fig. 2 suggests the implementation may already use an effectively corrected antisymmetric rotation; this should be checked. I recommend major revision rather than rejection, provided the authors correct the generator and either restrict or extend the claim appropriately."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: the paper's interpretation of DFT+U+V as DFT+U with density-dependent hybridized projectors is probably right, and the explicit generator is a nice piece of work. But the central derivation has a conjugation error for complex off-diagonal density-matrix elements. It's fixable, but it has to be fixed before the equivalence claim is safe.\n\nWhat's genuinely new: the explicit first-order generator S in Eq. (5), the proof that DFT+U+V reproduces DFT+U on rotated projectors to order V/U, and the argument that the +V correction is really an on-site curvature penalty applied to projectors that already contain covalency. The frozen-projector obstruction in S7/S8 is also new and useful, and the paper is honest about the limits: first-order validity, failure near lambda=1/2, ortho-atomic projectors only. The NiO illustration is nice and the data are archived.\n\nThe soft spot is not the disclosed limitations — it's an actual error. Eq. (11) defines the inter-site potential with Lambda^{JI}_{ji}, but Eq. (18), the functional derivative of E^V, gives Lambda^{IJ}_{ij}. For real density matrices these are the same; for complex ones they differ by conjugation. This propagates into the generator Eq. (5), and the identity Eq. (19) that connects Tr[rho V_V] to 2E^V only holds when Lambda is real. So the first-order equivalence as written fails for a complex inter-site density matrix. The stress-test note shows the fix is straightforward — conjugate the numerator in S — and once that's done the proof goes through. But the current text is not correct in the general case, and the numerical test at Gamma in NiO probably doesn't expose it because those matrix elements are real or effectively so.\n\nThat's a load-bearing flaw in the formal proof, though not in the physical message. If the authors repair the conjugation, the paper would be solid. I'd send it to review — the topic matters, the idea is important, and the rest of the algebra is careful — but the reviewers should insist on a corrected derivation and a check that the potential equivalence holds for complex Lambda. For my own work, I wouldn't cite this version; I'd wait for the revision.","headline":"Good physical idea and mostly clean proof, but the central algebra has a conjugation error for complex density matrices that needs fixing before the equivalence claim is safe.","tokens_in":23008,"tokens_out":8439,"would_cite":false,"duration_ms":77733,"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":"To first order in V/U, the inter-site +V correction of DFT+U+V is exactly an on-site DFT+U correction evaluated on a density-dependent set of Hubbard projectors hybridized with the orbitals of neighboring sites.","keywords":["DFT+U+V","inter-site Hubbard corrections","Hubbard projectors","self-interaction error","covalency","strongly correlated materials","ortho-atomic scheme"],"falsifier":"In a plane-wave PAW code (which uses non-orthogonal projectors), construct S^ from Eq. (5) and compare the DFT+U+V energy and potential against DFT+U with the rotated projectors over a range of V/U. If the difference scales as V/U rather than (V/U)^2, the ortho-atomic assumption is load-bearing. Alternatively, in an ortho-atomic code, drive an orbital occupation toward 1/2 and check whether the residual ceases to be O(V^2/U).","tokens_in":22071,"feed_emoji":"⚛️","tokens_out":8546,"duration_ms":82477,"temperature":0.7,"pith_summary":"This paper proves that the inter-site Hubbard correction +V, used in the DFT+U+V method, does not introduce a new kind of interaction: to first order in the ratio V/U, a DFT+U+V energy is exactly the standard DFT+U energy evaluated on a rotated set of Hubbard projectors that are hybridized with orbitals of neighboring sites. The rotation is generated by an explicit anti-Hermitian operator built from the inter-site density matrix and the Hubbard parameters, and the projectors are necessarily density-dependent because the generator is a functional of the density. The reason a reader should care is that this explains why +V works: it moves the on-site occupation-curvature penalty (the self-interaction correction) onto orbitals that already contain the covalent character of the bond, counteracting the over-localization and broken bonds that plain +U causes. For practitioners it also reshapes the theory: the choice of projectors, Hubbard parameters, and whether to include +V are not independent decisions.","feed_headline":"Inter-site Hubbard +V is +U on rotated projectors","feed_subtitle":"The inter-site correction adds no new physics—it shifts where the on-site penalty is measured, restoring covalency.","key_machinery":"The load-bearing object is the anti-Hermitian generator S^ (Eq. 5) that rotates the 'block-diagonalizing' projectors |φ^{Iσ}_i⟩—the natural orbitals that diagonalize each site's occupation matrix—into hybridized projectors e^{S^}|φ^{Iσ}_i⟩ that mix in the orbitals of V-coupled neighbors. The proof reduces the equivalence to the commutator identity [S^, V^U] = (1/2)V^V, which is solvable precisely because both the inter-site energy and potential are quadratic in the density, so 2E_V = Tr(ρ V^V). This single identity forces the form of S^, makes the projectors density-dependent, and yields the energy and potential equalities order by order.","core_discovery":"Central result: to first order in V/U, E^U[S^] = E^U + E^V + O(V^2/U), and the same for the potential. E^U[S^] is the DFT+U energy with projectors e^{S^}|φ^{Iσ}_i⟩, where S^ is the anti-Hermitian generator whose matrix elements mix site I and J through V^{IJ}, the inter-site occupancy Λ^{JIσ}_{ji}, and the denominator U^I(1-2λ^{Iσ}_i)-U^J(1-2λ^{Jσ}_j). Because S^ is built from the density matrix, the projectors are density-dependent. The corollary is that +V is the same self-interaction penalty as +U, but measured on projectors that carry the bond's covalent character—which is why +V restores covalency. With frozen projectors, the equality holds for the energy with e^{S^} and for the potenti","pith_inferences":["If the all-orders diagonalization of the corrective potential remains of +U type (the paper shows the first-order term is, while the second-order term is not), then extended Hubbard functionals beyond Dudarev-type +V might also reduce to some projector redefinition, unifying the extended-Hubbard family.","The divergence of the generator as an occupation approaches 1/2 delimits the practical range of the equivalence; comparing energies in systems with a nearly half-filled strongly hybridized orbital (e.g., mixed-valence oxides) would reveal where the first-order statement breaks.","The density-dependent projectors suggest a self-consistent scheme in which the rotation is updated during the SCF cycle; one testable prediction is that such a scheme exactly matches DFT+U+V energies and potentials while sidestepping the frozen-projector mismatch."],"forward_implications":["Any DFT+U+V calculation can be replaced, to first order in V/U, by a DFT+U calculation on hybridized density-dependent projectors, bringing the full DFT+U implementation ecosystem (response properties, extensions) to bear on +V problems.","The +V correction targets the same local occupation-curvature self-interaction as +U; it should be parametrized by linearizing the energy in the occupation of the hybridized projectors, a condition the usual off-diagonal response-matrix recipe does not satisfy.","With frozen projectors, no single rotation reproduces the DFT+U+V potential: the spectrum requires the e^{2S^} rotation while the energy requires e^{S^}, so frozen-projector DFT+U approximates DFT+U+V only up to the residual W^ term.","Applying +U to the hybridized projectors enhances covalency for less-than-half-filled antibonding orbitals (the minority-spin e_g states in NiO), directly counteracting the over-localization caused by +U on bare projectors.","Projector choice, U values, and the decision to include +V are coupled degrees of freedom, not independent choices; V is, at heart, a choice of projectors."],"fun_headline_variants":["+V is +U on density-rotated projectors","Inter-site +V equals on-site +U to first order","Hubbard +V is +U with covalent projectors","Density-dependent projectors unify +V and +U","The +V correction is +U on a rotated basis"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The proof assumes the Hubbard projectors are orthonormal across sites (the ortho-atomic scheme); with non-orthogonal projectors, inter-site overlap matrices enter the operator algebra and the simple generator form cannot be assumed to hold.","fun_headline_variants_meta":{"raw":{"variants":["+V is +U on density-rotated projectors","Inter-site +V equals on-site +U to first order","Hubbard +V is +U with covalent projectors","Density-dependent projectors unify +V and +U","The +V correction is +U on a rotated basis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000406,"raw_usage":{"total_tokens":1959,"prompt_tokens":769,"completion_tokens":1190,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":1119}},"tokens_in":513,"tokens_out":1190,"duration_ms":11057,"temperature":1.0,"reasoning_tokens":1119,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T16:10:50.379984+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a plane-wave PAW code (which uses non-orthogonal projectors), construct S^ from Eq. (5) and compare the DFT+U+V energy and potential against DFT+U with the rotated projectors over a range of V/U. If the difference scales as V/U rather than (V/U)^2, the ortho-atomic assumption is load-bearing. Alternatively, in an ortho-atomic code, drive an orbital occupation toward 1/2 and check whether the residual ceases to be O(V^2/U).","supporting_citations":[],"review_version":1}