Radial Distribution Function in a Two Dimensional Core-Shoulder Particle System

Andrew J Archer; Gerhard Kahl; Michael Wassermair; Roland Roth

arxiv: 2603.24537 · v2 · submitted 2026-03-25 · ❄️ cond-mat.soft · cond-mat.stat-mech

Radial Distribution Function in a Two Dimensional Core-Shoulder Particle System

Michael Wassermair , Gerhard Kahl , Andrew J Archer , Roland Roth This is my paper

Pith reviewed 2026-05-15 00:18 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech

keywords radial distribution functiondensity functional theorycore-shoulder potentialOrnstein-Zernike equationtest-particle methodtwo dimensionspair correlations

0 comments

The pith

In a 2D core-shoulder system, the Ornstein-Zernike route to the radial distribution function can outperform the test-particle route.

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

The paper investigates two methods for obtaining the radial distribution function from classical density functional theory in a two-dimensional system of particles with core-shoulder repulsions. One method inserts a test particle and minimizes the grand potential functional to find the surrounding density profile. The other extracts the pair direct correlation function from the second derivative of the excess free energy and solves the Ornstein-Zernike equation. Conventional wisdom holds that the test-particle method should be more accurate for approximate functionals since it uses only first derivatives. The calculations demonstrate that this is not the case here, with the Ornstein-Zernike results sometimes agreeing better with simulation data.

Core claim

For the two-dimensional core-shoulder particle system studied, the radial distribution function g(r) computed via the test-particle route in density functional theory is not always more accurate than the version obtained from the Ornstein-Zernike equation using the direct correlation function.

What carries the argument

The two alternative routes to the radial distribution function g(r) within classical density functional theory: the test-particle route and the Ornstein-Zernike route based on the pair direct correlation function.

If this is right

Accuracy of approximate DFT functionals for pair structure can reverse between the two routes in systems with shoulder potentials.
Both routes should be evaluated when applying DFT to compute g(r) in two-dimensional soft-matter models.
Core-shoulder interactions in 2D may expose limitations in the expected superiority of single-derivative methods.
The choice of route affects how well DFT matches exact g(r) depending on the state point.

Where Pith is reading between the lines

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

Similar reversals might occur in other low-dimensional systems with competing length scales in the potential.
Future functional approximations could be designed to balance performance across both routes rather than optimizing for one.
Testing in three dimensions would clarify if the two-dimensional geometry is essential to the observation.

Load-bearing premise

The specific choice of core-shoulder potential and the restriction to two dimensions are representative enough that the observed reversal challenges the general preference for the test-particle route.

What would settle it

Exact simulation data or results from a more accurate functional for the same 2D core-shoulder system where the test-particle g(r) is consistently closer to the true values than the Ornstein-Zernike g(r) would contradict the paper's finding.

read the original abstract

An important quantity in liquid state theory is the radial distribution function $g(r)$. It can be calculated within the framework of classical density functional theory in two very distinct ways. In the test-particle route, one fixes a single fluid particle, turning it into an external potential in which the inhomogeneous structure of the fluid is calculated by minimising the functional. The second route to $g(r)$ in density functional theory employs the Ornstein-Zernike equation and the pair direct correlation function, that can be obtained from the second functional derivatives of the excess free energy functional. Since typically an approximate excess free energy functional is employed, one generally expects that the test-particle route, which requires only one functional derivative, to be more accurate than the Ornstein-Zernike route. Here we study a two dimensional core-shoulder particle system and present results that challenge this expectation. Our results show that in this system test-particle results for $g(r)$ are not always better than results obtained via the Ornstein-Zernike route.

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.

Desk Editor's Note private letter to a colleague

The abstract claims a 2D core-shoulder system where the Ornstein-Zernike route beats the test-particle route for g(r), flipping the usual DFT expectation, but without any numbers or details this stays unverified.

read the letter

