Dielektryk Hopfielda

Dielektryk Hopfielda – w mechanice kwantowej model dielektryka składającego się z kwantowych oscylatorów harmonicznych oddziałujących z modami kwantowego pola elektromagnetycznego.

Oddziaływanie kolektywne modów polaryzacji ładunku ze wzbudzeniami próżni, fotonami prowadzi do zaburzenia liniowej relacji dyspersji fotonów oraz stałej dyspersji fal ładunku poprzez uniknięcie przecięcia między dwiema liniami dyspersji polarytonów[1]. Podobnie do fononów akustycznych i optycznych daleko od rezonansu jedna gałąź dyspersji zachowuje się jak fotony a druga jak fale ładunku. Matematycznie dielektryk Hopfielda dla jednego modu wzbudzeń jest równoważny paczce Trojanskiej w przybliżeniu harmonicznym. Model Hopielda dielektryka przewiduje istnienie wiecznie związanych fotonów podobnych do promieniowania Hawkinga wewnątrz materii o gęstości proporcjonalnej do siły sprzężenia pomiędzy polem i materią.

Teoria

Hamiltonian skwantowanego dielektryka Lorentza składającego się z N {\displaystyle N} oscylatorów harmonicznych oddziałujących z kwantowym polem elektromagnetycznym może być zapisany w przybliżeniu dipolowym jako:

H = A = 1 N p A 2 2 m + m ω 2 2 x A 2 e x A E ( r A ) + λ = 1 2 d 3 k a λ k + a λ k c k , {\displaystyle H=\sum \limits _{A=1}^{N}{\frac {{p_{A}}^{2}}{2m}}+{\frac {m\omega ^{2}}{2}}{x_{A}}^{2}-e{x_{A}}\cdot E(r_{A})+\sum \limits _{\lambda =1}^{2}\int d^{3}ka_{\lambda k}^{+}a_{\lambda k}\hbar ck,}

gdzie:

E ( r A ) = i L 3 λ = 1 2 d 3 k [ c k 2 ϵ 0 ] 1 2 [ e λ ( k ) a λ ( k ) exp ( i k r A ) H . C . ] {\displaystyle E(r_{A})={\frac {i}{L^{3}}}\sum \limits _{\lambda =1}^{2}\int d^{3}k\left[{\frac {ck}{2\epsilon _{0}}}\right]^{\frac {1}{2}}[e_{\lambda }(k)a_{\lambda }(k)\exp(ikr_{A})-H.C.]}

jest operatorem pola elektrycznego działającym w punkcie położenia r A . {\displaystyle r_{A}.}

Wyrażając go przy pomocy operatorów kreacji i anihilacji oscylatorów harmonicznych otrzymujemy

H = A = 1 N ( a A + a A ) ω e 2 β ( a A + a A + ) E ( r A ) + λ k a λ k + a λ k c k . {\displaystyle H=\sum \limits _{A=1}^{N}(a_{A}^{+}\cdot a_{A})\hbar \omega -{\frac {e}{{\sqrt {2}}\beta }}(a_{A}+{a_{A}}^{+})\cdot E(r_{A})+\sum _{\lambda }\sum _{k}a_{\lambda k}^{+}a_{\lambda k}\hbar ck.}

Zakładając, że oscylatory położone są na węzłach jakiejś regularnej sieci krystalicznej ciała stałego i stosując polarytonową transformatę Fouriera

B k + = 1 N A = 1 N exp ( i k r A ) a A + , {\displaystyle B_{k}^{+}={\frac {1}{\sqrt {N}}}\sum \limits _{A=1}^{N}\exp(ikr_{A})a_{A}^{+},}
B k = 1 N A = 1 N exp ( i k r A ) a A {\displaystyle B_{k}={\frac {1}{\sqrt {N}}}\sum \limits _{A=1}^{N}\exp(-ikr_{A})a_{A}}

oraz definiując rzuty fal ładunku oscylatorów na kierunki polaryzacji pola elektromagnetycznego

B λ k + = e λ ( k ) B k + , {\displaystyle B_{\lambda k}^{+}=e_{\lambda }(k)\cdot B_{k}^{+},}
B λ k = e λ ( k ) B k , {\displaystyle B_{\lambda k}=e_{\lambda }(k)\cdot B_{k},}

po pominięciu wkładów podłużnych nie oddziałujących z polem elektromagnetycznym możemy uzyskać Hamiltonian Hopfielda

