Actually, that splitting is caused by both spin-orbit (fine structure) and spin-spin (hyperfine) tensor effects. Unlike the case with chemical elements, spin-orbit and spin-spin effects are similar in size, because both particles have the same mass.
That's also true of positronium:
Spectroscopy of excited state positronium -- Hyperfine Interactions, Volume 89, Number 1 - SpringerLink
Including
[1109.4476] Observation of the h_b states and Beyond, the bottomonium list becomes:
eta-b(n
1S
0): L=0, S=0, J=0, n=1
upsilon(n
3S
1): L=0, S=1, J=1, n=1,2,3,4,?,?
h-b(n
1P
1): L=1, S=0, J=1, n=1,2
chi-b(n
3P
1,2,3): L=1, S=1, J=0,1,2, n=1,2,3 (1,2 known previously, 3 just discovered)
upsilon(1
3D
2): L=2, S=1, J=2, n=1
The new one's mass: 10.539 +/- 0.004 (stat.) +/- 0.008 (syst.) GeV
I once calculated the sizes of the spin-orbit and spin-spin splittings of charmonium and bottomonium states. For the chi-q states, I used
{J=L-1 value, J=L value, J=L+1 value}
The spin-orbit effect I found as a multiple of <L.S> = (J(J+1) - L(L+1) - S(S+1))/2
{-(L+1), -1, L}
The tensor spin-spin effect I found as a multiple of <(3(n.S)
2 - (S.S))/2>
{-(L+1)/(2(2L-1)), 1/2, -1/((2(2L+3))}
Charmonium:
1S - L=0, n=1: Scalar spin-spin: 116.6 MeV
2S - L=0, n=2: Scalar spin-spin: 49.1 MeV
1P - L=1, n=1: Scalar spin-spin: -0.6 MeV
1
3P - L=1, S=1, n=1: Spin-orbit: 35.0 MeV, Tensor spin-spin: 40.6 MeV
Bottomonium:
1S - L=0, n=1: Scalar spin-spin: 71.4 MeV
1P - L=1, n=1: Scalar spin-spin: 1.6 MeV
2P - L=1, n=2: Scalar spin-spin: 0.5 MeV
1
3P - L=1, S=1, n=1: Spin-orbit: 13.7 MeV, Tensor spin-spin: 13.1 MeV
2
3P - L=1, S=1, n=2: Spin-orbit: 9.3 MeV, Tensor spin-spin: 9.1 MeV
The scalar spin-spin term is sizable only for L = 0, much as for the inverse-square case.
However, the n,L energy levels have a noticeable departure from the inverse-square case.
Charmonium:
2S - 1S: 606 MeV
2S - 1P: 148 MeV
Bottomonium (S = 1 only):
2S - 1S: 563 MeV
2S - 1P: 123 MeV
3S - 2S: 332 MeV
3S - 2P: 95 MeV
2P - 1D: 99 MeV
4S - 3S: 224 MeV
4S - 3P: 40 MeV
The main potential is clearly not quite inverse square.