Blog: Update 2/22: column density.nb

File column density.nb, 5.6 KB (added by adebrech, 9 years ago)
Line 
1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 11.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 158, 7]
13NotebookDataLength[ 5575, 165]
14NotebookOptionsPosition[ 5046, 142]
15NotebookOutlinePosition[ 5380, 157]
16CellTagsIndexPosition[ 5337, 154]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[{
22 RowBox[{
23 RowBox[{"mH", "=",
24 RowBox[{"1.6733", "*",
25 SuperscriptBox["10",
26 RowBox[{"-", "24"}]]}]}], ";"}], "\[IndentingNewLine]",
27 RowBox[{
28 RowBox[{"lScale", "=",
29 RowBox[{"3.427", "*",
30 SuperscriptBox["10", "11"]}]}], ";"}], "\[IndentingNewLine]",
31 RowBox[{
32 RowBox[{"velScale", "=", "3636589"}], ";"}], "\[IndentingNewLine]",
33 RowBox[{
34 RowBox[{"\[Rho]Scale", "=",
35 RowBox[{"3.2112", "*",
36 SuperscriptBox["10",
37 RowBox[{"-", "15"}]]}]}], ";"}], "\[IndentingNewLine]",
38 RowBox[{
39 RowBox[{"timeScale", "=", "94234.69"}], ";"}], "\[IndentingNewLine]",
40 RowBox[{
41 RowBox[{"Rs", "=",
42 RowBox[{"0.32", "*", "lScale"}]}], ";"}], "\[IndentingNewLine]",
43 RowBox[{
44 RowBox[{"MgIIFrac", "=",
45 SuperscriptBox["10",
46 RowBox[{"-", "5"}]]}], ";"}], "\[IndentingNewLine]",
47 RowBox[{
48 RowBox[{"columnLength", "=",
49 RowBox[{"2.678", "*", "lScale"}]}], ";"}]}], "Input",
50 CellChangeTimes->{{3.695123264855075*^9, 3.695123266549275*^9}, {
51 3.695123397014388*^9, 3.695123397802359*^9}, {3.695123429692853*^9,
52 3.695123432030191*^9}, {3.6951235227564917`*^9, 3.695123523631727*^9}, {
53 3.695123597060893*^9, 3.69512360581527*^9}, 3.6951236467953587`*^9, {
54 3.6951236966551733`*^9, 3.695123731215706*^9}}],
55
56Cell[BoxData[
57 RowBox[{
58 RowBox[{"MgIIColDens", "=",
59 RowBox[{
60 FractionBox[
61 RowBox[{"\[Rho]avg", " ", "\[Rho]Scale"}], "mH"], "columnLength", " ",
62 "MgIIFrac"}]}], ";"}]], "Input",
63 CellChangeTimes->{{3.695123270283332*^9, 3.695123283532696*^9}, {
64 3.695123385162396*^9, 3.6951235105502367`*^9}, {3.695123552678598*^9,
65 3.695123585625667*^9}, {3.695123616595901*^9, 3.695123651805203*^9}}],
66
67Cell[CellGroupData[{
68
69Cell[BoxData[{
70 RowBox[{
71 RowBox[{"\[Rho]avg", "=",
72 RowBox[{"2.12", "*",
73 SuperscriptBox["10",
74 RowBox[{"-", "6"}]]}]}],
75 ";"}], "\[IndentingNewLine]", "MgIIColDens", "\[IndentingNewLine]",
76 RowBox[{"Clear", "[", "\[Rho]avg", "]"}]}], "Input",
77 CellChangeTimes->{{3.695123655311982*^9, 3.695123672164901*^9}, {
78 3.695123757727964*^9, 3.695123765276033*^9}, {3.6951260502954483`*^9,
79 3.695126050602833*^9}, {3.695126147236927*^9, 3.695126152419984*^9}}],
80
81Cell[BoxData["3.733826056682245`*^10"], "Output",
82 CellChangeTimes->{3.695123672572246*^9, 3.6951237666513977`*^9,
83 3.695126152705969*^9, 3.695475045781752*^9, 3.6967818983695*^9}]
84}, Open ]],
85
86Cell[CellGroupData[{
87
88Cell[BoxData[{
89 RowBox[{"density", "=",
90 RowBox[{"Solve", "[",
91 RowBox[{
92 RowBox[{"MgIIColDens", "\[Equal]",
93 RowBox[{"2", "*",
94 SuperscriptBox["10", "17"]}]}], ",", "\[Rho]avg"}],
95 "]"}]}], "\[IndentingNewLine]",
96 RowBox[{
97 RowBox[{"\[Rho]slope", "=",
98 RowBox[{"1.46804", "*",
99 SuperscriptBox["10",
100 RowBox[{"-", "8"}]]}]}], ";"}], "\[IndentingNewLine]",
101 RowBox[{
102 RowBox[{"time", "=",
103 RowBox[{"Solve", "[",
104 RowBox[{
105 RowBox[{
106 RowBox[{
107 RowBox[{"\[Rho]slope", " ", "cycles"}], "\[Equal]", "\[Rho]avg"}], "/.",
108 RowBox[{"density", "[",
109 RowBox[{"[", "1", "]"}], "]"}]}], ",", "cycles"}], "]"}]}],
110 ";"}], "\[IndentingNewLine]",
111 RowBox[{"years", "=",
112 RowBox[{
113 FractionBox[
114 RowBox[{"cycles", "/.",
115 RowBox[{"time", "[",
116 RowBox[{"[", "1", "]"}], "]"}]}], "10"],
117 FractionBox["timeScale",
118 RowBox[{"3.154", "*",
119 SuperscriptBox["10", "7"]}]]}]}]}], "Input",
120 CellChangeTimes->{{3.695123967538151*^9, 3.695124018436553*^9}, {
121 3.695125692478299*^9, 3.6951258414082327`*^9}, {3.6951260171382923`*^9,
122 3.695126066960822*^9}, {3.695126098978674*^9, 3.6951261395523777`*^9}, {
123 3.695475037043419*^9, 3.695475037312537*^9}, 3.696780748832787*^9,
124 3.696780779028981*^9, 3.696781891812632*^9}],
125
126Cell[BoxData[
127 RowBox[{"{",
128 RowBox[{"{",
129 RowBox[{"\[Rho]avg", "\[Rule]", "11.355644145264563`"}], "}"}],
130 "}"}]], "Output",
131 CellChangeTimes->{{3.695125812074988*^9, 3.69512581817631*^9},
132 3.695125848213621*^9, {3.6951260595902243`*^9, 3.6951261397915163`*^9}, {
133 3.695475039416697*^9, 3.695475045861493*^9}, {3.6967807448635883`*^9,
134 3.6967807492456284`*^9}, 3.6967807793213043`*^9, 3.696781899901896*^9}],
135
136Cell[BoxData["231112.2696699705`"], "Output",
137 CellChangeTimes->{{3.695125812074988*^9, 3.69512581817631*^9},
138 3.695125848213621*^9, {3.6951260595902243`*^9, 3.6951261397915163`*^9}, {
139 3.695475039416697*^9, 3.695475045861493*^9}, {3.6967807448635883`*^9,
140 3.6967807492456284`*^9}, 3.6967807793213043`*^9, 3.696781899908024*^9}]
141}, Open ]]
142},
143WindowSize->{1920, 1056},
144WindowMargins->{{0, Automatic}, {Automatic, 0}},
145FrontEndVersion->"11.0 for Linux x86 (64-bit) (July 28, 2016)",
146StyleDefinitions->"Default.nb"
147]
148(* End of Notebook Content *)
149
150(* Internal cache information *)
151(*CellTagsOutline
152CellTagsIndex->{}
153*)
154(*CellTagsIndex
155CellTagsIndex->{}
156*)
157(*NotebookFileOutline
158Notebook[{
159Cell[558, 20, 1270, 33, 185, "Input"],
160Cell[1831, 55, 409, 9, 59, "Input"],
161Cell[CellGroupData[{
162Cell[2265, 68, 472, 10, 77, "Input"],
163Cell[2740, 80, 182, 2, 34, "Output"]
164}, Open ]],
165Cell[CellGroupData[{
166Cell[2959, 87, 1306, 36, 144, "Input"],
167Cell[4268, 125, 423, 8, 32, "Output"],
168Cell[4694, 135, 336, 4, 30, "Output"]
169}, Open ]]
170}
171]
172*)
173