float autoScaleOrbital() {
		float w = (psi_n * (psi_n + 3) - 5f) / psi_Znuc;
		if (w < 1)
			w = 1;
		if (psi_n < 3)
			w += 1;
		float aMax = 0;
		if (!isAnisotropic)
			return w;
		for (int i = 3; --i >= 0;)
			if (anisotropy[i] > aMax)
				aMax = anisotropy[i];
		return w * aMax;
	}