🔬 mesoscale physics

Continuum Limit of Nonlocal Electrostatics in Random Media

This paper derives a deterministic two-scale continuum limit for the electrostatic energy of stationary, ergodic random charge density fields, demonstrating that while microscopic fluctuations can alter local energy contributions, the total energy converges to a sum of the ensemble-averaged microscopic cell energy and the macroscopic electrostatic field energy.

Original authors: Prashant K. Jha, Kaushik Dayal

Published 2026-09-09

📖 8 min read🧠 Deep dive

CC BY 4.0

Original authors: Prashant K. Jha, Kaushik Dayal

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

Electricity is often thought of as a force that acts between distinct objects, like the snap of a static shock or the pull of a magnet. But at the deepest level, electricity is a story of countless tiny charges, positive and negative, jostling within the atoms that make up every solid object. In a perfect crystal, these charges sit in neat, repeating rows, creating a predictable pattern. However, in the real world, especially at temperatures above absolute zero, atoms are never perfectly still. They vibrate, shift, and fluctuate around their average positions. These tiny movements change how the charges interact with one another, creating a complex, random landscape of electrical forces that is difficult to describe with the smooth, continuous equations scientists usually use to model materials.

For decades, physicists have struggled to bridge the gap between this chaotic microscopic reality and the smooth, predictable behavior we see in the macroscopic world. The challenge is that the electrical force between charges is long-range; a charge here can feel the influence of a charge far away, and the sum of all these interactions depends heavily on the specific arrangement of every single atom. When researchers try to average out the microscopic chaos to find a simple rule for the whole material, they often lose critical information about how those tiny fluctuations affect the material's energy. This is particularly important for understanding how materials respond to electric fields, how they store energy, and how they might behave in advanced technologies like sensors or energy storage devices.

In a new study, researchers Prashant K. Jha and Kaushik Dayal have developed a rigorous mathematical framework to solve this problem. They asked a fundamental question: if you have a material where the microscopic charges are arranged randomly but follow certain statistical rules, what does the total electrical energy look like when you zoom out to see the material as a whole? They did not just guess or simulate; they proved that under specific conditions, the chaotic, random energy of the microscopic world settles into two distinct, predictable parts as you look at larger and larger scales. Their work shows that the total energy of such a material is not a single, simple number, but a combination of two very different contributions that tell different stories about the material's internal state.

The researchers began by imagining a material divided into tiny blocks, or cells. Inside each cell, the charges are arranged in a way that is perfectly balanced, meaning the total positive charge equals the total negative charge, so the cell has no net electric charge. However, the charges inside might not be perfectly centered; they could be slightly shifted, giving the cell a tiny magnetic-like property called a dipole moment. In a random material, these shifts vary from cell to cell. Some cells might have their charges shifted one way, others another, and the size of the shift changes randomly. The researchers assumed that while the specific arrangement in any single cell is random, the overall statistical pattern of these shifts remains the same throughout the material. This is a concept known as stationarity, where the rules of the game do not change as you move from one spot to another.

To find the total energy, the team had to separate the interactions into two types. The first type involves charges interacting with other charges within the same tiny cell. The second type involves charges in one cell interacting with charges in a different cell, potentially far away. This separation is crucial because the rules for how these interactions behave change as you look at larger scales. The researchers used a powerful mathematical tool called ergodic averaging. In simple terms, this means that if you look at a large enough volume of the material, the average behavior of the random fluctuations will match the average behavior you would get if you looked at every possible version of the material. This allows them to replace the messy, random details with a clean, deterministic average.

Their analysis revealed a surprising and precise result. The total energy of the material splits into a "local" part and a "nonlocal" part. The local part represents the energy of the charges interacting within their own tiny cells. This energy retains the memory of the microscopic randomness. Even though the material looks smooth from a distance, the energy stored in these tiny cells depends on the specific statistics of how the charges fluctuate. If the charges jitter more or less, or if their distribution changes shape, the local energy changes, even if the overall average position of the charges stays the same. This part of the energy is essentially the sum of the average microscopic energy of every single cell in the material.

