In this paper, we address various optical soliton solutions and demonstrate the different dynamics of solitary waves to a (3+1)-dimensional nonlinear Schrödinger equation (NLSE) with parabolic law (NLSE) using a newly created powerful and effective method named as the extended generalized Riccati equation mapping method. This technique presents an organized manner to reveal the essential dynamics. There is great significance of the nonlinear Schrödinger equations and their several formulations in numerous fields of science, particularly in nonlinear optics, optical fibers, quantum electronics, and plasma physics. Through the use of numerical simulations and mathematical analysis, we explore the characteristics and behavior of these solitary wave solutions in a variety of scientific contexts. These results demonstrate the essential complexity of the governing equation and yield novel derived solutions. These solutions contribute to a better understanding of nonlinear wave phenomena by highlighting the fundamental dynamics establishing solitary waves in the NLSE. To enhance our wider knowledge, we provide effective graphic representations of the nonlinear wave structures in the derived solutions utilizing a variety of graphs, including 3D, 2D, and density plots. Moreover, a specific transformation has been applied to transform the system into a planar dynamical system, and several phase portraits have been presented to examine its behavior. Furthermore, upon introducing a perturbed term, chaotic behavior has been observed across different parameter values through both two-dimensional and three-dimensional graphics.