Lennard-Jones Work: Calculating Equilibrium Position

  • Thread starter Thread starter Nacho Verdugo
  • Start date Start date
  • Tags Tags
    lennard-jones work
Click For Summary
SUMMARY

The discussion focuses on calculating the work done by the Lennard-Jones force when bringing two atoms from infinity to their equilibrium position. The force is defined as F(r) = F0 [2(σ/r)¹³ - (σ/r)⁷]. The key steps include integrating this force from A = -∞ to B = x0, where x0 represents the equilibrium position. The final expression for work is W = F0σ⁷/6 (1/x0⁶) - F0σ¹³/6(1/x0¹²). The equilibrium position is identified as the distance r0 where the force equals zero.

PREREQUISITES
  • Understanding of the Lennard-Jones potential and its mathematical representation.
  • Knowledge of integral calculus, specifically definite integrals.
  • Familiarity with concepts of force and equilibrium in physics.
  • Ability to manipulate algebraic expressions and solve equations.
NEXT STEPS
  • Learn how to derive the equilibrium position r0 from the Lennard-Jones force equation.
  • Study the implications of the Lennard-Jones potential in molecular dynamics simulations.
  • Explore numerical integration techniques for calculating work done in varying force fields.
  • Investigate the physical significance of attractive and repulsive forces in diatomic molecules.
USEFUL FOR

Students in physics or chemistry, particularly those studying intermolecular forces, molecular dynamics, or thermodynamics, will benefit from this discussion.

Nacho Verdugo
Messages
15
Reaction score
0

Homework Statement


I need to calculate the work donde by the Lennard-Jones Law, considering:

F(r)=F0 [2(σ/r)13-(σ/r)7]

when approximating two atoms from infinity to the equilibrium position between both atoms

Homework Equations



First thing I don't know how to calculate is the equilibrium position (x0) between two arbitrary atoms

The Attempt at a Solution



I just integrated this force, considering

WAB=∫F⋅dr ; r=A to r=B

As I don't know the x0, I just integrate from A to B, to later analyze the result. Assuming A=-∞ and B=x0

W= ∫ F0 [2(σ/r)13-(σ/r)7] dr ; r=A to r=B
= integrating...
= 2F0σ13 ( r-12 /-12) - F0σ7 (r-6/-6) ; still need to evaluate

Considering that A=-∞, I finally obtain:

W= F0σ7/6 (1/x06)-F0σ13/6(1/x012)

Is it right until here?
 
Physics news on Phys.org
What is the force at the equilibrium position. What does "equilibrium" mean?
What are your limits of integration? Consider particle A already at the origin. Then calculate the work done by the force on particle B as it is brought from infinity to x0.
 
kuruman said:
What is the force at the equilibrium position. What does "equilibrium" mean?
What are your limits of integration? Consider particle A already at the origin. Then calculate the work done by the force on particle B as it is brought from infinity to x0.
The statement of this problems is:

In diatomic molecules, the constituent atoms exert attractive forces between themselves at great distances and repulsive forces at short distances. For many molecules the Lennard-Jones law is a good approximation for the modulus of these forces.

$$ F(r)=F_0 \cdot [2(\frac{\sigma}{r})^{13}-(\frac{\sigma}{r})^7] $$

where r is the distance between the center of the nucleus of the atoms, ## \sigma ## a length parameter and ## F_0 ## a constant.

Determine the work to be performed by this force or by an external agent to approximate two atoms from infinity to the position of balance between the two.
- So, what I mean to equilibrium position is the position of balance between the two atoms, but there is no more information about that.

- The limits of integration, as I don't know the position of balance, are ## A=-\infty ## to ## B=x_0 ## (the so called position of balance).

- Didn't catch up the last advice
 
Last edited:
Nacho Verdugo said:
So, what I mean to equilibrium position is the position of balance between the two atoms, but there is no more information about that.
How big is the force on molecule A when it is at equilibrium or, as you say, in balance with molecule B?
 
kuruman said:
How big is the force on molecule A when it is at equilibrium or, as you say, in balance with molecule B?

I guess it is zero right? I that way it is in equilibrium.
 
Exactly. So at what inter-molecular distance r0 is the force zero?
 
kuruman said:
Exactly. So at what inter-molecular distance r0 is the force zero?

I don't know :C I took a while to answer cause of that. I think it should be at a huge distance, right?
 
It is simpler than you think. The force in general is
$$F(r)=F_0 \cdot [2(\frac{\sigma}{r})^{13}-(\frac{\sigma}{r})^7]$$
At r = r0 you have
$$0=F(r_0)=F_0 \cdot [2(\frac{\sigma}{r_0})^{13}-(\frac{\sigma}{r_0})^7]$$
Do you see what you have to do?
 
kuruman said:
It is simpler than you think. The force in general is
$$F(r)=F_0 \cdot [2(\frac{\sigma}{r})^{13}-(\frac{\sigma}{r})^7]$$
At r = r0 you have
$$0=F(r_0)=F_0 \cdot [2(\frac{\sigma}{r_0})^{13}-(\frac{\sigma}{r_0})^7]$$
Do you see what you have to do?
yes! with $$F(r_0)=0$$ I can actually find $$r_0$$ and later evaluate that value in the integral.
 
  • #10
Go for it! :smile:
 
  • Like
Likes   Reactions: Nacho Verdugo
  • #11
kuruman said:
Go for it! :smile:
thank you very much!:biggrin:
 
  • #12
Nacho Verdugo said:
yes! with $$F(r_0)=0$$ I can actually find $$r_0$$ and later evaluate that value in the integral.
 
  • #13
Nacho Verdugo said:
The statement of this problems is:

In diatomic molecules, the constituent atoms exert attractive forces between themselves at great distances and repulsive forces at short distances. For many molecules the Lennard-Jones law is a good approximation for the modulus of these forces.

$$ F(r)=F_0 \cdot [2(\frac{\sigma}{r})^{13}-(\frac{\sigma}{r})^7] $$

where r is the distance between the center of the nucleus of the atoms, ## \sigma ## a length parameter and ## F_0 ## a constant.

Determine the work to be performed by this force or by an external agent to approximate two atoms from infinity to the position of balance between the two.
- So, what I mean to equilibrium position is the position of balance between the two atoms, but there is no more information about that.

- The limits of integration, as I don't know the position of balance, are ## A=-\infty ## to ## B=x_0 ## (the so called position of balance).

- Didn't catch up the last advice
 
  • #14
Is there a question you wish to ask?
 
  • #15
kuruman said:
Is there a question you wish to ask?
No thanks! I just made a mistake editing this.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
695
  • · Replies 12 ·
Replies
12
Views
4K
Replies
15
Views
7K
Replies
1
Views
3K
Replies
1
Views
5K
  • · Replies 4 ·
Replies
4
Views
3K
Replies
3
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 9 ·
Replies
9
Views
3K