🔬 condensed matter

Thermodynamics of hard-sphere fluids in polydisperse random porous media: Extended scaled particle theory

This paper extends scaled particle theory to provide a fully analytical description of the thermodynamic properties of hard-sphere fluids confined in size-polydisperse random porous media, demonstrating excellent agreement with Monte Carlo simulations and integral-equation theories across a broad range of matrix packing fractions and polydispersity degrees.

Original authors: T. Hvozd, M. Hvozd, M. Holovko

Published 2026-03-24

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: T. Hvozd, M. Hvozd, M. Holovko

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 organize a massive party in a room that is already crowded with furniture. But here's the twist: the furniture isn't just chairs and tables; it's a chaotic mix of giant beanbags, tiny stools, and everything in between, all frozen in place. Now, you want to invite guests (the fluid) to move around in the empty spaces left by this furniture.

This paper is about figuring out exactly how crowded the room feels to your guests and how much "effort" (energy) it takes to squeeze them into the space, depending on how messy and varied the furniture is.

Here is the breakdown of the science, translated into everyday concepts:

The Setting: The "Frozen Chaos" Room

In the world of physics, scientists often study Hard-Sphere Fluids. Think of these as perfect, bouncy balls that can't overlap. They are the simplest model for liquids, gases, and even complex things like proteins in your body.

Usually, scientists study these balls in an empty room. But in the real world (like in gels, rocks, or biological cells), these balls are trapped inside porous media—materials full of holes and obstacles.

The problem this paper tackles is that most previous theories assumed the "furniture" (the obstacles) was all the same size. But in reality, the obstacles are polydisperse. This is a fancy word for "mixed sizes." Some obstacles are huge boulders; others are tiny pebbles.

The Problem: The Old Maps Were Wrong

Scientists have a tool called Scaled Particle Theory (SPT). Imagine SPT as a map that predicts how hard it is to squeeze a new guest into the room based on the size of the furniture.

The Old Map: Previous versions of this map worked great when the furniture was all the same size or when the room wasn't too crowded.
The Glitch: When the room got very crowded (high packing fraction) and the furniture was a wild mix of sizes, the old maps started to break. They predicted things that didn't match reality, like saying a guest could fit in a space that was actually too small.

The Solution: A New, Smarter Map

The authors of this paper, T. Hvozd, M. Hvozd, and M. Holovko, decided to fix the map. They developed a new, more sophisticated version of the theory (which they call SPT2b3 and a new improved version, SPT2b3).

Here is how they did it, using an analogy:

Measuring the "Empty Space": They realized that just knowing the total volume of the room isn't enough. You need to know the shape of the empty spaces. If you have a room filled with giant boulders, the empty space is huge but hard to reach. If you have a room filled with tiny pebbles, the empty space is small but everywhere.
The "Probe" Test: They imagined sending in a "probe" guest of different sizes to test the room.
- A tiny probe can wiggle into the smallest cracks between the pebbles.
- A giant probe can only fit in the big gaps between the boulders.
The New Formula: By mathematically tracking how these probes of different sizes fit into the mixed-size furniture, they created a new equation. This equation accounts for the fact that mixing big and small obstacles creates a more efficient "packing" than having them all the same size.

The Results: Why It Matters

The authors tested their new map against Monte Carlo simulations (which are basically super-computer games where they randomly move millions of balls to see what happens).

The Verdict: Their new "SPT2b3**" map was incredibly accurate. It matched the computer simulations perfectly, even when the room was very crowded and the furniture was very mixed up.
The "Sweet Spot": They found that when you have a mix of sizes, the fluid behaves differently than if everything was uniform. The new theory captures this nuance perfectly.

The Real-World Impact

Why should you care about bouncy balls in a crowded room? Because this isn't just about balls. This math applies to:

Medicine: Understanding how proteins move and clump together inside your cells, which are packed with other molecules (a "crowded" environment).
Materials Science: Designing better filters, batteries, or catalysts where chemicals need to flow through tiny, messy pores.
Soft Matter: Creating better gels, paints, or colloids where particle size varies.

The Bottom Line

Think of this paper as the difference between using a generic "Do Not Enter" sign on a door versus having a detailed, 3D blueprint of the hallway behind it. The authors took a complex, messy situation (mixed sizes in a crowded space) and built a precise, analytical tool to predict exactly how things will behave. They didn't just fix the math; they gave scientists a reliable way to understand how fluids behave in the messy, crowded, real world.

1. Problem Statement

The study addresses the thermodynamic description of Hard-Sphere (HS) fluids confined within size-polydisperse random porous media. While the thermodynamics of fluids in monodisperse porous matrices have been extensively studied using analytical methods like Scaled Particle Theory (SPT) and integral equations, there is a significant gap in the theoretical understanding of systems where the porous matrix itself is polydisperse (i.e., composed of obstacles with a distribution of sizes).

Existing approaches for polydisperse matrices often rely on:

Numerical calculations or Monte Carlo (MC) simulations, which lack analytical transparency.
Assumptions of narrow size distributions.
Limitations in accuracy at high matrix packing fractions ( $\eta_0$ ).

The authors aim to develop a fully analytical framework that simultaneously accounts for the size-polydispersity of the matrix obstacles and the packing fraction of the fluid, providing accurate predictions for chemical potentials and pressures across a broad range of conditions.

