- 156

- 4

If I have a 3 atom lattice laid out at positions ##x_1 = 1, x_2 = 2, x_3 = 3##

##x_1## and ##x_3## are stationary, but ##x_2## interacts with ##x_1## and ##x_3## via the Lennard-Jones interatomic potential.

I therefore worked out that the total acceleration for ##x_2## for ##m = \sigma = \epsilon = 1## is:

[tex]a_2 = 24(2r_{21}^{-13} - r_{21}^{-7}) + 24(-2r_{32}^{-13} + r_{32}^{-7});[/tex]

where ##r_{21} = |x_2 - x_1| , r_{32} = |x_3 - x_2|##

Is this correct?

##x_1## and ##x_3## are stationary, but ##x_2## interacts with ##x_1## and ##x_3## via the Lennard-Jones interatomic potential.

I therefore worked out that the total acceleration for ##x_2## for ##m = \sigma = \epsilon = 1## is:

[tex]a_2 = 24(2r_{21}^{-13} - r_{21}^{-7}) + 24(-2r_{32}^{-13} + r_{32}^{-7});[/tex]

where ##r_{21} = |x_2 - x_1| , r_{32} = |x_3 - x_2|##

Is this correct?

Last edited by a moderator: