-
Notifications
You must be signed in to change notification settings - Fork 1
/
homonuclear_levels.nb
6387 lines (6225 loc) · 281 KB
/
homonuclear_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[ 287200, 6378]
NotebookOptionsPosition[ 280459, 6154]
NotebookOutlinePosition[ 280880, 6170]
CellTagsIndexPosition[ 280837, 6167]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[TextData[{
"Electronic Energy Levels of ",
Cell[BoxData[
FormBox[
SubsuperscriptBox["H", "2", "+"], TraditionalForm]]],
"-like Atomic Dimers"
}], "Title"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "diatomic_base.m"}], ";"}]}], "Input"],
Cell[CellGroupData[{
Cell[TextData[{
"Calculate energy ",
Cell[BoxData[
FormBox[
RowBox[{
RowBox[{"E", "/",
SuperscriptBox["Z", "2"]}], "\[Equal]",
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{"p", "/",
RowBox[{"(",
RowBox[{"Z", " ", "R"}], ")"}]}], ")"}], "2"]}]}], TraditionalForm]]],
" in dependence of ",
Cell[BoxData[
FormBox[
RowBox[{"Z", " ", "R"}], TraditionalForm]]]
}], "Section"],
Cell[BoxData[
RowBox[{
RowBox[{"CalcP", "[",
RowBox[{
RowBox[{"{",
RowBox[{"l_", ",", "m_"}], "}"}], ",", "pstart_", ",", "ZR_List"}],
"]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{"p0", "=", "pstart"}], "}"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{"p0", "=",
RowBox[{
RowBox[{"FindRoot", "[",
RowBox[{
RowBox[{"RadialPZero", "[",
RowBox[{
RowBox[{"{",
RowBox[{"l", ",", "m", ",", "p"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"2", " ", "#"}], ",", "0"}],
RowBox[{"(*", "homonuclear", "*)"}], "}"}], ",", "10"}], "]"}],
",",
RowBox[{"{",
RowBox[{"p", ",", "p0"}], "}"}]}], "]"}], "\[LeftDoubleBracket]",
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}]}], "}"}],
"&"}], "/@", "ZR"}]}], "]"}]}]], "Input"],
Cell[BoxData[
RowBox[{
RowBox[{"CalcE", "[",
RowBox[{"lm_", ",", "pstart_", ",", "ZR_List"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"First", "[", "#", "]"}], ",",
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"Last", "[", "#", "]"}], "/",
RowBox[{"First", "[", "#", "]"}]}], ")"}], "2"]}]}], "}"}], "&"}], "/@",
RowBox[{"CalcP", "[",
RowBox[{"lm", ",", "pstart", ",", "ZR"}], "]"}]}]}]], "Input"],
Cell[BoxData[
RowBox[{
RowBox[{"CalcErepul", "[",
RowBox[{"lm_", ",", "pstart_", ",", "ZR_List"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"First", "[", "#", "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"Last", "[", "#", "]"}], "/",
RowBox[{"First", "[", "#", "]"}]}], ")"}], "2"]}], "+",
FractionBox["1",
RowBox[{"First", "[", "#", "]"}]]}]}], "}"}], "&"}], "/@",
RowBox[{"CalcP", "[",
RowBox[{"lm", ",", "pstart", ",", "ZR"}], "]"}]}]}]], "Input"],
Cell[BoxData[
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Use", " ", "energy", " ", "limit",
FormBox[
RowBox[{" ",
RowBox[{
RowBox[{
RowBox[{"E", "/",
SuperscriptBox["Z", "2"]}], "\[Equal]",
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{"p", "/",
RowBox[{"(",
RowBox[{"Z", " ", "R"}], ")"}]}], ")"}], "2"]}]}], "\[Rule]",
RowBox[{
RowBox[{"-", "1"}], "/",
RowBox[{"(",
RowBox[{"2",
SuperscriptBox["n", "2"]}], ")"}]}]}]}],
TraditionalForm], " ", "for", " ", "R"}], "\[Rule]",
RowBox[{"\[Infinity]", " ", "as", " ", "starting", " ", "point"}]}], " ",
"*)"}]], "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[CellGroupData[{
Cell[TextData[{
Cell[BoxData[
FormBox[
RowBox[{
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{"p", "/",
RowBox[{"(",
RowBox[{"Z", " ", "R"}], ")"}]}], ")"}], "2"]}], "\[Rule]",
RowBox[{
RowBox[{"-", "1"}], "/",
RowBox[{"(",
RowBox[{"2", " ",
SuperscriptBox["1", "2"]}], ")"}]}]}], TraditionalForm]]],
" for R\[Rule]\[Infinity]"
}], "Subsection"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"p1s\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"CalcP", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[", "ZRrange", "]"}]}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+", "#"}], ")"}], "/", "2"}]}], "}"}], "&"}], "/@",
"ZRrange"}]}], "}"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<p(1\!\(\*SubscriptBox[\(s\[Sigma]\), \(g\)]\))\>\""}], "}"}]}],
",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{"Blue", ",",
RowBox[{"{",
RowBox[{"Gray", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwtxQk02wccwPE/6mgcVcdcCSJIWou77u330zKWTlfGW1dle17jWHVdh7qv
0rWdlnpbjUc3Fdq5qec5t5RoaIsQRx21MIxFba64j63v7fPe932pQVd82NIE
QZz8r7fXlFvpf6D1PvpXBYvT1S3QO9PC5YaTM+5EHKfL2L2LsRkCFUuOAzLu
nd53+ukYnsl98us2yw6/p9eY33KnI/ZExEyb2aLvSKrTvIkJjgt8+HFojWTe
vM6CIw2b21TzhtIs0exRWlZJAhWTQ2jBUtNMvAbekawBA0yVXk0Rk8wwiZnq
17NFRv6urhIlhoHa9XPZ+WxdLKaPvCNlbIqRhJ4AVbUxqjzQ+2cFY0xj8zpl
1jTRhVPByzY0wkffVCeNKWrg4pGnm5JgQ/y8wMrB9xM1TPFJkTytouDXiwNJ
EXxV7OpQIzjjuojrBWur7SpIc7U65+eljWfVO1JCXZUwZ/7YD6f3NFGl/KZ3
I0FCYehlYe6sOh4kn+C5bsnjbW29jElCDb2eicYLaHJYHsAe8GSposmBhLwd
cwgtqh1E43HK6DLQcJ69Ko0twUqMqgoSzuxS517flcKD9tArI0YK6CHOjnF3
I9Cz68Zvo8Oy2E8zj8wx3QfCcTBL0C6DURdXifilHeiJ1WgJnpLCM9PF2aNT
W5BP4/b1MAnUYl/SJ29swNid+2neHntA8cpc8GCuA1/kH6qYtA0dk3aMsdQ1
MKbs95ZsbkAl2XiuYmUFgnfzDMprJOCcqH/nj8RlOMpP0z9VuArJnuq6w9Ql
sMptWqG1L8On4eZn9ZsXwZFumpd+dAm47bEmVy8sQKjdximqxiKQwgwTrBl/
weW+WKvmj8WQdF1lhKw3B9a31CStvXOw9aV8gYrNLGzmTLtyv5sF3x0hVear
adj7k66xmjgNGRdTwpRfTME/Dwf7tThTMJEZruPrPgl6Y9rlo2IRpFiVWFJa
J8CGNXCdNfUaanOGGXIm45DcNbrWRB2DQ9hQSGofgRVW2ReUklfwo18/7U3W
MHSSK2s02UOwzc81/DZjEC6dM5UX+Q+Az756Wm+tEKjbrzZd0/uBMBkRDB/0
Qfxnm4FKdQL48Dwo58/2gt3QE053VTdcE1fbHhzphsh788aktefg28MsrYt9
DoHM8mx+fCeILhyWH6R0goG+8mYlPAPh4dv13I0OKIt+kCnnxIN6e6cVXRke
tKhNKBSx2yBAItYROLfB1l1mUWE8F8gzG0EaLlyIvmqp+DihFUYlEB3/USs8
jiKdKOM3QZdS5qRjURP4/z2m0OjTAOLj1N9t7RtASzhTfFO3HhpTnVVfytbD
3EHQyTDtOjBy2wvvVaiDOE8Ox96jFnhvHr5MdK6FpDj6i/tu1aCmvhWya10N
LSHr+qOsCnBcZtn88l4F5Oe9VQqDRfsfNGWVAvG/AGeGrEgo5P4LqxDtkg==
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVxX1Q03UAx/FBCQuQYAxZEA8DYhJyPE9Ej88XpYNbabI7BaVIvGOQSVZG
wlCYoYkgkgYcHHAh08MYThN64EF5UAfSkGd5ymM/BgzpkCcHAwGL4/dP77v3
vbhHTwhj9BkMxu7/Xvf/qVBYsN4vaFOup0Jt7ILdgKACgZWLLQdHVEgR81pz
g2+jZ2eY3HlZBXGoVLo95Fek+waIT7IpaF4f3f05pxJrHi45P3hRsOoavX7B
+jdw9+4T3wyjEPlikPmn8A+My7MkbgkUyhKM/MoV1Thctl96rIjCqa89jctO
18HKsPKr1mYKS1nupSXJ9VD1OxrVLFKoZT1jlsY0IrzUJqhq2wjKTxVfNgh4
AK28u5PEjsDebrPuFh7hoVSvXls2gih32RVFcjMqYvy6EqdH8O2PE85GLx/D
lFUcfDZQDX7vXalSrkTZ4e9yBnLVSD6kizKpbIeH8XBm9IIa3OU+XdC5TqSv
+aWYRY/iiwgXw+HIbtjkPW5cfDqK5ndv3bGM6YV3iNvxDw6NYU5QfsT2Rh/q
rr4fEfJ8DKktAy+ruYNw+IZ/5sv0cfgIur8XUH8jq6uqROivgc0gRzYwOYwL
fMWx35c1mL7W02klpaDbPzmR2TmB1XEee/6MGinzFhlL959Dl6cOqs8Yw+m/
Zks+ezQJ73SWtu6JBu8EXxRHaf5BfEeSV83Hk+BP1HGe8aYQx1/cw2VPwbfw
1TWzjBfYwXMpOGc+gyKLnz5qYs7AK796zqlpFhzlJ0vm1AzMFWl2e0rmwVQH
ZXM7ZiFaKbCX3dHCZ02ZfaB/Ds62a09u6BZx1TRVI1yZh2I4Ms44ZRnqafeq
TdBi8FJuWljIKuIVrY03pQsodKrvaHNnkA11aEti14ooPbLhEhg7erLbm94g
G75CaMv5+wNPN5ENV/G6Ke5EvyOTSKKrth83Z5BakclWeYURifIxYIcf1CMe
t/2Hh8SbSX540XT8PX0i+zSmO1RgRvsmucixyVQxWLQGpCsuvit/zILWkORN
uOZ8uGpJ+xZxCvKKOLCXQ2tCWh6yGNIha1pTIhFKtA1yW0JcG8RH+szI1NsN
Oq3IgZZFdkkrHlxxcKRlkwRZVNjPTGfaLeQ6r3+LnrMLLYcoVqxNbBO30lqT
s/rzkkkjN1pbkhrrJNJTu9Pak5pGs4LeNE9aLhlqFyrExJvWiZC2k4lqN1/a
98i+/Lv3lgV8Wh5Jymw39ZT607qSsMseu84H7KTdRiwN5jqLrQJpPci/aRz2
sw==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"p(1\\!\\(\\*SubscriptBox[\\(s\[Sigma]\\), \\(g\\)]\\))\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-0.2941205549586153,
15.546125946802059`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E1s\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"CalcE", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[",
RowBox[{"Rest", "[", "ZRrange", "]"}], "]"}]}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<E(1\!\(\*SubscriptBox[\(s\[Sigma]\), \
\(g\)]\))/\!\(\*SuperscriptBox[\(Z\), \(2\)]\)\>\""}], "}"}]}], ",",
RowBox[{"AxesOrigin", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Blue"}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwBIQTe+yFib1JlAgAAAEEAAAACAAAAFQby0Z0XPkDCOk6LQBDhv1GOLzyF
OjtA3aDI1uQs4b95i88OMKM4QEukCmSNTOG/TpiyvfpKNkCVp0z6i2/hv0DN
c3jjKzRAHm+lSTuW4b/bz1LGekAyQGFTLiAAweG/usEQmdWDMEDDGyI2S/Dh
v39nI2YB4y1ADANJdpwk4r+BAvOA6gorQOSgL7OJXuK/xvwaDB14KEDxNL4P
zJ7iv6Qn2BYBJCZAwN1FCVjm4r91q15RnwgkQEyGUiqCNuO/PKOsw5EgIkC/
mU7qK5Hjv+0PwPj1ZiBAHJtxwej447+AUoD1wK4dQDgasmD/cOS/yMQSAKPb
GkDw7F3NKf3kvyNBMVhVTRhAIOxocg6h5b+W6CmUS/0VQDgV3eeYX+a/0mht
0pjlE0BiMvDkeDrnv4kYG4vgABJAi+b8svsx6L+rXWXTSEoQQOGYHMs8Rem/
L68439t6DUDnfHuzgHLqv7tmDCiurApAK+lomIy367//wmhv2CIIQLn+sSrk
Ee2/SMiFvtnWBUAe1lUn7H7uvwA51I/PwgNAmgUDRPn777/NeRS7ZuEBQMdo
tRoqw/C/myO/0M0tAEDUyWGjnI3xv9qMlYNRR/0/vV9Puutb8r/G31toC376
Pxd23JCoLPO/JB3+zqX49z8rJBGpX/7zv2b8mR+rsPU/QBnN1JvP9L8RxoMe
Q6DzP2YoIYHqnvW/MZi58iPC8T/Pz9cX4Wr2vzknJpmEEfA/UjsGKyMy979o
NvdDIRTtPzaC5Q5p8/e/bdB5MbpP6j8b7Tt8hq34vzKIEvW8zuc/4W/C0XBf
+b/4luNBv4rlP02MFo5ECPq//eYnFPN94z+N2YqmSaf6v++h1NEXo+E/S45v
cfY7+78b2hir2erfPxJfLPPwxfu/M0rTgkrh3D9forl6DkX8v3/I2fS5Ido/
odLmmFG5/L/ySqpgHaXXP1a5O5vmIv2/xxytsBVl1T+OpA3MHoL9v2tYJgff
W9M/qAMiy2rX/b8h+tf4QYTRPwi8b1tUI/6/e1n4Xgyzzz9zNvQCeGb+vzbV
sqPMrsw/ofv70X6h/r90kOgkCvTJP509dpkY1f6/Xi2rkcZ7xz+WrT3G9gH/
vx8eDfitP8U/ZL4bBMgo/7+qd52OBjrDP4V82Lg0Sv+/uxLdCKJlwT8OWo9X
3Gb/v/mNLaKge78/XMPFf1N//793cTALp3y8P2RWydIilP+/p3UKNarGuT/8
t4tkxqX/v1vu2gi4Urc/0rcPq620/78X0uSkhxq1Pw6v2s87wf+/5/9iQmkY
sz9oNANJyMv/v3pIo6M3R7E/PrG0op/U/79+eifKlUSvP10yw2UE3P+/u2f2
HtlKrD/KxMMLMOL/v5uZmZmZmak/jEG+9FPn/79/tStd
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"E(1\\!\\(\\*SubscriptBox[\\(s\[Sigma]\\), \
\\(g\\)]\\))/\\!\\(\\*SuperscriptBox[\\(Z\\), \\(2\\)]\\)\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-1.9939765510826843`, 0}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E1s\[Sigma]Repul", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"CalcErepul", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[",
RowBox[{"Rest", "[", "ZRrange", "]"}], "]"}]}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<E(1\!\(\*SubscriptBox[\(s\[Sigma]\), \
\(g\)]\))/\!\(\*SuperscriptBox[\(Z\), \(2\)]\)+1/(Z R)\>\""}], "}"}]}], ",",
RowBox[{"AxesOrigin", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", "0.5"}], "}"}]}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwBIQTe+yFib1JlAgAAAEEAAAACAAAAFQby0Z0XPkCfDe/HBQDgv1GOLzyF
OjtAK2tgoggA4L95i88OMKM4QDa2guYMAOC/TpiyvfpKNkA2z4lJEwDgv0DN
c3jjKzRAfeLP4hwA4L/bz1LGekAyQP+bansrAOC/usEQmdWDMECIb5ZcQgDg
v39nI2YB4y1AFXLAf2gA4L+BAvOA6gorQFwe3g2uAOC/xvwaDB14KEBTko5x
OAHgv6Qn2BYBJCZAhWGt5lgC4L91q15RnwgkQMkQ3IqvBOC/PKOsw5EgIkCR
mN9uVgngv+0PwPj1ZiBA5qKqqwQS4L+AUoD1wK4dQC4FsGUNIeC/yMQSAKPb
GkD0uHl4HTngvyNBMVhVTRhAYMpV1rFc4L+W6CmUS/0VQD0vnVVtjeC/0mht
0pjlE0Ay5ndYlMvgv4kYG4vgABJARBnunOIV4b+rXWXTSEoQQGSkstW3aeG/
L68439t6DUDbSXQWbsPhv7tmDCiurApAIppWhKwe4r//wmhv2CIIQLYCjcef
duK/SMiFvtnWBUAGEEGUFsbivwA51I/PwgNAOlbX9okH47/NeRS7ZuEBQIpf
pQcaNeO/myO/0M0tAEB+Ka9qhUjjv9qMlYNRR/0/OnHUeB8747/G31toC376
P4/MKRLHBeO/JB3+zqX49z9NkYjs3qDiv2b8mR+rsPU/XgQzf0cE4r8RxoMe
Q6DzP7VWwkFaJ+G/MZi58iPC8T/YPjWw5QDgvzknJpmEEfA/mDMUrFIO3b9o
NvdDIRTtP0SMRf6hX9m/bdB5MbpP6j+oorgp3N/UvzKIEvW8zuc/eHfLhLjv
zr/4luNBv4rlP0h71Mo1HsK//eYnFPN94z8AGCp0nbiYv++h1NEXo+E/YBbK
2266vD8b2hir2erfPxCD7ZUNPdE/M0rTgkrh3D/E62WwwL/cP3/I2fS5Ido/
OoiGrJ/s5D/ySqpgHaXXP0wWRvyHV+w/xxytsBVl1T+odwh+ilryP2tYJgff
W9M/1pWWJd4N9z8h+tf4QYTRPwIdONIaUvw/e1n4Xgyzzz9gBkBtchoBQDbV
sqPMrsw/VDXg5aliBEB0kOgkCvTJPwhq6EcQCghAXi2rkcZ7xz91X0GW6xkM
QB8eDfitP8U/+6cRG0NOEECqd52OBjrDP8pDf/elzhJAuxLdCKJlwT+CVNSt
dJQVQPmNLaKge78/zxgZELemGEB3cTALp3y8P838LyI5DRxAp3UKNarGuT+P
vDGPn9AfQFvu2gi4Urc/II8ZET/9IUAX0uSkhxq1P6olkEG4SiRA5/9iQmkY
sz+bpu6+mtYmQHpIo6M3R7E/2sTdj2anKUB+eifKlUSvP2CT4/5MxCxAu2f2
HtlKrD+2/mh/oRowQJuZmZmZmak/5hu0wIoBMkD/HPqk
"]]},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"E(1\\!\\(\\*SubscriptBox[\\(s\[Sigma]\\), \
\\(g\\)]\\))/\\!\\(\\*SuperscriptBox[\\(Z\\), \\(2\\)]\\)+1/(Z R)\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{-1, 0.5},
PlotRangeClipping->True]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"p1p\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"CalcP", "[",
RowBox[{
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[", "ZRrange", "]"}]}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+", "#"}], ")"}], "/", "2"}]}], "}"}], "&"}], "/@",
"ZRrange"}]}], "}"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<p(1\!\(\*SubscriptBox[\(p\[Sigma]\), \(u\)]\))\>\""}], "}"}]}],
",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{"Purple", ",",
RowBox[{"{",
RowBox[{"Gray", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0.5, 0, 0.5], LineBox[CompressedData["
1:eJwtxQk0GwYYwPGkRtUtpMlUgjrCUHfmar+Pdqqm3Xhbsay21/dc69RTR11F
62prpZ2xhB40xSpu88xRcaxxTMR9HyUzlves5pjwHFvf2++9//vrXQ3zDjhC
IBBc/+vdyXLrg08oZ7DBJlCSqmGBXlkWzmmOTvgVW58hwzTD2EyRiiXXHr8c
JRw4PjXBS+zaV7seTFyVtTW7+xEDURgRIza1xT72jyYrhoY4LfIWxKE1rlFn
myUO+tjUrsYZTbHErBjrGW6CHiYF6QcSxeZ4TrmXfnpYB28f2UiWKJhiWvwp
VuGONgr2tJRoMcZ4rDbUmhKghS8YE8eJBkZo7M/LWVGhYhTP3+uZvAGa+TiL
xGtkdOaWdz7SPYnCb9TtPWU1cVW1TboVqIssmo+4z42Eyd7JW22VNCyx4ZyJ
rlbD7l9JBO60FspvMy9kxKugvouV7+cXqTgr3YhwH1XEvBWTHz7eJ+OUiWBG
JugYDgWHDrGXNHDOUKU2/YOjeI96IvMNgYQdW5k3srVkkXclYNjdQw2rLXcy
Wqxl0KLKfn46ThmdNUlWUWFEbA5UMq4sV0DmWAKpsOAQDjuCwyZOymPPJj1p
PW8P3LvTWifHZJF/0SjKWrgDBIeRbFGHDF53mmuOdtoGYaxmc+ACEX3V3WyC
5zehQJ8/IDQnoHe84p5CwzpMfZeb4nV+Hx761QpT29ZAMM8KVkzchbYk6k9d
G6tgQDvoL5ZuQ/h7Ot9PtUsgcI+jw6veglvpb6Z3l5ZBXZBCP1u4AfMURynx
iyWwYjeu63f8DcSl8RwNihgcGEacVPU1KGEJ2W/JCxDM3D6rp7kKCz2SsZKn
cxA6EGvV9IkEVDmtj4tipsH6LmmrpX8ZgjKXXMgVEyDNE7vw7y9BiFPLSqnl
GOz/wdDcuCUGQkS6UejOMLwtGhmkcBcgLKu3fuhgEE5MUXmTknm4QVNlJTgM
gI3H8B2PhRlYfJV+ydZNCEndk5uNelNA6xJFJH3bC+seZV/TisfBju4akzTT
BV3aFdXkgFEY++y2z/tZr+Gar9HRedYwDMrU9J+O6wS93XGpS+ogrIdIRK35
7RDvJ/VXqhPB8p/MSr8RPjBHa7l9lX2gbOjIiupqgciHKwYKmz2QY1xe4SnX
BP7mvEeC+C4wSIyx8LzXADp0ZWkFvAZiflByzYV6KLv5JEvOsRNUTSOPh7vU
QTNpVv55QDv42ZXp6l6rgZ0H5s8L4/kgiFCL60ypgpvhloqlCS0w1+jVVnGn
HEqjFOzKBI1AzsFF3d9eAuuvKflfvBvgVFXKdbnLpUAZ+v1FhlY9KH7awN7V
Loblw6uuIdQ6MP1ZyHPQ4kKcO5f74fka8KF3l2a7F0FiHKM391wV5IZRonPc
nkFz0D/0SY9yGMu9XO7t+RgKOO+8hAQ9vntNPgcI/yuLfDBodWWR/y97n+vu
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVxX1Q03UAx/FBCQuQYAxZEA8DYhJyPE9Ej88XpYNbabI7BaVIvGOQSVZG
wlCYoYkgkgYcHHAh08MYThN64EF5UAfSkGd5ymM/BgzpkCcHAwGL4/dP77v3
vbhHTwhj9BkMxu7/Xvf/qVBYsN4vaFOup0Jt7ILdgKACgZWLLQdHVEgR81pz
g2+jZ2eY3HlZBXGoVLo95Fek+waIT7IpaF4f3f05pxJrHi45P3hRsOoavX7B
+jdw9+4T3wyjEPlikPmn8A+My7MkbgkUyhKM/MoV1Thctl96rIjCqa89jctO
18HKsPKr1mYKS1nupSXJ9VD1OxrVLFKoZT1jlsY0IrzUJqhq2wjKTxVfNgh4
AK28u5PEjsDebrPuFh7hoVSvXls2gih32RVFcjMqYvy6EqdH8O2PE85GLx/D
lFUcfDZQDX7vXalSrkTZ4e9yBnLVSD6kizKpbIeH8XBm9IIa3OU+XdC5TqSv
+aWYRY/iiwgXw+HIbtjkPW5cfDqK5ndv3bGM6YV3iNvxDw6NYU5QfsT2Rh/q
rr4fEfJ8DKktAy+ruYNw+IZ/5sv0cfgIur8XUH8jq6uqROivgc0gRzYwOYwL
fMWx35c1mL7W02klpaDbPzmR2TmB1XEee/6MGinzFhlL959Dl6cOqs8Yw+m/
Zks+ezQJ73SWtu6JBu8EXxRHaf5BfEeSV83Hk+BP1HGe8aYQx1/cw2VPwbfw
1TWzjBfYwXMpOGc+gyKLnz5qYs7AK796zqlpFhzlJ0vm1AzMFWl2e0rmwVQH
ZXM7ZiFaKbCX3dHCZ02ZfaB/Ds62a09u6BZx1TRVI1yZh2I4Ms44ZRnqafeq
TdBi8FJuWljIKuIVrY03pQsodKrvaHNnkA11aEti14ooPbLhEhg7erLbm94g
G75CaMv5+wNPN5ENV/G6Ke5EvyOTSKKrth83Z5BakclWeYURifIxYIcf1CMe
t/2Hh8SbSX540XT8PX0i+zSmO1RgRvsmucixyVQxWLQGpCsuvit/zILWkORN
uOZ8uGpJ+xZxCvKKOLCXQ2tCWh6yGNIha1pTIhFKtA1yW0JcG8RH+szI1NsN
Oq3IgZZFdkkrHlxxcKRlkwRZVNjPTGfaLeQ6r3+LnrMLLYcoVqxNbBO30lqT
s/rzkkkjN1pbkhrrJNJTu9Pak5pGs4LeNE9aLhlqFyrExJvWiZC2k4lqN1/a
98i+/Lv3lgV8Wh5Jymw39ZT607qSsMseu84H7KTdRiwN5jqLrQJpPci/aRz2
sw==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"p(1\\!\\(\\*SubscriptBox[\\(p\[Sigma]\\), \\(u\\)]\\))\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-0.573876294405314,
15.546125946802059`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output",ImageCache->GraphicsData["CompressedBitmap", \
"\<\
eJztWwtYVVUWvsIVsWx6TA9npm+iKU2bsvqayiktKyvTptFsRmvSdErCVz4i
BdHwci+XS7xf8pALorxMEBW4gCL4AG9CeBlQSZAYQoNJChVHUmwY1l5nn8e9
Z3HRbz5n5vv6Pjj33nP2v/Zaa/9rnbX32ecP8729Fiyd773wvfker6ycv9xr
4Xsfeby8bGX/KdchGs0QL41G8zcPDXzv6/8qHPr6//oWwgHOXN3359iHVmMz
e2o8PfFfZ2kfy0679h/bNRadWWPr/+Xi0FSnsbT3jR5UUyb1fqGp/LQLO7rB
ZXbVTvooEuKqgNifg//7VLTwNNsIO55EhzzLPtw0fTazQli7Racx2/rGsB83
KKxlTfqbe+osmnbl91FcmkrzfukugnWCcPbrNwp/erLzrjKlVHq53wlExQ70
zQjVy+y3ek9PoJNuww77wH0c+TBXQibxAcJbYi/st01j9mTXRinswOZDVcTe
L4hlI2i2OSqvFMsZKm+uJhbbuYtiCHED+GCsfJw8lfYOQB758MuujRqIa47S
Uf3heFkwk5AsDIyyLSFWw9Xg8cJdQsh26p6h6Exdv3OloBQ0GoGSLLbBUgZ5
fCMJU+tqjGycwQBCNDYbJrOT0nq4HEaJ0zg6Zn25f7nMMZz4LB1pHuSWsWHx
VAwNEZ1CYLirYyTZo2XBL48fGGFm6UBpZjgNk2mv3oXkjEHnl9Qd83awXw8K
5FFkc0zdwoi6Knpweruyb6ouWkEwnew6V77P4Z4lRIy7Mwy/Z3EYdc/q9+3A
9yyZu8o3TdqEzobjM7Isxyl7qyBKUo6JGuIgypo8PlldFGf3bSq0dZTTEukR
Kcl5QIHxFF0ANYDcP+y/3aLRmW2usvzNu0aHjFTev6C5bPjMZnngctHyITGz
Nuh7+zaEV3oM7gbJmjFyBws0d+NH1F6yBttoZVEkT6VSypCc6SoXJWoHxGXX
XTWk/2+Cjz5QeCwc3IWW0JHyFmivp+gGLW+kUiGoMdeFMNtmdqYrelcLes6u
z6pviB0Tu9waYXWTnC4KUzs37D/c7nr08ZMu/z1dfqZxTCdazk5ZChFJLiOz
i0rwKglPBO81BIQb1/qWnq6eiOLlxRAYU5qKmrT/h07/SZf/fV14YMinkVou
wf5m5kB0xzuFfF7uGETXHhjuCq37/8ecaThTlD4lHWrFkd3t3WNUjHvoOpy7
WdRcvEULzhmMA6+lthLFCSLm2VJtkCbgU/Pn62DyCJnq9ks6ngpzxfrMlaTE
tXhgUktLC/uFTNVCsgQeAB+AC0z+4/ZaDrIOZjVJe2hoaFZOTs7I7u5uYdkA
srDoZtbDNDlsKCu5ent715eXl3eZTKblVqtVqLdEBbNyZufAdwZ/1QHOXQO9
jm9ra+NRzCuk8W3WNux6Cs0KQV2Prq4uO1drwTdQ4QeV+pS69/b0sosv06IE
94JV4JD4goICUE24F/kcDDoIUw+PrpYu1OolmSgmSTS/p/9mVlwMQubZbEI5
7Arm2BIeTZjekNeAAibTugyTzGKusRv+vOzp2XC/FM16gRaF/B0Gpsi5CvUn
eAcyC6rzPC3jRuwa1IHBBuvASoHTj7bb2kGUxJVJztTRgjrgYRA3vaGhwUVJ
HTBQpM6ztLSfcQ3a28ecOXNGnijAz0CjSS3lLajVRFrOz9nRjdMIqAODraDO
MzRcyNqgAox6UGlpKfgHzgJvuJ9Z46edaeHOfQOipjQ1NcFZ8DGwB8Sxxr+n
xdwBRxdQxoVTuqcXTAGTxIzxFC3gl9IYgfd4HmAX5RNyljxGIr+ApjCYkBDu
6u7slmcodP/vnHXIOAp5qCUyMhKGlCchyAgaWZpTlXG3pAfnOjAIkpBSj8do
GR5SCnirbledwvBHncHcIFGDC0x7ivZs3j5/O7hbHPZxNPxeiXx52dnZMOxz
awtrFcP9MA2/T+odFOeDrSDdb2k4LldoJzcfbD664akNCsqPdQYbCrF71hRk
ys8wZPDqXb5woIrFS27A0KL09PSTURFRB1JmpEDQi/QcTcOFdTIIDTCYBzr4
CrseZc/Rh6TsElHsU1ydNCEJAoo1vo/uaJxkY7fRYGwLnx7O3YMd3UtjH0Es
0NianJzcGhEQ8XX42HAp+3vQ2MekfsEf4BegsJSqf01jMcq0MA4wnGCuOJx3
0zBcltR6Vx6sPG/UGTuDXwkGVdnFX9EwzAZDIei+SIxPPGtabao0v2yWNP0F
jR0vWcmzDHgXBgmxd9FYYa0Ohvxk1ONRc2ozaxFzJ42ZgCYm5efmXzD6GU/E
vB4jkuB2GjZRYltjTFhMt3GJMcay3CKpeZsMy6DP8SFoaoK7LBjmVR1fzRrf
SneE919WtsQXeBUcSXg+AWoGdvFmGvYiwrJzUnMuGj42FGQsyxCD6CYaNllS
sj3UEHrB+IHRuzKikl0cQcOwftL+pe6LOihAIEfB4LGLN9AwrOC04G/IioUZ
cyUl3WnYVIRVmsPNPYblhqRdgbtEMrvRsGmSkmdN/qZzQX8NEsk8lIa9hiMN
OWV2vbUewg6UBaVxpLGy0KCdDs9+1RYt/ohNwUXgqvTcxbliWSPU4cPZhxv1
GNNVReoMPFebEJTwg2GxIaIotAgVFFZF7aYZrApgJfd0VGdRVWXVhUD/wO+C
3w6GRPGjoyYqj8nUNJmJLoNiZ1Z9RT1kKXAbK5J75VLVH4mpuQzndFqe1ZPy
ffNhzC/JpSmnTGoz/VkoBZCQensMXobVBzcd/EE+guqPQ9TMfAulrarYW3FJ
/4n+dNjcMGDFRZpLb0uOWVRVWAWDD0HNqpALNGwOdgTpAowP3mPcAyacpxFz
EQGs6gjxDzkXtChoQc32mnM0Yh4iwoq3F18OWBfQEDsvFvJ1F42Yj8ZAWbWo
KrcKNIPgZynwOxr2Pp7jOdDvQOwBTaesOWu9ALW5p+tU1/mgdUHfmxaY3jhe
cvxbWqynFKZRls8sl/Rr9NVJXknMtR00zEsyYt3+lP3wzAaeTrKb+Tc0bJFE
JGgOhiw+nH34NI1YgogJrfWtPYa1htNhnmETWm2tbTRiqUSWjNyU3N4AnwBL
xscZzL+tNGwZdgR1o7E0rBQsAl//nUaskMjCs/fs+qL6r2jER4iY2ljZ2Buw
NuBk9OLocR0nOk7SCG/JmF2ZcZlXdH663K1rtjJjGu1HfxWKf6TjVEdYsbGY
F/snaPE+iIDRhtL+8ManN85oKGtooBFrEPF+TVHNZb2v/njc0rjRna2dx2iE
n0SVfakRqVcCVgdsyvPPY1Spp2HrsKMnT7WeCi3xLwFjIIvV0Qh/PoItXZA2
96e+kPpi8+fNtTRCh4hFVTurrgT4BVQneSdB6ByhEQHcYd3d1UmmJCDXhoLA
ghoagU//+mcEDc3RFl8LL+uraYQREXyVsnTztM1A+cM0woQIXVlmGTDEmuyT
DLNGqz1BgiWOH0kMTIShWF8eU36IlhyCiJnHao8l7Vq4i1coFTQiDBG8OCnI
mJkBHD9AIyIQEVS6qRS0L03zSwNy7KMRkRLRm6MCoi7rV+i9qrdUl9OIaES8
VVdVB0SH6RyUg3tpRCwioBW0zst+JxuIvodGbEBE8s74nWAHTOLAjhIaEY+I
ia2Nra0Rn0Rc0q/Uz7NtsxXTiEREvFtbUSuvqCw0YiMieNGUlL8iH7hRYM8N
M2Lzsk3ZwIu07Z9u1+TTYlMkWn8bsi7kYuDSwNe/LPxyJ43AfQbaBTVlNUBr
PgfLoxGbEcEXLjYU+BZA4sulEel4riJFnwLOjLFEWjQ5dPMMiRVdprUmuN2/
fHLvyc9oRBYiVlgtVvA9L8+yacRWRNhXYJk0ArtnN8naeH38Jf2H+oAycxkq
K+1Ec77VUq0gRG9o59msNpi3dX76wadA8S1y4YPerqhW3eHosIICajWIHYg4
VlCkDdwLsROK9cIq7u2oPF/T5gVdqnOxKhvP1JRH+rJb+vG4dXHA6lUVGRVI
dqGYveo9Z2od5UtxfN7oa2wP9QqFOMYNNEIlf3V7udR6KcRe4P4Gsc/LzCR7
JjndBaXGpCI8B0USFEtQNAnpSZXVxZJr28LXhoNrl1u3WeNpxG5ErDxUcugH
vZ8eqiWItTgasUfiHgybYhEkhobtlQIO1kB4eR1NI/BOw+rs02FrwrpMS0yz
6kvrI2nEfkQE7/5sN9RNRxIWJkBIhMsQDID3SVbwrTyUdkhRWIfS4iskE6A5
wKCwDqERhyQTYKHpu+APg6c3HGwIphGfIyK2cEshVEuV5pVmMCGIRhzmY97T
Aw7lO3WMNKJaYgnfwgCTJAONwNJLCxXSxUDfwI6QZSHwPYBGHJFYkp67MRfu
1QUZfjhBWE/DsI5kpR8kHr4Q5U8j6jgC16AgBCEU19GIo4iAsrXH4GP4JmxF
2BOnG0/72bPkmGRC9rb4bT/q1uhytvpvZSzxpcU38DHv7IKaGlgCudSHRpyQ
TOAThDm1u2tX0YgmREArIPqxuFVx4zqaO7xpxEmJ8sXpIelQU6fsCNzBjFlJ
w3CepR3X0drBC2u4ya6gETiXEycIJVumbpnaeLhxGY34GhGrKrZWAE2OJPom
3nO24+xSGtEmUX5falAqjExsYUjhEhqBk1/tE6e+OgVLoXyCsIhGCHuQ+ASh
KP3NdJgfe9mz5B/YTle2uQycak1emwzae9KSv5UI8nmyPhlQhr0xexfQiE7O
2KPN8gnCezTie0TwCQJoDwSZTyPOIiK0JLUE/Lkv1T/1lp7unndpxDnJjvoN
6zdAQQvZdC6N6JbskE8Q3qER/0QEnyDsypyTCXa8TSN6EBFfkFAAGpWm+aeB
HbNpBK6wsYlOU/T6aD5BmEUjLiNiamNNo3yC8Cd7blzBdvK1VJjezKQl/wsR
eVmhWRALeVmGLAjRGTQCj+L05mLgCjZHmO4EIdQgfAmUTxRedwITnujz2UJa
3so8YPtrTmDCemd+hikD6J6ZE5ajmeYEgxsKtNMaaxs7Qv1CLwQuC4R5w6tO
YMMFw47uO8onD5BYX0HY4Ap7cnGcHYWtBZAHIbnHFvoXAsNekvfgplGr7om3
i1yonoSdB8AAa/L6ZOBztCXRMpm2Re0VGmkxXS70jeM1x783+ZqgrALyviAX
Ouj3W9QKYna8BXuBJSz5SrKwzQJhV/eGCDkeSAVxBTrKYrSAx4Q9GAhz+vqF
2oo8O+IOB1YsnYj5JKbH4G1Yuz9r/0T7MRjcmwukv+7AbmDifsG42ggVFsTj
ADs42BEfOGrt17gH2GrBjviM05XvqGFVwHgnmDevz3X2/RpfnKS+803EwuY2
5cZiu8CRtbtBhTrSuwl9fLMy+/HSgD3ZZCi6p+GyyS3QyVP1NTF2RkgB8HCJ
iTGz5jxUpBDCaau4Q9o+zFSwatlQfNllqMOLbKCtWk5jxxcFLd0VruauE99q
4m+T3KmuoBpMcV56riZi5HN3/nap8o0QNQXtXgzh71z2f4r7xlU9KL6+KqLk
pwd44UrxCokoTTHCA/tOVO3WwftOroAdUYapcFPhqNtpvtnvuFdL645Yvu9Y
/hKMTENRIWLk3HiIKl4V6rPP9vw5sTIli5fdnBAYtRimCGP7sbyda0qEqsK1
AmakmB6ktx7kGDKfDfk3zycGIw==\
\>"]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E1p\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"CalcE", "[",
RowBox[{
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[",
RowBox[{"Rest", "[", "ZRrange", "]"}], "]"}]}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<E(1\!\(\*SubscriptBox[\(p\[Sigma]\), \
\(u\)]\))/\!\(\*SuperscriptBox[\(Z\), \(2\)]\)\>\""}], "}"}]}], ",",
RowBox[{"AxesOrigin", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Purple"}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0.5, 0, 0.5], LineBox[CompressedData["
1:eJwBIQTe+yFib1JlAgAAAEEAAAACAAAAFQby0Z0XPkAUtE2LQBDhv1GOLzyF
OjtALEfA1uQs4b95i88OMKM4QIe3pWONTOG/TpiyvfpKNkC3cJP2i2/hv0DN
c3jjKzRAv1KDLTuW4b/bz1LGekAyQN5nfnL/wOG/usEQmdWDMED3GwK5R/Dh
v39nI2YB4y1ACI9TQo0k4r+BAvOA6gorQMKYR6xQXuK/xvwaDB14KEB9b1BR
EZ7iv6Qn2BYBJCZAh6vEXTvk4r91q15RnwgkQG1dgFIKMeO/PKOsw5EgIkBV
l6iyYYTjv+0PwPj1ZiBAcUeH+6Ld47+AUoD1wK4dQBoAHHiLO+S/yMQSAKPb
GkAjAZRPI5zkvyNBMVhVTRhAaBQQ68H85L+W6CmUS/0VQPFqh/4lWuW/0mht
0pjlE0CUyN5lnLDlv4kYG4vgABJAH4A3RjH85b+rXWXTSEoQQPd3Nu3oOOa/
L68439t6DUDTUyrI/2Lmv7tmDCiurApAccZwTjF35r//wmhv2CIIQPIyDnkE
c+a/SMiFvtnWBUB2zM9xGVXmvwA51I/PwgNAW+Sr6W8d5r/NeRS7ZuEBQH1a
Tt6YzeW/myO/0M0tAECKJHk4w2jlv9qMlYNRR/0/ZiOfUJPz5L/G31toC376
P3SMOQ69c+S/JB3+zqX49z9TjRoWaO/jv2b8mR+rsPU//5cMc3ls478RxoMe
Q6DzP5ri4aPs7+K/MZi58iPC8T/X/yiAZH3ivzknJpmEEfA/x8BEvwgX4r9o
NvdDIRTtP1+25W2pveG/bdB5MbpP6j9/a9anDXHhvzKIEvW8zuc/PRGStE4w
4b/4luNBv4rlPxVCiPYo+uC//eYnFPN94z+48oJROM3gv++h1NEXo+E/zD2Y
Mh+o4L8b2hir2erfP9OgPUKcieC/M0rTgkrh3D8RUNADlHDgv3/I2fS5Ido/
Vf/B/BJc4L/ySqpgHaXXP7sN271LS+C/xxytsBVl1T/1M7D7kj3gv2tYJgff
W9M/hESSAFsy4L8h+tf4QYTRPzCJESsvKeC/e1n4Xgyzzz87rGvPryHgvzbV
sqPMrsw/x2UKnY4b4L90kOgkCvTJP8nhxIyLFuC/Xi2rkcZ7xz/UFrdOchLg
vx8eDfitP8U/ewXPKBgP4L+qd52OBjrDP+b+0DVaDOC/uxLdCKJlwT9ZcVX1
Gwrgv/mNLaKge78/XibtHkYI4L93cTALp3y8P1VniKzFBuC/p3UKNarGuT+S
LCESiwXgv1vu2gi4Urc/I7FcmokE4L8X0uSkhxq1P6upUeG2A+C/5/9iQmkY
sz/d5t1nCgPgv3pIo6M3R7E/tp79On0C4L9+eifKlUSvP6gYb6sJAuC/u2f2
HtlKrD/EX54SqwHgv5uZmZmZmak/D6lkol0B4L+EORiA
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"E(1\\!\\(\\*SubscriptBox[\\(p\[Sigma]\\), \
\\(u\\)]\\))/\\!\\(\\*SuperscriptBox[\\(Z\\), \\(2\\)]\\)\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-0.7020498783151209, 0}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E1p\[Sigma]Repul", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"CalcErepul", "[",
RowBox[{
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "2"}], ",",
RowBox[{"Reverse", "[",
RowBox[{"Rest", "[", "ZRrange", "]"}], "]"}]}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<E(1\!\(\*SubscriptBox[\(p\[Sigma]\), \
\(u\)]\))/\!\(\*SuperscriptBox[\(Z\), \(2\)]\)+1/(Z R)\>\""}], "}"}]}], ",",
RowBox[{"AxesOrigin", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Purple"}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"-", "0.6"}], ",", "0.2"}], "}"}]}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[
{RGBColor[0.5, 0, 0.5], LineBox[CompressedData["
1:eJwBIQTe+yFib1JlAgAAAEEAAAACAAAAFQby0Z0XPkDxhu7HBQDgv1GOLzyF
OjtAehFYoggA4L95i88OMKM4QHLJHeYMAOC/TpiyvfpKNkBYmNBFEwDgv0DN
c3jjKzRAHsatxhwA4L/bz1LGekAyQHywus0qAOC/usEQmdWDMEC8b3bfPgDg
v39nI2YB4y1AEf7KS1kA4L+BAvOA6gorQDoW9gZ1AOC/xvwaDB14KEDfzCCz
fQDgv6Qn2BYBJCZATC8sOzwA4L91q15RnwgkQNTPE2Zv/t+/PKOsw5EgIkBP
LHNuGPnfv+0PwPj1ZiBAdp6Ay33t37+AUoD1wK4dQCDWM/oy19+/yMQSAKPb
GkBPml/1LbDfvyNBMVhVTRhAT+X5ncpw37+W6CmUS/0VQOwJj9j0D9+/0mht
0pjlE0DI+Myyb4Pev4kYG4vgABJAsGVRYDDA3b+rXWXTSEoQQPQGme/Huty/
L68439t6DUCOQUZW2mfbv7tmDCiurApA0O68dKK82b//wmhv2CIIQN1t0iuA
r9e/SMiFvtnWBUC8DHa9hzjVvwA51I/PwgNA92kAOQFS0r/NeRS7ZuEBQOKh
I8J68c2/myO/0M0tAECC6ZJxPVjGv9qMlYNRR/0/MKmmptq2u7/G31toC376
P1DNpuovA6K/JB3+zqX49z9gm/L3hMemP2b8mR+rsPU/jFhq3dq5wD8RxoMe
Q6DzP/RdeHI4msw/MZi58iPC8T/ewaL+767UPzknJpmEEfA/Ijh7gSiM2z9o
NvdDIRTtP+uHwrBXeeE/bdB5MbpP6j9jHcU7EXrlPzKIEvW8zuc/p/C/zaTS
6T/4luNBv4rlP7O377LSju4//eYnFPN94z/Rt/gHy93xP++h1NEXo+E/yxDg
xY2z9D8b2hir2erfP2zviDdm0Pc/M0rTgkrh3D9I9eqktDz7P3/I2fS5Ido/
FBfJ8JcB/z/ySqpgHaXXP8+eOF2ClAFAxxytsBVl1T8eAR/m794DQGtYJgff
W9M/nrs3uA1mBkAh+tf4QYTRPzmKD8xrMAlAe1n4Xgyzzz+LNt96QkUMQDbV
sqPMrsw/spmbp4WsD0B0kOgkCvTJPzIoudh8txFAXi2rkcZ7xz9GONkyJcsT
QB8eDfitP8U/5bb+FpIWFkCqd52OBjrDPw5D+96nnxhAuxLdCKJlwT/afA1F
6GwbQPmNLaKge78/2uQsLIOFHkB3cTALp3y8P76iqJC0+CBAp3UKNarGuT9+
Qgjjr9wiQFvu2gi4Urc/CLvVLIzzJEAX0uSkhxq1P/Fgdk2kQidA5/9iQmkY
sz+6LpEh488pQHpIo6M3R7E/F4GEsNKhLEB+eifKlUSvPyII5fCsvy9Au2f2
HtlKrD8FSJAnt5gxQJuZmZmZmak/t9rsEvV/M0AJrAwU
"]]},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"E(1\\!\\(\\*SubscriptBox[\\(p\[Sigma]\\), \
\\(u\\)]\\))/\\!\\(\\*SuperscriptBox[\\(Z\\), \\(2\\)]\\)+1/(Z R)\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{-0.6, 0.2},
PlotRangeClipping->True]], "Output"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[TextData[{
Cell[BoxData[
FormBox[
RowBox[{
RowBox[{
RowBox[{"-", "2"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{"p", "/",
RowBox[{"(",
RowBox[{"Z", " ", "R"}], ")"}]}], ")"}], "2"]}], "\[Rule]",
RowBox[{
RowBox[{"-", "1"}], "/",
RowBox[{"(",
RowBox[{"2", " ",
SuperscriptBox["2", "2"]}], ")"}]}]}], TraditionalForm]]],
" for R\[Rule]\[Infinity]"
}], "Subsection"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"p2s\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"CalcP", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"3", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "4"}], ",",
RowBox[{"Reverse", "[", "ZRrange", "]"}]}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"#", ",",
RowBox[{
RowBox[{"(",
RowBox[{"3", "+", "#"}], ")"}], "/", "4"}]}], "}"}], "&"}], "/@",
"ZRrange"}]}], "}"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<p(2\!\(\*SubscriptBox[\(s\[Sigma]\), \(g\)]\))\>\""}], "}"}]}],
",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{"Green", ",",
RowBox[{"{",
RowBox[{"Gray", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 1, 0], LineBox[CompressedData["
1:eJwtxQk0G3YcwPG477rFsaSu0T3MsdUUb79/isei1bHqWN7MM+Jei3XOqi7t
09EZ1vYllNJU0yYRV1UjKogmrO77XKpY0KZTZ5pRW9/b573v+1pFnQ2JUcRg
MMf/6/2NVTeGK7CfIw5+LmZJYImCi5y9r3h6oYC8eySvSjzKLBw85EL3QO7D
ke6MtxYoiNr4RE50R1qO3oXnfjdDqD8tY9HhU5SGFwEjBotmB0OEWcgNTQ++
8pMnGKPWTj3aOMUFNVFP2DErDdHFWBuywqITipRT4jP29NElxc28NU0HVNxz
2rM0Tw8J98y1cRlHkMbXLSJp0CF0137KRMHWDrWPnNF/MquFzrMigm+r26Lh
pNKPg6s1kDedLSixtEZ6/PHNepoakup2yLbJlkjSK+cjvgrKC8nb7uDgkOWI
MnVJWxn1dBtg6LPmSPcL8ZJ5riKyIbiGhZ40RR0BZJ6yrgK6ufLR9cB9Y1Ri
6fXiccYBjMQlj1CXDVH+eIW9QHUffjG1KHyOMUCa+WqHKAI5sL6NGQ0g6qFy
v+KXiiwZONd5iGezdNCUEpY0yd0BHln7CIetiXbNdmkK0i046Io7O2Wtjv48
s/psirAJAT1X2qcnVBBlo3Fe1vwGMMfGfhvsUkIJ9KL2o/7r0J9pxCMvKKBv
5PgmG4kUym34Q/1OGJSuESpJTn0JM9duUIL994HJCS9ZcFwFoZgUp5Urh6er
QdPBhhKwxb0bqJHtAtN7zUVkvQzkPdphVv024AiBhKywRdAXUvA+VZuAJcQN
DTYsgCuVu2HT9Qaiu3xIoY7P4Zi9He2y/jqQZioScgrnIc5918fKSAoTywUm
XbszkDyU6dp6ag0Skxq/T6RNgdtVg+22AQk4vyB/dy9pAmQ3Fwn8gmXoNha+
nk0eg/2/7I02LyzCqQ/LS9m0Efi7emwYS1+AxgLdc9JXQ2AxY8qaXhPD+FqN
Q2DvAHxCHP2ZuDAHK7FbUpO9Z3CxZ3qLazUDOVE4B5XkXtggMiNxNZOQMrFe
64YVgeiD2nrjmHFIC+LKdHa6ITHMTk1MGoWgXWnKVyoCsJJPygiXhyF0Qq1T
g9AJ2eGyCO2mQTDSL6vaI/LBfbyR3sfpA7N6j1FKRBv8WLxiq7nVC2WtmD7i
Iy5EOLFKhNkigPK50aQTLXAYryOrhaegP8o7nWHRDMz0iiJVTwGEe51MjcM3
Ac9gXv1OTCcwI60aK79sgLe/Ot2pyuaDBoPzA5dUB+kpLlqMnDagiatccqLZ
wDiveZQp5ALVq4wd/fABkF7PqD8OaQFxFeud1JcB2JGlu/nmzSDws1ex1asB
yUHU8XjTJkilDuVQ9eiQFUCnf+bfADh/3q0wQjXkZtn/ccO3Dv65FP6T0Os2
8GJ38NNENly7cF8JfG9BOe29B5AjL4y/f50GmP+JalIf7kU94v8LS237iw==
"]]},
{GrayLevel[0.5], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxVznlM03ccxvGCirUWRwsF1gqCkBbGiCDjUvTzBTqFmjmBbKjEDlFKiRA1
SDiKUgUilzq2YcqQwKhEBwUc7OKyCKMwJoVyiMBMZRz+Wu62SEGO4fj9s3fy
5PXvYxtxOTjSkEAg+G3uvf8Pg8KC9/0IvP/CoCHqrfUQRwLFQ4KLTlcxuJHM
6sxnV4O6yCGUeh2D5ACx2PP4T8CmWzV5Z2PwZiPCL9qyFjYsTDmtIgwsescf
3qb/AnlM1pH+xxiEzQ4Tfw/+DeIGjzXn12PwKJ7kXi6rg8qQx0HucgwSrrrs
fpTSCLlpZ7wmxjBYvuNcWiKQwpUYyo5X7zb/UF8RSyOfwceKOzZnzVVQnlB0
1+hQKwgGVfwSNxXsszbWV0Ib/KMQKldCVMB1rsiTCdrhRQ09ujFBBde+xuxJ
uj9hZ9bJUlSsAo+BGvHzqucQ4ny4c7pTBYIzei65thtEGXNtLu9UYLsyqPdN
VwB7gw6ZLmq4dJq5UxnWB8Z/jb8WXFJD+97KJ7TIASgQ034tlahBwykPtyob
BG4gTSLWqSG1Y0hXZzsMJt6FRz3ZU+DG6bvFGf0bjkWfPeX4YAoYw5YVQ2ol
+DMuuIetTcHcD/0KC/EopDjE69r407A2yTLTXh8D6fKuthTlNOjvj/lKsyfg
6URiT3XEDBzMpC42yt8AuZR9pFo7A7E9Sa71n6uBcY2bO/vtLPA9lvxtzWbg
+661U8RP58CbxSxIp8yDPCuQeX77PLiK6jR2LQtA5lWtZl6ZB4oszdq/RAvp
CyXtJtg88FYL9lU8WQR66mT7RzELYG+1Li/TL0HgbWmc4/oCyJRh/N03VkCz
ZD58q1gDw7n5aUHH18DjwkR4f5AWCu2kPV3OBCTP+wqbNNdBV5JZA2/UAE36
ZJ+rndcBwbv/XnfLNiTcq26PGl2EgI6Mp0MvdiD7cT9lzNRb2GjhX365n4iE
53/2jKHooYFHdqiSkBDXzcgs9MtlOFDtpRxJNkai0AdzsU0rUHEusi+AY4K2
XIUsS0bOawIVbbkOvfzYXtGEKdpyA+5jjt+dWKOhEzqzzrJ1ArLzdT39xWeW
6JsSS0k50xB1/EEliEfo6GLELqEidvNnsHCxucoKIcfm5PDB7Wjmg2b9Is8G
1wj5iCWteTb7cYkovoIbVEy0xyWhh6yX5gb2TFwykq3SyVaJDrh70E1DrVBN
ckLoZhKKzzBBqVF2PIMxZ9RcS6kVLlBQ/TOTgoE0F1xTNNIdLEtGB3FpCHXF
JY45fYJrgU6KappWOB64H6KknO49LmIvXAYKunvAJ+PQYVxrRDPSKIosjiLS
tmVMnmOD/gUWNfnW
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"p(2\\!\\(\\*SubscriptBox[\\(s\[Sigma]\\), \\(g\\)]\\))\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-0.08744792371239561,
8.382601912613008}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E2s\[Sigma]", "=",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"CalcE", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"3", "+",
RowBox[{"Max", "[", "ZRrange", "]"}]}], ")"}], "/", "4"}], ",",
RowBox[{"Reverse", "[",
RowBox[{"Rest", "[", "ZRrange", "]"}], "]"}]}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<Z R [a.u.]\>\"", ",",
"\"\<E(2\!\(\*SubscriptBox[\(s\[Sigma]\), \
\(g\)]\))/\!\(\*SuperscriptBox[\(Z\), \(2\)]\)\>\""}], "}"}]}], ",",
RowBox[{"AxesOrigin", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Green"}]}], "]"}]}]], "Input"],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 1, 0], LineBox[CompressedData["
1:eJwBIQTe+yFib1JlAgAAAEEAAAACAAAAFQby0Z0XPkAfUgfAb93Dv1GOLzyF
OjtAEfSjkDw7xL95i88OMKM4QDALJWD8oMS/TpiyvfpKNkDZdye+Qg/Fv0DN
c3jjKzRANeO4Y7yGxb/bz1LGekAyQGi6PA4/CMa/usEQmdWDMEC4qNzs25TG
v39nI2YB4y1AXgGolu0tx7+BAvOA6gorQByChsYV1ce/xvwaDB14KECYxqvq
HYzIv6Qn2BYBJCZALceDTbVUyb91q15RnwgkQETIj+wZMMq/PKOsw5EgIkCk
zy79yx7Lv+0PwPj1ZiBAX+EYH20gzL+AUoD1wK4dQCvo21XVM82/yMQSAKPb
GkABnIffS1fOvyNBMVhVTRhAqy7zXceIz7+W6CmUS/0VQHHXLJoQY9C/0mht
0pjlE0AhNKrxmgbRv4kYG4vgABJAVO22nPqt0b+rXWXTSEoQQAz8jfwzWNK/
L68439t6DUBoCfTxVwTTv7tmDCiurApAu5VKuoKx07//wmhv2CIIQLMkxpHa
XtS/SMiFvtnWBUAua4SPjgvVvwA51I/PwgNAUT5k69W21b/NeRS7ZuEBQI1c
FLrvX9a/myO/0M0tAEAkfAgbIwbXv9qMlYNRR/0/ArSTzL+o17/G31toC376
P2A+GRcfR9i/JB3+zqX49z8GGHb+pODYv2b8mR+rsPU/sHpJp8F02b8RxoMe
Q6DzP6McNNzyAtq/MZi58iPC8T8xoY2dxYravzknJpmEEfA/Y9WqpNcL279o
NvdDIRTtP1oD9sLYhdu/bdB5MbpP6j9sYysFjPjbvzKIEvW8zuc/jlhchMhj
3L/4luNBv4rlP7g5K9F5x9y//eYnFPN94z9IGRjqnyPdv++h1NEXo+E/2Wud
s0543b8b2hir2erfP3FZ8u+sxd2/M0rTgkrh3D+a4wq88gvev3/I2fS5Ido/
5wMtn2dL3r/ySqpgHaXXP6l2akFghN6/xxytsBVl1T+OlcbiO7fev2tYJgff
W9M/XbAysGHk3r8h+tf4QYTRP2HYhRQ+DN+/e1n4Xgyzzz8qHSEiQC/fvzbV
sqPMrsw/j4JWLtdN3790kOgkCvTJP+NDbrJwaN+/Xi2rkcZ7xz8sMPh/dn/f
vx8eDfitP8U/WberT02T37+qd52OBjrDPxXqE6xTpN+/uxLdCKJlwT8sNCI0
4bLfv/mNLaKge78/8KnPLUa/3793cTALp3y8PzR2TV7Lyd+/p3UKNarGuT83
HtwbstLfv1vu2gi4Urc/Mh4MjTTa378X0uSkhxq1P7xYrQmG4N+/5/9iQmkY
sz+FQ82S0+Xfv3pIo6M3R7E/VLqYWETq379+eifKlUSvP3r6kUb67d+/u2f2
HtlKrD/NYiaQEvHfv5uZmZmZmak/ZR8tOabz37/OtR0N
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["\"Z R [a.u.]\"", TraditionalForm],
FormBox[
"\"E(2\\!\\(\\*SubscriptBox[\\(s\[Sigma]\\), \
\\(g\\)]\\))/\\!\\(\\*SuperscriptBox[\\(Z\\), \\(2\\)]\\)\"",
TraditionalForm]},
AxesOrigin->{0, 0},
PlotRange->{{0., 30.092251893604118`}, {-0.4992461737024228, 0}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"E2s\[Sigma]Repul", "=",
RowBox[{"ListPlot", "[",