---
title: 'Method of Planes (MoP): Simulations & Analysis'
url: https://www.emergentmind.com/topics/method-of-planes-mop
type: topic
---

# Method of Planes (MoP): Simulations & Analysis

Method of Planes (MoP) denotes a plane-based methodology that appears in two distinct technical lineages. In molecular simulation and non-equilibrium molecular dynamics, it is a surface-based procedure for computing density, velocity, momentum flux, and local stress on a finite-area control plane from molecular crossing events and pair interactions [2107.07785]. In analysis and geometric PDE, closely related literature uses the reflection-based method of moving planes to derive symmetry, monotonicity, and rigidity by comparing a function or hypersurface with its reflection across a moving hyperplane [1811.05202]. The shared geometric object is a plane, but the operative role of the plane differs: in microscopic transport it is a control surface carrying fluxes, whereas in PDE it is a moving comparison interface.

## 1. Terminological scope and lineages

Two usages dominate the literature.

| Usage | Plane object | Primary purpose |
|---|---|---|
| Surface-based MoP in microscopic systems | Finite-area control plane \(S_k\) | Surface density, fluxes, stress, and conservation-law-consistent local mechanics |
| Method of moving planes in analysis | Hyperplane \(T_\lambda\) moving through space | Symmetry, monotonicity, rigidity, and quantitative approximate symmetry |

In the microscopic setting, a control plane \(S_k\) is a finite area oriented with outward normal along Cartesian direction \(k\). Surface-averaged quantities are defined by counting molecular passages across \(S_k\) during a sampling interval and by summing pairwise interaction forces whose bond lines intersect the plane [2107.07785]. In the analytical setting, one fixes a direction \(e\in S^{n-1}\), defines hyperplanes \(T_\lambda=\{x\in\mathbb{R}^n:x\cdot e=\lambda\}\), reflects points by \(x^\lambda=x-2(x\cdot e-\lambda)e\), and compares a solution \(u\) with its reflection \(u^\lambda(x)=u(x^\lambda)\) as the plane moves [1811.05202].

The historical trajectories are separate. The reflection method originates in Alexandrov’s reflection principle and soap bubble theorem, Serrin’s symmetry result in potential theory, and the Gidas–Ni–Nirenberg symmetry theory for semilinear elliptic equations, with later quantitative developments and nonlocal extensions [1811.05202]. The surface-flux formulation in molecular simulation is rooted in the Irving–Kirkwood framework and was reformulated for non-equilibrium systems with macroscopic flow through the velocity distribution function, then extended to machine-learning potentials through an explicit pairwise force decomposition compatible with MoP control-volume balances [2107.07785].

## 2. Surface-based MoP in molecular simulation

In molecular dynamics, the Method of Planes computes local mechanical stress on a surface element rather than assigning stress to atom positions. The stress tensor on a surface with outward normal in direction \(k\) is decomposed into kinetic and interaction parts,
\[
\tau_{kl}=\tau_{kl}^{\mathrm{kin}}+\tau_{kl}^{\mathrm{int}},
\]
where \(\tau_{kl}\) is the \(l\)-direction stress acting on that surface [2107.07785]. In the standard equilibrium formulation, the kinetic part is obtained from momentum transfer by molecules crossing the plane, and the interaction part is obtained from pairwise forces whose bond lines cross the plane.

