Numerical simulation of coupled Klein–Gordon–Schrödinger equations: RBF partition of unity

Abstract

The coupled Klein–Gordon–Schrödinger equations have significant implications in quantum field theory, particle physics, cosmology, and nonlinear dynamics. In this study, we propose an efficient method for numerically simulating this system. The proposed approach involves employing the radial basis function partition of unity for spatial discretization. This method utilizes scaled Lagrange basis functions with polyharmonic spline kernels, taking advantage of the scalability property of the polyharmonic kernel to ensure stability in the approximation process. The resulting spatially discretized system yields a nonlinear time-dependent set of differential equations. To solve this system, we combine an implicit central finite difference scheme with a predictor–corrector procedure to overcome the nonlinearity. In the numerical results section, we assess the efficiency, accuracy, and versatility of the proposed method by conducting several simulations in large domains over extended periods. We present and compare these results with those obtained using alternative methods, demonstrating the effectiveness of the proposed approach in accurately solving the coupled Klein–Gordon–Schrödinger equations.