The nonlocal part, on the other hand, behaves very differently. This part describes the energy of the interactions between different cells. The researchers found that all the random fluctuations in the individual cells cancel out in this long-range interaction. Instead of seeing a chaotic mess of random dipoles, the material behaves as if it has a single, smooth, and predictable "effective polarization." This is a measure of the average shift of charges across the entire material. The nonlocal energy depends only on this average shift, completely ignoring the specific, random jitters of the individual cells. It is as if the material has forgotten the details of the microscopic chaos and is only responding to the average direction in which the charges are pushed.

A critical piece of the puzzle was a term called the depolarization energy. When a cell is polarized, it creates an internal electric field that opposes the polarization, much like a spring pushing back when you stretch it. In the mathematical breakdown, this depolarization energy appeared in both the local and nonlocal parts, but with opposite signs. In the local calculation, it added a positive amount of energy, while in the nonlocal calculation, it subtracted the exact same amount. When the researchers added the two parts together to get the total energy, these terms canceled each other out perfectly. This cancellation is a profound result: it means that the messy, cell-specific depolarization effects do not appear in the final, total energy of the material. The total energy is simply the sum of the average microscopic energy of the cells plus the energy of the smooth, macroscopic electric field generated by the average polarization.

The study also explored what happens when the random fluctuations are "centered," meaning the charges shift randomly but their average position remains exactly where it was in a perfect, non-random crystal. In this specific case, the researchers showed that the effective polarization remains unchanged. The material looks exactly the same from a distance, and the long-range electric field is identical to that of a perfect crystal. However, the local energy is different. Because the charges are jittering, the average microscopic energy inside each cell increases. This means that two materials could look identical from the outside, have the same overall electric field, and respond to external forces in the same way, yet one could store more internal energy than the other simply because its atoms are vibrating more.

This finding has important implications for how we understand and model materials. It proves that a description based only on the average polarization is incomplete. Such a description can predict the long-range electric fields and how the material interacts with external objects, but it cannot predict the internal energy stored in the microscopic fluctuations. The researchers demonstrated that their mathematical framework works for both perfectly repeating crystals and for materials with random, independent shifts in their atomic positions. They verified that their assumptions hold true for a model where charges move independently from cell to cell, a scenario that mimics the thermal motion of atoms in a solid.

The work provides a rigorous foundation for separating the microscopic and macroscopic worlds in electrostatics. It confirms that while the long-range behavior of a random material is governed by simple, smooth averages, the short-range energy retains a fingerprint of the underlying randomness. This distinction is vital for accurately calculating the energy of materials at finite temperatures, where thermal motion is always present. The researchers' method allows scientists to calculate the total energy of a complex, random material by breaking it down into manageable parts: the average energy of the tiny cells and the energy of the smooth field they create together.

In the end, the study offers a clear, mathematically proven picture of how electricity works in disordered materials. It shows that nature does not simply average out all the chaos; instead, it preserves the microscopic details in the local energy while presenting a smooth, predictable face to the macroscopic world. This dual nature means that to fully understand a material's energy, one must look at both the forest and the trees. The forest, represented by the effective polarization, dictates the long-range forces, but the trees, with their individual, random movements, determine the internal energy stored within the material. This insight refines our understanding of dielectric materials and provides a new tool for predicting how they will behave in real-world applications where temperature and disorder play a significant role.

Technical Summary: Continuum Limit of Nonlocal Electrostatics in Random Media

Problem Statement
Classical continuum electrostatic theories describe the energetic consequences of polarization while suppressing microscopic charge distributions. However, at finite temperatures, charged atoms fluctuate around mean positions, creating random charge density fields. The long-range nature of Coulomb interactions makes the transition from these random microscopic fields to a macroscopic energy nontrivial. Previous analyses have addressed thermodynamic limits for Coulomb matter, homogenization in stochastic lattices, and mean-field models of disordered crystals. This paper addresses a specific gap: deriving a three-dimensional continuum-energy decomposition for a prescribed stationary random charge density that distinguishes the effect of microscopic fluctuations on local electrostatic interactions from their effect on the macroscopic polarization field.

Methodology
The authors employ a two-scale continuum limit approach, separating the domain into three length scales: the microscopic scale $\ell$ (atomic lattice), the material point scale $\varepsilon$ (containing many atoms), and the macroscopic scale $L$ (the body), where $\ell \ll \varepsilon \ll L$ .

