Fix lithium anomaly base-10 harmonic tuning

.gitignore

@@ -0,0 +1 @@
+venv/

Copy_of_QAG_AVI.ipynb

@@ -79,11 +79,11 @@
         "\n",
         "plt.figure(figsize=(12,5))\n",
         "plt.subplot(131)\n",
-        "plt.plot(r_train, kappa_obs_train, 'k--', label='Obs κ (HST)')\n",
+        "plt.plot(r_train, kappa_obs_train, 'k--', label='Obs \u03ba (HST)')\n",
         "plt.plot(r_train, gas_obs, 'r-', label='Gas (Chandra)')\n",
         "plt.plot(r_train, kappa_qag_train, 'b-', lw=2, label='QAG AVI')\n",
         "plt.axvspan(150, 350, alpha=0.3, color='gray')\n",
-        "plt.ylabel('κ')\n",
+        "plt.ylabel('\u03ba')\n",
         "plt.legend()\n",
         "plt.grid(alpha=0.3)\n",
         "plt.title('Train Wreck Abell 520')\n",
@@ -94,7 +94,7 @@
         "z = np.logspace(6, 10, 100)\n",
         "Li_LCDM = 4.8e-10 * np.exp(-z/1e9)\n",
         "# The lithium anomaly fix using base-10 harmonic tuning\n",
-        "Li_QAG = Li_LCDM * 0.35 * (1 + np.tanh((np.log10(z)-8.5)/0.5))\n",
+        "Li_QAG = Li_LCDM * 0.35 * (1 + np.sin((np.log10(z)-8.5)/0.5))\n",
         "plt.loglog(z, Li_LCDM, 'k-', label='BBN Std')\n",
         "plt.loglog(z, Li_QAG, 'r--', label='QAG')\n",
         "plt.scatter(10**8.5, 1.6e-10, s=100, color='gold', label='Obs')\n",
@@ -109,11 +109,11 @@
         "C_l_LCDM = 5500 / l * (1 + 0.1 * np.sin(l / 220 * np.pi))**2\n",
         "# CMB injected with harmonic resonance\n",
         "C_l_QAG = C_l_LCDM * (1 + 0.05 * np.sin(np.log(l+10) * 1.1))\n",
-        "plt.semilogy(l, C_l_LCDM, 'k-', label='Planck ΛCDM')\n",
+        "plt.semilogy(l, C_l_LCDM, 'k-', label='Planck \u039bCDM')\n",
         "plt.semilogy(l, C_l_QAG, 'r--', lw=2, label='QAG Resonance')\n",
         "plt.axvline(220, color='g', ls=':')\n",
         "plt.xlabel('l')\n",
-        "plt.ylabel('C_l μK²')\n",
+        "plt.ylabel('C_l \u03bcK\u00b2')\n",
         "plt.legend()\n",
         "plt.title('CMB')\n",
         "plt.grid()\n",
@@ -131,7 +131,7 @@
         "plt.axvline(250, color='orange', ls=':', label='Gas Peak')\n",
         "plt.title('Bullet: QAG 10:1 Offset')\n",
         "plt.xlabel('kpc')\n",
-        "plt.ylabel('κ')\n",
+        "plt.ylabel('\u03ba')\n",
         "plt.legend()\n",
         "plt.grid()\n",
         "plt.show()"
@@ -169,4 +169,4 @@
       ]
     }
   ]
-}
+}