A central construction is the plane-average relation connecting a microscopic particle quantity \(\xi^i\) with a macroscopic field \(\xi(S_k,t)\) defined on \(S_k\):
\[
\int_{-\infty}^{\infty} d\mathbf{v}\, f(S_k,\mathbf{v},t)\,\xi(S_k,t)\,|v_k|
=
\lim_{\delta t\to 0}
\frac{1}{S_k\delta t}
\Bigg\langle
\sum_{i}^{\text{crossing }S_k}
m\,\xi^i
\Bigg\rangle.
\]
Here \(f(\mathbf{x},\mathbf{v},t)\) is the local velocity distribution function. By setting \(\xi=1/|v_k|\), one obtains the surface density,
\[
\rho(S_k,t)=
\lim_{\delta t\to 0}
\frac{1}{S_k\delta t}
\Bigg\langle
\sum_i^{\text{crossing }S_k}
\frac{m}{|v_k^i|}
\Bigg\rangle,
\]
and by setting \(\xi=v_l/|v_k|\), one obtains the surface mass flux,
\[
\rho u_l(S_k,t)=
\lim_{\delta t\to 0}
\frac{1}{S_k\delta t}
\Bigg\langle
\sum_i^{\text{crossing }S_k}
\frac{m\,v_l^i}{|v_k^i|}
\Bigg\rangle.
\]
The streaming velocity on the plane follows as
\[
u_l(S_k,t)=\frac{\rho u_l(S_k,t)}{\rho(S_k,t)}.
\]
Because this estimator is defined for any component \(l\), including \(l\neq k\), any velocity component can be measured on the plane. That point is operationally important: it enables the separation of kinetic momentum flux into advection and stress in non-equilibrium systems with macroscopic flow [2107.07785].

The kinetic contribution to stress is written in terms of peculiar velocity \(c=v-u\), but the plane formulation makes the advection–stress split explicit. In particular,
\[
\tau_{kl}^{\mathrm{kin}}(S_k,t)
=
-\lim_{\delta t\to 0}\frac{1}{S_k\delta t}
\Bigg\langle
\sum_i^{\text{crossing }S_k}
\frac{m\,v_k^i v_l^i}{|v_k^i|}
\Bigg\rangle
+\rho u_lu_k(S_k,t),
\]
with the advective term reconstructed by
\[
\rho u_lu_k(S_k,t)=
\frac{\rho u_l(S_k,t)\cdot \rho u_k(S_k,t)}{\rho(S_k,t)}.
\]
The interaction part is accumulated from forces \(F_l^{ij}\) for pairs \((i,j)\) whose bond line crosses \(S_k\). The resulting formulation defines density, mass flux, momentum flux, and stress directly on a surface. The molecular-dynamics papers emphasize that this surface definition is essential for consistency with the mass and momentum conservation laws in dynamic systems [2107.07785].

## 3. Conservation laws, implementation, and validation in molecular and ML settings

The surface-based MoP is organized around the macroscopic balance laws. Mass conservation is written as
\[
\frac{\partial \rho}{\partial t}+\frac{\partial (\rho u_k)}{\partial x_k}=0,
\]
and momentum conservation as
\[
\frac{\partial (\rho u_l)}{\partial t}
+
\frac{\partial (\rho u_lu_k)}{\partial x_k}
=
\frac{\partial \tau_{kl}}{\partial x_k}
+\rho F_l.
\]
Because the convective momentum flux \(\rho u_lu_k\) and the stress \(\tau_{kl}\) appear as separate surface terms, the microscopic passage statistics across \(S_k\) must be split accordingly. This is the distinctive role of MoP in non-equilibrium molecular dynamics [2107.07785].

In practice, the limit \(\delta t\to 0\) is replaced by the MD timestep \(\Delta t\). Crossings are detected from particle segments between \(x^i(t)\) and \(x^i(t+\Delta t)\), and the paper using velocity Verlet evaluates
\[
\mathbf{v}^i=\frac{\mathbf{x}^i(t+\Delta t)-\mathbf{x}^i(t)}{\Delta t}
\]
to maintain consistency between flux calculations and position updates and to avoid divisions by \(|v_k^i|\approx 0\) [2107.07785]. The same source states two practical assumptions: \(f(\mathbf{x},\mathbf{v},t)\) does not vary appreciably within distances \(|v|\Delta t\) across the oblique pillar used in the derivation, and the crossing velocities are evaluated consistently with the integration scheme.

