{
"cells": [
{
"cell_type": "markdown",
"id": "086f9d1b",
"metadata": {},
"source": [
"Copyright 2025 Marc Kamionkowski and Kris Sigurdson.\n",
"Licensed under Apache-2.0 (http://www.apache.org/licenses/LICENSE-2.0).\n",
"\n",
"This notebook provides some simple python scripts to make plots and movies from the output of NSphere-SIDM. \n",
"\n",
"The scripts use data produced by `nsphere --sidm --ntimesteps 12500000 --dtwrite 25000 --nout 250 --tfinal 1100 --sort 3` (the parameters used for the plots in the paper). This is then run with 100000 particles, the default halo profile, an NFW profile with scale radius $r_s = 1.18$ kpc, mass $= 1.15 \\times 10^{-9}$ M$_\\odot$, and a large-radius cutoff function $[1+(r/r_{\\text{vir}})^{10}]^{-1}$ with $r_{\\text{vir}} = 19$ (concentration parameter). The SIDM cross section (per unit mass) is isotropic and velocity-independent, and taken to be 50 cm$^2$/g. All of these can be changed using command line options, run `nsphere --help` for more information.\n",
"\n",
"The jupyter notebook is set up to be run from the NSphere-SIDM directory. It was run with Python 3.12.4"
]
},
{
"cell_type": "markdown",
"id": "e55c1626",
"metadata": {},
"source": [
"## SIDM Core Collapse Analysis, Visualization and Plots"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "bbf95eff",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Using 154 x 153 bins\n",
"Histogram data shape: (251, 154, 153)\n"
]
},
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"