Abstract:We analyse the numbers of closed paths of length $k\in\mathbb{N}$ on two important regular lattices: the hexagonal lattice (also called $\textit{graphene}$ in chemistry) and its dual triangular lattice. These numbers form a moment sequence of specific random variables connected to the distance of a position of a planar random flight (in three steps) from the origin. Here, we refer to such a random variable as a $\textit{random eigenvalue}$ of the underlying lattice. Explicit formulas for the probability density and characteristic functions of these random eigenvalues are given for both the hexagonal and the triangular lattice. Furthermore, it is proven that both probability distributions can be approximated by a functional of the random variable uniformly distributed on increasing intervals $[0,b]$ as $b\to\infty$. This yields a straightforward method to simulate these random eigenvalues without generating graphene and triangular lattice graphs. To demonstrate this approximation, we first prove a key integral identity for a specific series containing the third powers of the modified Bessel functions $I_n$ of $n$th order, $n\in\mathbb{Z}$. Such series play a crucial role in various contexts, in particular, in analysis, combinatorics, and theoretical physics.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to study the spectral properties of graphene (i.e., the hexagonal lattice) and its dual triangular lattice. Specifically, the author analyzes the number of closed paths of length \(k\) on these lattices, and these path numbers form the moment sequence of specific random variables. These random variables are related to the position distance of a planar random flight (within three steps) starting from the origin. The author gives explicit formulas for the probability density function and the characteristic function of these random variables and proves that these probability distributions can be approximated by the functions of random variables uniformly distributed on the interval \([0, b]\) as \(b \to \infty\). This provides a simple method for simulating these random variables without generating graphene and triangular lattice graphs.
To verify this approximation method, the author first proves a key integral identity for a specific series of the modified Bessel function \(I_n\) that contains the cube of \(I_n\). This type of series plays an important role in analysis, combinatorics, and theoretical physics.
### Main contributions
1. **Analysis of the number of closed paths**: The author analyzes the number of closed paths of length \(k\) on the hexagonal lattice and the triangular lattice, and these path numbers form the moment sequence of specific random variables.
2. **Probability density and characteristic function of random variables**: Gives explicit formulas for the probability density function and the characteristic function of these random variables.
3. **Approximation method**: Proves that these probability distributions can be approximated by the functions of random variables uniformly distributed on the interval \([0, b]\) as \(b \to \infty\).
4. **Integral identity of the modified Bessel function**: Proves a key integral identity for a specific series of the modified Bessel function \(I_n\), which provides a theoretical basis for the approximation method.
### Key formulas
- **Integral representation of the number of closed paths**:
\[
\mu_{2k}(H) = \frac{1}{(2\pi)^3} \int_0^{2\pi} \int_0^{2\pi} \int_0^{2\pi} \left| e^{i\phi_1} + e^{i\phi_2} + e^{i\phi_3} \right|^{2k} d\phi_1 d\phi_2 d\phi_3
\]
- **Probability density function of random variable \(X\)**:
\[
f_X(x) = \int_0^\infty t x J_0(t x) J_0^3(t) dt, \quad x \geq 0
\]
- **Probability density functions of random variables \(H\) and \(T\)**:
\[
f_H(x) = \begin{cases}
\frac{\sqrt{3} |x|}{\pi (3 + x^2)} \, _2F_1 \left( \frac{1}{3}, \frac{2}{3}; 1; \frac{x^2 (9 - x^2)^2}{(3 + x^2)^3} \right), & \text{if } x \in [-3, 3] \\
0, & \text{else}
\end{cases}
\]
\[
f_T(x) = \begin{cases}
\frac{\sqrt{3}}{\pi (3 + x)} \, _2F_1 \left( \frac{1}{3}, \frac{2}{3}; 1; \frac{x (9 - x)^2}{(3 + x)^3} \right), & \text{if } x \in [0, 9] \\
0, & \text{else}
\end{cases}
\]
- **Approximation of random variable \(Y_b\)**:
\[
Y_b = \cos(X_b) + \cos(\beta X_b) + \cos((1 + \beta) X_b)