1 |
gezelter |
4195 |
(*CacheID: 234*) |
2 |
|
|
(* Internal cache information: |
3 |
|
|
NotebookFileLineBreakTest |
4 |
|
|
NotebookFileLineBreakTest |
5 |
|
|
NotebookDataPosition[ 0, 0] |
6 |
|
|
NotebookDataLength[ 19818, 604] |
7 |
|
|
NotebookOptionsPosition[ 18290, 543] |
8 |
|
|
NotebookOutlinePosition[ 18645, 559] |
9 |
|
|
CellTagsIndexPosition[ 18602, 556] |
10 |
|
|
WindowFrame->Normal*) |
11 |
|
|
|
12 |
|
|
(* Beginning of Notebook Content *) |
13 |
|
|
Notebook[{ |
14 |
|
|
Cell[BoxData[{ |
15 |
|
|
RowBox[{ |
16 |
|
|
RowBox[{"rc", "=."}], ";"}], "\[IndentingNewLine]", |
17 |
|
|
RowBox[{ |
18 |
|
|
RowBox[{"\[Kappa]", "=."}], ";"}], "\[IndentingNewLine]", |
19 |
|
|
RowBox[{ |
20 |
|
|
RowBox[{ |
21 |
|
|
RowBox[{"f", "[", |
22 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
23 |
|
|
RowBox[{ |
24 |
|
|
RowBox[{"Erfc", "[", |
25 |
|
|
RowBox[{"\[Kappa]", " ", "r"}], "]"}], "/", "r"}]}], |
26 |
|
|
";"}], "\[IndentingNewLine]", |
27 |
|
|
RowBox[{ |
28 |
|
|
RowBox[{ |
29 |
|
|
RowBox[{"g", "[", |
30 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
31 |
|
|
RowBox[{ |
32 |
|
|
RowBox[{ |
33 |
|
|
RowBox[{"(", |
34 |
|
|
RowBox[{ |
35 |
|
|
RowBox[{"-", "2"}], "*", "\[Kappa]"}], ")"}], "/", |
36 |
|
|
RowBox[{"(", |
37 |
|
|
RowBox[{ |
38 |
|
|
RowBox[{"E", "^", |
39 |
|
|
RowBox[{"(", |
40 |
|
|
RowBox[{ |
41 |
|
|
RowBox[{"r", "^", "2"}], "*", |
42 |
|
|
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
43 |
|
|
RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-", |
44 |
|
|
RowBox[{ |
45 |
|
|
RowBox[{"Erfc", "[", |
46 |
|
|
RowBox[{"r", "*", "\[Kappa]"}], "]"}], "/", |
47 |
|
|
RowBox[{"r", "^", "2"}]}]}]}], " ", ";"}], "\[IndentingNewLine]", |
48 |
|
|
RowBox[{ |
49 |
|
|
RowBox[{ |
50 |
|
|
RowBox[{"h", "[", |
51 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", |
52 |
|
|
RowBox[{ |
53 |
|
|
RowBox[{ |
54 |
|
|
RowBox[{"(", |
55 |
|
|
RowBox[{"4", "*", "\[Kappa]"}], ")"}], "/", |
56 |
|
|
RowBox[{"(", |
57 |
|
|
RowBox[{ |
58 |
|
|
RowBox[{"E", "^", |
59 |
|
|
RowBox[{"(", |
60 |
|
|
RowBox[{ |
61 |
|
|
RowBox[{"r", "^", "2"}], "*", |
62 |
|
|
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
63 |
|
|
RowBox[{"Sqrt", "[", "Pi", "]"}], "*", |
64 |
|
|
RowBox[{"r", "^", "2"}]}], ")"}]}], "+", |
65 |
|
|
RowBox[{ |
66 |
|
|
RowBox[{"(", |
67 |
|
|
RowBox[{"4", "*", |
68 |
|
|
RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/", |
69 |
|
|
RowBox[{"(", |
70 |
|
|
RowBox[{ |
71 |
|
|
RowBox[{"E", "^", |
72 |
|
|
RowBox[{"(", |
73 |
|
|
RowBox[{ |
74 |
|
|
RowBox[{"r", "^", "2"}], "*", |
75 |
|
|
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
76 |
|
|
RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+", |
77 |
|
|
RowBox[{ |
78 |
|
|
RowBox[{"(", |
79 |
|
|
RowBox[{"2", "*", |
80 |
|
|
RowBox[{"Erfc", "[", |
81 |
|
|
RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/", |
82 |
|
|
RowBox[{"r", "^", "3"}]}]}]}], ";"}], "\[IndentingNewLine]", |
83 |
|
|
RowBox[{ |
84 |
|
|
RowBox[{ |
85 |
|
|
RowBox[{"f0", "[", |
86 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
87 |
|
|
RowBox[{ |
88 |
|
|
RowBox[{"f", "[", |
89 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", |
90 |
|
|
RowBox[{"f", "[", |
91 |
|
|
RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-", |
92 |
|
|
RowBox[{ |
93 |
|
|
RowBox[{"(", |
94 |
|
|
RowBox[{"r", "-", "rc"}], ")"}], |
95 |
|
|
RowBox[{"g", "[", |
96 |
|
|
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], " ", |
97 |
|
|
";"}], "\[IndentingNewLine]", |
98 |
|
|
RowBox[{ |
99 |
|
|
RowBox[{ |
100 |
|
|
RowBox[{"g0", "[", |
101 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
102 |
|
|
RowBox[{ |
103 |
|
|
RowBox[{"g", "[", |
104 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", |
105 |
|
|
RowBox[{"g", "[", |
106 |
|
|
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", |
107 |
|
|
RowBox[{ |
108 |
|
|
RowBox[{ |
109 |
|
|
RowBox[{"h0", "[", |
110 |
|
|
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", |
111 |
|
|
RowBox[{"h", "[", |
112 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], ";"}], "\[IndentingNewLine]", |
113 |
|
|
RowBox[{ |
114 |
|
|
RowBox[{ |
115 |
|
|
RowBox[{"vprime", "[", |
116 |
|
|
RowBox[{"r_", ",", "method_"}], "]"}], ":=", |
117 |
|
|
RowBox[{"Switch", "[", |
118 |
|
|
RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", |
119 |
|
|
RowBox[{"g", "[", |
120 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
121 |
|
|
"TSF", ",", " ", |
122 |
|
|
RowBox[{"g0", "[", |
123 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", "SP", |
124 |
|
|
",", " ", |
125 |
|
|
RowBox[{"g", "[", |
126 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", |
127 |
|
|
"GSF", ",", " ", |
128 |
|
|
RowBox[{ |
129 |
|
|
RowBox[{"g", "[", |
130 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], " ", "-", " ", |
131 |
|
|
RowBox[{"g", "[", |
132 |
|
|
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], "\[IndentingNewLine]", |
133 |
|
|
"]"}]}], ";"}], "\[IndentingNewLine]", |
134 |
|
|
RowBox[{ |
135 |
|
|
RowBox[{ |
136 |
|
|
RowBox[{"vprimeprime", "[", |
137 |
|
|
RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ", |
138 |
|
|
RowBox[{"Switch", "[", |
139 |
|
|
RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",", |
140 |
|
|
RowBox[{"h", "[", |
141 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
142 |
|
|
"TSF", ",", " ", |
143 |
|
|
RowBox[{"h0", "[", |
144 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
145 |
|
|
"SP", ",", " ", |
146 |
|
|
RowBox[{"h", "[", |
147 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
148 |
|
|
"GSF", ",", " ", |
149 |
|
|
RowBox[{"h", "[", |
150 |
|
|
RowBox[{"r", ",", "\[Kappa]"}], "]"}]}], "\[IndentingNewLine]", "]"}]}], |
151 |
|
|
";"}], "\[IndentingNewLine]", |
152 |
|
|
RowBox[{ |
153 |
|
|
RowBox[{ |
154 |
|
|
RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=", |
155 |
|
|
RowBox[{"r", " ", |
156 |
|
|
RowBox[{"Switch", "[", |
157 |
|
|
RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",", |
158 |
|
|
RowBox[{ |
159 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
160 |
|
|
RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]", |
161 |
|
|
"2", ",", |
162 |
|
|
RowBox[{ |
163 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
164 |
|
|
RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]", |
165 |
|
|
"3", ",", |
166 |
|
|
RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}], |
167 |
|
|
";"}], "\[IndentingNewLine]", |
168 |
|
|
RowBox[{ |
169 |
|
|
RowBox[{ |
170 |
|
|
RowBox[{"T", "[", |
171 |
|
|
RowBox[{"\[Alpha]_", ",", "\[Beta]_", ",", " ", "method_"}], "]"}], " ", ":=", |
172 |
|
|
" ", |
173 |
|
|
RowBox[{ |
174 |
|
|
RowBox[{ |
175 |
|
|
RowBox[{"KroneckerDelta", "[", |
176 |
|
|
RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*", " ", |
177 |
|
|
RowBox[{"(", |
178 |
|
|
RowBox[{ |
179 |
|
|
RowBox[{"vprime", "[", |
180 |
|
|
RowBox[{"r", ",", " ", "method"}], "]"}], "/", "r"}], ")"}]}], "+", |
181 |
|
|
" ", |
182 |
|
|
RowBox[{ |
183 |
|
|
RowBox[{"R", "[", "\[Alpha]", "]"}], |
184 |
|
|
RowBox[{"R", "[", "\[Beta]", "]"}], |
185 |
|
|
RowBox[{ |
186 |
|
|
RowBox[{"(", |
187 |
|
|
RowBox[{ |
188 |
|
|
RowBox[{"vprimeprime", "[", |
189 |
|
|
RowBox[{"r", ",", "method"}], "]"}], " ", "-", " ", |
190 |
|
|
RowBox[{ |
191 |
|
|
RowBox[{"vprime", "[", |
192 |
|
|
RowBox[{"r", ",", "method"}], "]"}], "/", "r"}]}], ")"}], "/", |
193 |
|
|
RowBox[{"r", "^", "2"}]}]}]}]}], ";"}], "\[IndentingNewLine]", |
194 |
|
|
RowBox[{ |
195 |
|
|
RowBox[{"TBare", " ", "=", " ", |
196 |
|
|
RowBox[{"MatrixForm", "[", |
197 |
|
|
RowBox[{"Table", "[", |
198 |
|
|
RowBox[{ |
199 |
|
|
RowBox[{"Table", "[", |
200 |
|
|
RowBox[{ |
201 |
|
|
RowBox[{"Integrate", "[", |
202 |
|
|
RowBox[{ |
203 |
|
|
RowBox[{ |
204 |
|
|
RowBox[{"r", "^", "2"}], |
205 |
|
|
RowBox[{"Integrate", "[", |
206 |
|
|
RowBox[{ |
207 |
|
|
RowBox[{"Integrate", "[", |
208 |
|
|
RowBox[{ |
209 |
|
|
RowBox[{ |
210 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
211 |
|
|
RowBox[{"T", "[", |
212 |
|
|
RowBox[{"i", ",", "j", ",", "BARE"}], "]"}]}], ",", |
213 |
|
|
RowBox[{"{", |
214 |
|
|
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
215 |
|
|
",", |
216 |
|
|
RowBox[{"{", |
217 |
|
|
RowBox[{"\[Phi]", ",", "0", ",", |
218 |
|
|
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
219 |
|
|
RowBox[{"{", |
220 |
|
|
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
221 |
|
|
RowBox[{"{", |
222 |
|
|
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
223 |
|
|
RowBox[{"{", |
224 |
|
|
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
225 |
|
|
";"}], "\[IndentingNewLine]", |
226 |
|
|
RowBox[{ |
227 |
|
|
RowBox[{"TSP", " ", "=", " ", |
228 |
|
|
RowBox[{"MatrixForm", "[", |
229 |
|
|
RowBox[{"Table", "[", |
230 |
|
|
RowBox[{ |
231 |
|
|
RowBox[{"Table", "[", |
232 |
|
|
RowBox[{ |
233 |
|
|
RowBox[{"Integrate", "[", |
234 |
|
|
RowBox[{ |
235 |
|
|
RowBox[{ |
236 |
|
|
RowBox[{"r", "^", "2"}], |
237 |
|
|
RowBox[{"Integrate", "[", |
238 |
|
|
RowBox[{ |
239 |
|
|
RowBox[{"Integrate", "[", |
240 |
|
|
RowBox[{ |
241 |
|
|
RowBox[{ |
242 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
243 |
|
|
RowBox[{"T", "[", |
244 |
|
|
RowBox[{"i", ",", "j", ",", "SP"}], "]"}]}], ",", |
245 |
|
|
RowBox[{"{", |
246 |
|
|
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
247 |
|
|
",", |
248 |
|
|
RowBox[{"{", |
249 |
|
|
RowBox[{"\[Phi]", ",", "0", ",", |
250 |
|
|
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
251 |
|
|
RowBox[{"{", |
252 |
|
|
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
253 |
|
|
RowBox[{"{", |
254 |
|
|
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
255 |
|
|
RowBox[{"{", |
256 |
|
|
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
257 |
|
|
";"}], "\[IndentingNewLine]", |
258 |
|
|
RowBox[{ |
259 |
|
|
RowBox[{"TGSF", " ", "=", " ", |
260 |
|
|
RowBox[{"MatrixForm", "[", |
261 |
|
|
RowBox[{"Table", "[", |
262 |
|
|
RowBox[{ |
263 |
|
|
RowBox[{"Table", "[", |
264 |
|
|
RowBox[{ |
265 |
|
|
RowBox[{"Integrate", "[", |
266 |
|
|
RowBox[{ |
267 |
|
|
RowBox[{ |
268 |
|
|
RowBox[{"r", "^", "2"}], |
269 |
|
|
RowBox[{"Integrate", "[", |
270 |
|
|
RowBox[{ |
271 |
|
|
RowBox[{"Integrate", "[", |
272 |
|
|
RowBox[{ |
273 |
|
|
RowBox[{ |
274 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
275 |
|
|
RowBox[{"T", "[", |
276 |
|
|
RowBox[{"i", ",", "j", ",", "GSF"}], "]"}]}], ",", |
277 |
|
|
RowBox[{"{", |
278 |
|
|
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
279 |
|
|
",", |
280 |
|
|
RowBox[{"{", |
281 |
|
|
RowBox[{"\[Phi]", ",", "0", ",", |
282 |
|
|
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
283 |
|
|
RowBox[{"{", |
284 |
|
|
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
285 |
|
|
RowBox[{"{", |
286 |
|
|
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
287 |
|
|
RowBox[{"{", |
288 |
|
|
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
289 |
|
|
";"}], "\[IndentingNewLine]", |
290 |
|
|
RowBox[{ |
291 |
|
|
RowBox[{"TTSF", " ", "=", " ", |
292 |
|
|
RowBox[{"MatrixForm", "[", |
293 |
|
|
RowBox[{"Table", "[", |
294 |
|
|
RowBox[{ |
295 |
|
|
RowBox[{"Table", "[", |
296 |
|
|
RowBox[{ |
297 |
|
|
RowBox[{"Integrate", "[", |
298 |
|
|
RowBox[{ |
299 |
|
|
RowBox[{ |
300 |
|
|
RowBox[{"r", "^", "2"}], |
301 |
|
|
RowBox[{"Integrate", "[", |
302 |
|
|
RowBox[{ |
303 |
|
|
RowBox[{"Integrate", "[", |
304 |
|
|
RowBox[{ |
305 |
|
|
RowBox[{ |
306 |
|
|
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
307 |
|
|
RowBox[{"T", "[", |
308 |
|
|
RowBox[{"i", ",", "j", ",", "TSF"}], "]"}]}], ",", |
309 |
|
|
RowBox[{"{", |
310 |
|
|
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
311 |
|
|
",", |
312 |
|
|
RowBox[{"{", |
313 |
|
|
RowBox[{"\[Phi]", ",", "0", ",", |
314 |
|
|
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
315 |
|
|
RowBox[{"{", |
316 |
|
|
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
317 |
|
|
RowBox[{"{", |
318 |
|
|
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
319 |
|
|
RowBox[{"{", |
320 |
|
|
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
321 |
|
|
";"}], "\[IndentingNewLine]", |
322 |
|
|
RowBox[{ |
323 |
|
|
RowBox[{"QBARE", " ", "=", |
324 |
|
|
RowBox[{"3", |
325 |
|
|
RowBox[{ |
326 |
|
|
RowBox[{"Part", "[", |
327 |
|
|
RowBox[{"TBare", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
328 |
|
|
RowBox[{"(", |
329 |
|
|
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
330 |
|
|
RowBox[{ |
331 |
|
|
RowBox[{"QTSP", " ", "=", " ", |
332 |
|
|
RowBox[{"3", |
333 |
|
|
RowBox[{ |
334 |
|
|
RowBox[{"Part", "[", |
335 |
|
|
RowBox[{"TSP", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
336 |
|
|
RowBox[{"(", |
337 |
|
|
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
338 |
|
|
RowBox[{ |
339 |
|
|
RowBox[{"QGSF", " ", "=", " ", |
340 |
|
|
RowBox[{"3", |
341 |
|
|
RowBox[{ |
342 |
|
|
RowBox[{"Part", "[", |
343 |
|
|
RowBox[{"TGSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
344 |
|
|
RowBox[{"(", |
345 |
|
|
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
346 |
|
|
RowBox[{ |
347 |
|
|
RowBox[{"QTSF", " ", "=", " ", |
348 |
|
|
RowBox[{"3", |
349 |
|
|
RowBox[{ |
350 |
|
|
RowBox[{"Part", "[", |
351 |
|
|
RowBox[{"TTSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
352 |
|
|
RowBox[{"(", |
353 |
|
|
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}]}], "Input", |
354 |
|
|
CellChangeTimes->{{3.613731558054296*^9, 3.613732200209785*^9}, |
355 |
|
|
3.6137322781282454`*^9, {3.61373233169027*^9, 3.613732397354877*^9}, { |
356 |
|
|
3.613732488348338*^9, 3.6137327713906393`*^9}, {3.613732810835518*^9, |
357 |
|
|
3.613732821957967*^9}, {3.613732862763562*^9, 3.613732913141625*^9}, { |
358 |
|
|
3.61373295595298*^9, 3.613733106904414*^9}, {3.613733177027177*^9, |
359 |
|
|
3.613733199740512*^9}, {3.61373326352699*^9, 3.6137332727080936`*^9}, { |
360 |
|
|
3.613733425731758*^9, 3.6137335719815893`*^9}, {3.613733605778666*^9, |
361 |
|
|
3.613733648768111*^9}, {3.613733776819376*^9, 3.6137338150569887`*^9}, { |
362 |
|
|
3.6137338699918118`*^9, 3.61373393057332*^9}, {3.613733985426153*^9, |
363 |
|
|
3.613733985628289*^9}, {3.613734054134171*^9, 3.6137340588514023`*^9}, { |
364 |
|
|
3.6137371511302147`*^9, 3.613737188118826*^9}, {3.613815304031159*^9, |
365 |
|
|
3.613815306899302*^9}, {3.613815415274508*^9, 3.613815450160325*^9}, { |
366 |
|
|
3.613815535717202*^9, 3.613815604127789*^9}, {3.613815665097225*^9, |
367 |
|
|
3.613815673383039*^9}, {3.613816023942026*^9, 3.6138160527230864`*^9}, { |
368 |
|
|
3.613825536655637*^9, 3.613825555245915*^9}, {3.613825660601698*^9, |
369 |
|
|
3.6138256692518787`*^9}, {3.613825968874196*^9, 3.613826000721004*^9}, { |
370 |
|
|
3.6138263115515537`*^9, 3.613826312408416*^9}, {3.613826919678602*^9, |
371 |
|
|
3.613826920319913*^9}, {3.613827178556149*^9, 3.6138272197490597`*^9}, { |
372 |
|
|
3.6138410090125504`*^9, 3.613841090576209*^9}, {3.613841142597171*^9, |
373 |
|
|
3.613841191624243*^9}, {3.6138412274048777`*^9, 3.613841229484283*^9}, { |
374 |
|
|
3.6138413897003736`*^9, 3.613841438224633*^9}, {3.61384150925097*^9, |
375 |
|
|
3.6138415133903227`*^9}, {3.613841547354535*^9, 3.613841573983226*^9}, { |
376 |
|
|
3.613898318732862*^9, 3.613898326606717*^9}, {3.614594975804741*^9, |
377 |
|
|
3.614595013097396*^9}, {3.61459504604961*^9, 3.614595067267371*^9}}], |
378 |
|
|
|
379 |
|
|
Cell[CellGroupData[{ |
380 |
|
|
|
381 |
|
|
Cell[BoxData["QBARE"], "Input", |
382 |
|
|
CellChangeTimes->{{3.61389845863691*^9, 3.6138984602794447`*^9}}], |
383 |
|
|
|
384 |
|
|
Cell[BoxData[ |
385 |
|
|
FractionBox[ |
386 |
|
|
RowBox[{"3", " ", |
387 |
|
|
RowBox[{"(", |
388 |
|
|
RowBox[{ |
389 |
|
|
RowBox[{ |
390 |
|
|
RowBox[{"-", |
391 |
|
|
FractionBox["8", "3"]}], " ", |
392 |
|
|
SuperscriptBox["\[ExponentialE]", |
393 |
|
|
RowBox[{ |
394 |
|
|
RowBox[{"-", |
395 |
|
|
SuperscriptBox["rc", "2"]}], " ", |
396 |
|
|
SuperscriptBox["\[Kappa]", "2"]}]], " ", |
397 |
|
|
SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+", |
398 |
|
|
RowBox[{ |
399 |
|
|
FractionBox["4", "3"], " ", "\[Pi]", " ", |
400 |
|
|
RowBox[{"Erf", "[", |
401 |
|
|
RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], |
402 |
|
|
RowBox[{"4", " ", "\[Pi]"}]]], "Output", |
403 |
|
|
CellChangeTimes->{3.61389846058971*^9, 3.61459514866705*^9}] |
404 |
|
|
}, Open ]], |
405 |
|
|
|
406 |
|
|
Cell[CellGroupData[{ |
407 |
|
|
|
408 |
|
|
Cell[BoxData["QTSP"], "Input", |
409 |
|
|
CellChangeTimes->{{3.613898461774685*^9, 3.61389846255107*^9}}], |
410 |
|
|
|
411 |
|
|
Cell[BoxData[ |
412 |
|
|
FractionBox[ |
413 |
|
|
RowBox[{"3", " ", |
414 |
|
|
RowBox[{"(", |
415 |
|
|
RowBox[{ |
416 |
|
|
RowBox[{ |
417 |
|
|
RowBox[{"-", |
418 |
|
|
FractionBox["8", "3"]}], " ", |
419 |
|
|
SuperscriptBox["\[ExponentialE]", |
420 |
|
|
RowBox[{ |
421 |
|
|
RowBox[{"-", |
422 |
|
|
SuperscriptBox["rc", "2"]}], " ", |
423 |
|
|
SuperscriptBox["\[Kappa]", "2"]}]], " ", |
424 |
|
|
SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+", |
425 |
|
|
RowBox[{ |
426 |
|
|
FractionBox["4", "3"], " ", "\[Pi]", " ", |
427 |
|
|
RowBox[{"Erf", "[", |
428 |
|
|
RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], |
429 |
|
|
RowBox[{"4", " ", "\[Pi]"}]]], "Output", |
430 |
|
|
CellChangeTimes->{3.6138984628297167`*^9, 3.614595151448133*^9}] |
431 |
|
|
}, Open ]], |
432 |
|
|
|
433 |
|
|
Cell[CellGroupData[{ |
434 |
|
|
|
435 |
|
|
Cell[BoxData["QGSF"], "Input", |
436 |
|
|
CellChangeTimes->{{3.613898464699469*^9, 3.613898465498917*^9}}], |
437 |
|
|
|
438 |
|
|
Cell[BoxData["1"], "Output", |
439 |
|
|
CellChangeTimes->{ |
440 |
|
|
3.613898465771669*^9, {3.614595155790647*^9, 3.61459516007798*^9}}] |
441 |
|
|
}, Open ]], |
442 |
|
|
|
443 |
|
|
Cell[CellGroupData[{ |
444 |
|
|
|
445 |
|
|
Cell[BoxData["QTSF"], "Input", |
446 |
|
|
CellChangeTimes->{{3.613898466555105*^9, 3.613898468075191*^9}}], |
447 |
|
|
|
448 |
|
|
Cell[BoxData["1"], "Output", |
449 |
|
|
CellChangeTimes->{3.613898468428026*^9, 3.614595163440785*^9}] |
450 |
|
|
}, Open ]], |
451 |
|
|
|
452 |
|
|
Cell[CellGroupData[{ |
453 |
|
|
|
454 |
|
|
Cell[BoxData[ |
455 |
|
|
RowBox[{"FullSimplify", "[", "QTSP", "]"}]], "Input", |
456 |
|
|
CellChangeTimes->{{3.613898477301433*^9, 3.613898481802129*^9}}], |
457 |
|
|
|
458 |
|
|
Cell[BoxData[ |
459 |
|
|
RowBox[{ |
460 |
|
|
RowBox[{"-", |
461 |
|
|
FractionBox[ |
462 |
|
|
RowBox[{"2", " ", |
463 |
|
|
SuperscriptBox["\[ExponentialE]", |
464 |
|
|
RowBox[{ |
465 |
|
|
RowBox[{"-", |
466 |
|
|
SuperscriptBox["rc", "2"]}], " ", |
467 |
|
|
SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}], |
468 |
|
|
SqrtBox["\[Pi]"]]}], "+", |
469 |
|
|
RowBox[{"Erf", "[", |
470 |
|
|
RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output", |
471 |
|
|
CellChangeTimes->{3.613898482310692*^9, 3.614595171545216*^9}] |
472 |
|
|
}, Open ]], |
473 |
|
|
|
474 |
|
|
Cell[CellGroupData[{ |
475 |
|
|
|
476 |
|
|
Cell[BoxData[ |
477 |
|
|
RowBox[{"FullSimplify", "[", "QBARE", "]"}]], "Input", |
478 |
|
|
CellChangeTimes->{{3.6145951775304317`*^9, 3.6145951844023333`*^9}}], |
479 |
|
|
|
480 |
|
|
Cell[BoxData[ |
481 |
|
|
RowBox[{ |
482 |
|
|
RowBox[{"-", |
483 |
|
|
FractionBox[ |
484 |
|
|
RowBox[{"2", " ", |
485 |
|
|
SuperscriptBox["\[ExponentialE]", |
486 |
|
|
RowBox[{ |
487 |
|
|
RowBox[{"-", |
488 |
|
|
SuperscriptBox["rc", "2"]}], " ", |
489 |
|
|
SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}], |
490 |
|
|
SqrtBox["\[Pi]"]]}], "+", |
491 |
|
|
RowBox[{"Erf", "[", |
492 |
|
|
RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output", |
493 |
|
|
CellChangeTimes->{3.614595184646675*^9}] |
494 |
|
|
}, Open ]], |
495 |
|
|
|
496 |
|
|
Cell[CellGroupData[{ |
497 |
|
|
|
498 |
|
|
Cell[BoxData[ |
499 |
|
|
RowBox[{"Limit", "[", |
500 |
|
|
RowBox[{"QBARE", ",", |
501 |
|
|
RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input", |
502 |
|
|
CellChangeTimes->{{3.6145952058026037`*^9, 3.6145952197375603`*^9}}], |
503 |
|
|
|
504 |
|
|
Cell[BoxData["0"], "Output", |
505 |
|
|
CellChangeTimes->{3.6145952203235397`*^9}] |
506 |
|
|
}, Open ]], |
507 |
|
|
|
508 |
|
|
Cell[CellGroupData[{ |
509 |
|
|
|
510 |
|
|
Cell[BoxData[ |
511 |
|
|
RowBox[{"Limit", "[", |
512 |
|
|
RowBox[{"QTSP", ",", |
513 |
|
|
RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input", |
514 |
|
|
CellChangeTimes->{{3.614595380663261*^9, 3.614595396297574*^9}}], |
515 |
|
|
|
516 |
|
|
Cell[BoxData["0"], "Output", |
517 |
|
|
CellChangeTimes->{{3.614595392877947*^9, 3.614595396693039*^9}}] |
518 |
|
|
}, Open ]], |
519 |
|
|
|
520 |
|
|
Cell[CellGroupData[{ |
521 |
|
|
|
522 |
|
|
Cell[BoxData[ |
523 |
|
|
RowBox[{"Limit", "[", |
524 |
|
|
RowBox[{"QGSF", ",", |
525 |
|
|
RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input", |
526 |
|
|
CellChangeTimes->{{3.6145953985545464`*^9, 3.614595405611649*^9}}], |
527 |
|
|
|
528 |
|
|
Cell[BoxData["1"], "Output", |
529 |
|
|
CellChangeTimes->{3.614595405937716*^9}] |
530 |
|
|
}, Open ]], |
531 |
|
|
|
532 |
|
|
Cell[CellGroupData[{ |
533 |
|
|
|
534 |
|
|
Cell[BoxData[ |
535 |
|
|
RowBox[{"Limit", "[", |
536 |
|
|
RowBox[{"QTSF", ",", |
537 |
|
|
RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input", |
538 |
|
|
CellChangeTimes->{{3.6145954076554337`*^9, 3.614595415455311*^9}}], |
539 |
|
|
|
540 |
|
|
Cell[BoxData["1"], "Output", |
541 |
|
|
CellChangeTimes->{3.6145954158135767`*^9}] |
542 |
|
|
}, Open ]] |
543 |
|
|
}, |
544 |
|
|
WindowSize->{884, 1337}, |
545 |
|
|
WindowMargins->{{Automatic, 292}, {Automatic, 0}}, |
546 |
|
|
FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \ |
547 |
|
|
2014)", |
548 |
|
|
StyleDefinitions->"Default.nb" |
549 |
|
|
] |
550 |
|
|
(* End of Notebook Content *) |
551 |
|
|
|
552 |
|
|
(* Internal cache information *) |
553 |
|
|
(*CellTagsOutline |
554 |
|
|
CellTagsIndex->{} |
555 |
|
|
*) |
556 |
|
|
(*CellTagsIndex |
557 |
|
|
CellTagsIndex->{} |
558 |
|
|
*) |
559 |
|
|
(*NotebookFileOutline |
560 |
|
|
Notebook[{ |
561 |
|
|
Cell[400, 13, 13403, 363, 811, "Input"], |
562 |
|
|
Cell[CellGroupData[{ |
563 |
|
|
Cell[13828, 380, 98, 1, 28, "Input"], |
564 |
|
|
Cell[13929, 383, 639, 19, 58, "Output"] |
565 |
|
|
}, Open ]], |
566 |
|
|
Cell[CellGroupData[{ |
567 |
|
|
Cell[14605, 407, 95, 1, 28, "Input"], |
568 |
|
|
Cell[14703, 410, 643, 19, 58, "Output"] |
569 |
|
|
}, Open ]], |
570 |
|
|
Cell[CellGroupData[{ |
571 |
|
|
Cell[15383, 434, 96, 1, 28, "Input"], |
572 |
|
|
Cell[15482, 437, 118, 2, 28, "Output"] |
573 |
|
|
}, Open ]], |
574 |
|
|
Cell[CellGroupData[{ |
575 |
|
|
Cell[15637, 444, 96, 1, 28, "Input"], |
576 |
|
|
Cell[15736, 447, 92, 1, 28, "Output"] |
577 |
|
|
}, Open ]], |
578 |
|
|
Cell[CellGroupData[{ |
579 |
|
|
Cell[15865, 453, 134, 2, 28, "Input"], |
580 |
|
|
Cell[16002, 457, 445, 13, 58, "Output"] |
581 |
|
|
}, Open ]], |
582 |
|
|
Cell[CellGroupData[{ |
583 |
|
|
Cell[16484, 475, 139, 2, 28, "Input"], |
584 |
|
|
Cell[16626, 479, 423, 13, 58, "Output"] |
585 |
|
|
}, Open ]], |
586 |
|
|
Cell[CellGroupData[{ |
587 |
|
|
Cell[17086, 497, 192, 4, 28, "Input"], |
588 |
|
|
Cell[17281, 503, 72, 1, 28, "Output"] |
589 |
|
|
}, Open ]], |
590 |
|
|
Cell[CellGroupData[{ |
591 |
|
|
Cell[17390, 509, 187, 4, 28, "Input"], |
592 |
|
|
Cell[17580, 515, 94, 1, 28, "Output"] |
593 |
|
|
}, Open ]], |
594 |
|
|
Cell[CellGroupData[{ |
595 |
|
|
Cell[17711, 521, 189, 4, 28, "Input"], |
596 |
|
|
Cell[17903, 527, 70, 1, 28, "Output"] |
597 |
|
|
}, Open ]], |
598 |
|
|
Cell[CellGroupData[{ |
599 |
|
|
Cell[18010, 533, 189, 4, 28, "Input"], |
600 |
|
|
Cell[18202, 539, 72, 1, 28, "Output"] |
601 |
|
|
}, Open ]] |
602 |
|
|
} |
603 |
|
|
] |
604 |
|
|
*) |
605 |
|
|
|