Метод Эйнштейна – Бриллюэна – Келлера

редактировать

Метод Эйнштейна – Бриллюэна – Келлера (EBK ) - это полуклассический метод (названный в честь Альберт Эйнштейн, Леон Бриллюэн и Джозеф Б. Келлер ), используемые для вычисления собственных значений в квантово-механических системах. Квантование EBK является улучшением по сравнению с квантованием Бора-Зоммерфельда, которое не учитывает скачки фазы каустики в классических поворотных точках. Эта процедура способна точно воспроизвести спектр трехмерного гармонического осциллятора, частицы в коробке и даже релятивистскую тонкую структуру водорода . atom.

В 1976–1977 гг. Берри и Табор вывели расширение для Gutzwiller формулы следа для плотности состояний интегрируемой системы, начиная с квантования EBK.

В последнее время был получен ряд результатов по вычислительным вопросам, связанным с этой темой, например, работа Эрик Дж. Хеллер и Эммануэль Дэвид Танненбаум с использованием метода градиентного спуска по уравнению в частных производных.

Содержание
  • 1 Процедура
    • 1.1 Пример: 2D атом водорода
  • 2 См. также
  • 3 Ссылки
Процедура

Для разделимой классической системы, определенной координатами (qi, pi); я ∈ {1, 2, ⋯, d} {\ displaystyle (q_ {i}, p_ {i}); я \ in \ {1,2, \ cdots, d \}}{\ displaystyle (q_ {i}, p_ {i}); i \ in \ {1,2, \ cdots, d \}} , дюйм каждая пара (qi, pi) {\ displaystyle (q_ {i}, p_ {i})}(q_i, p_i) описывает замкнутую функцию или периодическую функцию в qi {\ displaystyle q_ { i}}q_ {i} , процедура EBK включает квантование интегралов по путям пи {\ displaystyle p_ {i}}p_ {i} по замкнутой орбите qi {\ displaystyle q_ {i}}q_ {i} :

I i = 1 2 π ∮ pidqi = ℏ (ni + μ i 4 + bi 2) {\ displaystyle I_ {i} = {\ frac {1} {2 \ pi}} \ oint p_ {i} dq_ {i} = \ hbar \ left (n_ {i} + {\ frac {\ mu _ {i}} {4}} + {\ frac {b_ {i}} {2}} \ right) }{\ displaystyle I_ {i} = {\ frac {1} {2 \ pi}} \ oint p_ {i} dq_ {i} = \ hbar \ left (n_ {i} + {\ frac {\ mu _ {i}} {4}} + {\ frac {b_ {i}} {2}} \ right)}

где I i {\ displaystyle I_ {i}}I_ {i} - это координата угла действия, ni {\ displaystyle n_ {i}}n_ {i} - положительное целое число, а μ i {\ displaystyle \ mu _ {i}}\ mu _ {i} и bi {\ displaystyle b_ {i}}b_ {i} - индексы Маслова. μ i {\ displaystyle \ mu _ {i}}\ mu _ {i} соответствует количеству классических поворотных точек на траектории qi {\ displaystyle q_ {i}}q_ {i} (границы Дирихле условие ), а bi {\ displaystyle b_ {i}}b_ {i} соответствует количеству отражений от твердой стенки (граничное условие Неймана ).

Пример: 2D атом водорода

Гамильтониан для нерелятивистского электрона (электрический заряд e {\ displaystyle e}e ) в атоме водорода:

H = pr 2 2 m + p φ 2 2 мр 2 - е 2 4 π ϵ 0 р {\ displaystyle H = {\ frac {p_ {r} ^ {2}} {2m}} + {\ frac {p _ {\ varphi} ^ {2}} {2mr ^ {2}}} - {\ frac {e ^ {2}} {4 \ pi \ epsilon _ {0} r}}}{\ displaystyle H = { \ frac {p_ {r} ^ {2}} {2m}} + {\ frac {p _ {\ varphi} ^ {2}} {2mr ^ {2}}} - {\ frac {e ^ {2}} {4 \ pi \ epsilon _ {0} r}}}

где pr {\ displaystyle p_ {r}}p_ {r} - канонический импульс к радиальному расстоянию r {\ displaystyle r}r , и p φ {\ displaystyle p _ {\ varphi}}{\ displaystyle п _ {\ varphi}} - канонический импульс азимутального угла φ {\ displaystyle \ varphi}\ varphi . Возьмите координаты угла действия:

I φ = constant = L {\ displaystyle I _ {\ varphi} = {\ text {constant}} = L}{\ displaystyle I _ {\ varphi} = {\ text {constant}} = L}

