Physics · Book 3 · Bachelor Year 1

University Physics — Year 1

University Physics — Year 1 · Bachelor Year 1

13Work, Energy, and Potential Energy

A cyclist stops pedaling at the top of a hill and arrives at the bottom at a speed that depends on the drop, not on the shape of the road. A bungee jumper falls sixty meters and stops, for an instant, exactly where the cord’s pull has eaten all the speed the fall had given. A hydrogen atom bound to a chlorine atom vibrates in a well whose depth is the energy it costs to tear the molecule apart. In each case a single number is exchanged between forms — motion, height, stretch, bond — and summed, it does not change. This chapter builds that bookkeeping from Newton’s second law: work and power, the kinetic energy theorem, potential energies for the forces that allow one, and the conservation of mechanical energy that turns many dynamical problems into one-line algebra and explains, through the shape of a potential well, which positions are stable.

The lecturer’s pendulum: released from the nose, the bowling ball returns to the nose and no further — mechanical energy is conserved, and never exceeded.
The lecturer’s pendulum: released from the nose, the bowling ball returns to the nose and no further — mechanical energy is conserved, and never exceeded.

13.1 Power and work

Definition 13.1 (Power and work of a force)

The power of a force F\vect F acting on a point moving at velocity v\vect v is

P=Fv(watts),\mathcal P = \vect F\cdot\vect v \qquad (\text{watts}),

and its work between times t1t_1 and t2t_2 (positions AA and BB) is the accumulated power,

WAB=t1t2Fv ⁣dt=ABF ⁣d(joules),W_{A\to B} = \int_{t_1}^{t_2}\vect F\cdot\vect v\,\dd t = \int_A^B \vect F\cdot\dd\vect\ell \qquad (\text{joules}),

the sum of the elementary works δW=F ⁣d\delta W = \vect F\cdot\dd\vect\ell along the elementary displacements  ⁣d=v ⁣dt\dd\vect\ell = \vect v\,\dd t of the trajectory. Work is positive if the force pushes along the motion (motive), negative if against it (resistive), zero if perpendicular.

Proposition 13.2 (Works one meets)

  • Constant force: WAB=FABW_{A\to B} = \vect F\cdot\vect{AB}, whatever the path. Weight: W=mg(zAzB)W = mg(z_A - z_B), zz counted upward.
  • Spring of stiffness kk, from elongation xAx_A to xBx_B: W=12k(xA2xB2)W = \tfrac12 k(x_A^2 - x_B^2).
  • Normal reaction of a fixed surface, string tension on a body moving perpendicular to it: W=0W = 0 (force \perp displacement).
  • Kinetic friction and drag: W<0W < 0 always (the force opposes the relative motion), and depends on the path length.

Proof. Constant F\vect F: F ⁣d=F ⁣d=FAB\int\vect F\cdot\dd\vect\ell = \vect F\cdot\int\dd\vect\ell = \vect F\cdot\vect{AB}; with F=mgez\vect F = -mg\vect e_z, FAB=mg(zBzA)\vect F\cdot\vect{AB} = -mg(z_B - z_A). Spring: F ⁣d=kx ⁣dx\vect F\cdot\dd\vect\ell = -kx\,\dd x, and xAxBkx ⁣dx=12k(xA2xB2)\int_{x_A}^{x_B}-kx\,\dd x = \tfrac12 k(x_A^2 - x_B^2). Friction: T=μNv/v\vect T = -\mu N\vect v/v, so Tv=μNv<0\vect T\cdot\vect v = -\mu Nv < 0.

The elementary work of a force along a small displacement of its point of application: only the component of F along the motion works. Summed along the path from A to B, it gives W_A B.
The elementary work of a force along a small displacement of its point of application: only the component of F\vect F along the motion works. Summed along the path from AA to BB, it gives WABW_{A\to B}.

13.2 The kinetic energy theorem

Theorem 13.3 (Kinetic energy and power theorems)

In an inertial frame, the kinetic energy Ek=12mv2E_k = \tfrac12 mv^2 of a point particle obeys

 ⁣dEk ⁣dt=Pi=(Fi)v,Ek(B)Ek(A)=iWi,AB:\frac{\dd E_k}{\dd t} = \sum\mathcal P_i = \Big(\sum\vect F_i\Big)\cdot\vect v, \qquad E_k(B) - E_k(A) = \sum_i W_{i,A\to B}:

the change of kinetic energy between two positions equals the total work of all the forces applied along the way.

Proof. Dot Newton’s second law with v\vect v: mav= ⁣d ⁣dt(12mvv)= ⁣dEk/ ⁣dtm\vect a\cdot\vect v = \tfrac{\dd}{\dd t}(\tfrac12 m\vect v\cdot\vect v) = \dd E_k/\dd t, while (Fi)v=Pi(\sum\vect F_i)\cdot\vect v = \sum\mathcal P_i. Integrate over time.

Example 13.4 (Braking distance)

A 1000kg1000\,\mathrm{kg} car at 100km/h100\,\mathrm{km}/\mathrm{h} (Ek=386kJE_k = 386\,\mathrm{kJ}) braked by a constant 7.0kN7.0\,\mathrm{kN} force stops after d=Ek/F=55md = E_k/F = 55\,\mathrm{m}: the distance grows as the square of the speed. Since the maximal braking force is μsmg\mu_smg (Chapter 12), the shortest stopping distance is v2/2μsgv^2/2\mu_sg whatever the car’s mass — and it doubles on a wet road where μs\mu_s halves.

13.3 Conservative forces and potential energy

Definition 13.5 (Conservative force; potential energy)

A force is conservative if its work between two points does not depend on the path followed. Equivalently, there exists a function of position, the potential energy EpE_p, such that

WAB=Ep(A)Ep(B)=ΔEp,δW= ⁣dEp,W_{A\to B} = E_p(A) - E_p(B) = -\Delta E_p , \qquad \delta W = -\dd E_p ,

EpE_p being defined up to an additive constant (a choice of origin).

Proposition 13.6 (Force from potential energy)

A conservative force derives from its potential energy:

F=gradEp=EpxexEpyeyEpzez,\vect F = -\overrightarrow{\operatorname{grad}}\,E_p = -\frac{\partial E_p}{\partial x}\,\vect e_x - \frac{\partial E_p}{\partial y}\,\vect e_y - \frac{\partial E_p}{\partial z}\,\vect e_z ,

and in one dimension Fx= ⁣dEp/ ⁣dxF_x = -\dd E_p/\dd x: the force points toward decreasing potential energy, “downhill” on the graph of EpE_p.

Proof. δW=Fx ⁣dx+Fy ⁣dy+Fz ⁣dz= ⁣dEp=(xEp ⁣dx+yEp ⁣dy+zEp ⁣dz)\delta W = F_x\dd x + F_y\dd y + F_z\dd z = -\dd E_p = -(\partial_xE_p\,\dd x + \partial_yE_p\,\dd y + \partial_zE_p\,\dd z) for every displacement: identify the coefficients (the differential of a function of several variables, from the mathematics course).

Proposition 13.7 (The usual potential energies)

  • Weight: Ep=mgz+constE_p = mgz + \text{const} (zz upward).
  • Spring: Ep=12k(0)2+constE_p = \tfrac12 k(\ell - \ell_0)^2 + \text{const}.
  • Newtonian gravitation of a mass MM: Ep=GMm/r+constE_p = -GMm/r + \text{const}, the constant usually chosen so that Ep0E_p \to 0 at infinity.
  • Electrostatic force on a charge qq in a potential VV: Ep=qVE_p = qV (Chapter 27).

Kinetic friction and fluid drag are not conservative (their work depends on the path, and is always negative).

Proof. Each is read off the works of Proposition 13.2: W=ΔEpW = -\Delta E_p. For gravitation, F=GMm/r2er\vect F = -GMm/r^2\,\vect e_r and δW=GMm ⁣dr/r2= ⁣d(GMm/r)\delta W = -GMm\,\dd r/r^2 = -\dd(-GMm/r).

13.4 Mechanical energy

Theorem 13.8 (Mechanical energy)

Let Em=Ek+EpE_m = E_k + E_p be the mechanical energy, where EpE_p is the sum of the potential energies of the conservative forces. Then

ΔEm=Wnc,\Delta E_m = W_{\text{nc}} ,

the total work of the non-conservative forces. If those do no work (no friction, or constraint forces perpendicular to the motion), EmE_m is conserved: the motion is conservative.

Proof. ΔEk=Wc+Wnc=ΔEp+Wnc\Delta E_k = W_{\text{c}} + W_{\text{nc}} = -\Delta E_p + W_{\text{nc}}.

Method 13.9 (Using energy)

When a problem asks for a speed at a position (not a time, not a force), try energy first:

  1. list the forces; check which work and which are conservative;
  2. write EmE_m at the two positions, with a clear origin for EpE_p;
  3. equate them (adding WncW_{\text{nc}} if friction acts) and solve for the unknown speed or position.

Energy never gives the time or the constraint forces; for those, go back to Newton’s law — often with the speed that energy just provided.

Example 13.10 (The loop’s speed law)

