Pulgon-tools: A toolkit for analysing and harnessing symmetries in quasi-1D systems

This paper introduces Pulgon-tools, an open-source software package that fills a gap in automated symmetry analysis for quasi-one-dimensional periodic systems by integrating structure generation, line-group symmetry detection, irreducible representation analysis, and harmonic interatomic force constant corrections within a unified framework.

The Gist

Imagine you are an architect trying to design a futuristic city made of tiny, infinitely long tubes (like nanotubes) and wires. In the world of 3D buildings (like skyscrapers), we have a very strict rulebook called "Space Groups" that tells us how to arrange bricks so the building is stable and symmetrical. We have excellent computer programs that can look at a 3D building and instantly say, "Ah, this is a perfect cube with a mirror on the side!"

But here's the problem: Nanotubes aren't 3D buildings. They are long, skinny, one-dimensional tubes. They have their own unique set of rules, called Line Groups.

Think of it like this:

• 3D Crystals are like a checkerboard. You can slide a tile left, right, up, or down, and it looks the same.
• Nanotubes are like a spiral staircase or a screw. To make it look the same, you can't just slide it; you have to twist it while you slide it.

Until now, scientists didn't have a good "rulebook" or a computer program to understand these spiral rules. They tried to use the 3D checkerboard rules, but it didn't work. The computer would look at a perfect spiral staircase and say, "This is a messy, broken building," because it couldn't see the twist.

Enter "Pulgon-tools."

The paper introduces a new, free software toolkit called Pulgon-tools (named after "Pulgon," a playful nod to the complexity of these shapes). It's like a specialized detective and a master builder rolled into one, designed specifically for these 1D tube-like structures.

Here is how it works, broken down into four simple jobs:

1. The Architect (Structure Generation)

Sometimes you want to build a new nanotube from scratch.

• The "Lego" Method: You tell the software, "I want a tube that twists 8 times and slides up a little bit." The software takes a few basic atoms (your "Lego bricks") and automatically snaps them together into a perfect, symmetrical tube based on those rules.
• The "Roll-Up" Method: Imagine taking a flat sheet of paper (like a graphene sheet) and rolling it into a tube. Depending on the angle you roll it, you get different types of tubes (some straight, some spiral). Pulgon-tools can take your "rolling instructions" (called chiral indices) and instantly build the 3D tube for you, ensuring the atoms are spaced perfectly.

2. The Detective (Symmetry Detection)

This is the most important part. You hand the software a picture of a messy, real-world nanotube (maybe it was built by a robot and has a few atoms slightly out of place).

• The Old Way: The software would look at it and say, "I can't find any symmetry. It's broken."
• The Pulgon Way: The detective looks closer. It ignores the tiny wobbles and asks: "If I twist this tube by 36 degrees and slide it up, does it look the same?"
• It figures out the Secret Code of the tube. Is it a simple slide? A screw twist? A mirror reflection? It identifies the exact "Line Group" family the tube belongs to. It's like looking at a fingerprint and instantly knowing which person it belongs to, even if the print is smudged.

3. The Translator (Irreps and Character Tables)

Once the detective knows the "Secret Code" (the symmetry), the software acts as a translator.

• In physics, we need to know how electrons or vibrations (sound) move through these tubes. The math for this is incredibly complex.
• Pulgon-tools translates the physical shape of the tube into a mathematical ID card (called Irreducible Representations).
• Think of this like a menu. Instead of giving you a giant, confusing list of every single atom's movement, it gives you a simplified menu of "flavors" of movement. It tells you: "Okay, this tube has 8 distinct 'flavors' of vibration. Here is exactly how they behave." This makes calculating properties (like heat or electricity) much faster and easier.

4. The Fixer (IFC Correction)

When scientists simulate these tubes on a computer, the numbers often get a little "sloppy" due to rounding errors.

• The Problem: A perfect tube should have a vibration that costs zero energy (like a sound wave moving through air). But because of computer errors, the simulation might say, "This tube is vibrating with a ghost energy," which is physically impossible.
• The Fix: Pulgon-tools acts as a spell-checker for physics. It looks at the messy numbers and gently nudges them back into line with the laws of physics (conservation of momentum). It ensures that if the tube is supposed to be perfectly smooth, the math says it is perfectly smooth, removing those "ghost" errors.

Why does this matter?

Nanotubes and nanowires are the future of technology. They are used in super-strong materials, tiny electronics, and medical sensors. But to design them, we need to understand their unique "spiral" symmetries.

Before Pulgon-tools, scientists were trying to fit a square peg (1D tubes) into a round hole (3D software). They had to do things manually, which was slow and prone to error.

Pulgon-tools is the new tool that finally gives scientists a dedicated, automated way to:

1. Build these tubes correctly.
2. Identify their hidden symmetries.
3. Translate their physics into simple math.
4. Fix any computer errors.

It's like giving a master carpenter a specialized saw for cutting curved wood, instead of forcing them to use a hammer. It makes the job of designing the materials of tomorrow much faster, more accurate, and much more fun.

Technical Summary

1. Problem Statement

Quasi-one-dimensional (quasi-1D) periodic systems, such as nanotubes and nanowires, possess intrinsic symmetries described by line groups rather than the three-dimensional space groups used for bulk crystals.

• The Gap: While mature software exists for 3D space-group detection (e.g., spglib, FINDSYM), there is a lack of widely available, automated tools for identifying line-group symmetries directly from atomistic structures. Existing 3D tools often fail to recognize unique quasi-1D operations (screw rotations, glide reflections along the axis, and arbitrarily subdivided rotational angles), resulting in incorrect low-symmetry classifications.
• Secondary Challenges: There are no standardized tools for:
  1. Generating quasi-1D structures consistent with line-group symmetry.
  2. Computing irreducible representations (irreps) and character tables specific to line groups.
  3. Enforcing physical invariance conditions (translational and rotational) on harmonic interatomic force constants (IFCs) in these systems, which is critical for accurate vibrational (phonon) calculations.

2. Methodology

The authors present Pulgon-tools, an open-source Python package designed to address these gaps through a modular, pipeline-oriented architecture. The software integrates four complementary components:

A. Structure Generation

The toolkit offers two distinct methods for constructing quasi-1D systems:

1. General Symmetry-Based Construction: Generates structures from a minimal atomic motif and specified line-group generators (axial point group and generalized translational group). It applies group expansion algorithms to build the full periodic structure.
2. Chiral Roll-Up Method: Specifically designed for $MX_2$-type nanotubes (e.g., $MoS_2$). It takes chiral indices $(n, m)$ of a 2D hexagonal parent lattice, calculates helical line-group parameters, maps the planar structure to cylindrical coordinates, and uses constrained optimization (Newton-Krylov method) to adjust bond lengths for structural stability.

B. Symmetry Detection

The detection algorithm decomposes the problem into two stages to identify the full line group $L = Z \rtimes P$:

1. Generalized Translational Group ($Z$) Identification:
   • Determines the primitive translation along the periodic axis.
   • Detects screw rotations ($C_{n_1} | f$) and glide reflections ($\sigma_V | L/2$) by testing candidate monomers and translation distances.
2. Axial Point Group ($P$) Identification:
   • Analyzes the primitive cell as a rigid cluster.
   • Determines the highest-order rotational symmetry ($C_{n_2}$), perpendicular twofold axes ($U_d$), and mirror planes to classify the group into one of the 13 line-group families (e.g., $L_{qp}$, $L_{nm}$, etc.).

C. Irreducible Representations (Irreps) and Character Tables

• Implements the analytical framework of Damnjanović and Milošević for all 13 line-group families.
• Constructs irreps and character tables on-the-fly (no precomputed lookup tables) using physically meaningful quantum numbers:
  • Axial wave vector $k \in (-\pi/L, \pi/L]$.
  • Angular momentum index $m$.
  • Parity labels ($\Pi_V, \Pi_H, \Pi_U$) corresponding to vertical mirrors, horizontal mirrors, and improper rotations.

D. Interatomic Force Constant (IFC) Correction

• Problem: Numerical noise or incomplete symmetry enforcement in first-principles calculations often leads to violations of translational and rotational invariance, causing unphysical phonon behaviors (e.g., imaginary frequencies at the Brillouin zone center).
• Solution: The module enforces acoustic sum rules, Born–Huang rotational sum rules, and Huang invariance conditions.
• Algorithm: Formulates the correction as a constrained quadratic optimization problem ($\min \|H - H_{raw}\|^2$ subject to $C \cdot H = 0$). It uses the CVXPY framework with the OSQP solver to efficiently handle large, sparse systems, ensuring the corrected IFCs deviate minimally from the raw data while satisfying physical laws.

3. Key Contributions

• First Automated Line-Group Detector: Pulgon-tools is the first widely available library to automatically detect line-group symmetries directly from atomistic coordinates, bridging the gap between 3D crystallography tools and quasi-1D physics.
• Comprehensive Toolkit: It unifies structure generation, symmetry analysis, representation theory, and force constant correction into a single ecosystem.
• Specialized $MX_2$ Generator: Provides a robust, dedicated tool for generating single- and multilayer $MoS_2$-type nanotubes from chiral indices, including bond-length relaxation.
• Physical Invariance Enforcement: Offers a rigorous method to correct IFCs for quasi-1D systems, eliminating spurious imaginary phonon modes and ensuring conservation of linear and angular momentum.
• Open-Source & Modular: Written in Python with an Apache 2.0 license, it features both a command-line interface and a Python API for integration into custom workflows (e.g., with Phonopy).

4. Results and Validation

The paper validates the software through several application examples:

• Structure Generation: Successfully generated $MoS_2$ nanotubes with various chiralities (zigzag, armchair, chiral) and general symmetry-based structures, verifying geometric consistency and bond-length stability.
• Symmetry Detection: Correctly identified the line group of a $(5,5)$ single-wall carbon nanotube (SWCNT) as family 4 ($T_{10}(L/2)C_5$), successfully navigating multiple monomer candidates and distinguishing between screw and glide symmetries.
• Irreps Generation: Generated character tables for a $MoS_2$ $(5,0)$ nanotube, correctly labeling 8 irreducible representations with quantum numbers $(k, m, \Pi)$ at $k=0$.
• IFC Correction: Applied to a $(12,12)$ $MoS_2$ nanotube. The raw IFCs (from machine-learning potentials) showed spurious imaginary frequencies near the $\Gamma$ point. After correction, these were eliminated, restoring the correct linear and quadratic acoustic branches in the phonon dispersion.

5. Significance

Pulgon-tools addresses a critical bottleneck in the computational study of low-dimensional materials. By providing a reliable method to identify line-group symmetries, it enables:

• Reduced Computational Cost: Symmetry-adapted sampling of the Brillouin zone.
• Physical Accuracy: Ensuring that vibrational and electronic calculations respect the fundamental symmetries of quasi-1D systems.
• Advanced Analysis: Facilitating symmetry-resolved studies of thermal transport, electronic band structures, and topological properties using irreducible representations.

The software fills a significant void in the computational materials science landscape, offering a standardized, automated approach for the modeling and analysis of nanotubes, nanowires, and other quasi-1D periodic systems.