13.3: Виведення Лагранжа
Мета цієї глави полягає в тому, щоб знайти напругуV(r) і щільність зарядуρch(r) навколо атома, і ми будемо використовувати обчислення варіацій для виконання цього завдання. Нам потрібно зробити деякі досить серйозні припущення, щоб зробити цю проблему керованою. Розглянемо ізольований нейтральний атом з безліччю електронів навколо нього. ПрипустимоT≈0 K, тому всі електрони займають мінімально можливі енергетичні рівні. Припустимо, атом сферично симетричний. Всі величини, з якими ми стикаємося, такі як напруга, щільність заряду та Лагранж, змінюються,r але не змінюються зθ абоϕ. Ми будемо використовувати сферичні координати з початком у ядра атома. Хоча величини змінюються залежно від положення, припускайте, що кількість не змінюється з часом. Щільність зарядуρch(r) говорить нам, де електрони, швидше за все, в середньому можна знайти. Це пов'язано з квантовою механічною хвильовою функцієюψ,
ρch=−q⋅|ψ|2
деq - величина заряду електрона. Припустимо, що всі електрони, що оточують атом, розподілені рівномірно і можуть розглядатися так, ніби вони були однорідною електронною хмарою деякої щільності заряду.
Виберіть один з електронів атома і подумайте, що відбувається, коли електрон переміщується радіально всередину і назовні. Малюнок13.3.1 ілюструє цю ситуацію. У міру руху електрона відбувається перетворення енергії. Мета цього розділу - записати гамільтоніан і лагранж для цього процесу перетворення енергії. Запишемо ці величини в одиницях енергії на одиницю об'єму на розглянутий валентний електрон.

