This work presents a theoretical approach to a highly multimode-guided wave-based fourth harmonic (FH) generation method by exploiting a zero-order Bessel–Gauss beam through the total internal reflection quasi phase-matching technique. An anisotropic crystal featuring a rectangular structure of beta barium borate has been used to carry out this analysis. In contrast to the conventional frequency conversion scheme using a Gaussian beam, the projected scheme yields a maximum six times superior efficiency of 6.9547% in a crystal length of nearly 6.3 mm, corresponding to a 266-nm FH wavelength. The analysis has also taken into account optical losses, such as absorption, Goos–Hänchen-shift, surface roughness, and the effect of the nonlinear law of reflection. Finally, a systematic study has been performed to describe the influence of parametric modification, such as variations in slab length and thickness, to provide a more comprehensive picture of the overall analysis.