如何使用gnuplot绘制像这样的Lennard Jones势函数？

Question

我开始使用gnuplot了，我试着画出一个像这样的LJ势，但只用12-6的势。另外，我想在图中显示使用的sigma和epsillon的值。

sigma = 2.6

epsillon =5

这就是我想要的样子：

​

​

这是我写的代码。

set xlabel "Distance r_0 ({\305})" 
set ylabel "Energy E (J/mol)
f(r,sigma,epsilon)=4*epsilon*((sigma/r)**12-(sigma/r)**6)  
plot[2:5.2] [-6:12] f(x,2.6,5) title '12-6 LJ'

如果有任何帮助和建议，我将不胜感激。我是一个自学成才的学生，如果这个问题微不足道，我深表歉意。

非常感谢你提前

Answer 1

这是第一次尝试：

f(r,sigma,epsilon)=4*epsilon*((sigma/r)**12-(sigma/r)**6)

sigma=2.6
epsilon=5

set arrow 1 from first 3,0 to first 3, f(3,sigma,epsilon)
set arrow 2 from first 3,0 to first 2.8, f(2.8,sigma,epsilon)
set label 1 "(r_0, {/Symbol e}_0)" at first 3,0 center offset 0,1

set xlabel "Distance r_0 ({\305})" 
set ylabel "Energy E (J/mol)

set xr [2:5.2]
set yr [-6:12]

set key reverse Left at 4,5

plot for [epsi= epsilon-1:epsilon+1:1 ] f(x,sigma,epsi) title sprintf("12-6 LJ {/Symbol s}=%.1f {/Symbol e}=%.1f",sigma,epsi)

set term pngcairo enhanced 
set out "LJ.png"
replot

结果如下：

​

​