On the frictionless loop of Problem 11.1 the track’s reaction is normal to the motion and works not: EmE_m is conserved. With z=R(1cosθ)z = R(1 - \cos\theta) above the bottom, 12mv2+mgR(1cosθ)=12mv02\tfrac12 mv^2 + mgR(1 - \cos\theta) = \tfrac12 mv_0^2: the speed law v2=v022gR(1cosθ)v^2 = v_0^2 - 2gR(1 - \cos\theta) used there, now derived in one line. The minimum entry speed 5gR\sqrt{5gR} then follows from the top condition v2gRv^2 \geq gR of Example 12.16.

Example 13.11 (Energy lost to friction)

A block slides down a slope of angle α\alpha and height hh with kinetic friction μd\mu_d: Wnc=μdmgcosα×h/sinα=μdmghcotαW_{\text{nc}} = -\mu_dmg\cos\alpha \times h/\sin\alpha = -\mu_dmgh\cot\alpha, so 12mv2=mgh(1μdcotα)\tfrac12 mv^2 = mgh(1 - \mu_d\cot\alpha). For h=2.0mh = 2.0\,\mathrm{m}, α=30\alpha = 30^\circ, μd=0.20\mu_d = 0.20: v=5.1m/sv = 5.1\,\mathrm{m}/\mathrm{s} instead of 6.3m/s6.3\,\mathrm{m}/\mathrm{s}; a third of the potential energy became heat.

13.5 Equilibrium, stability, and phase portraits

Proposition 13.12 (Equilibrium and stability in one dimension)

A particle moving on an axis under a conservative force of potential energy Ep(x)E_p(x):

  • is in equilibrium at x0x_0 iff F(x0)=Ep(x0)=0F(x_0) = -E_p'(x_0) = 0: the extrema of EpE_p;
  • the equilibrium is stable if EpE_p has a local minimum there (Ep(x0)>0E_p''(x_0) > 0): a small displacement produces a restoring force; unstable at a maximum;
  • near a stable equilibrium, small motions are harmonic, with

    ω02=Ep(x0)m.\omega_0^2 = \frac{E_p''(x_0)}{m} .

Proof. Taylor’s formula (mathematics course): Ep(x)Ep(x0)+12Ep(x0)(xx0)2E_p(x) \approx E_p(x_0) + \tfrac12 E_p''(x_0)(x - x_0)^2, so FEp(x0)(xx0)F \approx -E_p''(x_0)(x - x_0) and mx¨=Ep(x0)(xx0)m\ddot x = -E_p''(x_0)(x - x_0) — a spring of stiffness Ep(x0)E_p''(x_0). For Ep<0E_p'' < 0 the “stiffness” is negative: exponential departure.

A one-dimensional potential landscape: minima are stable equilibria, the maximum an unstable one. A particle of energy E_1 below the barrier is trapped in one well between two turning points (where E_p = E_1, arrows) and oscillates; with E_2 above the barrier it passes freely from one well to the other.
A one-dimensional potential landscape: minima are stable equilibria, the maximum an unstable one. A particle of energy E1E_1 below the barrier is trapped in one well between two turning points (where Ep=E1E_p = E_1, arrows) and oscillates; with E2E_2 above the barrier it passes freely from one well to the other.

Definition 13.13 (Turning points; bound and free motion)

For a conservative one-dimensional motion of energy EmE_m, the particle can only be where Ep(x)EmE_p(x) \leq E_m (since Ek0E_k \geq 0); the points where Ep=EmE_p = E_m are turning points, where the velocity vanishes and reverses. Trapped between two turning points in a potential well, the motion is bound and periodic; able to reach infinity, it is free; a potential barrier higher than EmE_m cannot be crossed.

Definition 13.14 (Phase portrait)

The phase portrait of a one-dimensional motion is the family of curves traced by the point (x,x˙)(x, \dot x) in the phase plane. For a conservative system each curve is a level line 12mx˙2+Ep(x)=Em\tfrac12 m\dot x^2 + E_p(x) = E_m: closed curves around stable equilibria (oscillations), open curves for free motion, and, through unstable equilibria, the separatrices that divide the two. Curves are traversed clockwise (x˙>0\dot x > 0 means xx increasing) and never cross.

Phase portrait of the simple pendulum, 1/2 2 + _02(1 - ) = const. Closed curves (blue): swings about the stable equilibrium = 0; open curves (orange): full rotations; between them the separatrix (red) through the unstable equilibria = ±π, the motion that just reaches the top.
Phase portrait of the simple pendulum, 12θ˙2+ω02(1cosθ)=const\tfrac12\dot\theta^2 + \omega_0^2(1 - \cos\theta) = \text{const}. Closed curves (blue): swings about the stable equilibrium θ=0\theta = 0; open curves (orange): full rotations; between them the separatrix (red) through the unstable equilibria θ=±π\theta = \pm\pi, the motion that just reaches the top.

