Abstract:The fundamental model of any periodic crystal is a periodic set of points at all atomic centres. Since crystal structures are determined in a rigid form, their strongest equivalence is rigid motion (composition of translations and rotations) or isometry (also including reflections). The recent isometry classification of periodic point sets used a complete invariant isoset whose size essentially depends on the bridge length, defined as the minimum `jump' that suffices to connect any points in the given set.
We propose a practical algorithm to compute the bridge length of any periodic point set given by a motif of points in a periodically translated unit cell. The algorithm has been tested on a large crystal dataset and is required for an efficient continuous classification of all periodic crystals. The exact computation of the bridge length is a key step to realising the inverse design of materials from new invariant values.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to calculate the bridge length of the set of periodic points, which is a crucial step in achieving the continuous classification of periodic crystals. Specifically, the author proposes a practical algorithm to calculate the bridge length of the motif of the point set within a given periodic translation unit cell. This algorithm has been tested on a large crystal data set and is necessary for the effective continuous classification of all periodic crystals. Accurate calculation of the bridge length is a key step in achieving the reverse design of materials from the invariant values of new materials.
### Background of the Paper and Problem Statement
At the atomic level, all solid crystalline materials can be modeled as periodic point sets, with these points located at the centers of all atoms (chemical properties can be included if necessary). Since the crystal structure is determined in a rigid form, their strongest equivalence relation is rigid motion (a combination of translation and rotation) or isometric transformation (also including reflection). The recent isometric classification of periodic point sets uses a complete invariant - isoset, whose size essentially depends on the bridge length, defined as the maximum "jump" required to connect any two points in a given point set.
### Definition of Bridge Length
For any finite or periodic point set \( S \subset \mathbb{R}^n \), the bridge length \( \beta(S) \) is the minimum distance such that any two points \( p, q \) in \( S \) can be connected by a finite sequence \( p = p_1, p_2, \ldots, p_k = q \), where each Euclidean distance \( |p_i - p_{i + 1}| \leq \beta(S) \) holds for all \( i = 1, \ldots, k - 1 \).
Equivalently, the bridge length is the minimum double - radius such that the union of closed balls with radius \( \frac{1}{2}\beta(S) \) around all points in \( S \) is connected.
### Main Problem
Design an algorithm to calculate the bridge length of an arbitrary periodic point set \( S \) with a fixed unit cell within polynomial time.
### Algorithm Idea
The core idea of this algorithm is to monitor the connectivity of the lifted graph by constructing a labelled quotient graph and gradually adding edges. When the lifted graph \( G \) is connected, the length of the last added edge is the bridge length \( \beta(S) \).
### Time Complexity
The main theorem shows that for any periodic point set \( S \) with \( m \) points in the unit cell \( U \), the bridge length \( \beta(S) \) can be calculated within time \( O(m^2 a(U) n N) \), where \( N \) is the time complexity of the Smith normal form and \( a(U) \) is the aspect ratio of the unit cell \( U \).
### Summary
The main contribution of this paper is to provide an efficient algorithm to calculate the bridge length of periodic point sets, which is of great significance for the continuous classification of periodic crystals and the reverse design of materials.