In this post, we check out the scientific power of Python by drawing upon the SciPy library to help use perform complex calculations.
Join the DZone community and get the full member experience.Join For Free
The -norm of a vector is defined to be the th root of the sum of its p th powers.
Such norms occur frequently in applications. Yoshio Koide discovered in 1981 that if you take the masses of the electron, muon, and tau particles, the ratio of the 1 norm to the 1/2 norm is very nearly 2/3. Explicitly,
to at least four decimal places. Since the ratio is tantalizingly close to 2/3, some believe there's something profound going on here and the value is exactly 2/3, but others believe it's just a coincidence.
The value of 2/3 is interesting for two reasons. Obviously it's a small integer ratio. But it's also exactly the midpoint between the smallest and largest possible value. More on that below.
Is the value 2/3 within the measure error of the constants? We'll investigate that with a little Python code.
The masses of particles are available in the
physical_constants dictionary in
scipy.constants. For each constant, the dictionary contains the best estimate of the value, the units, and the uncertainty (standard deviation) in the measurement .
from scipy.constants import physical_constants as pc from scipy.stats import norm def pnorm(v, p): return sum([x**p for x in v])**(1/p) def f(v): return pnorm(v, 1) / pnorm(v, 0.5) m_e = pc["electron mass"] m_m = pc["muon mass"] m_t = pc["tau mass"] v0 = [m_e, m_m, m_t] print(f(v0))
This says that the ratio of the 1 norm and 1/2 norm is 0.666658, slightly less than 2/3. Could the value be exactly 2/3 within the resolution of the measurements? How could we find out?
The function f above is minimized when its arguments are all equal, and maximized when its arguments are maximally unequal. To see this, note that f(1, 1, 1) = 1/3 and f(1, 0, 0) = 1. You can prove that those are indeed the minimum and maximum values. To see if we can make f larger, we want to increase the largest value, the mass of tau, and decrease the others. If we move each value one standard deviation in the desired direction, we get
v1 = [m_e - m_e, m_m - m_m, m_t + m_t] print(f(v1))
which returns 0.6666674, just slightly bigger than 2/3. Since the value can be bigger than 2/3, and less than 2/3, the intermediate value theorem says there are values of the constants within one standard deviation of their mean for which we get exactly 2/3.
Now that we've shown that it's possible to get a value above 2/3, how likely is it? We can do a simulation, assuming each measurement is normally distributed.
N = 1000 above = 0 for _ in range(N): r_e = norm(m_e, m_e).rvs() r_m = norm(m_m, m_m).rvs() r_t = norm(m_t, m_t).rvs() t = f([r_e, r_m, r_t]) if t > 2/3: above += 1 print(above)
When we I ran this, I got 168 values above 2/3 and the rest below. So based solely on our calculations here, not taking into account any other information that may be important, it's plausible that Koide's ratio is exactly 2/3.
 Strictly speaking, we should take the absolute values of the vector components. Since we're talking about masses here, I simplified slightly by assuming the components are non-negative.
Also, what I'm calling the p "norm" is only a norm if p is at least 1. Values of p less than 1 do occur in application, even though the functions they define are not norms. I'm pretty sure I've blogged about such an application, but I haven't found the post.
 The SciPy library obtained its values for the constants and their uncertainties from the CODATA recommended values, published in 2014.
Published at DZone with permission of John Cook , DZone MVB. See the original article here.
Opinions expressed by DZone contributors are their own.