High-performance implementation of a Runge-Kutta finite-difference scheme for the Higgs boson equation in the de Sitter spacetime

Abstract High performance computations are presented for the Higgs boson equation in the de Sitter Spacetime using explicit fourth order Runge–Kutta scheme on the temporal discretization and fourth order finite difference discretization in space. In addition to the full, ( 3 + 1 ) − dimensional equation we also examine the ( 1 + 1 ) − dimensional radial solutions. The numerical code for the ( 3 + 1 ) − dimensional equation is programmed in CUDA Fortran and is performed on NVIDIA Tesla K40c GPU Accelerators. The radial form of the equation is simulated in MATLAB. The numerical results demonstrate the existing theoretical result that under certain conditions bubbles form in the scalar field. We also demonstrate the known blow-up phenomena for the solutions of the related semilinear Klein–Gordon equation with imaginary mass. Our numerical studies suggest several previously not known properties of the solution for the Higgs boson equation in the de Sitter spacetime for which theoretical proofs do not exist yet: 1. smooth solution exists for all time if the initial conditions are compactly supported and smooth; 2. under some conditions no bubbles form; 3. solutions converge to step functions related to unforced, damped Duffing equations.

[1]  N. Konyukhova,et al.  On the stability of a self-similar spherical bubble of a scalar Higgs field in de Sitter space , 2005 .

[2]  Jalil Rashidinia,et al.  Numerical solution of the nonlinear Klein-Gordon equation , 2010, J. Comput. Appl. Math..

[3]  S. Lipson,et al.  The impossibility of free tachyons , 1971 .

[4]  K. Yagdjian.  On the Global Solutions of the Higgs Boson Equation , 2010, 1009.3076.

[5]  S. Jerez-Galiano,et al.  Two finite-difference schemes that preserve the dissipation of energy in a system of modified wave equations , 2010 .

[6]  Salah M. El-Sayed,et al.  A numerical solution of the Klein-Gordon equation and convergence of the decomposition method , 2004, Appl. Math. Comput..

[7]  Patrick Joly,et al.  Construction and Analysis of Fourth-Order Finite Difference Schemes for the Acoustic Wave Equation in Nonhomogeneous Media , 1996 .

[8]  P. J. van der Houwen,et al.  The development of Runge-Kutta methods for partial differential equations , 1996 .

[9]  L. Vázquez,et al.  Analysis of Four Numerical Schemes for a Nonlinear Klein-Gordon Equation , 1990 .

[10]  H. Epstein,et al.  de Sitter Tachyons and Related Topics , 2014, 1403.3319.

[11]  Wilhelm Schlag,et al.  Numerical study of the blowup/global existence dichotomy for the focusing cubic nonlinear Klein–Gordon equation , 2010, 1011.2015.

[12]  Ricardo Carretero-González,et al.  Numerical Stability of Explicit Runge-Kutta Finite Difference Schemes for the Nonlinear Schrödinger Equation , 2011, ArXiv.

[13]  Muhammet Yazıcı,et al.  Approximate solutions to the nonlinear Klein-Gordon equation in de Sitter spacetime , 2016 .

[14]  K. Yagdjian.  The semilinear Klein-Gordon equation in de Sitter spacetime , 2009, 0903.0089.

[15]  Mehdi Dehghan,et al.  Numerical solution of the nonlinear Klein-Gordon equation using radial basis functions , 2009 .

[16]  C. Bogey,et al.  A family of low dispersive and low dissipative explicit schemes for flow and noise computations , 2004 .

[17]  S. Coleman.  Aspects of symmetry: Selected Erice lectures: Contents , 1985 .

[18]  Victor Pereyra,et al.  Symbolic generation of finite difference formulas , 1978 .

[19]  P. C. Ray,et al.  SOLUTION OF NON-LINEAR KLEIN-GORDON EQUATION WITH A QUADRATIC NON-LINEAR TERM BY ADOMIAN DECOMPOSITION METHOD , 2009 .

[20]  Mehdi Dehghan,et al.  Fourth-order compact solution of the nonlinear Klein-Gordon equation , 2009, Numerical Algorithms.

[21]  R. Fry,et al.  Smooth bump functions and the geometry of banach spaces , 2002 .

[22]  L. Hörmander,et al.  Lectures on Nonlinear Hyperbolic Differential Equations , 1997 .