Static vector solitons in a topological mechanical lattice

The paper addresses how multiple topological solitons can interact and propagate in macroscopic mechanical systems, a capability largely unexplored in prior mechanical metamaterials that typically support a single scalar kink. The authors motivate the work by highlighting topological solitons (kinks, vortices, skyrmions) as robust, particle-like excitations relevant in condensed matter, soft matter, and optics, and by recent interest in mechanical metamaterials and topological Maxwell lattices that exhibit boundary-localized zero modes. Building on analogies to the SSH model and to coupled double Peierls chains with sublattice symmetry breaking and chiral solitons, the authors propose a mechanical lattice that can host multiple, interacting solitons with vector character. The purpose is to realize static vector solitons through coupled fields, uncover their domain-wall dynamics connecting multiple degenerate ground states, and connect topological polarization with multistability, thereby enabling robust information processing concepts in mechanical platforms.

The paper surveys foundational studies of topological solitons in polymers, dislocation dynamics, magnetic domain walls, DNA, and liquid crystals, emphasizing their stability and collision behaviors. It reviews mechanical metamaterials demonstrating programmable localized deformations, transition fronts, and topological mechanical lattices that enable mechanical signal transmission, logic, and actuation. Prior one-dimensional topological Maxwell lattices analogous to SSH chains have shown boundary modes and conductive solitons (single scalar kinks). More recent atomic-scale studies of coupled double chains (indium wires) with sublattice symmetry breaking revealed chiral solitons protected by a nontrivial configuration space with multiple ground states, though long-range interactions complicate mechanical analogues. In other fields (ferroelectric phase transitions, nonlinear optics, multispecies condensates), coupled field equations can yield vector solitons with superposed kink components. The authors position their work to bring these multi-soliton/vector-soliton phenomena into macroscopic mechanical lattices.

• Lattice design: A one-dimensional topological mechanical (Maxwell) lattice is constructed from hinged rhombic linkages (rigid, gray) connected by elastic bonds (red). Each unit cell has a common node with two translational DoFs (u, v). Geometry is parameterized by rhombus side length r, height 2h, and lattice spacing a, with glide-reflection symmetry. Interconnected bonds impose two length constraints between neighboring cells; admissible linkage shaping region satisfies u^2 + (v + h)^2 ≤ 4r^2.
• Topological characterization: The compatibility matrix C maps infinitesimal nodal displacements to bond elongations. For periodic lattices, the Fourier-transformed C(k) yields the dynamic matrix D(k) = m^{-1} C(k) k_s C^†(k). Distorted lattices show a gapped phonon spectrum; boundary-localized zero modes are predicted by the winding number of arg det C(k): w = (1/2πi) ∮ ∂k ln det C(k) dk. The topological polarization R_w a identifies the edge localization of zero modes. Finite lattices under open boundaries have two zero modes (via Maxwell-Calladine index theorem) that localize according to polarization, with distinct penetration depths extracted from complex roots of det C(k) = 0.
• Nonlinear map iteration (zero-energy constraints): Finite deformations are computed by enforcing constant bond lengths (Eqs. (2) and (3), not explicitly printed here) forming a recursive map (u_{n+1}, v_{n+1}) = f(u_n, v_n). Given the first node perturbation along a boundary zero mode, subsequent node positions are solved iteratively (Newton’s method), producing domain walls and homogeneous end states (fixed points) corresponding to multiple ground states.
• Continuum model: In the long-wavelength limit (lattice spacing a small relative to deformation variation scale), nodal coordinates are approximated by fields u(x), v(x), x = n + 1/2. Expanding to second order in fields and first order in derivatives yields coupled first-order ODEs for normalized fields U = u/u_m, V = v/v_m with coefficient α: a[(U+V)^2 − (U+V)^2] + U' + V' = 0 and a[(U−V)^2 − (U−V)^2] + U' − V' = 0. These are solved using variables U±V to obtain analytic vector-soliton solutions that are linear superpositions of two kinks centered at x1 and x2 (Eqs. (7), (8)).
• Effective potential: An effective potential energy density P(u, v) is derived showing four degenerate minima (ground states) and identifying soliton trajectories as minimum-energy pathways on a rugged energy surface; polarization changes correspond to traversing domain walls between minima.
• Experimental prototype: A physical lattice with parameters (r, h, a) = (20 mm, 31.8 mm, 40 mm) was fabricated (photosensitive resin linkages, metal screws). Boundary displacements were applied to drive transformations along chosen pathways; edge stiffness asymmetry was measured by applying 2 mm displacements at edges and recording forces.
• Generalization: Analytical criteria for topological transitions and existence of multiple ground states are derived (e.g., boundary lines v = ± u_u/h; invariance P(u, v) = P(−u, v)). A distorted lattice example (u, v) = (1.23 u_0, 0.25 v_0) demonstrates persistence of four ground states and vector solitons with asymmetric amplitudes.

