Two-Scale Analysis of the Electrostatics of Dielectric Crystals: Emergence of Polarization Density and Boundary Charges
Shoham Sen, Yang Wang, Timothy Breitzman, Kaushik Dayal
TL;DR
This work develops a rigorous two-scale homogenization framework for electrostatics in dielectric crystals with locally periodic charge distributions. By employing $2$-scale convergence and a dipole scaling $\rho_l=\frac{1}{l}\tilde{\rho}_l$, it derives a homogenized potential $\Phi_0$ governed by $\Delta_x\Phi_0 = \mathbb{1}_{\Omega}\operatorname{div}(\boldsymbol{p}_0)$ together with a boundary condition that includes a bulk polarization $\boldsymbol{p}_0$ and a surface density $\sigma$ arising from partial unit cells; the analysis demonstrates that unit-cell choices alter $\boldsymbol{p}_0$ and $\sigma$ but leave the energy and the homogenized potential invariant. The paper also shows the necessity of a surface charge contribution on $\partial\Omega\setminus\Gamma_d$ and clarifies the non-uniqueness of the dipole moment while preserving a unique macroscopic potential, thereby bridging microscopic charge distributions to macroscopic electrostatic behavior in crystals. These results provide a mathematically rigorous basis for incorporating polarization and boundary effects into continuum models of dielectric crystals and solid electrolytes, with implications for energy storage and electromechanics.
Abstract
Ionic crystals, such as solid electrolytes and complex oxides, are central to modern technologies for energy storage, sensing, actuation, and other functional applications. An important fundamental issue in the atomic and quantum-scale modeling of these materials is defining the macroscopic polarization. In a periodic crystal, the usual definition of the polarization as the first moment of the charge density in a unit cell is found to depend qualitatively - allowing even a change in the sign - and quantitatively on the choice of unit cell. We examine this issue using a rigorous approach based on the framework of 2-scale convergence. By examining the continuum limit of when the lattice spacing is much smaller than the characteristic dimensions of the body, we show that the 2-scale limit provides both a bulk polarization as well as a surface charge density supported on the boundary of the body. Further, different choices of the periodic unit cell of the body lead to correspondingly different partial unit cells at the boundary; these choices give to different bulk polarization and surface charges but compensate such that the electric field and energy are independent of the choice of unit cell.
