beta0 = -k0 * np.sin(theta0) * np.sin(phi0)