Intramolecular Force and Potential Energy

SECTION 1

Historical Context & Motivation

The question of what holds atoms together within a molecule has occupied chemists and physicists for well over a century. Early in the nineteenth century, chemists recognized that certain substances resisted decomposition far more stubbornly than others, suggesting that the forces binding atoms together varied enormously in strength and character. The concept of intramolecular forces — the attractive and repulsive interactions operating within a molecule — gradually crystallized as a distinct category, separate from the weaker forces acting between molecules. The marriage of this concept with the language of potential energy provided a quantitative framework that could predict bond lengths, vibrational frequencies, and dissociation thresholds — transforming chemistry from a qualitative taxonomy into a predictive physical science.

1812

Berzelius and Electrochemical Dualism

Jöns Jacob Berzelius proposed that chemical bonds arise from the electrostatic attraction between oppositely charged atoms, laying the groundwork for understanding intramolecular forces as fundamentally electrical in nature.

1916

Lewis Electron-Pair Bond

Gilbert N. Lewis introduced the shared electron-pair model of the covalent bond, establishing that atoms could be held together by mutual sharing of electrons rather than outright transfer, expanding the taxonomy of intramolecular forces.

1927

Heitler–London Treatment of H₂

Walter Heitler and Fritz London applied quantum mechanics to the hydrogen molecule, deriving the first rigorous potential energy curve for a covalent bond and demonstrating that bond formation corresponds to an energy minimum.

1929

Morse Potential Function

Philip Morse proposed an analytically tractable potential energy function for diatomic molecules that accurately captures both the attractive and repulsive regimes, becoming a standard model for interpreting spectroscopic data.

1951

Lennard-Jones & Molecular Mechanics

Computational approaches began parameterizing intramolecular potential energy terms — stretches, bends, and torsions — into force fields, enabling the simulation of complex molecular behavior and laying the basis for modern computational chemistry.

The central question that this lesson addresses is deceptively simple: why do atoms form bonds at particular distances, and how does the potential energy of a molecular system encode all the information about bond strength, equilibrium geometry, and vibrational behavior? Answering this question requires us to dissect the competing forces at work inside a molecule and translate them into the universal language of energy.

SECTION 2

Core Principles & Definitions

Before diving into the mathematical treatment, it is essential to establish a clear conceptual vocabulary. Intramolecular forces are the net result of several simultaneous interactions: nucleus–electron attraction, electron–electron repulsion, and nucleus–nucleus repulsion. The potential energy of the system is a scalar function that encapsulates the cumulative effect of all these forces at any given internuclear separation. When this function reaches a minimum, the system achieves its most stable configuration — the equilibrium bond length.

Intramolecular Force

Any force that operates within a molecule — including covalent bonds, ionic bonds, and metallic bonds. These forces are typically 10–100 times stronger than intermolecular forces and determine molecular geometry, bond energy, and reactivity.

Potential Energy Curve

A plot of the system's potential energy V(r) as a function of internuclear distance r. The shape of this curve — its depth, width, and asymmetry — encodes the bond's strength, equilibrium length, and anharmonicity.

Bond Dissociation Energy (Dₑ)

The energy difference between the minimum of the potential energy curve and the dissociation asymptote (r → ∞). It represents the energy required to completely separate two bonded atoms from their equilibrium configuration, ignoring zero-point energy corrections.

Equilibrium Bond Length (rₑ)

The internuclear distance at which the net force between two bonded atoms is zero and the potential energy is at its minimum. At r < rₑ, repulsion dominates; at r > rₑ, attraction dominates.

Force–Energy Relationship

The force F(r) is the negative gradient of the potential energy: F = −dV/dr. This fundamental relationship connects the observable force between atoms to the underlying energy surface, and implies that the force vanishes precisely at the potential energy minimum.

✦ KEY TAKEAWAY