Для радиальной координаты r {\ displaystyle r}r :

pr = 2 m E - L 2 r 2 - e 2 4 π ϵ 0 р {\ displaystyle p_ {r} = {\ sqrt {2mE - {\ frac {L ^ {2}} {r ^ {2}}} - {\ frac {e ^ {2}} {4 \ pi \ epsilon _ {0} r}}}}}{\ displaystyle p_ {r} = {\ sqrt {2mE - {\ frac { L ^ {2}} {r ^ {2}}} - {\ frac {e ^ {2}} {4 \ pi \ epsilon _ {0} r}}}}}
I r = 1 π ∫ r 1 r 2 prdr = me 2 4 π ϵ 0 - 2 m E - | L | {\ displaystyle I_ {r} = {\ frac {1} {\ pi}} \ int _ {r_ {1}} ^ {r_ {2}} p_ {r} dr = {\ frac {me ^ {2} } {4 \ pi \ epsilon _ {0} {\ sqrt {-2mE}}}} - | L |}{\ displaystyle I_ {r} = {\ frac {1} {\ pi}} \ int _ {r_ {1}} ^ {r_ {2}} p_ {r} dr = {\ frac {me ^ {2}} {4 \ pi \ epsilon _ {0} {\ sqrt {-2mE}}}} - | L |}

, где мы интегрируем две классические точки поворота r 1, r 2 {\ displaystyle r_ {1}, r_ {2}}{\ displaystyle r_ {1}, r_ {2}} (μ r = 2 {\ displaystyle \ mu _ {r} = 2}{\ displaystyle \ mu _ {r} = 2} )

E = - me 4 32 π 2 ϵ 0 2 (I r + L) 2 {\ displaystyle E = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ epsilon _ {0} ^ {2} (I_ {r} + L) ^ {2}}}}{\ displaystyle E = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ epsilon _ {0} ^ {2} ( I_ {r} + L) ^ {2}}}}

Использование квантования EBK br = μ φ = b φ = 0, n φ = m {\ displaystyle b_ {r} = \ mu _ {\ varphi} = b _ {\ varphi} = 0, n _ {\ varphi } = m}{\ displaystyle b_ {r} = \ mu _ {\ varphi} = b _ {\ varphi} = 0, n _ {\ varphi} = m} :

L = ℏ м; м знак равно 0, 1, 2, ⋯ {\ displaystyle L = \ hbar m \ quad; \ quad m = 0,1,2, \ cdots}{\ displaystyle L = \ hbar m \ quad; \ quad m = 0,1,2, \ cdots}
I r = ℏ (n r + 1/2); nr знак равно 0, 1, 2, ⋯ {\ displaystyle I_ {r} = \ hbar (n_ {r} +1/2) \ quad; \ quad n_ {r} = 0,1,2, \ cdots}{\ displaystyle I_ {r} = \ hbar (n_ {r} +1/2) \ quad; \ quad n_ {r} = 0,1,2, \ cdots}
E = - меня 4 32 π 2 ϵ 0 2 ℏ 2 (nr + m + 1/2) 2 {\ displaystyle E = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ epsilon _ {0} ^ {2} \ hbar ^ {2} (n_ {r} + m + 1/2) ^ {2}}}}{\ displaystyle E = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ эпсилон _ {0} ^ {2} \ hbar ^ {2} (n_ {r} + m + 1/2) ^ {2}}}}

и сделав n = nr + m + 1 { \ displaystyle n = n_ {r} + m + 1}{\ displaystyle n = n_ {r} + m + 1} восстанавливается спектр двумерного атома водорода:

E n = - me 4 32 π 2 ϵ 0 2 ℏ 2 (n - 1 / 2) 2; п = 1, 2, 3, ⋯ {\ displaystyle E_ {n} = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ epsilon _ {0} ^ {2} \ hbar ^ { 2} (n-1/2) ^ {2}}} \ quad; \ quad n = 1,2,3, \ cdots}{\ displaystyle E_ {n} = - {\ frac {me ^ {4}} {32 \ pi ^ {2} \ epsilon _ {0} ^ {2} \ hbar ^ {2} (n-1/2) ^ {2}}} \ quad; \ quad n = 1,2,3, \ cdots}

Обратите внимание, что в этом случае L {\ displaystyle L}L почти совпадает с обычным квантованием оператора углового момента на плоскости L z {\ displaystyle L_ {z}}L_ {z} . Для трехмерного случая метод EBK для полного углового момента эквивалентен поправке Лангера.

См. Также
  • значок Физический портал
Ссылки
Последняя правка сделана 2021-05-18 09:42:53
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте