In this thesis we consider harmonic maps and barycentric maps in the context of higher Teichmüller theory. We are particularly interested in how these maps can be used to study Hitchin representations. The main results of this work are as follows. Our first result states that equivariant harmonic maps into non-compact symmetric spaces that satisfy suitable non-degeneracy conditions depend in a real analytic fashion on the metric of the domain manifold and the representations they are associated to. For our second result we consider the energy functional on Teichmüller space that is associated to a Hitchin representation. We prove that this functional is strictly plurisubharmonic for Hitchin representations into either PSL(n, R), PSp(2n, R),...