A set of novel, grid-robust, higher order vector basis functions is propose
d for the method-of-moments (MoM) solution of integral equations for three-
dimensional (3-D) electromagnetic (EM) problems, These basis functions are
defined over curvilinear triangular patches and represent the unknown elect
ric current density within each patch using the Lagrange interpolation poly
nomials, The highlight of these basis functions is that the Lagrange interp
olation points are chosen to be the same as the nodes of the well-developed
Gaussian quadratures, As a result, the evaluation of the integrals in the
MoM is greatly simplified, Additionally, the surface of an object to he ana
lyzed can be easily meshed because the new basis functions do not require t
he side of a triangular patch to be entirely shared by another triangular p
atch, which is a very stringent requirement for traditional vector basis fu
nctions. The proposed basis functions are implemented with point matching f
or the MoM solution of the electric-field integral equation, the magnetic-f
ield integral equation, and the combined-field integral equation. Numerical
examples are presented to demonstrate the higher order convergence and the
grid robustness for defective meshes using the new basis functions.