- #1
barnflakes
- 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: