[PoM] Lennard-Jones potential parameters

Click For Summary
SUMMARY

The discussion focuses on calculating the Lennard-Jones potential parameters, ε and σ, for the interaction between 11B boron and 14N nitrogen nuclei. The user successfully derived the moment of inertia (I) and the equilibrium separation (R0) but struggled with the relationship between the force constant (k) and the vibrational frequency (ν). The correct expression for ε was ultimately determined to be ε = [μ(2πcν)²σ²(6√2)¹⁴]/72, yielding a value of 1.5192 x 10-18 J after correcting units. This solution provides a clear path for future calculations involving molecular interactions.

PREREQUISITES
  • Understanding of Lennard-Jones potential and its parameters (ε and σ).
  • Familiarity with quantum mechanics concepts, particularly vibrational frequencies.
  • Knowledge of classical mechanics, specifically moment of inertia and force constants.
  • Ability to perform dimensional analysis in physical equations.
NEXT STEPS
  • Explore advanced topics in quantum mechanics related to molecular vibrations.
  • Learn about dimensional analysis techniques in physical chemistry.
  • Investigate the application of Lennard-Jones potential in computational chemistry simulations.
  • Study the derivation and implications of force constants in molecular dynamics.
USEFUL FOR

Students and researchers in physical chemistry, molecular physics, and computational chemistry who are involved in modeling molecular interactions and potential energy surfaces.

BRN
Messages
107
Reaction score
10
Hi guys! I need your help!

1. Homework Statement


Use the function of Lennard-Jones V (R) = ε [(σ / R)12 - (σ / R)6] as model for the adiabatic potential energy in function of the separation between the 11B boron and nitrogen nuclei 14N. You determine the parameters ε and σ to reproduce the spectroscopic values of the vibrational quantum, ν = 1514.6 cm-1, and the separation of 1,666 cm-1 between the lines of rotational BN molecule.

The Attempt at a Solution


[/B]
I have

ΔErot2/I=1.666 cm-1=3.3091*10-23 J I=ħ2/ΔErot=3.3603*10-46 Kgm2

and

R0=√(I/μ)=1.8125*10-10 m

At this point, I calculating the minimum oh the potential:

V(0)=∂V(R)/∂R =0 ⇒ R06=2σ6 ⇒ σ=6√(R06/2)=1.6147*10-10 m


For ε, I calculating k by:

k=∂2V(R)/∂R2=6σ6εR-14(26σ6-7R6)

but here I'm lost...

Who can help me? Please!
 
Physics news on Phys.org
BRN said:
For ε, I calculating k by:

k=∂2V(R)/∂R2=6σ6εR-14(26σ6-7R6)

but here I'm lost...
You need an equation to relate k to ##\nu##, then the only unknown in there will be ε.
 
  • Like
Likes   Reactions: BRN
Tanks for your answer,
the only equation i know to relate k to ν is:

ν=1/(2πc)√(k/μ)

so, i have (using R=R0):

ε=[μ(2πcν)2σ7(6√2)14]/(72σ6)

but is dimensional incorrect...
 
BRN said:
For ε, I calculating k by:

k=∂2V(R)/∂R2=6σ6εR-14(26σ6-7R6)
When you let R = Ro, you should then be able to simplify to get a nice expression relating ##k## to ##\varepsilon## and ##\sigma##. This might make it a little easier to get the correct expression for ##\varepsilon## in terms of ##\nu##.
 
Last edited:
BRN said:
Tanks for your answer,
the only equation i know to relate k to ν is:

ν=1/(2πc)√(k/μ)

so, i have (using R=R0):

ε=[μ(2πcν)2σ7(6√2)14]/(72σ6)

but is dimensional incorrect...
I think your expression for ##\varepsilon## is correct except that you appear to be off by a factor of ##\sigma##. Try to simplify the expression as much as possible.
 
  • Like
Likes   Reactions: BRN
Oh Damn!
yes, the equation is correct (apart the simplifications), but υ must be convert to m-1 and not to J!

So:

ε=[μ(2πcν)2σ2(6√2)14]/72=1.5192*10-18 J

Now it's ok!

Tanks at all!
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
Replies
5
Views
5K
  • · Replies 1 ·
Replies
1
Views
686
  • · Replies 3 ·
Replies
3
Views
819
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 14 ·
Replies
14
Views
2K
Replies
11
Views
6K