Método das Relações de Recorrência aplicado no problema do elétron em campo magnético constante e uniforme

Problem of a single electron in an external magnetic field solved via Method of Recurrence Relations

Resumo

Nesse trabalho utilizamos o Método das Relações de Recorrência (MRR) para estudar o comportamento de um elétron na presença de um campo magnético uniforme aplicado ao longo do eixo z de um sistema de referência estabelecido. Encontramos que as funções de autocorrelação temporais nas direções x e y, respectivamente C_x(t) e C_y(t), oscilam periodicamente no tempo, enquanto na direção z, a função de autocorrelação temporal, C_z(t), tem valor constante e igual a um, corroborando o resultado esperado de que o spin do elétron precessa em torno da direção de aplicação do campo magnético uniforme e independente do tempo.

Palavras-chave:
Dinâmica de Spins; Método das Relações de Recorrência

Abstract

We use the Method of Recurrence Relations (MRR) to study the behavior of a single electron in the presence of a uniform magnetic field along the z-direction. We find that the time-dependent autocorrelation functions in the x- and y-directions, C_x(t) and C_y(t), respectively, oscillate periodically in time, while in the z-direction the temporal autocorrelation function C_z(t) has a constant value, corroborating the expected result that the spin precesses about the direction of the applied magnetic field.

Keywords:
Spin Dynamics; Method of Recurrence Relations

1. Introdução

Materiais magnéticos estão presentes no nosso dia a dia desde as aplicações mais simples como, por exemplo, os imãs de geladeira até as mais sofisticadas como armazenamento de informações em computadores. Atualmente, uma quantidade significativa de informações está gravada magneticamente em discos rígidos de computadores. A gravação das informações é feita utilizando-se o código binário (0 ou 1). Magneticamente, 0 poderia ser representado como spin orientado para baixo e 1 por spin orientado para cima em relação a algum sistema de referência. Hoje em dia, muitos pesquisadores trabalham na spintrônica (ou eletrônica de spins) ^[1] que visa controlar a corrente num dispositivo eletrônico, não somente pela carga do elétron, mas também pelo seu spin. Os objetivos principais são o de se aumentar a velocidade de processamento de informações e diminuir as perdas de energia nas transmissões de informação. Além disso, em Medicina, por exemplo, temos métodos de diagnóstico baseados em ressonância de spins, tanto eletrônicos (Espectroscopia de Ressonância Paramagnética Eletrônica) quanto nucleares (Ressonância Magnética Nuclear) ^[2,3]. Vê-se, por tudo isso, a importância do estudo de sistemas de spins.

Existem muitas técnicas analíticas ou numéricas ^[4] para a obtenção de soluções exatas ou aproximadas para modelos de spins. Dentre essas técnicas temos o Método das Relações de Recorrência (MRR) ^[5,6,7], que é bastante útil na obtenção da dinâmica de sistemas de muitos corpos quânticos. Esse método permite obter informações apuradas para propriedades do sistema como, funções correlação, funções espectrais e correntes de spins, dentre outras. O MRR utiliza o processo Grarn-Schmidt de ortogonalização, o qual é bastante geral e permite a construção de um conjunto completo de vetores ortogonalizados do espaço de Hilbert. O método de relações de recorrência foi formulado por Howard Lee ^[8,9], como uma generalização do formalismo desenvolvido, dentre outros, por Mori ^[10,11].

2. Spin do elétron e elétron em campo magnético constante e uniforme

2.1. Spin do elétron