Think of two bonded atoms as a ball resting in a bowl. The bottom of the bowl is the equilibrium bond length — the ball naturally settles there because any displacement raises its potential energy. Pushing the ball up either side of the bowl requires energy input: squeezing atoms too close (climbing the steep inner wall) costs energy due to electron–electron and nuclear repulsion, while pulling them apart (climbing the shallow outer wall) costs energy because you work against the attractive forces. The depth of the bowl corresponds to the bond dissociation energy — deeper bowls mean stronger bonds.

SECTION 3

Visualizing the Potential Energy Curve

The potential energy curve for a diatomic molecule is arguably the single most informative diagram in all of chemical bonding. It reveals the equilibrium bond length, the bond dissociation energy, the curvature that governs vibrational frequency, and the asymmetry that produces anharmonic effects. The following diagram illustrates the canonical shape of such a curve and labels its critical features.

The curve shows three regimes: at short distances (left, repulsive region), electron cloud overlap and nuclear repulsion drive V(r) steeply upward. At large distances (right, attractive region), the atoms are essentially independent. At the equilibrium distance rₑ, attraction and repulsion balance exactly, and V(r) reaches its minimum. The well depth D_e equals the bond dissociation energy.

Notice the pronounced asymmetry of the curve. The repulsive wall at small r rises far more steeply than the attractive tail at large r — a feature that has profound spectroscopic consequences. Because the restoring force is not symmetric about r_e, molecular vibrations are anharmonic: the average bond length increases with increasing vibrational quantum number, and overtone transition frequencies are not exact integer multiples of the fundamental frequency. This anharmonicity is directly encoded in the shape of V(r).

SECTION 4

Mathematical Framework

Two analytical potential energy functions dominate introductory treatments of diatomic bonding. The simpler of the two is the harmonic oscillator approximation, which models the potential well as a perfect parabola. The more realistic model is the Morse potential, which captures the asymmetry of the true curve and allows for bond dissociation. Both connect to the force–energy relationship F = −dV/dr, which is the thread linking potential energy surfaces to observable forces.

Force–Energy Relationship

FORCE AS NEGATIVE GRADIENT

F(r) = −dV(r)/dr

F(r) = net intramolecular force along the bond axis; V(r) = potential energy as a function of internuclear distance r. At r = r_e, dV/dr = 0, so F = 0 — the hallmark of equilibrium.

Harmonic Oscillator Approximation

HARMONIC POTENTIAL

V(r) = ½ k (r − rₑ)²

k = force constant (N/m or equivalently, bond stiffness); r_e = equilibrium bond length. This parabolic approximation is valid only for small displacements from r_e. Differentiating gives F(r) = −k(r − r_e), Hooke's law.

Morse Potential

MORSE POTENTIAL

V(r) = Dₑ [1 − e^(−a(r − rₑ))]²

D_e = well depth (bond dissociation energy from the bottom of the well); a = √(k / 2D_e), a parameter controlling the width of the well; r_e = equilibrium bond length. As r → ∞, V → D_e, correctly predicting dissociation. Near r_e, the Morse potential reduces to the harmonic form.

COULOMB POTENTIAL (IONIC BOND)

V(r) = − (z₊ z₋ e²) / (4π ε₀ r) + B / rⁿ

For ionic bonds, the attractive term is the Coulomb potential between charges z₊e and z₋e. The repulsive term B/rⁿ (Born repulsion, n ≈ 9–12) accounts for electron cloud overlap at short range. ε₀ = permittivity of free space, e = elementary charge.

The force constant k is experimentally accessible through infrared and Raman spectroscopy. The vibrational frequency of a diatomic molecule is ν = (1/2π)√(k/μ), where μ is the reduced mass m₁m₂/(m₁ + m₂). A stiffer bond (larger k) means a higher vibrational frequency and a narrower, deeper potential energy well. Thus spectroscopic measurements provide a direct window into the curvature of the potential energy surface near equilibrium.

SECTION 5

Classification of Intramolecular Forces

