---
title: Pseudoperiodic SBC in Spherical Systems
url: https://www.emergentmind.com/topics/pseudoperiodic-spherical-boundary-condition-sbc
type: topic
---

# Pseudoperiodic SBC in Spherical Systems

Searching arXiv for papers on "Pseudoperiodic Spherical Boundary Condition" and closely related usage.
Pseudoperiodic Spherical Boundary Condition (SBC) is a context-dependent boundary-construction concept used to impose spherical structure without simulating a full sphere or without inheriting the artifacts of conventional translational periodicity. In the literature represented here, the term denotes at least four distinct but related ideas: an adaptation of revised periodic boundary conditions for curved membranes based on small-angle rotations [1010.0067]; a rotationally invariant spherical-domain protocol with antipodal ghosts for isotropic 3D particle simulations [2507.07888]; an angle-independent phase-closure condition for moving spherical-harmonic cavities [2604.27525]; and an extrapolated interpretation of periodic boundary forcing for the spherical KdV–Burgers equation [2011.14189]. Across these settings, “pseudoperiodic” indicates that periodicity is not realized by ordinary translation on a flat lattice, but by an alternative closure rule tailored to spherical geometry, spherical symmetry, or spherical propagation.

## 1. Terminological scope and core idea

The term SBC does not denote a single universal formalism. In membrane simulation, SBC replaces translational periodicity with small-angle rotations that act as approximate spherical isometries, so that a local patch can represent a curved membrane of radius $R$ [1010.0067]. In isotropic particle simulation, SBC denotes a spherical simulation domain with a ghost-halo protocol based on the antipodal map $G(r) = r - 2R \hat r$, eliminating cubic lattice artifacts while preserving a permeable boundary [2507.07888]. In moving cavity theory, SBC is the requirement that every center-to-boundary-to-center ray accumulates the same round-trip phase $\Phi(\theta)$, thereby preserving spherical-harmonic eigenstructure under motion [2604.27525]. In the KdV–Burgers setting, the phrase is not introduced in the original paper; rather, the provided synthesis treats the conservation-law-based handling of periodic or slowly modulated spherical boundary forcing as informative for a “Pseudoperiodic Spherical Boundary Condition” [2011.14189].

A common structural theme is the substitution of ordinary translational closure by a geometry-aware surrogate. In the membrane case, closure is local and rotational; in the particle case, it is topological and antipodal; in the cavity case, it is phase-theoretic; in the PDE case, it is boundary-forcing-based. This suggests that SBC is best understood as a family of non-translational closure prescriptions adapted to spherical or quasi-spherical settings, rather than as a single algorithmic standard.

## 2. Rotational SBC for curved membranes

In "Approximate Modeling of Spherical Membrane" [1010.0067], SBC is introduced as an adaptation of revised periodic boundary conditions that enforces spherical symmetry locally by replacing translational periodicity with a set of small-angle rotations. Two symmetry operations, $S_1$ and $S_2$, rotate about fixed axes $\hat a_1$ and $\hat a_2$ by angles $\delta\theta_1$ and $\delta\theta_2$, and a pseudoperiodic image is generated by
$$
S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.
$$
For small angles, the rotations commute to linear order, so
$$
S(n_1,n_2)\approx R(\theta,\hat a), \qquad \theta \hat a \approx n_1 \delta\theta_1 \hat a_1 + n_2 \delta\theta_2 \hat a_2.
$$

The geometry is formulated through a local tangent plane at a reference point $r_0 = R n_0$ on a sphere of radius $R$. With tangent vectors $e_1$ and $e_2$, local coordinates $(u,v)$ are mapped to the sphere via the exponential map. Writing
$$
s=\sqrt{u^2+v^2}, \qquad \hat w = \frac{u e_1 + v e_2}{s},
$$
the mapped normal and position are
$$
n(u,v)=n_0\cos(s/R)+\hat w \sin(s/R), \qquad r(u,v)=R\,n(u,v).
$$
Using Rodrigues’ formula, a boundary crossing in the local patch is implemented by a rotation rather than a translation:
$$
T = R(\hat a_1,\pm \delta\theta_1)
$$
for left/right crossings and
$$
T = R(\hat a_2,\pm \delta\theta_2)
$$
for bottom/top crossings.

