In this study, complex dynamics of Briggs–Rauscher reaction system is investigated analytically and numerically. First, the Briggs–Rauscher reaction system is reduced into a new nonlinear parametric oscillator. The Melnikov method is used to derive the condition of the appearance of horseshoe chaos in the cases ω = Ω and ω ≠ Ω. The performed numerical simulations confirm the obtained analytical predictions. Second, the prediction of coexisting attractors is investigated by solving numerically the new nonlinear parametric ordinary differential equation via the (...) fourth-order Runge–Kutta algorithm. As results, it is found that the new nonlinear chemical system displays various coexisting behaviors of symmetric and asymmetric attractors. In addition, the system presents a rich variety of bifurcations phenomena such as symmetry breaking, symmetry restoring, period doubling, reverse period doubling, period-m bubbles, reverse period-m bubbles, intermittency, and antimonotonicity. On the contrary, emerging chaotic band attractors and period-1, period-3, period-9, and period-m bubbles routes to chaos occur in this system. (shrink)