Skip to content

Instantly share code, notes, and snippets.

@syrte

syrte/FFB_equations.ipynb

Created November 12, 2023 07:29
Show Gist options
  • Save syrte/d92b7bd63c83d6943bca95ac621f863d to your computer and use it in GitHub Desktop.
Save syrte/d92b7bd63c83d6943bca95ac621f863d to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
FFB_equations.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "47a41ec9",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 1;\n",
" var nbb_unformatted_code = \"# code format\\n%load_ext nb_black\";\n",
" var nbb_formatted_code = \"# code format\\n%load_ext nb_black\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# code format\n",
"%load_ext nb_black"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "662d6b2c",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 2;\n",
" var nbb_unformatted_code = \"from astropy.units import *\\nfrom astropy import units as u\\nfrom astropy import constants as c\";\n",
" var nbb_formatted_code = \"from astropy.units import *\\nfrom astropy import units as u\\nfrom astropy import constants as c\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from astropy.units import *\n",
"from astropy import units as u\n",
"from astropy import constants as c"
]
},
{
"cell_type": "markdown",
"id": "1b653997",
"metadata": {},
"source": [
"# common variables\n",
"\n",
"define commonly used variables, which serve as ireducible units untill substitution"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "a9985af2",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 3;\n",
" var nbb_unformatted_code = \"# scaled variables\\nMv_sc = u.def_unit(\\\"[Mv/10^10.8Msun]\\\", format={\\\"latex\\\": r\\\"M_{v,10.8}\\\"})\\nzp1_sc = u.def_unit(\\\"[(1+z)/10]\\\", format={\\\"latex\\\": r\\\"(1+z)_{10}\\\"})\\nSFR_sc = u.def_unit(\\\"[SFR/(Msun/yr)]\\\", format={\\\"latex\\\": r\\\"SFR_{M_\\\\odot/yr}\\\"})\\n\\nlambdas_sc = u.def_unit(\\\"lambda_s025\\\", format={\\\"latex\\\": r\\\"\\\\lambda_{s,.025}\\\"})\\nlambdad_sc = u.def_unit(\\\"lambda_d025\\\", format={\\\"latex\\\": r\\\"\\\\lambda_{d,.025}\\\"})\\n\\n# dimentionless numbers\\neta = u.def_unit(\\\"eta\\\", format={\\\"latex\\\": r\\\"\\\\eta\\\"})\\neps = u.def_unit(\\\"eps\\\", format={\\\"latex\\\": r\\\"\\\\epsilon\\\"})\\n\\n# bases\\nbases = [Mv_sc, zp1_sc, SFR_sc, lambdas_sc, lambdad_sc, eta, eps]\";\n",
" var nbb_formatted_code = \"# scaled variables\\nMv_sc = u.def_unit(\\\"[Mv/10^10.8Msun]\\\", format={\\\"latex\\\": r\\\"M_{v,10.8}\\\"})\\nzp1_sc = u.def_unit(\\\"[(1+z)/10]\\\", format={\\\"latex\\\": r\\\"(1+z)_{10}\\\"})\\nSFR_sc = u.def_unit(\\\"[SFR/(Msun/yr)]\\\", format={\\\"latex\\\": r\\\"SFR_{M_\\\\odot/yr}\\\"})\\n\\nlambdas_sc = u.def_unit(\\\"lambda_s025\\\", format={\\\"latex\\\": r\\\"\\\\lambda_{s,.025}\\\"})\\nlambdad_sc = u.def_unit(\\\"lambda_d025\\\", format={\\\"latex\\\": r\\\"\\\\lambda_{d,.025}\\\"})\\n\\n# dimentionless numbers\\neta = u.def_unit(\\\"eta\\\", format={\\\"latex\\\": r\\\"\\\\eta\\\"})\\neps = u.def_unit(\\\"eps\\\", format={\\\"latex\\\": r\\\"\\\\epsilon\\\"})\\n\\n# bases\\nbases = [Mv_sc, zp1_sc, SFR_sc, lambdas_sc, lambdad_sc, eta, eps]\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# scaled variables\n",
"Mv_sc = u.def_unit(\"[Mv/10^10.8Msun]\", format={\"latex\": r\"M_{v,10.8}\"})\n",
"zp1_sc = u.def_unit(\"[(1+z)/10]\", format={\"latex\": r\"(1+z)_{10}\"})\n",
"SFR_sc = u.def_unit(\"[SFR/(Msun/yr)]\", format={\"latex\": r\"SFR_{M_\\odot/yr}\"})\n",
"\n",
"lambdas_sc = u.def_unit(\"lambda_s025\", format={\"latex\": r\"\\lambda_{s,.025}\"})\n",
"lambdad_sc = u.def_unit(\"lambda_d025\", format={\"latex\": r\"\\lambda_{d,.025}\"})\n",
"\n",
"# dimentionless numbers\n",
"eta = u.def_unit(\"eta\", format={\"latex\": r\"\\eta\"})\n",
"eps = u.def_unit(\"eps\", format={\"latex\": r\"\\epsilon\"})\n",
"\n",
"# bases\n",
"bases = [Mv_sc, zp1_sc, SFR_sc, lambdas_sc, lambdad_sc, eta, eps]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "5f5bb6c6",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 4;\n",
" var nbb_unformatted_code = \"def subs_unit(quantity, substitute_dict):\\n \\\"substitute unit by replacement\\\"\\n \\\"\\\"\\\"\\n quantity : astropy.units.Quantity\\n \\n substitute_dict : dict\\n Dict for substitution, e.g. {km:100*m, Gyr:1e9*yr}\\n \\\"\\\"\\\"\\n # iterative substitution for 10 times or untill all items repaced\\n for i in range(10):\\n count = 0\\n for unit, replace in substitute_dict.items():\\n unit_dict = dict(zip(quantity.unit._bases, quantity.unit._powers))\\n if unit in unit_dict:\\n power = unit_dict[unit]\\n if power != 0:\\n quantity = quantity * (replace / unit) ** power\\n count += 1\\n if count == 0:\\n break\\n return quantity\\n\\n\\ndef convert_ffb_threhold(quantity):\\n return subs_unit(quantity, {Mv_sc: zp1_sc**-6.2000001})\\n # use 6.2000001 instead of 6.2 to avoid displying fraction\";\n",
" var nbb_formatted_code = \"def subs_unit(quantity, substitute_dict):\\n \\\"substitute unit by replacement\\\"\\n \\\"\\\"\\\"\\n quantity : astropy.units.Quantity\\n \\n substitute_dict : dict\\n Dict for substitution, e.g. {km:100*m, Gyr:1e9*yr}\\n \\\"\\\"\\\"\\n # iterative substitution for 10 times or untill all items repaced\\n for i in range(10):\\n count = 0\\n for unit, replace in substitute_dict.items():\\n unit_dict = dict(zip(quantity.unit._bases, quantity.unit._powers))\\n if unit in unit_dict:\\n power = unit_dict[unit]\\n if power != 0:\\n quantity = quantity * (replace / unit) ** power\\n count += 1\\n if count == 0:\\n break\\n return quantity\\n\\n\\ndef convert_ffb_threhold(quantity):\\n return subs_unit(quantity, {Mv_sc: zp1_sc**-6.2000001})\\n # use 6.2000001 instead of 6.2 to avoid displying fraction\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"def subs_unit(quantity, substitute_dict):\n",
" \"substitute unit by replacement\"\n",
" \"\"\"\n",
" quantity : astropy.units.Quantity\n",
" \n",
" substitute_dict : dict\n",
" Dict for substitution, e.g. {km:100*m, Gyr:1e9*yr}\n",
" \"\"\"\n",
" # iterative substitution for 10 times or untill all items repaced\n",
" for i in range(10):\n",
" count = 0\n",
" for unit, replace in substitute_dict.items():\n",
" unit_dict = dict(zip(quantity.unit._bases, quantity.unit._powers))\n",
" if unit in unit_dict:\n",
" power = unit_dict[unit]\n",
" if power != 0:\n",
" quantity = quantity * (replace / unit) ** power\n",
" count += 1\n",
" if count == 0:\n",
" break\n",
" return quantity\n",
"\n",
"\n",
"def convert_ffb_threhold(quantity):\n",
" return subs_unit(quantity, {Mv_sc: zp1_sc**-6.2000001})\n",
" # use 6.2000001 instead of 6.2 to avoid displying fraction"
]
},
{
"cell_type": "markdown",
"id": "210fc519",
"metadata": {},
"source": [
"# Steady wind"
]
},
{
"cell_type": "markdown",
"id": "9dcf9e2b",
"metadata": {},
"source": [
"## wind speed"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "8238ffcb",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$3333.0977 \\; \\mathrm{\\frac{km}{\\eta^{1/2}\\,s}}$"
],
"text/plain": [
"<Quantity 3333.09773416 km / (eta(1/2) s)>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 5;\n",
" var nbb_unformatted_code = \"M_sn = 0.2 * (SFR_sc * Msun / yr)\\nE_sn = 7e41 * erg / s * SFR_sc\\nM_out = eta * M_sn\\nM_in = eps**-1 * (SFR_sc * Msun / yr)\\n\\nVwind = (2 * E_sn / M_out) ** 0.5\\nVwind = Vwind.decompose([km, s, kg] + bases)\\nVwind\";\n",
" var nbb_formatted_code = \"M_sn = 0.2 * (SFR_sc * Msun / yr)\\nE_sn = 7e41 * erg / s * SFR_sc\\nM_out = eta * M_sn\\nM_in = eps**-1 * (SFR_sc * Msun / yr)\\n\\nVwind = (2 * E_sn / M_out) ** 0.5\\nVwind = Vwind.decompose([km, s, kg] + bases)\\nVwind\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"M_sn = 0.2 * (SFR_sc * Msun / yr)\n",
"E_sn = 7e41 * erg / s * SFR_sc\n",
"M_out = eta * M_sn\n",
"M_in = eps**-1 * (SFR_sc * Msun / yr)\n",
"\n",
"Vwind = (2 * E_sn / M_out) ** 0.5\n",
"Vwind = Vwind.decompose([km, s, kg] + bases)\n",
"Vwind"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "8dd61110",
"metadata": {
},
"outputs": [
{
"data": {
"text/plain": [
"1.988409870698051e+30"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 6;\n",
" var nbb_unformatted_code = \"Msun.to(kg)\";\n",
" var nbb_formatted_code = \"Msun.to(kg)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Msun.to(kg)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "dbbf183a",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/latex": [
"$3323.4259 \\; \\mathrm{\\frac{km}{s}}$"
],
"text/plain": [
"<Quantity 3323.42594321 km / s>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 7;\n",
" var nbb_unformatted_code = \"((2 * 7e41 * erg / s / (0.2 * 2e30 * kg / yr)) ** 0.5).to(km / s)\";\n",
" var nbb_formatted_code = \"((2 * 7e41 * erg / s / (0.2 * 2e30 * kg / yr)) ** 0.5).to(km / s)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"((2 * 7e41 * erg / s / (0.2 * 2e30 * kg / yr)) ** 0.5).to(km / s)"
]
},
{
"cell_type": "markdown",
"id": "cdd78906",
"metadata": {},
"source": [
"## density"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "565281b0",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 8;\n",
" var nbb_unformatted_code = \"R_sc = u.def_unit(\\\"(R/kpc)\\\", format={\\\"latex\\\": r\\\"R_{kpc}\\\"})\";\n",
" var nbb_formatted_code = \"R_sc = u.def_unit(\\\"(R/kpc)\\\", format={\\\"latex\\\": r\\\"R_{kpc}\\\"})\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"R_sc = u.def_unit(\"(R/kpc)\", format={\"latex\": r\"R_{kpc}\"})"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "b896c91d",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 9;\n",
" var nbb_unformatted_code = \"rho_inner = 0.3 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2\\nrho_outer = 0.05 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2\";\n",
" var nbb_formatted_code = \"rho_inner = 0.3 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2\\nrho_outer = 0.05 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rho_inner = 0.3 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2\n",
"rho_outer = 0.05 * (M_out / E_sn) ** 0.5 * M_out / (R_sc * kpc) ** 2"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "25c892c0",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1.0072129 \\; \\mathrm{\\frac{\\eta^{3/2}\\,M_{p}\\,SFR_{M_\\odot/yr}}{R_{kpc}^{2}\\,cm^{3}}}$"
],
"text/plain": [
"<Quantity 1.00721287 [SFR/(Msun/yr)] eta(3/2) M_p / ((R/kpc)2 cm3)>"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 10;\n",
" var nbb_unformatted_code = \"rho_inner.decompose([M_p, cm, s, R_sc] + bases) / 1e-3\";\n",
" var nbb_formatted_code = \"rho_inner.decompose([M_p, cm, s, R_sc] + bases) / 1e-3\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rho_inner.decompose([M_p, cm, s, R_sc] + bases) / 1e-3"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "8a805a26",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1.6786881 \\; \\mathrm{\\frac{\\eta^{3/2}\\,M_{p}\\,SFR_{M_\\odot/yr}}{R_{kpc}^{2}\\,cm^{3}}}$"
],
"text/plain": [
"<Quantity 1.67868812 [SFR/(Msun/yr)] eta(3/2) M_p / ((R/kpc)2 cm3)>"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 11;\n",
" var nbb_unformatted_code = \"rho_outer.decompose([M_p, cm, s, R_sc] + bases) / 1e-4\";\n",
" var nbb_formatted_code = \"rho_outer.decompose([M_p, cm, s, R_sc] + bases) / 1e-4\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rho_outer.decompose([M_p, cm, s, R_sc] + bases) / 1e-4"
]
},
{
"cell_type": "markdown",
"id": "604a5c30",
"metadata": {},
"source": [
"# Radius"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "4cb8673d",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 12;\n",
" var nbb_unformatted_code = \"Vvir = 148 * km / s * Mv_sc ** (1 / 3) * zp1_sc**0.5\\nRvir = 12.3 * kpc * Mv_sc ** (1 / 3) * zp1_sc**-1\";\n",
" var nbb_formatted_code = \"Vvir = 148 * km / s * Mv_sc ** (1 / 3) * zp1_sc**0.5\\nRvir = 12.3 * kpc * Mv_sc ** (1 / 3) * zp1_sc**-1\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Vvir = 148 * km / s * Mv_sc ** (1 / 3) * zp1_sc**0.5\n",
"Rvir = 12.3 * kpc * Mv_sc ** (1 / 3) * zp1_sc**-1"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "28a2236e",
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/latex": [
"$81.262462 \\; \\mathrm{\\frac{Myr}{(1+z)_{10}^{3/2}}}$"
],
"text/plain": [
"<Quantity 81.26246167 Myr / [(1+z)/10](3/2)>"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 13;\n",
" var nbb_unformatted_code = \"tvir = Rvir / Vvir\\ntvir.decompose([Myr, kpc] + bases)\";\n",
" var nbb_formatted_code = \"tvir = Rvir / Vvir\\ntvir.decompose([Myr, kpc] + bases)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tvir = Rvir / Vvir\n",
"tvir.decompose([Myr, kpc] + bases)"
]
},
{
"cell_type": "markdown",
"id": "1703f95e",
"metadata": {},
"source": [
"## Shell"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "6aefc3cf",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1.1260465 \\; \\mathrm{\\frac{\\epsilon\\,\\eta^{1/2}}{M_{v,10.8}^{1/3}\\,(1+z)_{10}^{1/2}}}$"
],
"text/plain": [
"<Quantity 1.12604653 eps eta(1/2) / ([(1+z)/10](1/2) [Mv/10^10.8Msun](1/3))>"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 14;\n",
" var nbb_unformatted_code = \"ratio_Rsh2_Rstr2 = Vwind * M_out / (4 * Vvir * M_in)\\nratio_Rsh2_Rstr2\";\n",
" var nbb_formatted_code = \"ratio_Rsh2_Rstr2 = Vwind * M_out / (4 * Vvir * M_in)\\nratio_Rsh2_Rstr2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ratio_Rsh2_Rstr2 = Vwind * M_out / (4 * Vvir * M_in)\n",
"ratio_Rsh2_Rstr2"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "3f6b0103",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.5289 \\; \\mathrm{\\frac{kpc\\,\\lambda_{s,.025}\\,M_{v,10.8}^{1/3}}{(1+z)_{10}^{1/2}}}$"
],
"text/plain": [
"<Quantity 0.5289 [Mv/10^10.8Msun](1/3) kpc lambda_s025 / [(1+z)/10](1/2)>"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 15;\n",
" var nbb_unformatted_code = \"Rstr = 0.043 * lambdas_sc * zp1_sc**0.5 * Rvir\\nRstr\";\n",
" var nbb_formatted_code = \"Rstr = 0.043 * lambdas_sc * zp1_sc**0.5 * Rvir\\nRstr\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rstr = 0.043 * lambdas_sc * zp1_sc**0.5 * Rvir\n",
"Rstr"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "ebba55ff",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/latex": [
"$0.56124403 \\; \\mathrm{\\frac{kpc\\,\\lambda_{s,.025}\\,\\epsilon^{1/2}\\,\\eta^{1/4}\\,M_{v,10.8}^{1/6}}{(1+z)_{10}^{3/4}}}$"
],
"text/plain": [
"<Quantity 0.56124403 [Mv/10^10.8Msun](1/6) eps(1/2) eta(1/4) kpc lambda_s025 / [(1+z)/10](3/4)>"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 16;\n",
" var nbb_unformatted_code = \"Rsh = Rstr * (ratio_Rsh2_Rstr2) ** 0.5\\nRsh\";\n",
" var nbb_formatted_code = \"Rsh = Rstr * (ratio_Rsh2_Rstr2) ** 0.5\\nRsh\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rsh = Rstr * (ratio_Rsh2_Rstr2) ** 0.5\n",
"Rsh"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "19aba0a6",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.56124403 \\; \\mathrm{\\frac{kpc\\,\\lambda_{s,.025}\\,\\epsilon^{1/2}\\,\\eta^{1/4}}{(1+z)_{10}^{1.78333335}}}$"
],
"text/plain": [
"<Quantity 0.56124403 eps(1/2) eta(1/4) kpc lambda_s025 / [(1+z)/10](1.78333335)>"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 17;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(Rsh)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(Rsh)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(Rsh)"
]
},
{
"cell_type": "markdown",
"id": "dcad5651",
"metadata": {},
"source": [
"## Disc"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "c033ba3c",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/latex": [
"$0.3075 \\; \\mathrm{\\frac{kpc\\,\\lambda_{d,.025}\\,M_{v,10.8}^{1/3}}{(1+z)_{10}}}$"
],
"text/plain": [
"<Quantity 0.3075 [Mv/10^10.8Msun](1/3) kpc lambda_d025 / [(1+z)/10]>"
]
},
"execution_count": 18,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 18;\n",
" var nbb_unformatted_code = \"Rd = (lambdad_sc * 0.025) * Rvir\\nRd\";\n",
" var nbb_formatted_code = \"Rd = (lambdad_sc * 0.025) * Rvir\\nRd\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rd = (lambdad_sc * 0.025) * Rvir\n",
"Rd"
]
},
{
"cell_type": "code",
"execution_count": 19,
"id": "efebe0f4",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.3075 \\; \\mathrm{\\frac{kpc\\,\\lambda_{d,.025}}{(1+z)_{10}^{3.0666667}}}$"
],
"text/plain": [
"<Quantity 0.3075 kpc lambda_d025 / [(1+z)/10](3.0666667)>"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 19;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(Rd)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(Rd)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(Rd)"
]
},
{
"cell_type": "markdown",
"id": "e1e9f03b",
"metadata": {},
"source": [
"## Galaxy radius, 2Re"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "fd8212e1",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.56124403 \\; \\mathrm{\\frac{kpc\\,\\lambda_{s,.025}\\,\\epsilon^{1/2}\\,\\eta^{1/4}\\,M_{v,10.8}^{1/6}}{(1+z)_{10}^{3/4}}}$"
],
"text/plain": [
"<Quantity 0.56124403 [Mv/10^10.8Msun](1/6) eps(1/2) eta(1/4) kpc lambda_s025 / [(1+z)/10](3/4)>"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 20;\n",
" var nbb_unformatted_code = \"Rsh\";\n",
" var nbb_formatted_code = \"Rsh\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rsh"
]
},
{
"cell_type": "code",
"execution_count": 21,
"id": "f34c3343",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.615 \\; \\mathrm{\\frac{kpc\\,\\lambda_{d,.025}\\,M_{v,10.8}^{1/3}}{(1+z)_{10}}}$"
],
"text/plain": [
"<Quantity 0.615 [Mv/10^10.8Msun](1/3) kpc lambda_d025 / [(1+z)/10]>"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 21;\n",
" var nbb_unformatted_code = \"2 * Rd\";\n",
" var nbb_formatted_code = \"2 * Rd\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"2 * Rd"
]
},
{
"cell_type": "markdown",
"id": "3cbc19e9",
"metadata": {},
"source": [
"# Gas fraction"
]
},
{
"cell_type": "code",
"execution_count": 22,
"id": "a192163c",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 22;\n",
" var nbb_unformatted_code = \"Mgas_sh = 1.04e5 * Msun * eta**1.5 * SFR_sc * (Rsh / kpc)\\nMgas_d = 1.04e5 * Msun * eta**1.5 * SFR_sc * (2 * Rd / kpc)\";\n",
" var nbb_formatted_code = \"Mgas_sh = 1.04e5 * Msun * eta**1.5 * SFR_sc * (Rsh / kpc)\\nMgas_d = 1.04e5 * Msun * eta**1.5 * SFR_sc * (2 * Rd / kpc)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Mgas_sh = 1.04e5 * Msun * eta**1.5 * SFR_sc * (Rsh / kpc)\n",
"Mgas_d = 1.04e5 * Msun * eta**1.5 * SFR_sc * (2 * Rd / kpc)"
]
},
{
"cell_type": "code",
"execution_count": 23,
"id": "f18b76a9",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 23;\n",
" var nbb_unformatted_code = \"SFR = 65 * Msun / yr * eps * Mv_sc**1.14000001 * zp1_sc**2.5\\nMstar = 0.16 * eps * Mv_sc * 10**10.8 * Msun\";\n",
" var nbb_formatted_code = \"SFR = 65 * Msun / yr * eps * Mv_sc**1.14000001 * zp1_sc**2.5\\nMstar = 0.16 * eps * Mv_sc * 10**10.8 * Msun\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"SFR = 65 * Msun / yr * eps * Mv_sc**1.14000001 * zp1_sc**2.5\n",
"Mstar = 0.16 * eps * Mv_sc * 10**10.8 * Msun"
]
},
{
"cell_type": "code",
"execution_count": 24,
"id": "8b8ef5cf",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 24;\n",
" var nbb_unformatted_code = \"eps_05 = u.def_unit(\\\"(eps/0.5)\\\", format={\\\"latex\\\": r\\\"\\\\epsilon_{0.5}\\\"})\\neta_6 = u.def_unit(\\\"(eta/6)\\\", format={\\\"latex\\\": r\\\"\\\\eta_6\\\"})\\nsubs = {SFR_sc: SFR / (Msun / yr), eps: eps_05 * 0.5, eta: eta_6 * 6}\\n\\nMgas_Mstar_sh = subs_unit(Mgas_sh / Mstar, subs)\\nMgas_Mstar_d = subs_unit(Mgas_d / Mstar, subs)\";\n",
" var nbb_formatted_code = \"eps_05 = u.def_unit(\\\"(eps/0.5)\\\", format={\\\"latex\\\": r\\\"\\\\epsilon_{0.5}\\\"})\\neta_6 = u.def_unit(\\\"(eta/6)\\\", format={\\\"latex\\\": r\\\"\\\\eta_6\\\"})\\nsubs = {SFR_sc: SFR / (Msun / yr), eps: eps_05 * 0.5, eta: eta_6 * 6}\\n\\nMgas_Mstar_sh = subs_unit(Mgas_sh / Mstar, subs)\\nMgas_Mstar_d = subs_unit(Mgas_d / Mstar, subs)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"eps_05 = u.def_unit(\"(eps/0.5)\", format={\"latex\": r\"\\epsilon_{0.5}\"})\n",
"eta_6 = u.def_unit(\"(eta/6)\", format={\"latex\": r\"\\eta_6\"})\n",
"subs = {SFR_sc: SFR / (Msun / yr), eps: eps_05 * 0.5, eta: eta_6 * 6}\n",
"\n",
"Mgas_Mstar_sh = subs_unit(Mgas_sh / Mstar, subs)\n",
"Mgas_Mstar_d = subs_unit(Mgas_d / Mstar, subs)"
]
},
{
"cell_type": "code",
"execution_count": 25,
"id": "a92c1feb",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.61126304 \\; \\mathrm{\\eta_6^{7/4}\\,(1+z)_{10}^{7/4}\\,\\lambda_{s,.025}\\,\\epsilon_{0.5}^{1/2}\\,M_{v,10.8}^{0.3066666766666667}}$"
],
"text/plain": [
"<Quantity 0.61126304 (eps/0.5)(1/2) (eta/6)(7/4) [(1+z)/10](7/4) [Mv/10^10.8Msun](0.3066666766666667) lambda_s025>"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 25;\n",
" var nbb_unformatted_code = \"Mgas_Mstar_sh / 1e-2\";\n",
" var nbb_formatted_code = \"Mgas_Mstar_sh / 1e-2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Mgas_Mstar_sh / 1e-2"
]
},
{
"cell_type": "code",
"execution_count": 26,
"id": "cbb16ee7",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/latex": [
"$0.60524151 \\; \\mathrm{\\eta_6^{3/2}\\,(1+z)_{10}^{3/2}\\,\\lambda_{d,.025}\\,M_{v,10.8}^{0.47333334333333343}}$"
],
"text/plain": [
"<Quantity 0.60524151 (eta/6)(3/2) [(1+z)/10](3/2) [Mv/10^10.8Msun](0.47333334333333343) lambda_d025>"
]
},
"execution_count": 26,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 26;\n",
" var nbb_unformatted_code = \"Mgas_Mstar_d / 1e-2\";\n",
" var nbb_formatted_code = \"Mgas_Mstar_d / 1e-2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Mgas_Mstar_d / 1e-2"
]
},
{
"cell_type": "markdown",
"id": "d22cf1fb",
"metadata": {},
"source": [
"## FFB"
]
},
{
"cell_type": "code",
"execution_count": 27,
"id": "eb70e42c",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.61126304 \\; \\mathrm{\\frac{\\eta_6^{7/4}\\,\\lambda_{s,.025}\\,\\epsilon_{0.5}^{1/2}}{(1+z)_{10}^{0.151333426000001}}}$"
],
"text/plain": [
"<Quantity 0.61126304 (eps/0.5)(1/2) (eta/6)(7/4) lambda_s025 / [(1+z)/10](0.151333426000001)>"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 27;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(Mgas_Mstar_sh) / 1e-2\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(Mgas_Mstar_sh) / 1e-2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(Mgas_Mstar_sh) / 1e-2"
]
},
{
"cell_type": "code",
"execution_count": 28,
"id": "f758eeda",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.60524151 \\; \\mathrm{\\frac{\\eta_6^{3/2}\\,\\lambda_{d,.025}}{(1+z)_{10}^{1.4346667760000016}}}$"
],
"text/plain": [
"<Quantity 0.60524151 (eta/6)(3/2) lambda_d025 / [(1+z)/10](1.4346667760000016)>"
]
},
"execution_count": 28,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 28;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(Mgas_Mstar_d) / 1e-2\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(Mgas_Mstar_d) / 1e-2\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(Mgas_Mstar_d) / 1e-2"
]
},
{
"cell_type": "markdown",
"id": "b561cd20",
"metadata": {},
"source": [
"## changing epsilon"
]
},
{
"cell_type": "code",
"execution_count": 29,
"id": "20fd98bc",
"metadata": {
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"5.664336862474418 6.547900426854397\n",
"0.06148232887224009 0.06804138174397717\n"
]
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 29;\n",
" var nbb_unformatted_code = \"subs_02 = {SFR_sc: SFR / (Msun / yr), eps: 0.2, eta: 21}\\nsubs_05 = {SFR_sc: SFR / (Msun / yr), eps: 0.5, eta: 6}\\nsubs_10 = {SFR_sc: SFR / (Msun / yr), eps: 1, eta: 1}\\n\\nMgas_Mstar_sh_02 = subs_unit(Mgas_sh / Mstar, subs_02)\\nMgas_Mstar_d_02 = subs_unit(Mgas_d / Mstar, subs_02)\\n\\nMgas_Mstar_sh_05 = subs_unit(Mgas_sh / Mstar, subs_05)\\nMgas_Mstar_d_05 = subs_unit(Mgas_d / Mstar, subs_05)\\n\\nMgas_Mstar_sh_10 = subs_unit(Mgas_sh / Mstar, subs_10)\\nMgas_Mstar_d_10 = subs_unit(Mgas_d / Mstar, subs_10)\\n\\nprint(Mgas_Mstar_sh_02 / Mgas_Mstar_sh_05, Mgas_Mstar_d_02 / Mgas_Mstar_d_05)\\n\\nprint(Mgas_Mstar_sh_10 / Mgas_Mstar_sh_05, Mgas_Mstar_d_10 / Mgas_Mstar_d_05)\";\n",
" var nbb_formatted_code = \"subs_02 = {SFR_sc: SFR / (Msun / yr), eps: 0.2, eta: 21}\\nsubs_05 = {SFR_sc: SFR / (Msun / yr), eps: 0.5, eta: 6}\\nsubs_10 = {SFR_sc: SFR / (Msun / yr), eps: 1, eta: 1}\\n\\nMgas_Mstar_sh_02 = subs_unit(Mgas_sh / Mstar, subs_02)\\nMgas_Mstar_d_02 = subs_unit(Mgas_d / Mstar, subs_02)\\n\\nMgas_Mstar_sh_05 = subs_unit(Mgas_sh / Mstar, subs_05)\\nMgas_Mstar_d_05 = subs_unit(Mgas_d / Mstar, subs_05)\\n\\nMgas_Mstar_sh_10 = subs_unit(Mgas_sh / Mstar, subs_10)\\nMgas_Mstar_d_10 = subs_unit(Mgas_d / Mstar, subs_10)\\n\\nprint(Mgas_Mstar_sh_02 / Mgas_Mstar_sh_05, Mgas_Mstar_d_02 / Mgas_Mstar_d_05)\\n\\nprint(Mgas_Mstar_sh_10 / Mgas_Mstar_sh_05, Mgas_Mstar_d_10 / Mgas_Mstar_d_05)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"subs_02 = {SFR_sc: SFR / (Msun / yr), eps: 0.2, eta: 21}\n",
"subs_05 = {SFR_sc: SFR / (Msun / yr), eps: 0.5, eta: 6}\n",
"subs_10 = {SFR_sc: SFR / (Msun / yr), eps: 1, eta: 1}\n",
"\n",
"Mgas_Mstar_sh_02 = subs_unit(Mgas_sh / Mstar, subs_02)\n",
"Mgas_Mstar_d_02 = subs_unit(Mgas_d / Mstar, subs_02)\n",
"\n",
"Mgas_Mstar_sh_05 = subs_unit(Mgas_sh / Mstar, subs_05)\n",
"Mgas_Mstar_d_05 = subs_unit(Mgas_d / Mstar, subs_05)\n",
"\n",
"Mgas_Mstar_sh_10 = subs_unit(Mgas_sh / Mstar, subs_10)\n",
"Mgas_Mstar_d_10 = subs_unit(Mgas_d / Mstar, subs_10)\n",
"\n",
"print(Mgas_Mstar_sh_02 / Mgas_Mstar_sh_05, Mgas_Mstar_d_02 / Mgas_Mstar_d_05)\n",
"\n",
"print(Mgas_Mstar_sh_10 / Mgas_Mstar_sh_05, Mgas_Mstar_d_10 / Mgas_Mstar_d_05)"
]
},
{
"cell_type": "markdown",
"id": "2fe59ea7",
"metadata": {},
"source": [
"# Metalicity"
]
},
{
"cell_type": "code",
"execution_count": 30,
"id": "9f0ce67a",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.22201569 \\; \\mathrm{\\frac{(1+z)_{10}^{1/2}\\,M_{v,10.8}^{1/3}}{\\eta^{1/2}\\,\\epsilon}}$"
],
"text/plain": [
"<Quantity 0.22201569 [(1+z)/10](1/2) [Mv/10^10.8Msun](1/3) / (eps eta(1/2))>"
]
},
"execution_count": 30,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 30;\n",
" var nbb_unformatted_code = \"f_omega = Rstr**2 / (4 * Rsh**2)\\nf_omega\";\n",
" var nbb_formatted_code = \"f_omega = Rstr**2 / (4 * Rsh**2)\\nf_omega\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"f_omega = Rstr**2 / (4 * Rsh**2)\n",
"f_omega"
]
},
{
"cell_type": "markdown",
"id": "95a74706",
"metadata": {
},
"source": [
"# Dust"
]
},
{
"cell_type": "code",
"execution_count": 31,
"id": "9b061a76",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 31;\n",
" var nbb_unformatted_code = \"NH_inner = 3.08e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc\\nNH_outer = 0.52e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc\";\n",
" var nbb_formatted_code = \"NH_inner = 3.08e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc\\nNH_outer = 0.52e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"NH_inner = 3.08e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc\n",
"NH_outer = 0.52e18 * cm**-2 * eta ** (3 / 2) * R_sc**-1 * SFR_sc"
]
},
{
"cell_type": "code",
"execution_count": 32,
"id": "e1fde876",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 32;\n",
" var nbb_unformatted_code = \"fdsn = u.def_unit(\\\"f_dsn\\\", format={\\\"latex\\\": r\\\"f_{d,sn}\\\"})\\nfdsn_65 = u.def_unit(\\\"(f_dsn/6.5)\\\", format={\\\"latex\\\": r\\\"f_{d,sn,6.5}\\\"})\";\n",
" var nbb_formatted_code = \"fdsn = u.def_unit(\\\"f_dsn\\\", format={\\\"latex\\\": r\\\"f_{d,sn}\\\"})\\nfdsn_65 = u.def_unit(\\\"(f_dsn/6.5)\\\", format={\\\"latex\\\": r\\\"f_{d,sn,6.5}\\\"})\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fdsn = u.def_unit(\"f_dsn\", format={\"latex\": r\"f_{d,sn}\"})\n",
"fdsn_65 = u.def_unit(\"(f_dsn/6.5)\", format={\"latex\": r\"f_{d,sn,6.5}\"})"
]
},
{
"cell_type": "code",
"execution_count": 33,
"id": "73a6f0fd",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 33;\n",
" var nbb_unformatted_code = \"sig_ext_mw = 1e-21 * cm**2\\nfd = eta**-1 * fdsn\\n\\ntau_inner = NH_inner * sig_ext_mw * fd\\ntau_outer = NH_outer * sig_ext_mw * fd\";\n",
" var nbb_formatted_code = \"sig_ext_mw = 1e-21 * cm**2\\nfd = eta**-1 * fdsn\\n\\ntau_inner = NH_inner * sig_ext_mw * fd\\ntau_outer = NH_outer * sig_ext_mw * fd\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sig_ext_mw = 1e-21 * cm**2\n",
"fd = eta**-1 * fdsn\n",
"\n",
"tau_inner = NH_inner * sig_ext_mw * fd\n",
"tau_outer = NH_outer * sig_ext_mw * fd"
]
},
{
"cell_type": "code",
"execution_count": 34,
"id": "d0e8fef5",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$3.08 \\; \\mathrm{\\frac{SFR_{M_\\odot/yr}\\,f_{d,sn}\\,\\eta^{1/2}}{R_{kpc}}}$"
],
"text/plain": [
"<Quantity 3.08 [SFR/(Msun/yr)] eta(1/2) f_dsn / (R/kpc)>"
]
},
"execution_count": 34,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 34;\n",
" var nbb_unformatted_code = \"tau_inner / 1e-3\";\n",
" var nbb_formatted_code = \"tau_inner / 1e-3\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_inner / 1e-3"
]
},
{
"cell_type": "code",
"execution_count": 35,
"id": "2dc197e2",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.52 \\; \\mathrm{\\frac{SFR_{M_\\odot/yr}\\,f_{d,sn}\\,\\eta^{1/2}}{R_{kpc}}}$"
],
"text/plain": [
"<Quantity 0.52 [SFR/(Msun/yr)] eta(1/2) f_dsn / (R/kpc)>"
]
},
"execution_count": 35,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 35;\n",
" var nbb_unformatted_code = \"tau_outer / 1e-3\";\n",
" var nbb_formatted_code = \"tau_outer / 1e-3\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_outer / 1e-3"
]
},
{
"cell_type": "markdown",
"id": "fe8acdff",
"metadata": {},
"source": [
"## FFB general"
]
},
{
"cell_type": "code",
"execution_count": 36,
"id": "10a90fb6",
"metadata": {
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/lizz/local/conda/envs/main/lib/python3.8/site-packages/astropy/units/quantity.py:614: RuntimeWarning: divide by zero encountered in reciprocal\n",
" result = super().__array_ufunc__(function, method, *arrays, **kwargs)\n"
]
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 36;\n",
" var nbb_unformatted_code = \"subs_sh = {SFR_sc: SFR / (Msun / yr), R_sc: Rsh / kpc, fdsn: 6.5 * fdsn_65}\\ntau_inner_shell = subs_unit(tau_inner, subs_sh)\\ntau_outer_shell = subs_unit(tau_outer, subs_sh)\\n\\nsubs_d = {SFR_sc: SFR / (Msun / yr), R_sc: 2 * Rd / kpc, fdsn: 6.5 * fdsn_65}\\ntau_inner_disc = subs_unit(tau_inner, subs_d)\\ntau_outer_disc = subs_unit(tau_outer, subs_d)\";\n",
" var nbb_formatted_code = \"subs_sh = {SFR_sc: SFR / (Msun / yr), R_sc: Rsh / kpc, fdsn: 6.5 * fdsn_65}\\ntau_inner_shell = subs_unit(tau_inner, subs_sh)\\ntau_outer_shell = subs_unit(tau_outer, subs_sh)\\n\\nsubs_d = {SFR_sc: SFR / (Msun / yr), R_sc: 2 * Rd / kpc, fdsn: 6.5 * fdsn_65}\\ntau_inner_disc = subs_unit(tau_inner, subs_d)\\ntau_outer_disc = subs_unit(tau_outer, subs_d)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"subs_sh = {SFR_sc: SFR / (Msun / yr), R_sc: Rsh / kpc, fdsn: 6.5 * fdsn_65}\n",
"tau_inner_shell = subs_unit(tau_inner, subs_sh)\n",
"tau_outer_shell = subs_unit(tau_outer, subs_sh)\n",
"\n",
"subs_d = {SFR_sc: SFR / (Msun / yr), R_sc: 2 * Rd / kpc, fdsn: 6.5 * fdsn_65}\n",
"tau_inner_disc = subs_unit(tau_inner, subs_d)\n",
"tau_outer_disc = subs_unit(tau_outer, subs_d)"
]
},
{
"cell_type": "code",
"execution_count": 37,
"id": "3f9010a6",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$2.3185993 \\; \\mathrm{\\frac{(1+z)_{10}^{13/4}\\,f_{d,sn,6.5}\\,M_{v,10.8}^{0.9733333433333334}\\,\\epsilon^{1/2}\\,\\eta^{1/4}}{\\lambda_{s,.025}}}$"
],
"text/plain": [
"<Quantity 2.31859927 (f_dsn/6.5) [(1+z)/10](13/4) [Mv/10^10.8Msun](0.9733333433333334) eps(1/2) eta(1/4) / lambda_s025>"
]
},
"execution_count": 37,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 37;\n",
" var nbb_unformatted_code = \"tau_inner_shell\";\n",
" var nbb_formatted_code = \"tau_inner_shell\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_inner_shell"
]
},
{
"cell_type": "code",
"execution_count": 38,
"id": "39f5d0dc",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.39145182 \\; \\mathrm{\\frac{(1+z)_{10}^{13/4}\\,f_{d,sn,6.5}\\,M_{v,10.8}^{0.9733333433333334}\\,\\epsilon^{1/2}\\,\\eta^{1/4}}{\\lambda_{s,.025}}}$"
],
"text/plain": [
"<Quantity 0.39145182 (f_dsn/6.5) [(1+z)/10](13/4) [Mv/10^10.8Msun](0.9733333433333334) eps(1/2) eta(1/4) / lambda_s025>"
]
},
"execution_count": 38,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 38;\n",
" var nbb_unformatted_code = \"tau_outer_shell\";\n",
" var nbb_formatted_code = \"tau_outer_shell\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_outer_shell"
]
},
{
"cell_type": "code",
"execution_count": 39,
"id": "8d966663",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$2.115935 \\; \\mathrm{\\frac{(1+z)_{10}^{7/2}\\,f_{d,sn,6.5}\\,\\epsilon\\,M_{v,10.8}^{0.8066666766666668}\\,\\eta^{1/2}}{\\lambda_{d,.025}}}$"
],
"text/plain": [
"<Quantity 2.11593496 (f_dsn/6.5) [(1+z)/10](7/2) [Mv/10^10.8Msun](0.8066666766666668) eps eta(1/2) / lambda_d025>"
]
},
"execution_count": 39,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 39;\n",
" var nbb_unformatted_code = \"tau_inner_disc\";\n",
" var nbb_formatted_code = \"tau_inner_disc\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_inner_disc"
]
},
{
"cell_type": "code",
"execution_count": 40,
"id": "11061904",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$0.35723577 \\; \\mathrm{\\frac{(1+z)_{10}^{7/2}\\,f_{d,sn,6.5}\\,\\epsilon\\,M_{v,10.8}^{0.8066666766666668}\\,\\eta^{1/2}}{\\lambda_{d,.025}}}$"
],
"text/plain": [
"<Quantity 0.35723577 (f_dsn/6.5) [(1+z)/10](7/2) [Mv/10^10.8Msun](0.8066666766666668) eps eta(1/2) / lambda_d025>"
]
},
"execution_count": 40,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 40;\n",
" var nbb_unformatted_code = \"tau_outer_disc\";\n",
" var nbb_formatted_code = \"tau_outer_disc\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tau_outer_disc"
]
},
{
"cell_type": "markdown",
"id": "30fdfa47",
"metadata": {},
"source": [
"## FFB threshold"
]
},
{
"cell_type": "code",
"execution_count": 41,
"id": "4d33911b",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$2.3185993 \\; \\mathrm{\\frac{f_{d,sn,6.5}\\,\\epsilon^{1/2}\\,\\eta^{1/4}}{\\lambda_{s,.025}\\,(1+z)_{10}^{2.7846668260000014}}}$"
],
"text/plain": [
"<Quantity 2.31859927 (f_dsn/6.5) eps(1/2) eta(1/4) / ([(1+z)/10](2.7846668260000014) lambda_s025)>"
]
},
"execution_count": 41,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 41;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(tau_inner_shell)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(tau_inner_shell)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(tau_inner_shell)"
]
},
{
"cell_type": "code",
"execution_count": 42,
"id": "8e7b82a6",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$2.115935 \\; \\mathrm{\\frac{f_{d,sn,6.5}\\,\\epsilon\\,\\eta^{1/2}}{\\lambda_{d,.025}\\,(1+z)_{10}^{1.5013334760000019}}}$"
],
"text/plain": [
"<Quantity 2.11593496 (f_dsn/6.5) eps eta(1/2) / ([(1+z)/10](1.5013334760000019) lambda_d025)>"
]
},
"execution_count": 42,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 42;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(tau_inner_disc)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(tau_inner_disc)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(tau_inner_disc)"
]
},
{
"cell_type": "markdown",
"id": "ed076e3f",
"metadata": {},
"source": [
"# Cooling radius"
]
},
{
"cell_type": "code",
"execution_count": 43,
"id": "165f251a",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 43;\n",
" var nbb_unformatted_code = \"Rcool = (\\n 4\\n * kpc\\n * (0.2 * eta) ** -2.9200001\\n * (Rsh / (0.3 * kpc)) ** 1.7900001\\n * (0.3 * SFR / (10 * Msun / yr)) ** -0.7900001\\n)\";\n",
" var nbb_formatted_code = \"Rcool = (\\n 4\\n * kpc\\n * (0.2 * eta) ** -2.9200001\\n * (Rsh / (0.3 * kpc)) ** 1.7900001\\n * (0.3 * SFR / (10 * Msun / yr)) ** -0.7900001\\n)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rcool = (\n",
" 4\n",
" * kpc\n",
" * (0.2 * eta) ** -2.9200001\n",
" * (Rsh / (0.3 * kpc)) ** 1.7900001\n",
" * (0.3 * SFR / (10 * Msun / yr)) ** -0.7900001\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 44,
"id": "1432ff66",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1418.1036 \\; \\mathrm{\\frac{\\lambda_{s,.025}^{0.7900001000000001}}{\\epsilon^{0.39500004999999994}\\,M_{v,10.8}^{0.7689334385666676}\\,(1+z)_{10}^{2.567500325}\\,\\eta^{2.722500075}}}$"
],
"text/plain": [
"<Quantity 1418.10358308 lambda_s025(0.7900001000000001) / ([(1+z)/10](2.567500325) [Mv/10^10.8Msun](0.7689334385666676) eps(0.39500004999999994) eta(2.722500075))>"
]
},
"execution_count": 44,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 44;\n",
" var nbb_unformatted_code = \"Rcool / Rsh\";\n",
" var nbb_formatted_code = \"Rcool / Rsh\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Rcool / Rsh"
]
},
{
"cell_type": "code",
"execution_count": 45,
"id": "5c0ea8de",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1418.1036 \\; \\mathrm{\\frac{(1+z)_{10}^{2.1998870710066827}\\,\\lambda_{s,.025}^{0.7900001000000001}}{\\epsilon^{0.39500004999999994}\\,\\eta^{2.722500075}}}$"
],
"text/plain": [
"<Quantity 1418.10358308 [(1+z)/10](2.1998870710066827) lambda_s025(0.7900001000000001) / (eps(0.39500004999999994) eta(2.722500075))>"
]
},
"execution_count": 45,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 45;\n",
" var nbb_unformatted_code = \"ratio_Rcool_Rsh_ffb = convert_ffb_threhold(Rcool / Rsh)\\nratio_Rcool_Rsh_ffb\";\n",
" var nbb_formatted_code = \"ratio_Rcool_Rsh_ffb = convert_ffb_threhold(Rcool / Rsh)\\nratio_Rcool_Rsh_ffb\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ratio_Rcool_Rsh_ffb = convert_ffb_threhold(Rcool / Rsh)\n",
"ratio_Rcool_Rsh_ffb"
]
},
{
"cell_type": "code",
"execution_count": 46,
"id": "c4242c28",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1418.1036 \\; \\mathrm{\\frac{(1+z)_{10}^{2.1998870710066827}\\,\\lambda_{s,.025}^{0.7900001000000001}}{\\epsilon^{0.39500004999999994}\\,\\eta^{2.722500075}}}$"
],
"text/plain": [
"<Quantity 1418.10358308 [(1+z)/10](2.1998870710066827) lambda_s025(0.7900001000000001) / (eps(0.39500004999999994) eta(2.722500075))>"
]
},
"execution_count": 46,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 46;\n",
" var nbb_unformatted_code = \"ratio_Rcool_Rsh_ffb\";\n",
" var nbb_formatted_code = \"ratio_Rcool_Rsh_ffb\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ratio_Rcool_Rsh_ffb"
]
},
{
"cell_type": "code",
"execution_count": 47,
"id": "7a417641",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1.0083403 \\; \\mathrm{\\lambda_{s,.025}^{0.7900001000000001}}$"
],
"text/plain": [
"<Quantity 1.00834033 lambda_s025(0.7900001000000001)>"
]
},
"execution_count": 47,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 47;\n",
" var nbb_unformatted_code = \"eps_val = 0.23\\neta_val = 5 / eps_val - 4\\nsubs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: 1})\";\n",
" var nbb_formatted_code = \"eps_val = 0.23\\neta_val = 5 / eps_val - 4\\nsubs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: 1})\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"eps_val = 0.23\n",
"eta_val = 5 / eps_val - 4\n",
"subs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: 1})"
]
},
{
"cell_type": "code",
"execution_count": 48,
"id": "0d1b1e01",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$1.0051915 \\; \\mathrm{\\lambda_{s,.025}^{0.7900001000000001}}$"
],
"text/plain": [
"<Quantity 1.00519151 lambda_s025(0.7900001000000001)>"
]
},
"execution_count": 48,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 48;\n",
" var nbb_unformatted_code = \"eps_val = 0.2\\neta_val = 5 / eps_val - 4\\nsubs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: (11 + 1) / 10})\";\n",
" var nbb_formatted_code = \"eps_val = 0.2\\neta_val = 5 / eps_val - 4\\nsubs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: (11 + 1) / 10})\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"eps_val = 0.2\n",
"eta_val = 5 / eps_val - 4\n",
"subs_unit(ratio_Rcool_Rsh_ffb, {eps: eps_val, eta: eta_val, zp1_sc: (11 + 1) / 10})"
]
},
{
"cell_type": "markdown",
"id": "14641811",
"metadata": {},
"source": [
"# Star formation history"
]
},
{
"cell_type": "code",
"execution_count": 49,
"id": "2253346b",
"metadata": {
},
"outputs": [
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 49;\n",
" var nbb_unformatted_code = \"N_gen_shell = 7.1 * Mv_sc ** (1 / 3) * zp1_sc\\nN_gen_disc = 10.6 * lambdad_sc ** (-5 / 6) * Mv_sc**-0.05 * zp1_sc ** (-1 / 3)\";\n",
" var nbb_formatted_code = \"N_gen_shell = 7.1 * Mv_sc ** (1 / 3) * zp1_sc\\nN_gen_disc = 10.6 * lambdad_sc ** (-5 / 6) * Mv_sc**-0.05 * zp1_sc ** (-1 / 3)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"N_gen_shell = 7.1 * Mv_sc ** (1 / 3) * zp1_sc\n",
"N_gen_disc = 10.6 * lambdad_sc ** (-5 / 6) * Mv_sc**-0.05 * zp1_sc ** (-1 / 3)"
]
},
{
"cell_type": "code",
"execution_count": 50,
"id": "61cd7982",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$7.1 \\; \\mathrm{\\frac{1}{(1+z)_{10}^{1.0666666999999999}}}$"
],
"text/plain": [
"<Quantity 7.1 1 / [(1+z)/10](1.0666666999999999)>"
]
},
"execution_count": 50,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 50;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(N_gen_shell)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(N_gen_shell)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(N_gen_shell)"
]
},
{
"cell_type": "code",
"execution_count": 51,
"id": "e74acfa6",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$10.6 \\; \\mathrm{\\frac{1}{(1+z)_{10}^{0.023333328333333292}\\,\\lambda_{d,.025}^{5/6}}}$"
],
"text/plain": [
"<Quantity 10.6 1 / ([(1+z)/10](0.023333328333333292) lambda_d025(5/6))>"
]
},
"execution_count": 51,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 51;\n",
" var nbb_unformatted_code = \"convert_ffb_threhold(N_gen_disc)\";\n",
" var nbb_formatted_code = \"convert_ffb_threhold(N_gen_disc)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"convert_ffb_threhold(N_gen_disc)"
]
},
{
"cell_type": "code",
"execution_count": 52,
"id": "65cf295c",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$11.445417 \\; \\mathrm{\\frac{Myr}{(1+z)_{10}^{0.43333330000000014}}}$"
],
"text/plain": [
"<Quantity 11.44541714 Myr / [(1+z)/10](0.43333330000000014)>"
]
},
"execution_count": 52,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 52;\n",
" var nbb_unformatted_code = \"(tvir / convert_ffb_threhold(N_gen_shell)).decompose([Myr, kpc, Msun] + bases)\";\n",
" var nbb_formatted_code = \"(tvir / convert_ffb_threhold(N_gen_shell)).decompose([Myr, kpc, Msun] + bases)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"(tvir / convert_ffb_threhold(N_gen_shell)).decompose([Myr, kpc, Msun] + bases)"
]
},
{
"cell_type": "code",
"execution_count": 53,
"id": "c0f425d3",
"metadata": {
},
"outputs": [
{
"data": {
"text/latex": [
"$7.66627 \\; \\mathrm{\\frac{Myr\\,\\lambda_{d,.025}^{5/6}}{(1+z)_{10}^{1.4766666716666668}}}$"
],
"text/plain": [
"<Quantity 7.66626997 lambda_d025(5/6) Myr / [(1+z)/10](1.4766666716666668)>"
]
},
"execution_count": 53,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"application/javascript": [
"\n",
" setTimeout(function() {\n",
" var nbb_cell_id = 53;\n",
" var nbb_unformatted_code = \"(tvir / convert_ffb_threhold(N_gen_disc)).decompose([Myr, kpc, Msun] + bases)\";\n",
" var nbb_formatted_code = \"(tvir / convert_ffb_threhold(N_gen_disc)).decompose([Myr, kpc, Msun] + bases)\";\n",
" var nbb_cells = Jupyter.notebook.get_cells();\n",
" for (var i = 0; i < nbb_cells.length; ++i) {\n",
" if (nbb_cells[i].input_prompt_number == nbb_cell_id) {\n",
" if (nbb_cells[i].get_text() == nbb_unformatted_code) {\n",
" nbb_cells[i].set_text(nbb_formatted_code);\n",
" }\n",
" break;\n",
" }\n",
" }\n",
" }, 500);\n",
" "
],
"text/plain": [
"<IPython.core.display.Javascript object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"(tvir / convert_ffb_threhold(N_gen_disc)).decompose([Myr, kpc, Msun] + bases)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.12"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": false,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment