🔬 condensed matter

Uniform distributions in nonuniform systems: Wall potentials generating constant density profiles in classical density functional theory

This paper solves the inverse problem in classical density functional theory by deriving explicit analytical and numerical expressions for wall potentials that generate perfectly flat equilibrium density profiles in planar, spherical, and cylindrical geometries for both hard-sphere and Lennard-Jones fluids within Rosenfeld's fundamental measure theory.

Original authors: Jiří Janek, Alexandr Malijevský

Published 2026-07-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jiří Janek, Alexandr Malijevský

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to arrange a crowd of people (fluid molecules) in a room next to a wall. Usually, when you put a wall in a room, the people naturally crowd up against it, creating layers, gaps, or specific patterns. This happens because the people bump into each other (repulsion) and sometimes stick together (attraction).

In the world of physics, this is described by something called Density Functional Theory (DFT). Usually, scientists use this theory to ask: "If I have a specific wall, how will the people arrange themselves?"

This paper asks the reverse question: "What kind of special wall do I need to build so that the people arrange themselves in a perfectly flat, uniform line, with no crowding or gaps at all?"

Here is a breakdown of how they solved this puzzle, using simple analogies:

1. The Goal: A Perfectly Flat Crowd

Normally, if you put a wall next to a fluid, the fluid creates "ripples" or "layers" near the wall, like water lapping against a dock. The authors wanted to find a "magic wall" that cancels out these ripples entirely, leaving the fluid density perfectly flat and uniform right up to the wall.

2. The Tools: Measuring the "Push" and "Pull"

To figure out what this magic wall looks like, the authors used a sophisticated mathematical toolkit called Rosenfeld's Fundamental Measure Theory.

