In this study, the Lengyel-Epstein system is under investigation analytically. This is the reaction–diffusion system leading to the concentration of the inhibitor chlorite and the activator iodide, respectively. These concentrations of the inhibitor chlorite and the activator iodide are shown in the form of wave solutions. This is a reaction†“diffusion model which considered for the first time analytically to explore the different abundant families of solitary wave structures. These exact solitary wave solutions are obtained by applying the generalized Riccati equation mapping method. The single and combined wave solutions are observed in shock, complex solitary-shock, shock singular, and periodic-singular forms. The rational solutions also emerged during the derivation. In the Lengyel-Epstein system, solitary waves can propagate at various rates. The harmony of the system’s diffusive and reactive effects frequently governs the speed of a single wave. Solitary waves can move at a variety of speeds depending on the factors and reaction kinetics. To show their physical behavior, the 3D and their corresponding contour plots are drawn for the different values of constants.