Validation in a quasi-1D steady-state Couette flow showed that MoP reproduced volume-average density and velocity profiles to within \(\approx 1\%\) in density and \(\approx 0.5\ \mathrm{m/s}\) in velocity. In bulk, \(\tau_{xx}=\tau_{zz}\) and \(\tau_{zx}=\tau_{xz}\) as expected for laminar Couette flow, and \(-\tau_{xx}\) (approximately \(-\tau_{zz}\)) equaled the imposed external pressure \(3.61\ \mathrm{MPa}\). The same study reports that \(\tau_{xx}-\rho u_xu_x\) differed from \(\tau_{zz}\) away from walls, highlighting the necessity of subtracting advection before interpreting \(\tau_{xx}\) as stress, and that the off-diagonal shear stress \(\tau_{zx}\equiv \tau_{xz}\) was constant in the liquid bulk and varied near walls due to solid–liquid friction [2107.07785]. In a quasi-2D moving-contact-line system, MoP revealed a clockwise “caterpillar-like” flow pattern from plane-based velocities, heterogeneous \(\tau_{zx}\) due to viscosity, and strong tensile stress near liquid–vapor interfaces due to interfacial tension.

A later extension derives MoP stress for the MACE potential within the Irving–Kirkwood framework [2509.16257]. For a plane with unit normal \(\mathbf{n}\) at point \(\mathbf{r}_p\), the MoP stress components are written as
\[
P^{\mathrm{MoP}}_{n\beta}(\mathbf{r}_p)
=
\frac{1}{A}
\Big\langle
\sum_i m_i v_{in}v_{i\beta}\,
\delta\!\big(\mathbf{n}\cdot(\mathbf{r}_i-\mathbf{r}_p)\big)
\Big\rangle
+
\frac{1}{4A}
\Big\langle
\sum_{i,j} f_{ij,\beta}\, d_{n;ij}(\mathbf{r}_p)
\Big\rangle,
\]
where \(d_{n;ij}\) is nonzero only if the interaction path crosses the plane. For MACE, the pairwise antisymmetric force decomposition
\[
\mathbf{f}_{ij}
\define
\frac{\partial U}{\partial \mathbf{r}_{ij}}
-
\frac{\partial U}{\partial \mathbf{r}_{ji}}
\]
is used so that \(\mathbf{f}_{ij}=-\mathbf{f}_{ji}\) and \(\mathbf{F}_i=\sum_{j\ne i}\mathbf{f}_{ij}\). The same work reports exact per-timestep momentum conservation for control volumes bounded by MoP planes, and for energy flux it distinguishes a preferred per-atom decomposition \(U=\sum_i U_i\) from a naive pairwise alternative that does not converge to exact energy conservation as \(\Delta t\to 0\) [2509.16257].

The water–ZrO\(_2\) testcase in that study uses a triclinic cell of dimensions \(13.6\times 13.6\times 40.3\) Å, two solid walls of approximately \(7\) Å at top and bottom, \(400\) slabs with \(\Delta z=0.1\) Å, and a temperature of \(500\) K. MoP yields a flat normal stress profile across \(401\) planes, while the local virial profile shows pronounced peaks near the solid–liquid interface; the paper further states that the control-volume momentum balance is verified to machine precision at every timestep [2509.16257].

## 4. Reflection-based moving planes in PDE and geometric analysis

