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≤˜r≤kf.
EkineticeV=143πk3f⋅kf∫˜r=0π∫˜θ=02π∫˜ϕ=0(ℏ2˜r22m)(−ρchq)(˜r2sin˜θd˜rd˜θd˜ϕ)
Інтеграл вище може бути оцінений безпосередньо. Першим кроком для його оцінки є витягування констант назовні. Як описано вище,ρch варіюється з,r але ні˜r, так що його можна витягнути за межі інтеграла теж.
EkineticeV=−143πk3f⋅ℏ2ρch2mqkf∫˜r=0π∫˜θ=02π∫˜ϕ=0˜r4sin˜θd˜rd˜θd˜ϕ
Інтеграл відокремлюється і може бути оцінений.
EkineticeV=−143πk3f⋅ℏ2ρch2mq(π∫˜θ=02π∫˜ϕ=0sin˜θd˜rd˜θd˜ϕ)(kf∫˜r=0˜r4d˜r)
EkineticeV=−143πk3f⋅ℏ2ρch2mq4π(kf55)
EkineticeV=−3ρchk2fℏ210mq
Щільність заряду - це функція положення в реальному просторіr, і ми знаходимося в процесі вирішення цієї функціїρch(r). Однак це також залежить від енергії ФерміEf, а отже, і хвильового вектора Ферміkf, для атома. Далі знаходимо взаємозв'язок міжρch іkf. Два електрони допускаються на рівень енергії (обертатися вгору і обертатися вниз), отже, на заповненийk стан. Кількість заповнених станів на атом у зворотному просторі пов'язано з щільністю заряду.
ρch=−2q(no. filled k statesunit vol. in k space)
У п. 6.3.1 ми побачили, що примітивна клітина у взаємному просторі була(2π)3 часом примітивною клітиною в реальному просторі, тому
(unit vol. k space)=(2π)3⋅(unit vol. real space)=(2π)3.
Ми знаємо дещо про хвильові вектори заповнених станів у зворотному просторі. ПриT=0 К найнижчі стани заповнюються, а всі інші порожні, і вони заповнюються до радіусаkf. Обсяг сфери радіусаkf задається43πk3f, і це являє собою кількістьk заповнених станів на об'єм зворотного простору. Тому ми можемо спростити вираз вище.
ρch=−2q⋅43πk3f⋅1(2π)3
ρch=−q3π2k3f
kf=(−3π2qρch)1/3
Ми хочемо написатиEkineticeV як функцію узагальненого шляхуV. Тепер ми можемо досягти цього завдання, об'єднавши рівняння\ ref {13.3.31} і\ ref {13.3.36}.
EkineticeV=−3ℏ210mqρch(−3π2qρch)2/3
EkineticeV=−3ℏ210mq(−3π2q)2/3ρ5/3ch
Електрична енергія - це твір заряду і напруги. Більш конкретно, з Рівняння 2.2.7, він дається
E=12QV.
Щільність електричної енергії тоді задається
EV=12ρchV.
Використовуйте рівняння\ ref {13.3.40}, щоб зв'язатиρchV і.
EkineticeV=12ρchV=−3ℏ210mq(−3π2q)2/3ρ5/3ch
Тепер ми пов'язували узагальнений шлях і узагальнений потенціал.
V=−3ℏ25mq(−3π2q)2/3ρ2/3ch
ρch=(−5mq3ℏ2⋅(−3π2q)−2/3)3/2V3/2
ρch=[(−5mq3ℏ2)3/2(−q3π2)]⋅V3/2
Нарешті, ми можемо писатиEkineticeV як функціюV.
EkineticeV=[(−5mq3ℏ2)3/2(−q3π2)]V5/2
Зверніть увагу, що кількість в дужках вище постійна. Коефіцієнтc0 визначається з терміну в дужках.
c0=(−5mq3ℏ2)3/2(−q3π2)
EkineticeV=c0V5/2
Тепер ми можемо описати всі терміни Лагранжа з точки зору нашого узагальненого шляху.
ECoulombenuclV+EeeinteractV=12ϵ|→∇V|2
EkineticeV=c0V5/2
Гамільтоніан являє собою загальну щільність енергії, а Лагранж являє собою різницю енергетичної щільності цих форм енергії. Гамільтоніан і Лагранж мають формуH=H(r,V,dVdr) іL=L(r,V,dVdr) деr знаходиться положення в сферичних координатах. Немаєθ абоϕ залежність відH абоL. Все сферично симетрично.
H=12ϵ|→∇V|2+c0V5/2
L=12ϵ|→∇V|2−c0V5/2
В якості осторонь розглянемо щеEf=μchem раз енергію Фермі. За допомогою деякої алгебри ми можемо записати її як функцію напруги. Використовуйте рівняння\ ref {13.3.21},\ ref {13.3.36} та\ ref {13.3.46}.
Ef=ℏ2k2f2m=ℏ22m(−3π2ρchq)2/3
Ef=ℏ22m(−3π2q)2/3[(−5mq3ℏ2⋅(−3π2q)−2/3)3/2V3/2]2/3
Ef=−5q6V
Зверніть увагу, що енергія Фермі - це лише масштабована версія напругиV щодо рівня землі наr=∞. Інженери-електрики часто використовують слово напруга синонімом потенціалу. Коли хіміки використовують термін хімічний потенціал, вони мають на увазі ту саму величину, просто масштабовану константою. Подібно до того, як напруга є фундаментальною кількістю електротехніки, яка представляє, як важко переміщати електрони навколо, хімічний потенціал є фундаментальною кількістю хімії, яка представляє, наскільки важко переміщати електрони навколо.