Plot[1-E^(-0.25 (DX)^2),{DX,0,4}]

ElementData["Carbon","CrystalStructure"]

HEnergy[r_]:=-(6/(r^11))+(5/(r^12))

Points=Values[Table[Solve[HEnergy[x]==i,x,Reals],{i,-1.3,-0.4,0.1}]]

L1=Show[Graphics[{Dashed,Line[{{First[Points[[5,1]]],HEnergy[First[Points[[5,1]]]]},{First[Points[[5,2]]],HEnergy[First[Points[[5,1]]]]}}]}],Graphics[{Disk[{Mean[{First[Points[[5,1]]],First[Points[[5,2]]]}],HEnergy[First[Points[[5,1]]]]},ImageScaled[0.01]]}]];

Clear[HL]

HL={0,0,0,0,0,0,0,0,0,0};For[i=1,i≤10,i++,
HL[[i]]={Graphics[{Dashed,Line[{{First[Points[[i,1]]],HEnergy[First[Points[[i,1]]]]},{First[Points[[i,2]]],HEnergy[First[Points[[i,1]]]]}}]}],Graphics[{Disk[{Mean[{First[Points[[i,1]]],First[Points[[i,2]]]}],HEnergy[First[Points[[i,1]]]]},ImageScaled[0.01]]}]}];

Dimensions[HL]

Show[HL]

HForce = -D[HEnergy[x],x]

NSolve[HForce==0,x,Reals]

Hr0=0.9090909090909091`;

HEMod=D[HForce,x]

HEModVal=HEMod/.x->Hr0

Hb=-HEModVal Hr0

HardPlot=Show[Plot[{HEnergy[x],0.1HForce},{x,0.8,3.0},PlotRange->{{0.8,1.5},{-1.5,1.5}},Axes->True,Frame->True,AspectRatio->1/2,PlotStyle->{,Hue[0.5]}],Plot[0.1(HEModVal x + Hb),{x,0.7,1.1},PlotStyle->{Dotted,RGBColor[{0,0,0}]}],Graphics[{Dashed,Line[{{Hr0,-1.5},{Hr0,1.5}}]}]]

ThermalHardPlot=Show[
Plot[{HEnergy[x]},{x,0.8,3.0},PlotRange->{{0.8,3.0},{-1.5,0.5}},Axes->True,Frame->True,AspectRatio->1],
Graphics[{Disk[{Hr0,HEnergy[Hr0]},ImageScaled[0.01]]}],HL
]

Hard=Plot[-(6/(r^11))+(5/(r^12)),{r,0.8,3},Frame->True,PlotRange->{{0.8,3.0},{-1.5,1.0}}]

SEnergy[r_]:=-(6/(r^3))+(100/(r^12))

SForce=-D[SEnergy[x],x]

NSolve[SForce==0,x,Reals]

Sr0=1.5946175031527288`

SEMod=D[SForce,x]

SEModVal=SEMod/.x->Sr0

Sb=-SEModVal Sr0

SEnergy[Sr0]

Clear[SPoints,SL]

SPoints=Values[Table[Solve[SEnergy[x]==i,x,Reals],{i,-1.0,-0.01,0.1}]];

Dimensions[SPoints]

First[SPoints[[11,1]]]

SL={0,0,0,0,0,0,0,0,0,0};For[i=1,i≤10,i++,
SL[[i]]={Graphics[{Dashed,Line[{{First[SPoints[[i,1]]],SEnergy[First[SPoints[[i,1]]]]},{First[SPoints[[i,2]]],SEnergy[First[SPoints[[i,1]]]]}}]}],Graphics[{Disk[{Mean[{First[SPoints[[i,1]]],First[SPoints[[i,2]]]}],SEnergy[First[SPoints[[i,1]]]]},ImageScaled[0.01]]}]}];

SoftPlot=Show[Plot[{SEnergy[x],0.1SForce},{x,0.8,3},PlotRange->{{0.8,3.0},{-1.5,1.5}},Axes->True,Frame->True,AspectRatio->1/2,PlotStyle->{Hue[0],Hue[0.1]}],Plot[0.1(SEModVal x+Sb),{x,1.3,1.8},PlotStyle->{Dotted,RGBColor[{0,0,0}]}],Graphics[{Dashed,Line[{{Sr0,-1.5},{Sr0,1.5}}]}]]

ThermalSoftPlot=Show[
Plot[{SEnergy[x]},{x,1.2,3.0},PlotRange->{{1.2,3.0},{-1.5,0.5}},Axes->False,Frame->True,AspectRatio->1,PlotStyle->Hue[0]],
Graphics[{Disk[{Sr0,SEnergy[Sr0]},ImageScaled[0.01]]}],SL
]

Show[ThermalHardPlot,ThermalSoftPlot,AspectRatio->1,FrameLabel->{"Separation (a.u.)","Energy (a.u.)"}]

Show[SoftPlot,HardPlot,PlotRange->{{0.8,2.0},{-1.5,1.5}},FrameLabel->{"Separation (a.u.)","Energy (a.u.)"}]

Soft=Plot[-(6/(r^3))+(100/(r^12)),{r,0.8,3},Frame->True,PlotRange->{{0.8,3.0},{-1.5,1.0}},PlotStyle->{Hue[0.01]}]

Show[Hard,Soft,FrameLabel->{"Separation (a.u.)","Energy (a.u.)"}]

Eng[n_,L_]:=(n^2)(0.377/(L^2))

P1=Plot[{Eng[1,1],Eng[2,1],Eng[3,1]},{x,0,1}]

P2=Plot[{Eng[1,2],Eng[2,2],Eng[3,2]},{x,0,1},PlotStyle->ColorData["Rainbow"]]

Show[P1,P2]