Final problem-solving workshop

Two integrals decide almost every normalization and expectation-value problem in PH11001: the Gaussian integral ∫e−ax2dx\int e^{-ax^2}dx, which shows up whenever a wavefunction has a bell-shaped Gaussian profile, and the sinusoidal-squared integral ∫sin⁡2(nπx/L) dx\int \sin^2(n\pi x/L)\,dx, which shows up whenever a wavefunction is confined between hard walls. This workshop is the course's final "integration and exam preparation" component — it introduces no new physics, only drills the toolkit you already have against representative problems spanning the wavefunction and particle-in-a-box material. Master these two families and you can normalize almost any bound-state wavefunction the course throws at you, compute ⟨x⟩\langle x\rangle and ⟨x2⟩\langle x^2\rangle quickly, and recognize the same half-angle trick wherever oscillatory averages appear — from box wavefunctions to wave intensities. The payoff is speed under exam pressure: quoting a standard result instead of re-deriving it from scratch.

The core insight

Every normalization or expectation-value problem in this course reduces to one of two standard integrals — the Gaussian ∫−∞∞e−ax2dx=π/a\int_{-\infty}^{\infty}e^{-ax^2}dx=\sqrt{\pi/a} or the sinusoidal-squared ∫0Lsin⁡2(nπx/L)dx=L/2\int_0^L\sin^2(n\pi x/L)dx=L/2 — which you should quote directly and combine with symmetry shortcuts rather than integrate from scratch.

The core idea

Gaussian-type integrals. These are the standard, quotable results for any Gaussian-shaped wavefunction such as a Gaussian wave packet ψ(x)∝e−x2/2σ2\psi(x)\propto e^{-x^2/2\sigma^2}:

∫−∞∞e−ax2 dx=πa,∫−∞∞x2e−ax2 dx=12aπa,∫−∞∞x e−ax2 dx=0\int_{-\infty}^{\infty} e^{-ax^2}\,dx = \sqrt{\frac{\pi}{a}}, \qquad \int_{-\infty}^{\infty} x^2 e^{-ax^2}\,dx = \frac{1}{2a}\sqrt{\frac{\pi}{a}}, \qquad \int_{-\infty}^{\infty} x\,e^{-ax^2}\,dx = 0

where: aa is a positive constant fixed by the physical width of the packet (for ψ∝e−x2/2σ2\psi\propto e^{-x^2/2\sigma^2}, a=1/σ2a=1/\sigma^2 in ∣ψ∣2|\psi|^2), and the third integral vanishes because x e−ax2x\,e^{-ax^2} is an odd function integrated over a symmetric interval. These three results are the complete toolkit needed to normalize a Gaussian wavefunction and to compute ⟨x⟩\langle x\rangle and ⟨x2⟩\langle x^2\rangle for it, without ever performing the integration by hand.

Sinusoidal-squared integrals. For the infinite-well wavefunction ϕn(x)=Asin⁡(nπx/L)\phi_n(x)=A\sin(n\pi x/L) on [0,L][0,L], normalization requires:

∫0Lsin⁡2 ⁣(nπxL)dx=L2\int_0^L \sin^2\!\left(\frac{n\pi x}{L}\right) dx = \frac{L}{2}

where: nn is the quantum number and LL is the well width; the result is independent of nn because it follows from the half-angle identity sin⁡2θ=12(1−cos⁡2θ)\sin^2\theta=\tfrac12(1-\cos2\theta) — the oscillating cos⁡(2nπx/L)\cos(2n\pi x/L) term completes a whole number of cycles over [0,L][0,L] and integrates to exactly zero, leaving precisely half of LL. The same half-angle trick governs time-averages of oscillatory quantities elsewhere in physics, e.g. ⟨cos⁡2(ωt)⟩=12\langle\cos^2(\omega t)\rangle=\tfrac12, which is why average power and intensity formulas (such as I=12cε0E02I=\tfrac12 c\varepsilon_0 E_0^2) always carry a factor of one-half.

Expectation values and symmetry. Once a wavefunction is normalized via ∫∣ψ∣2dx=1\int|\psi|^2dx=1, expectation values take the form ⟨x⟩=∫ψ∗x ψ dx\langle x\rangle=\int\psi^*x\,\psi\,dx and ⟨x2⟩=∫ψ∗x2 ψ dx\langle x^2\rangle=\int\psi^*x^2\,\psi\,dx. For a real, symmetric probability density — such as ∣ϕn(x)∣2|\phi_n(x)|^2 about the midpoint of a box, or ∣ψ(x)∣2|\psi(x)|^2 about the origin for a centred Gaussian — you rarely need to grind through the integral by hand: symmetry alone fixes ⟨x⟩\langle x\rangle, and only ⟨x2⟩\langle x^2\rangle needs the standard-integral machinery above. This division of labour — symmetry for the first moment, standard integrals for the second — is exactly the pattern examiners test: a quick qualitative argument for ⟨x⟩\langle x\rangle, followed by a short calculation for ⟨x2⟩\langle x^2\rangle using whichever of the two integral families matches the wavefunction's shape.

Choosing the right family. Before reaching for either toolkit, identify the wavefunction's shape and domain: a bell-shaped exponential extending over all space signals the Gaussian family, while a sine confined between hard walls on a finite interval signals the sinusoidal family. The two never mix within a single problem in this course — a wavefunction is either an unbounded Gaussian packet or a bounded box eigenstate, never both at once — so this identification step alone resolves most of the ambiguity about which formula to quote.

Figure
Two probability density curves side by side, a Gaussian bell centred at the origin and a sine-squared hump confined between two walls, each shaded symmetric about its own midpoint
Two families of integrals, two symmetric densities Gaussian packet, symmetric about x = 0 midpoint particle in a box, symmetric about x = L over 2 midpoint 0 L Symmetry about the midpoint fixes the mean position without integration
Derivation ANormalizing a particle-in-a-box wavefunction with the sin-squared integral

The goal is to show, using only the standard sinusoidal-squared integral above and no fresh integration by parts, that ϕn(x)=2/Lsin⁡(nπx/L)\phi_n(x)=\sqrt{2/L}\sin(n\pi x/L) is a correctly normalized wavefunction for an infinite well of width LL, for every quantum number nn — the single most common normalization check in this course, and the one that recurs in the worked example below.

Step 1. Write the normalization condition to be verified: total probability over the well must equal one, ∫0L∣ϕn(x)∣2 dx=1\int_0^L |\phi_n(x)|^2\,dx = 1. Substitute the trial wavefunction ϕn(x)=Asin⁡(nπx/L)\phi_n(x)=A\sin(n\pi x/L) with AA left as an unknown constant, so the condition becomes A2∫0Lsin⁡2(nπx/L) dx=1A^2\int_0^L \sin^2(n\pi x/L)\,dx = 1.

Step 2. Evaluate the integral using the standard sinusoidal-squared result quoted above rather than re-deriving it: ∫0Lsin⁡2(nπx/L) dx=L/2\int_0^L \sin^2(n\pi x/L)\,dx = L/2, valid for every positive integer nn because the cos⁡(2nπx/L)\cos(2n\pi x/L) term always completes a whole number of cycles over [0,L][0,L].

Step 3. Substitute back: A2⋅L2=1A^2\cdot\dfrac{L}{2}=1, so A2=2/LA^2=2/L and therefore A=2/LA=\sqrt{2/L}. This confirms the standard normalization constant used throughout the particle-in-a-box material, for any well width LL and any quantum number nn.

Step 4. Interpret the result: because the integral L/2L/2 does not depend on nn, the same normalization constant A=2/LA=\sqrt{2/L} works for the ground state and every excited state of the box — only the number of oscillations changes with nn, not the total probability.

Worked example

Normalize a Gaussian wave packet and find ⟨x2⟩\langle x^2\rangle. A particle's wavefunction is ψ(x)=Ne−x2/2σ2\psi(x)=Ne^{-x^2/2\sigma^2} with σ=2.0×10−10 m\sigma=2.0\times10^{-10}\ \text{m}. Find the normalization constant NN and the expectation value ⟨x2⟩\langle x^2\rangle.

Step 1. Set up the normalization integral. Requiring ∫−∞∞∣ψ∣2dx=1\int_{-\infty}^{\infty}|\psi|^2dx=1 gives N2∫−∞∞e−x2/σ2dx=N2πσ2=1N^2\int_{-\infty}^{\infty}e^{-x^2/\sigma^2}dx=N^2\sqrt{\pi\sigma^2}=1, using the standard Gaussian result with a=1/σ2a=1/\sigma^2.

Step 2. Solve for NN. This gives N=(πσ2)−1/4N=(\pi\sigma^2)^{-1/4}. Numerically, σπ=(2.0×10−10)(1.772)=3.545×10−10\sigma\sqrt\pi=(2.0\times10^{-10})(1.772)=3.545\times10^{-10}, so N≈5.31×104 m−1/2N\approx 5.31\times10^{4}\ \text{m}^{-1/2}.

Step 3. Compute ⟨x2⟩\langle x^2\rangle using the second standard Gaussian integral: ⟨x2⟩=N2∫−∞∞x2e−x2/σ2dx=N2⋅σ22πσ2\langle x^2\rangle=N^2\int_{-\infty}^{\infty}x^2e^{-x^2/\sigma^2}dx=N^2\cdot\dfrac{\sigma^2}{2}\sqrt{\pi\sigma^2}. Since N2πσ2=1N^2\sqrt{\pi\sigma^2}=1 from Step 1, this simplifies immediately.