This construction is approximate but controlled. Its accuracy depends on the separation of scales
$$
L/R \ll 1,
$$
with $L$ the linear size of the patch, and on sufficiently short-ranged interactions [1010.0067]. Non-commutativity of $S_1$ and $S_2$ introduces errors of order $O(\delta\theta^2)$, giving a practical error scaling $O((L/R)^2)$ when $\delta\theta_i \sim L_i/R$. A practical error estimator is to compare energies computed with $S_1^{n_1}S_2^{n_2}$ versus $S_2^{n_2}S_1^{n_1}$.

The principal significance of this SBC is computational. Instead of simulating a complete sphere with millions of interacting particles, one can simulate a small patch while controlling curvature and extracting elastic properties. The method was demonstrated for single- and multilayer graphene, using patches as small as $2$–$32$ atoms, with orders-of-magnitude lower computational cost than full-sphere simulation [1010.0067].

## 3. Curvature energetics and elastic moduli extraction

The membrane SBC of [1010.0067] is closely tied to Helfrich elasticity. For a membrane with spontaneous curvature $c_0$,
$$
E = \int \left[\frac{\kappa}{2}(2H+c_0)^2 + \kappa_G K\right] dA.
$$
For $c_0=0$, a sphere has $H=1/R$, $K=1/R^2$, and $A=4\pi R^2$, yielding
$$
E_{\text{sphere}} = 8\pi \kappa + 4\pi \kappa_G,
$$
or equivalently an energy density
$$
g = \frac{2\kappa+\kappa_G}{R^2}.
$$
For a cylinder of radius $R$,
$$
H = \frac{1}{2R}, \qquad K=0, \qquad g = \frac{\kappa}{2R^2}.
$$
By fitting the energy density versus $R^{-2}$ for both cylinder and sphere, one obtains
$$
\kappa = 2R^2 g_{\text{cyl}}, \qquad \kappa_G = R^2 g_{\text{sphere}} - 2\kappa.
$$

The separation of $\kappa_G$ is topological, via Gauss–Bonnet:
$$
\int K\,dA = 2\pi \chi.
$$
For a sphere, $\chi=2$ and $\int K\,dA=4\pi$; for a cylinder, $\chi=0$. SBC can also be applied to negative Gaussian curvature regions such as $R_1=-R_2=R$, where
$$
g = -\kappa_G/R^2,
$$
providing an independent check on $\kappa_G$ [1010.0067].

The graphene demonstrations reported a fitted mean bending modulus
$$
\kappa \approx 1.61 \ \text{eV}
$$
from cylindrical RPBC and a Gaussian curvature modulus
$$
\kappa_G \approx -0.70 \ \text{eV}
$$
from spherical SBC [1010.0067]. The spherical setup used a 2-atom skewed unit cell with axes $\hat a_1,\hat a_2$ separated by $150^\circ$ and $\delta\theta_i = 2.5 \ \text{\AA}/R'$. Accuracy was checked by small non-commutativity error, agreement between $N=8$ and $N=32$ cells, and independent negative-curvature calculations. Radii down to
$$
R_{\min} \approx 10 \ \text{\AA}
$$
with
$$
\delta\theta_{\max} \approx 15^\circ
$$
were tractable [1010.0067].

For multilayer AB-stacked graphene, the reported values were
$$
\kappa_2 \approx 180 \ \text{eV}, \qquad \kappa_{G2} \approx -140 \ \text{eV},
$$
and
$$
\kappa_3 \approx 690 \ \text{eV}, \qquad \kappa_{G3} \approx -600 \ \text{eV},
$$
in reasonable agreement with the plate-bending estimates
$$
\kappa_n = n\kappa_1 + E h^3 (n^3-n)/12, \qquad
\kappa_{Gn} = n\kappa_{G1} - E h^3 (n^3-n)/12
$$
for layer separation $h=3.4 \ \text{\AA}$ [1010.0067].

A crucial limitation is the strain criterion for solid membranes:
$$
E_s/E_c \sim \frac{[E h \rho^4 /(108R^2)]}{2\kappa+\kappa_G} \ll 1.
$$
If this is violated, nonlocal stress fields dominate and SBC becomes ill-defined [1010.0067]. The method is therefore especially suitable for liquid membranes, and conditionally applicable to solid membranes provided strain energy remains subdominant.