Example 13.15 (The pendulum’s landscape)

Ep=mgcosθE_p = -mg\ell\cos\theta: minimum at θ=0\theta = 0, Ep=mgE_p'' = mg\ell, so ω02=g/(m2)=g/\omega_0^2 = g\ell/(m\ell^2) = g/\ell — the small-oscillation result of Example 12.12 read off the well’s curvature; maximum at θ=π\theta = \pi, the inverted pendulum, unstable. Energy Em<mgE_m < mg\ell: oscillation between ±θmax\pm\theta_{\max}; Em>mgE_m > mg\ell: it goes over the top and rotates; Em=mgE_m = mg\ell: the separatrix, an infinitely slow approach to the top.

13.6 Exercises

Exercise 13.1

A 20kg20\,\mathrm{kg} crate is lifted 3.0m3.0\,\mathrm{m} at constant speed in 5.0s5.0\,\mathrm{s}. Work of the lifting force, of the weight, and power supplied. The same crate pushed 3.0m3.0\,\mathrm{m} along a floor with μd=0.30\mu_d = 0.30: work of friction.

Solution

Solution of Exercise 13.1.

Lifting force +mgh=20×9.81×3.0=589J+mgh = 20 \times 9.81 \times 3.0 = 589\,\mathrm{J}; weight 589J-589\,\mathrm{J}; power 589/5.0=118W589/5.0 = 118\,\mathrm{W}. Friction: μdmgd=0.30×196×3.0=177J-\mu_dmgd = -0.30 \times 196 \times 3.0 = -177\,\mathrm{J}.

Exercise 13.2

Kinetic energy of a 1000kg1000\,\mathrm{kg} car at 100km/h100\,\mathrm{km}/\mathrm{h}; stopping distance with a constant 7.0kN7.0\,\mathrm{kN} braking force; what happens to the distance if the speed is doubled?

Solution

Solution of Exercise 13.2.

Ek=12×1000×27.82=386kJE_k = \tfrac12 \times 1000 \times 27.8^2 = 386\,\mathrm{kJ}; d=Ek/F=55md = E_k/F = 55\,\mathrm{m}; doubling vv quadruples EkE_k and dd: 220m220\,\mathrm{m}.

Exercise 13.3

A spring of stiffness 500N/m500\,\mathrm{N}/\mathrm{m} is compressed by 10cm10\,\mathrm{cm} and released against a 50g50\,\mathrm{g} ball on a frictionless horizontal track. Energy stored; speed of the ball.

Solution

Solution of Exercise 13.3.

Ep=12×500×0.102=2.5JE_p = \tfrac12 \times 500 \times 0.10^2 = 2.5\,\mathrm{J}; v=2×2.5/0.050=10m/sv = \sqrt{2 \times 2.5/0.050} = 10\,\mathrm{m}/\mathrm{s}.

Exercise 13.4

A pendulum of length 1.0m1.0\,\mathrm{m} is released from rest at 6060^\circ: speed at the bottom by energy; then, with Newton’s law, the tension there (in units of mgmg).

Solution

Solution of Exercise 13.4.

12mv2=mg(1cos60)\tfrac12 mv^2 = mg\ell(1 - \cos 60^\circ): v=g=3.1m/sv = \sqrt{g\ell} = 3.1\,\mathrm{m}/\mathrm{s}. Tension: Tmg=mv2/=mgT - mg = mv^2/\ell = mg, T=2mgT = 2mg.

Exercise 13.5 ★★

A block slides from rest down a 3030^\circ slope of height 2.0m2.0\,\mathrm{m}, μd=0.20\mu_d = 0.20. Speed at the bottom by the energy theorem; fraction of the initial potential energy lost to friction.

Solution

Solution of Exercise 13.5.

Friction works μdmgcosα(h/sinα)-\mu_dmg\cos\alpha\,(h/\sin\alpha), so

12mv2=mgh(1μdcotα)=mgh(10.35),\tfrac12 mv^2 = mgh(1 - \mu_d\cot\alpha) = mgh(1 - 0.35),

v=2×9.81×2.0×0.65=5.1m/sv = \sqrt{2 \times 9.81 \times 2.0 \times 0.65} = 5.1\,\mathrm{m}/\mathrm{s}; 35%35\% lost.

Exercise 13.6 ★★

Using Ep=GMm/rE_p = -GMm/r, show that the speed needed to escape the Earth’s attraction from its surface is ve=2GM/Rv_e = \sqrt{2GM/R}, and compute it (M=5.97×1024kgM = 5.97 \times 10^{24}\,\mathrm{kg}, R=6.37×106mR = 6.37 \times 10^{6}\,\mathrm{m}). What is the escape speed from a body with the same density but twice the radius?

