(* 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[      9163,        217]
NotebookOptionsPosition[      8433,        189]
NotebookOutlinePosition[      8856,        205]
CellTagsIndexPosition[      8813,        202]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["(* Al atom Kimball, Ne-tbp core 27.12.2012 *)", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.565604726881878*^9, 3.5656047336366897`*^9}}],

Cell[CellGroupData[{

Cell["\<\
Clear[k1,k2,k3,k4, sig1,sig2,c,z,R1,R2,R3,R4,a,b,f,g,d];
c = {k1 -> 1.0, k2 -> 1.0, k3 -> 1.0, k4 -> 1.0, sig1 -> 0.28, sig2 -> 0.3, \
sig3 -> 0.3};
z=13.0; 
ad = Sqrt[3]/2;
T = 2.25*k1/R1^2+6.75*k2/R2^2+3.375*k3/R3^2+2.25*k4/R4^2 /. c;  
b = (R1+R2)/2+Sqrt[(R2+R3)^2-3*(R1+R2)^2/4];
a = (R1+R2)*Sqrt[3];d = b*Sqrt[3];
f = Sqrt[(R1+R4)^2+(R1+R2)^2];g=Sqrt[b^2+(R1+R4)^2];Vee=3.0*sig1/R1+3.0*3.0*\
sig2/R2 /. c;
Vee = Vee + \
12/(R1+R2)+6/b+12/a+4/(R1+R4)+12/f+3/d+12/(R2+R3)+6/(b+R1+R2)+6/g+0.5/(R1+R4);
Vne=-3.0*z/R1-6*z/(R1+R2)-2*z/(R1+R4)-3*z/b;

func = T + Vee + Vne /. c;

t = FindMinimum[func, {R1,0.109}, {R2,0.4198}, \
{R3,1.668},{R4,0.416},{MaxIterations -> 500},{Method-> \"Newton\"}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{
  3.565087777574009*^9, {3.5650878706125727`*^9, 3.565087885728999*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "241.5006007013421`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.1091130118246442`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.41981752153157376`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.6681677324825779`"}], ",", 
     RowBox[{"R4", "\[Rule]", "0.4159724167065372`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5656038983488226`*^9, 3.5656041694772987`*^9, 3.5656042701910753`*^9, {
   3.5656058260564437`*^9, 3.5656058910773582`*^9}, {3.5656059446322517`*^9, 
   3.565605974459504*^9}, 3.5656060262515955`*^9, 3.5656060804772906`*^9, 
   3.565606140802596*^9, 3.5656062555719976`*^9, 3.5656063876262302`*^9, {
   3.565606504938436*^9, 3.5656065166228566`*^9}, {3.5656065745769587`*^9, 
   3.5656066233114443`*^9}, 3.5656066600651083`*^9, 3.565606799232953*^9, {
   3.565606853396248*^9, 3.565606912925953*^9}, 3.565606974514861*^9, 
   3.5656071553191786`*^9, 3.565607198281654*^9, 3.577420398434599*^9, 
   3.577420466778319*^9, 3.577420527774426*^9, 3.5774205692860985`*^9, 
   3.5774206054069624`*^9, 3.5774206429874287`*^9, 3.577420741931602*^9, 
   3.5774208106341224`*^9, 3.5774208561706023`*^9, 3.577420909054695*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["-(Vne+Vee)/T /. t[[2]]", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.565087777574009*^9, 3.565087788384828*^9}}],

Cell[BoxData["1.9999999999505236`"], "Output",
 CellChangeTimes->{
  3.5656038983644223`*^9, 3.5656041694772987`*^9, 3.5656042701910753`*^9, {
   3.5656058260720434`*^9, 3.565605891092958*^9}, {3.565605944647852*^9, 
   3.565605974459504*^9}, 3.5656060262671957`*^9, 3.5656060804928904`*^9, 
   3.5656061408181963`*^9, 3.565606255587598*^9, 3.56560638764183*^9, {
   3.5656065049540358`*^9, 3.5656065166384563`*^9}, {3.5656065745769587`*^9, 
   3.5656066233114443`*^9}, 3.5656066600807085`*^9, 3.565606799232953*^9, {
   3.565606853396248*^9, 3.565606912941553*^9}, 3.565606974530461*^9, 
   3.5656071553191786`*^9, 3.565607198281654*^9, 3.577420398450199*^9, 
   3.5774204667939186`*^9, 3.5774205277900257`*^9, 3.5774205693016987`*^9, 
   3.5774206054381623`*^9, 3.5774206430030284`*^9, 3.5774207419472017`*^9, 
   3.5774208106497226`*^9, 3.5774208561862025`*^9, 3.5774209090702953`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot1=Graphics3D[{Opacity[0.5],{Lighter[Green,0.2],
Sphere[{-(R1+R2)/2,ad*(R1+R2),0},R2],
Sphere[{-(R1+R2)/2,-ad*(R1+R2),0},R2],
Sphere[{R1+R2,0,0},R2]},
{Darker[Red,1],Sphere[{0,0,0},R1]},
{Lighter[Blue,0.2],
Sphere[{0,0,R1+R4},R4],
Sphere[{0,0,-(R1+R4)},R4]},
Sphere[{-b,0,0},R3],
Sphere[{b/2,ad*b,0},R3], 
Sphere[{b/2,-ad*b,0},R3]}] /. t[[2]];

Show[plot1,{AspectRatio -> Automatic, Axes -> True}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.565087777574009*^9, 3.565087795872841*^9}, {
   3.565589646182991*^9, 3.565589647352993*^9}, {3.56558975929879*^9, 
   3.5655897682376056`*^9}, {3.5655898688577824`*^9, 
   3.5655898708545856`*^9}, {3.565589913848261*^9, 3.565589914300662*^9}, {
   3.5655899927376003`*^9, 3.565590027790861*^9}, {3.565595979449672*^9, 
   3.5655959871248856`*^9}, {3.565596027872157*^9, 3.565596042504983*^9}, {
   3.5655961175723147`*^9, 3.565596126932331*^9}, {3.565603588516678*^9, 
   3.565603634224758*^9}, 3.565603666548015*^9, {3.5656041066715884`*^9, 
   3.5656041591968803`*^9}, {3.565604260441058*^9, 3.5656042655422673`*^9}, {
   3.565605784513571*^9, 3.565605887224151*^9}, {3.565605919079407*^9, 
   3.5656059687030945`*^9}, 3.565606022008388*^9, {3.565606065423264*^9, 
   3.5656060712264743`*^9}, {3.5656061304441786`*^9, 
   3.5656061491018114`*^9}, {3.565606249316387*^9, 3.5656062500963883`*^9}, {
   3.5656063058196864`*^9, 3.5656063265989227`*^9}, {3.565606363804988*^9, 
   3.5656063705598*^9}, {3.5656064915848126`*^9, 3.5656065122548485`*^9}, {
   3.5656065602873335`*^9, 3.5656066562587023`*^9}, {3.565606782244523*^9, 
   3.5656067934921427`*^9}, {3.565606839699424*^9, 3.5656068438490314`*^9}, {
   3.5656068753454866`*^9, 3.5656069090259457`*^9}, {3.565606962222039*^9, 
   3.565606966886448*^9}, {3.565607142168355*^9, 3.565607149781169*^9}, {
   3.5656071939760466`*^9, 3.5656071947404475`*^9}, {3.57742025672395*^9, 
   3.5774202899676085`*^9}, {3.577420375892559*^9, 3.5774204532998953`*^9}, {
   3.5774205064803886`*^9, 3.577420522688817*^9}, 3.5774205639976892`*^9, 
   3.577420601397755*^9, {3.5774206374650187`*^9, 3.57742063810462*^9}, 
   3.577420731947584*^9, {3.577420768170848*^9, 3.5774207998233037`*^9}, {
   3.577420830227757*^9, 3.577420849680991*^9}, {3.577420904125087*^9, 
   3.5774209050454884`*^9}}],

Cell[BoxData[
 Graphics3DBox[
  {Opacity[0.5], 
   {RGBColor[0.2, 1., 0.2], 
    SphereBox[{-0.264465266678109, 0.45806727872373715`, 0}, 
     0.41981752153157376`], 
    SphereBox[{-0.264465266678109, -0.45806727872373715`, 0}, 
     0.41981752153157376`], 
    SphereBox[{0.528930533356218, 0, 0}, 0.41981752153157376`]}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.1091130118246442]}, 
   {RGBColor[0.2, 0.2, 1.], 
    SphereBox[{0, 0, 0.5250854285311813}, 0.4159724167065372], 
    SphereBox[{0, 0, -0.5250854285311813}, 0.4159724167065372]}, 
   SphereBox[{-2.301584995398827, 0, 0}, 1.6681677324825779`], 
   SphereBox[{1.1507924976994135`, 1.9932310749844744`, 0}, 
    1.6681677324825779`], 
   SphereBox[{1.1507924976994135`, -1.9932310749844744`, 0}, 
    1.6681677324825779`]},
  AspectRatio->Automatic,
  Axes->True,
  ImageSize->{499.6962794051012, 454.},
  ViewPoint->{2.9193980507035877`, 0.6055004361950559, -1.600151318880725},
  ViewVertical->{-0.20007553888185597`, -0.03779228891592261, \
-2.1551803583411595`}]], "Output",
 CellChangeTimes->{
  3.5656038983644223`*^9, 3.565604169492899*^9, 3.5656042702066755`*^9, {
   3.5656058263216443`*^9, 3.5656058911085577`*^9}, {3.565605944663452*^9, 
   3.5656059744751043`*^9}, 3.5656060262671957`*^9, 3.5656060805084906`*^9, 
   3.5656061408181963`*^9, 3.565606255587598*^9, 3.56560638767303*^9, {
   3.565606504969636*^9, 3.5656065166540565`*^9}, {3.5656065745925584`*^9, 
   3.565606623327044*^9}, 3.5656066600963087`*^9, 3.565606799264153*^9, {
   3.565606853411848*^9, 3.565606912941553*^9}, 3.565606974530461*^9, 
   3.5656071553347783`*^9, 3.565607198297254*^9, 3.577420398465799*^9, 
   3.5774204667939186`*^9, 3.577420527805626*^9, 3.577420569317299*^9, 
   3.5774206054381623`*^9, 3.577420643018628*^9, 3.577420741962802*^9, 
   3.5774208106497226`*^9, 3.5774208561862025`*^9, 3.5774209090858955`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.565603744173752*^9, 3.5656037681665936`*^9}, 
   3.5656038774603853`*^9},
 NumberMarks->False],

Cell["", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.565087777574009*^9, 3.5650878034232545`*^9}, {
   3.565087934182684*^9, 3.5650879365538883`*^9}, 3.565603832766307*^9}]
},
WindowSize->{1010, 740},
WindowMargins->{{Automatic, 46}, {Automatic, 73}},
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[557, 20, 155, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[737, 26, 837, 22, 252, "Input"],
Cell[1577, 50, 1225, 22, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2839, 77, 130, 2, 31, "Input"],
Cell[2972, 81, 888, 12, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3897, 98, 2305, 41, 235, "Input"],
Cell[6205, 141, 1879, 35, 469, "Output"]
}, Open  ]],
Cell[8099, 179, 143, 3, 31, "Input"],
Cell[8245, 184, 184, 3, 31, "Input"]
}
]
*)

(* End of internal cache information *)