# Day 238: Week 34 Review - Banach and Hilbert Spaces ## Schedule Overview (8 hours) | Session | Duration | Focus | |---------|----------|-------| | Morning | 3 hours | Comprehensive review and concept map | | Afternoon | 3 hours | Problem-solving marathon | | Evening | 2 hours | Quantum applications and synthesis | ## Review Objectives By the end of today, you will: 1. **Synthesize** all concepts from Week 34 into a coherent framework 2. **Master** the relationships between normed, Banach, inner product, and Hilbert spaces 3. **Solve** comprehensive problems integrating multiple concepts 4. **Connect** the mathematical framework to quantum mechanics 5. **Prepare** for Week 35's study of bounded operators --- ## 1. Week 34 Concept Map ### 1.1 The Hierarchy of Spaces ``` Vector Space │ │ + norm ▼ Normed Space ──────────────────────────┐ │ │ │ + completeness │ + inner product ▼ ▼ Banach Space Inner Product Space │ │ │ │ + completeness │ ▼ └───────────────────────────► Hilbert Space │ │ + separability ▼ Separable Hilbert Space (≅ ℓ² or ℂⁿ) ``` ### 1.2 Key Relationships | Property | Normed | Banach | Inner Product | Hilbert | |----------|--------|--------|---------------|---------| | Vector space | ✓ | ✓ | ✓ | ✓ | | Norm defined | ✓ | ✓ | ✓ (induced) | ✓ (induced) | | Complete | ✗ | ✓ | ✗ | ✓ | | Inner product | ✗ | ✗ | ✓ | ✓ | | Parallelogram law | ✗ | ✗ | ✓ | ✓ | | Self-dual | ✗ | ✗ | ✗ | ✓ | --- ## 2. Master Theorem Summary ### 2.1 Definitions $$\boxed{\begin{aligned} &\textbf{Norm:} && \|x\| \geq 0, \|\alpha x\| = |\alpha|\|x\|, \|x+y\| \leq \|x\| + \|y\| \\[5pt] &\textbf{Banach Space:} && \text{Complete normed space} \\[5pt] &\textbf{Inner Product:} && \langle x,y\rangle = \overline{\langle y,x\rangle}, \text{ linear in 1st arg, } \langle x,x\rangle \geq 0 \\[5pt] &\textbf{Hilbert Space:} && \text{Complete inner product space} \\[5pt] &\textbf{Separable:} && \text{Has countable dense subset (equiv: countable ONB)} \end{aligned}}$$ ### 2.2 Fundamental Inequalities $$\boxed{\begin{aligned} &\textbf{Cauchy-Schwarz:} && |\langle x,y\rangle| \leq \|x\| \|y\| \\[5pt] &\textbf{Triangle:} && \|x+y\| \leq \|x\| + \|y\| \\[5pt] &\textbf{Parallelogram:} && \|x+y\|^2 + \|x-y\|^2 = 2\|x\|^2 + 2\|y\|^2 \\[5pt] &\textbf{Bessel:} && \sum_n |\langle e_n,x\rangle|^2 \leq \|x\|^2 \\[5pt] &\textbf{Parseval:} && \sum_n |\langle e_n,x\rangle|^2 = \|x\|^2 \text{ (for ONB)} \end{aligned}}$$ ### 2.3 Key Theorems | Theorem | Statement | |---------|-----------| | **Polarization** | $$\langle x,y\rangle = \frac{1}{4}\sum_{k=0}^3 i^k \|x + i^k y\|^2$$ | | **Jordan-von Neumann** | Parallelogram law ⟺ norm from inner product | | **Riesz-Fischer** | $$L^2$$ is complete | | **Gram-Schmidt** | LI set → orthonormal set (same span) | | **Riesz Representation** | $$f \in \mathcal{H}^* \Rightarrow \exists! y: f(x) = \langle y, x\rangle$$ | | **Classification** | All sep. inf-dim Hilbert spaces ≅ $$\ell^2$$ | --- ## 3. Quantum Mechanics Dictionary ### 3.1 Mathematical Structure ↔ Physical Interpretation | Mathematics | Quantum Mechanics | |-------------|-------------------| | Hilbert space $$\mathcal{H}$$ | State space | | Vector $$\psi \in \mathcal{H}$$, $$\|\psi\| = 1$$ | Pure quantum state $$\|\psi\rangle$$ | | Inner product $$\langle\phi,\psi\rangle$$ | Probability amplitude $$\langle\phi\|\psi\rangle$$ | | $$\|\langle\phi,\psi\rangle\|^2$$ | Transition probability | | Orthogonality $$\langle\phi,\psi\rangle = 0$$ | Distinguishable states | | Orthonormal basis $$\{e_n\}$$ | Complete set of states (CSCO) | | Fourier coefficients $$c_n = \langle e_n, \psi\rangle$$ | Probability amplitudes | | Parseval: $$\sum \|c_n\|^2 = 1$$ | Probability conservation | | Completeness: $$\sum \|e_n\rangle\langle e_n\| = I$$ | Resolution of identity | | Riesz representation | Bra-ket duality | | Isomorphism $$L^2 \cong \ell^2$$ | Position ↔ energy representation | ### 3.2 The Postulates in Hilbert Space Language 1. **States**: Rays in a separable Hilbert space $$\mathcal{H}$$ 2. **Observables**: Self-adjoint operators on $$\mathcal{H}$$ (Week 35) 3. **Measurement outcomes**: Eigenvalues, with probabilities $$|c_n|^2$$ 4. **Dynamics**: Unitary evolution preserving $$\|\psi\|$$ (Schrödinger equation) --- ## 4. Comprehensive Problem Set ### Problem 1: Space Classification **Determine whether each space is a Banach space, Hilbert space, both, or neither.** (a) $$(\mathbb{R}^3, \|\cdot\|_2)$$ (b) $$(\ell^1, \|\cdot\|_1)$$ (c) $$(C[0,1], \|\cdot\|_\infty)$$ (d) $$(C[0,1], \|\cdot\|_2)$$ where $$\|f\|_2 = \sqrt{\int_0^1 |f|^2}$$ (e) $$(L^2[0,1], \|\cdot\|_2)$$ **Solution**: (a) **Hilbert space** (hence also Banach). Finite-dimensional, norm from standard inner product. (b) **Banach but not Hilbert**. Complete but violates parallelogram law: $$x = (1,0,0,\ldots), y = (0,1,0,\ldots)$$ $$\|x+y\|_1^2 + \|x-y\|_1^2 = 4 + 4 = 8 \neq 4 = 2(1+1)$$ ✗ (c) **Banach but not Hilbert**. Complete (uniform limit of continuous functions is continuous). Violates parallelogram law with similar example. (d) **Neither Banach nor Hilbert**. Not complete (Day 232 counterexample). However, if it were complete, it would be Hilbert (norm satisfies parallelogram law). (e) **Hilbert space** (hence also Banach). Complete by Riesz-Fischer. $$\square$$ --- ### Problem 2: Inner Product Verification **Prove or disprove: $$\langle f, g\rangle = \int_0^1 f(x)g(x) + f'(x)g'(x) \, dx$$ defines an inner product on $$C^1[0,1]$$ (continuously differentiable functions).** **Solution**: **(IP1) Conjugate symmetry**: For real functions, $$\langle g, f\rangle = \int_0^1 gf + g'f' = \int_0^1 fg + f'g' = \langle f, g\rangle$$ ✓ **(IP2) Linearity**: Clear from linearity of integral and derivative. ✓ **(IP3) Positive definiteness**: $$\langle f, f\rangle = \int_0^1 f(x)^2 + f'(x)^2 \, dx \geq 0$$ ✓ If $$\langle f, f\rangle = 0$$, then $$\int_0^1 f^2 = 0$$ and $$\int_0^1 (f')^2 = 0$$. Since $$f$$ is continuous and $$f^2 \geq 0$$, we have $$f(x) = 0$$ for all $$x$$. ✓ **Conclusion**: Yes, this is an inner product on $$C^1[0,1]$$. $$\square$$ **Note**: This is the $$H^1$$ Sobolev inner product, important in PDEs and quantum mechanics. --- ### Problem 3: Gram-Schmidt in L² **Apply Gram-Schmidt to $$\{1, e^x, e^{2x}\}$$ in $$L^2[0,1]$$ (first two steps).** **Solution**: **Step 1**: $$p_0 = 1$$ $$\|1\|^2 = \int_0^1 1 \, dx = 1$$ $$e_0 = 1$$ **Step 2**: $$\langle e_0, e^x\rangle = \int_0^1 e^x \, dx = e - 1$$ $$u_1 = e^x - (e-1) \cdot 1 = e^x - e + 1$$ $$\|u_1\|^2 = \int_0^1 (e^x - e + 1)^2 \, dx$$ Let me compute this: $$= \int_0^1 e^{2x} - 2(e-1)e^x + (e-1)^2 \, dx$$ $$= \frac{e^2 - 1}{2} - 2(e-1)(e-1) + (e-1)^2 = \frac{e^2-1}{2} - (e-1)^2$$ $$= \frac{e^2-1}{2} - e^2 + 2e - 1 = \frac{e^2-1-2e^2+4e-2}{2} = \frac{-e^2+4e-3}{2} = \frac{-(e-1)(e-3)}{2}$$ $$= \frac{(e-1)(3-e)}{2}$$ Numerically: $$\approx \frac{1.718 \times 0.282}{2} \approx 0.242$$ $$e_1 = \frac{e^x - e + 1}{\|u_1\|} \approx \frac{e^x - e + 1}{0.492}$$ $$\square$$ --- ### Problem 4: Parseval Application **Use Parseval to prove: $$\sum_{n=1}^\infty \frac{1}{(2n-1)^2} = \frac{\pi^2}{8}$$** **Solution**: Consider $$f(x) = 1$$ on $$[0, \pi]$$, extended as an odd function to $$[-\pi, \pi]$$: $$f(x) = \text{sign}(x) = \begin{cases} 1 & x > 0 \\ -1 & x < 0 \end{cases}$$ The Fourier sine series is: $$f(x) = \sum_{n=1}^\infty b_n \sin(nx)$$ where $$b_n = \frac{2}{\pi}\int_0^\pi \sin(nx) \, dx = \frac{2}{\pi} \cdot \frac{1 - (-1)^n}{n}$$ For odd $$n = 2k-1$$: $$b_{2k-1} = \frac{4}{\pi(2k-1)}$$ For even $$n$$: $$b_n = 0$$ Parseval (for sine series on $$[0, \pi]$$): $$\frac{2}{\pi}\int_0^\pi f(x)^2 \, dx = \sum_{n=1}^\infty b_n^2$$ Left side: $$\frac{2}{\pi} \cdot \pi = 2$$ Right side: $$\sum_{k=1}^\infty \frac{16}{\pi^2(2k-1)^2}$$ Therefore: $$2 = \frac{16}{\pi^2} \sum_{k=1}^\infty \frac{1}{(2k-1)^2}$$ $$\boxed{\sum_{n=1}^\infty \frac{1}{(2n-1)^2} = \frac{\pi^2}{8}}$$ $$\square$$ --- ### Problem 5: Completeness and Approximation **Prove that continuous functions are dense in $$L^2[0,1]$$, then explain why this doesn't contradict the fact that $$C[0,1]$$ with $$L^2$$ norm is incomplete.** **Solution**: **Part 1: Density** Let $$f \in L^2[0,1]$$ and $$\varepsilon > 0$$. We need $$g \in C[0,1]$$ with $$\|f - g\|_2 < \varepsilon$$. Step functions are dense in $$L^2$$ (can approximate any measurable function). Any step function can be approximated in $$L^2$$ by a continuous function (smooth out the jumps over intervals of width $$\delta$$). Therefore, continuous functions are dense in $$L^2$$. **Part 2: No contradiction** $$(C[0,1], \|\cdot\|_2)$$ being incomplete means: there exist Cauchy sequences in $$C[0,1]$$ whose limits (in the $$L^2$$ metric) are NOT continuous functions. These limits ARE in $$L^2[0,1]$$ (which is complete), just not in the subset $$C[0,1]$$. $$C[0,1]$$ is dense but not closed in $$L^2[0,1]$$. The completion of $$(C[0,1], \|\cdot\|_2)$$ IS $$L^2[0,1]$$. Analogy: $$\mathbb{Q}$$ is dense in $$\mathbb{R}$$ but not complete. The completion of $$\mathbb{Q}$$ is $$\mathbb{R}$$. $$\square$$ --- ### Problem 6: Quantum Mechanics Application **A quantum harmonic oscillator is prepared in the state:** $$|\psi\rangle = \frac{1}{2}|0\rangle + \frac{1}{2}|1\rangle + \frac{1}{\sqrt{2}}|2\rangle$$ **(a) Verify normalization using Parseval.** **(b) What is the probability of measuring energy $$E_n = (n + 1/2)\hbar\omega$$?** **(c) Compute $$\langle\psi|0\rangle\langle 0|\psi\rangle$$ and interpret.** **Solution**: **(a)** Coefficients: $$c_0 = 1/2$$, $$c_1 = 1/2$$, $$c_2 = 1/\sqrt{2}$$ Parseval: $$\||\psi\rangle\|^2 = \sum_n |c_n|^2 = \frac{1}{4} + \frac{1}{4} + \frac{1}{2} = 1$$ ✓ **(b)** Probabilities: - $$P(E_0 = \hbar\omega/2) = |c_0|^2 = 1/4 = 25\%$$ - $$P(E_1 = 3\hbar\omega/2) = |c_1|^2 = 1/4 = 25\%$$ - $$P(E_2 = 5\hbar\omega/2) = |c_2|^2 = 1/2 = 50\%$$ - $$P(E_n) = 0$$ for $$n \geq 3$$ **(c)** $$\langle\psi|0\rangle = \overline{c_0} = 1/2$$ $$\langle 0|\psi\rangle = c_0 = 1/2$$ $$\langle\psi|0\rangle\langle 0|\psi\rangle = \frac{1}{2} \cdot \frac{1}{2} = \frac{1}{4}$$ This equals $$|c_0|^2 = P(E_0)$$, the probability of measuring the ground state energy. The operator $$|0\rangle\langle 0|$$ is the projector onto the ground state. $$\square$$ --- ## 5. Synthesis: The Big Picture ### 5.1 Why This Week Matters Week 34 established the **mathematical language of quantum mechanics**: 1. **Hilbert spaces** provide the arena where quantum states live 2. **Inner products** give probability amplitudes 3. **Orthonormal bases** correspond to measurement outcomes 4. **Parseval's identity** ensures probability conservation 5. **Completeness** guarantees physical limits exist 6. **Separability** allows us to work with countable bases ### 5.2 The Road Ahead **Week 35: Bounded Linear Operators** - Operators are the mathematical representation of observables - The spectral theorem connects operators to measurement outcomes - Adjoint operators lead to self-adjoint (Hermitian) observables **Future topics building on this foundation**: - Unbounded operators (momentum, position) - Spectral theory - Quantum dynamics (unitary evolution) - Tensor products (composite systems) --- ## 6. Self-Assessment Checklist ### Definitions (Can you state precisely?) - [ ] Norm and its three axioms - [ ] Banach space - [ ] Inner product and its axioms - [ ] Hilbert space - [ ] Orthonormal set vs. orthonormal basis - [ ] Separability ### Theorems (Can you state and prove?) - [ ] Every inner product induces a norm - [ ] Cauchy-Schwarz inequality - [ ] Parallelogram law - [ ] Polarization identity - [ ] Bessel's inequality - [ ] Parseval's identity - [ ] Riesz-Fischer ($$L^2$$ is complete) - [ ] Gram-Schmidt process - [ ] Riesz representation theorem - [ ] Classification of separable Hilbert spaces ### Computations (Can you perform?) - [ ] Verify norm/inner product axioms - [ ] Compute $$\ell^p$$ norms - [ ] Apply Gram-Schmidt - [ ] Compute Fourier coefficients - [ ] Use Parseval to evaluate sums - [ ] Construct explicit isomorphisms ### Connections (Can you explain?) - [ ] Why quantum states live in Hilbert space - [ ] Why Parseval = probability conservation - [ ] Why completeness relation = resolution of identity - [ ] Why $$L^2 \cong \ell^2$$ (all representations equivalent) --- ## 7. Computational Lab: Week Review ```python """ Day 238 Computational Lab: Week 34 Review Comprehensive exercises synthesizing Banach and Hilbert space concepts """ import numpy as np import matplotlib.pyplot as plt from scipy.integrate import quad from scipy.special import hermite # ============================================================ # Part 1: Space Classification Visual # ============================================================ def visualize_space_hierarchy(): """ Create a visual representation of the space hierarchy. """ fig, ax = plt.subplots(figsize=(12, 8)) # Define positions for boxes positions = { 'Vector Space': (0.5, 0.9), 'Normed Space': (0.5, 0.7), 'Banach Space': (0.25, 0.5), 'Inner Product Space': (0.75, 0.5), 'Hilbert Space': (0.5, 0.3), 'Sep. Hilbert': (0.5, 0.1) } # Define arrows arrows = [ ('Vector Space', 'Normed Space', '+ norm'), ('Normed Space', 'Banach Space', '+ complete'), ('Normed Space', 'Inner Product Space', '+ inner prod'), ('Banach Space', 'Hilbert Space', '+ inner prod'), ('Inner Product Space', 'Hilbert Space', '+ complete'), ('Hilbert Space', 'Sep. Hilbert', '+ separable') ] # Draw boxes for name, (x, y) in positions.items(): bbox = dict(boxstyle='round,pad=0.3', facecolor='lightblue', alpha=0.7) ax.text(x, y, name, ha='center', va='center', fontsize=12, bbox=bbox, transform=ax.transAxes) # Draw arrows for start, end, label in arrows: x1, y1 = positions[start] x2, y2 = positions[end] ax.annotate('', xy=(x2, y2+0.04), xytext=(x1, y1-0.04), arrowprops=dict(arrowstyle='->', color='black'), transform=ax.transAxes) ax.text((x1+x2)/2 + 0.05, (y1+y2)/2, label, fontsize=9, transform=ax.transAxes) # Add examples examples = { (0.05, 0.7): r'Ex: $\mathbb{R}^n$', (0.05, 0.5): r'Ex: $\ell^1$, $C[0,1]$', (0.95, 0.5): r'Ex: incomplete $L^2$', (0.05, 0.3): r'Ex: $L^2$, $\ell^2$', (0.05, 0.1): r'All $\cong \ell^2$!' } for (x, y), text in examples.items(): ax.text(x, y, text, fontsize=10, style='italic', transform=ax.transAxes) ax.axis('off') ax.set_title('Hierarchy of Spaces (Week 34)', fontsize=14, fontweight='bold') plt.tight_layout() plt.savefig('space_hierarchy.png', dpi=150) plt.show() # ============================================================ # Part 2: Inequality Verification Suite # ============================================================ def verify_all_inequalities(): """ Numerically verify all key inequalities from the week. """ np.random.seed(42) print("Inequality Verification Suite") print("=" * 60) # Generate random vectors n = 100 x = np.random.randn(n) + 1j * np.random.randn(n) y = np.random.randn(n) + 1j * np.random.randn(n) # 1. Cauchy-Schwarz lhs_cs = np.abs(np.vdot(x, y)) rhs_cs = np.linalg.norm(x) * np.linalg.norm(y) print(f"1. Cauchy-Schwarz: || = {lhs_cs:.6f} ≤ ||x|| ||y|| = {rhs_cs:.6f}") print(f" Satisfied: {lhs_cs <= rhs_cs + 1e-10}") # 2. Triangle inequality lhs_tri = np.linalg.norm(x + y) rhs_tri = np.linalg.norm(x) + np.linalg.norm(y) print(f"\n2. Triangle: ||x+y|| = {lhs_tri:.6f} ≤ ||x|| + ||y|| = {rhs_tri:.6f}") print(f" Satisfied: {lhs_tri <= rhs_tri + 1e-10}") # 3. Parallelogram law lhs_para = np.linalg.norm(x+y)**2 + np.linalg.norm(x-y)**2 rhs_para = 2 * (np.linalg.norm(x)**2 + np.linalg.norm(y)**2) print(f"\n3. Parallelogram: ||x+y||² + ||x-y||² = {lhs_para:.6f}") print(f" 2(||x||² + ||y||²) = {rhs_para:.6f}") print(f" Equal: {np.abs(lhs_para - rhs_para) < 1e-10}") # 4. Bessel's inequality (using standard basis as ONB) # x is already in ℓ^2 (finite vector), standard basis x_for_bessel = np.random.randn(1000) + 1j * np.random.randn(1000) # Use first 100 basis vectors coeffs = x_for_bessel[:100] lhs_bessel = np.sum(np.abs(coeffs)**2) rhs_bessel = np.linalg.norm(x_for_bessel)**2 print(f"\n4. Bessel (first 100 coeffs): Σ|c_n|² = {lhs_bessel:.6f}") print(f" ||x||² = {rhs_bessel:.6f}") print(f" Satisfied: {lhs_bessel <= rhs_bessel + 1e-10}") # 5. Parseval (full basis) lhs_parseval = np.sum(np.abs(x_for_bessel)**2) print(f"\n5. Parseval (full): Σ|c_n|² = {lhs_parseval:.6f}") print(f" ||x||² = {rhs_bessel:.6f}") print(f" Equal: {np.abs(lhs_parseval - rhs_bessel) < 1e-10}") # ============================================================ # Part 3: Quantum Mechanics Synthesis # ============================================================ def quantum_synthesis(): """ Comprehensive quantum mechanics example using all concepts. """ print("\n" + "=" * 60) print("Quantum Mechanics Synthesis") print("=" * 60) x = np.linspace(-10, 10, 1000) # Harmonic oscillator eigenstates def psi_n(x, n): Hn = hermite(n) norm = 1 / np.sqrt(2**n * np.math.factorial(n) * np.sqrt(np.pi)) return norm * Hn(x) * np.exp(-x**2 / 2) # Create a quantum state coeffs = np.array([0.5, 0.5, 1/np.sqrt(2), 0, 0]) # c_0, c_1, c_2, c_3, c_4 print("\n1. State: |ψ⟩ = 0.5|0⟩ + 0.5|1⟩ + (1/√2)|2⟩") # Verify normalization (Parseval) norm_sq = np.sum(np.abs(coeffs)**2) print(f"\n2. Normalization (Parseval): Σ|c_n|² = {norm_sq:.6f}") print(f" State is {'normalized' if np.abs(norm_sq - 1) < 1e-10 else 'NOT normalized'}") # Probabilities print("\n3. Measurement probabilities:") for n, c in enumerate(coeffs): if np.abs(c) > 1e-10: print(f" P(E_{n} = {n+0.5}ℏω) = |c_{n}|² = {np.abs(c)**2:.4f}") # Inner products (transition amplitudes) print("\n4. Transition amplitudes:") # Inner product with |+⟩ = (|0⟩ + |1⟩)/√2 plus_coeffs = np.array([1/np.sqrt(2), 1/np.sqrt(2), 0, 0, 0]) amplitude = np.vdot(coeffs, plus_coeffs) # <+|ψ> print(f" ⟨+|ψ⟩ = {amplitude:.4f}") print(f" P(+|ψ) = |⟨+|ψ⟩|² = {np.abs(amplitude)**2:.4f}") # Completeness relation print("\n5. Completeness relation: Σ|n⟩⟨n| = I") print(" Verification: Σ|⟨n|ψ⟩|² = Σ|c_n|² = 1 ✓") # Position representation psi_x = sum(c * psi_n(x, n) for n, c in enumerate(coeffs)) norm_x = np.trapz(np.abs(psi_x)**2, x) print(f"\n6. Position representation:") print(f" ∫|ψ(x)|²dx = {norm_x:.6f} (should equal 1)") # Visualize fig, axes = plt.subplots(2, 2, figsize=(14, 10)) # Energy representation ax = axes[0, 0] ax.bar(range(len(coeffs)), np.abs(coeffs)**2, color='steelblue', alpha=0.7) ax.set_xlabel('$n$') ax.set_ylabel('$|c_n|^2$') ax.set_title('Energy Representation') ax.set_xticks(range(len(coeffs))) # Position representation ax = axes[0, 1] ax.fill_between(x, np.abs(psi_x)**2, alpha=0.5, color='green') ax.plot(x, np.abs(psi_x)**2, 'g-', linewidth=2) ax.set_xlabel('$x$') ax.set_ylabel('$|\\psi(x)|^2$') ax.set_title('Position Representation') ax.set_xlim(-6, 6) # Individual eigenstates ax = axes[1, 0] for n in range(3): ax.plot(x, psi_n(x, n) + n*0.5, linewidth=2, label=f'$\\psi_{n}(x)$') ax.set_xlabel('$x$') ax.set_ylabel('$\\psi_n(x)$ (offset)') ax.set_title('Harmonic Oscillator Eigenstates') ax.legend() ax.set_xlim(-6, 6) # Orthonormality verification ax = axes[1, 1] n_max = 5 overlap = np.zeros((n_max, n_max)) for i in range(n_max): for j in range(n_max): overlap[i, j] = np.trapz(psi_n(x, i) * psi_n(x, j), x) im = ax.imshow(overlap, cmap='RdBu', vmin=-0.5, vmax=1.5) ax.set_xlabel('$n$') ax.set_ylabel('$m$') ax.set_title('Orthonormality: $\\langle m|n \\rangle = \\delta_{mn}$') plt.colorbar(im, ax=ax) for i in range(n_max): for j in range(n_max): ax.text(j, i, f'{overlap[i,j]:.2f}', ha='center', va='center', fontsize=9) plt.suptitle('Quantum Harmonic Oscillator: Hilbert Space Structure', fontsize=14) plt.tight_layout() plt.savefig('quantum_synthesis.png', dpi=150) plt.show() # ============================================================ # Part 4: Week Summary Statistics # ============================================================ def week_summary(): """ Display week summary with key formulas. """ print("\n" + "=" * 60) print("WEEK 34 SUMMARY: BANACH AND HILBERT SPACES") print("=" * 60) summary = """ DAY 232: Normed Spaces and Banach Spaces ---------------------------------------- • Norm axioms: positivity, homogeneity, triangle inequality • Banach space = complete normed space • ℓ^p spaces (p ≥ 1) are Banach spaces • Hölder and Minkowski inequalities DAY 233: Inner Product Spaces ---------------------------- • Inner product axioms: conjugate symmetry, linearity, positive definiteness • Cauchy-Schwarz: |⟨x,y⟩| ≤ ||x|| ||y|| • Polarization identity: recover inner product from norm • Parallelogram law: characterizes inner product spaces DAY 234: Hilbert Spaces and L² ----------------------------- • Hilbert space = complete inner product space • L² = square-integrable functions (Lebesgue) • Riesz-Fischer: L² is complete • Wave functions ψ ∈ L² with ||ψ|| = 1 DAY 235: Orthonormal Sets and Gram-Schmidt ----------------------------------------- • Orthonormal: ⟨e_i, e_j⟩ = δ_ij • Gram-Schmidt: construct ONS from LI set • Fourier coefficients: c_n = ⟨e_n, x⟩ • Bessel's inequality: Σ|c_n|² ≤ ||x||² DAY 236: Orthonormal Bases and Parseval -------------------------------------- • ONB = complete orthonormal system • Parseval: Σ|c_n|² = ||x||² (equality for ONB) • Basel problem: Σ1/n² = π²/6 via Parseval • Completeness: Σ|e_n⟩⟨e_n| = I DAY 237: Separable Hilbert Spaces -------------------------------- • Separable = countable dense subset = countable ONB • Classification: all sep. inf-dim Hilbert spaces ≅ ℓ² • Riesz representation: f(x) = ⟨y, x⟩ for unique y • L² ≅ ℓ² (position ↔ energy representation) KEY QUANTUM CONNECTIONS: ----------------------- • State space = separable Hilbert space • ⟨φ|ψ⟩ = probability amplitude • |⟨φ|ψ⟩|² = transition probability • Σ|c_n|² = 1 = probability conservation """ print(summary) # ============================================================ # Main Execution # ============================================================ if __name__ == "__main__": print("=" * 60) print("Day 238: Week 34 Review - Computational Lab") print("=" * 60) print("\n1. Visualizing space hierarchy...") visualize_space_hierarchy() print("\n2. Verifying all inequalities...") verify_all_inequalities() print("\n3. Quantum mechanics synthesis...") quantum_synthesis() print("\n4. Week summary...") week_summary() print("\n" + "=" * 60) print("Week 34 Review Complete!") print("Next: Week 35 - Bounded Linear Operators") print("=" * 60) ``` --- ## 8. Looking Ahead: Week 35 ### Topics to Come **Bounded Linear Operators** - Definition and operator norm - Space $$B(\mathcal{H})$$ of bounded operators - Adjoint operators - Self-adjoint (Hermitian) operators - Unitary operators - Projection operators ### Why Operators Matter In quantum mechanics: - **Observables** are self-adjoint operators - **Symmetries** are unitary operators - **Measurement** involves projection operators - **Dynamics** follow $$i\hbar \frac{d|\psi\rangle}{dt} = H|\psi\rangle$$ where $$H$$ is an operator --- ## 9. Final Summary ### Week 34 Achievements This week, we built the complete mathematical foundation for quantum mechanics: 1. **Day 232**: Established norms, Banach spaces, and $$\ell^p$$ spaces 2. **Day 233**: Introduced inner products, Cauchy-Schwarz, parallelogram law 3. **Day 234**: Defined Hilbert spaces, proved $$L^2$$ is complete 4. **Day 235**: Developed Gram-Schmidt and Bessel's inequality 5. **Day 236**: Characterized orthonormal bases via Parseval's identity 6. **Day 237**: Classified separable Hilbert spaces as isomorphic to $$\ell^2$$ 7. **Day 238**: Synthesized everything into a coherent framework ### The Quantum Physics Takeaway $$\boxed{\text{Quantum mechanics = Linear algebra in infinite-dimensional Hilbert space}}$$ Every concept from finite-dimensional linear algebra (vectors, inner products, bases, projections) extends beautifully to Hilbert spaces, with completeness ensuring that limits and infinite series behave properly. --- *"The theory of Hilbert spaces is the natural setting for quantum mechanics. This week's material is not just mathematics—it is the language in which Nature speaks at the quantum level."*

Day 238: Week 34 Review - Banach and Hilbert Spaces

Today's Schedule (7 hours)

Day 238: Week 34 Review - Banach and Hilbert Spaces

Today's Schedule (7 hours)