Solution

Solution of Exercise 13.6.

Em=12mv2GMm/R0E_m = \tfrac12 mv^2 - GMm/R \geq 0 for the body to reach infinity (where Ep=0E_p = 0 and v0v \to 0): ve=2GM/Rv_e = \sqrt{2GM/R}, i.e.

ve=2×6.67×1011×5.97×10246.37×106=11.2km/s.v_e = \sqrt{\frac{2 \times 6.67\times10^{-11} \times 5.97\times10^{24}}{6.37\times10^6}} = 11.2\,\mathrm{km}/\mathrm{s}.

Same density, MR3M \propto R^3: veRv_e \propto R, so 22.4km/s22.4\,\mathrm{km}/\mathrm{s}.

Exercise 13.7 ★★

Bungee: a 70kg70\,\mathrm{kg} jumper, cord of natural length 20m20\,\mathrm{m} and stiffness 50N/m50\,\mathrm{N}/\mathrm{m}, jumps from rest. Find the maximal stretch of the cord, the total drop, and the maximal tension (in units of the weight). Where is the speed greatest?

Solution

Solution of Exercise 13.7.

Lowest point: mg(0+x)=12kx2mg(\ell_0 + x) = \tfrac12 kx^2, 25x2687x13734=025x^2 - 687x - 13734 = 0, x=41mx = 41\,\mathrm{m}; total drop 61m61\,\mathrm{m}; Tmax=kx=2.0kN=3.0mgT_{\max} = kx = 2.0\,\mathrm{kN} = 3.0\,mg (net 2g2\,\mathrm{g} upward). Speed is greatest where the net force vanishes, kx=mgkx = mg: x=14mx = 14\,\mathrm{m}, 34m34\,\mathrm{m} below the jump.

Exercise 13.8 ★★

A cyclist (80kg80\,\mathrm{kg} with bike) climbs a 5%5\% slope at 5.0m/s5.0\,\mathrm{m}/\mathrm{s} against a drag 12ρCxSv2\tfrac12\rho C_xSv^2 with CxS=0.50m2C_xS = 0.50\,\mathrm{m}^{2}. Power against gravity, against drag, total.

Solution

Solution of Exercise 13.8.

Gravity: mgvsinα=80×9.81×5.0×0.05=196Wmgv\sin\alpha = 80 \times 9.81 \times 5.0 \times 0.05 = 196\,\mathrm{W}; drag: 12ρCxSv3=0.5×1.2×0.50×125=38W\tfrac12\rho C_xSv^3 = 0.5 \times 1.2 \times 0.50 \times 125 = 38\,\mathrm{W}; total 234W234\,\mathrm{W}.

Exercise 13.9 ★★

A particle of mass mm has Ep(x)=E0[(x/a)42(x/a)2]E_p(x) = E_0[(x/a)^4 - 2(x/a)^2]. Find the equilibria and their stability, the depth of the wells, and the angular frequency of small oscillations about a stable one.

Solution

Solution of Exercise 13.9.

Ep=4E0(x3/a4x/a2)=0E_p' = 4E_0(x^3/a^4 - x/a^2) = 0: x=0,±ax = 0, \pm a. Ep=E0(12x2/a44/a2)E_p'' = E_0(12x^2/a^4 - 4/a^2): 4E0/a2-4E_0/a^2 at 00 (unstable maximum, Ep=0E_p = 0), +8E0/a2+8E_0/a^2 at ±a\pm a (stable minima, Ep=E0E_p = -E_0): wells of depth E0E_0; ω02=8E0/ma2\omega_0^2 = 8E_0/ma^2.

Exercise 13.10 ★★★

For the harmonic oscillator Ep=12kx2E_p = \tfrac12 kx^2, show that the phase curves are ellipses, give their semi-axes for energy EE, and explain why all are traversed in the same time.

Solution

Solution of Exercise 13.10.

12mx˙2+12kx2=E\tfrac12 m\dot x^2 + \tfrac12 kx^2 = E: an ellipse of semi-axes xmax=2E/kx_{\max} = \sqrt{2E/k} and x˙max=2E/m\dot x_{\max} = \sqrt{2E/m}. The motion is harmonic with period 2πm/k2\pi\sqrt{m/k} whatever the amplitude (isochronism), so every ellipse takes the same time.

Exercise 13.11 ★★★

A small object slides from rest at the top of a frictionless hemisphere of radius RR. Using energy and Newton’s law, find the angle from the vertical at which it leaves the surface, the height, and its speed then.

