Pseudoperiodic SBC in Spherical Systems

• Pseudoperiodic SBCs are innovative boundary prescriptions that replace standard translational periodicity with geometry-aware closures tailored for spherical domains.
• They enable efficient simulation of curved membranes, isotropic particle domains, and moving cavities by preserving spherical symmetry and critical energy characteristics.
• These methods reduce computational cost and lattice artifacts while providing actionable insights into elasticity, shock dynamics, and boundary forcing 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 (Koskinen et al., 2010); a rotationally invariant spherical-domain protocol with antipodal ghosts for isotropic 3D particle simulations (Dedola et al., 10 Jul 2025); an angle-independent phase-closure condition for moving spherical-harmonic cavities (Meucci, 30 Apr 2026); and an extrapolated interpretation of periodic boundary forcing for the spherical KdV–Burgers equation (Samokhin, 2020). 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$ (Koskinen et al., 2010). 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 (Dedola et al., 10 Jul 2025). 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 (Meucci, 30 Apr 2026). 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” (Samokhin, 2020).

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" (Koskinen et al., 2010), 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

$G(r) = r - 2R \hat r$0

The geometry is formulated through a local tangent plane at a reference point $G(r) = r - 2R \hat r$1 on a sphere of radius $G(r) = r - 2R \hat r$2. With tangent vectors $G(r) = r - 2R \hat r$3 and $G(r) = r - 2R \hat r$4, local coordinates $G(r) = r - 2R \hat r$5 are mapped to the sphere via the exponential map. Writing

$G(r) = r - 2R \hat r$6

the mapped normal and position are

$G(r) = r - 2R \hat r$7

Using Rodrigues’ formula, a boundary crossing in the local patch is implemented by a rotation rather than a translation:

$G(r) = r - 2R \hat r$8

for left/right crossings and

$G(r) = r - 2R \hat r$9

for bottom/top crossings.

This construction is approximate but controlled. Its accuracy depends on the separation of scales

$\Phi(\theta)$0

with $\Phi(\theta)$1 the linear size of the patch, and on sufficiently short-ranged interactions (Koskinen et al., 2010). Non-commutativity of $\Phi(\theta)$2 and $\Phi(\theta)$3 introduces errors of order $\Phi(\theta)$4, giving a practical error scaling $\Phi(\theta)$5 when $\Phi(\theta)$6. A practical error estimator is to compare energies computed with $\Phi(\theta)$7 versus $\Phi(\theta)$8.

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 $\Phi(\theta)$9–$S_1$0 atoms, with orders-of-magnitude lower computational cost than full-sphere simulation (Koskinen et al., 2010).

3. Curvature energetics and elastic moduli extraction

The membrane SBC of (Koskinen et al., 2010) is closely tied to Helfrich elasticity. For a membrane with spontaneous curvature $S_1$1,

$S_1$2

For $S_1$3, a sphere has $S_1$4, $S_1$5, and $S_1$6, yielding

$S_1$7

or equivalently an energy density

$S_1$8

For a cylinder of radius $S_1$9,

$S_2$0

By fitting the energy density versus $S_2$1 for both cylinder and sphere, one obtains

$S_2$2

The separation of $S_2$3 is topological, via Gauss–Bonnet:

$S_2$4

For a sphere, $S_2$5 and $S_2$6; for a cylinder, $S_2$7. SBC can also be applied to negative Gaussian curvature regions such as $S_2$8, where

$S_2$9

providing an independent check on $\hat a_1$0 (Koskinen et al., 2010).

The graphene demonstrations reported a fitted mean bending modulus

$\hat a_1$1

from cylindrical RPBC and a Gaussian curvature modulus

$\hat a_1$2

from spherical SBC (Koskinen et al., 2010). The spherical setup used a 2-atom skewed unit cell with axes $\hat a_1$3 separated by $\hat a_1$4 and $\hat a_1$5. Accuracy was checked by small non-commutativity error, agreement between $\hat a_1$6 and $\hat a_1$7 cells, and independent negative-curvature calculations. Radii down to

$\hat a_1$8

with

$\hat a_1$9

were tractable (Koskinen et al., 2010).