## 4. Antipodal-ghost SBC for isotropic 3D particle simulations

A distinct SBC is introduced in "Pseudoperiodic Spherical Boundary Conditions: Efficient And Isotropic 3D Particle Simulations Without Lattice Artifacts" [2507.07888]. Here the simulation domain is the closed ball
$$
S = \{ r \in \mathbb{R}^3 : \|r\| \le R \},
$$
with a ghost-halo shell
$$
H = \{ r : R-r_c \le \|r\| \le R \},
$$
where $r_c$ is the interaction cutoff. The key map is antipodal:
$$
G(r)=r-2R\hat r, \qquad \hat r = r/\|r\|.
$$

When a particle enters the shell $H$, a ghost copy is created at $G(r_{\text{COM}})$, carrying the same velocities or Brownian increment statistics and the same internal state. While the real particle remains inside $S$, both real and ghost coexist; when the real particle completely exits $S$, it is deleted and the ghost is promoted to real. No coordinate wrapping or discontinuous reimaging occurs [2507.07888].

This is “pseudoperiodic” because it preserves a permeable boundary, mass or particle conservation, and correct cross-boundary interactions, but without a global lattice. Distances are standard Euclidean distances in $\mathbb{R}^3$, not geodesic distances on the sphere. A minimum-image-like effect is achieved by including realized ghost candidates and taking
$$
r_{ij,\text{SBC}} = \min \{\|r_i-r_j\|,\ \|r_i-G(r_j)\|,\ \|G(r_i)-r_j\|\},
$$
with self pairs optionally excluded to avoid an artifact at $r=2R$ [2507.07888].

The formal appeal of this SBC is exact rotational equivariance. The domain is invariant under $Q \in SO(3)$, and the antipodal map satisfies
$$
G(Qr)=QG(r).
$$
Shell-membership decisions depend only on $\|r\|$, and Euclidean distances are rotationally invariant. Therefore the algorithm commutes with rotation, and no lattice orientation exists to imprint preferred directions [2507.07888]. The paper reports that nearest-neighbor-shell azimuth and elevation distributions remain isotropic up to $\phi = 40\%$, whereas conventional PBC exhibit clear fourfold symmetry already at $\phi = 1\%$ [2507.07888].

The paper isolates three PBC artifact sources in crowded regimes: finite-size truncation of long-wavelength fluctuations in a cubic box, artificial recurrence and coordinate-wrapping impulses, and boundary-shape imprint from the cubic lattice [2507.07888]. SBC is designed to remove the last two while maintaining the practical advantages of a permeable boundary.

## 5. Algorithmic properties, validation, and limits of the antipodal formulation

The antipodal-ghost SBC is implemented with standard cell lists or Verlet lists over a cubic bounding volume, culled by the spherical predicate $\|r\| \le R + r_c$ and augmented by ghost particles [2507.07888]. Forces are computed with standard Euclidean pair distances; for general short-range potentials one may use truncated and shifted interactions
$$
u_{\text{tr}}(r)=
\begin{cases}
u(r)-u(r_c), & r<r_c,\\
0, & r\ge r_c.
\end{cases}
$$
No special force-wrapping is required.

The kinematics of ghost motion follow directly from the antipodal map. If a real particle moves from $r_{\ell,r}$ to $r_{f,r}=r_{\ell,r}+dr$, then
$$
r_{\ell,g}=r_{\ell,r}-2R\hat r_{\ell,r}, \qquad
r_{f,g}=r_{f,r}-2R\hat r_{f,r},
$$
and the ghost displacement is
$$
d_g = dr + 2R(\hat r_{\ell,r} - \hat r_{f,r}).
$$
The angle $\xi$ between $dr$ and $d_g$ is explicitly given in the paper’s supporting information [2507.07888]. Special cases include $\xi=0$ for purely radial motion and $\xi \to \pi$ for purely tangential motion at the boundary.

The following reported validation results characterize the method [2507.07888]:

