  Using data taken at 23 center-of-mass energies between 4.0 and 4.6 GeV with
the BESIII detector at the BEPCII collider and with a total integrated
luminosity of approximately $15$ fb$^{-1}$, the process $e^+e^-\to 2(p
\bar{p})$ is studied for the first time. The Born cross sections for $e^+e^-\to
2(p \bar{p})$ are measured, and no significant structure is observed in the
lineshape. The baryon pair ($pp$ and $\bar{p}\bar{p}$) invariant mass spectra
are consistent with phase space, therefore no hexaquark or di-baryon state is
found.
