Visualizing the topological pentagon states of a giant C₅₄₀ metamaterial

The study explores how topological defects, particularly pentagonal disclinations inherent to closed-cage fullerenes, can be engineered and directly visualized using macroscopic topological metamaterials. While buckyballs such as C60 have long been known and their topology recognized for influencing electronic properties, realizing and probing defect-bound states in real molecules is challenging. Topological metamaterials provide a controllable platform to emulate Dirac materials and real-space topological defects. The research question centers on whether a giant, 3D-printed C540 acoustic metamaterial can host and reveal pentagon-localized topological states predicted by gauge-theory treatments of fullerenes, and how these states manifest in real space when probed by sound. The work aims to connect a continuum gauge description and index theorem to experimentally observable defect states, highlighting a topology rooted in real-space lattice defects rather than Bloch band invariants.

The paper situates itself within the broader context of topological phases in photonics, acoustics, and mechanics, where metamaterials have emulated quantum Hall, spin Hall, and valley Hall effects and enabled studies of higher-order states, non-Hermitian phase transitions, and non-Abelian physics. Prior research has established that strain and curvature in graphene introduce pseudo-magnetic gauge fields that markedly affect electronic properties, and that topological defects like disclinations and dislocations can bind states. Gauge-theory and Dirac-equation treatments have been applied to fullerenes to describe their low-energy spectra via effective monopole flux models. Wave-based realizations of topological vortices, disclinations, and dislocations have been demonstrated in photonic and acoustic platforms, laying groundwork for probing defect-bound states in engineered lattices. This work extends those ideas to a 3D printed spherical fullerene metamaterial to directly observe pentagon-localized states.

Theory and modeling: The authors adopt a gauge-theory framework for a spherical fullerene. Starting from a honeycomb lattice, they linearize the tight-binding model to a low-energy Dirac Hamiltonian H = vF σ·∇Ψ. Constructing a truncated icosahedron introduces twelve pentagonal disclinations, which are modeled as effective gauge fluxes. The Dirac operator on the sphere incorporates a monopole field at the center with charge g = 3/2, smearing the total pentagon flux uniformly over the sphere. The spherical Dirac Hamiltonian is H_sphere = vF[σ·(∇ − A) + σ·Δ]Ψ = εΨ, where A is the monopole field and Δ encodes a Kekulé mass term used to tune confinement (Δ = 0 for index analysis). An index theorem after Jackiw-Rossi-Weinberg, index = ∫∇×A·dS = 4g, predicts the number of zero modes. The spectrum depends on angular momentum and Kekulé distortion through ε = √(1 + (l/r)^2) − g^2/r^2 Δ^2 (as presented). The theory is benchmarked against tight-binding and finite-element simulations. Design and fabrication: A giant C540 acoustic metamaterial was 3D printed, forming a closed-cage network of hollow spherical cavities (meta-atoms) connected by cylindrical waveguides, mimicking single and double bonds. The sample comprises 12 printed parts assembled into a sphere of overall diameter 76.7 cm. Each hollow sphere has diameter d0 = 3.6 cm; center-to-center spacing l = 4.5 cm. Intra-cell tube diameter d1 = 0.75 cm is fixed, while inter-cell diameter d2 is varied to implement a Kekulé distortion (design target d2 = 0.5 d1; experimental choice d2 = 0.7 d1 for practicality). Numerical simulations: Finite-element-method (FEM) simulations (COMSOL, pressure acoustics module) compute eigenspectra and eigenfields. Air parameters: density ρ0 = 1.21 kg/m^3, sound speed c0 = 343 m/s at 20 °C. Hard-wall boundary conditions are used on cavity and waveguide surfaces. Intrinsic losses are modeled by introducing c = c0(1 + 0.02i). Comparisons are made among gauge theory (GT), tight-binding method (TBM), and FEM. Experiments: Sound is launched by a loudspeaker inserted into selected cavities; pressure fields are measured using condenser microphones inserted into desired cavities. Signals are amplified (BSWA TECH PA300) and acquired by a digitizer (NI PXI-4499) with LabVIEW processing. Spectra are probed at two sets of locations: (i) 15 selected cavities within the hexagonal network (bulk-like) and (ii) averaged over corners of a selected pentagon (defect-localized). Symmetry-resolved excitations are tested by driving opposing pentagons either in-phase or out-of-phase to distinguish T1g (even) and T1u (odd) triplet character under inversion. Triplet labeling follows Mulliken notation for the Ih icosahedral symmetry.