In analysis, the method of moving planes proceeds by fixing a direction \(e\in S^{n-1}\), defining the hyperplane
\[
T_\lambda=\{x\in\mathbb{R}^n:x\cdot e=\lambda\},
\]
the reflected point
\[
x^\lambda=x-2(x\cdot e-\lambda)e,
\]
the cap
\[
\Sigma_\lambda=\{x\in\Omega:x\cdot e>\lambda\},
\]
and the reflected function \(u^\lambda(x)=u(x^\lambda)\). One starts with \(\lambda\) so large that \(T_\lambda\) does not intersect \(\Omega\), decreases \(\lambda\), and stops at the critical position
\[
\lambda^*=\inf\{\lambda:\Sigma'_t\subset \Omega \text{ for all } t\in(\lambda,\infty)\}.
\]
At \(\lambda=\lambda^*\), either an interior tangency or a boundary tangency occurs. The comparison function \(w_\lambda=u-u^\lambda\) is then analyzed in the reflected region. In the model Serrin problem \(\Delta u=-1\) in \(\Omega\), \(u=0\) on \(\partial\Omega\), one has \(\Delta w_\lambda=0\) in the reflected region; the strong maximum principle and Hopf lemma handle interior tangency, while Serrin’s corner lemma handles boundary tangency [1811.05202].

The method yields classical rigidity and symmetry theorems. For Serrin’s overdetermined problem, if \(\Omega\subset\mathbb{R}^n\) is bounded with \(\partial\Omega\in C^{2,\alpha}\) and
\[
\begin{cases}
\Delta u=-1 & \text{in }\Omega,\\
u=0 & \text{on }\partial\Omega,\\
\partial_\nu u=c & \text{on }\partial\Omega,
\end{cases}
\]
then \(\Omega\) is a ball and \(u\) is radial; if \(\Omega=B_R(x_0)\), then
\[
u(x)=\frac{R^2-|x-x_0|^2}{2n},
\qquad
\partial_\nu u|_{\partial B_R}=-\frac{R}{n}.
\]
The same reflection mechanism underlies Gidas–Ni–Nirenberg-type radial symmetry for positive solutions of \(\Delta u+f(u)=0\) in balls or \(\mathbb{R}^n\), and Alexandrov’s theorem that a bounded connected \(C^2\) hypersurface with constant mean curvature is a sphere [1811.05202].

The geometric version rewrites the hypersurface locally as graphs over a tangent plane. Constant mean curvature gives
\[
\operatorname{div}\!\left((1+|\nabla u_i|^2)^{-1/2}\nabla u_i\right)=H,
\]
and the graph difference satisfies a linear uniformly elliptic equation \(Lw=0\), so maximum principles again force coincidence at the critical plane. The survey also records a nonlocal Alexandrov theorem: for \(s\in(0,1/2)\), if the \(s\)-nonlocal mean curvature \(H^s\) is constant on \(\partial\Omega\), then \(\Omega\) is a ball [1811.05202].

## 5. Quantitative, nonlocal, and direct nonlocal variants

A major development is the quantitative method of moving planes, where exact symmetry is replaced by stability estimates. The survey defines the Fraenkel asymmetry
\[
\alpha(\Omega)=\inf_{x\in\mathbb{R}^n,r>0}\frac{|\Omega\Delta B_r(x)|}{|\Omega|},
\]
the Hausdorff distance between \(\partial\Omega\) and a sphere, and the oscillation of boundary data \(\operatorname{osc}_{\partial\Omega}(\partial_\nu u)\). For the Serrin problem, Aftalion–Busca–Reichel define
\[
\operatorname{def}(\Omega)=\|\partial_\nu u-c\|_{C^1(\partial\Omega)}
\]
and obtain concentric balls \(B_r\subset\Omega\subset B_R\) with
\[
R-r\le C |\log \operatorname{def}(\Omega)|^{-1/n}
\]
when the deficit is small. For the discrete Serrin problem, Ciraolo–Magnanini–Sakaguchi obtain the linear bound \(R-r\le C\,\operatorname{def}(\Omega)\). For quantitative Alexandrov, Ciraolo–Vezzoni use \(\operatorname{def}(\Omega)=\operatorname{osc}(H)\) and prove \(R-r\le C\,\operatorname{def}(\Omega)\), together with \(C^{1,\alpha}\)-smallness of \(\partial\Omega\) relative to a sphere. For quantitative nonlocal Alexandrov, Ciraolo–Figalli–Maggi–Novaga define a Lipschitz deficit of \(H^s\) and prove
\[
p(\Omega)\le C(n)\,\operatorname{diam}(\Omega)^{2n+2s+1}\operatorname{def}(\Omega),
\]
with no touching ball or connectedness assumptions [1811.05202].

The analytical toolkit behind these estimates includes quantitative Hopf boundary point lemmas, quantitative Serrin corner lemmas, interior and boundary Harnack inequalities, Carleson-type estimates, interior elliptic regularity, and the use of touching ball geometry to convert smallness of \(w_\lambda\) into small geometric asymmetry. The survey organizes the quantitative argument into approximate symmetry in each direction, localization of critical hyperplanes near an approximate center, control of inner and outer radii, and, for Alexandrov, upgrade to \(C^{1,\alpha}\)-closeness [1811.05202].

The same reflection framework has been adapted to genuinely nonlocal operators. For the uniformly elliptic nonlocal Bellman operator
\[
F_su(x):=\inf_{A\in\mathcal{M}_{\vartheta,\Theta}}
\mathrm{PV}\int_{\mathbb{R}^n}\frac{u(y)-u(x)}{|A^{-1}(y-x)|^{n+2s}}\,dy,
\]
Dai and Qin establish strong maximum principles for antisymmetric functions, maximum principles in bounded and unbounded sets, narrow region principles, decay-at-infinity principles, and sliding methods. These tools yield symmetry, monotonicity, uniqueness, asymptotic results in bounded and unbounded domains, and applications to the uniformly elliptic nonlocal Monge–Ampère operator \(D_s\) of Caffarelli–Charro [2004.02879].

For the fractional Laplacian, a direct method of moving planes works with the integral operator itself rather than the Caffarelli–Silvestre extension. The key ingredients are again strong maximum principles for antisymmetric functions, narrow region principles, and decay at infinity, leading to radial symmetry and monotonicity in critical cases and nonexistence in subcritical cases for equations such as \((-\Delta)^s u=u^p\) [1411.1697]. An analogous direct method has also been developed for the logarithmic Schrödinger operator \((I-\Delta)^{\log}\), with kernel
\[
(I-\Delta)^{\log}u(x)=
c_N\,\mathrm{P.V.}\int_{\mathbb{R}^N}
\frac{u(x)-u(y)}{|x-y|^N\kappa(|x-y|)}\,dy,
\]
where \(\kappa(r)=2^{1-N/2}r^{N/2}\mathcal{K}_{N/2}(r)\). Under the asymptotic condition \(\lim_{|x|\to\infty}u(x)=a\) and \(m>pa^{p-1}\), nonnegative bounded Dini continuous solutions of
\[
(I-\Delta)^{\log}u+mu=u^p
\]
are radially symmetric and radially nonincreasing [2210.09811].

## 6. Limitations, comparisons, and recurring points of confusion

A recurrent source of confusion is terminological. In molecular simulation, MoP is a surface-flux construction tied to control surfaces and conservation laws. In PDE, the moving-planes method is a reflection and comparison method whose decisive ingredients are maximum principles and the behavior of reflected caps. The acronym may overlap in informal usage, but the technical content is distinct.

Within molecular simulation, MoP is repeatedly contrasted with volume-average and local virial approaches. The 2021 formulation states that density, mass, and momentum fluxes are defined on a surface rather than in a volume, and presents this as essential for consistency with mass and momentum conservation in dynamic systems [2107.07785]. The MACE study makes the comparison sharper: the local virial formula is described as strictly valid for whole periodic domains, not for local inhomogeneous bins, and near walls it produces peaks aligned with density oscillations, whereas plane-based MoP maintains a flat normal stress profile under mechanical equilibrium [2509.16257]. The same paper also identifies limitations: long-range electrostatics of Ewald type are not included in the presented MoP, and accurate energy flux for many-body ML models requires a per-atom energy decomposition rather than a naive pairwise one.

Within analysis, the scope of moving planes is broad but conditional. The survey states that typical PDE formulations require \(\partial\Omega\in C^{2,\alpha}\) or comparable touching-ball assumptions, together with validity of maximum principles for the operator. Convexity is not required for the sharp quantitative Alexandrov result, although in qualitative settings it prevents bubbling; for quantitative results, touching-ball conditions are used except in the nonlocal Alexandrov theorem, where nonlocality precludes bubbling [1811.05202]. Related techniques include moving spheres, which is especially useful in conformal or scale-invariant problems, and the sliding method, which uses translations rather than reflections and is employed for monotonicity and symmetry in unbounded domains [1811.05202].

Taken together, these lineages show how a plane can function either as a geometric boundary carrying fluxes or as a moving mirror enforcing comparison. In molecular mechanics, the plane localizes traction and momentum transfer; in PDE and geometry, the plane organizes symmetry and rigidity arguments. The persistence of the term across both areas reflects a shared dependence on planar structure, but the mathematical and physical content is determined by the surrounding theory rather than by the name alone.

Source: https://www.emergentmind.com/topics/method-of-planes-mop