Solution

Solution of Exercise 13.11.

Energy: v2=2gR(1cosθ)v^2 = 2gR(1 - \cos\theta). Radial: mgcosθN=mv2/Rmg\cos\theta - N = mv^2/R, so N=mg(3cosθ2)N = mg(3\cos\theta - 2), zero at cosθ=2/3\cos\theta = 2/3 (θ=48\theta = 48^\circ): height 2R/32R/3, speed 2gR/3\sqrt{2gR/3}.

Exercise 13.12 ★★★

A car at 90km/h90\,\mathrm{km}/\mathrm{h} skids to rest in 80m80\,\mathrm{m} with locked wheels. Find μd\mu_d from energy. Where did the energy go? Why does an anti-lock system (wheels kept rolling) stop shorter?

Solution

Solution of Exercise 13.12.

12mv2=μdmgd\tfrac12 mv^2 = \mu_dmgd: μd=625/(2×9.81×80)=0.40\mu_d = 625/(2 \times 9.81 \times 80) = 0.40. Heat in the tire–road contact (and a skid mark). Rolling wheels use static friction, μs>μd\mu_s > \mu_d: a larger braking force, a shorter distance, and steering kept.

13.7 Problem: The energy well of a chemical bond

Problem 13.1

Weekend problem — two atoms, one curve: the potential energy of a bond read as a landscape — its floor, its walls, the small vibrations at the bottom, and the energy it takes to climb out

The interaction of the two atoms of a diatomic molecule is modeled by the potential energy

Ep(r)=ϵ[(σr)122(σr)6],E_p(r) = \epsilon\left[\left(\frac{\sigma}{r}\right)^{12} - 2\left(\frac{\sigma}{r}\right)^{6}\right],

rr the distance between the nuclei. For hydrogen chloride take ϵ=4.4eV\epsilon = 4.4\,\mathrm{eV}, σ=0.127nm\sigma = 0.127\,\mathrm{nm}, and let the light hydrogen atom (m=1.67×1027kgm = 1.67 \times 10^{-27}\,\mathrm{kg}) move while the heavy chlorine stays put. 1eV=1.60×1019J1\,\mathrm{eV} = 1.60 \times 10^{-19}\,\mathrm{J}; kB=1.38×1023J/Kk_B = 1.38 \times 10^{-23}\,\mathrm{J}/\mathrm{K}; NA=6.02×1023N_A = 6.02 \times 10^{23}.

Part I — The landscape.

  1. Give the limits of EpE_p as r0r \to 0 and rr \to \infty, and sketch the curve.
  2. Show that EpE_p has a single minimum, at r=σr = \sigma, of value ϵ-\epsilon.
  3. With u=(σ/r)6u = (\sigma/r)^6, show that Ep=ϵ(u22u)E_p = \epsilon(u^2 - 2u): a parabola in uu. Recover the minimum from it.
  4. Express the force F(r)F(r) (positive if repulsive) and give its sign on each side of σ\sigma.
  5. Compute EpE_p and FF at r=1.2σr = 1.2\sigma (in eV and in nN).
  6. What is the dissociation energy of the molecule, in eV and in kJ/mol\mathrm{kJ}/\mathrm{mol}? (Measured bond energy of HCl: 431kJ/mol431\,\mathrm{kJ}/\mathrm{mol}.)
  7. Check the stability criterion at r=σr = \sigma: compute Ep(σ)E_p''(\sigma).

Part II — Small vibrations.

  1. Write the second-order Taylor expansion of EpE_p around σ\sigma and identify an effective stiffness kk.
  2. Compute kk (in N/m\mathrm{N}/\mathrm{m}).
  3. Deduce the angular frequency, the frequency and the wavelength of the vibration, and its wavenumber 1/λ1/\lambda in cm1\mathrm{cm}^{-1}.
  4. The measured vibration of HCl is at 2990cm12990\,\mathrm{cm}^{-1}. Comment on the agreement and on what the model gets right.
  5. At 300K300\,\mathrm{K} the mean vibrational energy is of order kBTk_BT. Estimate the amplitude of vibration and compare with σ\sigma.
  6. Deduce the maximal speed of the hydrogen atom in that vibration.

Part III — Large oscillations and the way out.

  1. Describe the phase portrait: which energies give closed curves, which open ones, and what separates them?
  2. For Em=ϵ/2E_m = -\epsilon/2, find the two turning points (in units of σ\sigma).
  3. Show that their midpoint lies beyond σ\sigma, and explain why heating (raising the vibrational energy) makes the mean bond length grow — the origin of thermal expansion.
  4. Is the period of large oscillations longer or shorter than 2π/ω02\pi/\omega_0? Argue from the shape of the walls.
  5. From Em=ϵ/2E_m = -\epsilon/2, what energy must be supplied to dissociate the molecule? If given as kinetic energy of the hydrogen atom at r=σr = \sigma, what speed is that?
  6. Two atoms approaching from far apart with positive total energy cannot form a bound molecule by themselves. Why? What is needed?

