पीयरल्स प्रतिस्थापन

पीयरल्स प्रतिस्थापन विधि, जिसका नाम रुडोल्फ पीयरल्स के मूल कार्य के नाम पर रखा गया है धीरे-धीरे बदलती चुंबकीय सदिश क्षमता की उपस्थिति में दृढ़ बंधन (टाइट बाइंडिंग) इलेक्ट्रॉनों का वर्णन करने के लिए एक व्यापक रूप से नियोजित अनुमान है।

बाहरी चुंबकीय सदिश क्षमता $$\mathbf{A}$$ की उपस्थिति में, अनुवाद ऑपरेटर, जो तंग-बाध्यकारी दृढ़ में हैमिल्टनियन के गतिज भाग का निर्माण करते हैं, बस हैं


 * $$\mathbf{T}_x = |m+1,n\rangle\langle m,n|e^{i\theta^x_{m,n}}, \quad \mathbf{T}_y = |m,n+1\rangle\langle m,n|e^{i\theta^y_{m,n}} $$और दूसरे परिमाणीकरण सूत्रीकरण में
 * $$\mathbf{T}_x = \boldsymbol{\psi}^\dagger_{m+1,n}\boldsymbol{\psi}_{m,n}e^{i\theta^x_{m,n}}, \quad \mathbf{T}_y = \boldsymbol{\psi}^\dagger_{m,n+1}\boldsymbol{\psi}_{m,n}e^{i\theta^y_{m,n}}. $$

चरणों को इस प्रकार परिभाषित किया गया है
 * $$ \theta^x_{m,n} = \frac{q}{\hbar}\int_m^{m+1} A_x(x,n)\text{d}x, \quad \theta^y_{m,n} = \frac{q}{\hbar}\int_n^{n+1} A_y(m,y) \text{d}y. $$

गुण
\begin{align} \boldsymbol{\nabla}\times\theta_{m,n}& = \Delta_x\theta^y_{m,n}-\Delta_y\theta^x_{m,n} = \left(\theta^y_{m+1,n}-\theta^y_{m,n}-\theta^x_{m,n+1}+\theta^x_{m,n}\right)\\ & = \frac{q}{\hbar}\int_{\text{unit cell}}\mathbf{A}\cdot \text{d}\mathbf{l} = 2\pi\frac{q}{h}\int \mathbf{B} \cdot \text{d}\mathbf{s} = 2\pi\phi_{m,n} \end{align}$$ और लैटिस के माध्यम से कुल प्रवाह है $ \Phi = \Phi_0\sum_{m,n}\phi_{m,n}$ साथ $$\Phi_0 = hc/e$$ गाऊसी इकाइयों में चुंबकीय प्रवाह क्वांटम होना।
 * 1) प्रति प्लैकेट $$\phi_{mn}$$ फ्लक्स क्वांटा की संख्या चरण कारक के लैटिस कर्ल से संबंधित है,
 * 1) फ्लक्स क्वांटा प्रति प्लैकेट $$\phi_{mn}$$ एकल कण अवस्था के संचित चरण से संबंधित है, $$ |\psi\rangle = \boldsymbol{\psi}_{i,j}|0\rangle $$ एक पट्टिका के आसपास:

औचित्य
यहां हम पियरल्स प्रतिस्थापन की तीन व्युत्पत्तियां देते हैं, जिनमें से प्रत्येक क्वांटम यांत्रिकी सिद्धांत के एक अलग सूत्रीकरण पर आधारित है।

स्वयंसिद्ध दृष्टिकोण
यहां हम पीयरल्स प्रतिस्थापन की एक सरल व्युत्पत्ति दे रहे हैं, जो द फेनमैन लेक्चर्स (खंड III, अध्याय 21) पर आधारित है। यह व्युत्पत्ति बताती है कि चुंबकीय क्षेत्र को हॉपिंग शर्तों में एक चरण जोड़कर टाइट-बाइंडिंग मॉडल में शामिल किया गया है और दिखाया गया है कि यह सातत्य हैमिल्टनियन के अनुरूप है। इस प्रकार, हमारा प्रारंभिक बिंदु हॉफस्टैटर हैमिल्टनियन है:

$$ H_0 = \sum_{m,n}\bigg(-te^{i\theta^x_{m,n}}\vert m\!+\!a,n \rangle \langle m,n\vert -te^{i\theta_{m,n}^y}\vert m,n\!+\!a\rangle\langle m,n\vert -\epsilon_0\vert m,n\rangle\langle m,n\vert\bigg)+ \text{h.c}. $$ अनुवाद संचालक $$\vert m+1\rangle\langle m\vert$$ इसके जनरेटर का उपयोग करके स्पष्ट रूप से लिखा जा सकता है, जो कि गति ऑपरेटर है। इस प्रतिनिधित्व के तहत इसे दूसरे क्रम तक विस्तारित करना आसान है,