| Observable | SBC result | Comparison |
|---|---|---|
| NNS orientational statistics | Isotropic up to $\phi = 40\%$ | PBC shows fourfold symmetry already at $\phi = 1\%$ |
| Angular momentum fluctuations | $\approx 4\%$ lower global $\Delta L/L$ fluctuation | Lower than PBC |
| Performance | Up to $\sim 60\%$ faster | Faster than optimized MIC at $\phi \gtrsim 10^{-3}$ and $N \in [10^3,10^4]$ |

The collision-rate tests show that SBC and PBC converge to the same collision-rate ground truth as $N$ increases, though PBC appears to converge faster because of spurious recurrence [2507.07888]. In pair statistics at $\phi = 40\%$, both methods have similar short-range structure, but SBC shows a many-body correlation peak near $2R-2r$, while PBC exhibits oscillations around unity before diverging at $L\sqrt{3}$ due to cubic symmetry [2507.07888].

The method currently targets short-range interactions. The paper does not treat Ewald-like long-range electrostatics; because such methods rely on periodicity, SBC would require alternative open-boundary treatments such as reaction-field, multipole or fast multipole, or specialized solvers [2507.07888]. This limitation is conceptually parallel to the membrane SBC of [1010.0067], where standard planar Ewald summation is likewise invalid and multipole-based summations over rotational images are recommended.

## 6. Phase-closure SBC for moving spherical-harmonic cavities

In "Lorentz-FitzGerald Contraction as the Unique Closure Condition for Moving Spherical-Harmonic Cavities" [2604.27525], SBC has a different meaning. It is the requirement that, for a cavity moving uniformly through a nondispersive mechanical wave medium at speed $v=\beta c$, the internal field acquires an angle-independent round-trip phase after a center-to-boundary-to-center ray cycle. For every polar angle $\theta$ relative to the direction of motion,
$$
\Phi(\theta) = \Phi^*,
$$
with $\Phi^*$ constant [2604.27525].

For a cavity boundary described by a radial function $r(\theta)$, pursuit geometry in the medium rest frame gives outward and return travel times
$$
\Delta t_{\text{out}}(\theta)=
\frac{r(\theta)\big(\beta\cos\theta+\sqrt{1-\beta^2\sin^2\theta}\big)}
{c(1-\beta^2)},
$$
$$
\Delta t_{\text{ret}}(\theta)=
\frac{r(\theta)\big(-\beta\cos\theta+\sqrt{1-\beta^2\sin^2\theta}\big)}
{c(1-\beta^2)},
$$
so that the two-way time is
$$
\Delta t_{\text{rt}}(\theta)=
\frac{2r(\theta)\sqrt{1-\beta^2\sin^2\theta}}
{c(1-\beta^2)}.
$$
With $k=\omega/c$ in a nondispersive medium, the round-trip phase is
$$
\Phi(\theta)=
\frac{2k\,r(\theta)\sqrt{1-\beta^2\sin^2\theta}}
{1-\beta^2}.
$$

Imposing $\Phi(\theta)=\Phi^*$ uniquely fixes the angular dependence of the boundary:
$$
r(\theta)=\frac{\Phi^*(1-\beta^2)}{2k}\cdot
\frac{1}{\sqrt{1-\beta^2\sin^2\theta}}.
$$
Writing the transverse scale as $a_\perp \equiv r(\pi/2)$ gives
$$
r(\theta)=\frac{a_\perp\sqrt{1-\beta^2}}{\sqrt{1-\beta^2\sin^2\theta}},
$$
which is an oblate ellipsoid of revolution with
$$
a_\parallel = r(0)=a_\perp\sqrt{1-\beta^2}=a_\perp/\gamma.
$$
Hence
$$
\frac{a_\parallel}{a_\perp}=\frac{1}{\gamma}=\sqrt{1-\beta^2}.
$$
According to the paper, no other angular function $r(\theta)$ satisfies the closure equation; within this framework, the deformation is unique up to overall scale [2604.27525].

Substituting the unique boundary into the round-trip time yields
$$
T=\Delta t_{\text{rt}}=\frac{2a_\perp}{c\sqrt{1-\beta^2}}.
$$
If $a_\perp=R_0$, then
$$
T=\gamma T_0, \qquad T_0=\frac{2R_0}{c}.
$$
The paper therefore presents contraction and period dilation as paired consequences of preserving spherical-harmonic eigenstructure by phase closure [2604.27525].