Part IV — Work, force, and temperature.

  1. Compute the work of the bond force when rr goes from 2σ2\sigma to σ\sigma. Is the force motive or resistive there?
  2. Compute the work an external agent must do to stretch the bond from σ\sigma to 1.5σ1.5\sigma at vanishing speed.
  3. Compute FF at r=0.9σr = 0.9\sigma and compare with r=1.2σr = 1.2\sigma: what does the asymmetry tell about compressing versus stretching?
  4. At what temperature would kBTk_BT equal ϵ\epsilon? Why do molecules dissociate at far lower temperatures in practice?
  5. A better model (Morse) changes the shape of the walls but keeps a minimum at the bond length. Which of the results above survive unchanged in method, and which numbers would change?
  6. Summarize the three ways energy methods were used: to find equilibrium, to find the vibration frequency, and to find what is needed to escape.
Solution

Solution of Problem 13.1.

1. r0r \to 0: ++\infty (the r12r^{-12} term); rr \to \infty: 00^-. A steep wall, a well, a long tail.

2. Ep=ϵ[12σ12/r13+12σ6/r7]=0E_p' = \epsilon[-12\sigma^{12}/r^{13} + 12\sigma^6/r^7] = 0 iff (σ/r)6=1(\sigma/r)^6 = 1: r=σr = \sigma, Ep=ϵ(12)=ϵE_p = \epsilon(1 - 2) = -\epsilon.

3. (σ/r)12=u2(\sigma/r)^{12} = u^2: Ep=ϵ(u22u)=ϵ[(u1)21]E_p = \epsilon(u^2 - 2u) = \epsilon[(u - 1)^2 - 1], minimal at u=1u = 1 (r=σr = \sigma) with value ϵ-\epsilon.

4. F=Ep=(12ϵ/σ)[(σ/r)13(σ/r)7]F = -E_p' = (12\epsilon/\sigma)[(\sigma/r)^{13} - (\sigma/r)^7]: positive (repulsive) for r<σr < \sigma, negative (attractive) beyond.

5. (1/1.2)6=0.335(1/1.2)^6 = 0.335, (1/1.2)12=0.112(1/1.2)^{12} = 0.112: Ep=ϵ(0.1120.670)=0.56ϵ=2.5eVE_p = \epsilon(0.112 - 0.670) = -0.56\epsilon = -2.5\,\mathrm{eV}. (1/1.2)13=0.0935(1/1.2)^{13} = 0.0935, (1/1.2)7=0.279(1/1.2)^7 = 0.279: F=(12ϵ/σ)(0.186)=2.2ϵ/σ=2.2×7.0×1019/1.27×1010=12nNF = (12\epsilon/\sigma)(-0.186) = -2.2\epsilon/\sigma = -2.2 \times 7.0\times10^{-19}/1.27\times10^{-10} = -12\,\mathrm{nN}.

6. ϵ=4.4eV=4.4×96.5=425kJ/mol\epsilon = 4.4\,\mathrm{eV} = 4.4 \times 96.5 = 425\,\mathrm{kJ}/\mathrm{mol}, close to the measured 431kJ/mol431\,\mathrm{kJ}/\mathrm{mol}.

7. Ep=ϵ[156σ12/r1484σ6/r8]E_p'' = \epsilon[156\sigma^{12}/r^{14} - 84\sigma^6/r^8]; at σ\sigma: 72ϵ/σ2>072\epsilon/\sigma^2 > 0: stable.

8. Epϵ+12(72ϵ/σ2)(rσ)2E_p \approx -\epsilon + \tfrac12(72\epsilon/\sigma^2)(r - \sigma)^2: k=72ϵ/σ2k = 72\epsilon/\sigma^2.

9. k=72×7.0×1019/(1.27×1010)2=3.1×103N/mk = 72 \times 7.0\times10^{-19}/(1.27\times10^{-10})^2 = 3.1 \times 10^{3}\,\mathrm{N}/\mathrm{m}.

10. ω0=k/m=3.1×103/1.67×1027=1.4×1015rad/s\omega_0 = \sqrt{k/m} = \sqrt{3.1\times10^3/1.67\times10^{-27}} = 1.4 \times 10^{15}\,\mathrm{rad}/\mathrm{s}; f=2.2×1014Hzf = 2.2 \times 10^{14}\,\mathrm{Hz}; λ=c/f=1.4µm\lambda = c/f = 1.4\,\text{µ}\mathrm{m}; 1/λ=7300cm11/\lambda = 7300\,\mathrm{cm}^{-1}.

