In this study, the \(\phi ^{6}\) -model expansion method is showed to be useful for finding solitary wave solutions to the Klein–Gordon (KG) equation. We develop a variety of solutions, including Jacobi elliptic functions, hyperbolic forms, and trigonometric forms, so greatly enhancing the range of exact solutions attainable. The 2D, 3D, and contour plots clearly show different types of solitary waves, like bright, dark, singular, and periodic solitons. This gives us a lot of information about how the KG equation doesn’t work in a straight line. Our findings highlight the \(\phi ^{6}\) model as a powerful tool to study nonlinear wave equations, improve our understanding of their complex dynamics, and increase the scope for theoretical exploration. The \(\phi ^{6}\) model expansion technique is exceptionally adaptable and may be utilised for a wide array of nonlinear partial differential equations. Despite its versatility, the technique may not be applicable to all nonlinear PDEs, especially those that do not meet the specified requirements or structures manageable by this technique. In theoretical physics, particularly in field theory and quantum mechanics, the Klein–Gordon equation is a classical model. By studying this model, we can illustrate the waves and particles movements at relativistic speeds. Among other areas, its significance in cosmology, quantum field theory, and the study of nonlinear optics are widely considered. Additionally, it provides exact solutions and nonlinear dynamics have various applications in applied mathematics and physics. The study is novel because it provides a new understanding of the complex behaviours and various waveforms of the controlling model by means of detailed evaluation. Future research could focus on further exploring the stability and physical implications of these solutions under different conditions, thereby advancing our knowledge of nonlinear wave phenomena and their applications in physics and beyond.
