(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      7305,        197]
NotebookOptionsPosition[      5964,        150]
NotebookOutlinePosition[      6955,        186]
CellTagsIndexPosition[      6912,        183]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 StyleBox["(* SeF6 als Se+6  + 6 F- in Oh 14.02.12 *)\n\
Clear[k1,k2,k3,k4,k5,k6,sig1,sig2,sig3,sig4,sig5,c,z,R1,R2,R3,R4,R5,S2,S3,S4];\
\nc = {k1 -> 0.88, k2 -> 1.15, k3 -> 3.35, k5 -> 0.88, k6 -> 2.73,sig1 -> \
0.3, \nsig2 -> 0.30, sig3 -> 0.30, sig4 -> 0.30, sig5 -> 0.302}; z=34.0; \
z1=9.0;\nad = Sqrt[3./8.]; sq2=Sqrt[2];\n\n(* Se+6 He+Ne shell *)\nT = \
2.25*k1/R1^2+9.*k2/R2^2 /. c;\n\
Vee=3.0*sig1/R1+12.*sig2/R2+16/(R1+R2)+24*ad/(R1+R2) /. c;\n\
Vne=-3.0*z/R1-8.65*z/(R1+R2);\nS2 = R2*4^(1/3);\n\n",
  FontSize->16],
 "(* Ar shell *)\nT = T + 9.*k3/R3^2 /. c;\nVee = \
Vee+12.*sig3/R3+80./(S2+R3)+24.*ad/(S2+R3) /. c;\nVne = Vne-8.45*z/(S2+R3);\n\
S3 = R3*4^(1/3);\n",
 StyleBox["\n(* F- He+Ne shell *)\nTf = 2.25*k5/R5^2+9.*k6/R6^2 /. c;\n\
Veef=3.0*sig1/R5+12*sig2/R6+16/(R5+R6)+24*ad/(R5+R6) /. c;\n\
Vnef=-3.0*z1/R5-8.32*z1/(R5+R6) /. c;\nS5 = R6*4^(1/3);\n\nT=T+6*Tf;\n\
Vee=Vee+6*Veef;\nVne=Vne+6*Vnef;\n\n(* Se+6/6 F- *)\nVnn = (6*z*z1 + \
12*z1^2/sq2+1.5*z1^2)/(S3+S5);\nVee = Vee+(6*18*10 + 12*10*10/sq2 + \
1.5*10*10)/(S3+S5);\nVne = Vne \
-(6*z*10+6*z1*18+24*z1*10/sq2+3*z1*10)/(S3+S5);\n\nfunc = T + Vnn + Vne + \
Vee;\n\nt = FindMinimum[func, {R1,0.0349}, {R2,0.126226},{R3,0.5635},\n    \
{R5,0.14932},{R6,1.4258},{Method -> Automatic}, {MaxIterations -> 500}]\n    \
\nN[Vne /. c /. t[[2]],10]\nN[Vnn /. c /. t[[2]],10]\nN[Vee /. c /. \
t[[2]],10]\nN[-(Vee+Vne+Vnn)/T /. c /. t[[2]],10]\nN[0.529177*(S3+S5) /. c /. \
t[[2]],10]\n",
  FontSize->16],
 StyleBox["N[",
  FontSize->16],
 "ratio=-(Vne+2*Vnn)/(Vee-Vnn)",
 StyleBox[" /. c /. t[[2]],8]",
  FontSize->16],
 StyleBox["\nd=S3+S5 /.c /.t[[2]]\n",
  FontSize->16],
 "(* xy plane *)\nplot2=Graphics[{\n      \
Circle[{0,0},R1],Circle[{0,0},S2],Circle[{0,0},S3],\n      Circle[{d,0},R5],\n\
      Circle[{d,0},S5],\n      Circle[{-d,0},S5],\n      Circle[{-d,0},R5],\n \
     Circle[{0,d},R5],\n      Circle[{0,d},S5],Circle[{0,-d},R5],\n      \
Circle[{0,-d},S5]} ] /. t[[2]] /. c\n\nShow[plot2,{AspectRatio -> \
Automatic,Axes -> True,\n    GridLines -> Automatic, PlotRange \\[Rule] \
{{-6,6},{-6,6}}, \n    Frame -> True}]\n"
}], "Input",
 PageWidth->Infinity,
 InitializationCell->True,
 CellChangeTimes->{{3.5802165486441154`*^9, 3.5802165568341303`*^9}, 
   3.580216596286599*^9},
 ShowSpecialCharacters->False,
 FontSize->16],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "3002.975337134227`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.03485135552623005`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.12496966139832615`"}], ",", 
     RowBox[{"R3", "\[Rule]", "0.5651896087657294`"}], ",", 
     RowBox[{"R5", "\[Rule]", "0.1491922681167669`"}], ",", 
     RowBox[{"R6", "\[Rule]", "1.340876953886245`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166184854383`*^9}},
 FontSize->16],

Cell[BoxData[
 RowBox[{"-", "8120.027196829662`"}]], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["874.1166427486546`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["1239.9598796803384`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["1.9999999999559719`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["1.6011266507669581`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["17.416732819150617`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData["3.025692066675154`"], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.03485135552623005], 
   CircleBox[{0, 0}, 0.1983769719678126], 
   CircleBox[{0, 0}, 0.8971825795162139], 
   CircleBox[{3.025692066675154, 0}, 0.1491922681167669], 
   CircleBox[{3.025692066675154, 0}, 2.12850948715894], 
   CircleBox[{-3.025692066675154, 0}, 2.12850948715894], 
   CircleBox[{-3.025692066675154, 0}, 0.1491922681167669], 
   CircleBox[{0, 3.025692066675154}, 0.1491922681167669], 
   CircleBox[{0, 3.025692066675154}, 2.12850948715894], 
   CircleBox[{0, -3.025692066675154}, 0.1491922681167669], 
   CircleBox[{0, -3.025692066675154}, 2.12850948715894]}]], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.03485135552623005], 
   CircleBox[{0, 0}, 0.1983769719678126], 
   CircleBox[{0, 0}, 0.8971825795162139], 
   CircleBox[{3.025692066675154, 0}, 0.1491922681167669], 
   CircleBox[{3.025692066675154, 0}, 2.12850948715894], 
   CircleBox[{-3.025692066675154, 0}, 2.12850948715894], 
   CircleBox[{-3.025692066675154, 0}, 0.1491922681167669], 
   CircleBox[{0, 3.025692066675154}, 0.1491922681167669], 
   CircleBox[{0, 3.025692066675154}, 2.12850948715894], 
   CircleBox[{0, -3.025692066675154}, 0.1491922681167669], 
   CircleBox[{0, -3.025692066675154}, 2.12850948715894]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-6, 6}, {-6, 6}}]], "Output",
 CellChangeTimes->{{3.5802166184386377`*^9, 3.5802166185010386`*^9}},
 FontSize->16]
}, Open  ]]
},
AutoGeneratedPackage->Automatic,
WindowSize->{988, 740},
WindowMargins->{{Automatic, 180}, {Automatic, 36}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PageHeaders->{{
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"], Inherited, 
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"]}, {
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"], Inherited, 
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"]}},
PrintingOptions->{"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->True,
"PrintingMargins"->{{36, 36}, {57.5625, 36}}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 2326, 46, 1154, "Input",
 InitializationCell->True],
Cell[2908, 70, 532, 13, 57, "Output"],
Cell[3443, 85, 147, 3, 33, "Output"],
Cell[3593, 90, 130, 2, 33, "Output"],
Cell[3726, 94, 131, 2, 33, "Output"],
Cell[3860, 98, 131, 2, 33, "Output"],
Cell[3994, 102, 131, 2, 33, "Output"],
Cell[4128, 106, 131, 2, 33, "Output"],
Cell[4262, 110, 130, 2, 33, "Output"],
Cell[4395, 114, 719, 13, 374, "Output"],
Cell[5117, 129, 831, 18, 379, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)