In the present paper a concept of vector nodal finite element has been introduced, algorithms of construction of the vector nodal basis functions with high approximate properties from special functional spaces are presented. Examples of high-order interpolation of harmonic, irrotational vector fields by the developed finite elements illustrate their approximate advantage in comparison with the standard Lagrange elements.