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

• Published in: Communications in nonlinear science & numerical simulation Vol. 68; pp. 15 - 30
• Main Authors: Balogh, Andras, Banda, Jacob, Yagdjian, Karen
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.03.2019; Elsevier Science Ltd
• Subjects: Accelerators; Bubbles; Computer simulation; Discretization; Finite difference method; FORTRAN; Higgs boson equation in the de Sitter spacetime; Higgs bosons; High-performance computations; Initial conditions; Klein-Gordon equation; Particle physics; Quantum physics; Runge-Kutta method; Simulation; Spacetime; Step functions; Theory; Unforced damped Duffing equations; Unforced damped Duffing equations; High-performance computations; Higgs boson equation in the de Sitter spacetime
• ISSN: 1007-5704; 1878-7274
• DOI: 10.1016/j.cnsns.2018.07.011

•High-performance computations shed light to the behavior of solutions.•Compactly supported smooth initial conditions lead to smooth global solution.•Under some conditions the solution forms no bubbles.•Solutions converge to step functions related to unforced, damped Duffing equations. 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.