Entre 1921 e 1922, Otto Stern e Walther Gerlach realizaram um experimento que visava observar se havia quantização espacial do momento angular orbital dos elétrons em átomos ^[12]. O experimento basicamente consistia em fazer um feixe de átomos, inicialmente de prata, passar por um campo magnético não homogêneo ao longo de um eixo (escolhido como eixo z em relação ao referencial do laboratório) produzido por um imã e então medir a distribuição desses átomos em um anteparo na saída do imã. Eles observaram que, aproximadamente metade dos átomos colidia com uma extremidade da placa e a outra metade na posição simetricamente oposta, não registrando praticamente nenhum átomo em qualquer posição intermediaria. A física clássica não explicava tal resultado, uma vez que, submetido a um campo magnético, um momento magnético clássico exibe um espectro continuo de valores. Em 1927 Phipps e Taylor repetiram o experimento de Stern e Gerlach com átomos neutros e com momento angular orbital nulo. Era esperado que tais átomos passassem pela região de campo magnético não homogêneo sem sofrerem nenhuma deflexão. No entanto, eles observaram que o feixe também se dividia em dois ao passar pela região com campo magnético. Esses experimentos indicavam que os átomos possuíam uma espécie de momento magnético intrínseco (M) que não estava relacionado ao momento orbital dos elétrons. Portanto, deveria existir também um momento angular intrínseco, que foi chamado de spin e representado pelo vetor S. A relação entre o momento magnético intrínseco e o momento angular intrínseco é dada por:

(1) $$\overrightarrow{M}=g\overrightarrow{S}.$$

A constante de proporcionalidade, g, é chamada razão giromagnética. O spin do elétron é um momento angular que existe independentemente do estado de movimento no qual a partícula se encontra. Cada partícula elementar tem um spin imutável associado a ela. Por exemplo, o spin do elétron é 1/2, o spin dos fótons é 1 e o spin dos grávitons é 2. Vale lembrar que partículas com spins semi-inteiros são ditas férmions ^[12] e obedecem ao princípio de exclusão de Pauli. Já partículas com valores de spins inteiros são chamadas de bósons ^[12] e não precisam obedecer ao princípio de exclusão de Pauli.

Wolfgang Pauli postulou que um elétron podia ser modelado por um vetor bidimensional complexo e criou uma teoria de momento angular que atuasse neste espaço. O estado mais geral que pode ser associando a uma partícula de spin 1/2, pode ser expresso como uma matriz coluna de dois elementos.

(2) $$\xi =\left({}_b^a\right)=\alpha {\xi }_{+}+\beta {\xi }_{-},$$

em que o índice + representa uma matriz que descreve um spin apontando para cima (↑) em relação a algum sistema de referência e o índice – representa o spin apontando para baixo (↓). Ou seja,

(3) $${\xi }_{+}=\left({}_0^1\right)=|\uparrow \rangle \text{e}{\xi }_{-}=\left({}_0^1\right)=|\downarrow \rangle ,$$

representam os estados acessíveis para o spin de uma partícula quântica. No que diz respeito ao estado de spin de uma partícula, são necessários dois números para caracterizar completamente o estado, são os números quânticos m (associado ao momento magnético) e o número quântico s (associado ao spin) que fornece a orientação espacial do momento magnético, ou seja, nos se informa se o spin aponta para cima ou para baixo em um dado sistema de referência. Para o elétron, s=1/2 e m₊ = +1/2, (↑), enquanto m₋ = -1/2, (↓).

Na descrição de Pauli, as componentes do momento angular de spin são dadas por

(4) $$S_x=\frac{\hslash }{2}{\sigma }_x,S_y=\frac{\hslash }{2}{\sigma }_{y}e{S}_z=\frac{\hslash }{2}{\sigma }_z,$$

onde ℏ é a constante de Planck, h, dividida por 2π.

As matrizes σ_x, σ_y e σ_z são ditas matrizes de Pauli com os seguintes valores:

(5) $$\begin{array}{l}{\sigma }_x=\left({}_1^0{}_0^1\right),{\sigma }_y=\left({}_i^0{}_0^{-i}\right)\hfill \\ \text{e}{\sigma }_z=\left({}_0^1{}_{-1}^0\right).\hfill \end{array}$$

Além disso, temos as seguintes relações de comutação para as componentes do momento angular:

(6) $$\begin{array}{l}\left[S_x,S_y\right]=i\hslash S_z\left[S_y,S_z\right]=i\hslash S_x\left[S_z,S_x\right]=i\hslash S_y\\ \left[S_y,S_x\right]=-i\hslash S_z\left[S_z,S_y\right]=i\hslash S_x\left[S_x,S_z\right]=-i\hslash S_y\end{array}$$

2.2. Eletron em campo magnético externo constante e uniforme

Vimos anteriormente que os elétrons possuem momento magnético intrínseco (M) relacionado ao seu momento angular intrínseco (S).

(7) $$\overrightarrow{M}=g\overrightarrow{S}=g\left(S_x\widehat{i}+S_y\widehat{j}+S_z\widehat{k}\right).$$

Quando um dipolo magnético é posicionado em um campo magnético B, o seu momento magnético tende a se alinhar com o campo e a menor energia de configuração do sistema ocorre quando o momento angular da partícula alinha-se com o campo externo.

No sistema estudado temos um elétron em repouso em um campo magnético uniforme que aponta para a direção z, ou seja,

(8) $$\overrightarrow{B}=B_0\widehat{k}.$$

A energia associada à interação do elétron com o campo é, então, dada pelo hamiltoniano:

(9) $$H=-\overrightarrow{M}\cdot \overrightarrow{B}=-g\left(S_x\widehat{i}+S_y\widehat{j}+S_z\widehat{k}\right)\cdot \text{}{B}_0\widehat{k}=-g{B}_0{S}_z.$$

A forma matricial do hamiltoniano é obtida via matrizes de Pauli,

(10) $$H=-g{B}_0{S}_z=-gB0\frac{\hslash }{2}{\sigma }_z=-g{B}_0\frac{\hslash }{2}\left({}_0^1{}_{-1}^0\right).$$

Sabemos que os autovetores (estados) e autovalores (no caso, energias) podem ser obtidos por

(11) $$H{\xi }_{+}=E_{+}{\xi }_{+}eH{\xi }_{-}=E_{-}{\xi }_{-}.$$

Os autovetores e autovalores são, então

(12) $${\xi }_{+}=\left({}_0^1\right)\text{com}{E}_{+}=-\left(\frac{g{B}_0\hslash }{2}\right),$$

(13) $${\xi }_{-}=\left({}_1^0\right)\text{com}{E}_{-}=+\left(\frac{g{B}_0\hslash }{2}\right).$$

A solução da equação de onda que descreve o sistema pode ser expressa em termos de estados estacionários:

(14) $$\begin{array}{l}\xi \left(t\right)=\alpha {\xi }_{+}{e}^{\frac{iE+t}{\hslash }}+\beta {\xi }_{-}e{-}^{\frac{iE-t}{\hslash }}=\alpha \left({}_0^1\right)e^{\frac{ig{B}_{0}t}{2}}+\\ \beta \left({}_1^0\right)e^{-\frac{ig{B}_{0}t}{2}}=\left({}_{\beta e^{-\frac{ig{B}_{0}t}{2}}}^{\alpha e^{\frac{ig{B}_{0}t}{2}}}\right).\end{array}$$

Para calcular as constantes α e β lembramos que a solução deve ser normalizada para qualquer tempo t e, em t=0, temos

(15) $$\xi \left(0\right)=\left({}_{\beta }^{\alpha }\right).$$

A condição de normalização nos diz que

(16) $$\langle \xi \left(0\right)|\xi \left(0\right)\rangle =1\Rightarrow \left(\alpha \beta \right)\left({}_{\beta }^{\alpha }\right)=1.$$

Logo |α|² + |β|² = 1 e

(17) $$\alpha =\cos \left(\frac{\theta }{2}\right)e\beta =sen\left(\frac{\theta }{2}\right),$$

em que θ é um ângulo fixo.

Logo, temos que a solução da equação de onda dependente do tempo, para uma partícula de spin 1/2 em um campo magnético uniforme:

(18) $$\xi \left(t\right)=\left({}_{sen\frac{\theta }{2}e{}^{-i\frac{g{B}_{0}t}{2}}}^{\cos \frac{\theta }{2}e{}^{i\frac{g{B}_{0}t}{2}}}\right).$$

Calculando o valor esperado do vetor momento angular de spin S, como uma função de tempo nós podemos obter informações úteis sobre a dinâmica desse sistema. O valor esperado de <S >é definido como

(19) $$\langle \overrightarrow{S}\rangle =\langle S_x\rangle \widehat{i}+\langle S_y\rangle \widehat{j}+\langle S_z\rangle \widehat{k},$$

com

(20) $$\begin{array}{l}\langle S_x\rangle =\langle \xi \left(t\right)|S_x|\xi \left(t\right)\rangle ,\langle S_y\rangle =\langle \xi \left(t\right)|S_y|\xi \left(t\right)\rangle \\ \text{e}\langle S_z\rangle =\langle \xi \left(t\right)|S_z|\xi \left(t\right)\rangle \end{array}$$

Para a componente x, teremos

(21) $$\begin{array}{l}\langle S_x\rangle =\left(\cos \left(\frac{\theta }{2}\right)e^{-i\frac{g{B}_{0}t}{2}}sen\left(\frac{\theta }{2}\right)e^{i\frac{g{B}_{0}t}{2}}\right).\\ \left[\frac{\hslash }{2}\left({}_1^0{}_0^1\right)\left({}_{en\left(\frac{\theta }{2}\right)e^{-i\frac{g{B}_{0}t}{2}}}^{\cos \left(\frac{\theta }{2}\right)e{}^{i\frac{g{B}_{0}t}{2}}}\right)\right].\end{array}$$

O que leva a

(22) $$\langle S_x\rangle =\frac{\hslash }{2}sen\left(\theta \right)\cos \left(g{B}_{0}t\right).$$

Analogamente,

(23) $$\langle S_y\rangle =-\frac{\hslash }{2}sen\left(\theta \right)sen\left(g{B}_{0}t\right)\text{e}\langle S_z\rangle =\frac{\hslash }{2}\cos \left(\theta \right).$$

Observamos então que as componentes x e y do momento angular de spin variam periodicamente no tempo, enquanto a componente z tem um valor constante. Representando o valor esperado de <S >em um sistema de coordenadas cartesiano, notamos que <S >está inclinado de um ângulo constante θ em relaηão ao eixo z e precessiona em torno desse eixo com uma frequência ω = g B₀. Essa frequência é a conhecida frequência de Larmor. Na Figura 1 apresentamos um diagrama do comportamento do vetor momento magnético de spin, S, do elétron numa região do espaço com campo magnético uniforme e independente do tempo.

Figura 1:
Precessão do momento angular total de spin S em torno da direção de aplicação do campo magnético externo uniforme. O ângulo formado entre o vetor S e o eixo z é igual a θ.

3. Método das relações de recorrência

Para estudarmos a dinâmica de um operador hermitiano (possui autovalores reais) em um determinado sistema físico, nós usaremos o Método das Relações de Recorrência (MRR). Esse método utiliza o processo Grarn-Schmidt de ortogonalização para a construção de uma base vetorial no espaço de Hilbert d dimensional.

Neste trabalho estamos interessados no cálculo das funções de autocorrelação temporais das variáveis dinâmicas σ_x, σ_y e σ_z, ou seja, queremos saber como as funções de autocorrelação das matrizes de Pauli comportam-se quando a partícula é submetida a um campo externo. Para obtermos essas informações utilizaremos o MRR, o qual é descrito abaixo.

Seja A um operador correspondente a um observável e H o hamiltoniano que descreve o sistema. Na representação de Heisenberg, a dependência temporal está no operador e, não na função de onda. Portanto, o operador A na descrição de Heisenberg toma a seguinte forma

(24) $$A\left(t\right)=e^{iHt}A\left(0\right)e^{-iHt},$$

na qual se fez ℏ = 1 por questões de simplicidade.

A evolução temporal de A é governada pela equação de Liouville

(25) $$\frac{d}{dt}A\left(t\right)=iLA\left(t\right).$$

Na Eq. 25, L é o operador de Liouvile, que é definido como

(26) $$LA\left(t\right)=\left[H,A\left(t\right)\right]=HA\left(t\right)-A\left(t\right)H.$$

No método das relações de recorrência, A(t) é considerado um vetor definido em um espaço de Hilbert, S, de dimensão d. A dimensionalidade do sistema dependerá do sistema tratado, bem como da variável dinâmica que se pretende estudar. O produto escalar (ou produto interno) em S é definido por meio do produto de Kubo,

(27) $$\left(X,Y\right)=\frac{1}{\beta }{\displaystyle \underset{0}{\overset{\beta }{\int }}d\lambda }\langle X\left(\lambda \right)Y^{†}\rangle -\langle X\rangle \langle Y^{†}\rangle ,$$

onde X e Y são vetores, definidos no espaço de Hilbert e os símbolos <> representam médias térmicas no ensemble canônico. Além disso, β = 1/(k_BT), em que T é a temperatura e k_B é a constante de Boltzmann. As médias térmicas no ensemble canônico são tomadas da seguinte maneira:

(28) $$\langle X\rangle =\frac{Tr{e}^{-\beta H}X}{Tr{e}^{-\beta H}}.$$

Na representação de Heisenberg,

(29) $$X\left(\lambda \right)=e^{\lambda H}X{e}^{-\lambda H}.$$

No limite de altas temperaturas (T → ∞), o produto escalar (X,Y) reduz-se a

(30) $$\left(X,Y\right)=\frac{TrX{Y}^{†}}{Tr1},$$

em que Tr1 é o número de estados do sistema.

A quantidade de interesse A(t) é dada como uma expansão em termos de vetores de uma base ortogonalizada, não normalizada, {f_ν} no espaço de Hilbert S,

(31) $$A\left(t\right)={\displaystyle \sum _{v=0}^{d-1}{a}_v\left(t\right)f_v.}$$

Sabemos que (Eq. 25),

(32) $$\frac{d}{dt}A\left(t\right)|{}_{t=0}^{}=iLA\left(t\right){|}_{t=0}=iLA\left(0\right).$$

Para verificar a relação de recorrência que gera a base de vetores f_ν, começamos definindo o primeiro vetor de base (ν = 0) como

(33) $$|f_1\rangle =iL|f_0\rangle =\Vert {\frac{d}{dt}{f}_0\left(t\right)\rangle |}_{t=0}.$$

Os vetores devem ser ortogonais, então os produtos escalares entre vetores diferentes deverão ser nulos. Portanto, <f₀|f₂>=< f₁|f₂>=0. O segundo vetor pode ser escrito como:

(34) $$|f_2\rangle =iL|f_1\rangle +|x\rangle .$$

Aplicando as condições de ortogonalidade teremos,

(35) $$\langle f_0|f_2\rangle =\langle f_0|iL|f_1\rangle +\langle f_0|x\rangle =0.$$

Como

(36) $$\langle f_0|iL=\langle f_1|,$$

obtemos,

(37) $$-\langle f_1|f_1\rangle +\langle f_0|x\rangle =0$$

De maneira análoga,

(38) $$\langle f_1|f_2\rangle =\langle f_0|iL|f_1\rangle +\langle f_1|x\rangle =0\Rightarrow \langle f_1|x\rangle =0.$$

Vamos supor que

(39) $$x=c{f}_0.$$

Levando a Eq. (39) na Eq.(37), vemos que,

(40) $$-\langle f_1|f_1\rangle +\langle f_0|c{f}_0\rangle =0$$

e ainda,

(41) $$c=\frac{\langle f_1|f_1\rangle }{\langle f_0|f_0\rangle }\equiv {\Delta }_1.$$

Portanto o segundo vetor da base (Eq. 34) fica dado por

(42) $$|f_2\rangle =iL|f_1\rangle +{\Delta }_1|f_0\rangle .$$

Os outros vetores da base são construídos maneira análoga, observando-se a ortogonalidade entre os mesmos. Os vetores f_ν podem ser construídos sistematicamente e satisfazem a seguinte relação de recorrência, dita Relação de Recorrência I (RRI),

(43) $$f_{v+1}=iL{f}_v+{\Delta }_v{f}_{v-1}0\le v\le d-1,$$

em que

(44) $${\Delta }_v=\frac{\left(f_v,f_v\right)}{\left(f_{v-1},f_{v-1}\right)},$$

com as seguintes definições,

(45) $$f_{-1}\equiv 0\text{e}{\Delta }_0\equiv 1.$$

Podemos fazer uma escolha arbitrária para o primeiro vetor de base f₀. Em geral, escolhe-se f₀ = A(0) e assim as condições iniciais impõem que a₀(0) = 1 e a_ν(0)=0 para valores de ν maiores ou iguais a 1.

Isso pode ser verificado facilmente, pois

(46) $$A\left(t\right)=a_0\left(t\right)f_0+a_1\left(t\right)f_1+\mathrm{...}+a_{d-1}\left(t\right)f_{d-1}.$$

E quando t = 0, temos,

(47) $$A\left(0\right)=a_0\left(0\right)f_0+a_1\left(0\right)f_1+\mathrm{...}+a_{d-1}\left(0\right)f_{d-1}=f_0.$$

A função a₀(t) representa a função de relaxação da teoria de resposta linear. No limite de altas temperaturas, a₀(t) é simplesmente a função de autocorrelação dependente do tempo a₀(t)=C(t). A função de autocorrelação temporal para o operador A(t) nos informa sobre a correlação do operador consigo mesmo em relação ao tempo t=0:

(48) $$C\left(t\right)=\langle \frac{A\left(t\right)A\left(0\right)}{A\left(0\right)A\left(0\right)}\rangle =\frac{\left(A\left(t\right),A\left(0\right)\right)}{\left(A\left(0\right),A\left(0\right)\right)}.$$

Expandindo A(t) na base de vetores apropriada, vemos que

(49) $$\begin{array}{l}C\left(t\right)=\frac{\left(A\left(t\right),A\left(0\right)\right)}{\left(A\left(0\right),A\left(0\right)\right)}\\ =\frac{\left(a_0\left(t\right)f_0+a_1\left(t\right)f_1+\mathrm{...},f_0\right)}{\left(f_0,f_0\right)}=a_0\left(t\right).\end{array}$$

Como a expansão de A(t) em termos dos vetores de base f_ν satisfaz a equação de Liouville, o uso da relação de recorrência dá origem a uma segunda relação de recorrência para os a_ν(t), dita Relação de Recorrência II (RRII):

(50) $${\Delta }_{v+1}{a}_{v+1}\left(t\right)=-\frac{d}{dt}{a}_v\left(t\right)+a_{v-1}\left(t\right)0\le v\le d-1,$$

em que a₋₁(t) ≡ 0.

Note que os Δ_ν são os únicos elementos que entram na construção dos a_ν(t), os quais determinam completamente a evolução temporal do operador.

Uma consequência de RRII é que a função de autocorrelação temporal a_o(t) tem derivada nula na origem, isto é,

(51) $$\frac{d}{dt}{a}_0\left(t\right)|{}_{t=0}=0,$$

o que impede, por exemplo, decaimento exponencial simples como solução exata para a_o(t).

Da relação de recorrência RRII, vemos que

(52) $$v=0{\Delta }_1{a}_1=-\frac{d}{dt}{a}_0\left(t\right).$$

Vamos aplicar agora a transformada de Laplace à Eq. 52

(53) $${\displaystyle \underset{0}{\overset{\infty }{\int }}{e}^{-zt}}{\Delta }_1{a}_1\left(t\right)dt=-{\displaystyle \underset{0}{\overset{\infty }{\int }}{e}^{-zt}}\frac{d{a}_0\left(t\right)}{dt}dt.$$

Integrando por partes o lado direito da equação, obtemos

(54) $${\Delta }_1{a}_1\left(z\right)=1-z{a}_0\left(z\right).$$

Repetindo o processo, ou seja, aplicando a transformada de Laplace aos termos com ν = 1, 2, 3, obtemos

(55) $${\Delta }_{v+1}{a}_{v+1}\left(z\right)=a_{v-1}\left(z\right)-z{a}_v\left(z\right)v\ge 1.$$

Partindo das relações de recorrência RRI e RRII, podemos escrever a₀(z) na forma de uma fração continuada:

