ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/multipole/DipoleTensor.nb
Revision: 4195
Committed: Thu Jul 17 18:24:41 2014 UTC (9 years, 11 months ago) by gezelter
File size: 26212 byte(s)
Log Message: 
Most recent version including dielectric stuff

File Contents

# Content
1 (*CacheID: 234*)
2 (* Internal cache information:
3 NotebookFileLineBreakTest
4 NotebookFileLineBreakTest
5 NotebookDataPosition[ 0, 0]
6 NotebookDataLength[ 26211, 815]
7 NotebookOptionsPosition[ 24416, 745]
8 NotebookOutlinePosition[ 24771, 761]
9 CellTagsIndexPosition[ 24728, 758]
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[{"s", "[",
86 RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
87 RowBox[{
88 RowBox[{
89 RowBox[{"(",
90 RowBox[{
91 RowBox[{"-", "12"}], "*", "\[Kappa]"}], ")"}], "/",
92 RowBox[{"(",
93 RowBox[{
94 RowBox[{"E", "^",
95 RowBox[{"(",
96 RowBox[{
97 RowBox[{"r", "^", "2"}], "*",
98 RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
99 RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
100 RowBox[{"r", "^", "3"}]}], ")"}]}], "-",
101 RowBox[{
102 RowBox[{"(",
103 RowBox[{"8", "*",
104 RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
105 RowBox[{"(",
106 RowBox[{
107 RowBox[{"E", "^",
108 RowBox[{"(",
109 RowBox[{
110 RowBox[{"r", "^", "2"}], "*",
111 RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
112 RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-",
113 RowBox[{
114 RowBox[{"(",
115 RowBox[{"8", "*", "r", "*",
116 RowBox[{"\[Kappa]", "^", "5"}]}], ")"}], "/",
117 RowBox[{"(",
118 RowBox[{
119 RowBox[{"E", "^",
120 RowBox[{"(",
121 RowBox[{
122 RowBox[{"r", "^", "2"}], "*",
123 RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
124 RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "-",
125 RowBox[{
126 RowBox[{"(",
127 RowBox[{"6", "*",
128 RowBox[{"Erfc", "[",
129 RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
130 RowBox[{"r", "^", "4"}]}]}]}], ";"}], "\[IndentingNewLine]",
131 RowBox[{
132 RowBox[{
133 RowBox[{"f2", "[",
134 RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
135 RowBox[{
136 RowBox[{"f", "[",
137 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
138 RowBox[{"f", "[",
139 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
140 RowBox[{
141 RowBox[{"(",
142 RowBox[{"r", "-", "rc"}], ")"}],
143 RowBox[{"g", "[",
144 RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
145 RowBox[{
146 RowBox[{
147 RowBox[{"(",
148 RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}], " ",
149 RowBox[{
150 RowBox[{"h", "[",
151 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}], " ", "-", " ",
152 RowBox[{
153 RowBox[{
154 RowBox[{"(",
155 RowBox[{"r", "-", "rc"}], ")"}], "^", "3"}], " ",
156 RowBox[{
157 RowBox[{"s", "[",
158 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "6"}]}]}]}],
159 ";"}], "\[IndentingNewLine]",
160 RowBox[{
161 RowBox[{
162 RowBox[{"g2", "[",
163 RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
164 RowBox[{
165 RowBox[{"g", "[",
166 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
167 RowBox[{"g", "[",
168 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
169 RowBox[{
170 RowBox[{"(",
171 RowBox[{"r", "-", "rc"}], ")"}],
172 RowBox[{"h", "[",
173 RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
174 RowBox[{
175 RowBox[{
176 RowBox[{"(",
177 RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}], " ",
178 RowBox[{
179 RowBox[{"s", "[",
180 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}]}]}], " ",
181 ";"}], "\[IndentingNewLine]",
182 RowBox[{
183 RowBox[{
184 RowBox[{"h2", "[",
185 RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
186 RowBox[{
187 RowBox[{"h", "[",
188 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
189 RowBox[{"h", "[",
190 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], " ", "-", " ",
191 RowBox[{
192 RowBox[{"(",
193 RowBox[{"r", "-", "rc"}], ")"}],
194 RowBox[{"s", "[",
195 RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}],
196 ";"}], "\[IndentingNewLine]",
197 RowBox[{
198 RowBox[{
199 RowBox[{"v21", "[",
200 RowBox[{"r_", ",", "method_"}], "]"}], ":=",
201 RowBox[{"Switch", "[",
202 RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", " ",
203 RowBox[{
204 RowBox[{"g", "[",
205 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], ",", " ",
206 "\[IndentingNewLine]", "TSF", ",", " ",
207 RowBox[{
208 RowBox[{"g2", "[",
209 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], ",",
210 "\[IndentingNewLine]", "SP", ",", " ",
211 RowBox[{
212 RowBox[{
213 RowBox[{"g", "[",
214 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-", " ",
215 RowBox[{
216 RowBox[{"g", "[",
217 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], ",",
218 "\[IndentingNewLine]", "GSF", ",", " ",
219 RowBox[{
220 RowBox[{
221 RowBox[{"g", "[",
222 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-", " ",
223 RowBox[{
224 RowBox[{"g", "[",
225 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}], " ", "-", " ",
226 RowBox[{
227 RowBox[{"(",
228 RowBox[{"r", "-", "rc"}], ")"}],
229 RowBox[{"(",
230 RowBox[{
231 RowBox[{
232 RowBox[{"h", "[",
233 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}], " ", "-", " ",
234 RowBox[{
235 RowBox[{"g", "[",
236 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
237 RowBox[{"rc", "^", "2"}]}]}], ")"}]}]}]}], "\[IndentingNewLine]",
238 "]"}]}], ";"}], "\[IndentingNewLine]",
239 RowBox[{
240 RowBox[{
241 RowBox[{"v22", "[",
242 RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ",
243 RowBox[{"Switch", "[",
244 RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",",
245 RowBox[{
246 RowBox[{"h", "[",
247 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
248 RowBox[{
249 RowBox[{"g", "[",
250 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], ",", " ",
251 "\[IndentingNewLine]", "TSF", ",", " ",
252 RowBox[{
253 RowBox[{"h2", "[",
254 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
255 RowBox[{
256 RowBox[{"g2", "[",
257 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], ",", " ",
258 "\[IndentingNewLine]", "SP", ",", " ",
259 RowBox[{
260 RowBox[{"h", "[",
261 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
262 RowBox[{
263 RowBox[{"g", "[",
264 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-",
265 RowBox[{"(",
266 RowBox[{
267 RowBox[{"h", "[",
268 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
269 RowBox[{
270 RowBox[{"g", "[",
271 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], ")"}]}], ",",
272 " ", "\[IndentingNewLine]", "GSF", ",", " ",
273 RowBox[{
274 RowBox[{"h", "[",
275 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
276 RowBox[{
277 RowBox[{"g", "[",
278 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-",
279 RowBox[{"(", " ",
280 RowBox[{
281 RowBox[{"h", "[",
282 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
283 RowBox[{
284 RowBox[{"g", "[",
285 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], ")"}], " ",
286 "-", " ",
287 RowBox[{
288 RowBox[{"(",
289 RowBox[{"r", "-", "rc"}], ")"}],
290 RowBox[{"(",
291 RowBox[{
292 RowBox[{"s", "[",
293 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
294 RowBox[{
295 RowBox[{"h", "[",
296 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}], " ", "+", " ",
297 RowBox[{
298 RowBox[{"g", "[",
299 RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
300 RowBox[{"rc", "^", "2"}]}]}], ")"}]}]}]}], "\[IndentingNewLine]",
301 "]"}]}], ";"}], "\[IndentingNewLine]",
302 RowBox[{
303 RowBox[{
304 RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=",
305 RowBox[{"r", " ",
306 RowBox[{"Switch", "[",
307 RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",",
308 RowBox[{
309 RowBox[{"Sin", "[", "\[Theta]", "]"}],
310 RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
311 "2", ",",
312 RowBox[{
313 RowBox[{"Sin", "[", "\[Theta]", "]"}],
314 RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
315 "3", ",",
316 RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}],
317 ";"}], "\[IndentingNewLine]",
318 RowBox[{
319 RowBox[{
320 RowBox[{"T", "[",
321 RowBox[{"\[Alpha]_", ",", "\[Beta]_", ",", " ", "method_"}], "]"}], " ", ":=",
322 " ",
323 RowBox[{
324 RowBox[{
325 RowBox[{"KroneckerDelta", "[",
326 RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*", " ",
327 RowBox[{"v21", "[",
328 RowBox[{"r", ",", " ", "method"}], "]"}]}], "+", " ",
329 RowBox[{
330 RowBox[{"R", "[", "\[Alpha]", "]"}],
331 RowBox[{"R", "[", "\[Beta]", "]"}],
332 RowBox[{"(",
333 RowBox[{
334 RowBox[{"v22", "[",
335 RowBox[{"r", ",", "method"}], "]"}], "/",
336 RowBox[{"r", "^", "2"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
337 RowBox[{
338 RowBox[{"TBare", " ", "=", " ",
339 RowBox[{"MatrixForm", "[",
340 RowBox[{"Table", "[",
341 RowBox[{
342 RowBox[{"Table", "[",
343 RowBox[{
344 RowBox[{"Integrate", "[",
345 RowBox[{
346 RowBox[{
347 RowBox[{"r", "^", "2"}],
348 RowBox[{"Integrate", "[",
349 RowBox[{
350 RowBox[{"Integrate", "[",
351 RowBox[{
352 RowBox[{
353 RowBox[{"Sin", "[", "\[Theta]", "]"}],
354 RowBox[{"T", "[",
355 RowBox[{"i", ",", "j", ",", "BARE"}], "]"}]}], ",",
356 RowBox[{"{",
357 RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
358 ",",
359 RowBox[{"{",
360 RowBox[{"\[Phi]", ",", "0", ",",
361 RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
362 RowBox[{"{",
363 RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
364 RowBox[{"{",
365 RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
366 RowBox[{"{",
367 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
368 ";"}], "\[IndentingNewLine]",
369 RowBox[{
370 RowBox[{"TSP", " ", "=", " ",
371 RowBox[{"MatrixForm", "[",
372 RowBox[{"Table", "[",
373 RowBox[{
374 RowBox[{"Table", "[",
375 RowBox[{
376 RowBox[{"Integrate", "[",
377 RowBox[{
378 RowBox[{
379 RowBox[{"r", "^", "2"}],
380 RowBox[{"Integrate", "[",
381 RowBox[{
382 RowBox[{"Integrate", "[",
383 RowBox[{
384 RowBox[{
385 RowBox[{"Sin", "[", "\[Theta]", "]"}],
386 RowBox[{"T", "[",
387 RowBox[{"i", ",", "j", ",", "SP"}], "]"}]}], ",",
388 RowBox[{"{",
389 RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
390 ",",
391 RowBox[{"{",
392 RowBox[{"\[Phi]", ",", "0", ",",
393 RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
394 RowBox[{"{",
395 RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
396 RowBox[{"{",
397 RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
398 RowBox[{"{",
399 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
400 ";"}], "\[IndentingNewLine]",
401 RowBox[{
402 RowBox[{"TGSF", " ", "=", " ",
403 RowBox[{"MatrixForm", "[",
404 RowBox[{"Table", "[",
405 RowBox[{
406 RowBox[{"Table", "[",
407 RowBox[{
408 RowBox[{"Integrate", "[",
409 RowBox[{
410 RowBox[{
411 RowBox[{"r", "^", "2"}],
412 RowBox[{"Integrate", "[",
413 RowBox[{
414 RowBox[{"Integrate", "[",
415 RowBox[{
416 RowBox[{
417 RowBox[{"Sin", "[", "\[Theta]", "]"}],
418 RowBox[{"T", "[",
419 RowBox[{"i", ",", "j", ",", "GSF"}], "]"}]}], ",",
420 RowBox[{"{",
421 RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
422 ",",
423 RowBox[{"{",
424 RowBox[{"\[Phi]", ",", "0", ",",
425 RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
426 RowBox[{"{",
427 RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
428 RowBox[{"{",
429 RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
430 RowBox[{"{",
431 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
432 ";"}], "\[IndentingNewLine]",
433 RowBox[{
434 RowBox[{"TTSF", " ", "=", " ",
435 RowBox[{"MatrixForm", "[",
436 RowBox[{"Table", "[",
437 RowBox[{
438 RowBox[{"Table", "[",
439 RowBox[{
440 RowBox[{"Integrate", "[",
441 RowBox[{
442 RowBox[{
443 RowBox[{"r", "^", "2"}],
444 RowBox[{"Integrate", "[",
445 RowBox[{
446 RowBox[{"Integrate", "[",
447 RowBox[{
448 RowBox[{
449 RowBox[{"Sin", "[", "\[Theta]", "]"}],
450 RowBox[{"T", "[",
451 RowBox[{"i", ",", "j", ",", "TSF"}], "]"}]}], ",",
452 RowBox[{"{",
453 RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
454 ",",
455 RowBox[{"{",
456 RowBox[{"\[Phi]", ",", "0", ",",
457 RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
458 RowBox[{"{",
459 RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
460 RowBox[{"{",
461 RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
462 RowBox[{"{",
463 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
464 ";"}], "\[IndentingNewLine]",
465 RowBox[{
466 RowBox[{"QBARE", " ", "=",
467 RowBox[{"3",
468 RowBox[{
469 RowBox[{"Part", "[",
470 RowBox[{"TBare", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
471 RowBox[{"(",
472 RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
473 RowBox[{
474 RowBox[{"QTSP", " ", "=", " ",
475 RowBox[{"3",
476 RowBox[{
477 RowBox[{"Part", "[",
478 RowBox[{"TSP", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
479 RowBox[{"(",
480 RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
481 RowBox[{
482 RowBox[{"QGSF", " ", "=", " ",
483 RowBox[{"3",
484 RowBox[{
485 RowBox[{"Part", "[",
486 RowBox[{"TGSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
487 RowBox[{"(",
488 RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
489 RowBox[{
490 RowBox[{"QTSF", " ", "=", " ",
491 RowBox[{"3",
492 RowBox[{
493 RowBox[{"Part", "[",
494 RowBox[{"TTSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/",
495 RowBox[{"(",
496 RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}]}], "Input",
497 CellChangeTimes->{{3.613731558054296*^9, 3.613732200209785*^9},
498 3.6137322781282454`*^9, {3.61373233169027*^9, 3.613732397354877*^9}, {
499 3.613732488348338*^9, 3.6137327713906393`*^9}, {3.613732810835518*^9,
500 3.613732821957967*^9}, {3.613732862763562*^9, 3.613732913141625*^9}, {
501 3.61373295595298*^9, 3.613733106904414*^9}, {3.613733177027177*^9,
502 3.613733199740512*^9}, {3.61373326352699*^9, 3.6137332727080936`*^9}, {
503 3.613733425731758*^9, 3.6137335719815893`*^9}, {3.613733605778666*^9,
504 3.613733648768111*^9}, {3.613733776819376*^9, 3.6137338150569887`*^9}, {
505 3.6137338699918118`*^9, 3.61373393057332*^9}, {3.613733985426153*^9,
506 3.613733985628289*^9}, {3.613734054134171*^9, 3.6137340588514023`*^9}, {
507 3.6137371511302147`*^9, 3.613737188118826*^9}, {3.613815304031159*^9,
508 3.613815306899302*^9}, {3.613815415274508*^9, 3.613815450160325*^9}, {
509 3.613815535717202*^9, 3.613815604127789*^9}, {3.613815665097225*^9,
510 3.613815673383039*^9}, {3.613816023942026*^9, 3.6138160527230864`*^9}, {
511 3.613825536655637*^9, 3.613825555245915*^9}, {3.613825660601698*^9,
512 3.6138256692518787`*^9}, {3.613825968874196*^9, 3.613826000721004*^9}, {
513 3.6138263115515537`*^9, 3.613826312408416*^9}, {3.613826919678602*^9,
514 3.613826920319913*^9}, {3.613827178556149*^9, 3.6138272197490597`*^9}, {
515 3.614594522437893*^9, 3.61459456139437*^9}}],
516
517 Cell[CellGroupData[{
518
519 Cell[BoxData[
520 RowBox[{"Expand", "[", "QBARE", "]"}]], "Input",
521 CellChangeTimes->{{3.6138263631875963`*^9, 3.613826368246337*^9},
522 3.6139024316559362`*^9}],
523
524 Cell[BoxData[
525 RowBox[{
526 RowBox[{"-",
527 FractionBox[
528 RowBox[{"2", " ",
529 SuperscriptBox["\[ExponentialE]",
530 RowBox[{
531 RowBox[{"-",
532 SuperscriptBox["rc", "2"]}], " ",
533 SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}],
534 SqrtBox["\[Pi]"]]}], "+",
535 RowBox[{"Erf", "[",
536 RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
537 CellChangeTimes->{3.613826419489653*^9, 3.613827358701158*^9,
538 3.613902416523858*^9, 3.6139057770507097`*^9, 3.613988354332242*^9,
539 3.614594677651219*^9}]
540 }, Open ]],
541
542 Cell[CellGroupData[{
543
544 Cell[BoxData[
545 RowBox[{"Expand", "[", "QTSP", "]"}]], "Input",
546 CellChangeTimes->{{3.613826369100939*^9, 3.613826374086033*^9}}],
547
548 Cell[BoxData[
549 RowBox[{
550 RowBox[{"-",
551 FractionBox[
552 RowBox[{"2", " ",
553 SuperscriptBox["\[ExponentialE]",
554 RowBox[{
555 RowBox[{"-",
556 SuperscriptBox["rc", "2"]}], " ",
557 SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}],
558 SqrtBox["\[Pi]"]]}], "-",
559 FractionBox[
560 RowBox[{"4", " ",
561 SuperscriptBox["\[ExponentialE]",
562 RowBox[{
563 RowBox[{"-",
564 SuperscriptBox["rc", "2"]}], " ",
565 SuperscriptBox["\[Kappa]", "2"]}]], " ",
566 SuperscriptBox["rc", "3"], " ",
567 SuperscriptBox["\[Kappa]", "3"]}],
568 RowBox[{"3", " ",
569 SqrtBox["\[Pi]"]}]], "+",
570 RowBox[{"Erf", "[",
571 RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
572 CellChangeTimes->{3.613988356366749*^9, 3.6145946819255857`*^9}]
573 }, Open ]],
574
575 Cell[BoxData["Q"], "Input",
576 CellChangeTimes->{{3.613907587432749*^9, 3.6139075875102386`*^9}}],
577
578 Cell[CellGroupData[{
579
580 Cell[BoxData[
581 RowBox[{"Expand", "[", "QGSF", "]"}]], "Input",
582 CellChangeTimes->{{3.6138263748828497`*^9, 3.613826377787304*^9}, {
583 3.613830939552072*^9, 3.613830971126011*^9}, {3.613905875167797*^9,
584 3.6139058760662394`*^9}}],
585
586 Cell[BoxData[
587 RowBox[{
588 RowBox[{"-",
589 FractionBox[
590 RowBox[{"2", " ",
591 SuperscriptBox["\[ExponentialE]",
592 RowBox[{
593 RowBox[{"-",
594 SuperscriptBox["rc", "2"]}], " ",
595 SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ", "\[Kappa]"}],
596 SqrtBox["\[Pi]"]]}], "-",
597 FractionBox[
598 RowBox[{"4", " ",
599 SuperscriptBox["\[ExponentialE]",
600 RowBox[{
601 RowBox[{"-",
602 SuperscriptBox["rc", "2"]}], " ",
603 SuperscriptBox["\[Kappa]", "2"]}]], " ",
604 SuperscriptBox["rc", "3"], " ",
605 SuperscriptBox["\[Kappa]", "3"]}],
606 RowBox[{"3", " ",
607 SqrtBox["\[Pi]"]}]], "-",
608 FractionBox[
609 RowBox[{"2", " ",
610 SuperscriptBox["\[ExponentialE]",
611 RowBox[{
612 RowBox[{"-",
613 SuperscriptBox["rc", "2"]}], " ",
614 SuperscriptBox["\[Kappa]", "2"]}]], " ",
615 SuperscriptBox["rc", "5"], " ",
616 SuperscriptBox["\[Kappa]", "5"]}],
617 RowBox[{"3", " ",
618 SqrtBox["\[Pi]"]}]], "+",
619 RowBox[{"Erf", "[",
620 RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]], "Output",
621 CellChangeTimes->{
622 3.613826419640369*^9, 3.6138270566935883`*^9, 3.6138273666891203`*^9, {
623 3.613830943408469*^9, 3.6138309716792517`*^9}, 3.613905783149455*^9,
624 3.613905876406279*^9, 3.613988358482745*^9, 3.614594685371746*^9}]
625 }, Open ]],
626
627 Cell[CellGroupData[{
628
629 Cell[BoxData[
630 RowBox[{"Expand", "[", "QTSF", "]"}]], "Input",
631 CellChangeTimes->{{3.61382637871056*^9, 3.61382638209606*^9}}],
632
633 Cell[BoxData["1"], "Output",
634 CellChangeTimes->{3.613826419723484*^9, 3.6138273774842777`*^9,
635 3.6139057853066607`*^9, 3.613988360907867*^9, 3.614594688339691*^9}]
636 }, Open ]],
637
638 Cell[CellGroupData[{
639
640 Cell[BoxData[
641 RowBox[{
642 RowBox[{"QGSF", "/.",
643 RowBox[{"rc", "\[Rule]", "11.5"}]}], " ", "/.",
644 RowBox[{"\[Kappa]", "\[Rule]", "0.2"}]}]], "Input",
645 CellChangeTimes->{{3.613826607232987*^9, 3.6138266182268*^9}}],
646
647 Cell[BoxData["0.8175716808656065`"], "Output",
648 CellChangeTimes->{3.613827062575906*^9, 3.613827398725811*^9,
649 3.613905789167375*^9, 3.6139075965611973`*^9, 3.613907713580661*^9,
650 3.6139883631950703`*^9}]
651 }, Open ]],
652
653 Cell[CellGroupData[{
654
655 Cell[BoxData[
656 RowBox[{"Limit", "[",
657 RowBox[{"QBARE", ",",
658 RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
659 CellChangeTimes->{{3.6145954324899607`*^9, 3.614595438979292*^9}}],
660
661 Cell[BoxData["0"], "Output",
662 CellChangeTimes->{3.614595439312615*^9}]
663 }, Open ]],
664
665 Cell[CellGroupData[{
666
667 Cell[BoxData[
668 RowBox[{"Limit", "[",
669 RowBox[{"QTSP", ",",
670 RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
671 CellChangeTimes->{{3.614595446545167*^9, 3.61459545529193*^9}}],
672
673 Cell[BoxData["0"], "Output",
674 CellChangeTimes->{3.614595456000258*^9}]
675 }, Open ]],
676
677 Cell[CellGroupData[{
678
679 Cell[BoxData[
680 RowBox[{"Limit", "[",
681 RowBox[{"QGSF", ",",
682 RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
683 CellChangeTimes->{{3.614595464309478*^9, 3.614595469692217*^9}}],
684
685 Cell[BoxData["1"], "Output",
686 CellChangeTimes->{3.6145954699957027`*^9}]
687 }, Open ]],
688
689 Cell[CellGroupData[{
690
691 Cell[BoxData[
692 RowBox[{"Limit", "[",
693 RowBox[{"QTSF", ",",
694 RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
695 CellChangeTimes->{{3.614595479813127*^9, 3.614595487321061*^9}}],
696
697 Cell[BoxData["1"], "Output",
698 CellChangeTimes->{3.614595487675913*^9}]
699 }, Open ]],
700
701 Cell[CellGroupData[{
702
703 Cell[BoxData[
704 RowBox[{"Expand", "[",
705 RowBox[{"h", "[",
706 RowBox[{"r", ",", "\[Kappa]"}], "]"}], "]"}]], "Input",
707 CellChangeTimes->{{3.614595726057316*^9, 3.6145957369270573`*^9}}],
708
709 Cell[BoxData[
710 RowBox[{"h", "[",
711 RowBox[{"r", ",", "\[Kappa]"}], "]"}]], "Output",
712 CellChangeTimes->{{3.6145957290402822`*^9, 3.614595737188229*^9},
713 3.614595821337948*^9}]
714 }, Open ]],
715
716 Cell[CellGroupData[{
717
718 Cell[BoxData[
719 RowBox[{"s", "[",
720 RowBox[{"r", ",", "\[Kappa]"}], "]"}]], "Input",
721 CellChangeTimes->{{3.6145957525304213`*^9, 3.614595789766198*^9}}],
722
723 Cell[BoxData[
724 RowBox[{"s", "[",
725 RowBox[{"r", ",", "\[Kappa]"}], "]"}]], "Output",
726 CellChangeTimes->{3.614595790131374*^9}]
727 }, Open ]],
728
729 Cell[CellGroupData[{
730
731 Cell[BoxData[
732 RowBox[{"Expand", "[",
733 RowBox[{
734 RowBox[{"s", "[",
735 RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",",
736 RowBox[{"{",
737 RowBox[{"r", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input",
738 CellChangeTimes->{{3.614595791670055*^9, 3.614595812670266*^9}}],
739
740 Cell[BoxData[
741 RowBox[{"s", "[",
742 RowBox[{"r", ",", "\[Kappa]"}], "]"}]], "Output",
743 CellChangeTimes->{{3.614595794966937*^9, 3.614595812909087*^9}}]
744 }, Open ]]
745 },
746 WindowSize->{884, 1337},
747 WindowMargins->{{Automatic, 355}, {Automatic, 0}},
748 FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \
749 2014)",
750 StyleDefinitions->"Default.nb"
751 ]
752 (* End of Notebook Content *)
753
754 (* Internal cache information *)
755 (*CellTagsOutline
756 CellTagsIndex->{}
757 *)
758 (*CellTagsIndex
759 CellTagsIndex->{}
760 *)
761 (*NotebookFileOutline
762 Notebook[{
763 Cell[400, 13, 17508, 501, 845, "Input"],
764 Cell[CellGroupData[{
765 Cell[17933, 518, 159, 3, 28, "Input"],
766 Cell[18095, 523, 541, 15, 58, "Output"]
767 }, Open ]],
768 Cell[CellGroupData[{
769 Cell[18673, 543, 128, 2, 28, "Input"],
770 Cell[18804, 547, 777, 24, 58, "Output"]
771 }, Open ]],
772 Cell[19596, 574, 95, 1, 28, "Input"],
773 Cell[CellGroupData[{
774 Cell[19716, 579, 230, 4, 28, "Input"],
775 Cell[19949, 585, 1278, 38, 58, "Output"]
776 }, Open ]],
777 Cell[CellGroupData[{
778 Cell[21264, 628, 126, 2, 28, "Input"],
779 Cell[21393, 632, 165, 2, 28, "Output"]
780 }, Open ]],
781 Cell[CellGroupData[{
782 Cell[21595, 639, 218, 5, 28, "Input"],
783 Cell[21816, 646, 208, 3, 28, "Output"]
784 }, Open ]],
785 Cell[CellGroupData[{
786 Cell[22061, 654, 190, 4, 28, "Input"],
787 Cell[22254, 660, 70, 1, 28, "Output"]
788 }, Open ]],
789 Cell[CellGroupData[{
790 Cell[22361, 666, 186, 4, 28, "Input"],
791 Cell[22550, 672, 70, 1, 28, "Output"]
792 }, Open ]],
793 Cell[CellGroupData[{
794 Cell[22657, 678, 187, 4, 28, "Input"],
795 Cell[22847, 684, 72, 1, 28, "Output"]
796 }, Open ]],
797 Cell[CellGroupData[{
798 Cell[22956, 690, 187, 4, 28, "Input"],
799 Cell[23146, 696, 70, 1, 28, "Output"]
800 }, Open ]],
801 Cell[CellGroupData[{
802 Cell[23253, 702, 186, 4, 28, "Input"],
803 Cell[23442, 708, 179, 4, 28, "Output"]
804 }, Open ]],
805 Cell[CellGroupData[{
806 Cell[23658, 717, 152, 3, 28, "Input"],
807 Cell[23813, 722, 127, 3, 28, "Output"]
808 }, Open ]],
809 Cell[CellGroupData[{
810 Cell[23977, 730, 269, 7, 28, "Input"],
811 Cell[24249, 739, 151, 3, 28, "Output"]
812 }, Open ]]
813 }
814 ]
815 *)
816