Processing math: 68%
Skip to main content
LibreTexts - Ukrayinska

13.3: Виведення Лагранжа

Мета цієї глави полягає в тому, щоб знайти напругуV(r) і щільність зарядуρch(r) навколо атома, і ми будемо використовувати обчислення варіацій для виконання цього завдання. Нам потрібно зробити деякі досить серйозні припущення, щоб зробити цю проблему керованою. Розглянемо ізольований нейтральний атом з безліччю електронів навколо нього. ПрипустимоT0 K, тому всі електрони займають мінімально можливі енергетичні рівні. Припустимо, атом сферично симетричний. Всі величини, з якими ми стикаємося, такі як напруга, щільність заряду та Лагранж, змінюються,r але не змінюються зθ абоϕ. Ми будемо використовувати сферичні координати з початком у ядра атома. Хоча величини змінюються залежно від положення, припускайте, що кількість не змінюється з часом. Щільність зарядуρch(r) говорить нам, де електрони, швидше за все, в середньому можна знайти. Це пов'язано з квантовою механічною хвильовою функцієюψ,

ρch=q|ψ|2

деq - величина заряду електрона. Припустимо, що всі електрони, що оточують атом, розподілені рівномірно і можуть розглядатися так, ніби вони були однорідною електронною хмарою деякої щільності заряду.

Виберіть один з електронів атома і подумайте, що відбувається, коли електрон переміщується радіально всередину і назовні. Малюнок13.3.1 ілюструє цю ситуацію. У міру руху електрона відбувається перетворення енергії. Мета цього розділу - записати гамільтоніан і лагранж для цього процесу перетворення енергії. Запишемо ці величини в одиницях енергії на одиницю об'єму на розглянутий валентний електрон.

13.3.1.png
Малюнок13.3.1: Ілюстрація атома.

Щоб зрозуміти, що відбувається при переміщенні електрона, розглянемо енергію атома більш детально. Закон Кулона, введений в рівнянні 1.6.2, говорить нам про те, що заряджені об'єкти чинять сили на інші заряджені об'єкти. Більш конкретно, напруженість електричного поля,E обумовлена точковим зарядомQ кулонів, відстань,r оточена матеріалом з діелектричною проникністю,ϵ задається

E=Qˆar4πϵr2.

Атом складається зN позитивно заряджених протонів. Розглянутий електрон відчуває привабливу кулонівську силу завдяки цим протонам. Крім того, атом маєN електрони, іN1 з них надають відштовхувальну кулонівську силу на розглянутий електрон. Оскільки існує поділ заряду і електричне поле, накопичується енергія. Називають складову енергії атома за рахунок кулонівського взаємодії між протонами ядра і розглянутим електрономECoulombenucl. Викликати кулонівську взаємодію між розглянутим електроном і всіма іншими електронамиEeeinteract. Атом також володіє кінетичною енергією. Називають кінетичну енергію ядраEkineticnucl і кінетичну енергію всіх електронівEkinetice. Енергія атома - це сума всіх цих термінів.

Eatom=ECoulombenucl.+Ekineticnucl+Eeeinteract+Ekinetice

Енергія, обумовлена спином електронів і протонів, ігнорується, як і енергія внаслідок взаємодії з будь-якими іншими поруч зарядженими об'єктами. ПриT0 K кінетична енергія ядра буде близька до нуля, тому ми можемо ігнорувати термін,Ekineticnucl0. КількістьEkinetice не може бути рівно нулем. У главі 6 ми побудували діаграми енергетичного рівня для електронів навколо атома. Навіть приT=0 К електрони мають деяку внутрішню енергію, і ця енергія позначається займаним енергетичним рівнем.

Якщо у нас є великий атом з безліччю електронів навколо нього, кулонівське взаємодія між будь-яким одним електроном і ядром екранується кулонівським взаємодією від усіх інших електронів. Більш конкретно, припустимо, що у нас є ізольований атом зN протонами в ядрі таN електронами навколо нього. Якщо ми виберемо один з електронів,ECoulombenucl для цього електрон описує енергію, що зберігається в електричному полі за рахунок поділу заряду між ядром позитивного зарядуNq і цим електроном. Однак є йN1 інші електрони, які мають негативний заряд. ТермінEeeinteract описує енергію, накопичену в електричному полі за рахунок поділу заряду міжN1 іншими електронами і розглянутим електроном. Ці терміни дещо скасовують один одного, оскільки розглянутий електрон взаємодіє зN протонами кожного позитивного заряду q іN1 електронів кожного з негативного заряду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=mv. (Зверніть увагу, що імпульсM і узагальнений імпульсM різні і мають різні одиниці.)

Ekineticee=|M|22m

Ми не знаємо, як змінюється енергія як функція положення r Замість цього, ми можемо записати енергію як функцію імпульсу кристалаMcrystal або хвильового вектораk, і ми знаємо щось про зміну цих величин. Кристалічний імпульс дорівнює хвильовому вектору, масштабованому константою Планка.

Mcrystal=k

Він має одиниці імпульсуkgms, і він був введений в п. 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 spacefilled k levels(Ekineticee)(evolume)d(vol. all k space)

Кількість електронів на одиницю об'єму задається

(evolume)=ρchq.

Обсяг, зайнятий у зворотному просторі43πk3f, дорівнює об'єму сфери радіусаkf.

EkineticeV=143πk3ffilled k levels(2˜r22m)(ρchq)d(vol. all k space)

Диференціальний елемент об'єму виражається як

d3|k|=˜r2sin˜θd˜rd˜θd˜ϕ.

EkineticeV=143πk3ffilled 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. Інженери-електрики часто використовують слово напруга синонімом потенціалу. Коли хіміки використовують термін хімічний потенціал, вони мають на увазі ту саму величину, просто масштабовану константою. Подібно до того, як напруга є фундаментальною кількістю електротехніки, яка представляє, як важко переміщати електрони навколо, хімічний потенціал є фундаментальною кількістю хімії, яка представляє, наскільки важко переміщати електрони навколо.