The answer by Martin is good, but I still want to continue along the thought path I had in mind. I hope I can give a different physical basis for confirmation of the number.
I'm sure there are better ways to do this, but I want to do it using only the information I have. I want to consider the transfer of some amount of mass ($m$ for now) from near the surface of the sphere to the inside of the sphere (fully integrating it). In both cases we can asses some amount of energy difference between the spherical blob of mass $M+m$ and the state of the sphere with mass $M$ with the $m$ mass hovering just above the surface. So the two states under consideration are:
- State 1: (M+m) big ball
- State 2: (M) ball next to (m) ball
I have no problem assuming $m \ll M$. Now, I want to write expression for both the gravitational binding energy as well as the surface binding energy of a ball. I'll do this for a generic sphere with a mass of $M$ and uniform density.
$$E_g(M) = - \frac{3 G M^2}{5 R(M)}$$
$$E_s(M) = - 4 \pi R(M)^2 \sigma$$
I'm leaving it in this form because we'll all agree that given the mass and the density, finding $R(M)$ isn't a problem. Now I want to write expressions for the difference in energy from state 2 to state 1. This is straight forward for the surface tension energy because the bodies are non-interacting. However, for the gravitational binding energy, there is still a binding energy between the large ball and the small ball that must be included. Keep in mind that state 1 is the lower energy state.
$$\Delta E_s = E_s(M) + E_s(m) - E_s(M+m)$$
$$\Delta E_g = E_g(M) + E_g(m) - \frac{G M m}{R(M)} - E_g(M+m)$$
Now, obviously, the idea would be to set these equal, assume that $m$ is small, and then find the $M$ as a solution to that equation. But that doesn't work! I think I have a major conceptual flaw in this approach, where the scaling of the surface area of the small $m$ blob just doesn't follow a scaling that works. I didn't know what to do, so I just removed the $E_s(m)$, abandoning all logical reasoning behind my work. But when I did this and used Martin's values, I obtained the following:
$$M=1.05 \times 10^6 kg$$
This was assuming $m=0.1 kg$. If I change that value it doesn't change $M$ very much, which is encouraging. This is a satisfying answer for me, because Martin's answer comes out to around 0.5 million kg.