\vert m\!+\!a\rangle\langle m\vert = \exp{\bigg(\!-\!\frac{i \mathbf{p}_xa}{\hbar}\bigg)}\vert m\rangle\langle m\vert = \left(1-\frac{i\mathbf{p}_x}{\hbar}a -\frac{\mathbf{p}_x^2}{2\hbar^2}a^2+\mathcal{O}(a^3) \right)\vert m\rangle\langle m\vert $$ और एक 2डी लैटिस में $$\vert m\!+\!a\rangle\langle m\vert \longrightarrow\vert m\!+\!a,n\rangle\langle m,n\vert$$. इसके बाद, हम चरण कारकों के दूसरे क्रम तक विस्तार करते हैं, यह मानते हुए कि सदिश क्षमता एक लैटिस रिक्ति (जिसे छोटा माना जाता है) पर महत्वपूर्ण रूप से भिन्न नहीं होती है।

\begin{align} e^{i\theta}& = 1+i\theta-\frac{1}{2}\theta^2 + \mathcal{O}(\theta^3),\\ \theta &\approx \frac{aqA_x}{\hbar},\\ e^{i\theta}& = 1 + \frac{iaqA_x}{\hbar} - \frac{a^2q^2A^2_x}{2\hbar^2} + \mathcal{O}(a^3). \end{align} $$ इन विस्तारों को हैमिल्टनियन यील्ड के प्रासंगिक हिस्से में प्रतिस्थापित करना

\begin{align} e^{i\theta}\vert m+a\rangle\langle m\vert +e^{-i\theta}\vert m\rangle\langle m+a\vert& = \bigg(1 + \frac{iaqA_x}{\hbar} - \frac{a^2q^2A^2_x}{2\hbar^2} + \mathcal{O}(a^3)\bigg)\bigg(1-\frac{i\mathbf{p}_x}{\hbar}a-\frac{\mathbf{p}_x^2}{2\hbar^2}a^2 +\mathcal{O}(a^3)\bigg)\vert m\rangle\langle m\vert+\text{h.c}\\ & = \bigg( 2-\frac{\mathbf{p}^2_x}{\hbar^2}a^2+\frac{q\lbrace \mathbf{p}_x,A_x \rbrace}{\hbar^2}a^2-\frac{q^2A_x^2}{\hbar^2}a^2+\mathcal{O}(a^3)\bigg) \vert m\rangle\langle m\vert\\ & = \bigg(-\frac{a^2}{\hbar^2}\big(\mathbf{p}_x-qA_x\big)^2+2+\mathcal{O}(a^3)\bigg) \vert m\rangle\langle m\vert. \end{align} $$ 2डी मामले में अंतिम परिणाम को सामान्यीकृत करते हुए, हम सातत्य सीमा पर हॉफस्टैटर हैमिल्टनियन पर पहुंचते हैं:
 * $$H_0 = \frac{1}{2m}\big(\mathbf{p}-q\mathbf{A}\big)^2+\tilde{\epsilon_0}$$

जहाँ प्रभावी द्रव्यमान है $$m = \hbar^2/2ta^2$$ और $$\tilde{\epsilon}_0 = \epsilon_0-4t$$.

अर्ध-शास्त्रीय दृष्टिकोण
यहां हम दिखाते हैं कि पीयरल्स चरण कारक गतिशील शब्द के कारण चुंबकीय क्षेत्र में एक इलेक्ट्रॉन के प्रसारक से उत्पन्न होता है $$q\mathbf{v}\cdot\mathbf{A} $$ लैग्रेंजियन में दिखाई दे रहा है। पथ अभिन्न सूत्रीकरण में, जो शास्त्रीय यांत्रिकी के क्रिया सिद्धांत को सामान्यीकृत करता है, साइट से संक्रमण आयाम $$j$$ समय पर $$t_j$$ साइट को $$i$$ समय पर $$t_i$$ द्वारा दिया गया है

