Ok, so looking at the formulas for ##\tau_C##, ##\tau_D##, and ##\tau_E## from my previous posts, the exact relationship between them is interesting.
It is evident that, contrary to intuition (at least mine), ##\tau_E > \tau_C##, i.e., the clock at rest at the center of the spherical mass runs faster than the clock in the circular orbit skimming the surface. This is easy to show; we just have to verify that
$$
\frac{3}{2} \sqrt{1 - \frac{2M}{R}} - \frac{1}{2} > \sqrt{1 - \frac{3M}{R}}
$$
which can be done by noting that both expressions go to 1 as ##M / R \rightarrow 0##, and showing that the derivative with respect to ##M / R## of the left-hand expression is always less negative than the derivative with respect to ##M / R## of the right-hand expression. This works out to
$$
- \frac{3}{2} \frac{1}{\sqrt{1 - 2M / R}} > - \frac{3}{2} \frac{1}{\sqrt{1 - 3M / R}}
$$
which confirms that, as ##M / R## gets larger (i.e., as the object gets more compact), ##\tau_E## decreases more slowly than ##\tau_C##, so we will always have ##\tau_E > \tau_C##.
The formula for ##\tau_D## is more complicated, but we can gain insight by looking at its behavior at some useful "test" values: ##R = 3M## and ##R = 6M##. I will also define the dimensionless radial coordinate ##\rho = r / R## to simplify the formulas.
For ##R = 6M##, we have ##\tau_C = 1 / \sqrt{2} \approx 0.7071## and ##\tau_E = \left( \sqrt{6} - 1 \right)/2 \approx 0.7247##. The formula for ##\tau_D## works out to
$$
\tau_D = \sqrt{\frac{7}{4} - \sqrt{\frac{3 - \rho^2}{2}} - \frac{1}{12} \rho^2 - \frac{3 v^2}{3 - \rho^2}}
$$
This is easy to evaluate at the endpoints, ##\rho = 1## (where ##v = 0##) and ##\rho = 0## (where ##v## has its maximum value, which we'll leave as ##v## in the formula). For ##\rho = 1##, we end up with ##\tau_D = \sqrt{2/3} \approx 0.8165## (note that this is equal to ##\tau_B##, the clock rate for the observer at rest on the surface at this radius), and for ##\rho = 0##, we end up with ##\tau_D = \sqrt{7/4 - \sqrt{3/2} - v^2} \approx \sqrt{0.5253 - v^2}##. If we assume that the maximum value of ##v## is approximately the same as the circular orbit velocity for ##R = 6M##, which is ##1/2##, then for ##\rho = 0## we obtain ##\tau_D \approx 0.5246##. This strongly suggests that a full numerical integration would show ##\tau_D < \tau_C < \tau_E## (but see further comments below on the possible behavior of the maximum value of ##v##).
For ##R = 3M##, it is even more interesting. We have ##\tau_C = 0## at this radius, corresponding to the fact that this radius is the photon sphere. We have ##\tau_E = \left( \sqrt{3} - 1 \right) / 2 \approx 0.3660##. The formula for ##\tau_D## becomes
$$
\tau_D = \sqrt{1 - \frac{1}{2} \sqrt{3 - 2 \rho^2} - \frac{1}{6} \rho^2 - \frac{3 v^2}{3 - 2 \rho^2}}
$$
This gives ##\tau_D = \sqrt{1/3} \approx 0.5774## for ##\rho = 1## (again, the same as ##\tau_B## for this radius). For ##\rho = 0## it gives ##\tau_D = \sqrt{1 - \sqrt{3} / 2 - v^2} \approx \sqrt{0.1340 - v^2}##. If ##v \rightarrow 1## for an object this compact, this obviously gives a negative quantity under the square root; but even if ##v## is anything larger than ##\sqrt{1 - \sqrt{3} / 2} \approx 0.3660##, the quantity under the square root will be negative (since ##v^2## is subtracted there). This indicates that, for an object this compact, either the "dropping the clock through the hole" experiment can't actually be done (unlikely since the value of ##\tau_E## shows us that a clock can sit at rest at the center of the object perfectly well), or the maximum velocity achieved in such an experiment, for objects this compact, is much less than the corresponding circular orbit velocity.
If the latter is the case, then we have to re-evaluate the ##R = 6M## case. If the maximum velocity achieved for that case is anything less than ##\sqrt{5/4 - \sqrt{3/2}} \approx 0.1589##, then the minimum value of ##\tau_D## (at ##\rho = 0##) will be greater than ##\tau_C##. Even a somewhat larger maximum velocity will still give an average value of ##\tau_D## greater than ##\tau_C##. This would be nice since we know we must have ##\tau_C < \tau_D## for the ##R = 3M## case (since ##\tau_C = 0## for this case), and we would expect the same inequality between the magnitudes to hold for all values of ##M / R## (at least, it seems to me that that ought to be the case).
That still leaves open the question of the relationship between ##\tau_D## and ##\tau_E##. Given the values above, I am inclined to think that ##\tau_D > \tau_E## holds; however, checking this would involve doing the full numerical integration, which is more than I want to tackle.
So, bottom line, it looks to me like the full set of inequalities is:
$$
\tau_C < \tau_E < \tau_D < \tau_B < \tau_A
$$
But numerical integration would be required to confirm whether ##\tau_D## belongs where it is, or between ##\tau_C## and ##\tau_E## (and to confirm that the same inequalities hold for all allowed values of ##M / R##).