c1={k1 -> 1.0, k2 -> 1.0, s1 -> 0.3, s2 -> 0.3};
Z = 3.0;
T=2.25*k1/P^2+1.125*k2/Q^2;
Vne=-3*Z/P-Z/(P+Q);
Vee=2/(P+Q)+3*s1/P;
func=T+Vne+Vee;
func = func /. c1;
N[t=FindMinimum[func, {P,0.55},{Q,3.1}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_1.gif]
[Graphics:Images/Li1_Ne1_gr_2.gif]
[Graphics:Images/Li1_Ne1_gr_3.gif]
[Graphics:Images/Li1_Ne1_gr_4.gif]
plot1=Graphics[{
Circle[{0,0},P],
Circle[{P+Q,0},Q],
Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_5.gif]

[Graphics:Images/Li1_Ne1_gr_6.gif]
Z = 4.0;
T=2.25*k1/P^2+2.25*k2/Q^2;
Vne=-3*Z/P-2*Z/(P+Q);
Vee=4/(P+Q)+3*s1/P+0.5/(P+Q);
func=T+Vne+Vee;
func = func /. c1;
N[t=FindMinimum[func, {P,0.40},{Q,1.9}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_7.gif]
[Graphics:Images/Li1_Ne1_gr_8.gif]
[Graphics:Images/Li1_Ne1_gr_9.gif]
[Graphics:Images/Li1_Ne1_gr_10.gif]
plot1=Graphics[{
Circle[{0,0},P],
Circle[{P+Q,0},Q],Circle[{-P-Q,0},Q],
Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_11.gif]

[Graphics:Images/Li1_Ne1_gr_12.gif]
Z=5;
T=2.25*k1/P^2+3*1.125*k2/Q^2;
Vne=-3*Z/P-3*Z/(P+Q);
Vee=6/(P+Q)+3*s1/P+3/(Sqrt[3]*(P+Q));
func=T+Vne+Vee;
func = func /. c1;
N[t=FindMinimum[func, {P,0.31},{Q,1.39}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_13.gif]
[Graphics:Images/Li1_Ne1_gr_14.gif]
[Graphics:Images/Li1_Ne1_gr_15.gif]
[Graphics:Images/Li1_Ne1_gr_16.gif]
plot1=Graphics[{
Circle[{0,0},P],
Circle[{P+Q,0},Q],Circle[{-0.5*(P+Q),0.866*(P+Q)},Q],
Circle[{-0.5*(P+Q),-0.866*(P+Q)},Q],Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_17.gif]

[Graphics:Images/Li1_Ne1_gr_18.gif]
Z=6;
T=2.25*k1/P^2+4*1.125*k2/Q^2;
Vne=-3*Z/P-4*Z/(P+Q);
Vee=8/(P+Q)+3*s1/P+6/(4*(P+Q)/Sqrt[6]);
func=T+Vne+Vee;
func = func /. c1;
N[t=FindMinimum[func, {P,0.26},{Q,1.39}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_19.gif]
[Graphics:Images/Li1_Ne1_gr_20.gif]
[Graphics:Images/Li1_Ne1_gr_21.gif]
[Graphics:Images/Li1_Ne1_gr_22.gif]
plot1=Graphics[{
Circle[{0,0},P],
Circle[{P+Q,0},Q],Circle[{-0.333*(P+Q),0.943*(P+Q)},Q],{Dashing[{0.02,0.02}],
Circle[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_23.gif]

[Graphics:Images/Li1_Ne1_gr_24.gif]
Z=7;
T=2.25*k1/P^2+5*1.125*k2/Q^2;
Vne=-3*Z/P-5*Z/(P+Q);
Vee=10/(P+Q)+3*s1/P+3*s2/Q+9/(4*(P+Q)/Sqrt[6]);
func=T+Vne+Vee;
func = func /. c1;
N[t=FindMinimum[func, {P,0.22},{Q,0.9}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_25.gif]
[Graphics:Images/Li1_Ne1_gr_26.gif]
[Graphics:Images/Li1_Ne1_gr_27.gif]
[Graphics:Images/Li1_Ne1_gr_28.gif]
plot1=Graphics[{
{GrayLevel[0.6],Disk[{P+Q,0},Q]},Circle[{P+Q,0},Q],Circle[{0,0},P],
Circle[{-0.333*(P+Q),0.943*(P+Q)},Q],{Dashing[{0.02,0.02}],
Circle[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_29.gif]

[Graphics:Images/Li1_Ne1_gr_30.gif]
Z=8;
T=2.25*k1/P^2+6*1.125*k2/Q^2;
Vne=-3*Z/P-6*Z/(P+Q);
Vee=12/(P+Q)+3*s1/P+2*3*s2/Q+13/(4*(P+Q)/Sqrt[6]);
func=T+Vne+Vee;func = func /. c1;
N[t=FindMinimum[func, {P,0.19},{Q,0.8}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_31.gif]
[Graphics:Images/Li1_Ne1_gr_32.gif]
[Graphics:Images/Li1_Ne1_gr_33.gif]
[Graphics:Images/Li1_Ne1_gr_34.gif]
plot1=Graphics[{
{GrayLevel[0.6],Disk[{P+Q,0},Q],Disk[{-0.333*(P+Q),0.943*(P+Q)},Q]},
Circle[{0,0},P],
{Dashing[{0.02,0.02}],
Circle[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Circle[{-0.333*(P+Q),0.943*(P+Q)},Q],Circle[{P+Q,0},Q],Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_35.gif]

[Graphics:Images/Li1_Ne1_gr_36.gif]
plot2=Graphics[{
{GrayLevel[0.6],Disk[{0,0.943*(P+Q)},Q],Disk[{0,0},Q]},
Circle[{0,0.943*(P+Q)},Q],Circle[{0,0},P],Circle[{-0.82*(P+Q),-0.471*(P+Q)},Q],
Circle[{0.82*(P+Q),-0.471*(P+Q)},Q],{Dashing[{0.02,0.02}],Circle[{0,0},Q]},Disk[{0,0},0.1]}] /. u;

Show[plot2,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_37.gif]

[Graphics:Images/Li1_Ne1_gr_38.gif]
Z=9;T=2.25*k1/P^2+7*1.125*k2/Q^2;
Vne=-3*Z/P-7*Z/(P+Q);
Vee=14/(P+Q)+3*s1/P+3*3*s2/Q+18/(4*(P+Q)/Sqrt[6]);
func=T+Vne+Vee;func = func /. c1;
N[t=FindMinimum[func, {P,0.16},{Q,0.7}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_39.gif]
[Graphics:Images/Li1_Ne1_gr_40.gif]
[Graphics:Images/Li1_Ne1_gr_41.gif]
[Graphics:Images/Li1_Ne1_gr_42.gif]
plot1=Graphics[{
{GrayLevel[0.6],Disk[{-0.333*(P+Q),0.943*(P+Q)},Q],
Disk[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Circle[{0,0},P],
{Dashing[{0.02,0.02}],
Circle[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Circle[{-0.333*(P+Q),0.943*(P+Q)},Q],Circle[{P+Q,0},Q],Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_43.gif]

[Graphics:Images/Li1_Ne1_gr_44.gif]
Z=10;
T=2.25*k1/P^2+4*2.25*k2/Q^2;
Vne=-3*Z/P-8*Z/(P+Q);
Vee=16/(P+Q)+3*s1/P+3*4*s2/Q+6*4/(4*(P+Q)/Sqrt[6]);
func=T+Vne+Vee;func = func /. c1;
N[t=FindMinimum[func, {P,0.14},{Q,0.62}],10]
u = t[[2]];
(* Vne *)
N[Vne /. u /. c1, 10]
(* Vee *)
N[Vee /. u /.c1, 10]
(* Virial ratio *)
N[-(Vne+Vee)/T /. u /. c1, 10]
[Graphics:Images/Li1_Ne1_gr_45.gif]
[Graphics:Images/Li1_Ne1_gr_46.gif]
[Graphics:Images/Li1_Ne1_gr_47.gif]
[Graphics:Images/Li1_Ne1_gr_48.gif]
plot1=Graphics[{{GrayLevel[0.6],
Disk[{P+Q,0},Q],Disk[{-0.333*(P+Q),0.943*(P+Q)},Q],{Dashing[{0.02,0.02}],
Disk[{-0.3333*(P+Q),-0.471*(P+Q)},Q]}},
Circle[{0,0},P],
Circle[{P+Q,0},Q],Circle[{-0.333*(P+Q),0.943*(P+Q)},Q],{Dashing[{0.02,0.02}],
Circle[{-0.3333*(P+Q),-0.471*(P+Q)},Q]},Disk[{0,0},0.1]}] /. u;

Show[plot1,{AspectRatio->Automatic, Frame -> True, Axes -> True, GridLines -> Automatic,
PlotRange -> {{-7,7},{-7,7}}}]

[Graphics:Images/Li1_Ne1_gr_49.gif]

[Graphics:Images/Li1_Ne1_gr_50.gif]

Converted by Mathematica      September 9, 2012