| 1 |
(*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 |
|
