Namaste! Chalo, aaj quantum mechanics ke ek dum mast topic, Rigid Rotor, ka poora post-mortem karte hain. Yeh B.Sc. wala basic nahi hai, yeh MSc/NET level ka detailed derivation hai. Hum har 'kyun' aur 'kaise' ko samjhenge, bilkul tutor style mein!
Sabse pehle, yeh 'Rigid Rotor' kya bala hai?
Ab MSc level ka twist: Yeh 'do atoms' (two-body problem) ko handle karna complicated hai. Hum isse simple 'one-body problem' mein convert karte hain.
Socho do dost (ek halka $m_1$, ek bhaari $m_2$) ek danda pakad ke ghoom rahe hain. Dono apne common 'center of mass' (C.M.) ke around ghoomenge.
Isse accha, hum ek naya, simpler system imagine karte hain:
"Ek 'imaginary' particle jiska mass $\mu$ (Reduced Mass) hai, woh C.M. se $r$ distance pe ghoom raha hai."
Yeh $\mu$ nikaalte kaise hain? Simple hai:
Moment of Inertia ($I$) hai: $I = m_1r_1^2 + m_2r_2^2$
Center of mass ke liye: $m_1r_1 = m_2r_2$ aur total distance $r = r_1 + r_2$.
Jab aap is maths ko solve karte ho, aapko milta hai:
$I = \left( \frac{m_1m_2}{m_1+m_2} \right) r^2$
Iss $ \frac{m_1m_2}{m_1+m_2} $ waali term ko hum kehte hain Reduced Mass ($\mu$).
Toh hamara final, simple formula ban gaya:
$$ I = \mu r^2 $$Lesson: Ab humein $m_1$ aur $m_2$ ko alag-alag track nahi karna. Hum bas ek particle (mass $\mu$) ko track karenge jo $r$ distance pe ghoom raha hai. Problem solved!
Ab humein is system ki total energy likhni hai. Quantum mein, hum energy ko Hamiltonian Operator ($\hat{H}$) kehte hain.
Total Energy = Kinetic Energy (ghoomne ki) + Potential Energy (position ki)
$\hat{H} = \hat{T} + \hat{V}$
Yahin hai sabse bada clue. Humne model ko kya bola? 'Rigid' Rotor.
Iska matlab distance $r$ fixed hai, badal nahi sakta. Potential energy tab aati hai jab distance change ho (jaise spring mein). Jab distance hi fix hai, toh koi potential energy nahi hai.
Toh, $V = 0$.
Analogy: Agar aap ek rubber band (spring) ko kheenchoge, usme potential energy store hogi. Par agar aap ek lohe ke dande ko ghumao, usme koi 'stretch' waali energy nahi aati. Hamara rotor lohe ka danda hai.
Toh saari energy sirf kinetic hai. Classical mechanics mein:
$E = \frac{1}{2}I\omega^2$
Humein Angular Momentum ($L = I\omega$) pasand hai, toh hum likhte hain:
$$ E = \frac{L^2}{2I} $$Quantum ka twist daalo (classical numbers ko operators bana do):
$$ \hat{H} = \frac{\hat{L}^2}{2I} $$Jahaan $\hat{L}^2$ hai angular momentum ka operator. Toh hamara final Hamiltonian hai:
$\hat{H} = \hat{T} + 0 = \frac{\hat{L}^2}{2I}$
Ab hum quantum ki master equation use karenge: The Time-Independent Schrödinger Equation.
$$ \hat{H}\psi = E\psi $$$\hat{H}$ ki value daal do jo humne upar nikali:
$$ \left( \frac{\hat{L}^2}{2I} \right) \psi = E\psi $$Thoda rearrange karo:
$$ \hat{L}^2\psi = (2IE)\psi $$Yeh hai hamari 'Problem Equation'. Agar humne $\psi$ (wavefunction) aur $E$ (energy) ke liye isse solve kar diya, toh humara kaam ho gaya!
Yahaan ek problem hai. $\hat{L}^2$ operator Cartesian coordinates (x, y, z) mein bohot messy dikhta hai.
Analogy: Agar aap ek merry-go-round pe ghoom rahe ho, toh aap apni position 'angle' (kitna degree ghoome) mein bataoge ya 'kitne feet left/right' mein? Obviously angle mein!
Waise hi, ghoomti hui cheez ke liye hum use karte hain Spherical Polar Coordinates ($r, \theta, \phi$).
Good News: Hamara rotor 'rigid' hai, yaani $r$ constant hai! Toh hamara wavefunction $\psi$ sirf $\theta$ aur $\phi$ pe depend karega.
$\psi = \psi(\theta, \phi)$
Jab hum $\hat{L}^2$ ko $r$ = constant maante hue spherical coordinates mein likhte hain, toh yeh aisa dikhta hai:
$$ \hat{L}^2 = -\hbar^2 \left[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2} \right] $$(Yahaan $\hbar = h/2\pi$, jise 'h-bar' kehte hain).
Ab iss poori cheez ko hamari Problem Equation mein daal do:
$$ -\hbar^2 \left[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial\psi}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2\psi}{\partial\phi^2} \right] = (2IE)\psi $$Yeh bohot daravna lag raha hai, hai na? Par isko solve karne ki ek standard trick hai.
Yeh ekdum 'Divide and Rule' policy hai. Hum maante hain ki hamara wavefunction $\psi$ do alag-alag, simpler functions ka product hai:
Analogy: Yeh 'Khichdi' $(\psi)$ ko solve karne jaisa hai. Hum 'Daal' ($\Theta$) aur 'Chawal' ($\Phi$) ko alag kar lete hain taaki unhe alag-alag handle kar sakein.
Jab hum yeh $\Theta\Phi$ ko uss badi equation mein daalte hain aur maths karte hain (LHS = $\theta$ part, RHS = $\phi$ part), toh equation do hisso mein toot jaati hai.
$\phi$ wala part bohot simple dikhta hai:
$$ \frac{d^2\Phi}{d\phi^2} + m^2\Phi = 0 $$(Yahaan $m^2$ ek 'separation constant' hai).
Yeh ek standard equation (Simple Harmonic) hai. Iska solution hai:
$\Phi(\phi) = A \cdot e^{im\phi}$
Ab aata hai Boundary Condition (Asli Quantum Twist):
Wavefunction $\psi$ ko 'single-valued' hona zaroori hai. Matlab, agar aap poora $360^\circ$ (ya $2\pi$ radians) ghoom ke waapis aao, toh aapko same value milni chahiye.
$\Phi(\phi) = \Phi(\phi + 2\pi)$
$A e^{im\phi} = A e^{im(\phi+2\pi)} = (A e^{im\phi}) \cdot (e^{im2\pi})$
Iska matlab $e^{im2\pi} = 1$ hona hi chahiye.
Euler's formula ($e^{ix} = \cos x + i \sin x$) se yeh tabhi possible hai jab:
$m = 0, \pm 1, \pm 2, \pm 3, \dots$
Yeh $m$ hamara pehla quantum number hai: Magnetic Quantum Number ($m_J$).
Ab bacha hua $\theta$ wala part aisa dikhta hai:
$$ \frac{1}{\sin\theta}\frac{d}{d\theta}\left(\sin\theta\frac{d\Theta}{d\theta}\right) + \left(\beta - \frac{m^2}{\sin^2\theta}\right)\Theta = 0 $$Yahaan $\beta = \frac{2IE}{\hbar^2}$ (bas ek naya naam diya simple karne ke liye).
Yeh equation $\sin\theta$ aur $\cos\theta$ mein hai, jo derivatives ke liye messy hain. Chalo ek naya, simpler variable introduce karte hain.
Substitution: Let $x = \cos\theta$
Toh $x$ ki range hogi $-1$ (south pole, $\theta=180^\circ$) se $+1$ (north pole, $\theta=0^\circ$).
Ab humein Chain Rule use karke $\frac{d}{d\theta}$ ko $\frac{d}{dx}$ mein badalna hai.
$\frac{dx}{d\theta} = -\sin\theta$
Toh, $\frac{d}{d\theta} = \frac{dx}{d\theta} \frac{d}{dx} = -\sin\theta \frac{d}{dx}$
Chalo equation ke pehle part ko convert karte hain:
$\frac{d\Theta}{d\theta} = -\sin\theta \frac{d\Theta}{dx}$
Toh $\sin\theta \frac{d\Theta}{d\theta} = \sin\theta \left(-\sin\theta \frac{d\Theta}{dx}\right) = -\sin^2\theta \frac{d\Theta}{dx}$
Kyunki $\sin^2\theta = 1 - \cos^2\theta = 1 - x^2$, yeh ban jaata hai:
$\sin\theta \frac{d\Theta}{d\theta} = -(1-x^2)\frac{d\Theta}{dx}$
Ab iska $\frac{d}{d\theta}$ lena hai... oh wait! Chain rule dobara!
$\frac{d}{d\theta} \left[ -(1-x^2)\frac{d\Theta}{dx} \right] = \left(-\sin\theta \frac{d}{dx}\right) \left[ -(1-x^2)\frac{d\Theta}{dx} \right]$
$\frac{d}{d\theta} [\dots] = \sin\theta \cdot \frac{d}{dx} \left[ (1-x^2)\frac{d\Theta}{dx} \right]$
Ab is poori cheez ko original equation ke $\frac{1}{\sin\theta}$ se multiply karo:
$\frac{1}{\sin\theta} \left[ \sin\theta \cdot \frac{d}{dx} \left[ (1-x^2)\frac{d\Theta}{dx} \right] \right] = \frac{d}{dx} \left[ (1-x^2)\frac{d\Theta}{dx} \right]$
Waah! Kitna clean ho gaya! Ab baaki parts ko substitute karte hain:
Hamari nayi equation bani:
$$ \frac{d}{dx} \left[ (1-x^2)\frac{dP(x)}{dx} \right] + \left(\beta - \frac{m^2}{1-x^2}\right)P(x) = 0 $$Yeh hai hamari famous Associated Legendre Differential Equation.
Isko solve karna thoda tricky hai. Hum 'Power Series Method' (Frobenius Method) use karte hain.
Case 1: Pehle $m=0$ (simplest case) solve karte hain.
Agar $m=0$, equation banti hai Legendre Equation:
$$ \frac{d}{dx} \left[ (1-x^2)\frac{dP(x)}{dx} \right] + \beta P(x) = 0 $$Product rule lagao:
$$ (1-x^2)\frac{d^2P}{dx^2} - 2x\frac{dP}{dx} + \beta P = 0 $$Hum assume karte hain ki solution ek infinite power series hai:
$P(x) = \sum_{k=0}^{\infty} a_k x^k$
Iske derivatives $P'$ aur $P''$ nikaalo aur equation mein daal do.
$P'(x) = \sum k a_k x^{k-1}$
$P''(x) = \sum k(k-1) a_k x^{k-2}$
Substitute karne par:
$(1-x^2)\sum k(k-1) a_k x^{k-2} - 2x\sum k a_k x^{k-1} + \beta \sum a_k x^k = 0$
$\sum k(k-1) a_k x^{k-2} - \sum k(k-1) a_k x^k - \sum 2k a_k x^k + \sum \beta a_k x^k = 0$
Ab humein $x$ ki saari powers ke coefficients ko ek saath lana hai. $x^k$ ke coefficient ko dekhte hain. Pehli term ($x^{k-2}$) ko $x^k$ banane ke liye, index $k$ ko $k+2$ se replace karo:
$\sum (k+2)(k+1) a_{k+2} x^k - \sum k(k-1) a_k x^k - \sum 2k a_k x^k + \sum \beta a_k x^k = 0$
Saare $x^k$ common le lo:
$\sum \left[ (k+2)(k+1)a_{k+2} - (k(k-1) + 2k - \beta)a_k \right] x^k = 0$
$\sum \left[ (k+2)(k+1)a_{k+2} - (k^2+k - \beta)a_k \right] x^k = 0$
Kyunki yeh har $x^k$ ke liye true hona chahiye, toh bracket ke andar ki poori cheez zero honi chahiye. Isse milta hai Recurrence Relation:
$$ (k+2)(k+1)a_{k+2} = (k(k+1) - \beta)a_k $$ $$ a_{k+2} = \frac{k(k+1) - \beta}{(k+2)(k+1)} a_k $$Yeh hamara 'recipe' hai. Agar $a_0$ pata hai, toh $a_2, a_4, \dots$ mil jaayenge. Agar $a_1$ pata hai, toh $a_3, a_5, \dots$ mil jaayenge.
The Quantum Twist (Boundary Condition)!
Hamara variable $x = \cos\theta$ hai. Iski value sirf $x = -1$ aur $x = +1$ (poles) ke beech ho sakti hai. Problem yeh hai ki yeh infinite series $x=\pm1$ par 'blow up' (diverge) kar jaati hai, matlab value infinity ho jaati hai.
Rule: Ek physical, 'well-behaved' wavefunction kabhi bhi infinity nahi ho sakta. Usse finite hona zaroori hai.
Solution: Hum is infinite series ko finite polynomial mein 'tod' (terminate) denge. Series terminate kaise hogi? Agar kisi point $k=J$ par, numerator zero ho jaaye, toh $a_{J+2}$ zero ho jaayega, aur uske baad $a_{J+4}, a_{J+6}$... sab zero ho jaayenge.
Toh, hum condition lagate hain:
$$ J(J+1) - \beta = 0 $$Yahaan se humein $\beta$ ki value mil gayi!
$$ \beta = J(J+1) \quad \text{jahaan } J=0, 1, 2, \dots $$Yeh jo finite polynomials bante hain, $P_J(x)$, inhein hi Legendre Polynomials kehte hain.
Case 2: Ab waapis $m \neq 0$ (Associated Legendre) par aate hain.
Poori equation yeh thi:
$$ \frac{d}{dx} \left[ (1-x^2)\frac{dP}{dx} \right] + \left(J(J+1) - \frac{m^2}{1-x^2}\right)P = 0 $$(Humne $\beta = J(J+1)$ daal diya).
Mathematicians ne bataya ki iska solution $P_J^m(x)$ (Associated Legendre) nikaalne ka ek 'jugaad' hai. Aap pehle simple $m=0$ wala $P_J(x)$ (Legendre Polynomial) lo aur usse $|m|$ baar differentiate karo, aur $(1-x^2)^{|m|/2}$ se multiply kar do.
$$ P_J^m(x) = (1-x^2)^{|m|/2} \frac{d^{|m|}}{dx^{|m|}} P_J(x) $$Yahaan ek aur condition aati hai. $P_J(x)$ ek $J$-degree polynomial hai (matlab $x^J$ uski highest power hai). Aap usse $J$ se zyada baar differentiate karoge toh answer zero aayega. Iska matlab, physically meaningful solution ke liye:
$|m| \le J$
Toh, $m$ ki allowed values hain:
$m = -J, -(J-1), \dots, 0, \dots, +(J-1), +J$
Toh humein dono quantum numbers aur unki conditions mil gayi, purely maths se!
Ab bas jod-tod karna hai. Humein $\beta$ ki do values pata hain:
Dono ko barabar rakh do:
$$ \frac{2IE}{\hbar^2} = J(J+1) $$Ab bas $E$ ke liye rearrange karo:
$$ E_J = \frac{\hbar^2}{2I} J(J+1) \quad \text{jahaan } J = 0, 1, 2, \dots $$
Yesss! Yeh hai hamare rigid rotor ke quantized energy levels.
Iska matlab kya hai?
Iska matlab yeh hai ki rigid rotor 'koi bhi' energy nahi le sakta. Woh sirf specific energy packets mein hi exist kar sakta hai, jo $J$ ki value pe depend karte hain.
Analogy (Energy ka Menu Card):
Aap hotel mein jaake yeh nahi bol sakte ki "mujhe 50 Rs. ki energy do". Menu fixed hai:
Aap sirf yahi 'dishes' (energy levels) order kar sakte ho.
Chemists (especially spectroscopists) ko Joules pasand nahi. Humein wavenumbers (cm⁻¹) pasand hain. Kyunki $E = hc\tilde{\nu}$.
Toh Energy (Joules) ko $hc$ se divide karke Wavenumber (cm⁻¹) banate hain:
$\tilde{E}_J = \frac{E_J}{hc} = \frac{1}{hc} \left[ \frac{\hbar^2}{2I} J(J+1) \right]$
Since $\hbar = h/2\pi$, $\hbar^2 = h^2/4\pi^2$
$\tilde{E}_J = \frac{1}{hc} \left[ \frac{h^2}{8\pi^2I} J(J+1) \right] = \left( \frac{h}{8\pi^2cI} \right) J(J+1)$
Yeh bracket waali poori term ($h, \pi, c, I$) ek specific molecule ke liye constant hai. Isse hum kehte hain Rotational Constant ($B$).
$$ B = \frac{h}{8\pi^2cI} \quad (\text{in cm}^{-1}) $$Toh hamara energy formula ban gaya super simple:
$\tilde{E}_J = B J(J+1)$
Ab menu card aur bhi simple hai: $E_0 = 0$, $E_1 = 2B$, $E_2 = 6B$, $E_3 = 12B$, ...
Hamari energy $E_J$ sirf $J$ pe depend karti hai. Usse $m$ se koi matlab nahi hai.
Par humne dekha ki har $J$ ke liye, $m$ ki kayi values possible hain (Chapter 4b ke solution se!):
$m = -J, -(J-1), \dots, 0, \dots, (J-1), +J$
Total kitni values hui? $J$ (negative) + $J$ (positive) + 1 (zero) = $(2J+1)$ values.
Iska matlab hai ki $(2J+1)$ alag-alag wavefunctions (states) hain jinki energy exactly same hai! Isse kehte hain Degeneracy.
Analogy (Hotel ka Floor):
Socho $J$ hai floor number aur $m$ hai room number.
- Ground Floor (J=0): $m=0$. Degeneracy = (2*0+1) = 1. Sirf ek room (Room 0) hai.
- First Floor (J=1): $m=-1, 0, +1$. Degeneracy = (2*1+1) = 3. Yahaan 3 rooms (Room -1, 0, 1) hain.
- Second Floor (J=2): $m=-2, -1, 0, +1, +2$. Degeneracy = (2*2+1) = 5. Yahaan 5 rooms hain.
Uss floor ke har room ka rent (Energy) same hai, par rooms (States/Wavefunctions) alag-alag hain.
Molecule light (photon) absorb karke $J=0$ se $J=1$ pe jump kar sakta hai. Par is 'party' mein entry ke liye kuch rules hain (Selection Rules):
Chalo dekhte hain absorption lines kahaan milengi. Molecule $J$ level se $J+1$ level pe jump karega.
Photon ki energy = Energy difference ($\Delta E$)
$\tilde{\nu}_{\text{line}} = \tilde{E}_{J+1} - \tilde{E}_{J}$
$\tilde{\nu}_{\text{line}} = [B(J+1)((J+1)+1)] - [B J(J+1)]$
$\tilde{\nu}_{\text{line}} = [B(J+1)(J+2)] - [B J(J+1)]$
$\tilde{\nu}_{\text{line}} = B(J+1) [ (J+2) - J ]$
$\tilde{\nu}_{\text{line}} = B(J+1) [ 2 ]$
$$ \tilde{\nu}_{\text{line}} = 2B(J+1) \quad (\text{jahaan } J=0, 1, 2, \dots) $$Yeh hai spectrum ka formula! Chalo lines plot karte hain (J yahaan *lower* state hai):
Result: Rotational spectrum mein humein equally spaced lines (barabar doori pe lines) milti hain, aur unke beech ka gap (spacing) hamesha $2B$ hota hai.
Mazaa aaya? Yeh hai Rigid Rotor ki poori kahani, zero se hero tak!
Time Remaining: 900s