The present research explores nematicons in liquid crystals (LCs) with quadruple power law nonlinearity utilizing the modified extended Fan sub-equation technique as an analytical tool to investigate the optical spatial soliton solutions. For the inaugural time, a novel version of nonlinearity is investigated in relation to LCs. There are distinct applications for the several wave solutions that have been created in optical handling data. The aforementioned modified extended Fan subequation approach offers novel, comprehensive solutions that are relatively easy to deploy in comparison to earlier, regular methodologies. This approach translates a coupled non-linear partial differential equation into a coupled ordinary differential equation through implementing a traveling wave conversion. This approach indicates that a large variety of traveling and solitary solutions that rely upon five parameters can be incorporated by the nematicons in LCs. In addition, the investigation yields solutions of the single and mixed non-degenerate Jacobi elliptic function form. Novel solutions, such as the periodic pattern, kink and anti-kink patterns, N-pattern, W-pattern, anti-Z-pattern, M-pattern, V-pattern, complexion pattern and anti-bell pattern, or dark soliton solutions of nematic LCs, have been constructed by means of modified extended Fan subequation technique through granting suitable values for the parameters. The computer software Mathematica 14 is used to illustrate several modulus, real and imaginary solutions visually in the form of contour, 2D, and 3D visualizations that help understand the concrete importance of the nematicons in LCs. This research additionally offers a physical comprehension of the obtained solutions and applications of model. The imposed approach is ultimately thought to be more potent and effective than alternative approaches, and the solutions found in this work could be beneficial in our understanding of soliton structures in LCs.