Abstract:This paper presents fast solvers for linear systems arising from the discretization of fractional nonlinear Schrödinger equations with Riesz derivatives and attractive nonlinearities. These systems are characterized by complex symmetry, indefiniteness, and a $d$-level Toeplitz-plus-diagonal structure. We propose a Toeplitz-based anti-symmetric and normal splitting iteration method for the equivalent real block linear systems, ensuring unconditional convergence. The derived optimal parameter is approximately equal to 1. By combining this iteration method with sine-transform-based preconditioning, we introduce a novel preconditioner that enhances the convergence rate of Krylov subspace methods. Both theoretical and numerical analyses demonstrate that the new preconditioner exhibits a parameter-free property (allowing the iteration parameter to be fixed at 1). The eigenvalues of the preconditioned system matrix are nearly clustered in a small neighborhood around 1, and the convergence rate of the corresponding preconditioned GMRES method is independent of the spatial mesh size and the fractional order of the Riesz derivatives.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the problem of rapidly solving the linear systems formed after the discretization of the Riesz fractional nonlinear Schrödinger equation (RFNSE). Specifically, these linear systems are characterized by complex symmetry, indefiniteness, and multi - level Toeplitz - plus - diagonal structures, which pose challenges to the solution.
#### 1. **Research background**
- **Standard Schrödinger equation (SSE)**: It is one of the most important equations in quantum mechanics and is a second - order partial differential equation.
- **Fractional Schrödinger equation (FSE)**: It is obtained by generalizing the standard diffusion operator to the fractional - order diffusion operator. There are two generalization methods: generalization to the fractional Laplacian operator or the Riesz fractional derivative.
- **Importance of numerical methods**: Due to the non - local nature of the fractional - order derivative, it is difficult to obtain the exact solution of the FSE, so numerical methods become crucial. Common numerical methods include the finite element method, the spectral method, the collocation method, and the finite difference method, etc.
#### 2. **Specific problems**
- **Riesz fractional - order nonlinear Schrödinger equation (RFNSE)**:
\[
iu_t+\frac{\partial^\alpha u(x,t)}{\partial |x|^\alpha}+\rho |u|^2u = 0,\quad x\in\mathbb{R}^d,\quad 0 < t\leq T
\]
where \(i=\sqrt{- 1}\), \(x=(x_1,x_2,\ldots,x_d)^\top\in\mathbb{R}^d\), \(\rho>0\) is a real constant (representing the attractive interaction between particles), \(u_0(x)\) is a complex - valued function, and the Riesz fractional derivative is defined as:
\[
\frac{\partial^\alpha u(x,t)}{\partial |x|^\alpha}=-\frac{1}{2\cos(\frac{\alpha\pi}{2})}\sum_{j = 1}^d\left[{}_-\!D_x^\alpha u(x,t)+x_jD_\infty^\alpha u(x,t)\right]
\]
where the left and right Riemann - Liouville fractional derivatives are respectively:
\[
{}_-\!D_x^\alpha u(x,t)=\frac{1}{\Gamma(2-\alpha)}\frac{\partial^2}{\partial x_j^2}\int_{-\infty}^{x_j}\frac{u(x_1,\ldots,\xi,\ldots,x_d,t)}{(x_j-\xi)^{\alpha - 1}}d\xi
\]
\[
x_jD_\infty^\alpha u(x,t)=\frac{1}{\Gamma(2-\alpha)}\frac{\partial^2}{\partial x_j^2}\int_{x_j}^{\infty}\frac{u(x_1,\ldots,\xi,\ldots,x_d,t)}{(\xi - x_j)^{\alpha - 1}}d\xi
\]
- **Linear systems after discretization**:
After discretizing the RFNSE using the linear implicit difference scheme, a series of complex linear systems in the form of \((D_d - T_d + iI)u = b\) will be obtained, where \(D_d\) is a diagonal matrix, \(T_d\) is a symmetric positive - definite multi - level Toeplitz matrix, and due to the parameter \(\rho>0\), \(D_d\)