Here “pseudoperiodic” means phase preservation modulo a single constant across angles rather than geometric periodicity on a fixed domain. This use differs sharply from both particle and membrane SBC, but it preserves the underlying idea that spherical organization is maintained by a nonstandard closure rule.

## 7. Periodic and “pseudoperiodic” spherical boundary forcing in KdV–Burgers dynamics

In "On periodic boundary solutions for cylindrical and spherical KdV-Burgers equations" [2011.14189], the original paper studies periodic perturbation at the boundary for the spherical KdV–Burgers equation
$$
u_t + \frac{1}{t}u = -2uu_x + \varepsilon^2 u_{xx} + \delta u_{xxx}.
$$
The source paper does not introduce the term “pseudoperiodic” or “quasi-periodic.” The provided synthesis states that interpreting slowly varying or modulated periodic forcing as a “Pseudoperiodic Spherical Boundary Condition” is an extrapolation of the paper’s conservation-law framework rather than a theorem stated in the paper [2011.14189].

The boundary forcing of interest is imposed at $x=0$, for example
$$
u(0,t)=A\sin(\omega t),
$$
with $u(b,t)=0$ on a sufficiently large computational domain [2011.14189]. A central identity is the conservation-law form
$$
[t^n u]_t = [t^n(-u^2 + \varepsilon^2 u_x + \delta u_{xx})]_x,
$$
which for spherical geometry ($n=1$) leads to the weighted boundary-flux relation
$$
\frac{1}{T}\int_0^{+\infty} u(x,T)\,dx
=
\frac{1}{T}\int_0^T \frac{t}{T}\,
\big(-\varepsilon^2 u_x(0,t) + u^2(0,t) - \delta u_{xx}(0,t)\big)\,dt.
$$
Defining
$$
J(t):=-\varepsilon^2 u_x(0,t) + u^2(0,t) - \delta u_{xx}(0,t),
$$
the right-hand side is a geometry-weighted time average of the boundary flux [2011.14189].

For suitable dispersion–dissipation balance, the asymptotic profile is a periodic sawtooth-like shock train with decreasing amplitude, preceded by a head shock of constant speed and height $V$ [2011.14189]. In the spherical case, the homothetic envelope is
$$
u_{\text{env}}(x,t)\approx V\sqrt{e}\,
\exp\!\left(-\frac{x}{2Vt\sqrt{e}}\right),
\qquad 0\le x\le Vt,
$$
and the area law gives
$$
\frac{e}{2}V^2 \approx \frac{1}{T}\int_0^T \frac{t}{T}\,J(t)\,dt.
$$
For constant boundary input $u(0,t)=M$ with $u_x(0,t)=0$ at the compression front,
$$
V = \frac{M}{\sqrt e}.
$$
For periodic $u(0,t)=A\sin t$ at $\omega=1$, the paper’s numerics give
$$
M \approx 0.467 A,
$$
hence
$$
V \approx \frac{0.467}{\sqrt e}A
$$
[2011.14189].

The extrapolated “pseudoperiodic” interpretation is that a slowly varying effective mean $M(t)$ could be tracked by a sliding, geometry-weighted average of $J(t)$, leading to a slowly varying head-shock speed
$$
V(t) \approx M(t)/\sqrt e
$$
and head-shock trajectory
$$
x \approx \int^t V(s)\,ds.
$$
Because the synthesis explicitly marks this as a natural extension rather than a theorem, it should be treated as suggestive rather than established [2011.14189].

Taken together, these works show that SBC is not a single boundary condition but a research motif: spherical closure imposed by rotations, antipodal ghosting, angle-independent phase closure, or geometry-weighted boundary forcing. What unifies them is the replacement of flat translational periodicity by a spherical or quasi-spherical prescription that preserves the structure relevant to the problem—local curvature energetics [1010.0067], isotropy and kinematic continuity [2507.07888], spherical-harmonic eigenstructure [2604.27525], or asymptotic shock organization under spherical boundary forcing [2011.14189].

Source: https://www.emergentmind.com/topics/pseudoperiodic-spherical-boundary-condition-sbc