From de889f688d37b16b3dd63ef38a3aa48c07b262c4 Mon Sep 17 00:00:00 2001 From: Luke Zoltan Kelley Date: Wed, 1 May 2024 12:36:56 -0700 Subject: [PATCH] Strip notebook --- notebooks/semi-analytic-models.ipynb | 98 ++++++++++++++-------------- 1 file changed, 49 insertions(+), 49 deletions(-) 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,