pith. sign in

arxiv: 2605.16666 · v1 · pith:YXACUU3Fnew · submitted 2026-05-15 · ⚛️ physics.chem-ph · cond-mat.dis-nn· cond-mat.mtrl-sci· cond-mat.str-el· physics.comp-ph

Reducing the Complexity of Density-Matrix Functionals in a Real-Space-Decomposed DF+RDMF Scheme with the Adaptive Cluster Approximation

Konstantin Tamoev , Robert Schade , Thomas D. K\"uhne This is my paper

Pith reviewed 2026-05-19 20:36 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.dis-nncond-mat.mtrl-scicond-mat.str-elphysics.comp-ph
keywords reduced density-matrix functional theorydensity functional theoryadaptive cluster approximationcarbon suboxidebending potentialreal-space decompositionelectronic correlationembedding methods
0
0 comments X

The pith

A real-space split of the Coulomb interaction lets reduced-density-matrix corrections be applied only where needed, fixing the geometry of carbon suboxide.

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

The paper formulates a hybrid DF+RDMF scheme that partitions the Coulomb interaction locally in real space so the expensive reduced-density-matrix functional applies only to the strongly correlated region while a standard density functional handles the rest. The local functionals are then compressed by the adaptive cluster approximation, which rotates the bath subspace unitarily before truncation to keep the local interactions intact but shrink the number of states treated explicitly. On the bending potential of C3O2 this correction reverses the linear preference given by plain PBE and produces a bent minimum in line with spectroscopic inferences. If the partitioning and truncation preserve accuracy, the method supplies a practical embedding route for adding explicit correlation to large or inhomogeneous systems without full exponential cost everywhere.

Core claim

By partitioning the Coulomb interaction locally in real space and evaluating the RDMF correction only for the strongly correlated part, then compressing the local density-matrix functionals via the adaptive cluster approximation that performs a unitary rotation of the bath subspace before truncation, the scheme reduces the number of explicitly correlated bath states while preserving the local interaction and stabilizes a bent configuration of C3O2 in qualitative agreement with spectroscopy, unlike the linear structure favored by semilocal PBE.

What carries the argument

The adaptive cluster approximation, which unitarily rotates the bath subspace before truncation to preserve local interactions while reducing the number of explicitly correlated bath states.

Load-bearing premise

The Coulomb interaction can be partitioned locally in real space such that the RDMF correction needs to be evaluated only for the strongly correlated part without significant loss of overall accuracy.

What would settle it

A full-system RDMF calculation on C3O2 that yields a linear energy minimum while the decomposed DF+RDMF/ACA version yields a bent minimum would demonstrate that the local partitioning introduces unacceptable error.

Figures

Figures reproduced from arXiv: 2605.16666 by Konstantin Tamoev, Robert Schade, Thomas D. K\"uhne.

