Modeling and design of high frequency electronic systems such as antennas and microwave devices require the rigorous numerical solution of Maxwell's equations. The frequency-domain (time-harmonic) tangential vector nite element method (TVFEM) for Maxwell equations results in a second-order dynamical electromagnetic model that must be repeatedly solved for multiple frequencies, excitation or material parameters each design loop. This leads to extremely long design turnaround that often is not optimal. This work will propose an accurate, error controllable and ef cient multi-parametric model order reduction scheme that signi cantly accelerate these parameters sweep. At the core of this work is the proper orthogonal decomposition (POD) sampling technique and balanced truncation (BT) algorithm that are used to reduce multi-parameter spaces that include frequency, material parameters and in nite array scan angles. The proposed methodology employs a novel computational scheme based on adaptive POD sampling and the singular value decomposition of the low-rank Hankel matrix. Numerical examples con rm the signi cant time savings and good accuracy of the method."