Щоб зрозуміти, що відбувається при переміщенні електрона, розглянемо енергію атома більш детально. Закон Кулона, введений в рівнянні 1.6.2, говорить нам про те, що заряджені об'єкти чинять сили на інші заряджені об'єкти. Більш конкретно, напруженість електричного поля,→E обумовлена точковим зарядомQ кулонів, відстань,r оточена матеріалом з діелектричною проникністю,ϵ задається
→E=Qˆar4πϵr2.
Атом складається зN позитивно заряджених протонів. Розглянутий електрон відчуває привабливу кулонівську силу завдяки цим протонам. Крім того, атом маєN електрони, іN−1 з них надають відштовхувальну кулонівську силу на розглянутий електрон. Оскільки існує поділ заряду і електричне поле, накопичується енергія. Називають складову енергії атома за рахунок кулонівського взаємодії між протонами ядра і розглянутим електрономECoulombenucl. Викликати кулонівську взаємодію між розглянутим електроном і всіма іншими електронамиEeeinteract. Атом також володіє кінетичною енергією. Називають кінетичну енергію ядраEkineticnucl і кінетичну енергію всіх електронівEkinetice. Енергія атома - це сума всіх цих термінів.
Eatom=ECoulombenucl.+Ekineticnucl+Eeeinteract+Ekinetice
Енергія, обумовлена спином електронів і протонів, ігнорується, як і енергія внаслідок взаємодії з будь-якими іншими поруч зарядженими об'єктами. ПриT≈0 K кінетична енергія ядра буде близька до нуля, тому ми можемо ігнорувати термін,Ekineticnucl≈0. КількістьEkinetice не може бути рівно нулем. У главі 6 ми побудували діаграми енергетичного рівня для електронів навколо атома. Навіть приT=0 К електрони мають деяку внутрішню енергію, і ця енергія позначається займаним енергетичним рівнем.
Якщо у нас є великий атом з безліччю електронів навколо нього, кулонівське взаємодія між будь-яким одним електроном і ядром екранується кулонівським взаємодією від усіх інших електронів. Більш конкретно, припустимо, що у нас є ізольований атом зN протонами в ядрі таN електронами навколо нього. Якщо ми виберемо один з електронів,ECoulombenucl для цього електрон описує енергію, що зберігається в електричному полі за рахунок поділу заряду між ядром позитивного зарядуNq і цим електроном. Однак є йN−1 інші електрони, які мають негативний заряд. ТермінEeeinteract описує енергію, накопичену в електричному полі за рахунок поділу заряду міжN−1 іншими електронами і розглянутим електроном. Ці терміни дещо скасовують один одного, оскільки розглянутий електрон взаємодіє зN протонами кожного позитивного заряду q іN−1 електронів кожного з негативного заряду−q. Однак терміни не йдуть повністю. Обчислення
ECoulombenucl+Eeeinteract
ускладнюється, тому що електрони знаходяться в русі, і ми насправді не знаємо, де вони знаходяться або навіть де їх найімовірніше знайти. Фактично, ми намагаємося вирішити, де вони, ймовірно, будуть знайдені.
Коли ми рухаємо розглянутий електрон всередину і назовні радіально, енергія передається між (ECoulombenucl+Eeeinteract) іEkinetice. Гамільтоніан - це сума цих двох форм енергії на одиницю об'єму, а Лагранж - різниця цих двох форм енергії на одиницю об'єму. Обидві величини мають одиниціJm3. Виберіть напругуV(r) як узагальнений шлях, а щільність зарядуρch(r) як узагальнений потенціал. Незалежною змінною цих величин є радіальне положенняr, а не час. Тепер ми можемо написати гамільтоніан і лагранж.
H(r,V,dVdr)=(ECoulombenuclV+EeeinteractV)+EkineticeV
L(r,V,dVdr)=(ECoulombenuclV+EeeinteractV)−EkineticeV
Наступний крок - написати
ECoulombenuclV+EeeinteractV
в плані шляхуV. Як детально описано в таблиці 12.2.3, щільність енергії, обумовлена електричним полем→E, задається
EV=12ϵ|→E|2.
Пам'ятайте, щоE представляє енергію, а→E представляє електричне поле. Електричне поле - негативний градієнт напругиV(r).
→E=−→∇V.
Ми можемо об'єднати ці вирази та рівняння 13.2.6, щоб записати перший член Гамільтоніана та Лагранжа через узагальнений шлях.
ECoulombenuclV+EeeinteractV=12ϵ|→∇V|2
H(r,V,dVdr)=(12ϵ|→∇V|2)+EkineticeV
L(r,V,dVdr)=(12ϵ|→∇V|2)−EkineticeV
Наступне завдання - описати термін, що залишивсяEkineticeV, як функцію узагальненого шляху теж. Це завдання трохи складніше. Ми продовжуємо застосовувати підхід до суворих наближень, поки це не буде керованим. Нам потрібно висловитиρch(r) як функціюV(r). Тоді з деякою алгеброю,EkineticeV може бути записана чисто як функціяV(r).
Ми хочемо узагальнити про кінетичну енергію електронів. Однак кожен електрон має свою швидкість→v і імпульс→M. Ці величини залежать від положення
→r=rˆar+θˆaθ+ϕˆaϕ
якимось невідомим способом. Крім того, розрахунокEkineticeV залежить від щільності зарядуρch(r), яка є невідомою величиною, яку ми намагаємося знайти. Нам більше пощастило, описуючи ці величини у взаємному просторі, введеному в п. 6.3. Положення позначається у зворотному просторі хвильовим вектором
→k=˜rˆar+˜θˆaθ+˜ϕˆaϕ
Ми можемо описати властивості матеріалу, описуючи, як вони змінюються залежно від положення в реальному просторі. Наприклад,ρch(r) представляє щільність заряду електронів як функцію відстані r від центру атома. Нам може бути цікаво, як інші величини, такі як енергія, необхідна для зривання електрона, або кінетична енергія, внутрішня для електрона, змінюються залежно від положення в реальному просторі. Замість того, щоб описувати, як величини змінюються залежно від положення в реальному просторі, ми можемо описати, як величини змінюються з просторовою частотою електронів. Це ідея, що стоїть за представленням величин у взаємному просторі. Нас може зацікавити, як змінюється щільність заряду електронів в залежності від просторової частоти зарядів у кристалі або іншому матеріалі, і це ідея, представлена функціями хвильового вектора, такими якρch(→k). Намагаємося вирішити для щільності зарядуρch(r). Ми очікуємо, що електрони частіше виявляються на певних відстаняхr від центру атома, ніж на інших відстанях. Однак немає закономірності до щільності заряду як функції хвильового вектора,ρch(→k). Припустимо, щоρch це приблизно постійна щодо|→k| до якогось рівня. При деякій більшій роботі це припущення дозволить нам вирішити для щільності зарядуρch(r).
Кінетична енергія одного електрона задається
Ekineticee−=12m|→v|2
деm - маса електрона. Ми можемо записати цю енергію з точки зору імпульсу,→M=m→v. (Зверніть увагу, що імпульс→M і узагальнений імпульсM різні і мають різні одиниці.)
Ekineticee−=|→M|22m
Ми не знаємо, як змінюється енергія як функція положення r Замість цього, ми можемо записати енергію як функцію імпульсу кристала→Mcrystal або хвильового вектора→k, і ми знаємо щось про зміну цих величин. Кристалічний імпульс дорівнює хвильовому вектору, масштабованому константою Планка.
→Mcrystal=ℏ→k
Він має одиниці імпульсуkg⋅ms, і він був введений в п. 6.3.2. Кінетична енергія одного електрона як функція імпульсу кристала задається
Ekineticee−=(→Mcrystal)22m=(ℏ|→k|)22m.
Вектор у зворотному просторі представлений Equation\ ref {13.3.14}, а Equation\ ref {13.3.18} можна спростити, оскільки ми передбачаємо сферичну симетрію˜θ=˜ϕ=0. Величина хвильового вектора стає|→k|=˜r, і ми можемо записати енергію як
Ekineticee−=ℏ2˜r22m.
Так само, як кожен електрон має свій імпульсm|→v|, кожен електрон має свій кришталевий імпульсℏ|→k|. Однак ми знаємо деяку інформацію про хвильовому|→k| векторі електронів в атомі. ПриT=0 К електрони займають найнижчі дозволені енергетичні стани. Енергетичні держави зайняті до якоїсь найвищої зайнятої держави, яка називається енергією ФерміEf. Хоча інженери-електрики використовують термін енергія Фермі, хіміки іноді використовують термін хімічний потенціалμchem. Найнижчі енергетичні стани, зайняті в той час як вищі порожні. Аналогічно хвильові вектори зайняті до деякого найвищого зайнятого хвильового вектора, який називається хвильовим вектором Ферміkf.
|→k|={filled state˜r<kfempty state˜r>kf
Енергія Фермі і хвильовий вектор Фермі пов'язані
Ef=ℏ2k2f2m.
Ми використовуємо ідею зворотного простору для написання виразу для кінетичної енергії електронів на одиницю об'єму [136, с. 49]. Кінетична енергія, обумовлена будь-яким електроном як функція положення у зворотному просторі, задається Equation\ ref {13.3.19}. Відзначимо, що при кожному|→k|=˜r значенні електрон має різну кінетичну енергію. Щоб знайти кінетичну енергію на одиницю об'єму за рахунок усіх електронів, ми інтегруємо по всьому|→k|=˜r в сферичні координати, які зайняті електронами, а потім ділимо на об'єм, зайнятий у→k просторі.
EkineticeV=1vol. occupied in k space⋅∫filled k levels(Ekineticee−)(e−volume)d(vol. all k space)
Кількість електронів на одиницю об'єму задається
(e−volume)=−ρchq.
Обсяг, зайнятий у зворотному просторі43πk3f, дорівнює об'єму сфери радіусаkf.
EkineticeV=143πk3f⋅∫filled k levels(ℏ2˜r22m)(−ρchq)d(vol. all k space)
Диференціальний елемент об'єму виражається як
d3|→k|=˜r2sin˜θd˜rd˜θd˜ϕ.
EkineticeV=143πk3f⋅∫filled k levels(ℏ2|→k|22m)(−ρchq)(˜r2sin˜θd˜rd˜θd˜ϕ)
Як описано вище, електрони займають стани у зворотному просторі тільки с0 \leq \tilde r \leq k_f.
\frac{E_{kinetic\, e}}{\mathbb{V}} = \frac{1}{\frac{4}{3}\pi k^3_f} \cdot \int\limits_{\tilde r = 0}^{k_f} \,\int\limits_{\tilde \theta = 0}^{\pi} \,\int\limits_{\tilde \phi = 0}^{2\pi} \left(\frac{\hbar^2 \tilde r^2}{2m}\right)\left(\frac{-\rho_{ch}}{q}\right) (\tilde r^2 \sin \tilde \theta\, d \tilde r \,d \tilde \theta \,d \tilde \phi) \nonumber
Інтеграл вище може бути оцінений безпосередньо. Першим кроком для його оцінки є витягування констант назовні. Як описано вище,\rho_{ch} варіюється з,r але ні\tilde r, так що його можна витягнути за межі інтеграла теж.
\frac{E_{kinetic\, e}}{\mathbb{V}} = \frac{-1}{\frac{4}{3}\pi k^3_f} \cdot \frac{\hbar^2 \rho_{ch}}{2mq} \int\limits_{\tilde r = 0}^{k_f} \,\int\limits_{\tilde \theta = 0}^{\pi} \,\int\limits_{\tilde \phi = 0}^{2\pi} \tilde r^4 \sin \tilde \theta\, d \tilde r \,d \tilde \theta \,d \tilde \phi \nonumber
Інтеграл відокремлюється і може бути оцінений.
\frac{E_{kinetic\, e}}{\mathbb{V}} = \frac{-1}{\frac{4}{3}\pi k^3_f} \cdot \frac{\hbar^2 \rho_{ch}}{2mq} \left( \int\limits_{\tilde \theta = 0}^{\pi} \,\int\limits_{\tilde \phi = 0}^{2\pi} \sin \tilde \theta\, d \tilde r \,d \tilde \theta \,d \tilde \phi \right) \left( \int\limits_{\tilde r = 0}^{k_f} \tilde r^4 d \tilde r \right) \nonumber
\frac{E_{kinetic \, e}}{\mathbb{V}} = \frac{-1}{\frac{4}{3} \pi k^3_f} \cdot \frac{\hbar^2 \rho_{ch}}{2mq}4\pi \left(\frac{{k_f}^5}{5}\right) \nonumber
\frac{E_{kinetic \, e}}{\mathbb{V}} =\frac{-3 \rho_{ch} k^2_f \hbar^2}{10mq} \label{13.3.31}
Щільність заряду - це функція положення в реальному просторіr, і ми знаходимося в процесі вирішення цієї функції\rho_{ch}(r). Однак це також залежить від енергії ФерміE_f, а отже, і хвильового вектора Ферміk_f, для атома. Далі знаходимо взаємозв'язок між\rho_{ch} іk_f. Два електрони допускаються на рівень енергії (обертатися вгору і обертатися вниз), отже, на заповненийk стан. Кількість заповнених станів на атом у зворотному просторі пов'язано з щільністю заряду.
\rho_{ch} = -2q \left( \frac{\text{no. filled } k \text{ states}}{\text{unit vol. in } k \text{ space}}\right) \nonumber
У п. 6.3.1 ми побачили, що примітивна клітина у взаємному просторі була(2 \pi )^3 часом примітивною клітиною в реальному просторі, тому
(\text{unit vol. } k \text{ space}) = (2 \pi )^3 \cdot (\text{unit vol. real space}) = (2 \pi )^3. \nonumber
Ми знаємо дещо про хвильові вектори заповнених станів у зворотному просторі. ПриT = 0 К найнижчі стани заповнюються, а всі інші порожні, і вони заповнюються до радіусаk_f. Обсяг сфери радіусаk_f задається\frac{4}{3}\pi k^3_f, і це являє собою кількістьk заповнених станів на об'єм зворотного простору. Тому ми можемо спростити вираз вище.
\rho_{ch} = -2q \cdot \frac{4}{3} \pi k_f^3 \cdot \frac{1}{(2 \pi )^3} \nonumber
\rho_{ch} = \frac{-q}{3 \pi^2}k_f^3 \nonumber
k_f = \left(\frac{-3 \pi^2}{q}\rho_{ch}\right)^{1/3} \label{13.3.36}
Ми хочемо написати\frac{E_{kinetic\, e}}{\mathbb{V}} як функцію узагальненого шляхуV. Тепер ми можемо досягти цього завдання, об'єднавши рівняння\ ref {13.3.31} і\ ref {13.3.36}.
\frac{E_{kinetic \, e}}{\mathbb{V}}= \frac{-3\hbar^2}{10mq}\rho_{ch} \left(\frac{-3\pi^2}{q}\rho_{ch}\right)^{2/3} \nonumber
\frac{E_{kinetic \, e}}{\mathbb{V}}= \frac{-3\hbar^2}{10mq} \left(\frac{-3\pi^2}{q}\right)^{2/3} \rho_{ch}^{5/3} \nonumber
Електрична енергія - це твір заряду і напруги. Більш конкретно, з Рівняння 2.2.7, він дається
E = \frac{1}{2}QV. \nonumber
Щільність електричної енергії тоді задається
\frac{E}{\mathbb{V}} = \frac{1}{2}\rho_{ch}V. \label{13.3.40}
Використовуйте рівняння\ ref {13.3.40}, щоб зв'язати\rho_{ch}V і.
\frac{E_{kinetic \, e}}{\mathbb{V}}= \frac{1}{2}\rho_{ch}V =\frac{-3\hbar^2}{10mq} \left(\frac{-3\pi^2}{q}\right)^{2/3} \rho_{ch}^{5/3} \nonumber
Тепер ми пов'язували узагальнений шлях і узагальнений потенціал.
V = \frac{-3\hbar^2}{5mq} \left(\frac{-3\pi^2}{q}\right)^{2/3} \rho_{ch}^{2/3} \nonumber
\rho_{ch} = \left(\frac{-5mq}{3\hbar^2} \cdot \left(\frac{-3\pi^2}{q}\right)^{-2/3} \right)^{3/2} V^{3/2} \nonumber
\rho_{ch} = \left[ \left(\frac{-5mq}{3\hbar^2}\right)^{3/2} \left(\frac{-q}{3\pi^2}\right)\right] \cdot V^{3/2} \nonumber
Нарешті, ми можемо писати\frac{E_{kinetic\, e}}{\mathbb{V}} як функціюV.
\frac{E_{kinetic \, e}}{\mathbb{V}}= \left[ \left(\frac{-5mq}{3\hbar^2}\right)^{3/2} \left(\frac{-q}{3\pi^2}\right)\right]V^{5/2} \nonumber
Зверніть увагу, що кількість в дужках вище постійна. Коефіцієнтc_0 визначається з терміну в дужках.
c_0 = \left(\frac{-5mq}{3\hbar^2}\right)^{3/2} \left(\frac{-q}{3\pi^2}\right) \label{13.3.46}
\frac{E_{kinetic \, e}}{\mathbb{V}}=c_{0} V^{5 / 2} \nonumber
Тепер ми можемо описати всі терміни Лагранжа з точки зору нашого узагальненого шляху.
\frac{E_{Coulomb\, e \,nucl}}{\mathbb{V}} + \frac{E_{e\, e \,interact}}{\mathbb{V}} = \frac{1}{2}\epsilon |\overrightarrow{\nabla}V|^2 \nonumber
\frac{E_{kinetic \, e}}{\mathbb{V}} = c_{0} V^{5 / 2} \nonumber
Гамільтоніан являє собою загальну щільність енергії, а Лагранж являє собою різницю енергетичної щільності цих форм енергії. Гамільтоніан і Лагранж мають формуH = H (r, V, \frac{dV}{dr}) і\mathcal{L} = \mathcal{L} (r, V, \frac{dV}{dr}) деr знаходиться положення в сферичних координатах. Немає\theta або\phi залежність відH або\mathcal{L}. Все сферично симетрично.
H = \frac{1}{2}\epsilon |\overrightarrow{\nabla}V|^2 + c_{0} V^{5 / 2} \nonumber
\mathcal{L} = \frac{1}{2}\epsilon |\overrightarrow{\nabla}V|^2 - c_{0} V^{5 / 2} \nonumber
В якості осторонь розглянемо щеE_f = \mu_{chem} раз енергію Фермі. За допомогою деякої алгебри ми можемо записати її як функцію напруги. Використовуйте рівняння\ ref {13.3.21},\ ref {13.3.36} та\ ref {13.3.46}.
E_f= \frac{\hbar^2k^2_f}{2m} = \frac{\hbar^2}{2m} \left( \frac{-3\pi^2 \rho_{ch}}{q}\right)^{2/3} \nonumber
E_f= \frac{\hbar^2}{2m} \left( \frac{-3\pi^2}{q}\right)^{2/3} \left[ \left(\frac{-5mq}{3\hbar^2} \cdot \left( \frac{-3\pi^2}{q}\right)^{-2/3}\right)^{3/2} V^{3/2}\right]^{2/3} \nonumber
E_f = \frac{-5q}{6}V \nonumber
Зверніть увагу, що енергія Фермі - це лише масштабована версія напругиV щодо рівня землі наr = \infty. Інженери-електрики часто використовують слово напруга синонімом потенціалу. Коли хіміки використовують термін хімічний потенціал, вони мають на увазі ту саму величину, просто масштабовану константою. Подібно до того, як напруга є фундаментальною кількістю електротехніки, яка представляє, як важко переміщати електрони навколо, хімічний потенціал є фундаментальною кількістю хімії, яка представляє, наскільки важко переміщати електрони навколо.