Newton trajectories for the tilted Frenkel–Kontorova model
by    W.Quapp   and J.M.Bofill
in:    Molecular Physics,117:9-12,1541-1558 (2019)
DOI: 10.1080/00268976.2019.1576930
To link to this article:
https://doi.org/10.1080/00268976.2019.1576930

(*  Figure 3 *)
(*  Make the side potential *)

sidePot=Plot[1.-Cos[x] ,{x,1.*Pi-0.3,3.*Pi+0.3},Axes->{True,False},AxesStyle->Directive[20],AxesLabel->{"x","SidePot"}]

einset={4.6886493185,7.8777212958};
lll=Table[ {einset[[i]],1.-Cos[einset[[i]]]},{i,2}]

arr1=Graphics[{Red,Arrow[{{4.0,1.0},{3.5,1.5}}]}];
arr2=Graphics[{Red,Arrow[{{8.4,0.95},{7.9,0.45}}]}];
arr3=Graphics[{Blue,Arrow[{{4.2,1.2},{3.7,1.7}}]}];
arr4=Graphics[{Blue,Arrow[{{8.6,1.2},{9.1,1.7}}]}];

yende=0.51-0.0051*124;

sspring[a0_,b0_,xb_,xe_,yb_]:=Module[{a=a0,x1=xb,x2=xe,b=b0,y1=yb,n=124},h=(x2-x1)/n;
xvalues=Table[k,{k,x1,x2,h}];
yvalues=Table[y1 -b  m + a Sin[m Pi/2],{m,0,n}];
Line[Transpose@{xvalues,yvalues}]];

feder=Graphics[{sspring[0.120,0.0,4.688649318,7.8777,1.0237],
Black,PointSize[0.02],Point[{4.6886,1.0237}],
Black,PointSize[0.02],Point[{7.8777,1.0237}]},PlotRange->All] 

x1=Show[Graphics[{Text[Style["x\*SubscriptBox[\(\),\(1\)]",Large],{4.341,1.1}],
Text[Style["x\*SubscriptBox[\(\),\(2\)]",Large],{8.512,1.1}],
Text[Style["3.19",Large],{6.40,0.71}],
Text[Style["1-cos(x)",Large],{5.0,1.744}],
PointSize[Large],Black,Point[{4.688649,1.-Cos[4.688649]}],PointSize[Large],Black,Point[{7.87772,1.-Cos[7.87772]}]},
FormatType->StandardForm,PlotRange->{{1.0,3.1},{1.0,3.1}}]];

kette=Show[sidePot,feder,x1,arr1,arr2,arr3,arr4]

Export["FKmodelx1x2SP.pdf",kette,ImageResolution->Automatic,
ImageSize->Automatic]