For multilayer AB-stacked graphene, the reported values were

$\hat a_2$0

and

$\hat a_2$1

in reasonable agreement with the plate-bending estimates

$\hat a_2$2

for layer separation $\hat a_2$3 (Koskinen et al., 2010).

A crucial limitation is the strain criterion for solid membranes:

$\hat a_2$4

If this is violated, nonlocal stress fields dominate and SBC becomes ill-defined (Koskinen et al., 2010). 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" (Dedola et al., 10 Jul 2025). Here the simulation domain is the closed ball

$\hat a_2$5

with a ghost-halo shell

$\hat a_2$6

where $\hat a_2$7 is the interaction cutoff. The key map is antipodal:

$\hat a_2$8

When a particle enters the shell $\hat a_2$9, a ghost copy is created at $\delta\theta_1$0, carrying the same velocities or Brownian increment statistics and the same internal state. While the real particle remains inside $\delta\theta_1$1, both real and ghost coexist; when the real particle completely exits $\delta\theta_1$2, it is deleted and the ghost is promoted to real. No coordinate wrapping or discontinuous reimaging occurs (Dedola et al., 10 Jul 2025).

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 $\delta\theta_1$3, not geodesic distances on the sphere. A minimum-image-like effect is achieved by including realized ghost candidates and taking

$\delta\theta_1$4

with self pairs optionally excluded to avoid an artifact at $\delta\theta_1$5 (Dedola et al., 10 Jul 2025).

The formal appeal of this SBC is exact rotational equivariance. The domain is invariant under $\delta\theta_1$6, and the antipodal map satisfies

$\delta\theta_1$7

Shell-membership decisions depend only on $\delta\theta_1$8, and Euclidean distances are rotationally invariant. Therefore the algorithm commutes with rotation, and no lattice orientation exists to imprint preferred directions (Dedola et al., 10 Jul 2025). The paper reports that nearest-neighbor-shell azimuth and elevation distributions remain isotropic up to $\delta\theta_1$9, whereas conventional PBC exhibit clear fourfold symmetry already at $\delta\theta_2$0 (Dedola et al., 10 Jul 2025).

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 (Dedola et al., 10 Jul 2025). 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 $\delta\theta_2$1 and augmented by ghost particles (Dedola et al., 10 Jul 2025). Forces are computed with standard Euclidean pair distances; for general short-range potentials one may use truncated and shifted interactions

$\delta\theta_2$2

No special force-wrapping is required.

The kinematics of ghost motion follow directly from the antipodal map. If a real particle moves from $\delta\theta_2$3 to $\delta\theta_2$4, then

$\delta\theta_2$5

and the ghost displacement is

$\delta\theta_2$6

The angle $\delta\theta_2$7 between $\delta\theta_2$8 and $\delta\theta_2$9 is explicitly given in the paper’s supporting information (Dedola et al., 10 Jul 2025). Special cases include $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$0 for purely radial motion and $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$1 for purely tangential motion at the boundary.

The following reported validation results characterize the method (Dedola et al., 10 Jul 2025):

Observable	SBC result	Comparison
NNS orientational statistics	Isotropic up to $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$2	PBC shows fourfold symmetry already at $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$3
Angular momentum fluctuations	$S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$4 lower global $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$5 fluctuation	Lower than PBC
Performance	Up to $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$6 faster	Faster than optimized MIC at $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$7 and $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$8

The collision-rate tests show that SBC and PBC converge to the same collision-rate ground truth as $S(n_1,n_2)=S_1^{n_1}S_2^{n_2}.$9 increases, though PBC appears to converge faster because of spurious recurrence (Dedola et al., 10 Jul 2025). In pair statistics at $G(r) = r - 2R \hat r$00, both methods have similar short-range structure, but SBC shows a many-body correlation peak near $G(r) = r - 2R \hat r$01, while PBC exhibits oscillations around unity before diverging at $G(r) = r - 2R \hat r$02 due to cubic symmetry (Dedola et al., 10 Jul 2025).

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 (Dedola et al., 10 Jul 2025). This limitation is conceptually parallel to the membrane SBC of (Koskinen et al., 2010), 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" (Meucci, 30 Apr 2026), SBC has a different meaning. It is the requirement that, for a cavity moving uniformly through a nondispersive mechanical wave medium at speed $G(r) = r - 2R \hat r$03, the internal field acquires an angle-independent round-trip phase after a center-to-boundary-to-center ray cycle. For every polar angle $G(r) = r - 2R \hat r$04 relative to the direction of motion,