• Discovery of static vector solitons in a macroscopic topological mechanical lattice: solitons arise from coupled fields and have components that are sums/differences of two kinks (superposed hyperbolic profiles), unlike classic scalar φ^4 or sine-Gordon kinks.
• Fourfold degenerate ground states: The lattice exhibits four homogeneous configurations (±u_0, 0) and (0, ±v_0), located in distinct topological phases in the phase diagram. Domain walls (solitons) interpolate between these states, reversing winding numbers and topological polarization.
• Independent boundary modes and branched motion pathways: Horizontal and vertical zero modes are independent with different penetration depths; finite deformations along these modes yield distinct domain walls, including kink-only and kink–antikink co-propagating structures.
• Continuum theory and analytic solutions: Coupled first-order ODEs admit vector-soliton solutions composed of two sub-solitons centered at x1 and x2. Special case x1 = x2 reduces to a single-component kink (SSH-like). Analytical profiles quantitatively agree with numerical iterations of the discrete zero-energy constraints across multiple cases.
• Effective potential energy perspective: The system’s energy landscape features four degenerate minima; vector solitons trace minimum-energy pathways on a rugged surface, explaining phase switching and multi-domain structures.
• Experimental validation: A physical prototype demonstrates reversible transformations along designed pathways, realizing transitions such as (0, −v_0) → (u_0, 0) → (−u_0, 0). Strong edge-stiffness asymmetry consistent with polarization: under 2 mm edge displacements, measured forces were ~0.05 N (soft edge) vs ~8.80 N (hard edge), a difference of over two orders of magnitude.
• Generality: Analytical criteria predict topological transition boundaries (e.g., v = ± u_u/h) and the existence of four ground states over a broad admissible geometric region. Distorted lattices retain inversion symmetry P(u, v) = P(−u, v) and support multi-phase domain structures; vector solitons with unequal component amplitudes persist.

The findings show that introducing coupled degrees of freedom into a topological Maxwell lattice enables static vector solitons that transfer zero modes from the edges into the bulk while switching between distinct topological phases. The superposition of sub-solitons provides a mechanism for multi-soliton interactions and programmable domain structures in macroscopic mechanical systems. By mapping the nonlinear deformations to a continuum theory with an effective multistable potential, the work bridges topological polarization and multistability: domain walls correspond to minimum-energy trajectories between degenerate minima. This perspective explains observed phenomena such as multiple homogeneous states, kink–antikink coexistence, and reversible microstructural evolution. The agreement among phonon/topological analysis, discrete zero-energy iterations, continuum analytic solutions, and experiments substantiates the theoretical framework and indicates robustness across parameter variations, suggesting routes toward topologically protected, low-energy mechanical information processing and computing.

The study demonstrates transformable topological mechanical lattices that host static vector solitons composed of superposed hyperbolic kinks. The finite isostatic lattice exhibits independent boundary modes, four degenerate ground states, and branched motion pathways; domain walls act as vector solitons that switch topological phases and reverse polarization. A continuum model yields closed-form vector-soliton solutions that match discrete simulations, and an effective potential explains multi-phase domain structures as minimum-energy pathways. Experiments validate reversible transformations and strong edge-stiffness asymmetry. Future work could explore: dynamic regimes and boosted vector solitons; integration into soliton-based logic gates, diodes, and algebraic mechanical computing; extensions to other transformable lattices (e.g., fixed boundaries, geared or rotating-block architectures); active or thermal actuation of polarized modes; and broader exploration of multistability–topology interplay in higher-dimensional or non-Hermitian mechanical systems.

• Continuum approximation: The theoretical model uses second-order field expansions and first-order derivatives; deviations from predicted topological transition lines (e.g., v = ± u_u/h) occur for large deformations due to truncation.
• Prototype nonidealities: Assembly tolerances and friction affect ease of deformation and pathway selection; horizontal displacements were more effective experimentally, indicating practical constraints.
• Geometric admissibility: Multiple ground states and vector solitons exist within an admissible geometric region; extreme distortions may violate connectivity or admissibility, limiting generality.
• Static focus: Emphasis is on static (zero-energy) vector solitons; while dynamic responses are mentioned, full dynamical characterization (e.g., dispersion, radiation losses) is beyond the present scope.
• Long-range interactions absent: Unlike atomic chains, the mechanical analogue neglects long-range interactions, which may alter behaviors in other realizations.