(*Interatomicpotential*)ϕ[r_,n_,m_,A_,B_]:=-;(*Interatomicforce-*)force[r_,n_,m_,A_,B_]:=-D[ϕ[r,n,m,A,B],r];Column[{With[{n=6,m=1,A=1,B=1},r0=SolveValues[{force[r,n,m,A,B]==0,r>0},r][[1]];Plot[{ϕ[r,n,m,A,B],Evaluate[force[r,n,m,A,B]]},{r,0,5},PlotRange{-1,2},Axes->False,Frame->True,FrameLabel{"Interatomic Distance (equilibrium distance)","Interatomic Potential(arb. units)"},FrameStyleDirective[12,Black],Epilog{Dashed,Line[{{0,0},{6,0}}],Line[{{r0,-1},{r0,2}}]}]],Style["Source: Hofmann, Philip. Solid state physics: an introduction. John Wiley & Sons, 2015."]}]
A
n
r
B
m
r
Out[]=
Source: Hofmann, Philip. Solid state physics: an introduction. John Wiley & Sons, 2015. |