Figure 3
Figure 3. Figure 3: FIG. 3. Convergence of the local density-matrix correction [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Bending potentials of C [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

Reduced density-matrix functional theory (RDMFT) provides a variational route to electronic correlations beyond conventional density-functional approximations, but explicit evaluations of density-matrix functionals still scale exponentially with the number of active one-particle states. We formulate and assess a real-space-decomposed density-functional plus reduced-density-matrix-functional (DF+RDMF) scheme in which the Coulomb interaction is partitioned locally in real space and the RDMF correction is evaluated only for the strongly correlated part of the interaction. The resulting local density-matrix functionals are further compressed using the adaptive cluster approximation (ACA), which performs a unitary rotation of the bath subspace before truncation and therefore preserves the local interaction while reducing the number of explicitly correlated bath states. As a molecular test case, we consider the bending potential of carbon suboxide, C$_3$O$_2$. While semilocal PBE favors a linear molecule, the DF+RDMF/ACA correction stabilizes a bent configuration in qualitative agreement with the quasilinear behavior inferred from spectroscopy. The approach provides a systematic embedding hierarchy for combining density functionals with explicitly correlated density-matrix corrections in extended or spatially inhomogeneous systems.

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

2 major / 1 minor

Summary. The manuscript formulates a real-space-decomposed DF+RDMF scheme in which the Coulomb interaction is partitioned locally so that the RDMF correction is applied only to the strongly correlated subsystem; the local functionals are then compressed via the adaptive cluster approximation (ACA), which performs a unitary bath rotation before truncation. As a test case the bending potential of C3O2 is examined, where the DF+RDMF/ACA correction is reported to stabilize a bent minimum in qualitative agreement with the quasilinear geometry inferred from spectroscopy, in contrast to the linear preference of semilocal PBE.

Significance. If the local partitioning and ACA truncation preserve the essential correlation physics, the approach supplies a systematic embedding hierarchy that could combine the efficiency of density functionals with explicit two-body corrections for spatially inhomogeneous or extended systems. The ACA construction that preserves the local interaction while reducing the number of bath states is a technically attractive feature.

major comments (2)
  1. [abstract / scheme formulation] Abstract and scheme-formulation paragraph: the central assertion that the real-space partitioning isolates the strongly correlated part 'without significant loss of overall accuracy' is load-bearing for the claim that the bent stabilization is a faithful capture of correlation rather than an embedding artifact. No tests varying the partitioning threshold, spatial cutoff, or bath size are described, leaving open the possibility that delocalized contributions across the central carbon or oxygen atoms are inadvertently excluded.
  2. [molecular test case] Molecular test case: only qualitative agreement with the quasilinear behavior is reported; no quantitative energy differences (e.g., barrier height relative to linear geometry), error bars, or comparisons to reference methods such as CCSD(T) or experiment are provided. This makes it impossible to assess whether the ACA truncation introduces geometry-dependent errors that could mimic the observed stabilization.
minor comments (1)
  1. [scheme formulation] Notation for the partitioned interaction and the ACA bath rotation should be defined explicitly with an equation or diagram early in the manuscript to avoid ambiguity when the scheme is applied to other systems.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and the positive assessment of the potential of the DF+RDMF/ACA scheme. We address each major comment below and describe the revisions planned for the manuscript.

read point-by-point responses

  1. Referee: [abstract / scheme formulation] Abstract and scheme-formulation paragraph: the central assertion that the real-space partitioning isolates the strongly correlated part 'without significant loss of overall accuracy' is load-bearing for the claim that the bent stabilization is a faithful capture of correlation rather than an embedding artifact. No tests varying the partitioning threshold, spatial cutoff, or bath size are described, leaving open the possibility that delocalized contributions across the central carbon or oxygen atoms are inadvertently excluded.

    Authors: We agree that explicit tests of the partitioning and truncation parameters are needed to substantiate the claim. In the revised manuscript we will add a dedicated subsection (or appendix) reporting calculations with varied partitioning thresholds, spatial cutoffs, and ACA bath sizes. These tests will show that the bent minimum remains stable and that the energy ordering is insensitive to reasonable changes in these parameters, thereby confirming that the real-space decomposition does not inadvertently exclude essential delocalized contributions. revision: yes

  2. Referee: [molecular test case] Molecular test case: only qualitative agreement with the quasilinear behavior is reported; no quantitative energy differences (e.g., barrier height relative to linear geometry), error bars, or comparisons to reference methods such as CCSD(T) or experiment are provided. This makes it impossible to assess whether the ACA truncation introduces geometry-dependent errors that could mimic the observed stabilization.

    Authors: The present manuscript focuses on the qualitative contrast with PBE and the agreement with the spectroscopically inferred quasilinear geometry. To strengthen the assessment, the revised results section will include quantitative energy differences (linear–bent gap and approximate barrier height), convergence data with respect to ACA bath size (providing an estimate of truncation error), and direct comparison with available experimental structural parameters. Where computationally feasible we will also reference CCSD(T) results for the relevant energy differences. These additions will allow a clearer evaluation of possible geometry-dependent truncation artifacts. revision: yes

Circularity Check

0 steps flagged

No significant circularity; DF+RDMF/ACA is an independent approximation hierarchy

full rationale

The paper presents a systematic embedding scheme that partitions the Coulomb interaction in real space, applies RDMF corrections only to the strongly correlated local part, and compresses the bath via ACA unitary rotation plus truncation. The bending potential result for C3O2 follows directly from applying this hierarchy to the molecule (PBE baseline plus local correction), without any equation reducing a claimed prediction to a fitted parameter or self-citation by construction. The central claim remains an independent numerical outcome of the stated approximations and is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review limits extraction; the scheme rests on domain assumptions about local interaction partitioning and ACA preservation of local terms, with no explicit free parameters or invented entities named.

axioms (2)
  • domain assumption The Coulomb interaction can be partitioned locally in real space allowing RDMF correction only on the strongly correlated part.
    Invoked in the formulation of the DF+RDMF scheme.
  • domain assumption Unitary rotation of the bath subspace before truncation in ACA preserves the local interaction.
    Central justification for the compression step.

pith-pipeline@v0.9.0 · 5761 in / 1329 out tokens · 55784 ms · 2026-05-19T20:36:12.856199+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.

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages

  1. [1]

    Semilocal PBE favors the linear structure, whereas the real-space-decomposed DF+RDMF/ACA selected-CI result stabilizes a bent geometry. With this partition the interaction contribution is ap- proximated as F W [ρ(1)]≈F WC1 DF [ρ(1)] +F WC2 DF [ρ(1)] +F WC3 RDMF[ρ(1)] +F WO1 DF [ρ(1)] +F WO2 DF [ρ(1)],(33) where the explicit RDMF correction is centered onC...

    work page

  2. [2]

    J. A. Pople, Rev. Mod. Phys.71, 1267 (1999)

    work page 1999

  3. [3]

    Kohn, Rev

    W. Kohn, Rev. Mod. Phys.71, 1253 (1999)

    work page 1999

  4. [4]

    Pernal and K

    K. Pernal and K. J. H. Giesbertz, inDensity-Functional Methods for Excited States(Springer International Pub- lishing, Cham, 2015) pp. 125–183

    work page 2015

  5. [5]

    T. L. Gilbert, Phys. Rev. B12, 2111 (1975)

    work page 1975

  6. [6]

    A. M. K. M¨ uller, Phys. Lett. A105, 446 (1984)

    work page 1984

  7. [7]

    Goedecker and C

    S. Goedecker and C. J. Umrigar, Phys. Rev. Lett.81, 866 (1998)

    work page 1998

  8. [8]

    Schade, C

    R. Schade, C. Bauer, K. Tamoev, L. Mazur, C. Plessl, and T. D. K¨ uhne, Phys. Rev. Research4, 033160 (2022)

    work page 2022

  9. [9]

    Gritsenko, K

    O. Gritsenko, K. Pernal, and E. J. Baerends, J. Chem. Phys.122, 204102 (2005)

    work page 2005

  10. [10]

    Piris, Phys

    M. Piris, Phys. Rev. Lett.127, 233001 (2021)

    work page 2021

  11. [11]

    N. N. Lathiotakis and M. A. L. Marques, J. Chem. Phys. 128, 184103 (2008)

    work page 2008

  12. [12]

    N. N. Lathiotakis, S. Sharma, J. K. Dewhurst, F. G. Eich, M. A. L. Marques, and E. K. U. Gross, Phys. Rev. A79, 7 040501 (2009)

    work page 2009

  13. [13]

    Pernal, K

    K. Pernal, K. J. H. Giesbertz, O. Gritsenko, and E. J. Baerends, J. Chem. Phys.127, 214101 (2007)

    work page 2007

  14. [14]

    K. J. H. Giesbertz, O. V. Gritsenko, and E. J. Baerends, Phys. Rev. Lett.105, 013002 (2010)

    work page 2010

  15. [15]

    Sharma, J

    S. Sharma, J. K. Dewhurst, N. N. Lathiotakis, and E. K. U. Gross, Phys. Rev. B78, 201103 (2008)

    work page 2008

  16. [16]

    Sharma, J

    S. Sharma, J. K. Dewhurst, S. Shallcross, and E. K. U. Gross, Phys. Rev. Lett.110, 116403 (2013)

    work page 2013

  17. [17]

    P. E. Bl¨ ochl, C. F. J. Walther, and T. Pruschke, Phys. Rev. B84, 205101 (2011)

    work page 2011

  18. [18]

    Schade,New methods for the ab-initio simula- tion of correlated systems, Ph.D

    R. Schade,New methods for the ab-initio simula- tion of correlated systems, Ph.D. thesis, Georg-August- Universit¨ at G¨ ottingen (2019)

    work page 2019

  19. [19]

    Schade and P

    R. Schade and P. E. Bl¨ ochl, Phys. Rev. B97, 245131 (2018)

    work page 2018

  20. [20]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

    work page 1996

  21. [21]

    P. E. Bl¨ ochl, R. Schade, L. Allen-Rump, S. Rajpuro- hit, A. Rathnakaran, K. Tamoev, M. Lokamani, and T. D. K¨ uhne, The CP-PAW code package for first- principles calculations from a user’s perspective (2026), arXiv:2601.12004 [cond-mat.mtrl-sci]

    work page arXiv 2026

  22. [22]

    A. J. Coleman, Rev. Mod. Phys.35, 668 (1963)

    work page 1963

  23. [23]

    Schilling, D

    C. Schilling, D. Gross, and M. Christandl, Phys. Rev. Lett.110, 040404 (2013)

    work page 2013

  24. [24]

    K. J. H. Giesbertz and M. Ruggenthaler, Phys. Rep.806, 1 (2019)

    work page 2019

  25. [25]

    D. A. Mazziotti, Phys. Rev. Lett.108, 263002 (2012)

    work page 2012

  26. [26]

    Levy, Proc

    M. Levy, Proc. Natl. Acad. Sci. U.S.A.76, 6062 (1979)

    work page 1979

  27. [27]

    S. M. Valone, J. Chem. Phys.73, 1344 (1980)

    work page 1980

  28. [28]

    Helbig,Orbital functionals in density-matrix- and current-density-functional theory, Ph.D

    N. Helbig,Orbital functionals in density-matrix- and current-density-functional theory, Ph.D. thesis, Freie Universit¨ at Berlin (2006)

    work page 2006

  29. [29]

    P. E. Bl¨ ochl, T. Pruschke, and M. Potthoff, Phys. Rev. B88, 205139 (2013)

    work page 2013

  30. [30]

    Baldsiefen, A

    T. Baldsiefen, A. Cangi, and E. K. U. Gross, Phys. Rev. A92, 052514 (2015)

    work page 2015

  31. [31]

    Diels and B

    O. Diels and B. Wolf, Ber. Dtsch. Chem. Ges.39, 689 (1906)

    work page 1906

  32. [32]

    Bereza and H

    S. Bereza and H. Staudinger, Ber. Dtsch. Chem. Ges.41, 4461 (1908)

    work page 1908

  33. [33]

    Staudinger, Ber

    H. Staudinger, Ber. Dtsch. Chem. Ges.53, 1073 (1920)

    work page 1920

  34. [34]

    L´ opez-Salas, J

    N. L´ opez-Salas, J. Kossmann, and M. Antonietti, Acc. Mater. Res.1, 117 (2020)

    work page 2020

  35. [35]

    Odziomek, P

    M. Odziomek, P. Giusto, J. Kossmann, N. V. Tarak- ina, J. Heske, S. M. Rivadeneira, W. Keil, C. Schmidt, S. Mazzanti, O. Savateev, L. Perdig´ on-Toro, D. Neher, T. D. K¨ uhne, M. Antonietti, and N. L´ opez-Salas, Adv. Mater.34, 2206405 (2022)

    work page 2022

  36. [36]

    R. B. Woodward and R. Hoffmann, J. Am. Chem. Soc. 87, 395 (1965)

    work page 1965

  37. [37]

    M. H. Garner, R. Hoffmann, S. Rettrup, and G. C. Solomon, ACS Cent. Sci.4, 688 (2018)

    work page 2018

  38. [38]

    L. A. Carreira, R. O. Carter, J. R. Durig, R. C. Lord, and C. C. Milionis, J. Chem. Phys.59, 1028 (1973)

    work page 1973

  39. [39]

    W. H. Weber and G. W. Ford, J. Mol. Spectrosc.63, 445 (1976)

    work page 1976

  40. [40]

    Fusina and I

    L. Fusina and I. M. Mills, J. Mol. Spectrosc.79, 123 (1980)

    work page 1980

  41. [41]

    P. R. Bunker, J. Mol. Spectrosc.80, 422 (1980)

    work page 1980

  42. [42]

    R. L. Lozes, J. R. Sabin, and J. Oddershede, J. Mol. Spectrosc.86, 357 (1981)

    work page 1981

  43. [43]

    Jensen and J

    P. Jensen and J. W. C. Johns, J. Mol. Spectrosc.118, 248 (1986)

    work page 1986

  44. [44]

    Masiello, S

    T. Masiello, S. L. Coy, B. C. West, B. C. Unruh, N. D. Zobov, O. L. Polyansky, J. Tennyson, and R. W. Field, J. Chem. Phys.123, 234308 (2005)

    work page 2005

  45. [45]

    Car and M

    R. Car and M. Parrinello, Phys. Rev. Lett.55, 2471 (1985)

    work page 1985

  46. [46]

    T. D. K¨ uhne, M. Krack, F. R. Mohamed, and M. Par- rinello, Phys. Rev. Lett.98, 066401 (2007)

    work page 2007