Think of this toolkit as a way to measure exactly how much the fluid molecules "push" against each other (because they are hard spheres that can't overlap) and how much they "pull" on each other (because of weak attractive forces, like a gentle magnet).
The authors calculated exactly how much "extra push" or "pull" the fluid creates naturally near a wall.
Then, they designed a wall potential (an invisible force field) that does the exact opposite. If the fluid pushes too hard, the wall pushes back just enough to cancel it out. If the fluid pulls together, the wall pulls apart just enough to flatten the line.

3. The Three Shapes: Flat, Round, and Tube

The researchers tested this idea in three different room shapes:

Flat Walls (Planar): Imagine a long, straight hallway. This was the easiest case. They found a neat, simple formula for the magic wall. It turned out to be a short-range force that gets stronger as the crowd gets denser.
Round Walls (Spherical): Imagine the fluid is surrounding a giant ball. Here, the curvature of the ball changes how the fluid behaves. The math became much more complicated (like a tangled knot of equations), but they still found an exact solution. They discovered that the "magic force" needed changes depending on how big the ball is.
Tube Walls (Cylindrical): Imagine the fluid is inside a long pipe. This was the hardest case. The math involved complex shapes called "elliptic integrals" that couldn't be written as a simple formula. So, they used powerful computers to calculate the answer numerically.

4. The Results: Curvature Matters

The most interesting finding is how the shape of the wall changes the "magic force" needed:

Hard-Sphere Fluid (Just Bumping): If the molecules only bump into each other (no sticking), the shape of the wall doesn't change the required force very much. A flat wall and a round ball need almost the same "cancelling" force.
Attractive Fluid (Bumping + Sticking): If the molecules also stick to each other, the shape matters a lot. For a round ball, the "cancelling" force needs to be weaker and shaped differently than for a flat wall. The curvature of the wall fights against the natural tendency of the fluid to clump together.

5. The Proof: Did it Work?

To prove their math wasn't just theory, they took their calculated "magic walls" and ran computer simulations.

The Test: They told the computer, "Here is the wall we designed. Now, show us the fluid."
The Result: The fluid arranged itself in a perfectly flat line, exactly as predicted. The tiny errors seen were just due to the computer's grid size, not a flaw in the theory.

Summary

In short, this paper is like a master architect who figured out exactly how to build a wall so that a crowd of people standing next to it would stand in a perfectly straight, uniform line, ignoring all their natural instincts to crowd or spread out.

They showed that:

It is possible to calculate this "perfect wall" explicitly.
The math is simple for flat walls, complex for round balls, and requires computers for tubes.
The shape of the wall (curvature) changes the force needed to keep the fluid flat, especially if the fluid molecules like to stick together.

This work provides a "recipe" for creating perfectly uniform fluids in specific geometries, which helps scientists test their computer models and understand how walls and fluids interact at a microscopic level.

Technical Summary: Uniform Distributions in Nonuniform Systems

Problem Statement
This work addresses the inverse problem within Classical Density Functional Theory (DFT) for inhomogeneous fluids. While standard DFT applications solve the "direct problem"—determining the equilibrium density profile $\rho(r)$ for a prescribed external potential $V(r)$ —this study seeks the converse: constructing the specific external wall potential $V(r)$ required to generate a perfectly flat, constant equilibrium density profile $\rho(r) = \rho_b$ in the fluid region adjacent to a substrate.

The motivation stems from the observation that realistic substrates typically induce packing oscillations, depletion layers, or adsorption. The authors aim to identify "structure-cancelling" wall fields that explicitly suppress these inhomogeneities. This problem is linked to interfacial theories of wetting and drying, where flat density profiles arise under special matching conditions between wall-fluid and fluid-fluid interactions. The study investigates how the required compensating field depends on substrate geometry (planar, spherical, and cylindrical) and fluid interactions.

Methodology
The authors employ Rosenfeld's Fundamental Measure Theory (FMT) to describe the hard-sphere reference system and a mean-field treatment for attractive interactions. Two one-component fluid models are analyzed:

A purely repulsive hard-sphere fluid.
A fluid with hard-sphere repulsion plus a truncated Lennard-Jones-like attractive tail.

The study considers three fundamental geometries:

Planar: A flat wall at $z=R$ .
Spherical: A solid sphere of radius $\varrho$ .
Cylindrical: An infinitely long solid cylinder of radius $\varrho$ .

The core methodology involves substituting the prescribed constant density profile into the Euler-Lagrange equation derived from the grand-potential functional. The external potential is recovered via:
$\beta V(\xi) = \beta\mu + c^{(1)}(\xi) - \ln[\Lambda^3 \rho(\xi)]$
where $c^{(1)}$ is the one-body direct correlation function. The authors derive explicit expressions for the weighted densities and $c^{(1)}$ for each geometry.

Analytical Treatment: For planar and spherical geometries, the convolution integrals of FMT are reduced to one-dimensional forms that allow for closed analytical solutions.
Numerical Treatment: For cylindrical geometry, the reduction of FMT kernels leads to complete elliptic integrals, precluding a simple closed form; thus, this case is evaluated via numerical quadrature.

Key Results

Planar Geometry: The inverse problem yields explicit, compact analytical expressions for the wall potential. For hard spheres, the potential is purely repulsive and short-ranged, with an amplitude increasing with packing fraction $\eta$ . When attractive interactions are included, the potential develops an attractive component, reflecting the competition between packing and cohesion. The contact value of the potential satisfies the exact sum rule $\beta p_{hs}/\rho_b = 1 + \beta V_{hs}^{pl}(R)$ .
Spherical Geometry: Analytical solutions exist but are substantially more complex, depending on both $\eta$ and the substrate radius $\varrho$ . The results show that curvature has a weak effect on the purely repulsive hard-sphere potential but a more pronounced effect on fluids with attractive interactions, where the attractive well becomes shallower as the sphere radius decreases.
Cylindrical Geometry: The wall potentials were determined numerically. The results demonstrate that cylindrical potentials lie intermediate between the planar and spherical limits for the same radius, reflecting the intermediate curvature of the cylinder.
Validation: Independent DFT calculations using the derived potentials confirm that the prescribed flat density profiles are recovered within numerical accuracy (relative deviations of order $10^{-4}$ ), validating the analytical and numerical constructions.

Significance and Claims
The paper claims to provide an explicit microscopic realization of structure-cancelling wall fields, moving beyond formal or purely numerical inversions. The primary contributions are:

Explicit Analytical Construction: The derivation of closed-form expressions for wall potentials in planar and spherical geometries that cancel interfacial structure.
Geometric Dependence: A clear demonstration of how substrate curvature modifies the compensating field, particularly highlighting the stronger curvature dependence in systems with attractive interactions compared to purely repulsive ones.
Reference Material: The work provides a compact collection of formulae for weighted densities and one-body direct correlation functions in the three fundamental geometries, serving as useful reference expressions for DFT implementations.
Theoretical Context: The study establishes that flat-profile conditions, often treated as asymptotic or phenomenological matching conditions in wetting theories, can be realized explicitly within a microscopic nonlocal density functional.

The authors note that the range of the compensating potential is finite, tied to the range of the FMT kernels and the truncated attractive interaction. They acknowledge that beyond mean-field theory, critical fluctuations would introduce long-ranged contributions, which are outside the scope of the current finite-ranged treatment. The work suggests future extensions to multicomponent fluids, systematic curvature expansions, and the application of this inverse viewpoint to dynamical and rheological contexts.