11. A factor 2.42.4 too high: the model’s walls are too steep (the r12r^{-12} core is a crude stand-in for a covalent bond), but the order of magnitude — an infrared vibration at 101410^{14} Hz — is right, and the method (curvature of the well) is exactly the one used with better potentials.

12. 12kA2kBT\tfrac12 kA^2 \approx k_BT: A=2kBT/k=8.3×1021/3.1×103=1.6×1012m=0.013σA = \sqrt{2k_BT/k} = \sqrt{8.3\times 10^{-21}/3.1\times10^3} = 1.6 \times 10^{-12}\,\mathrm{m} = 0.013\sigma: the bond barely shivers.

13. vmax=Aω0=1.6×1012×1.4×1015=2.2km/sv_{\max} = A\omega_0 = 1.6\times10^{-12} \times 1.4\times10^{15} = 2.2\,\mathrm{km}/\mathrm{s}.

14. ϵ<Em<0-\epsilon < E_m < 0: closed curves around (σ,0)(\sigma, 0) — vibration, bound; Em>0E_m > 0: open curves — the atoms separate; Em=0E_m = 0: the separatrix, an atom just able to reach infinity.

15. u22u=12u^2 - 2u = -\tfrac12 with u=(σ/r)6u = (\sigma/r)^6: u=1±1/2=1.71u = 1 \pm 1/\sqrt2 = 1.71 or 0.2930.293; r=σu1/6r = \sigma u^{-1/6}: 0.915σ0.915\sigma and 1.23σ1.23\sigma.

16. Midpoint 1.07σ>σ1.07\sigma > \sigma: the outer wall is softer than the inner one, so the atom spends more room outside; the time-averaged distance grows with energy — heated bonds lengthen, solids expand.

17. Longer: the restoring force on the soft outer side is weaker than the harmonic approximation assumes, so the return takes more time.

18. ϵ/2=2.2eV=3.5×1019J\epsilon/2 = 2.2\,\mathrm{eV} = 3.5 \times 10^{-19}\,\mathrm{J}; as kinetic energy at σ\sigma: v=2×3.5×1019/1.67×1027=2.0×104m/sv = \sqrt{2 \times 3.5\times10^{-19}/1.67\times10^{-27}} = 2.0 \times 10^{4}\,\mathrm{m}/\mathrm{s}.

19. Energy conservation: Em>0E_m > 0 stays >0> 0, the curve is open, the atoms fly apart after one encounter. A third body (another molecule, a wall, a photon) must carry the excess away.

20. W=Ep(2σ)Ep(σ)=ϵ(2122×26)+ϵ=0.97ϵ=4.3eVW = E_p(2\sigma) - E_p(\sigma) = \epsilon(2^{-12} - 2 \times 2^{-6}) + \epsilon = 0.97\epsilon = 4.3\,\mathrm{eV}, positive: the attraction pulls the atom in, motive.

21. Wext=Ep(1.5σ)Ep(σ)=ϵ(0.00770.176)+ϵ=0.83ϵ=3.7eVW_{\text{ext}} = E_p(1.5\sigma) - E_p(\sigma) = \epsilon(0.0077 - 0.176) + \epsilon = 0.83\epsilon = 3.7\,\mathrm{eV}.

22. (1/0.9)13=3.93(1/0.9)^{13} = 3.93, (1/0.9)7=2.09(1/0.9)^7 = 2.09: F=+22ϵ/σ=+120nNF = +22\epsilon/\sigma = +120\,\mathrm{nN}, ten times the attraction at 1.2σ1.2\sigma: compressing a bond by 10%10\% costs far more than stretching it by 20%20\%.

23. T=ϵ/kB=7.0×1019/1.38×1023=5.1×104KT = \epsilon/k_B = 7.0\times10^{-19}/1.38\times10^{-23} = 5.1 \times 10^{4}\,\mathrm{K}. Real gases dissociate well below because energies are distributed: a fraction of molecules carries many times kBTk_BT (Chapter 20).

24. Survive: equilibrium at the minimum, k=Epk = E_p'' at the minimum, dissociation energy == depth, turning points from Ep=EmE_p = E_m, asymmetry \Rightarrow expansion. Change: the values of kk, of the frequency and of the turning points.

25. Equilibrium: Ep=0E_p' = 0; vibration: ω02=Ep/m\omega_0^2 = E_p''/m; escape: EmE_m compared with the well’s rim.

Terms defined in this chapter

See all 393 terms in the glossary