Intramolecular forces are broadly classified into three categories — covalent, ionic, and metallic — each with a distinct potential energy profile. The differences arise from how electrons are distributed between the bonding atoms: shared (covalent), transferred (ionic), or delocalized over a lattice (metallic). The following diagram and table summarize these distinctions.

Comparison of potential energy profiles for three types of intramolecular forces. The covalent curve (solid yellow) typically has the shortest equilibrium distance. The ionic curve (dashed violet) features a deeper well for strongly electropositive/electronegative pairs and a long-range Coulombic tail. The metallic curve (dotted cyan) is characteristically shallower, reflecting the delocalized nature of metallic bonding.

Comparison of intramolecular force types and their potential energy characteristics

Property	Covalent Bond	Ionic Bond	Metallic Bond
Electron behavior	Shared between two atoms	Transferred from cation to anion	Delocalized across lattice (electron sea)
Typical Dₑ range	150–1100 kJ/mol	400–4000 kJ/mol (lattice energy)	100–400 kJ/mol
Directionality	Highly directional (σ, π overlap)	Non-directional (Coulombic)	Non-directional
V(r) dominant terms	Exchange + correlation energy	Coulomb attraction + Born repulsion	Electron kinetic energy + ion–electron Coulomb
Example	H–H in H₂ (436 kJ/mol)	Na⁺Cl⁻ (787 kJ/mol lattice)	Na metal (107 kJ/mol cohesive)

It is worth emphasizing that many real bonds exhibit intermediate character. The bond in HCl, for instance, has significant polar covalent character — neither purely covalent nor purely ionic — and its potential energy curve reflects this mixed parentage. Electronegativity differences between the bonding atoms determine where a given bond falls on the covalent-to-ionic continuum, and this positioning is directly visible in the shape and depth of V(r).

SECTION 6

Worked Example: Morse Potential for HCl

Let us apply the Morse potential to the HCl molecule. Given spectroscopic data: D_e = 431.62 kJ/mol, r_e = 1.275 × 10⁻¹⁰ m, and the force constant k = 516.3 N/m. We will (a) calculate the Morse parameter a, (b) determine V(r) at r = 1.50 × 10⁻¹⁰ m, and (c) find the force at that distance.

Morse Potential Analysis of HCl

