Thermoelectric generation systems (TEGSs) are used to convert temperature difference and heat flow into DC power based on the Seebeck theorem. The basic unit of TEGS is the thermoelectric module (TEM). TEGSs have gained increasing interest in the research fields of sustainable energy. The output power from TEM is mostly reliant on differential temperature between the hot and cold sides of the TEM added to the value of the load. As such, a robust MPPT strategy (MPPTS) is required to ensure that the TEGS is operating near to the MPP while varying the operating conditions. Two main drawbacks may occur in the conventional MPPTSs: low dynamic response, such as in the incremental resistance (INR) method, and oscillations around MPP at steady state, such as in the hill climbing (HC) method. In the current research work, an optimized fractional MPPTS is developed to improve the tracking performance of the TEGS, and remove the two drawbacks of the conventional MPPTSs. The proposed strategy is based on fractional order control (FOC). The main advantage of FOC is that it offers extra flexible time and frequency responses of the control system consent for better and robust performance. The optimal parameters of the optimized fractional MPPTS are identified by a manta ray foraging optimization (MRFO). To verify the robustness of the MRFO, the obtained results are compared with ten other algorithms: particle swarm optimization; whale optimization algorithm; Harris hawks optimization; heap-based optimizer; gradient-based optimizer; grey wolf optimizer; slime mould algorithm; genetic algorithm; seagull optimization algorithm (SOA); and tunicate swarm algorithm. The maximum average cost function of 4.92934 kWh has been achieved by MRFO, followed by SOA (4.5721 kWh). The lowest STD of 0.04867 was also accomplished by MRFO. The maximum efficiency of 99.46% has been obtained by MRFO, whereas the lowest efficiency of 74.01% was obtained by GA. Finally, the main findings proved the superiority of optimized fractional MPPTS compared with conventional methods for both steady-state and dynamic responses.