  We study symmetry of solutions of the mean field equation \[ \Delta u
+\rho(\frac{Ke^u}{\int_{T_\epsilon} Ke^u} -\frac{1}{|T_\epsilon|} )=0\] on the
flat torus $T_\epsilon=[-\frac{1}{2\epsilon}, \frac{1}{2\epsilon}] \times
[-\frac{1}{2}, \frac{1}{2}]$ with $0<\epsilon \leq 1$, where $K\in
C^2({T}_\epsilon)$ is a positive function with $-\Delta \ln K \leq
\frac{\rho}{|T_\epsilon|}$ and $\rho \leq 8\pi$. We prove that if $(x_0,y_0)$
is a critical point of the function $u+ln(K)$, then $u$ is evenly symmetric
about the lines $x=x_0$ and $y=y_0$, provided $K$ is evenly symmetric about
these lines. In particular we show that all solutions are one-dimensional if
$K\equiv 1$ and $\rho \leq 8\pi$. The results are sharp and answer a conjecture
of Lin and Lucia affirmatively. We also prove some symmetry results for mean
field equations on annulus.

Analysis of PDEs