जहाँ एकीकरण ऑपरेटर, $$ \int_{\mathbf{r}(t_i)}^{\mathbf{r}(t_j)} \mathcal{D}[\mathbf{r}(t)] $$ से सभी संभावित पथों के योग को दर्शाता है $$\mathbf{r}(t_i)$$ को $$\mathbf{r}(t_j)$$ और $$\mathcal{S}[\mathbf{r}_{ij}] = \int_{t_i}^{t_j} L[\mathbf{r}(t),\dot{\mathbf{r}}(t),t] \mathrm{d}t$$ शास्त्रीय क्रिया (भौतिकी) है, जो एक कार्यात्मक है जो एक प्रक्षेपवक्र को अपने तर्क के रूप में लेती है। हम उपयोग करते हैं $$\mathbf{r}_{ij}$$ अंतबिंदुओं के साथ एक प्रक्षेपवक्र को दर्शाने के लिए $$r(t_i), r(t_j)$$. प्रणाली के लैग्रेंजियन को इस प्रकार लिखा जा सकता है जहाँ $$ L^{(0)} $$ चुंबकीय क्षेत्र की अनुपस्थिति में लैग्रेंजियन है। संबंधित क्रिया पढ़ती है अब, यह मानते हुए कि केवल एक ही मार्ग दृढ़ता में योगदान देता है, हमारे पास है इसलिए, एक चुंबकीय क्षेत्र के अधीन एक इलेक्ट्रॉन का संक्रमण आयाम एक चरण में चुंबकीय क्षेत्र की अनुपस्थिति में होता है।

एक और व्युत्पत्ति
हैमिल्टनियन द्वारा दिया गया है
 * जहाँ $$ U\left(\mathbf{r}\right)$$ क्रिस्टल लैटिस के कारण संभावित परिदृश्य है। बलोच प्रमेय का दावा है कि समस्या का समाधान:$$H\Psi_{\mathbf{k}}(\mathbf{r}) = E\left(\mathbf{k}\right)\Psi_{\mathbf{k}}(\mathbf{r})$$, बलोच योग प्रपत्र में मांगा जाना है

\Psi_{\mathbf{k}}(\mathbf{r}) = \frac{1}{\sqrt{N}}\sum_{\mathbf{R}}e^{i\mathbf{k}\cdot\mathbf{R}}\phi_\mathbf{R}\left(\mathbf{r}\right), $$ जहाँ $$N$$ इकाई सेल्स की संख्या है, और $$\phi_\mathbf{R}$$ वानियर फलन के रूप में जाने जाते हैं। संगत आइगेन मान $$E\left(\mathbf{k}\right)$$, जो क्रिस्टल गति के आधार पर बैंड बनाते हैं $$\mathbf{k}$$, आव्यूह तत्व की गणना करके प्राप्त किए जाते हैं

E\left(\mathbf{k}\right) = \int d\mathbf{r}\ \Psi_\mathbf{k}^*(\mathbf{r})H\Psi_\mathbf{k}(\mathbf{r}) = \frac{1}{N}\sum_{\mathbf{R}\mathbf{R}^{\prime}}e^{i\mathbf{k}\left(\mathbf{R}^{\prime}-\mathbf{R}\right)} \int d\mathbf{r}\ \phi^*_\mathbf{R}\left(\mathbf{r}\right)H\phi_{\mathbf{R}^{\prime}}\left(\mathbf{r}\right) $$ और अंततः सामग्री-निर्भर होपिंग इंटीग्रल्स पर निर्भर होते हैं
 * $$t_{12} = -\int

d\mathbf{r}\ \phi^*_{\mathbf{R}_1}\left(\mathbf{r}\right)H\phi_{\mathbf{R}_2}\left(\mathbf{r}\right).$$ चुंबकीय क्षेत्र की उपस्थिति में हैमिल्टनियन में परिवर्तन होता है जहाँ $$q$$ कण का आवेश है. इसमें संशोधन करने के लिए, वानियर फलन को बदलने पर विचार करें