Step 4. Read off the answer: ⟨x2⟩=σ22=(2.0×10−10)22=2.0×10−20 m2\langle x^2\rangle=\dfrac{\sigma^2}{2}=\dfrac{(2.0\times10^{-10})^2}{2}=\mathbf{2.0\times10^{-20}\ m^2}. The spread ⟨x2⟩=σ/2≈1.4×10−10 m\sqrt{\langle x^2\rangle}=\sigma/\sqrt2\approx1.4\times10^{-10}\ \text{m} is directly the Δx\Delta x that would enter the uncertainty relation Δx Δp≥ℏ/2\Delta x\,\Delta p\geq\hbar/2 from the wave-mechanics block.

Competitive stretch

The same packet's ground-state cousin, ϕ1(x)=2/Lsin⁡(πx/L)\phi_1(x)=\sqrt{2/L}\sin(\pi x/L) on a box of width LL, gives ⟨x⟩=L/2\langle x\rangle=L/2 purely by symmetry — no integration needed, since ∣ϕ1(x)∣2|\phi_1(x)|^2 is symmetric about the midpoint x=L/2x=L/2 (check ∣ϕ1(L−x)∣2=∣ϕ1(x)∣2|\phi_1(L-x)|^2=|\phi_1(x)|^2 using sin⁡(π−θ)=sin⁡θ\sin(\pi-\theta)=\sin\theta). Recognizing this shortcut in an exam saves the entire integration-by-parts calculation that a brute-force approach would require.

Common traps

  • Using ∫−∞∞e−ax2dx=π/a\int_{-\infty}^{\infty}e^{-ax^2}dx=\sqrt{\pi/a} or ∫0Lsin⁡2(nπx/L)dx=L/2\int_0^L\sin^2(n\pi x/L)dx=L/2 outside their valid domain — the Gaussian result needs the full real line, the sine-squared result needs a full-period interval like [0,L][0,L]; applying L/2L/2 to a partial interval such as [0,L/2][0,L/2] gives the wrong answer because the cosine term no longer completes whole cycles.
  • Grinding through integration by parts when a symmetry check would give ⟨x⟩\langle x\rangle instantly — always check whether ∣ψ∣2|\psi|^2 is even, odd, or symmetric about the interval's midpoint first, since a fully worked-out integral and the symmetry shortcut must agree, and the shortcut is far faster under exam conditions.
  • Treating the normalization constant NN or AA as a free parameter chosen for convenience — it is fixed uniquely by requiring total probability equal one, and re-using a different value later in the same problem (say, to simplify an unrelated expression) silently breaks the normalization.
  • Dropping the factor of one-half that comes from ⟨sin⁡2⟩=⟨cos⁡2⟩=12\langle\sin^2\rangle=\langle\cos^2\rangle=\tfrac12 when computing average power, intensity, or energy — a frequent slip under time pressure, and one that produces an answer exactly double the correct value.
  • Forgetting that ∫0Lsin⁡2(nπx/L)dx=L/2\int_0^L\sin^2(n\pi x/L)dx=L/2 holds for every quantum number nn — students sometimes assume the normalization constant A=2/LA=\sqrt{2/L} changes between the ground state and excited states, when in fact it is identical for all nn.

Ask the tutor: want a few more practice problems mixing Gaussian and sinusoidal integrals, or a walkthrough applying this toolkit to a specific past-exam question?

Why do I need to memorize these integrals instead of just integrating them each time?

Because exam time is the real constraint, not mathematical ability. Deriving ∫e−ax2dx=π/a\int e^{-ax^2}dx=\sqrt{\pi/a} from scratch requires a trick (squaring the integral and switching to polar coordinates) that eats several minutes you do not have during a problem set or exam. Quoting the standard result and moving straight to the physics — solving for a normalization constant or an expectation value — is what the course expects, and what distinguishes a fluent answer from a slow one.

How do I know whether to use the Gaussian integral or the sinusoidal one?

Look at the wavefunction's shape and domain. If the wavefunction is a bell-shaped exponential like e−x2/2σ2e^{-x^2/2\sigma^2} extending over all space (a free-particle wave packet), you need the Gaussian family. If the wavefunction is a sine confined between hard walls, like sin⁡(nπx/L)\sin(n\pi x/L) on [0,L][0,L] (a particle in a box), you need the sinusoidal-squared family. The two never mix within a single problem in this course, so identifying the wavefunction type immediately tells you which toolkit applies.

Does the symmetry shortcut for ⟨x⟩\langle x\rangle always work, or only sometimes?

It only works when the probability density ∣ψ(x)∣2|\psi(x)|^2 is symmetric about the centre of the interval you are integrating over — true for every particle-in-a-box eigenstate about x=L/2x=L/2, and true for a centred Gaussian packet about x=0x=0. If the wavefunction were shifted off-centre, or were a superposition of states with different symmetries, the shortcut would fail and you would need to integrate directly. Always verify the symmetry explicitly (as in the worked example) before invoking it, rather than assuming it.

Before this