• Gauge-theory prediction and index theorem: Modeling the twelve pentagon disclinations as a central monopole with charge g = 3/2 yields an index of 4g = 6, predicting six zero modes (three per Dirac point, with valley folding to Γ forming a double Dirac cone). The Kekulé mass term tunes confinement without adding gauge flux.
• Agreement across models: Eigenfrequency spectra from gauge theory (GT), tight-binding (TBM), and finite-element method (FEM) show consistent evolution of isolated pentagonal and hexagonal resonances versus angular momentum and Kekulé distortion.
• Real-space localization: FEM and experiments reveal strong spatial confinement of acoustic fields to pentagons, directly visualizing defect-bound states on the spherical C540 metamaterial.
• Kekulé-enhanced confinement: Implementing a Kekulé distortion (via tube-hopping ratio d2/d1) reduces bulk leakage of defect states and slightly shifts them from zero energy, improving spectral distinguishability.
• Experimental spectra and frequencies: A gapped Dirac spectrum is observed. Fullerene ground-state pentagon modes produce a peak near 850 Hz (c1/e1), with a higher-order excitation near 1225 Hz (c2/e2). Zero-mode triplets split into T1g and T1u components at approximately 869 Hz and 865 Hz, respectively, distinguished by even/odd inversion symmetry in measured and simulated pressure distributions.
• Symmetry-resolved excitation: In-phase/out-of-phase driving of opposing pentagons selectively excites T1g (even under inversion) and T1u (odd) triplets, confirming their symmetry properties.
• Scale and construct: The macroscopic C540 structure features 12 pentagons and 260 hexagons, with overall diameter 76.7 cm, enabling direct probing of phenomena analogous to electronic fullerenes using sound.

The findings validate that real-space topological defects in a spherical fullerene geometry host robust, pentagon-localized zero modes whose count is fixed by a real-space index theorem rather than Bloch band topology. By mapping these modes in a macroscopic acoustic platform, the work bridges gauge-theory predictions (monopole flux model) with directly observable field patterns, confirming the topological origin and symmetry character (T1g/T1u) of the defect states. The Kekulé distortion provides a practical control knob to enhance confinement and spectral isolation, improving defect-state detection. This approach offers insights into how topological defects govern low-energy excitations in fullerenes and demonstrates the utility of metamaterials to emulate and visualize nanoscale emergent phenomena in accessible wave systems.

The authors 3D-printed a giant C540 acoustic metamaterial that emulates fullerene topology and directly visualizes pentagon-localized topological states. A continuum gauge-theory description with an effective monopole field, corroborated by TBM and FEM, predicts and explains six zero modes tied to pentagonal disclinations. Experiments confirm strong defect confinement, symmetry-resolved triplets (T1g/T1u), and tunability via a Kekulé distortion. The platform opens avenues to investigate real-space topological phenomena in carbon allotropes and other complex networks using wave-based probes. Future directions include scaling and adapting the metamaterial framework to graphullerene sheets and other carbon architectures, exploring disorder robustness and non-Hermitian effects, and leveraging symmetry control for targeted state engineering.

• Practical fabrication constraints led to choosing an inter-cell tube diameter ratio d2 = 0.7 d1 (rather than the tighter theoretical design d2 = 0.5 d1), which may slightly affect optimal confinement and spectral placement.
• Intrinsic acoustic losses are present and were modeled (c = c0(1 + 0.02i)); losses broaden spectral features and can influence measured Q-factors and peak amplitudes.
• The gauge-theory treatment homogenizes discrete pentagon fluxes into an effective central monopole on the sphere, an approximation that may omit fine lattice-scale details (tight-binding and FEM were used to cross-validate).