Step 1 — Identify Given ValuesD_e = 431.62 kJ/mol = 7.168 × 10⁻¹⁹ J/molecule (dividing by Avogadro's number). r_e = 1.275 × 10⁻¹⁰ m. k = 516.3 N/m. Target distance: r = 1.50 × 10⁻¹⁰ m.

Step 2 — Calculate the Morse Parameter aThe Morse parameter is a = √(k / 2D_e). Substituting: a = √(516.3 / (2 × 7.168 × 10⁻¹⁹)) = √(3.603 × 10²⁰) = 1.898 × 10¹⁰ m⁻¹.

a = 1.898 × 10¹⁰ m⁻¹

Step 3 — Calculate V(r) at r = 1.50 × 10⁻¹⁰ mFirst compute the exponent: a(r − r_e) = 1.898 × 10¹⁰ × (1.50 − 1.275) × 10⁻¹⁰ = 1.898 × 0.225 = 0.4271. Then: V(r) = D_e[1 − e⁻⁰·⁴²⁷¹]² = 7.168 × 10⁻¹⁹ × [1 − 0.6524]² = 7.168 × 10⁻¹⁹ × (0.3476)² = 7.168 × 10⁻¹⁹ × 0.1208 = 8.66 × 10⁻²⁰ J. Converting: V = 52.2 kJ/mol.

V(1.50 Å) = 8.66 × 10⁻²⁰ J ≈ 52.2 kJ/mol above the minimum

Step 4 — Calculate the Force at r = 1.50 ÅDifferentiating the Morse potential: F(r) = −dV/dr = −2aD_e e⁻ᵃ⁽ʳ⁻ʳₑ⁾[1 − e⁻ᵃ⁽ʳ⁻ʳₑ⁾]. Since r > r_e, the bond is stretched. F = −2(1.898 × 10¹⁰)(7.168 × 10⁻¹⁹)(0.6524)(0.3476) = −2 × 1.898 × 10¹⁰ × 7.168 × 10⁻¹⁹ × 0.2268 = −6.17 × 10⁻⁹ N. The negative sign indicates the force is restoring — directed back toward r_e.

F(1.50 Å) ≈ −6.17 nN (restoring, toward equilibrium)

Step 5 — Interpret the ResultsAt 1.50 Å the HCl bond is stretched about 17.6% beyond equilibrium. The potential energy is 52.2 kJ/mol above the well minimum — about 12% of D_e. The restoring force of roughly 6 nN illustrates that even modest stretches generate significant intramolecular forces. The fact that V(r) at this displacement is well below D_e confirms the bond is far from dissociation.

SECTION 7

Harmonic vs. Morse: Strengths & Limitations

Choosing between the harmonic oscillator and the Morse potential involves a trade-off between mathematical simplicity and physical realism. Each model has a well-defined domain of validity, and recognizing the boundaries of each is essential for selecting the appropriate tool for a given problem.

Comparison of the harmonic oscillator and Morse potential models for diatomic molecules

Feature	Harmonic Oscillator	Morse Potential
Functional form	V = ½k(r − rₑ)²	V = Dₑ[1 − e⁻ᵃ⁽ʳ⁻ʳₑ⁾]²
Symmetry about rₑ	Perfectly symmetric (parabola)	Asymmetric — steep repulsive wall, shallow attractive tail
Bond dissociation	Cannot model (V → ∞ as r → ∞)	Correctly predicts V → Dₑ as r → ∞
Number of parameters	2 (k, rₑ)	3 (Dₑ, a, rₑ)
Vibrational energy levels	Equally spaced: Eᵥ = ℏω(v + ½)	Converging: Eᵥ = ℏω(v + ½) − ℏωχₑ(v + ½)²
Best use case	Low-v fundamentals, quick estimates, pedagogical clarity	Overtone analysis, dissociation thresholds, high-v spectroscopy
Key limitation	Fails badly for \|r − rₑ\| > ~10% of rₑ	Still approximate; ignores coupling between vibration & rotation

✦ KEY TAKEAWAY

The harmonic oscillator is like modeling a mountain valley as a perfect V-shaped trench — adequate for describing motion at the bottom, but it predicts infinitely high walls that prevent escape. The Morse potential, by contrast, captures the fact that the valley eventually flattens out: if you give a molecule enough energy, the atoms can fly apart. In research-grade computational chemistry, neither model is typically used directly; instead, ab initio methods compute V(r) point by point on a grid, and analytical functions like Morse are fitted to these data for convenience.

SECTION 8

Connection to Potential Energy Surfaces & Advanced Theory

The one-dimensional V(r) curve for a diatomic is a special case of the far richer concept of a potential energy surface (PES). For a polyatomic molecule with N atoms, the potential energy depends on 3N − 6 internal coordinates (3N − 5 for a linear molecule), generating a high-dimensional hypersurface whose minima, saddle points, and valleys dictate molecular structure, conformational dynamics, and reaction pathways. Understanding the diatomic V(r) curve is thus the essential first step toward navigating the full PES.

From diatomic V(r) to the full potential energy surface

Concept	Diatomic V(r)	Polyatomic PES
Dimensionality	1D — single internuclear coordinate r	(3N − 6)D — bond lengths, angles, dihedrals
Equilibrium	Single minimum at rₑ	Multiple minima (conformers, isomers)
Transition states	Not applicable (no reaction path)	Saddle points connecting minima; govern reaction rates via transition state theory
Computation	Exact QM solution for H₂⁺; analytic fits (Morse)	Ab initio (DFT, CCSD(T)); molecular mechanics force fields
Key insight	Curvature → vibrational frequency	Gradient → forces in molecular dynamics; Hessian → normal modes

In courses on physical chemistry and quantum mechanics, you will encounter the Born–Oppenheimer approximation, which underpins the entire concept of a potential energy surface. This approximation separates electronic and nuclear motion, allowing one to compute V(r) by solving the electronic Schrödinger equation at fixed nuclear positions. Without it, the notion of a smooth, well-defined V(r) curve would not exist — making it one of the most consequential ideas in theoretical chemistry. Furthermore, the concepts of molecular orbital theory and valence bond theory provide complementary frameworks for understanding why particular combinations of atomic orbitals yield bonding (lower energy) or antibonding (higher energy) molecular states, directly shaping the V(r) curve you have studied in this lesson.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

Explain why the potential energy curve for a diatomic molecule is asymmetric — steep on the left (short r) and shallow on the right (long r). What distinct physical interactions are responsible for each side of the curve, and how does this asymmetry manifest in the molecule's spectroscopic behavior?

PROBLEM 2 — BASIC CALCULATION

The force constant for the C=O bond in carbon monoxide is k = 1902 N/m and the equilibrium bond length is r_e = 1.128 × 10⁻¹⁰ m. Using the harmonic approximation, calculate the potential energy when the bond is compressed to r = 1.08 × 10⁻¹⁰ m.

PROBLEM 3 — INTERMEDIATE

For HF, the Morse parameters are D_e = 570 kJ/mol, r_e = 0.917 × 10⁻¹⁰ m, and a = 2.218 × 10¹⁰ m⁻¹. (a) Determine the force constant k from these Morse parameters. (b) Compare with the experimentally reported k = 966 N/m and comment on the agreement.

PROBLEM 4 — APPLIED

A researcher studying a novel diatomic molecule AB observes a fundamental vibrational frequency of ν̃ = 2145 cm⁻¹ and determines the bond dissociation energy to be D₀ = 1070 kJ/mol. The reduced mass is μ = 1.139 × 10⁻²⁶ kg. Using these data and the Morse model, estimate the anharmonicity constant χ_e. (Hint: for the Morse oscillator, χ_e = ℏω / 4D_e, and D₀ ≈ D_e − ½ℏω + ¼ℏωχ_e.)

PROBLEM 5 — CRITICAL THINKING

Consider two hypothetical diatomic molecules X₂ and Y₂. Both have the same equilibrium bond length r_e = 1.50 Å, but X₂ has D_e = 200 kJ/mol while Y₂ has D_e = 800 kJ/mol. Both have the same reduced mass. Compare and contrast the potential energy curves of X₂ and Y₂ in terms of: (a) well depth, (b) force constant, (c) Morse parameter a, (d) vibrational frequency, and (e) anharmonicity constant χ_e. Discuss which molecule would be easier to dissociate photolytically and why.

SUMMARY

Summary

Intramolecular forces — covalent, ionic, and metallic — are the strong interactions that hold atoms together within molecules and extended solids. These forces arise from the balance between attractive interactions (nucleus–electron Coulomb attraction, electron sharing, or charge transfer) and repulsive interactions (electron–electron and nucleus–nucleus repulsion). The potential energy curve V(r) encodes this balance as a function of internuclear distance, reaching a minimum at the equilibrium bond length rₑ where the net force is zero. The depth of this minimum defines the bond dissociation energy Dₑ.

Quantitatively, the harmonic oscillator approximation V = ½k(r − rₑ)² captures behavior near equilibrium but fails at large displacements. The Morse potential extends the model to include anharmonicity and correctly predicts bond dissociation. The fundamental relation F = −dV/dr connects force and energy, and the curvature of V(r) at the minimum gives the force constant k, which determines vibrational frequency. These one-dimensional concepts generalize to the multidimensional potential energy surface that governs all of molecular structure and reactivity.