The energy of a graph was introduced by Gutman in 1978 as the sum of the absolute values of the eigenvalues of its adjacency matrix. We study the energy of integral circulant graphs, also called gcd graphs, which can be characterized by their vertex count n and a set D of divisors of n in such a way that they have vertex set Z/nZ and edge set {{a, b} : a, b ∈ Z/nZ, gcd(a − b, n) ∈ D}.Given an arbitrary prime power p s , we determine all divisor sets maximising the energy of an integral circulant graph of order p s . This enables us to compute the maximal energy E max (p s ) among all integral circulant graphs of order p s .