(56) $$a_0\left(z\right)=\frac{1}{z+\frac{{\Delta }_1}{z+\frac{{\Delta }_2}{z+\mathrm{...}}}}.$$

A função de autocorrelação temporal, a₀(t), é determinada aplicando-se a transformada de Laplace inversa em a₀(z).

(57) $$a_0\left(t\right)=L^{-1}\left(a_0\left(z\right)\right).$$

Quando utilizamos o método das relações de recorrência, calculamos primeiro os vetores da base e, a seguir, os deltas. Vimos acima que os deltas determinam a₀(z) e, por conseguinte, a função de autocorrelação temporal.

4. Aplicação do Método das Relações de Recorrência no problema do elétron em um campo magnético externo uniforme

Pretendemos obter o comportamento dinâmico das componentes do momento angular do spin de um elétron imerso em um campo magnético uniforme que aponta na direção z de um sistema de coordenadas.

4.1. Função de autocorrelação temporal na direção y

A princípio obteremos os vetores da base e os deltas para posteriormente calcularmos as funções de autocorrelação temporais. Para o cálculo da função de autocorrelação temporal na direção y, escolhemos

(58) $$f_0={\sigma }_y.$$

A relação de recorrência RRI nos diz como proceder para encontrar o segundo vetor da base f₁:

(59) $$f_1=iL{f}_0+{\Delta }_0{f}_{-1}=iL{f}_0=i\left[H,f_0\right]=i\left[H,{\sigma }_y\right],$$

já que f₋₁ = 0 e Δ₀ = 0.

Como vimos na seção 2, o hamiltoniano para este problema é dado por:

(60) $$H=-\frac{g{B}_0}{2}{\sigma }_z.$$

Na qual se fez ℏ = 1 por uma questão de coerência com o MRR.

Logo,

(61) $$f_1=i\left[-\frac{g{B}_0}{2}{\sigma }_{z,}{\sigma }_y\right]=-i\frac{g{B}_0}{2}\left[{\sigma }_{z,}{\sigma }_y\right]=-i\frac{g{B}_0}{2}\left(-2i{\sigma }_x\right)=-g{B}_0{\sigma }_x.$$

Para o cálculo de Δ₁, utilizamos

(62) $${\Delta }_1=\frac{\left(f_1,f_1\right)}{\left(f_0,f_0\right)}=\frac{\left(-g{B}_0{\sigma }_x,-g{B}_0{\sigma }_x\right)}{\left({\sigma }_y,{\sigma }_y\right)}={\left(g{B}_0\right)}^2.$$

Para calcular o próximo vetor da base, devemos fazer ν = 1 na RRI:

(63) $$f_2=iL{f}_1+{\Delta }_1{f}_0=i\left[-\frac{g{B}_0}{2}{\sigma }_z,-g{B}_0{\sigma }_y\right]+{\left(g{B}_0\right)}^2{\sigma }_y=-{\left(g{B}_0\right)}^2{\sigma }_y+{\left(g{B}_0\right)}^2{\sigma }_y=0.$$

Como consequência,

(64) $${\Delta }_2=\frac{\left(f_2{f}_2\right)}{\left(f_{1,}{f}_1\right)}=0.$$

Pode-se mostrar utilizando-se RRI que f₃ = f₄ =...= f_d-1 = 0. Além disso, Δ₃ = Δ_{4 =...=} Δ_d-1 = 0.

Logo podemos calcular o valor de a₀(z),

(65) $$a_0\left(z\right)=\frac{1}{z+\frac{{\Delta }_1}{z+0}}=\frac{1}{z+\frac{{\left(g{B}_0\right)}^2}{z}}=\frac{z}{z^2+{\left(g{B}_0\right)}^2}.$$

A função de autocorrelação temporal é então obtida tomando-se a transforma de Laplace inversa de a₀(z),

(66) $$a_0\left(t\right)=L^{-1}{a}_0(z)=\cos \left(\frac{g{B}_0}{2}t\right)=C_y\left(t\right).$$

Portanto, a₀(t) pode assumir valores no intervalo de -1 a 1 na direção do eixo y.

