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:
|
---|
10 | NotebookFileLineBreakTest
|
---|
11 | NotebookFileLineBreakTest
|
---|
12 | NotebookDataPosition[ 158, 7]
|
---|
13 | NotebookDataLength[ 5575, 165]
|
---|
14 | NotebookOptionsPosition[ 5046, 142]
|
---|
15 | NotebookOutlinePosition[ 5380, 157]
|
---|
16 | CellTagsIndexPosition[ 5337, 154]
|
---|
17 | WindowFrame->Normal*)
|
---|
18 |
|
---|
19 | (* Beginning of Notebook Content *)
|
---|
20 | Notebook[{
|
---|
21 | Cell[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 |
|
---|
56 | Cell[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 |
|
---|
67 | Cell[CellGroupData[{
|
---|
68 |
|
---|
69 | Cell[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 |
|
---|
81 | Cell[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 |
|
---|
86 | Cell[CellGroupData[{
|
---|
87 |
|
---|
88 | Cell[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 |
|
---|
126 | Cell[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 |
|
---|
136 | Cell[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 | },
|
---|
143 | WindowSize->{1920, 1056},
|
---|
144 | WindowMargins->{{0, Automatic}, {Automatic, 0}},
|
---|
145 | FrontEndVersion->"11.0 for Linux x86 (64-bit) (July 28, 2016)",
|
---|
146 | StyleDefinitions->"Default.nb"
|
---|
147 | ]
|
---|
148 | (* End of Notebook Content *)
|
---|
149 |
|
---|
150 | (* Internal cache information *)
|
---|
151 | (*CellTagsOutline
|
---|
152 | CellTagsIndex->{}
|
---|
153 | *)
|
---|
154 | (*CellTagsIndex
|
---|
155 | CellTagsIndex->{}
|
---|
156 | *)
|
---|
157 | (*NotebookFileOutline
|
---|
158 | Notebook[{
|
---|
159 | Cell[558, 20, 1270, 33, 185, "Input"],
|
---|
160 | Cell[1831, 55, 409, 9, 59, "Input"],
|
---|
161 | Cell[CellGroupData[{
|
---|
162 | Cell[2265, 68, 472, 10, 77, "Input"],
|
---|
163 | Cell[2740, 80, 182, 2, 34, "Output"]
|
---|
164 | }, Open ]],
|
---|
165 | Cell[CellGroupData[{
|
---|
166 | Cell[2959, 87, 1306, 36, 144, "Input"],
|
---|
167 | Cell[4268, 125, 423, 8, 32, "Output"],
|
---|
168 | Cell[4694, 135, 336, 4, 30, "Output"]
|
---|
169 | }, Open ]]
|
---|
170 | }
|
---|
171 | ]
|
---|
172 | *)
|
---|
173 |
|
---|