The paper's main finding is that in this two-dimensional core-shoulder particle system the radial distribution function from the Ornstein-Zernike route is sometimes more accurate than the one from the test-particle route. This runs counter to the standard expectation that the test-particle method should win when an approximate excess free-energy functional is in use, since it only requires a single functional derivative rather than the second derivatives needed for the direct correlation function in the OZ equation. They reach this by applying the usual DFT routes to the specific interaction and reporting the reversal in accuracy ordering. That is the concrete new piece: a case where the hierarchy does not hold for this model in 2D. The abstract sets the comparison out plainly and identifies the system clearly, so the claim itself is easy to understand. The limitation is obvious from the start. Only the abstract is available, so there are no numerical values, no specification of the excess free-energy approximation, no density or temperature ranges, no error bars, and no description of how accuracy was measured between the two routes. Without those, it is impossible to judge whether the reported reversal is robust or tied to a narrow set of parameters. The work sits squarely in classical liquid-state theory and soft-matter DFT. A reader who already knows the test-particle versus OZ comparison would see the potential value as a possible exception, but anyone outside that niche would get little from it. I would send the paper to peer review. The claim is specific and falsifiable, so referees can check the actual calculations and decide whether the evidence supports the reversal. If the numerics hold up, it supplies a useful reference point; if they do not, the paper can be corrected or rejected on that basis.

Referee Report

1 major / 0 minor

Summary. The manuscript examines the radial distribution function g(r) for a two-dimensional core-shoulder particle system within classical density functional theory. It contrasts the test-particle route (minimizing the grand-potential functional with one particle fixed as an external potential) against the Ornstein-Zernike route (using the pair direct correlation function obtained from the second functional derivative of an approximate excess free-energy functional). The central claim, based on numerical results for this specific system, is that the test-particle route does not always yield more accurate g(r) than the Ornstein-Zernike route, contrary to the usual expectation for approximate functionals.

Significance. If substantiated by the full calculations, the result supplies a concrete counter-example to the prevailing heuristic that the test-particle route is systematically superior when an approximate excess free-energy functional is employed. This could inform method selection in liquid-state theory for soft-core potentials in two dimensions. The work is grounded in a direct numerical comparison rather than additional fitted parameters, which strengthens its potential utility as a diagnostic case.

major comments (1)

[Abstract] Abstract: the claim that test-particle results 'are not always better' rests on unspecified numerical comparisons. No information is given on the approximate excess free-energy functional, the core-shoulder potential parameters, the density range, or quantitative error measures (e.g., integrated deviations from simulation benchmarks or direct comparison of the two routes), preventing verification that the reported counter-example is load-bearing rather than an artifact of a particular approximation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the major comment below.

read point-by-point responses

Referee: [Abstract] Abstract: the claim that test-particle results 'are not always better' rests on unspecified numerical comparisons. No information is given on the approximate excess free-energy functional, the core-shoulder potential parameters, the density range, or quantitative error measures (e.g., integrated deviations from simulation benchmarks or direct comparison of the two routes), preventing verification that the reported counter-example is load-bearing rather than an artifact of a particular approximation.

Authors: We agree that the abstract would benefit from greater specificity to allow readers to immediately assess the strength of the counter-example. The body of the manuscript contains the full numerical comparisons, including the approximate excess free-energy functional employed, the core-shoulder potential parameters, the density range, and quantitative error measures (integrated deviations from simulation benchmarks) for both routes. To address the referee's concern, we will revise the abstract to incorporate concise statements of these details and the direct comparison of the two routes. This change will make the central claim more transparent without altering the manuscript's conclusions. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The abstract presents a direct numerical comparison of two routes to g(r) (test-particle vs. Ornstein-Zernike) using the same approximate excess free-energy functional in one specific 2D core-shoulder system. No derivation, ansatz, fitted parameter, or self-citation chain is described that reduces the central claim to its own inputs by construction. The result is a falsifiable computational finding rather than an algebraic identity or renamed input.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard liquid-state assumptions; no new free parameters, invented entities, or ad-hoc axioms are introduced in the abstract.

axioms (2)

domain assumption An approximate excess free-energy functional is employed
Standard in classical DFT calculations of structure
standard math The Ornstein-Zernike equation relates the pair correlation function to the direct correlation function
Fundamental relation in liquid theory invoked for the second route

pith-pipeline@v0.9.0 · 5457 in / 1220 out tokens · 57898 ms · 2026-05-15T00:18:13.821891+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

test-particle route, which requires only one functional derivative, to be more accurate than the Ornstein-Zernike route

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.