4.2. Função de autocorrelação temporal na direção x

O problema na direção x é simétrico ao problema em y. Para se fazer os cálculos basta seguir a mesma sequência apresentada na seção 4.1. Pode-se mostrar que, para f₀ = σ_x,

(67) $$a_0\left(t\right)=L^{-1}{a}_0(z)=\cos \left(\frac{g{B}_0}{2}t\right)=C_x\left(t\right).$$

4.3. Função de autocorrelação temporal na direção z

Para direção z, escolhemos o primeiro vetor da base como

(68) $$f_0={\sigma }_z,$$

o que leva a

(69) $$f_1=i\left[-\frac{g{B}_0}{2}{\sigma }_z,{\sigma }_z\right]=-i\frac{g{B}_0}{2}\left[{\sigma }_z,{\sigma }_z\right]=i\frac{g{B}_0}{2}\left(0\right)=0$$

e Δ₁ = 0. Além disso, f₂ =f₃ =...= f_d−1 =0 e Δ₂ = Δ₃=...=Δ_d−1 =0.

Portanto,

(70) $$a_0\left(z\right)=\frac{1}{z},$$

com tranfomada de Laplace inversa,

(71) $$a_0\left(t\right)=L^{-1}\left(\frac{1}{z}\right)=1=C_z\left(t\right),$$

ou seja, a função de autocorrelação temporal na direção de aplicação do campo apresenta valor constante. O momento angular de spin S_z está sempre positivamente correlacionado consigo mesmo, para todo tempo t.

5. Considerações finais

Com a aplicação do método de relação de recorrência (MRR), conseguimos obter o comportamento dinâmico em três dimensões do spin do elétron quando submetido ao campo magnético uniforme. Verificamos que na direção do campo (eixo z) a componente do spin, S_z, possui uma função de autocorrelação temporal com valor constante e igual a 1, ou seja, C_z(t)=1. Porém, nos eixos x e y, verificou-se que essas funções oscilam periodicamente no tempo e podem assumir valores que vão desde -1 a 1. Este resultado era esperado, pois como vimos no cálculo dos valores esperados dos spins, S_x e S_y são funções periódicas do tempo, enquanto S_z apresenta um valor constante, ou seja, a componente z do spin está sempre positiva e integralmente correlacionada consigo mesma.

Referências

• ^[1]
  C. Chappert, A. Fert, and F.N. Van Dau, Nature Materials 6, 813 (2007).
• ^[2]
  A.P. Guimarães, Revista Brasileira de Ensino de Física 22, 360 (2000).
• ^[3]
  C.M.T. Dieguez, L.V. Montanheiro, L.B. Cleto, L.B.M.J.C. Bonfim and C.A. Dartora, Revista Brasileira de Ensino de Física 40, e1310 (2018).
• ^[4]
  J.B. Parkinson and D.J.J. Farnell, An Intoduction to Quantum Spin Systems Lec. Notes Phys. 816 (Springer, Berlin, 2010).
• ^[5]
  G. Muller, and V.S. Viswanath, The Recursion Method: Application to Many-Body Dynamics (Springer-Verlag, Berlin, 1994).
• ^[6]
  J. Florencio and M.H. Lee, Physical Review B 35, 1835 (1987).
• ^[7]
  M.E.S. Nunes, J.A. Plascak and J. Florencio, Physica A 332, 1 (2004).
• ^[8]
  M.H. Lee, Phys. Rev. Lett. 49, 1072 (1982).
• ^[9]
  M.H. Lee, Phys. Rev. B 26, 2547 (1982).
• ^[10]
  H. Mori, Progress of Theoretical Physics 34, 399 (1965).
• ^[11]
  H. Mori, Progress of Theoretical Physics 33, 423 (1965).
• ^[12]
  D.J. Griffiths, Mecânica Quântica (Pearson Prentice Hall, São Paulo, 2011).

Datas de Publicação

• Publicação nesta coleção
  2018

Histórico

• Recebido
  16 Nov 2017
• Revisado
  21 Dez 2017
• Aceito
  09 Jan 2018