\begin{align} \tilde{\phi}_\mathbf{R}(\mathbf{r}) = e^{i \frac{q}{\hbar} \int_\mathbf{R}^\mathbf{r} \mathbf{A}(\mathbf{r}',t) \cdot dr'} \phi_\mathbf{R}(\mathbf{r}), \end{align} $$ जहाँ $$\phi_\mathbf{R} \equiv \tilde{\phi}_\mathbf{R}(\mathbf{A}\to 0)$$. यह नई बलोच तरंग को कार्यशील बनाता है

\tilde{\Psi}_\mathbf{k}(\mathbf{r}) = \frac{1}{\sqrt{N}} \sum_{\mathbf{R}} e^{i \mathbf{k}\cdot\mathbf{R}} \tilde{\phi}_\mathbf{R}(\mathbf{r}), $$ समय पर पूर्ण हैमिल्टनियन के मूल अवस्था में $$t$$, पहले जैसी ही ऊर्जा के साथ है। इसे देखने के लिए हम सबसे पहले प्रयोग करते हैं $$\mathbf{p} = -i \hbar \nabla$$ लिखना
 * $$\begin{align}

\tilde{H}(t) {\tilde{\phi}_\mathbf{R}(\mathbf{r})} & = \left[ \frac{(\mathbf{p} - q\mathbf{A}(\mathbf{r},t))^2}{2m} + U(\mathbf{r}) \right] e^{i\frac{q}{\hbar} \int_\mathbf{R}^\mathbf{r} \mathbf{A}(\mathbf{r}',t) \cdot d\mathbf{r}'} \phi_\mathbf{R}(\mathbf{r}) \\ & = e^{i\frac{q}{\hbar} \int_\mathbf{R}^\mathbf{r} A(\mathbf{r}',t) \cdot d\mathbf{r}'} \left[\frac{(\mathbf{p} - q\mathbf{A}(\mathbf{r},t) + q \mathbf{A}(\mathbf{r},t))^2}{2m} + U(\mathbf{r}) \right] \phi_\mathbf{R}(\mathbf{r}) \\ & = e^{i\frac{q}{\hbar} \int_\mathbf{R}^\mathbf{r} A(\mathbf{r}',t) \cdot d\mathbf{r}'} H \phi_\mathbf{R}(\mathbf{r}). \end{align}$$ फिर जब हम अर्ध-संतुलन में होपिंग इंटीग्रल की गणना करते हैं (यह मानते हुए कि सदिश क्षमता धीरे-धीरे बदलती है)

\begin{align} \tilde{t}_{\mathbf{R}\mathbf{R}'}(t)& = -\int d\mathbf{r}\ \tilde{\phi}_\mathbf{R}^*(\mathbf{r})\tilde{H}(t)\tilde{\phi}_{\mathbf{R}'}(\mathbf{r}) \\ & = - \int d\mathbf{r}\ \phi_\mathbf{R}^*(\mathbf{r})e^{i\frac{q}{\hbar} \left[-\int_{\mathbf{R}}^\mathbf{r} \mathbf{A}(\mathbf{r}',t)\cdot d\mathbf{r}'+\int_{\mathbf{R}'}^\mathbf{r} \mathbf{A}(\mathbf{r}',t)\cdot d\mathbf{r}'\right] } H \phi_{\mathbf{R}'}(\mathbf{r}) \\ & = - e^{i\frac{q}{\hbar}\int_{\mathbf{R}'}^{\mathbf{R}} \mathbf{A}(\mathbf{r}',t)\cdot d\mathbf{r}' } \int d\mathbf{r}\ \phi_\mathbf{R}^*(\mathbf{r})e^{i\frac{q}{\hbar}\Phi_{\mathbf{R}',\mathbf{r},\mathbf{R}}} H \phi_{\mathbf{R}'}(\mathbf{r}), \end{align} $$ जहाँ हमने परिभाषित किया है $$\Phi_{\mathbf{R}',\mathbf{r},\mathbf{R}} = \oint_{\mathbf{R}'\to \mathbf{r} \to \mathbf{R} \to \mathbf{R}'}\mathbf{A}(\mathbf{r}',t)\cdot d\mathbf{r}'$$, तीन स्थिति तर्कों द्वारा बनाए गए त्रिभुज के माध्यम से प्रवाह है। चूंकि हम मान लेते हैं $$\mathbf{A}(\mathbf{r},t)$$ लैटिस पैमाने पर लगभग एक समान है - वह पैमाना जिस पर वानियर अवस्था को पदों पर स्थानीयकृत किया जाता है $$\mathbf{R}$$ - हम अनुमान लगा सकते हैं $$\Phi_{\mathbf{R},\mathbf{r},\mathbf{R}'} \approx 0$$, वांछित परिणाम दे रहा है,

इसलिए, उठाए गए चरण कारक के अलावा, आव्यूह तत्व चुंबकीय क्षेत्र के बिना मामले के समान हैं, जिसे पीयरल्स चरण कारक दर्शाया गया है। यह अत्यधिक सुविधाजनक है, तब से हमें चुंबकीय क्षेत्र मान की परवाह किए बिना समान सामग्री मापदंडों का उपयोग करने को मिलता है, और संबंधित चरण को ध्यान में रखना संगणनात्मक रूप से तुच्छ है। इलेक्ट्रॉनों के लिए ($$ q = -e $$) यह हॉपिंग शब्द को प्रतिस्थापित करने के समान है $$t_{ij}$$ साथ $$t_{ij}e^{-i\frac{e}{\hbar}\int_i^j\mathbf{A}\cdot d\mathbf{l}}$$