Efficient determination of the critical parameters and the statistical quantities for Klein-Gordon and sine-Gordon equations with a singular potential using generalized polynomial chaos methods

TL;DR

This paper introduces a generalized polynomial chaos method combined with spectral collocation techniques to accurately and efficiently determine critical parameters and mean solutions for Klein-Gordon and sine-Gordon equations with singular potentials, addressing challenges posed by singularities.

Contribution

The work develops a novel gPC-based approach with spectral collocation to precisely identify critical parameters in wave equations with singular potentials, overcoming convergence issues caused by singularities.

Findings

The method accurately determines critical potential strength for Klein-Gordon equations.

It effectively finds initial velocities leading to different solution behaviors in sine-Gordon equations.

Numerical results demonstrate high efficiency and precision of the proposed approach.

Abstract

We consider the Klein-Gordon and sine-Gordon type equations with a point-like potential, which describes the wave phenomenon in disordered media with a defect. The singular potential term yields a critical phenomenon--that is, the solution behavior around the critical parameter value bifurcates into two extreme cases. Pinpointing the critical value with arbitrary accuracy is even more challenging. In this work, we adopt the generalized polynomial chaos (gPC) method to determine the critical values and the mean solutions around such values. First, we consider the critical value associated with the strength of the singular potential for the Klein-Gordon equation. We expand the solution in the random variable associated with the parameter. The obtained partial differential equations are solved using the Chebyshev collocation method. Due to the existence of the singularity, the Gibbs phenomenon appears in the solution, yielding a slow convergence of the numerically computed critical value. To deal with the singularity, we adopt the consistent spectral collocation method. The gPC method, along with the consistent Chebyshev method, determines the critical value and the mean solution highly efficiently. We then consider the sine-Gordon equation, for which the critical value is associated with the initial velocity of the kink soliton solution. The critical behavior in this case is that the solution passes through (particle-pass), is trapped by (particle-capture), or is reflected by (particle-reflection) the singular potential if the initial velocity of the soliton solution is greater than, equal to, or less than the critical value, respectively. We use the gPC mean value rather than reconstructing the solution to find the critical parameter. Numerical results show that the critical value can be determined efficiently and accurately by using the proposed method.