diff --git a/notebooks/semi-analytic-models.ipynb b/notebooks/semi-analytic-models.ipynb
index 95d54de1..a93b6acd 100644
--- a/notebooks/semi-analytic-models.ipynb
+++ b/notebooks/semi-analytic-models.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "markdown",
-   "id": "8d47f1a3",
+   "id": "0",
    "metadata": {},
    "source": [
     "# `holodeck` - Semi-Analytic Models"
@@ -10,7 +10,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "af59b140",
+   "id": "1",
    "metadata": {},
    "source": [
     "For more information on the holodeck SAMs, see the [holodeck getting started guide](https://holodeck-gw.readthedocs.io/en/main/getting_started/index.html)."
@@ -19,7 +19,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "f9a9d9bc",
+   "id": "2",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -41,7 +41,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "630dd2a1",
+   "id": "3",
    "metadata": {},
    "source": [
     "## Quick Start"
@@ -50,7 +50,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "f691ee46",
+   "id": "4",
    "metadata": {},
    "source": [
     "Construct a Semi-Analytic Model (SAM) using all of the default components"
@@ -59,7 +59,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "4f80bc42",
+   "id": "5",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -71,7 +71,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "5bb4f96a",
+   "id": "6",
    "metadata": {},
    "source": [
     "Choose the edges of the frequency bins at which to calculate the GWB"
@@ -80,7 +80,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "a4e14d8e",
+   "id": "7",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -95,7 +95,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "34f11e96",
+   "id": "8",
    "metadata": {},
    "source": [
     "Calculate GWB from this SAM.  \n",
@@ -107,7 +107,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "abde3ae2",
+   "id": "9",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -119,7 +119,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "ebc65eca",
+   "id": "10",
    "metadata": {},
    "source": [
     "Plot GWB over multiple realizations"
@@ -128,7 +128,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "d836901c",
+   "id": "11",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -139,7 +139,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "0e256c02",
+   "id": "12",
    "metadata": {},
    "source": [
     "Slightly fancier plot:"
@@ -148,7 +148,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b69aed2a",
+   "id": "13",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -185,7 +185,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "4356677c",
+   "id": "14",
    "metadata": {},
    "source": [
     "## Specifics"
@@ -193,7 +193,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "3dbf39b8",
+   "id": "15",
    "metadata": {},
    "source": [
     "### Constructing a SAM"
@@ -201,7 +201,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "d736718a",
+   "id": "16",
    "metadata": {},
    "source": [
     "SAMs are built from simple analytic models to derive the number-density of MBH binaries.\n",
@@ -222,7 +222,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "8e8f3189",
+   "id": "17",
    "metadata": {},
    "source": [
     "The SAMs are initialized over a 3-dimensional parameter space of total MBH mass ($M = m_1 + m_2$), MBH mass ratio ($q = m_2 / m_1 \\leq 1$), and redshift ($z$).  The `holodeck` code typically refers to the number of bins in each of these dimensions as `M`, `Q`, and `Z`; for example, the shape of the number-density of galaxy mergers will be `(M, Q, Z)`."
@@ -231,7 +231,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "7dacb18d",
+   "id": "18",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -261,7 +261,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "7a971041",
+   "id": "19",
    "metadata": {},
    "source": [
     "### number density and the SAM grid"
@@ -269,7 +269,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "ddd1b513",
+   "id": "20",
    "metadata": {},
    "source": [
     "The formation rate of MBH-MBH 'binaries' is calculated in `Semi_Analytic_Model.static_binary_density`, evaluated at the edges of the grid so that it's shape is the number of bins in each dimension, plus one, i.e. `(M+1, Q+1, Z+1)`.  `static_binary_density` is implemented as a `@property` so that the first time the value is accessed it is calculated and cached, and then returned immediately on subsequent calls. \n",
@@ -280,7 +280,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "3cefecee",
+   "id": "21",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -301,7 +301,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "4689a2eb",
+   "id": "22",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -335,7 +335,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "dc7b2575",
+   "id": "23",
    "metadata": {},
    "source": [
     "### total number of binaries in a universe"
@@ -343,7 +343,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "6a9d679e",
+   "id": "24",
    "metadata": {},
    "source": [
     "Above, we calculated the volumetric number-density rate of binary mergers.  Here, we calculate the total number of binaries in a simulated universe at particular GW frequencies of interest.  The SAM models currently assume circular binary orbits, so that the GW frequency is exactly twice the orbital frequency.\n",
@@ -366,7 +366,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "31caefdd",
+   "id": "25",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -381,7 +381,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "875f617c",
+   "id": "26",
    "metadata": {},
    "source": [
     "For each point in the 3-dimensional SAM grid, we will be calculating the number of binaries at each frequency.  So the returned values will be 4-dimensional with an additional axis with `F` frequency bins added: `(M, Q, Z, F)`."
@@ -389,7 +389,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "894b20da",
+   "id": "27",
    "metadata": {},
    "source": [
     "**GW-Only Binary Evolution**"
@@ -398,7 +398,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "d0cd03af",
+   "id": "28",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -417,7 +417,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "63496676",
+   "id": "29",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -442,7 +442,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "38215e23",
+   "id": "30",
    "metadata": {},
    "source": [
     "Convert from differential number of binaries to actual number of binaries"
@@ -451,7 +451,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "d3c9d78f",
+   "id": "31",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -467,7 +467,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "79774740",
+   "id": "32",
    "metadata": {},
    "source": [
     "Calculate the total number of binaries in certain ranges of parameter space"
@@ -476,7 +476,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "986c695b",
+   "id": "33",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -507,7 +507,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "0cbf07d3",
+   "id": "34",
    "metadata": {},
    "source": [
     "### Self-Consistent Binary Evolution (Phenomenological Model)"
@@ -516,7 +516,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "d704a3e7",
+   "id": "35",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -532,7 +532,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "3212639a",
+   "id": "36",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -566,7 +566,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "1cb3479a",
+   "id": "37",
    "metadata": {},
    "source": [
     "## Gravitational Waves"
@@ -574,7 +574,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "857fab5c",
+   "id": "38",
    "metadata": {},
    "source": [
     "### Compare GWB and CW"
@@ -583,7 +583,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b8098cc1",
+   "id": "39",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -604,7 +604,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "b224fa4c",
+   "id": "40",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -628,7 +628,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "526f1ed3",
+   "id": "41",
    "metadata": {},
    "source": [
     "Calculate the distribution of GWB Amplitudes at 1/yr"
@@ -637,7 +637,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "e9fc9e99",
+   "id": "42",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -680,7 +680,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "4b567f98",
+   "id": "43",
    "metadata": {},
    "source": [
     "## Plot GWB Amplitude Distribution vs. M-MBulge parameters"
@@ -689,7 +689,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "20f4383b",
+   "id": "44",
    "metadata": {},
    "source": [
     "Calculate GWB amplitudes at $f = 1/yr$ over a grid of M-Mbulge parameters, specifically the amplitude and power-law."
@@ -698,7 +698,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "55cb90ca",
+   "id": "45",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -727,7 +727,7 @@
   {
    "attachments": {},
    "cell_type": "markdown",
-   "id": "cc13e780",
+   "id": "46",
    "metadata": {},
    "source": [
     "Plot the interquartile ranges for each power-law, as a function of normalization"
@@ -736,7 +736,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "22a7c9c1",
+   "id": "47",
    "metadata": {},
    "outputs": [],
    "source": [
@@ -774,7 +774,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.4"
+   "version": "3.11.7"
   },
   "toc": {
    "base_numbering": 1,