Abstract This research explores new soliton solutions to the Atangana–Baleanu derivative (ABD) fractional system of equations for ion sound and Langmuir waves (FISLW). We utilize the fractional ABD operator to convert our system into an ordinary differential equations. In recent years, machine learning (ML) evolves significantly in the context of data analysis and computing different solutions, which typically enables systems to operate wisely. Now, we are going to use numerous ML tools including matplotlib. pyplot as plt, scipy.integrate, mpl toolkits.mplot3d, and Axes3D to generate various types of optical solutions by using complete discriminant of the polynomial method. We will also analyze solutions for the hyperbolic function, trigonometric function, Jacobian elliptic function (JEF), and other solitary wave solutions. Solitons have extensive uses in pure and applied mathematics, including nonlinear partial differential equations: the Boussinesq equation, the nonlinear Schrödinger equation, and the sine-Gordon equation, Lie groups, Lie algebras, and differential and algebraic geometry. In addition, we study the chaotic behaviour, i.e., 2D, 3D, time series, Poincarè maps, and sensitivity analysis of our governing model. Sensitivity analysis explores how changes in a system’s variables affects its behaviour.