Tiexiang Li

Jingkai Zhang, Xing-Long Lyu, Tiexiang Li et al.

Band-structure calculations for three-dimensional anisotropic phononic crystals require the repeated solution of large elastic generalized eigenvalue problems along Bloch paths. In standard staggered-grid discretizations, anisotropic coupling may involve derivative components located at incompatible grid positions, so additional interpolation or averaging closures are often introduced. This paper proposes a parity-blocked rotated staggered discretization based on four Bloch-periodic body-diagonal differences. The directional derivatives are reconstructed from these diagonal differences, leading to a Hermitian $B_hC_hB_h^H$ generalized eigenvalue formulation that incorporates anisotropic derivative coupling without separate interpolation closures. On even grids, when the stiffness and mass matrices are nodewise local multiplication matrices, the body-diagonal shifts preserve two independent parity invariants. The discrete velocity space is then decomposed exactly into four mutually independent block subspaces, and the full discrete spectrum can be recovered by solving the four smaller eigenvalue problems and merging their spectra. The full and block formulations are further organized in a unified Fourier SVD framework, which supports $Γ$-point zero-mode treatment, shift-invert Krylov iteration, inner PCG solves, and GPU matrix-vector products. Numerical experiments for a three-dimensional two-phase anisotropic phononic crystal show that the block implementation preserves the full-space spectrum while substantially reducing the wall-clock time. The results demonstrate that the proposed method provides a structured and efficient solver for large-scale band-structure computations of three-dimensional anisotropic phononic crystals.

Teng-Chao Sun, Tiexiang Li, Wen-Wei Lin et al.

In this paper, we propose a spectral framework that embeds 1D and 2D quasiperiodic Helmholtz eigenvalue problems into higher-dimensional (2D and 4D) periodic spaces via the projection method \cite{jiang2014numerical, jiang2024numerical}. To effectively map the elevated high-dimensional states back to the original physical space, we establish a novel validation framework based on the weighted expectation of pointwise Rayleigh quotients. Supported by comprehensive error and spectral analysis, we demonstrate that the eigenvalues derived from this expectation align more authentically with the original quasiperiodic model, ultimately yielding a more appropriate and reliable eigenpair solution. Numerical experiments on continuous media demonstrate that our approach offers an accurate, robust, and scalable tool for solving quasiperiodic Helmholtz eigenvalue problems.

Tsung-Ming Huang, Tiexiang Li, Wei-De Li et al.

Calculation of band structure of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalue. In this paper, we try to address this problem in such a broad context that band structure of three dimensional isotropic photonic crystals with all 14 Bravais lattices can be efficiently computed in a unified framework. We uncover the delicate machinery behind several key results of our work and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments in our framework. In this work particular effort is made on reformulating the Bloch boundary condition for all 14 Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by systematic use of commutativity among them, which has been overlooked previously, and reducing eigen-decomposition of double-curl operator to the canonical form of a 3x3 complex skew-symmetric matrix under unitary congruence. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm.

Zhen-Chen Guo, Tiexiang Li, Ying-Ying Zhou

To solve the Bethe-Salpeter eigenvalue problem with distinct sizes, two efficient methods, called ΓQR algorithm and Γ-Lanczos algorithm, are proposed in this paper. Both algorithms preserve the special structure of the initial matrix $H=\begin{bmatrix}A & B-\overline{B} & -\overline{A}\end{bmatrix}$, resulting the computed eigenvalues and the associated eigenvectors still hold the properties similar to those of $H$. Theorems are given to demonstrate the validity of the proposed two algorithms in theory. Numerical results are presented to illustrate the superiorities of our methods.