diff --git a/CHANGELOG.md b/CHANGELOG.md
index 9ca8fce021..210e3e6914 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -6,6 +6,7 @@
 - [[PR#336]](https://github.com/lanl/singularity-eos/pull/336) Included code and documentation for a full, temperature consistent, Mie-Gruneisen EOS based on a pressure power law expansion in eta = 1-V/V0. PowerMG.
 
 ### Fixed (Repair bugs, etc)
+- [[PR343]](https://github.com/lanl/singularity-eos/pull/343) Add chemical potentials to stellar collapse gold files
 - [[PR342]](https://github.com/lanl/singularity-eos/pull/342) Fix missing using statement in stellar collapse root finding routines
 - [[PR341]](https://github.com/lanl/singularity-eos/pull/341) Short-circuit HDF5 machinery when cray-wrappers used in-tree
 - [[PR340]](https://github.com/lanl/singularity-eos/pull/335) Fix in-tree builds with plugin infrastructure
diff --git a/CMakeLists.txt b/CMakeLists.txt
index f7039e5027..4780fcf7ba 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -29,9 +29,9 @@ project(
 
 list(APPEND CMAKE_MODULE_PATH "${PROJECT_SOURCE_DIR}/cmake")
 
-set(SINGULARITY_GOLDFILES_VERSION "goldfiles-1.6.2")
+set(SINGULARITY_GOLDFILES_VERSION "goldfiles-1.8.0")
 set(SINGULARITY_GOLDFILE_HASH
-    33ef74b29937cc1347b525f72662933dd2b0556550f6541f97fc1de8a01c3c2a)
+    249772f3314c4b6b9386aa08895280698ae905ef188f4383b819f28c1484603b)
 
 # ------------------------------------------------------------------------------#
 # Options
diff --git a/utils/scripts/make_tabulated_ideal_sc.py b/utils/scripts/make_tabulated_ideal_sc.py
index 850b93bbe2..dd5b4d7668 100755
--- a/utils/scripts/make_tabulated_ideal_sc.py
+++ b/utils/scripts/make_tabulated_ideal_sc.py
@@ -95,6 +95,11 @@
 Xn = 1. - Xp
 Abar = np.ones_like(PRS)
 Zbar = np.ones_like(PRS)
+mu_e = np.ones_like(PRS)
+mu_n = np.ones_like(PRS)
+mu_p = np.ones_like(PRS)
+muhat = np.ones_like(PRS)
+munu = np.ones_like(PRS)
 
 # Output
 if filename is None:
@@ -123,3 +128,8 @@
   f.create_dataset("Xn", data = Xn)
   f.create_dataset("Abar", data = Abar)
   f.create_dataset("Zbar", data = Zbar)
+  f.create_dataset("mu_e", data = mu_e)
+  f.create_dataset("mu_n", data = mu_n)
+  f.create_dataset("mu_p", data = mu_p)
+  f.create_dataset("muhat", data = muhat)
+  f.create_dataset("munu", data = munu)