The concept of cubic–quartic (CQ) optical solitons emerged a couple of years ago out of extreme necessity when the chromatic dispersion (CD) effect ran low. This led to replenishing the low CD count with the CQ dispersive effect so that the necessary balance between CD and self-phase modulation (SPF) was sustained for the existence and propel of solitons for long distances through optical fibers. Later, this concept was applied to the Lakshmanan–Porsezian–Daniel (LPD) model [1,2,3,4], not including the well-known and most visible one: the nonlinear Schrödinger’s equation. The vector-coupled LPD equation in birefringent fibers was considered with the aid of the extended version of Jacobi’s elliptic function expansion scheme in [1], where Jacobi’s doubly periodic wave solutions are found. These solutions, in the limiting case, give rise to dark solitons, singular solitons or periodic solutions. The generalized LPD model with arbitrary refractive index is considered via the Jacobi and Weierstrass elliptic functions in [2], where solitary waves corresponding to optical solitons are recovered for an arbitrary refractive index. Periodic waves are also revealed for the classical case atn = 1 and for the new casen = 1/2. The LPD model is investigated, with the aid of the Ablowitz–Kaup–Newell–Segur (AKNS) hierarchy, to discover novel soliton breathers on the zero background in [3], where the specific conditions for the appearance of the standing and moving soliton breathers are also derived. Additionally, simple formulas for the soliton-breather periods are structured to understand the possibilities for controlling breather dynamics. The LPD equation in polarization-preserving fibers is discussed by virtue of the complete discrimination system for the polynomial method to obtain a complete list of all envelope patterns and to show a variety of dynamical properties of the patterns in [4], where solitons, singularity, periodicity, and double periodicity emerged from the exact parameter condition of existence for each pattern. Several results have been recovered using the CQ–LPD model with its scalar version [5,6,7]. The CQ–LPD equation is adopted using the Lie symmetry analysis in [5], where straddled singular–bright optical solitons are reported using the reduced ordinary differential equations that are handled using Kudryashov’s method. The perturbed CQ–LPD model is addressed by the method of undetermined coefficients in [6], where polarization-preserving fibers and birefringent fibers are studied. The conservation laws for polarization-preserving fibers are also retrieved and enumerated. The existence criteria for the displayed solitons are also presented. The perturbed CQ–LPD model is employed for both with and without polarization in [7], where a full spectrum of soliton solutions are retrieved using the sine-Gordon equation approach. Powerful Lie symmetry analysis has been successfully applied to several models that arise in mathematical physics [1,2,3,4,6,7,8,9,10,11,12,13]. In particular, a successful application of Lie symmetry analysis to the scalar version of the CQ–LPD model, yielded dark and singular solitons [5]. Lie symmetry is one of the most powerful methods for obtaining optical solitons with the governing models in nonlinear optics. In the last few decades, Lie’s method has been described in a number of excellent textbooks and has been applied to a number of physical and engineering models. The study of the group of infinitesimal transformations, in other words, Lie group point transformations, has an important place in this method. The Lie symmetry method, developed in the 19th century (1842–1899) by the Norwegian mathematician Sophus Lie, is exceptionally algorithmic. This method systematically combines and expands famous ad hoc methodologies for constructing optical solitons with model equations in optical fiber communications.

Turning the page, and moving forward, it is about time to address the CQ–LPD model in birefringent fibers when polarization-mode dispersion kicks in. The perturbation terms are also taken into account in the vector CQ–LPD model. Lie symmetry analysis, when applied, leads to coupled ordinary differential equations (ODEs) [14,15]. These ODEs would be subsequently addressed using two powerful integration schemes: Kudryashov’s method and the improved F-expansion scheme. These lead to dark and singular soliton solutions to the model. The details of Lie symmetry analysis followed by the derivation of the soliton solutions from the derived ODEs are detailed and exhibited in the rest of the paper after a succinct introduction to the model.

The Lie symmetry method [5,14,15,16,17] is implemented in the coupled CQ–LPD model with perturbation terms for birefringent fibers (1) and (2) in this section. To solve model (1) and (2), the wave transformations are structured as follows: where the real-valued functions$u=u\left(x,t\right)$ and$w=w\left(x,t\right)$ are the amplitude components of the soliton, while the real-valued function$v=v\left(x,t\right)$ is the phase component of the soliton. Inserting Equation (3) into Equation (1), we have the real and imaginary parts: and

Similarly, inserting Equation (3) into Equation (2), we have the real and imaginary portions, which due to lengthy expressions, we have not written here. For the above system of equations, let us consider a one-parameter ($ϵ$) Lie group of transformations as follows: where$\xi , \tau$ and${\eta }_l$ for$l=1, 2, 3$ are infinitesimals and$ϵ<<1$ is a very small parameter. The associated vector field for Equation (6) is indicated below:

Hence, the fourth prolongation formula [14,15] for Equation (7) is considered as where${\eta }_1^t$,${\eta }_2^t$,${\eta }_3^t$,${\eta }_1^x$,${\eta }_2^x$,${\eta }_3^x$,${\eta }_1^{xxx}$,${\eta }_2^{xxx}$,${\eta }_3^{xxx}$,${\eta }_1^{xxxx}$,${\eta }_2^{xxxx},$ and${\eta }_3^{xxxx}$ represent the extended infinitesimals. By using$p{r}^{\left(4\right)}V(\Delta )=0$ whenever$\Delta =0$ in Equations (4) and (5), and other equations obtained from Equation (2), we obtain the invariance conditions. Plugging the extended infinitesimals into the invariance condition and equating the coefficients of different derivative terms that are equal to zero leads to the following: where$C_l$,$l=1,2,3$ are arbitrary constants. Thus, we obtain only trivial symmetries.