Obligatory Covid-19 blogging

SARS-CoV-2 spike ectodomain structure (open state)
SARS-CoV-2 spike ectodomain structure (open state) https://3dprint.nih.gov/discover/3DPX-013160
Over at Practical Ethics I have blogged a bit:

The Unilateralist Curse and Covid-19, or Why You Should Stay Home: why we are in a unilateralist curse situation in regards to staying home, making it rational to stay home even when it seems irrational.

Taleb and Norman had a short letter Ethics of Precaution: Individual and Systemic Risk making a similar point, noting that recognizing the situation type and taking contagion dynamics into account is a reason to be more cautious. It is different from our version in the effect of individual action: we had single actor causing full consequences, the letter has exponential scale-up. Also, way more actors: everyone rather than just epistemic peers, and incentives that are not aligned since actors do not bear the full costs of their actions. The problem is finding strategies robust to stupid, selfish actors. Institutional channeling of collective rationality and coordination is likely the only way for robustness here.

Never again – will we make Covid-19 a warning shot or a dud? deals with the fact that we are surprisingly good at forgetting harsh lessons (1918, 1962, Y2K, 2006…), giving us a moral duty to try to ensure appropriate collective memory of what needs to be recalled.

This is why Our World In Data, the Oxford COVID-19 Government Response Tracker and IMF’s policy responses to Covid-19 are so important. The disjointed international responses act as natural experiments that will tell us important things about best practices and the strengths/weaknesses of various strategies.

When the inverse square stops working

In physics inverse square forces are among the most reliable things. You can trust that electric and gravitational fields from monopole charges decay like 1/r^2. Sure, dipoles and multipoles may add higher order terms, and extended conductors like wires and planes produce other behaviour. But most of us think we can trust the 1/r^2 behaviour for spherical objects.

I was surprised to learn that this is not at all true recently when a question at the Astronomy Stack Exchange asked about whether gravity changes near the surface of dense objects.

Electromagnetism does not quite obey the inverse square law

The cause was this paper by John Lekner, that showed that there can be attraction between conductive spheres even when they have the same charge! (Popular summary in Nature by Philip Ball) The “trick” here is that when the charged spheres approach each other the charges on the surface redistribute themselves which leads to a polarization. Near the the other sphere charges are pushed away, and if one sphere has a different radius from the other the “image charge” can be opposite, leading to a net attraction.

Charge distribution on two spherical conductors with the same net charge.
Charge distribution on two spherical conductors with the same net charge.

The formulas in the paper are fairly non-intuitive, so I decided to make an approximate numeric model. I put 500 charges on two spheres (radius 1.0 and 2.0) and calculated the mutual electrostatic repulsion/attraction moving them along the surface. Iterate until it stabilizes, then calculate the overall force on one of the spheres.

Force between two equally charged (blue) and two oppositely charged (red) spheres, and force time squared distance.
Force between two equally charged (blue) and two oppositely charged (red) spheres, and force time squared distance.

The result is indeed that the 1/r^2 law fails as the spheres approach each other. The force times squared distance is constant until they get within a few radii, and then the equally charged sphere begins to experience less repulsion and the oppositely charged spheres more attraction than expected. My numerical method is too sloppy to do a really good job of modelling the near-touching phenomena, but it is enough to show that that the inverse square effect is not true for conductors close enough.

Gravity doesn’t obey the inverse square law either

My answer to the question was along the same lines: if two spherical bodies get close to each other they are going to deform, and this will also affect the force between them. In this case it is not charges moving on the surface, but rather gravitational and tidal distortion turning them into ellipsoids. Strictly speaking, they will turn into more general teardrop shapes, but we can use the ellipsoid as an approximation. If they have fixed centres of mass they will be prolate ellipsoids, while if they are orbiting each other they would be general three-axis ellipsoids.

Calculating the gravitational field of an ellipsoid has been done and has a somewhat elegant answer that unfortunately needs to be expressed terms of special functions. The gravitational potential in the system is just the sum of the potentials from both ellipsoids. The equilibrium shapes would correspond to ellipsoids with the same potential along their entire surface; maybe there is an explicit solution, but it does look likely to be an algebraic mess of special functions.

I did a numeric exploration instead. To find the shape I started with spheres and adjusted the semi-major axis (while preserving volume) so the potential along the surface became more equal at the poles. After a few iterations this gives a self-consistent shape. Then I calculated the force (the derivative of the potential) due to this shape on the other mass.

Semi-major axis, force, and force times distance squared for two self-gravitating unit volume ellipsoids at different center-of-mass distances.
Semi-major axis, force, and force times distance squared for two self-gravitating unit volume ellipsoids at different center-of-mass distances.

The result is indeed that the force increases faster than 1/r^2 as the bodies approach each other, since they elongate and eventually merge (a bit before this they will deviate from my ellipsoidal assumption).

This was the Newtonian case. General relativity is even messier. In intense gravitational fields space-time is curved and expanded, making even the meaning of the distance in the inverse square law problematic. For black holes the Paczyński–Wiita potential is an approximation of the potential, and it deviates from the U(r)=-GM/r potential as U_{PW}(r)=-GM/(r-R_S) (where R_S is the Schwarzschild radius). It makes the force increase faster than the classical potential as we approach r=R_S.

Normally we assume that charges and masses stay where they are supposed to be, just as we prefer to reason as if objects are perfectly rigid or well described by point masses. In many situations this stops being true and then the effective forces can shift as the objects and their charges shift around.