2. Methodology

The authors extend the Scaled Particle Theory (SPT) framework, specifically building upon the SPT2 approach previously developed for monodisperse matrices.

Theoretical Framework

System Definition: The system consists of a monodisperse HS fluid (annealed) confined in a quenched, size-polydisperse HS matrix.
Core Concept: SPT calculates thermodynamic properties by determining the reversible work required to insert a scaled test particle of radius $R_s$ into the fluid/matrix system.
Key Parameters:
- Geometric Porosity ( $\phi_0$ ): The fraction of free volume not occupied by the matrix.
- Thermodynamic Porosity ( $\phi_1$ ): Related to the excess chemical potential of a fluid particle at infinite dilution in the matrix.
- Polydispersity: Handled by replacing discrete summations over matrix species with integrals over a continuous size distribution function $f(\sigma_0)$ , utilizing moments ( $m_l$ ) of the distribution.

Derivation of Approximations

The authors derive a hierarchy of analytical approximations to resolve divergences found in earlier SPT formulations at high densities:

SPT2a: A baseline approximation replacing $\phi_1$ with $\phi_0$ to avoid divergence, but known to underestimate porous media effects.
SPT2b: Includes a logarithmic divergence at $\eta_1 = \phi_1$ .
SPT2b1: Corrects SPT2b by expanding logarithmic terms to remove the divergence.
SPT2b3:* Introduces a third porosity parameter ( $\phi^*$ ) to describe the maximum allowable packing fraction, providing accuracy across low to medium densities.
SPT2b3 (New Contribution): The authors identify that SPT2b3* exhibits deviations at high matrix packing fractions ( $\eta_0 \geq 0.25$ ). They propose a new approximation, SPT2b3,** by neglecting the last term in the expansion of the logarithmic term used in SPT2b3*. This modification significantly improves accuracy at high densities.

Corrections

To further enhance accuracy, the authors incorporate the Carnahan–Starling (CS) correction into the pressure and chemical potential expressions, which accounts for the non-ideal behavior of hard spheres at high fluid densities.

3. Key Contributions

Generalization to Polydisperse Matrices: The work successfully generalizes the SPT2 framework to handle arbitrary continuous size distributions of matrix obstacles (e.g., rectangular and Schultz–Zimm distributions), moving beyond the limitations of monodisperse or narrow-distribution models.
Development of SPT2b3:** The introduction of the SPT2b3 approximation is the primary theoretical innovation. It resolves the accuracy degradation observed in previous SPT variants at high matrix packing fractions ( $\eta_0 \approx 0.3$ ).
Fully Analytical Expressions: The paper provides closed-form analytical expressions for the excess chemical potential and pressure of the confined fluid, eliminating the need for numerical integration for standard distributions.
Extension of Applicability: The framework extends the valid range of SPT theories to higher matrix packing fractions ( $0.1 \leq \eta_0 \leq 0.3$ ) and broader degrees of polydispersity.

4. Results

The theoretical predictions were validated against Monte Carlo (MC) simulation data (specifically from Leon et al., Ref. [45]) and compared with other integral equation theories.

Rectangular Size Distribution:
- At low matrix packing fractions ( $\eta_0 = 0.1 - 0.2$ ), existing approximations (SPT2b1, SPT2b3*) showed good agreement with MC data.
- As $\eta_0$ increased to 0.25 and 0.3, previous approximations (SPT2a, SPT2b1, SPT2b3*) began to deviate significantly from simulation results, often bracketing the data incorrectly or diverging.
- The new SPT2b3 approximation demonstrated excellent agreement with MC simulations across the entire range of $\eta_0$ (up to 0.3) and various polydispersity widths.
Schultz–Zimm Distribution:
- The SPT2b3** approach was successfully applied to continuous Schultz–Zimm distributions, confirming its versatility for arbitrary continuous size distributions.
- Results confirmed that increasing matrix polydispersity generally lowers the chemical potential at fixed fluid densities, an effect that becomes more pronounced at higher matrix packing fractions.
Trends:
- Increasing the matrix packing fraction ( $\eta_0$ ) significantly decreases the chemical potential of the confined fluid.
- The "pore size distribution" is shown to be sensitive to the width of the matrix particle size distribution, particularly at high $\eta_0$ .

5. Significance and Impact

Reference System for Complex Fluids: The derived analytical expressions provide a robust reference system for studying more complex confined fluids, such as ionic liquids, polymers, and proteins in crowded environments.
Phase Behavior Prediction: By adding attractive or associative interactions to this HS reference, the model can predict phase separation, aggregation, and diffusion in soft matter systems (e.g., protein solutions in porous gels).
Computational Efficiency: The fully analytical nature of the SPT2b3** approach allows for rapid calculation of thermodynamic properties without the computational cost of MC simulations, making it highly valuable for screening and modeling applications.
Foundation for Future Extensions: The authors note that this framework can be extended to describe size-polydisperse fluids in polydisperse matrices and to calculate transport properties like self-diffusion coefficients by adapting the porosity parameters.

In summary, this paper bridges a critical gap in statistical mechanics by providing a highly accurate, analytical, and versatile theory for hard-sphere fluids in realistic, polydisperse porous environments, validated against high-fidelity simulations.