$G(r) = r - 2R \hat r$05

with $G(r) = r - 2R \hat r$06 constant (Meucci, 30 Apr 2026).

For a cavity boundary described by a radial function $G(r) = r - 2R \hat r$07, pursuit geometry in the medium rest frame gives outward and return travel times

$G(r) = r - 2R \hat r$08

$G(r) = r - 2R \hat r$09

so that the two-way time is

$G(r) = r - 2R \hat r$10

With $G(r) = r - 2R \hat r$11 in a nondispersive medium, the round-trip phase is

$G(r) = r - 2R \hat r$12

Imposing $G(r) = r - 2R \hat r$13 uniquely fixes the angular dependence of the boundary:

$G(r) = r - 2R \hat r$14

Writing the transverse scale as $G(r) = r - 2R \hat r$15 gives

$G(r) = r - 2R \hat r$16

which is an oblate ellipsoid of revolution with

$G(r) = r - 2R \hat r$17

Hence

$G(r) = r - 2R \hat r$18

According to the paper, no other angular function $G(r) = r - 2R \hat r$19 satisfies the closure equation; within this framework, the deformation is unique up to overall scale (Meucci, 30 Apr 2026).

Substituting the unique boundary into the round-trip time yields

$G(r) = r - 2R \hat r$20

If $G(r) = r - 2R \hat r$21, then

$G(r) = r - 2R \hat r$22

The paper therefore presents contraction and period dilation as paired consequences of preserving spherical-harmonic eigenstructure by phase closure (Meucci, 30 Apr 2026).

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" (Samokhin, 2020), the original paper studies periodic perturbation at the boundary for the spherical KdV–Burgers equation

$G(r) = r - 2R \hat r$23

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 (Samokhin, 2020).

The boundary forcing of interest is imposed at $G(r) = r - 2R \hat r$24, for example

$G(r) = r - 2R \hat r$25

with $G(r) = r - 2R \hat r$26 on a sufficiently large computational domain (Samokhin, 2020). A central identity is the conservation-law form

$G(r) = r - 2R \hat r$27

which for spherical geometry ($G(r) = r - 2R \hat r$28) leads to the weighted boundary-flux relation

$G(r) = r - 2R \hat r$29

Defining

$G(r) = r - 2R \hat r$30

the right-hand side is a geometry-weighted time average of the boundary flux (Samokhin, 2020).

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 $G(r) = r - 2R \hat r$31 (Samokhin, 2020). In the spherical case, the homothetic envelope is

$G(r) = r - 2R \hat r$32

and the area law gives

$G(r) = r - 2R \hat r$33

For constant boundary input $G(r) = r - 2R \hat r$34 with $G(r) = r - 2R \hat r$35 at the compression front,

$G(r) = r - 2R \hat r$36

For periodic $G(r) = r - 2R \hat r$37 at $G(r) = r - 2R \hat r$38, the paper’s numerics give

$G(r) = r - 2R \hat r$39

hence

$G(r) = r - 2R \hat r$40

(Samokhin, 2020).

The extrapolated “pseudoperiodic” interpretation is that a slowly varying effective mean $G(r) = r - 2R \hat r$41 could be tracked by a sliding, geometry-weighted average of $G(r) = r - 2R \hat r$42, leading to a slowly varying head-shock speed

$G(r) = r - 2R \hat r$43

and head-shock trajectory

$G(r) = r - 2R \hat r$44

Because the synthesis explicitly marks this as a natural extension rather than a theorem, it should be treated as suggestive rather than established (Samokhin, 2020).

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 (Koskinen et al., 2010), isotropy and kinematic continuity (Dedola et al., 10 Jul 2025), spherical-harmonic eigenstructure (Meucci, 30 Apr 2026), or asymptotic shock organization under spherical boundary forcing (Samokhin, 2020).