Homoclinic puzzles and chaos in a nonlinear laser model

Abstract: We present a case study elaborating on the multiplicity and self-similarity of homoclinic and heteroclinic bifurcation structures in the 2D and 3D parameter spaces of a nonlinear laser model with a Lorenz-like chaotic attractor. In a symbiotic approach combining the traditional parameter continuation methods using MatCont and a newly developed technique called the Deterministic Chaos Prospector (DCP) utilizing symbolic dynamics on fast parallel computing hardware with graphics processing units (GPUs), we exhibit how specific codimension-two bifurcations originate and pattern regions of chaotic and simple dynamics in this classical model. We show detailed computational reconstructions of key bifurcation structures such as Bykov T-point spirals and inclination flips in 2D parameter space, as well as the spatial organization and 3D embedding of bifurcation surfaces, parametric saddles, and isolated closed curves (isolas).