In this paper, we primarily focus on the nonlinear paraxial wave equation in Kerr media, a generalization of the nonlinear Schrödinger equation utilized to characterize the dynamics of optical beam propagation. By using three potent analytical methods, namely, the Sine-Gordon expansion method, the functional variable method, and the Bernoulli ([Formula: see text])-expansion method, numerous novel soliton solutions are derived. These solutions, comprising hyperbolic, trigonometric and exponential functions, represent significant additions to the field of optics. Furthermore, we elucidate the physical characteristics of these novel optical soliton solutions by presenting a series of three-dimensional (3D) and two-dimensional (2D) graphs with appropriate parameter values.