Mathematical Framework: The microscopic charge density $\rho(x, y, \omega)$ is modeled as a stationary and ergodic random field under lattice translations. The field is generated by a reference cell charge density $\bar{\rho}$ and a measure-preserving, ergodic lattice action $\{T_k\}$ .
Energy Decomposition: Following Marshall and Dayal [24], the total discrete electrostatic energy is decomposed into:
- Local Energy ( $E_{loc}$ ): Interactions between charges within the same material point.
- Nonlocal Energy ( $E_{nloc}$ ): Interactions between charges in distinct material points.
Averaging and Convergence: The authors utilize lattice ergodic theorems to show that spatial averages over large numbers of microscopic cells converge almost surely to deterministic ensemble averages. They analyze the convergence of finite-volume potentials to infinite-volume limits, distinguishing between the full-space potential and boundary-corrected potentials arising from the finite material point regions.
Assumptions: The derivation relies on three key assumptions:
- (A1) Continuity of the charge density map and exact charge neutrality within every reference cell realization.
- (A2) Existence and $L^2$ -integrability of the microscopic Coulomb potential.
- (A3) Strong convergence of finite-volume field corrections to the gradient of the potential of a uniformly polarized cube.

Key Contributions and Results
The paper establishes that under the stated assumptions, the discrete local and nonlocal energies converge almost surely to deterministic continuum limits.

Local Limit: The local energy converges to the sum of two terms:
- The ensemble-averaged microscopic Coulomb energy: $\frac{1}{2} \int_\Omega \mathbb{E}[\int_Q \bar{\rho}\bar{h} \, ds] \, dx$ .
- A cell-depolarization term: $\frac{1}{6} \int_\Omega |p(x)|^2 \, dx$ , where $p(x)$ is the effective polarization.
Nonlocal Limit: The nonlocal energy converges to the macroscopic electrostatic field energy minus a depolarization term:
- $\frac{1}{2} \int_{\mathbb{R}^3} |\nabla \phi|^2 \, dx - \frac{1}{6} \int_\Omega |p(x)|^2 \, dx$ , where $\phi$ is the potential generated by the effective polarization $p$ .
Total Energy Cancellation: A critical result is the cancellation of the cell-depolarization terms ( $+\frac{1}{6}|p|^2$ and $-\frac{1}{6}|p|^2$ ) in the total energy. Consequently, the total continuum energy is:
$E_{total} = \frac{1}{2} \int_\Omega \mathbb{E}\left[\int_Q \bar{\rho}\bar{h} \, ds\right] dx + \frac{1}{2} \int_{\mathbb{R}^3} |\nabla \phi|^2 dx$
This consists of the ensemble-averaged microscopic cell energy and the macroscopic field energy.
Information Retention: The local and nonlocal limits retain different information. The nonlocal limit depends solely on the effective polarization $p(x)$ , losing the specific statistics of the microscopic fluctuations. The local limit retains the ensemble-averaged statistics of the charge density.
Verification Cases:
- Periodic Case: The framework recovers the known periodic two-scale limit when the probability space is a single point.
- Random Displacement Model: The authors verify the assumptions for a model of independent cellwise random displacements. They demonstrate that centered microscopic fluctuations can alter the local energy (by changing the ensemble-averaged microscopic term) without changing the effective polarization or the nonlocal energy.

Significance and Claims
The paper claims that ergodic averaging does not eliminate all dependence of electrostatic energy on microscopic randomness. While the long-range (nonlocal) interactions are fully determined by the deterministic effective polarization, the short-range (local) interactions retain the statistical signature of the microscopic charge distribution.

The authors emphasize that continuum descriptions based only on effective polarization are insufficient to capture the energetic effects of microscopic fluctuations. Such descriptions correctly predict the macroscopic field energy and its dependence on boundaries but fail to account for the ensemble-averaged microscopic cell energy. This distinction is particularly relevant for finite-temperature systems where polarization fluctuations interact with depolarizing fields and surface screening. The work provides a rigorous mathematical foundation for separating these effects in random media, extending previous periodic analyses to stationary and ergodic random settings.

More like this