I read an article saying that the inner product between divergence-free current and a gradient field is zero.
A divergence-free surface current is $\nabla\cdot\vec{J}=0$, and $\vec{J}$ could be represented as $\vec{J}=\nabla\times(\psi\hat{n})$, where $\hat{n}$ is the normal vector of the surface. So the statement becomes: $\left( \nabla\times(\psi\hat{n}) \right) \cdot \nabla \varphi=0$.
I think according to the identity: $$\nabla\cdot(\vec{A}\times\vec{B})=\vec{B}\cdot(\nabla\times\vec{A})-\vec{A}\cdot(\nabla\times\vec{B})$$ we have $$\nabla\times(\psi\hat{n})\cdot \nabla \varphi=\nabla\cdot(\psi\hat{n}\times\nabla\varphi)+\psi\hat{n}\cdot\nabla\times\nabla\varphi=\nabla\cdot(\psi\hat{n}\times\nabla\varphi),$$ but what next?
Update Thank you Luboš Motl. I suppose I now understand why, but I don't have enough reputation to reply below, so just update here my answer.
The goal is to prove $\int_s \vec{J}\cdot\nabla\varphi ds=0$ The whole process is as follows:
First, $\vec{J}$cannot go across the surface edge, so $\vec{J}\cdot\hat{t}=0$, where $\hat{l}$ is the surface edge direction and $\hat{t}=\hat{l}\times\hat{n}$ is the edge out direction.
Second, according to the identity $$\nabla\cdot(\vec{A}\times\vec{B})=\vec{B}\cdot(\nabla\times\vec{A})-\vec{A}\cdot(\nabla\times\vec{B}) \, ,$$ we have \begin{align} \vec{J}\cdot\nabla\varphi &=\nabla\times(\psi\hat{n})\cdot \nabla \varphi \\ &=\nabla\cdot(\psi\hat{n}\times\nabla\varphi)+\psi\hat{n}\cdot\nabla\times\nabla\varphi \\ &=\nabla\cdot(\psi\hat{n}\times\nabla\varphi) \end{align} since $$\nabla\times(f\vec{A})=\nabla{f}\times\vec{A}+f(\nabla\times A)$$ $$\psi\hat{n}\times\nabla\varphi= -\nabla \times (\varphi\psi\hat{n}) + \varphi \nabla \times(\psi\hat{n}) \, .$$ Then \begin{align} \nabla\cdot(\psi\hat{n}\times\nabla\varphi) &=\nabla\cdot(-\nabla\times(\varphi\psi\hat{n})+\varphi\nabla\times(\psi\hat{n})) \\ &=\nabla\cdot(\varphi\nabla\times(\psi\hat{n})) \end{align} Finally, \begin{align} \int_s \vec{J}\cdot\nabla\varphi ds &=\int_s\nabla\times(\psi\hat{n})\cdot \nabla \varphi ds \\ &= \int_s \nabla\cdot(\varphi\nabla\times(\psi\hat{n}))ds \\ &=\oint_l \varphi\nabla\times(\psi\hat{n})\cdot\hat{t}dl \\ &=\oint_l \varphi\vec{J}\cdot\hat{t}dl \\ &=0 \, . \end{align}
I think here the important things are:
Generally speaking, divergence-free current usually can be expressed as $\vec{J}=\nabla\times\vec{T}$, and $\vec{J}=\nabla\times(\psi\hat{n})$ is specially for surface current.
the $\hat{n}$ is only valid on the surface(there is no meaning of $\hat{n}$ for point in side of a body). the integral is on the surface rather than on the body. According to the original article, it is just talking about PEC and surface current.