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.