H = λ k ( B λ k + B λ k + 1 2 ) ω + c k a λ k + a λ k + i e ϵ 0 m ω N V c k [ B λ k a λ k + B λ k + a λ k B λ k + a λ k + B λ k a λ k + ] . {\displaystyle H=\sum _{\lambda }\sum _{k}(B_{\lambda k}^{+}B_{\lambda k}+{\frac {1}{2}})\hbar \omega +\hbar cka_{\lambda k}^{+}a_{\lambda k}+{\frac {ie\hbar }{\sqrt {\epsilon _{0}m\omega }}}{\sqrt {\frac {N}{V}}}{\sqrt {ck}}[B_{\lambda k}a_{\lambda -k}+B_{\lambda k}^{+}a_{\lambda k}-B_{\lambda k}^{+}a_{\lambda -k}^{+}-B_{\lambda k}a_{\lambda k}^{+}].}

Ponieważ oddziaływanie nie miesza polaryzacji ten może być przekształcony do postaci normalnej z częstościami własnymi gałęzi polarytonowych

H = λ k [ Ω + ( k ) C λ + k + C λ + k + Ω ( k ) C λ k + C λ k ] + c o n s t {\displaystyle H=\sum _{\lambda }\sum _{k}\left[\Omega _{+}(k)C_{\lambda +k}^{+}C_{\lambda +k}+\Omega _{-}(k)C_{\lambda -k}^{+}C_{\lambda -k}\right]+const}

z równaniem własnym

[ C λ ± k , H ] = Ω ± ( k ) C λ ± k , {\displaystyle [C_{\lambda \pm k},H]=\Omega _{\pm }(k)C_{\lambda \pm k},}
C λ ± k = c 1 a λ k + c 2 a λ k + c 3 a λ k + + c 4 a λ k + + c 5 B λ k + c 6 B λ k + c 7 B λ k + + c 8 B λ k + , {\displaystyle C_{\lambda \pm k}=c_{1}a_{\lambda k}+c_{2}a_{\lambda -k}+c_{3}a_{\lambda k}^{+}+c_{4}a_{\lambda -k}^{+}+c_{5}B_{\lambda k}+c_{6}B_{\lambda -k}+c_{7}B_{\lambda k}^{+}+c_{8}B_{\lambda -k}^{+},}

gdzie:

Ω ( k ) 2 = ω 2 + Ω 2 ( ω 2 Ω 2 ) 2 + 4 g ω 2 Ω 2 2 , {\displaystyle \Omega _{-}(k)^{2}={\frac {\omega ^{2}+\Omega ^{2}-{\sqrt {{(\omega ^{2}-\Omega ^{2})}^{2}+4g\omega ^{2}\Omega ^{2}}}}{2}},}
Ω + ( k ) 2 = ω 2 + Ω 2 + ( ω 2 Ω 2 ) 2 + 4 g ω 2 Ω 2 2 , {\displaystyle \Omega _{+}(k)^{2}={\frac {\omega ^{2}+\Omega ^{2}+{\sqrt {{(\omega ^{2}-\Omega ^{2})}^{2}+4g\omega ^{2}\Omega ^{2}}}}{2}},}

z

Ω ( k ) = c k , {\displaystyle \Omega (k)=ck,}

(dyspersja fotonów w próżni)

i

g = N e 2 V m ϵ 0 ω 2 {\displaystyle g={\frac {Ne^{2}}{Vm\epsilon _{0}\omega ^{2}}}}

jest bezwymiarową stałą sprzężenia proporcjonalną do gęstości N / V {\displaystyle N/V} (ilości oscylatorów na jednostkę objętości) dielektryka z częstością Lorentza ω {\displaystyle \omega } (dyspersja fal ładunku w przybliżeniu ciasnego wiązania). Można zauważyć, że w odróżnieniu od próżni pola elektromagnetycznego bez materii wartość oczekiwana średniej liczby fotonów < a λ k + a λ k > {\displaystyle <a_{\lambda k}^{+}a_{\lambda k}>} jest różna od zera w stanie podstawowym Hamlitoniamu polarytonowego C k ± | 0 >= 0 {\displaystyle C_{k\pm }|\mathbf {0} >=0} podobnie do promieniowania Hawkinga w pobliżu czarnej dziury z powodu efektu Unruha-Daviesa. Można z łatwością zauważyć, że mniejsza częstość Ω {\displaystyle \Omega _{-}} staje się urojona, kiedy stała sprzężenia przekracza wartość krytyczna g > 1 , {\displaystyle g>1,} co sugeruje, że dielektryk Hopfielda ulega nadpromienistej przemianie fazowej.

Przypisy

  1. J.JJ.J. Hopfield J.JJ.J., Theory of the Contribution of Excitons to the Complex Dielectric Constant of Crystals, wyd. 5, t. 112, Physical Review, 1958, s. 1555–1567, DOI: 10.1103/PhysRev.112.1555, Bibcode: 1958PhRv..112.1555H .