/
heteronuclear_levels.nb
2336 lines (2289 loc) · 116 KB
/
heteronuclear_levels.nb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 7.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 145, 7]
NotebookDataLength[ 118866, 2327]
NotebookOptionsPosition[ 116930, 2263]
NotebookOutlinePosition[ 117352, 2279]
CellTagsIndexPosition[ 117309, 2276]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["\<\
Single-Electron Energy Levels of Heteronuclear Diatomic Molecules\
\>", "Title"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "diatomic_base.m"}], ";"}]}], "Input"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"relative", " ", "nuclear", " ", "charge", " ", "difference"}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"\[CapitalDelta]Zrelval", "=",
FractionBox[
RowBox[{"5", "-", "4"}],
RowBox[{
RowBox[{"(",
RowBox[{"5", "+", "4"}], ")"}], "/", "2"}]]}]}]], "Input"],
Cell[BoxData[
FractionBox["2", "9"]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[TextData[{
"Calculate single-electron energies in dependence of ",
Cell[BoxData[
FormBox[
RowBox[{"Z", " ", "R"}], TraditionalForm]]]
}], "Section"],
Cell[BoxData[
RowBox[{
RowBox[{"CalculatePlist", "[",
RowBox[{
RowBox[{"{",
RowBox[{"l_", ",", "m_"}], "}"}], ",", "pstart_", ",", "ZR_List", ",",
"\[CapitalDelta]Zrel_"}], "]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{"p0", "=", "pstart"}], "}"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{"p0", "=",
RowBox[{"CalculateP", "[",
RowBox[{
RowBox[{"{",
RowBox[{"l", ",", "m"}], "}"}], ",", "p0", ",",
RowBox[{
RowBox[{"{",
RowBox[{"2", ",", "\[CapitalDelta]Zrel"}], "}"}], "#"}]}],
"]"}]}]}], "}"}], "&"}], "/@", "ZR"}]}], "]"}]}]], "Input"],
Cell[BoxData[
RowBox[{
RowBox[{"ZRrange", "=",
RowBox[{"Prepend", "[",
RowBox[{
RowBox[{"Exp", "[",
RowBox[{"Range", "[",
RowBox[{
RowBox[{"Log", "[", "0.05", "]"}], ",",
RowBox[{"Log", "[", "32", "]"}], ",", "0.1"}], "]"}], "]"}], ",",
"0"}], "]"}]}], ";"}]], "Input"],
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
"n", ",", "l", ",", "m", ",", "pFitSlopeOffset", ",", "numRadCoeffs", ",",
"color"}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"paraInfo", "=",
RowBox[{"ParaInfo", "[", "\[CapitalDelta]Zrelval", "]"}]}],
";"}]}]], "Input"],
Cell[BoxData[
RowBox[{
RowBox[{"colors", "=",
RowBox[{"{",
RowBox[{"(*",
RowBox[{"1", "s\[Sigma]"}], "*)"}],
RowBox[{"Blue", ",",
RowBox[{"(*",
RowBox[{"1", "p\[Sigma]"}], "*)"}], "Purple", ",",
RowBox[{"(*",
RowBox[{"2", "s\[Sigma]"}], "*)"}], "Green", ",",
RowBox[{"(*",
RowBox[{"2", "p\[Sigma]"}], "*)"}], "Red", ",",
RowBox[{"(*",
RowBox[{"1", "d\[Sigma]"}], "*)"}], "Black", ",",
RowBox[{"(*",
RowBox[{"1", "p\[Pi]"}], "*)"}], "Brown", ",",
RowBox[{"(*",
RowBox[{"1", "p",
RowBox[{"(",
RowBox[{"-", "\[Pi]"}], ")"}]}], "*)"}], "Brown", ",",
RowBox[{"(*",
RowBox[{"1", "d\[Pi]"}], "*)"}], "Orange", ",",
RowBox[{"(*",
RowBox[{"1", "d",
RowBox[{"(",
RowBox[{"-", "\[Pi]"}], ")"}]}], "*)"}], "Orange", ",",
RowBox[{"(*",
RowBox[{"1", "f\[Sigma]"}], "*)"}], "Magenta", ",",
RowBox[{"(*",
RowBox[{"2", "p\[Pi]"}], "*)"}],
RowBox[{"Darker", "[", "Yellow", "]"}], ",",
RowBox[{"(*",
RowBox[{"2", "p",
RowBox[{"(",
RowBox[{"-", "\[Pi]"}], ")"}]}], "*)"}],
RowBox[{"Darker", "[", "Yellow", "]"}], ",",
RowBox[{"(*",
RowBox[{"3", "s\[Sigma]"}], "*)"}], "Gray"}], "}"}]}], ";"}]], "Input"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"QuantNumbersToString", "[",
RowBox[{"#", "\[LeftDoubleBracket]",
RowBox[{"1", ";;", "3"}], "\[RightDoubleBracket]"}], "]"}], ",",
RowBox[{"#", "\[LeftDoubleBracket]",
RowBox[{"4", ",", "1"}], "\[RightDoubleBracket]"}]}], "}"}], "&"}], "/@",
"paraInfo"}]], "Input"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"\<\"1s\[Sigma]\"\>", ",", "0.5555555555555556`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1p\[Sigma]\"\>", ",", "0.4444444444444444`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"2s\[Sigma]\"\>", ",", "0.2777777777777778`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"2p\[Sigma]\"\>", ",", "0.2222222222222222`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1d\[Sigma]\"\>", ",", "0.2777777777777778`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1p(+\[Pi])\"\>", ",", "0.2777777777777778`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1p(-\[Pi])\"\>", ",", "0.2777777777777778`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1d(+\[Pi])\"\>", ",", "0.2222222222222222`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1d(-\[Pi])\"\>", ",", "0.2222222222222222`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"1f\[Sigma]\"\>", ",", "0.2222222222222222`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"2p(+\[Pi])\"\>", ",", "0.18518518518518517`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"2p(-\[Pi])\"\>", ",", "0.18518518518518517`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\<\"3s\[Sigma]\"\>", ",", "0.18518518518518517`"}], "}"}]}],
"}"}]], "Output"]
}, Open ]],
Cell[BoxData[
RowBox[{"(*", " ",
RowBox[{"2", "p",
RowBox[{"(",
RowBox[{"\[PlusMinus]", "\[Pi]"}],
SubscriptBox[")", "u"]}], " ", "3",
SubscriptBox["s\[Sigma]", "g"], " ", "3",
SubscriptBox["p\[Sigma]", "u"], " ", "3", "p",
RowBox[{"(",
RowBox[{"\[PlusMinus]", "\[Pi]"}],
SubscriptBox[")", "u"]}]}], " ", "*)"}]], "Input"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "\[LessEqual]",
RowBox[{"Length", "[", "paraInfo", "]"}]}], ",",
RowBox[{"i", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"curInfo", "=",
RowBox[{
"paraInfo", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}]}],
";", "\[IndentingNewLine]",
RowBox[{
SubscriptBox["pPlot", "i"], "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"CalculatePlist", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{
"curInfo", "\[LeftDoubleBracket]", "2",
"\[RightDoubleBracket]"}], ",",
RowBox[{"Abs", "[",
RowBox[{
"curInfo", "\[LeftDoubleBracket]", "3",
"\[RightDoubleBracket]"}], "]"}]}], "}"}], ",",
RowBox[{
RowBox[{
RowBox[{"curInfo", "\[LeftDoubleBracket]",
RowBox[{"4", ",", "1"}], "\[RightDoubleBracket]"}],
RowBox[{"Max", "[", "ZRrange", "]"}]}], "+",
RowBox[{"curInfo", "\[LeftDoubleBracket]",
RowBox[{"4", ",", "2"}], "\[RightDoubleBracket]"}]}], ",",
RowBox[{"Reverse", "[", "ZRrange", "]"}], ",",
"\[CapitalDelta]Zrelval"}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{
RowBox[{
RowBox[{"curInfo", "\[LeftDoubleBracket]",
RowBox[{"4", ",", "1"}], "\[RightDoubleBracket]"}], "#"}],
"+",
RowBox[{"curInfo", "\[LeftDoubleBracket]",
RowBox[{"4", ",", "2"}], "\[RightDoubleBracket]"}]}]}], "}"}],
"&"}], "/@", "ZRrange"}]}], "}"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<Z R [a.u.]\>\"", ",",
RowBox[{"\"\<p(\>\"", "<>",
RowBox[{"QuantNumbersToString", "[",
RowBox[{"curInfo", "\[LeftDoubleBracket]",
RowBox[{"1", ";;", "3"}], "\[RightDoubleBracket]"}], "]"}],
"<>", "\"\<)\>\""}]}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{
"colors", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}],
",",
RowBox[{"{",
RowBox[{"Gray", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"Print", "[",
SubscriptBox["pPlot", "i"], "]"}], ";"}]}], "]"}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwtk4k3FAgcgAc5tkSOJlEYU8ZEGee46vdzbDKpRImkYh2hnqTaduMpVwel
rHrE5nwy7iNnMfMaTZJj5MgxJJPJsRVtE+Nq67393vve9xd8JN9QV39JAoHg
9cOfDUuLVZLptEfOjtBIUaA5Vsvai9JP2mHMhZUdHJ4ppqraMovkbXFKaf+f
C+4mGJIX/fZ+DiD52rT25Xkj5FBo/CDPnag6893UuYqGjNuGNnFW1qjTRI41
XmOICWPJw4N2lsilPj10YsgA/Qflde6F0rEHF4Pmu7dh2h9JA8pNZqg2UJKs
9EUPsT38kkDfFP35xsyvxhT0jKv+i1NjjE2d0cWbU7YiXdf4YdoxI6xj1JDo
Klvwjrnl2DkNGuLJD3XEUh2MCiQHSAi24wdKJa3Bl4SkU915+g36KLW8HO5g
qI2vu9+u2VNPxY3eyBeRNdHQ2qIqfJiCUYQ6crzyJnxI6SdKbNHFB/612khR
RxctA+vzN7egDdE3wueEGuYrdKtFKZCxN9Lc2auOiMmqKYyVEhIKrE6ndhiu
x0lF9rwoQBs7zk243OCoYFDy5sXNoZp4xs/tRUyYMlp5VEXo2WxCGz137ghd
CVWPNoeMGqmjh5sVY0x1HW62NfI4vE8N4UW4k3itAkaVmtDaE4jolvilO3qj
PJ53KHFlT6mik9mlf9qdVmNGFmmqN1AFie4OgqsZcnhDTSNhlKCMbvRK4hlF
WZyjaigNVa1Dbrj6dFCWNJ5yJ4e8WlbAR6zpZvt9q/A6x7U+ni2PvzfNBRao
SOGTAHm90uLV6O5t27NjVgJ7jCuj7FlySNxUvXxWQMCcK+/Go0QyKMnLrHHn
fQdFtyPDnc7SOC7lZ1C4sAwrFj1Jnc+kUH1gItjdbgmE3/OPOR6WRH4RY+hk
3gK4uU7k60lKoHKDkYuzphg+7XzQuD7pO6y0i7nSpXPQl3gv5qDjMoRrySzz
Dn6DGIUd6U+oi3DWvX0XXU4EFq9Kmh1NxKDNT9UTd/4LNTcpJB2fOaizmt8l
VfAFfltK0yoqF0FzRZjIL2kWEgLOVIdrfYX3nhOFTnEzkLgywLcs+gIBtJPe
QqXP0MehHho/OAu6pyImg00+giVF98cRM+AcyvQXnpmG7JlyixTeR/j41tLX
nD0JsKLnosKcBunnhXudqRPQeF3zcwhzEjb05S1SC4RgeF1Z9LTjA4gynb+N
WI9D1GJ+7FGiEGjFL43iBQIw7vbb1hr5HsRWtP0OD8fgjl0lgyUpADaLsG13
yDv4mN3TtSH3HdQd7+CwnEZh1sSCNHBsFELbr9b3hY1AUKt2yNFjIxAx+TxX
6SwfVmSCpdjb+JBRFuPYkDgIUS0DX+tJg2CaJ8Pw4/bDtWHLTM9d/WB2nJtc
rPEGMloomcLIPhgvddK0TOiFyL7OscHhHhjVkFqsU+mBYA9d2bde3VA2sj7g
fdVrGFL0neeKuuDJ2VcrEcFd0NY4xLvN5EGRyFduvzkPyo9kC8qmOmCw8JG/
p0EHmPdW5raVtsGskN4yo9QGW1kBC9YZrcAOuyQ202iFOLlr1yXKW2Apnaqr
vLsFPLO+UT0+cOENo8Eh9g4XNDXXzpfAcyjlvT+cPd8Mh04338uq4cDUFbu4
lAgOWM/4ssVOz0DxAlsvU/0ZpHuUaRXPsaH0nn2mWhcbxLe252RdZoF3dL8W
34YFC+PSbg5mjSB7N0XNQbIRbvkMuC2pPQF9HyEt9XMDFNQSdjeR6+G1foEF
TbIevD4NytW51oIL+0FMAr0WvC6HM/fmVsOebJ/4A7eqgXrx0N9rVB7D/sZf
vXZLPAayJ9G+NrMSQgnHzbzvVsKlPbm5dMcKyFll3pJjXQHX5K+GZsiWw3mX
QdeEpTJ42HE/0LipBLw+2/huLSuBqgMVcdKpRZC2K7hBmFAE6Wk/YYLLS5eL
/UlMIPzPdlPbX55K1rP+AwhAp8Y=
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk2k0FAoYhoeiudmyOxFNhEk1jGIm6fsmMrYSkrG0kCxJjlS3bgsqcRNK
pUT2220sEdeWRGm7EsmSNbPYpWxDY736eZ9z3vOe5/d7XpJnoMNRUQKB4LaU
X/1/OJCY8As20GP4tBg+BwrtnoaL3c8Gox8Ru5kCDiTXxftQX+TCl1srqzyI
XIiQDAtMWpEPKpRpRaEaF85aZmQYM5/CHyLxwg2GXNB0UTIrSSmASvZjL2Ub
LpDP7HsoIf8PdFPP13G8uOB2Pphtk1EE229HHcgJXfLv7cRShxJQ2yKzJj2Z
C49LCBYvNMtg5kGG0cgLLkR7tDnOqZRDkP1F/edcLsz0ijmab62AqrQT+owV
PBBGb0pPPV8JiqXa5gUUHiSy8jRypqtAaYdsnpUbD0xGPauEVq/gYJrFDovr
PNh3/PXd1OJqeDAuljH8nAfq6lI/c+ENTBVta7o8wQOX1Ckyq/8tXNJK/N11
Mx/CiRGRIvnvwdduyuFdAB/WV3rPmCTVANvD6a1WAR+Mmgsyap/Uwvh4il3z
LB/yndP4eUN1YC0fpmJm0wO1FR2fYtifoPm8SywntQc6ZDx/vhU0gI2pqiR5
vgeOsbRXdLs1Qk/AXFja4V642FLPa+9qAvFWArR96IWk9zopfRdbwDLxxLNj
0AcRXfQUlx2tsJfiHDxX3gch79smy0jt4O55msPY2Q8L4seWVW3ohNDespNz
n/vBr2atv6v7V+jhzWTPnRiAMUMaqc2dA51qLBtN5UEYSWtqUM7gwpTL8lue
tYNwc2eBdaUoH17FWyRWxg4BtdFrQ83FHsiqYX919ByGkNlHV12V+uCI0bny
Pbu+ASVSTvC8rh9W+zZQp+kjUBGp/sOfPQi2aeUVuYzvAAu6e+XZw+C+b8by
+4EfkDaaT7vzaQS8tQ02ki1Hga6jnXBVdhTiIo5lWY2MQks1eV+v/RgMGipf
tUgZgxsLbZ307HHoZLuvX3toHKK8A4qCNSYhnb9sYmjTBByZS9DIzhfADW5W
6yfJSSi+rkNa5zENk+FCgdz8JNA+5L5mGgrhlqp5gARhCq5Ib04sJ8/CKqLR
foryNLTcuHvFnjkPboqOX1aZ/4Tvpg8qFGMXIS3fyc8jUgiODgOPdEVFED/m
MCd7Z6Bv8ZE700kUndQeyr48OAcLtKbY+lfL8I1bfUrf5DzIODp31duKoUDc
Q+5J1iKkh3J7QwTieFv6hxTzJQGbqAUhZpVEzKONWdyrFMFyb0ndJzkr0fYd
KBU3iGJktUPZtSpJtFUWeG8WLkPf/Zr+H+alMVZev49tLIbTZFXZjsJVSJFz
9p2IEcc/VVSjOAQ5NFNVvyQkEDEplTTU7COPrGGiaUPsb3jKPNehakgBkyn+
UGYqgSFPDPU/Rilh9HE6R01KCtcwDFhOu1VwQM+x5B5BBhVcX/tzDFajrvCv
3AvRq3Abq/CC7nY1PLNbuiQzTxb94tbMrglUR2d7q/mjA3I4KFP1U+C9FvtT
jihcoStgnMId64VcEl57SX3MzVTER9KNKiHSmtg8IU6L3qiMezU2mpy6roXL
n3YQn31WwWSdViURLW10C9vl8D5hNVJMaIXBXTpYzcjPYl5Vw8+N3RKWZWR8
ip0lM3fUkeTbmKn3TA8rTp6+E+CzFkN8NL1F+Juw7F37YswDEt40ovNOqupj
W81H5aTedWisTU1OcDfAvFLj2iBrLXQJL7pdXUzF0EFnqtyH9Uu7Bp/l621B
L7qX399+OphwLrZN7sVWpO9pLb2mS8aj7ZLr7gYaI61l7JvhSj2M4sV1te+k
46UAgykv2U1oHUPZHr7NBI2OqjSIMShYraPf6ediilu6uAfU6/TRP/Nyd3w6
IMvt2SFjeSreV2CwsyUZ+NbO0jo9yBCLVpgJEg/vxH/jZSedh7dg0NIbxOvN
sHHRNPzsZSP8D3T7qGc=
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(1s\[Sigma])\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.09769797540612109,
22.509651496540595`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"FindRoot", "::", "lstol"}], "MessageName"], ":",
" ", "\<\"The line search decreased the step size to within tolerance \
specified by AccuracyGoal and PrecisionGoal but was unable to find a \
sufficient decrease in the merit function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindRoot/lstol\\\", ButtonNote -> \
\\\"FindRoot::lstol\\\"]\\)\"\>"}]], "Message", "MSG"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0.5, 0, 0.5], LineBox[CompressedData["
1:eJwtk3k0FAgcxwc5apAriRwTGZMVpsNI9ftJbyeqrUS5OjfEaG1ppUOTNJFJ
2k4TYUYz2+TMREmOiC0zriZknDkSlUQTRdPWe/t57/s+/37/+ZD2hnsGKhMI
BP8f++mDnDO6ag1uuK+nd6SVRsVCdTd5yu41yJY/Xjo65oDJBq7CLE1X5N/x
Mjc5Z4+MW6e7r/EA7R4/eze53Q6ryA4dIb6rsH3FR8GWTbboccF+JWuFC3Z3
dDGdwijI7r3UKVvjjKdFeY1jWWQMlGkuuBruhBeqdKgKojVyjia16ZUtwxAt
SadyghViXURUn+1SLOjw9mq2tERfVuHlqiIqDueLXQ+1ktDJmprGCXDEoxu8
0qZ4FnhxuXPvIRMH3BFcphk20xyZwZZBSn12OMyqS9L8PB9J+6W3bB/aojhl
YG2Xhgk+l3YT1xVTUL5QPdwM5qG9C00U0UnGV3WyDauvzcU08ktDJStrlPwx
ze4mGuJm819cDv/41fOpo6jxpgEKtKVGTG1LjDRRkGZ66OMlgyseihwSVp9N
52dq6+HQ7IpJeZAFru8XHW4d0sGQS6ZTpuFmSMzVknaKZuMKH9EJm5XzMYRv
f+7YQy008HvC6HE0xgMlVqG2/UQ0dXX08d5ohDpUvaBEu1nIzF3iUMc2xE0N
d1RplzXw8Nocz4phA9xBFbc26KpjagZpuDlYH0W6mk1WAlU8Z2TC7iHo4fiY
uYXPxhk4QTHRbRfpYI8WdZmKqgru32bJEH/TxhVUSuI2sRLGV3kWn63QxGj1
k1yrdAKWBGna5GbPwtSoU+OZ3xXwglrAdCvXwMtLc4/slE0D79SrAaZcDSta
ogjJ9V9h9tbtnQ0bVDEyLJHF7ZgEBe1FUkOlCsYqJ3tJZkzA6++CALq3Mq7O
nIzYTJfDVs83AhtlJZR9Tq12vTUOI6tulM5J+g7Di3flsk3GoOX81dgt9G8Q
I3fwHckehVjtxSkllCko4fmxJvJGgCbOeUJf8gXEkt8SKLR3UJRAJi3YMwEG
K1Vq60eH4PdpjnlWvhz07uqkdbcOAjvoQGGE+Sdw4LJnfR0egPOKtg7nrDEQ
Pjmdd4jcDy1VFK+BLR/Bz5/PZyT0gjPZ+kcRo5CnOOs8ovMKuKP5tCuN76Gv
KZqk8Wc3gMJms77wLfgxFQLOlw4ojTf7wBAOge3wW2fTYhnYx+vJH9UPwjzx
FCNG+BKYU4IzfoavgSCIdVhQ0wJU6b5FtdH9cMw9LkE8qxkurinwKFfug8rG
uA1FEVJ4z33RNDfzFZx220mK+dYEH5fQSG0BPRCZtotH5jdCSK0Fwy+gCwb9
xk98bKgHhVqoSsWiDriUFBlvHCcB5tO2T8UkGRgjXfckoxbiOp3TfVe/BHXO
/DrXo08h9Sk5/XV0CyRfT/b6kF0D0S0NvbLOFzDDU/osWKUaQn2s1bv9pbBT
EmTbdrQK2mfvnayRN8ExO+7CQq1KkJS2N14QNsKaGQIH4aMKyN/O7csbrgee
6ZFk+5nlsLy5IFOSKwFjtfyAhf8+goXlQV9dUmvh5INz0Vl5D4GlERevlP8U
4vSJTFb5A/DN+EzxGawBlzfbtUrGi8DMTGsyB6ohxjUs8G/3QvAKe3I1o6gK
ej0l51llInAZ3Vvxxb0SIDSBRt9cACk+eebZExUQVXrQ/d50PnxJtONlHC+H
4zPpScLYPPg6oLp17bJSaH638XquSw4k7mnbOm1UAhEab9kmZllw+z7h1zLL
YnA00nVsJgvBf0Sm8cDzPlxMevthyvsf8D8eIVyfWQifq3mGg3w+UCK9bhL1
74HR4+dnLebcAktfQ7f76QXwXWBUs5vLg6h1mZlO9LuQyhgLZK/jQpxmTHiq
ej5UPr9bnKWRAWn114KpZTkQdE+L31d+E0Sb7rJUk7NAYhP41y5OCqRwfiKE
mqnb0pwbHCD8j5TY3EWv6Sn/D1PMssk=
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk/k7FAgch4cKGyW3HMPknOTuINrv1xFFSigm2g2ZcWyp2NSDZ0KWEkkp
Imcp64yQxDhjO1CONaPJHI5SiTK52fpx3+f5PO/z+QNekk+wi58wgUDw/LGf
/j8cyEj/SSF4i1cfE6/mQOWBh7Fr0oqgTatnrXkrB7I6b9JMGkpAZvAkNaqH
A3ESUcGZouVw/uts2DyfA+f25OfvsH8IqbR383UCDmhQ5G1qsiuASnoSMyTG
BfJZtzviMo9g1NPw3gVVLniGhxQ65leB0ds8hXrTH3+CJfbYpQbKyRLDjxy5
8KCGYNegUQsL4Sopl/24kOjNdF1UrIMm1hLFN4oL8yNrXG231QNB+SuFls2F
uUT9vJxwBvw9R+xuYnAhw6NMrXimEaTSWB65PC5YTPo0zu1thodKb4zNxHjg
9kdrak51C4TrBkU/NeIBkbhutgTaIOVaR6G/Fw8oOd/JHmPPIC2GlxKewINY
sbh4ofIOqD2jJ2newAMtBnXeIvM5qO+3NRuf5sH2vor8l6UvwbN5H5NtyIdy
91x+2Xgn5FnTOP7BfHhZP9idVNgNUr6yzNFKPgxK+sw+E7yGvOdyp54s8SHQ
Q1t0yLMHnBtuCu13GobI/i4ei90LZjI0E9n8Ycjs0MkejewHkb1aqfErwxDH
Ns+m/DoAIYzdXka+I0DvYE7XkljAj5rZmNE5AssigasaN78FAnEhydd6FAKe
qwcd8XoHztGpl4wYozBlakZienFAXoHTV7V7DD7n9r5WyOfCrH03VvWPQbJ1
hQNDmA+dRUFlK2feg0nP8c3PI4chXqzO8azyB6AvFFw8Ij8KKm36/4i//gCG
8dKCp51jIKtu1dx2Yxzq44lfggo/wKkxH+MG2keAZV1nmcKPwC80ENVy/AS5
k+VmN7o/Q3/+fo4VfgZzHe30i1KTcADC2H/umYD+FrLbyMEpYKlp5ezy+wJX
lplvzYu+QrlbZM3GfZOQQD1RFaI2De+jA5rjpybBdzFdrahcAP1KX+xY+VNQ
fVmHtMl7BppXG+hJHv8KZi9KWu1N5+B+WYC7pek3iFlvkFFHXoAQkVnbK9LT
0H8lNeag/RJEhFkX4yoBTOy6XS93dQWcEm3dHol8B1eX9wW6wkJI/VdD5Chx
BkZXCrzsDwljGL/8/ivHWVg2673a1bwKd8gz6Y+T50DS1Z3dtW8NkrOdHkxN
zEPeBe4IXSCCxOHhstP+i9BrUkG3YYihw3TUaY+VJaijSuiWFq/FoV/TNSar
ViC+xaX2r0YJtCw+6uLRTkD/wxpBL5bW44D1t4ukZ0I4Q1aWGqzcgAyrodrm
AWG8pKicwCFIY5ecYlo7YTVm5pDG+2gyyDm08xYL1mCobYlL47gsOlC++LXf
EkF6qanRqwR5PD2pfyFyrRiqWhl7HHJSRJfDknPSt39B2SOtQRxjJWQ2WUa0
7xHHnR6VEbqWKhgv55AgunEdBqSoLqgGE1Ez+sFC8DpJ/CDZOCugquMnugrb
/8YGTJG94bBcQsIsmuCTUo0UFqzvUaSv18AwcvbcqUlpdFbbYhF6WRMzEyNO
GFjLYpbOgLyQpjZyXwUOmpfKoaGFWWUIWwcztjbWZW9XwDc9Q+J7asn4jcr/
PZWtiCT/nrt6T/SwY7xOPeKeEtJpGlQhvj4mwrWXw8kqmLzdnHdG2QiTbPdS
rucRcYe2SVa6lzEyrm851XRSHSmxVddbqk3QhCQjCM0jIb4KOcfX24oq3jFx
ahObMP38VaZ0wzacu2s3n+WqiX4siU2pwTsw5Z2O+/deLUzgpbBZ1uaY5j7d
6hyqgw5JhpaxOy3wtxq7+6u3krFFx+htAGUXqqVmLt6U18Ogu9FDN/MAtVOS
VxqI+pgma1VYJGGFfRUhGwIPGGKVqI0g45g15p4L/T7UZ4Snf9Qg0mWDkqYv
7lxXNsH/AFjtr2I=
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(1p\[Sigma])\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.5242041923705706,
18.2220448783777}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"FindRoot", "::", "lstol"}], "MessageName"], ":",
" ", "\<\"The line search decreased the step size to within tolerance \
specified by AccuracyGoal and PrecisionGoal but was unable to find a \
sufficient decrease in the merit function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindRoot/lstol\\\", ButtonNote -> \
\\\"FindRoot::lstol\\\"]\\)\"\>"}]], "Message", "MSG"],
Cell[CellGroupData[{
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 1, 0], LineBox[CompressedData["
1:eJwtk4k3FAgcgJkKhVFuYZgcQzWJDoT9/aSSylaohqHWbKjUztaobCVtaESh
Gr0mHa5qJw1iGcZZo5I7gxpy507OCSltvbffe9/7/oKPzGC6+RFkZGToP/zZ
49zwJXK1TsjUVG0/cpKCOfJO0oTfNmKxGvfV7RYTvKXuyEtTckS69nGG0gFj
DEy92H4zGXDJ91b5e9+WoYiy+v1hTwfc/3o6eE8OGbfFWNhHbLDD7vB6ok+4
IUZ3XW9t3miLkdkeDoR6Evo1Ky2LZ1rjtUciEl2kh9y/YiWqxetQGO9k+1ay
FLGaFdy9Yi2aXWiMLFHVQc+InBuiXCssG1iXFxuohdamVve43pYYl5rKGe/U
wLj1tl0ndFcji6bzqpiljqEBRv6y3VRUeMuunSCpIfmQOHWFcAVu/4UNgs4l
WC9uV9yab47je5mWbsLFaGFnk81qpeBGjYxhn19V8B7lnaassSlGnhf3812V
cZfBSrugKGMUD6tYdTIV8SFRrB1KNEI9o2jemayFeF2ds22OT0ZBVL9klZoC
DqiUTkv9DZEuJ9Jqj5XDw9f1Z/WZJOyPiZ65YboAN9Cyz5nZ6yHRiUM73TgP
1b3KAjssl2LpiE9IzW0C6jta0va4amOPYdsZaZAshqavWV0drYkOvIPF3gwZ
DNrEdysdVMcBBu+qqmgO7iSSBxsD1JDKiZ6ldn2Fy9q60R0yqjjfzuHZaYNZ
mDLXXdKSvRj7IuYWxQfPwKG9RoGV34gYNp2i5z80BZEit/xLpUo4TV1aMP/U
ZyjwVzJLf7IIU4OXtpZoSqHBKivUqUQBifsJbzSqJiD5QmdPqFQOYy2VBvic
cVBx39dau2MBelaFHXZjjsGcTUNs7fN5WC2R6g95jULv94feznsISFXuChp7
8gnc3fofmhFkce9wBnHZ84/wyeF2kUbsd/hsonNJMDEITVfiw3Y7f4PQnDg7
QxyAMOKqhALzWQhZ4RfQ86APbCr5Zc5rZsAj2/7dI6NeyI2ikJf5TsHjtsu7
9fM+wO9fuQZpmVJ4vjIygMjohmj/Yzksg0mIKfcp8jDsgitzkve2aePAG/gY
KxnrgCaRuUfP7jHIrafNmKW1gy3F9McRo7CPQh9Rj2+FpNFMG07dMAi2t1Vz
77cAzJntUuMNwYvNb3PHyiVQFEkaCeQNQMj5VWVHVd6BRaSqtLCmD9Im8+uY
x5sgdPZhuJdmL4Sp+/myhhvASnxweUXIB2gz2MQjXBRD3MasbSWEbngAZ1OF
1HoYTmp4o5XSCTqNRVtcRupgbI0NWeLdAeeab56URNXC4QrDQC/vNiC4XvLt
8amGObkj80qXv4fr0pETzD2VEFoumcwnN4OWvbuXK/M1sFtt73v+8g7Sd3rb
/Ml7BXfKKfd7Q5rAlX0srlv2JYQ01XY1tzbAV/IHl6pTZXCEZirfTheDXlNF
UYCcCFpUGNMvpW9A0bqlZlHGM6gqaqmL4dWBnk3XQjKzFDL3JXVnDNaAWupj
7r6oYljfmJVSlV4F8jkuvRcPFIJJif8XuzsV0KxpraRDF0KEAjtSNrMc/kh+
1hF/Jg88Ez+b0/pegoX5WMDtglwgkZSn+fAChm2Em4Z0c8DjaFl8Yq4IdAq0
qExuNtiNMkpnXJ6DsrEJg2qVBQm0DIMnU6VQKHy8cseHTJi5Sk1OPFsCzyCp
iOSdAV96FrhvWlcEQe15Ar4uH676Sty/ahfApOqE64b5afCPQGZLsVE+DM1S
J4dUeUD/1KyQ5yaAGE5Zpc7mR0A/y+JtT8mBRJMnU863HoD5KY+7imr/wuK3
J+VMFqWCkaemk+B+Fuy1f8pi30qG4K0pKdbOT+EBv8cx0DEJ2Ep/M+/IZ8Ij
gXL6CCER7tXcDLAq5kMiKfALL/cuZO98GrHgVhrkrV/5evJaAiRwf8KDMSFb
O4/DBZn/sV3ruLCQkF/yHwsEsUA=
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk/k7FAgcxgeRwsix2JVjjHKUa1RjUn2/UZGyOXPvdsiR2hzVk3JFpIhW
jibWncy6lauSKZJy51hnrqWZccsQotWP+3ne9/k87x/wks5csjjHSyAQHNb7
w/+HBYn0HzBg6zuv9bDg6YmiUP6HOZC4WJ2VcZ0FyU3xrpRXeTBhzD+ofYsF
t4VvXkraWAiRCSt3/aJYcM04I4NqVASVhHq/ODoLyHZShmUpxTBgZVUc8ZgF
alet/hKSeAae9P2unsUscLjhwziWUQLJ799qnmSu76kewXKLMiA+uk00b2ZB
dhnhyCtyBfg3mideHWDBvdPdlt9kXkB59gGZTzMsWB7ltzy0uxIGzagPU/nY
sHRPIz31RhVsdpRQ6pZmQ6JtgULuIhOEUsL/ydRgg/7MGebS0TdQYEZWIB9m
g9WFmrjU0mowXQpuvfIbG+TlRb7mwVtYsZc7/9qXDXapC2q2n2sBYhr1NRLY
ECp4O5ynsA6eqdxgjJSyYVuVy7J+0gdoL1Q/K97Nhj0dxRkN+Q3gu9T8hLPK
hkKbtJECThM8voORads40FDZ2xLFaIHxylwjTzMO9Iqe+VrLbYWsZJGw2AAO
nLfdvnHAoQ3m9qwxzQo54N/ZPNzT3w58F4JSJkY5kFSnkjLm3wmZplSnEoVx
uN1PS7E70AXZCTm0NadxCKzrnq8g9YBm0d4ZxbRxWBM4z8dU7wPmnCMtgD0O
7h8UPewdP4G316p1AHUCZnX1SN2Og/AHMffxnxETMJnW3iqdMQQrO4ONdcYm
4L5BsUkV7wicPpbROGA0CZQ2Z/UP/v+C9ajVtFrRJASuZN2ylxoDSuJLpydK
U6AVLs592fQZRJ650VeSpqAyXH7ag8GGrJb+T52K0wBrqmYSjHHQVPdXXsyf
hrSZQr3YlknQfSATL6k8AzSV7fRbYjPw3nB3R+bVGeisVrMaNZ+FebqQQUbz
DESudffRcuZgyipf/L7WLES4XCzxUZiH1tmbvtz4WTj7ja6QU8iFWk/rC24C
c1B6V4WkdHoRlt98u34nYA706vNqjHSXoE4kvZSf8AVCiJqJL9RWQHCDbE5C
xBfojIwLMTdahfy7DA9BpXmY2v+o8qfo72DW9fcz7Zp5sLRgZany8mCqSyd5
2ZsLY9+zHI2seVE9qin8lNYCrOm1Rze/4cPhvM8nfJYXQNTSpr/5OD9eD9He
ZNOxCOlBQ6OBXAF8QJwWMXr9FdopxYGGVYJYoDd7JKFqCV64CKvm527G4+9A
qrR1GcKrLSrCmMJ4XJrrorm0Am4nyR71q0SMltAeY1BXYVFNVqz36RbUErdx
+xK1BndkZCMGCeKYLiJa5UcgYFIqidPhKoG5V24ls94R8PKhPAsmRxKbE0pS
+pJ5MDBfV7sxQgqPlG71jArjRbmDOrbWpjI4H8foPRbEh5L2NR6DOr+g6tLj
PL97G3Cv7VM/1X1b8aopsSyzgB/dY+RW5C7Jo4350dVzLAFkizK/cl0U8XPK
WckQmiDGSMaarOWRMOw1JXsocxNmEdtkAolk7PgioHdvpxCaKezUv3xXGTcU
9Qo+/yiMySpdUjzK29Hh5mGLOjoRtfT1nvr0q+Cwom0H2WALfmwbEDKuUENb
xSkHfxMxJLm1Ze54vgMrva/EXnQVx0BXsgvPiAZWvOv5HvVIAu/voQ17y2pj
94dG6aRRSaRupyTTHXWwoJza4GUihXahJQ+qSykYxLahiNdLIzb6XBvZsQud
ac7uT9x/RrpvdLf4q91I+7WrPExVFs/1CCvFXaKiXufshO5mOYwYjunvMaBh
wEWdBWcxBTSJ0toXulcfEymRrOBURaxW0e5zt9uPu/qHnOSbSOiRGTwQnw5o
6/D8d6oEGR9KHmTkCB/E2hPGJuleyliy0ZCbeMoA38eLzduMb0Ov9TcINBti
2/f9odeCVfA/frq4hQ==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(2s\[Sigma])\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.09769797540612125,
11.754825748270298`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwtk3kwFAocx9dGxLIREjk256YIuVb1+2mVUqJdlaNepZKop2J6NWpWckb0
XLlyty85W1lJsSFPiJy5yZUUZbM506uZ95n5zuf/78yH4urFOEMkEAguv/bb
lxICpFc20ZG9ebfyUI06FovSBUknduKqQw5Eju8GjJe1zM4hWeKUmiaj2pqC
nln+A3EZgDpcFakzVDWs0t7Se85pO5KSR40IO1TQJkJ/WyDNAgdZA1I7DqzH
sKGovu6d5tgWaqfa6KeIZ7pJG2K9TJEeIOml16qACdciu2TKjVHf1cVckr4W
8Y331WHdrdhrT3m1q1EOnQKLo6u4hkijRTj3XZZFUy3DlISjBpg+jWSC/hq8
a2I+dFlpC85bKYn+IMgg66y6m9DwZpQZXNVg+WE1Utxbs3Sf6eKIe/+hbB4Z
W1oHJPaUUtFEL2Lf/XZJ1LcwK/Lu08Z/KP0xueIkTNHulBfS0EKqW48a7Zg4
2qtusvC5rYHagTbhvDoxZEu1KrCk1PGqMPdUt70oRsnG2CznUfBaxQzLdlIE
P5J5cwI3NVzvcT1ELkUYz0UpLyp7qaDGk63zHidXIM2x6LrOtvX4iG/w+ZUx
EWWdqz0HDRRxyTmR1LdOCJUtDRwP2Srg1VDGrnQSAVn5RlvehMmjxycjygXF
ZfCxymPwJmTRv3Nvb5vxEiSnUSbaz67BxPih9QyPBQhVUAob/PXDVLvLjVHu
HMxSlaR7ilajf/WYBltxFtwPq3vW/5DCgB9BB4tjBBBSxSgN4pFQg0e0MNWY
gTI3kk5+rjheFGq0ta/lQ5shh0WvEMOcgONMTb9pyPB7P8oSrETO9/xtfTZf
gcw80te0XwS9a2OMre2mYNmsLbKpcgVGJ+23Paz7GcZ+so9aHyIijaM3flxn
ApiMcbYOUQiNKiXl+FbjMLU98YVc5E8I6t/F5LDGoCM89tZB6x+wP0BN8ti7
EbglpZdURl0EQ1GRskqbYTCrz6u2NpoHSs7zgsrO98C9rU3ZcHIW6CqxK638
BuHUUoJqTqEAciObrrn49UOY24Vib9UZuKR8ZaRcrxfCl7t6zXP4MMevL3Fc
1Q0dVVSH0YPT0DlRIV0k2Qnm2lq/ivgKz9JnroF5B6R/LTSLeTsJr2sDgk+E
tgEs69ivyf4EpkdeLlxZbIEXISpfPLM/QoZQjv+X8GbQD5ERPG/8AD3jZp8O
WLwF1iI7wFl+DIJDTV9eEWsEw9bTG+tujMDAvY8vhevr4e5Ojk0FcRhG43os
TZ+9hsn0tua1me9BaqlOiNnyL0wbmVG6jg6CSMPpyC/kGjhXp+bpfLQfDqyz
jeB7VcPySo8VvI29gJaMztAvlcCq7ZoppXSDj42/nmb4SwjuM0912tEJApW8
GS86D5JrtVPHbnTAKZYaK06rHG50NA1197WB8Xnk3v1WBh6OWqIDLq2QFk0u
sJ4ohR6y61yNoBmyZEZp/qJPoeFFz9uI7LfwnZQjbLuHC4VH0ocLJhqBrCuV
+S37CZi0czIb8hvA8bCy7wK1CDQr3BYskutg/Do4r6t5DIFiwSFChbUgvkm4
0OR6ITilfac6fqiBgemDD2pi8kFFRXIuD17BkRKfv/KO5YLD+erYNG4VsL67
O72zewQWX11583sr4Sed7UJzfwhJjgWqubM88JKL4opmsWH+zuaMNN8KmF14
MskmPICFURGmlfELmD/M8Xjsmwl3TnYxlxTKYIQZdnNKJgMelhB2l6uXQm9g
d7NddRq4THWLPWWUAJlEFPQopYKLr3f2vsxi4If5bpboTQbqFYf7EmuegNYf
zfPjrxJB3UmeXpLKAW8R/Yt/dsTD1T2ZmabWjyG1hX3hjfg9CCbd9EoWLYQ8
VhDz/IlYSGmMO2tYngchEclKQ23RUGT3OFAkPge65APK+GeiICnhN9lwu5nf
FCb9NxD+p0Wivd+6ZrDiP6Tqpss=
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk3s01AkcxT0SaSJGsuXRkMdQicR41PdrFEp70qgMaUNFpfKajU46U0ls
FEdUIu8seUbIljHlkVBqTTKYKFEJ/WYYrzRWf+4995577p+fPy7FN5BxVEZK
SurAQn71/0VAasovFUDXpfULJqBi94MouduFEPziIV3JjID0Vzf9zTnF0FiV
6tFvQUA06WJgmnwZpNtJ7kpoBIQ75+RYOT2AwZRqyektBOh5qDtUZ5SDR8xY
vxWdAOqZvXeXkh/CTKcL1duJgAPnQgtcciqBVP1l+/ddC3usW+ERoxpOKspa
jzIIyK+WcuTo1cC8dnSurwcB13z4bnMaj6HsDLuT6U3A7KCc27bNtcC9EtHJ
9ydg5tr67MxzdRDAkWcRgQs8zFKdoikuLMqz3J8fToAt4cud2fEMiq6wSEqX
CNh7siE5s6oedFrurDGNI0Bbe9n0Ah0cWRtivOoWAR6Zk1Tm5yYYynz06X02
AVEK0THSZc1Q0tvicL2UAP06v1nbtBaIIatct6klwPJteU5bSRvc2uh2ZqqV
gDL3rIHS4VeQ//o3WncPAW21Pa+vF7yGBo3tVTMjBPQo+043id9AQYhc2ul5
Ak4wDeT7DnRAUunRC25kIZzvbP/YLeDBq4YUUqORENKaDTOGzncCnbFD/yMI
IVpgneGxtQsEgSor25hCYDfzJ2oo3XC2/3xTXqgQJItPyHKNe+HCc87A7QQh
HG9ZE+Dp9R6SvI/qtpUKQbiJRuF79UOjXPifp94IYTSL92ZlzgeQqFi4FouF
kEAv31knMwATw93f6jRFYN5xxLjl/CfQfTaWKHAUAftH3mVP9SFwCnYR0Vki
MI1RFT959RnknBOCzO6JoDZG+3tAwVforWHGfeeLACRGruSCb3D2fqv9C9Vx
yCLKaEmvR0GcTw5f7joO1oYGKZdVCOC1HwzTShyHznrq3sE9QphyWnJsN38c
4iT8XutCETDvc3pIBhMQ63eqMlRnAlhPLXwawifg8FyKTmGZGAKuhemNvpmA
qquGFF2fKfCRuZUsby4GWmtxg9OmGeC9e2iSniKGSKUNqY+pP8BGNkLXcckk
dMYlR+5x+gkcW8r+zIuTMLblTu2K+HnIc9Q19pSdAjfGlzwjGWms+KYSuSF+
Cobm87yc9sngk+CujiH9aZDQePHtz2RRM3YEvJqnQdnNXdC+Sw6T2g144WEz
kH3hwyBbvBiVQ0TvaBazwDMvZzvUKaD+LJdWND8Lj/1IRiVFihiU0aYy3f0D
YuoZNVe4JLQrOshgPp+DY/v1Alp/KmEXffwypeknTFFXq/RULMc6+76aZ10S
+EtjdWy/lCp+USS1HZaSwrRMyvBbfzISxnddRjOkkLWtmMEdVsN7RVPlAy7S
yC7ZtPFlrDrSptW5W+VkUMvejLnvdw2kFkZHLHkpg2qeDQH9ZqtwxFubczhX
Fm2YFRFGdprYoHivMTJuER5P1PqhFaiNqva8n9VX5PCrMnda7LcGR9iagmNJ
izFRLWmnpJiC6f7ikVXV8pin1KHBVtLDMGrGTBChgK4662xZV9di2rWIUxvo
iphu2KUuvdYAP7w80WNdshRNbWkVoQJDTLXgPs6wXIb/dvQtda6h4rjfwKFk
gRJSjnXkmvxjgkfEB5dzGMuR7a/nJz2wHpsynP1CnFUwwdL6Y8jqjeiUnVu0
Y48qWhmYp6d4mWHdjXVBT0+T0SOq8kZ9lTmaU8hiVrYa4svQ8AETC9T0iYzW
GVuBKWfj+aqczTiT6zib7rYSj3aTdJMDrTDxvaH7JE8DYz8mCrrp1njbfaLB
lbUKd143tYuyscU/qh3/XmShifWGG3uPe2xBneS0uZvq2hiQe6nvZjYgeUcz
+5DyGrytZl9QSLLHE3xrhUpDClbKO4hTvemYFc6a7Nuni8ELb1jc7oDKm1rv
3kjVw/8A77a5FA==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(2p\[Sigma])\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.5242041923705705,
9.888443865601243}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
GraphicsBox[{{}, {},
{GrayLevel[0], LineBox[CompressedData["
1:eJwtk4k3FAgcgB25Irlz5JgIQyJKjvj9RLnaEiVXKdeQ3RAdbJoKUSolyjE5
S02OHCFrY3KU1LjGKIMVQpMtRyb3tPXefu997/sLPoJPiLM/Hw8Pj+dPfzUs
PVZSsN0aFaZ5WnuoWlgpZM3JPLoTxaY/iY5Ja2KajBW1UMwKO+w2LpAzNTD4
/qXBO3mA8TTZ2CuW6tioZdAf5G6BFhfJJySWCOhwQ39HnJk5rn3wxNK+XQ0T
h5MHWDtNcVVUtOn6cFX0Z4ltSA3ZjkHLVofrQpUxPTKpV6puG2o0UtivU5QQ
6eFnR3S34lFm9dbxfgV0j6u83VhliNf8Vs1I2crjdk3DrHSvLXhu4nzo2U45
vGlsOnxSyQBjk7+fkY6SRTJJPYB3RA+TJkvYZpYySAhk3Nf9SxcvnCLORSpK
YxdjUNSuhohGu21eTq6WQn1zk4rwAS30jSm6NioqiVla7+V4NTRRdJ7B1FeU
QCfVTeYRVzXQevML8pUAcSwQZ8iTxdVR9lOjlctjMUyWSXHgFhPQzLBPqHx5
NbLX0uY5AWpokiJqcNRTBIOSlZeUQ1SQFZmgSG8SQjO3inPaO9bj5NThoPUm
gijj0RT8YYsiTtqJNOhUrkJlqy1uB3+TR0e1TsfrO/iRXGJkQE+UQ6q+FD2m
nRcjbIqdaZ9l8IrJ/pXeEB6k5BA+M0nSaNw5C2HAhSvySokfeKQwdd9W6c6G
JZgjKkn2VUhgYbZ3vcyJBQh0VQ9+syKOjMCo2mSbOUhodK65TBNDd/1n7Xp2
HKgNENMuKVqNS4qBdIHIb9BtWE62rhfGvaSeXReZ05B3YWiUzBHEf4MiHxZ7
TMFal0MD7XsEcPwr258i8hW4Jt1J7Q38uDihdIuePQFjPwq8bA/y4fkCxjur
UDa4OH8q0ObjRU6lCFcwYhy+WmQ8l036AU5mOqQ1D0ah51pqzH7bFVDMULmk
wR2BGPHNmbXEJTAPLi4bjhoGkzfFTbZGC2C0aDZ3ed0QVF3VImw4NgfBqcYW
1bmD4LucrlpYygEZD2Kuqd0AJAb8URmuOgvppCjvV8p9cI3b229aOAP0sAZH
XZVe6GkkHhjdPw1XV0WQSna9A1MtzZ9HTEGpr50S/y0m5E6VmqR0fAFQcaBv
XWEAcLWdpKkTIJBxoCsjvgueJ6hMBlPZ4ONfO0PX6wT9BCnO323jwK5L7IK+
diAvFcR6yI1B1iGFLu2rdDBk+Om0Rn+EjL2Ugiq/N3BzZ7lDPd8IdDeeqr5L
eg1fcrs71+UPQVaLN4OW9AqmjUwIvV4fAE4dIQSwmiGoVS3Yw+sfCL4ZIsix
aQKu4HF+mk4/nEk4xhCnNwC5pXe2hsCCAWM5+djQFxA/YJrtbvkeslc+XnDS
owGlRSt7LLoHimTZTpv56yC6p32YNdANm84dT7d8VwvH3TSFBj0ZMHN2qiyP
UQN9a33mX3I64ds8KTyQUw1vn/d13KB2wP5N8kNM0yooPZQ78uRzG/A3337h
l/UUjJnl+W9L3sJjjZDMOrUK2FgfsGhOaQUFzYLVNXVlECccn8Bb2gKZyeDv
fLoU3HO+E93GX0KJvfrew4kloKKyZr4YmmHd0C5zomsRHPi9KTWnqhG8+mUo
PfaPwXzKh7Zg3wDHnT+qhPg+gky3J6pFczRwfcBfaJRTAAvX9fJy/qyH+vPc
PWkr92FxVMDFZttzYFLbTmJkPlw/1uuyLF8L+bmOYUsSefCommd3nXoNMI07
97xoyAHPryzhZ87V4OnqTvZWyAbPP8OpjvmVwEprnhFnUYB4+sA9UemnILxE
1XBpygB1dznr6uxyUPX28dZhpsFZu/z87bZl0C03/rBV5C7Ei10MoQiVgt92
9bJ73qmQ1XaHZFhXDGHLR5Qsum9Dxb6yOIG0QhhCvsUT/smQmf4LKkicDMrY
IHkLeP4n9WHM4omns/X/AY9oqFk=
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk/k7FAgcxgeRwsix2JVjjHKUa1RjUn2/UZGyOXPvdsiR2hzVk3JFpIhW
jibWncy6lauSKZJy51hnrqWZccsQotWP+3ne9/k87x/wks5csjjHSyAQHNb7
w/+HBYn0HzBg6zuv9bDg6YmiUP6HOZC4WJ2VcZ0FyU3xrpRXeTBhzD+ofYsF
t4VvXkraWAiRCSt3/aJYcM04I4NqVASVhHq/ODoLyHZShmUpxTBgZVUc8ZgF
alet/hKSeAae9P2unsUscLjhwziWUQLJ799qnmSu76kewXKLMiA+uk00b2ZB
dhnhyCtyBfg3mideHWDBvdPdlt9kXkB59gGZTzMsWB7ltzy0uxIGzagPU/nY
sHRPIz31RhVsdpRQ6pZmQ6JtgULuIhOEUsL/ydRgg/7MGebS0TdQYEZWIB9m
g9WFmrjU0mowXQpuvfIbG+TlRb7mwVtYsZc7/9qXDXapC2q2n2sBYhr1NRLY
ECp4O5ynsA6eqdxgjJSyYVuVy7J+0gdoL1Q/K97Nhj0dxRkN+Q3gu9T8hLPK
hkKbtJECThM8voORads40FDZ2xLFaIHxylwjTzMO9Iqe+VrLbYWsZJGw2AAO
nLfdvnHAoQ3m9qwxzQo54N/ZPNzT3w58F4JSJkY5kFSnkjLm3wmZplSnEoVx
uN1PS7E70AXZCTm0NadxCKzrnq8g9YBm0d4ZxbRxWBM4z8dU7wPmnCMtgD0O
7h8UPewdP4G316p1AHUCZnX1SN2Og/AHMffxnxETMJnW3iqdMQQrO4ONdcYm
4L5BsUkV7wicPpbROGA0CZQ2Z/UP/v+C9ajVtFrRJASuZN2ylxoDSuJLpydK
U6AVLs592fQZRJ650VeSpqAyXH7ag8GGrJb+T52K0wBrqmYSjHHQVPdXXsyf
hrSZQr3YlknQfSATL6k8AzSV7fRbYjPw3nB3R+bVGeisVrMaNZ+FebqQQUbz
DESudffRcuZgyipf/L7WLES4XCzxUZiH1tmbvtz4WTj7ja6QU8iFWk/rC24C
c1B6V4WkdHoRlt98u34nYA706vNqjHSXoE4kvZSf8AVCiJqJL9RWQHCDbE5C
xBfojIwLMTdahfy7DA9BpXmY2v+o8qfo72DW9fcz7Zp5sLRgZany8mCqSyd5
2ZsLY9+zHI2seVE9qin8lNYCrOm1Rze/4cPhvM8nfJYXQNTSpr/5OD9eD9He
ZNOxCOlBQ6OBXAF8QJwWMXr9FdopxYGGVYJYoDd7JKFqCV64CKvm527G4+9A
qrR1GcKrLSrCmMJ4XJrrorm0Am4nyR71q0SMltAeY1BXYVFNVqz36RbUErdx
+xK1BndkZCMGCeKYLiJa5UcgYFIqidPhKoG5V24ls94R8PKhPAsmRxKbE0pS
+pJ5MDBfV7sxQgqPlG71jArjRbmDOrbWpjI4H8foPRbEh5L2NR6DOr+g6tLj
PL97G3Cv7VM/1X1b8aopsSyzgB/dY+RW5C7Jo4350dVzLAFkizK/cl0U8XPK
WckQmiDGSMaarOWRMOw1JXsocxNmEdtkAolk7PgioHdvpxCaKezUv3xXGTcU
9Qo+/yiMySpdUjzK29Hh5mGLOjoRtfT1nvr0q+Cwom0H2WALfmwbEDKuUENb
xSkHfxMxJLm1Ze54vgMrva/EXnQVx0BXsgvPiAZWvOv5HvVIAu/voQ17y2pj
94dG6aRRSaRupyTTHXWwoJza4GUihXahJQ+qSykYxLahiNdLIzb6XBvZsQud
ac7uT9x/RrpvdLf4q91I+7WrPExVFs/1CCvFXaKiXufshO5mOYwYjunvMaBh
wEWdBWcxBTSJ0toXulcfEymRrOBURaxW0e5zt9uPu/qHnOSbSOiRGTwQnw5o
6/D8d6oEGR9KHmTkCB/E2hPGJuleyliy0ZCbeMoA38eLzduMb0Ov9TcINBti
2/f9odeCVfA/frq4hQ==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(1d\[Sigma])\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-1.2940511499132061`,
11.831716835608448`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"]
}, Open ]],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"FindRoot", "::", "lstol"}], "MessageName"], ":",
" ", "\<\"The line search decreased the step size to within tolerance \
specified by AccuracyGoal and PrecisionGoal but was unable to find a \
sufficient decrease in the merit function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindRoot/lstol\\\", ButtonNote -> \
\\\"FindRoot::lstol\\\"]\\)\"\>"}]], "Message", "MSG"],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"General", "::", "stop"}], "MessageName"], ":",
" ", "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(FindRoot :: \
lstol\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG"],
Cell[CellGroupData[{
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0.6, 0.4, 0.2], LineBox[CompressedData["
1:eJwtk4s3EwocgD3yqBZhSZbHUmZJiR6Efj+9uD3Io9qibtwiuSLc6uqxrrgp
pKTyinlEu6yYpKUhUmrLaLgxi8hrlYSlpezWOfc75zvfX/CR/UM9D6goKSn5
/PBnj6TF6KgL1yMJ9oRvSKRgucZ6Wca+dehmooXTZIswlejMKiI4I+NCZbFh
1EIMzo/uupYLmJIvjAswNMM6inVnEN0Jt9S5qQ2LyLj54jLH2DUO6JS4ZECc
b4rxPcmSjnX2qCmW9sv1TfBAB2HB1dDVKC9wMz9BMMK0P5PadatWorD98j53
CgnxRcTxXssVSKjamP9o/zykx5Zfqbtng9Nb5eoetXNxtblNVprvcsSTw/Nv
O+njpVX2PeEka+x52kXkthCREWgWoNxrhV5iVZ0zsXpIPijKt3xgiQfcYpmO
brr4UtQ105VLxQHTUT83Kx1c5mBXFiGhYP2to1nW5NmYRXmlr7zQHOXRPFnL
di3cbrLEIfLCQtSJlosmowlYoCUyYGiZYZ7CeM/D1hmYTEzZPMUm43uz2VXS
jdNxSLvmiyzAFN+op4YebtLAoGSjSaNQYzxZalQREq6Oa2hlJy0c5+ORWsbg
IaoaEnc/Du5ebojdjCXa28ZU0ch5OW3HNgOMC4iYQ29UQcZtW+sX8foYYbgp
1oyrjJEb2J41UiImuNZ3n+EoYSaTLG0N1EP5KrnVx3AFnDcgxXcr6aLXgwek
oPTvMEEl6YjLZmMTQxAyrWMSDu40C+Z/10IxX3LhrPVXiKvz5P5dQ8DP7Hfz
BZlfoDKAYHG7eAbOJQk8uowmoMWGw1hfrYmlw34TIo4Mcs+86WPI1HFm4fuJ
5t3joO21SyLcqoaK71GdWnPGYMquJUlYq4rrIkVEfvcn6FcU+LrsUMHi9OQ2
d94IeHkOFlioKCPJmv+2w+sjDDul8+YkKUBANNAe3/IB2hKunvVw+Q6TX6M0
c/zewVmtpRmV1Emo5Ur52VeGwI7PfuxiK4c+ov2v5d0DcO8ChbzAbwIusiWG
bJd++O1bmklRiQzsrtGB9vQtxAeElEeYjEMTt5afuLcXEqbaO+2LRkHDWHx+
MaEH2uqo3n0en+BUVOEGHWE32FPMfxwxArOlYqedW7sgZ6TELqXpAwiBpBCa
SgCmLLbrsd6Bv28u8fQ8MfDijD8Gs4ZAVTPGfu+KdlgWpyt72DgAW5nmSadD
/wXGZEHMbv1+qHw/dnKqoRVsRPsXPz/1FkbTLCO0oQUureNsrlbphZ31L8/L
hC/hQ05L89y8NzBu6T6hEdUMn2ztyO2+3RA0ct231KEJgp6bBu/2fQ13n797
0WXeCFPqh1RrFndC8a7Ry/8o+MBoaB/nkjsgNWUpaZbiGZyT2GfT174Ce8PM
5GMLGiCzgZLdf6oNUp1DfBiBT+BUm7CnQ9IClmOHNV/zH8MhmrlGl48I+uYT
qGHb6kCs7f/liawZ+M/GP9tKH4GAJ266yGqCsgbBexdmDZTsyum9I20Ea9Jg
+ktxFaxq5eQJbguAxaNs7Mt9CIuqA746ZD6HK4a7wu4mP4BYzXNxyiUNoLej
zdbv1n2gMz9TaQNP4I+9N5asfXMPjI1nfWFDPZA1HK8K1pSD9++PrzLv1UF9
bkwkl1MGDiP+NfJfaiEsKba/YxMHMmh3TIonasB5xf5oK1kJyBOtcpknqkF0
OYWXd+wOfO1T89qwkgdD3neCiTZsSPRr9/pmUAnOtH6PMP0iuFWhtKnKjAsm
TRzXcjILfIY7NO97VkDhEVHdXI9C8DkRwdqSVw5UVjNRkXMTqEe9b8zUuwuB
FeGDyrr5YEbXX1+RzYGQSP63nqxcOO6al7fapRRCXh+bdXFTDpwj/BWaqVEC
umH0MaIGE7IarwXaVLHhwMd9UjrvBpS5l8aqpRZBpGt9Au16BmSk/YQFN9cm
+WunpYHS/7Qy1TpHGpnV/wHUF7d/
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk/k7FAgcxgeRwsix2JVjjHKUa1RjUn2/UZGyOXPvdsiR2hzVk3JFpIhW
jibWncy6lauSKZJy51hnrqWZccsQotWP+3ne9/k87x/wks5csjjHSyAQHNb7
w/+HBYn0HzBg6zuv9bDg6YmiUP6HOZC4WJ2VcZ0FyU3xrpRXeTBhzD+ofYsF
t4VvXkraWAiRCSt3/aJYcM04I4NqVASVhHq/ODoLyHZShmUpxTBgZVUc8ZgF
alet/hKSeAae9P2unsUscLjhwziWUQLJ799qnmSu76kewXKLMiA+uk00b2ZB
dhnhyCtyBfg3mideHWDBvdPdlt9kXkB59gGZTzMsWB7ltzy0uxIGzagPU/nY
sHRPIz31RhVsdpRQ6pZmQ6JtgULuIhOEUsL/ydRgg/7MGebS0TdQYEZWIB9m
g9WFmrjU0mowXQpuvfIbG+TlRb7mwVtYsZc7/9qXDXapC2q2n2sBYhr1NRLY
ECp4O5ynsA6eqdxgjJSyYVuVy7J+0gdoL1Q/K97Nhj0dxRkN+Q3gu9T8hLPK
hkKbtJECThM8voORads40FDZ2xLFaIHxylwjTzMO9Iqe+VrLbYWsZJGw2AAO
nLfdvnHAoQ3m9qwxzQo54N/ZPNzT3w58F4JSJkY5kFSnkjLm3wmZplSnEoVx
uN1PS7E70AXZCTm0NadxCKzrnq8g9YBm0d4ZxbRxWBM4z8dU7wPmnCMtgD0O
7h8UPewdP4G316p1AHUCZnX1SN2Og/AHMffxnxETMJnW3iqdMQQrO4ONdcYm
4L5BsUkV7wicPpbROGA0CZQ2Z/UP/v+C9ajVtFrRJASuZN2ylxoDSuJLpydK
U6AVLs592fQZRJ650VeSpqAyXH7ag8GGrJb+T52K0wBrqmYSjHHQVPdXXsyf
hrSZQr3YlknQfSATL6k8AzSV7fRbYjPw3nB3R+bVGeisVrMaNZ+FebqQQUbz
DESudffRcuZgyipf/L7WLES4XCzxUZiH1tmbvtz4WTj7ja6QU8iFWk/rC24C
c1B6V4WkdHoRlt98u34nYA706vNqjHSXoE4kvZSf8AVCiJqJL9RWQHCDbE5C
xBfojIwLMTdahfy7DA9BpXmY2v+o8qfo72DW9fcz7Zp5sLRgZany8mCqSyd5
2ZsLY9+zHI2seVE9qin8lNYCrOm1Rze/4cPhvM8nfJYXQNTSpr/5OD9eD9He
ZNOxCOlBQ6OBXAF8QJwWMXr9FdopxYGGVYJYoDd7JKFqCV64CKvm527G4+9A
qrR1GcKrLSrCmMJ4XJrrorm0Am4nyR71q0SMltAeY1BXYVFNVqz36RbUErdx
+xK1BndkZCMGCeKYLiJa5UcgYFIqidPhKoG5V24ls94R8PKhPAsmRxKbE0pS
+pJ5MDBfV7sxQgqPlG71jArjRbmDOrbWpjI4H8foPRbEh5L2NR6DOr+g6tLj
PL97G3Cv7VM/1X1b8aopsSyzgB/dY+RW5C7Jo4350dVzLAFkizK/cl0U8XPK
WckQmiDGSMaarOWRMOw1JXsocxNmEdtkAolk7PgioHdvpxCaKezUv3xXGTcU
9Qo+/yiMySpdUjzK29Hh5mGLOjoRtfT1nvr0q+Cwom0H2WALfmwbEDKuUENb
xSkHfxMxJLm1Ze54vgMrva/EXnQVx0BXsgvPiAZWvOv5HvVIAu/voQ17y2pj
94dG6aRRSaRupyTTHXWwoJza4GUihXahJQ+qSykYxLahiNdLIzb6XBvZsQud
ac7uT9x/RrpvdLf4q91I+7WrPExVFs/1CCvFXaKiXufshO5mOYwYjunvMaBh
wEWdBWcxBTSJ0toXulcfEymRrOBURaxW0e5zt9uPu/qHnOSbSOiRGTwQnw5o
6/D8d6oEGR9KHmTkCB/E2hPGJuleyliy0ZCbeMoA38eLzduMb0Ov9TcINBti
2/f9odeCVfA/frq4hQ==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(1p(+\[Pi]))\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.03282883344474078,
11.7759128321482}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0.6, 0.4, 0.2], LineBox[CompressedData["
1:eJwtk4s3EwocgD3yqBZhSZbHUmZJiR6Efj+9uD3Io9qibtwiuSLc6uqxrrgp
pKTyinlEu6yYpKUhUmrLaLgxi8hrlYSlpezWOfc75zvfX/CR/UM9D6goKSn5
/PBnj6TF6KgL1yMJ9oRvSKRgucZ6Wca+dehmooXTZIswlejMKiI4I+NCZbFh
1EIMzo/uupYLmJIvjAswNMM6inVnEN0Jt9S5qQ2LyLj54jLH2DUO6JS4ZECc
b4rxPcmSjnX2qCmW9sv1TfBAB2HB1dDVKC9wMz9BMMK0P5PadatWorD98j53
CgnxRcTxXssVSKjamP9o/zykx5Zfqbtng9Nb5eoetXNxtblNVprvcsSTw/Nv
O+njpVX2PeEka+x52kXkthCREWgWoNxrhV5iVZ0zsXpIPijKt3xgiQfcYpmO
brr4UtQ105VLxQHTUT83Kx1c5mBXFiGhYP2to1nW5NmYRXmlr7zQHOXRPFnL
di3cbrLEIfLCQtSJlosmowlYoCUyYGiZYZ7CeM/D1hmYTEzZPMUm43uz2VXS
jdNxSLvmiyzAFN+op4YebtLAoGSjSaNQYzxZalQREq6Oa2hlJy0c5+ORWsbg
IaoaEnc/Du5ebojdjCXa28ZU0ch5OW3HNgOMC4iYQ29UQcZtW+sX8foYYbgp
1oyrjJEb2J41UiImuNZ3n+EoYSaTLG0N1EP5KrnVx3AFnDcgxXcr6aLXgwek
oPTvMEEl6YjLZmMTQxAyrWMSDu40C+Z/10IxX3LhrPVXiKvz5P5dQ8DP7Hfz
BZlfoDKAYHG7eAbOJQk8uowmoMWGw1hfrYmlw34TIo4Mcs+86WPI1HFm4fuJ
5t3joO21SyLcqoaK71GdWnPGYMquJUlYq4rrIkVEfvcn6FcU+LrsUMHi9OQ2
d94IeHkOFlioKCPJmv+2w+sjDDul8+YkKUBANNAe3/IB2hKunvVw+Q6TX6M0
c/zewVmtpRmV1Emo5Ur52VeGwI7PfuxiK4c+ov2v5d0DcO8ChbzAbwIusiWG
bJd++O1bmklRiQzsrtGB9vQtxAeElEeYjEMTt5afuLcXEqbaO+2LRkHDWHx+
MaEH2uqo3n0en+BUVOEGHWE32FPMfxwxArOlYqedW7sgZ6TELqXpAwiBpBCa
SgCmLLbrsd6Bv28u8fQ8MfDijD8Gs4ZAVTPGfu+KdlgWpyt72DgAW5nmSadD
/wXGZEHMbv1+qHw/dnKqoRVsRPsXPz/1FkbTLCO0oQUureNsrlbphZ31L8/L
hC/hQ05L89y8NzBu6T6hEdUMn2ztyO2+3RA0ct231KEJgp6bBu/2fQ13n797
0WXeCFPqh1RrFndC8a7Ry/8o+MBoaB/nkjsgNWUpaZbiGZyT2GfT174Ce8PM
5GMLGiCzgZLdf6oNUp1DfBiBT+BUm7CnQ9IClmOHNV/zH8MhmrlGl48I+uYT
qGHb6kCs7f/liawZ+M/GP9tKH4GAJ266yGqCsgbBexdmDZTsyum9I20Ea9Jg
+ktxFaxq5eQJbguAxaNs7Mt9CIuqA746ZD6HK4a7wu4mP4BYzXNxyiUNoLej
zdbv1n2gMz9TaQNP4I+9N5asfXMPjI1nfWFDPZA1HK8K1pSD9++PrzLv1UF9
bkwkl1MGDiP+NfJfaiEsKba/YxMHMmh3TIonasB5xf5oK1kJyBOtcpknqkF0
OYWXd+wOfO1T89qwkgdD3neCiTZsSPRr9/pmUAnOtH6PMP0iuFWhtKnKjAsm
TRzXcjILfIY7NO97VkDhEVHdXI9C8DkRwdqSVw5UVjNRkXMTqEe9b8zUuwuB
FeGDyrr5YEbXX1+RzYGQSP63nqxcOO6al7fapRRCXh+bdXFTDpwj/BWaqVEC
umH0MaIGE7IarwXaVLHhwMd9UjrvBpS5l8aqpRZBpGt9Au16BmSk/YQFN9cm
+WunpYHS/7Qy1TpHGpnV/wHUF7d/
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk/k7FAgcxgeRwsix2JVjjHKUa1RjUn2/UZGyOXPvdsiR2hzVk3JFpIhW
jibWncy6lauSKZJy51hnrqWZccsQotWP+3ne9/k87x/wks5csjjHSyAQHNb7
w/+HBYn0HzBg6zuv9bDg6YmiUP6HOZC4WJ2VcZ0FyU3xrpRXeTBhzD+ofYsF
t4VvXkraWAiRCSt3/aJYcM04I4NqVASVhHq/ODoLyHZShmUpxTBgZVUc8ZgF
alet/hKSeAae9P2unsUscLjhwziWUQLJ799qnmSu76kewXKLMiA+uk00b2ZB
dhnhyCtyBfg3mideHWDBvdPdlt9kXkB59gGZTzMsWB7ltzy0uxIGzagPU/nY
sHRPIz31RhVsdpRQ6pZmQ6JtgULuIhOEUsL/ydRgg/7MGebS0TdQYEZWIB9m
g9WFmrjU0mowXQpuvfIbG+TlRb7mwVtYsZc7/9qXDXapC2q2n2sBYhr1NRLY
ECp4O5ynsA6eqdxgjJSyYVuVy7J+0gdoL1Q/K97Nhj0dxRkN+Q3gu9T8hLPK
hkKbtJECThM8voORads40FDZ2xLFaIHxylwjTzMO9Iqe+VrLbYWsZJGw2AAO
nLfdvnHAoQ3m9qwxzQo54N/ZPNzT3w58F4JSJkY5kFSnkjLm3wmZplSnEoVx
uN1PS7E70AXZCTm0NadxCKzrnq8g9YBm0d4ZxbRxWBM4z8dU7wPmnCMtgD0O
7h8UPewdP4G316p1AHUCZnX1SN2Og/AHMffxnxETMJnW3iqdMQQrO4ONdcYm
4L5BsUkV7wicPpbROGA0CZQ2Z/UP/v+C9ajVtFrRJASuZN2ylxoDSuJLpydK
U6AVLs592fQZRJ650VeSpqAyXH7ag8GGrJb+T52K0wBrqmYSjHHQVPdXXsyf
hrSZQr3YlknQfSATL6k8AzSV7fRbYjPw3nB3R+bVGeisVrMaNZ+FebqQQUbz
DESudffRcuZgyipf/L7WLES4XCzxUZiH1tmbvtz4WTj7ja6QU8iFWk/rC24C
c1B6V4WkdHoRlt98u34nYA706vNqjHSXoE4kvZSf8AVCiJqJL9RWQHCDbE5C
xBfojIwLMTdahfy7DA9BpXmY2v+o8qfo72DW9fcz7Zp5sLRgZany8mCqSyd5
2ZsLY9+zHI2seVE9qin8lNYCrOm1Rze/4cPhvM8nfJYXQNTSpr/5OD9eD9He
ZNOxCOlBQ6OBXAF8QJwWMXr9FdopxYGGVYJYoDd7JKFqCV64CKvm527G4+9A
qrR1GcKrLSrCmMJ4XJrrorm0Am4nyR71q0SMltAeY1BXYVFNVqz36RbUErdx
+xK1BndkZCMGCeKYLiJa5UcgYFIqidPhKoG5V24ls94R8PKhPAsmRxKbE0pS
+pJ5MDBfV7sxQgqPlG71jArjRbmDOrbWpjI4H8foPRbEh5L2NR6DOr+g6tLj
PL97G3Cv7VM/1X1b8aopsSyzgB/dY+RW5C7Jo4350dVzLAFkizK/cl0U8XPK
WckQmiDGSMaarOWRMOw1JXsocxNmEdtkAolk7PgioHdvpxCaKezUv3xXGTcU
9Qo+/yiMySpdUjzK29Hh5mGLOjoRtfT1nvr0q+Cwom0H2WALfmwbEDKuUENb
xSkHfxMxJLm1Ze54vgMrva/EXnQVx0BXsgvPiAZWvOv5HvVIAu/voQ17y2pj
94dG6aRRSaRupyTTHXWwoJza4GUihXahJQ+qSykYxLahiNdLIzb6XBvZsQud
ac7uT9x/RrpvdLf4q91I+7WrPExVFs/1CCvFXaKiXufshO5mOYwYjunvMaBh
wEWdBWcxBTSJ0toXulcfEymRrOBURaxW0e5zt9uPu/qHnOSbSOiRGTwQnw5o
6/D8d6oEGR9KHmTkCB/E2hPGJuleyliy0ZCbeMoA38eLzduMb0Ov9TcINBti
2/f9odeCVfA/frq4hQ==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox["\"p(1p(-\[Pi]))\"", TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 39.61737269377307}, {-0.03282883344474078,
11.7759128321482}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Print"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[1, 0.5, 0], LineBox[CompressedData["
1:eJwtk4k3FAgcgJErBjkaIyWDjHGUs4j29xtqlWorVK7KVEaaSjmWVhq2REOS
pZLc8rLT1CCkcczS6cwxFRKZTUSSqHE0bb2333vf+/6Cj7w/2D1ARkpKyveH
P3si46y6fJsLbtD8Q2Q8aYjlCi7Tmf7OOKJPK8oqNsCrWrRiDoGGe7tVYyIi
ycgs/LP/cj4gxcxwpcwefWygWL4K8l6Hb+K/OXpx9NAteZVT3FpHzOUFB8nV
L8XEwdS+HmcHzNIknqXPLcGAHoJBevAapIqTCDOeOphx8mK3Rq0dkh5z5uOb
tRFbQiNFZrZox35KH91HRO+48r8aKqzxgl/GDiWVxbjG2Do7w88KxWwLM367
JqasdhgM0bVE1sFeo1KOBrICDRnSIgucpr/YeO6qOpIPdRaa3TfD9ixiV0v6
Iuzo7FfeWEVFG2mqkO2shqsc7ctC+yjo1NtCGj6hgtmUl0RpI2P8ELAwjF2h
jNuXmzuGsY3Qo7Q5baWqEhapdpJYqoZ44iXHOTxCEVO10twkXDLG+XaM13yS
xxE1gXiaoY/HjdJpQ6fkMCh12dyyYD2cyjs1KKcji2u9yk6ZOC1Fh/YqhuMj
GdTyecAcsFqCVyxtG1POSuMympXXzq0kzHQPWDLpIYWs2zaWLYlEFDGjn/oe
kkDYeq674L0WXlnPx8lN83A9l/xeGKiJcWO7eEObZ+E8STdxQEoD7fTZUx5H
xfCVqqveW7YIvx/VOm/H/QKHdhkym76pYjXPdRdTaRoSGtyrzgkIeIlbTas5
8xn4DILJ7VtKWKPcrNOvPQld1qUslzpFnCxhkMMfTkB+zJu3rGl5bNQ6EmRN
/AhqHrv72rbIoUHnsSRB6xhI7LsuttUvQLrPC25dyXsY+l7k57pTBsd3OGkx
qofBw324yERGGoVdfmqS0SEYX3etZvHF77AtiDE7tu4tPE9KP7PD9Rt0JIl9
7t0RwRnVlZl86hzo7AuxXAiDYN/EfeBqMwN9Cdoc7sgAVLApZAP6V+Ao3Ta1
8e2HA/MZyzm8abBcMhnVpN4HiYyj5aHLp8CU6eGd9bkHkiTdrxw4k+A35VIX
8O0lPG+ger7d8Qn2pCmfJ5u/AAeK8Y8jJsB8y16jhdFCyJvg2ac9+wCtmO6U
9KETQGKyXbN4FGTpFbFbozqgJkHvI7N4BMT8iTw9o3ZYlaAxXd36DhSt02Vn
mtqANVd01oc4BP461yK1T7eAdedB08bofyGRUthg6tsEKc6lbnUyIhDxNysp
+j+FD3ld7doFbyCF172g/vxj+GRjT+72GwDZNo5sofAhBDXqM338XsORqLgI
C9oDkMgfXiAwfQX9BqRqlcZ6YD3pnqoi90D+hhC6/rF/IL7PIcf7l5cwR8wr
CzITwPUnlJyh6OdQukiTSpGqhejnbYM9fV0wdeuklbCLD4e9jBX6fTuB/Fty
Q0V7FfSq7Rc/mm4Hzl43/4HPldBc0/ssufgZjJlzes1NYoG3O090530rkEtC
xvnKsbBaWFrQfLsZHBPdXptIYmBFHWPW8Xoj0LwtQvcoxEKcYnyCNO8JXArP
OVxxlQXeuV+oXu8egV4/ycJhezTo6amIufAQRJTYh2KrKPA88iA9t6IBpJNN
yt+tPQmOE/sFM5vqoVObILE8FAGZXneW3/oqgATbTQTbynCYuWCRnxtVB4Hm
SqMxK8Jg9q2cx3q7GvjjhqqgixcCF+jdHvMkPrBe30gg7j4BNyulfq01rIKR
rA4NT+Jx8B3vUbznXgmVKxL2nO4+Cr5RocWbC8qhfjjFZFcJE6i/e2Ypa94F
di3hk2ZJEBh6E10qc0ohS/dKrlAYCJEbCwrWuJaA0d+nI8eXMiCeEBt8XYEH
9+/arhadOQjZrZcDrWu5wG6vZjEUD0DZtpI4uasc0MhZmbrzJh0yM35SDALS
DYtyuj9I/Q9t8TBpJju77j98g6Ef
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVk3s01Akcxb1Jk9dIymMMwlCJpJlU3y8KpT2JHkPapQdKRczWtHGmksfm
1RFF5J2SZ4SsjCm0Qsma1HiUUnohdX4/Y9BY++fec++55/75+ePSD4R4HpaT
kZHZN5//+v8iIDPjPxXDqwsr501A9Y670YrpJXDyyT0nNRsCsp9dDbTll0Fr
bab3kB0BsZTzIVnKlZC9QXpDyiSA61ZQsM71LnzIqJOe2EiAibeOc11OFXjH
jQ+tcyKAcWrXjYXUeyDpdWf4uRKw72x4sXtBDVDqPm35tn1+j/ep3Pesg2Oq
8qwxTwJu18m48E3qYc4wtvCANwGJ/iKvWd0GqDzF62X7ETD9QdFr89pGEMRE
9IoCCZAkrszPPdsEwXxlzkTIPA+7glYqFoBCkf2e21wCHCYOCCRbH0FpDIei
doGAXcda0nJrm4HWft3IOoEAQ8NFU/N0cMg0zHLZNQK8cycZ7I+PYST3/vvX
+QREq8TGyVa2QflAu3NSBQHLmwKmHbLaIY6qmbS+kQD7F1UFneWdcG211ylx
BwGVe/OGK748g9vPlzL7+gnobOx/nlT8HFp0t9RKRgnoVz8w9ZjshuIwxawT
cwQcZZspv9nXA6kVh895UUmI7O161zcohGctGZRWCxKy2sxzRiJ7wclz6/J3
QELsICvHe9MrGAzRXNLJJoHXJiLq6X1wZijycVE4CVKlo/ICywE49zd/OP0y
CUfajYJ9fF9Dqt9h484KEr6vYdJFvkPQqsj9/Xg3CWN5wu4lBW9BqmnnUUaS
cNmpaluT3DAQX/q+NulPgm3PIcv2yPdg/Gg8ZdBlEngzRRd9dEbA9aT7DyfO
JFjHaZEPnn0ERbfLoTY3J6ExzvBbcPFnGKhnJ3wTTQJILTyoxV/hzJ0Oxyda
YsibqGSmPh8D8jaVq+EhBpa5WcZFzQkQdu0/bZAiht5mxq4PO7+D2HVB0A6R
GBKkogFWyQ9g3+H3U8ymID7geE04jQDOQzv/Fu4UHJzNoJVUkhCceNpkrHsK
ai+Z0439xeAvdy1N2VYCzI6yFtc1EhC+vGeVnSGBKLVVmQ2MGVgvH2HssmAa
ehPSona6/gS+A31P7vlpGN94vXFx8hwUuRhb+sjPgJfnpyILOVms/qoZtSp5
Bkbminxdd8vhg5OvekaWz4KUKUzueiSP+vGj4Ns2C+peewe7titiapeZkHv6
J+Sfe/uBRyqhetiPl0w7KQhtq3jOTSq4fFrALJ2TQkMAxaK8VBVDczo1p/rm
IK7Zsz5GQEHD3lWM0oMyGLTHJLjjpxoeyzvzSGNYBsUMPc3+ag0cyv3m639M
Fv/U1YsfktHCT6qUzoMycpiVS//yIpCKE5Y33Mdy5JCzucxT8EUbb5aKq4bd
5ZFXvmb103gdZE7pCDYpKqCBow179y+6yCiJjVjwVAG1fVqCh2yW4aifIf9g
oSKuZ1dHWGzQxxbVm61RCUp4JMVgxiDEELUchT/rYpTxs7pgigwwwlGe/mBQ
qgqmaKduk5bRMTuQHF1WtwCL1Hp0eWomeJqRIwmdUEUP2goHziVTzEqMOL7K
iYLZ5q90ZE3N8O3To/2s8kVo7cCsDh80x0w7QUOOvTr+0/NmoVs9A/+o/l7R
fkID6UE9hVZ/WeEhcr8G31MTeYEmAbLDK/FxjltAmJsWXrZnvQvTW42u+YWl
W3dScZ2ZbXaGrw02XVkR+vCENnpH11xprrVFWzqV5OQvRnwazh22skN9/6hY
2rgOZpxJFmnx16Kk0GU620sXD/dRjNNC1mHKa/O9k8KlGP8uZbDPiYXpe4kW
D44ebkuy3hC93gF/rXO5pWBngM3mqweOeG9EWlrW7FUdGgYXXnhzNR+QurWN
91uyEaZrOxaXUBzxqIilUnOLjjXKzmSmnxPmcTmTb14Y48n5Nyh1OaP6mo4b