由于我对所有这些东西都不熟悉，所以我想验证我的程序是否正确。

我需要计算系统的 Lennard-Jones 势、势能和力。

程序如下：

1）读取输入数据：正方形的边长-a，粒子数-N，参数$\sigma$,$\epsilon$

2) 对于$i=1,n$在 x 和 y 轴上随机插入 n 个点。

• 2a) 对于$j=i+1,n+1$计算 x 和 y 轴上粒子之间的平方差。

• 2b) 如果$r_{ij} \leq 2^{1/6}\sigma$返回 2a) 步骤

• 2c) 两个循环结束

3) 对于$i=1,n$

• 3a) 对于$j=i+1,+1$计算伦纳德·琼斯势 $$u_{ij}=4\epsilon[(\sigma/r_{ij})^{12}-(\sigma/r_{ij})^6]$$ 在哪里$r_{ij}$在 2a) 步骤中计算
• 3b) 计算系统的势能 $$U_N=\sum_{i=1}^N\sum_{j=i+1}^N u_{ij}(r_{ij})$$
• 3c) 计算力 $$F=\sum_{j=1,j\neq i}^NF_{ij}$$ 在哪里 $$F_{ij}=(-\nabla_ju_{ij}(r_{ij}))$$
• 3d) 两个循环结束

我真的很感激检查它是否正确。我也很感激任何与这个问题相关的书籍——我目前正在使用 Allend Tildesley，但一开始对我来说相当困难。