Accessibility / Report Error

COEFICIENTES DE TRANSMISSÃO E REFLEXÃO PELO MÉTODO DA AMPLITUDE VARIÁVEL

TRANSMISSION AND REFLECTION COEFFICIENTS BY THE VARIABLE AMPLITUDE METHOD

Resumo

In this work, a simple derivation of the variable amplitude method using the variation of parameters to solve a differential equation is presented. The variable amplitude method was originally devised by Tikochinsky in 1977, using the quantum theory of scattering. The method is applied to two model potentials, the rectangular potential barrier and the Eckart potential, both with analytical solutions for the reflection coefficient. Numerical results will be compared with the exact values for several energies. The problem of calculating the reflection coefficient, usually involving extensive algebra as described in several textbooks, is reduced to solving a first order differential equation with initial condition. The method is very simple to apply, representing an attractive tool for teaching introductory quantum mechanics. A simple computer code is available from which reflection coefficients for the Eckart potential can be calculated.

Keywords:
transmission coefficient; quantum tunneling; variation of parameters


Keywords:
transmission coefficient; quantum tunneling; variation of parameters

INTRODUÇÃO

No estudo do movimento unidimensional de uma partícula de massa m e energia E, existem duas perspectivas: tratar o problema do ponto de vista da mecânica clássica ou de acordo com a mecânica quântica. No tratamento clássico, quando não ocorre perdas de energia de qualquer forma, a partícula fica presa a uma região do espaço x, tal que sua energia potencial V(x) é sempre menor ou igual a E. Em outras palavras, a energia total é conservada, V(x)+mv2/2=E. Este é um dos aspectos mais fascinantes que distingue a mecânica quântica da mecânica clássica. Na teoria quântica o estado da partícula é representado por uma função ψ(x) obtida pela solução da equação de Schrödinger, cujo módulo ao quadrado (x)|2dx é interpretado como a probabilidade de encontrar a partícula na região x+dx. Neste caso, como a função é espalhada além do limite x0, em que x0 satisfaz a relação V(x0)=E, existe uma probabilidade finita de que uma partícula com energia total E penetre em uma região do espaço onde a energia total é menor que a energia potencial V(x).

O processo, conhecido como tunelamento quântico, é de extrema relevância para explicação de várias fenômenos físicos e químicos, como o decaimento alfa,11 Gamow, G.; Z. Phys. 1928, 51, 204.,22 Gurney, R. W.; Condon, E. U.; Phys. Rev. 1929, 33, 127. a emissão de elétrons por metais frios sob a influência de um campo forte,33 Fowler, R. W.; Nordheim, L.; Proc. R. Soc. London, Ser. A 1928, 119, 173. reações de transferência de hidrogênio, fundamentais para a compreensão do efeito isotópico cinético H/D,44 Bell, R. P.; The Tunnel Effect in Chemistry, Chapman and Hall: London, 1980.,55 Johnston, H. S.; Rapp, D.; J. Am. Chem. Soc. 1961, 83, 1. mecanismos de transferência de elétrons66 Levich, V. G.; Dogonadze, R. R.; Dokl. Akad. Nauk SSSR 1960, 133, 158.,77 Marcus, R. A.; J. Chem. Phys. 1966, 45, 4493. e inversão espectral.88 Khairutdinov, R. F.; Zamaraev, K. I.; Zhdanov, V. P. Em Comprehensive Chemical Kinetics; Compton, R. G., eds.; Elsevier: Amsterdam, 1989, vol. 30.,99 Tserkis, S. T.; Moustakidis, Ch. C.; Massen, S. E.; Panos, C. P.; Phys. Lett. A 2014, 378, 497. Além disso, o tunelamento quântico é de importância para a fabricação de novos semicondutores e dispositivos eletrônicos.1010 Johnson, E. A. Em Low-Dimensional semiconductor structures; Barnham, K., Vvedensky, D., eds.; Cambridge Univ. Press: Cambridge, 2001, cap. 2.

A compreensão das propriedades físicas dos sistemas requer em muitos casos a avaliação da equação unidimensional de Schrödinger para um potencial específico. No entanto, soluções analíticas são limitadas a potenciais modelos, que embora sejam úteis para a análise da situação física em questão, nem sempre representam sistemas reais. Podemos destacar o potencial degrau e o potencial barreira retangular apresentado em diversos livros-textos para introdução do assunto.1111 Powell, J. L.; Crasemann, B.; Quantum Mechanics, Addison-Wesley: Massachusetts, 1961.,1212 Cohen-Tannoudji, C.; Diu, B.; Laloë, F.; Quantum Mechanics, Wiley: New York, 1977, vol. 1. Exemplos de potenciais para os quais se tem soluções analíticas podem ser encontrados na referência Razavy.1313 Razavy, M.; Quantum Theory of Tunneling, World Scientific: New Jersey, 2003.

Dentre os métodos numéricos destaca-se o método da matriz de transferência que é amplamente utilizado para cálculos de coeficientes de transmissão/reflexão.1414 Vatannia S.; Gildenblat, G.; IEEE J. Quantum Electron. 1996, 32, 1093.

15 Sanchez-Soto, L. L.; Monzona, J. J.; Barriuso, A. G.; Carinena, J. F.; Phys. Rep. 2012, 513, 191.

16 Mukherjee, S.; Karmakar, R.; Deyasi, A.; International Journal of Soft Computing and Engineering 2011, 1, 41.
-1717 Das, S.; Am. J. Phys. 2015, 83, 590. O trabalho de Rozman et al.1818 Rozman, M. G.; Reineker, P.; Tehver, R.; Phys. Rev. A: At., Mol., Opt. Phys. 1994, 49, 3310. destaca algumas deficiências do método, como o aumento do esforço computacional decorrente do aumento de segmentos em que o potencial é dividido e erros numéricos de arredondamento ao lidar com barreiras opacas.

Outro método, originalmente apresentado por Tikochinsky em 1977,1919 Tikochinsky, Y.; Ann. Phys. 1977, 103, 185. consiste em estabelecer uma equação diferencial de primeira ordem para a amplitude de reflexão similar a equação de Calogero.2020 Viterbo, V. D.; Lemes, N. H. T.; Braga, J. P.; Rev. Bras. Ensino Fis. 2014, 36, 1310.,2121 Calogero, F.; Variable phase approach to potential scattering, Academic Press: New York, 1967. A equação de Calogero relaciona o deslocamento de fase ao potencial intermolecular e tem sido utilizada na identificação de estados ligados2222 Braga, J. P.; Murrell, J. N.; Mol. Phys. 1984, 53, 295. e no estudo do problema inverso da obtenção do potencial intermolecular a partir de dados de deslocamento de fase.2323 Lemes, N. H. T.; Braga, J. P.; Alves, M. O.; Costa, É. D'M.; J. Mol. Model. 2014, 20, 2317. O método de Tikochinsky, método da amplitude variável, se aplica de forma geral a potenciais localizados e consiste de uma ferramenta simples e de fácil implementação computacional.

Nesse trabalho busca-se enfatizar a simplicidade do método proposto por Tikochinsky para cálculos de coeficientes de transmissão/reflexão. Embora, na literatura, tenha sido aplicado na maioria dos casos para potenciais com limites assintóticos V(x±∞) = 0,2424 Kidun, O.; Fominykh, N.; Berakdar, J.; Phys. Rev. A: At., Mol., Opt. Phys. 2005, 71, 022703. o método também é aplicável a potenciais com condições de contorno diferentes, como será discutido. Uma derivação da equação diferencial será feita por um novo caminho utilizando o método da variação dos parâmetros, uma poderosa ferramenta matemática para resolver equações diferenciais lineares não homogêneas de ordem superior. Soluções numéricas serão obtidas para coeficientes de transmissão/reflexão e comparadas com resultados analíticos para potenciais modelos.

EQUAÇÃO DIFERENCIAL PARA A AMPLITUDE DE REFLEXÃO

Uma derivação da equação diferencial para a amplitude de reflexão, com o enfoque principalmente voltado para a física de espalhamento,2525 Pliego Jr., J. R.; Braga, J. P.; Bosanac, S. D.; Phys. Rev. A: At., Mol., Opt. Phys. 1995, 53, 342. pode ser encontrada no livro de Calogero e Degasperis.2626 Calogero, F.; Degasperis, A. Em Studies in Mathematics and its Applications; Lions, J. L., Papanicolaou, G., Rockafeliar, R. T., Fujita, H., eds.; North-Holland: Amsterdam, 1982, vol. 13. Nesta seção, diferentemente da abordagem adotada no livro citado, a equação diferencial é obtida utilizando o método da variação dos parâmetros conforme método geral apresentado no trabalho de Ronveaux.2727 Ronveaux, A.; Am. J. Phys. 1969, 37, 135.

Considerando a equação unidimensional de Schrodinger

(1)

busca-se a solução em termos das soluções particulares para o problema homogêneo associado, isto é, ψ1 = eikx e ψ2 = e-ikx. A quantidade U(x) se relaciona com o potencial V(x) por U(x) = 2mV(x)/ħ2 e o número de onda k se relaciona com a energia E por m é a massa da partícula e ħ a constante reduzida de Planck., em que

A solução para o problema original, Equação 1, é colocada na forma

(2)

com a suposição de que

(3)

que é a idéia básica do método da variação dos parâmetros.2828 Zill, D. G.; Cullen, M. S.; Equações Diferenciais, Makron Books: São Paulo, 2007. A substituição de ψ e sua segunda derivada na Equação 1 fornece

(4)

As Equações 3 e 4 compõe o sistema de equações

(5)

ou na forma matricial WX ' = UWX, na qual

O vetor X ' é encontrado por X ' = MX com M = W -1UW, isto é,

(6)

As funções A(x) e B(x) modulam as amplitudes das funções de onda incidente e refletida, portanto, segundo a interpretação de Born, o quadrado da razão destas amplitudes, em x → -∞, definem a probabilidade de reflexão ou coeficiente de reflexão desejado. Assim, assumindo r(x) como r(x) = B(x) / A(x), obtém-se a equação diferencial

ou

(7)

que é uma equação de Riccati,2727 Ronveaux, A.; Am. J. Phys. 1969, 37, 135. cuja solução em x → -∞ fornece o coefiente de reflexão R = |r(x → -∞)|2. A probabilidade de transmissão ou coeficiente de transmissão pode ser facilmente obtida da relação R+T = 1, que reflete a conservação de partículas no processo.

Neste último caso, para avaliar os elementos da matriz M, calcula-se a inversa de W analiticamente. A Equação 7 aplica-se na condição em que o potencial para x → -∞ tem o valor zero. Essa condição pode ser sempre alcançada somando-se uma constante à energia e ao potencial.

A condição inicial para resolução numérica do problema, desde que as partículas incidam sobre a barreira da esquerda para a direita, é dada por r0 = r(k,x0) = x0 é um valor grande da coordenada. Com a substituição de ψ1, ψ2 e do determinante det(W) = -2ik na Equação 7, tem-se o problema de valor inicial a ser resolvido, e , em que

(8)

A quantidade r0 é a amplitude de reflexão para o potencial degrau V(x0)H(x - x0) com H sendo a função de Heaviside. O termo é de importância para potenciais com limites assintóticos VD = V(x → +∞) ≠ 0 para levar em conta a contribuição do potencial a direita de x0 ao se iniciar a integração nesse ponto. Para o caso em que V(x → +∞) = 0, não existem partículas que se deslocam para a esquerda à direita da barreira, consequentemente, a condição inicial r(x → ∞) deve ser nula, portanto, r0 se reduz a zero para valores grandes de x.

RESULTADOS E DISCUSSÃO

Devido à importância pedagógica do potencial barreira retangular, esse potencial será considerado como um primeiro exemplo. O segundo potencial abordado será o potencial de Eckart.2929 Eckart, C.; Phys. Rev. 1930, 35, 1303. O potencial apresenta o limite assintótico, VD ≠ 0, o que possibilita mostrar a eficiência do método para uma condição de contorno diferente. Para os dois casos, considerou-se que a partícula incide da esquerda para a direita. A integração numérica foi realizada para valores de coordenada de 20 a -20 utilizando o método de Runge-Kutta quarta e quinta ordem com passo adaptativo.3030 Forsythe, G. E.; Malcolm, M. A.; Moler, C. B.; Computer Methods for Mathematical Computations, Prentice Hall: New Jersey, 1977. Os cálculos foram feitos em unidades atômicas.

Potencial barreira retangular

Um potencial contínuo que representa uma barreira pode ser estabelecido com a função,

(9)

No limite de α grande, o potencial tende ao potencial barreira retangular (Figura 1), em que Vmax é o máximo do potencial e L a largura da barreira.

Figura 1
Aproximação para o potencial barreira retangular, α = 100, Vmax = 1 e L = 5

As soluções exatas do caso limite para E > Vmax e E < Vmax, são dadas, respectivamente, por

(10)

(11)

com κ = .1515 Sanchez-Soto, L. L.; Monzona, J. J.; Barriuso, A. G.; Carinena, J. F.; Phys. Rep. 2012, 513, 191. e =

A solução numérica para o coeficiente de transmissão, T = 1 - R, em que R = |r|2, com parâmetros α = 100, Vmax = 1 e L = 5, é apresentada na Figura 2 em comparação com as soluções exatas.

Figura 2
Coeficiente de transmissão para o potencial barreira retangular, (asteriscos) solução numérica da Eq. (8), (linha tracejada) solução analítica

O erro médio cometido é menor que 0,05%, neste caso, o erro além de ser função dos parâmetros ajustáveis do método de resolução da equação diferencial que controla a precisão, é também função do valor α considerado para a aproximação. A solução correta é obtida tanto para E > Vmax, quanto para E < Vmax. Os valores de T > 0 para E < Vmax indicam o tunelamento das partículas, isto é, mesmo que as partículas tenham energia menor que a da barreira, ainda há uma probabilidade finita de essas a atravessarem, o que classicamente não seria permitido. Também, mesmo que as partículas possuam energia maior que a da barreira, E > Vmax, ainda se observa T < 1, quando classicamente se esperaria T = 1.

O potencial considerado se anula assintoticamente, tanto para valores negativos, quanto para valores positivos de x, um caso com condições de contorno diferentes pode ser ilustrado com o potencial de Eckart.2929 Eckart, C.; Phys. Rev. 1930, 35, 1303.

Potencial de Eckart

O potencial de Eckart é dado pela função,1313 Razavy, M.; Quantum Theory of Tunneling, World Scientific: New Jersey, 2003.,2929 Eckart, C.; Phys. Rev. 1930, 35, 1303.

(12)

em que a e b são parâmetros arbitrários. O potencial possui limites diferentes como função dos parâmetros e possui solução analítica para o tunelamento o que permite a comparação com os resultados numéricos. Na Figura 3 é apresentado o potencial de Eckart onde se verifica a modificação no limite VD em função da modificação dos parâmetros.

Figura 3
Potencial de Eckart para diferentes valores do parâmetro a com b = 8a

A Figura 4 apresenta os coeficientes de reflexão, R, para valores de energia tomados no intervalo, E < 4,5. Os resultados estão em comparação com a solução analítica,<

(13)

em que β = b > 1/4. e δ = , apropriada para descrever o coeficiente de reflexão para

Figura 4
Coeficiente de reflexão R para o potencial de Eckart para diferentes valores de a com b = 8a, (linha contínua) exato, (círculos) solução numérica

Como se verifica, a resolução da Equação 8 leva a valores de coeficientes de reflexão corretos. O erro médio em relação à solução analítica foi de 1,54 × 10-6%. Uma das vantagens de se utilizar o método da amplitude variável é a possibilidade de utilizar técnicas de integração com passo adaptativo, o que contorna os problemas citados por Rozman et al.1818 Rozman, M. G.; Reineker, P.; Tehver, R.; Phys. Rev. A: At., Mol., Opt. Phys. 1994, 49, 3310. no método da matriz de transferência.

CONCLUSÃO

A abordagem pelo método da amplitude variável consiste em um caminho elementar para reduzir os cálculos dos coeficientes de transmissão ou reflexão a um problema de valor inicial, utilizando uma equação diferencial de primeira ordem. A equação diferencial pode ser resolvida pelas técnicas com passo adaptativo, o que leva a uma solução precisa com um mínimo esforço computacional. Utilizando uma condição inicial adequada o método torna-se aplicável de maneira geral a potenciais localizados com diferentes condições de contorno. O método aqui apresentado pode ser adaptado para o cálculo da matriz de espalhamento e versões em mais de uma dimensão devem ser investigadas.

De simples implementação computacional, o método representa uma ferramenta atrativa para compreender os processos de tunelamento para os mais diversos potenciais. A exemplo da equação de Calogero, a equação diferencial para a amplitude de reflexão pode ser utilizada para estudar o problema inverso de recuperação do potencial a partir de dados de coeficiente de transmissão ou reflexão.

MATERIAL SUPLEMENTAR

Um exemplo de código para cálculos de coeficiente de reflexão desenvolvido em MATLAB está disponível como material suplementar em http://quimicanova.sbq.org.br, com acesso livre.

https://minio.scielo.br/documentstore/1678-7064/XPmcTr5QWpRJ4ddHvDgsLjf/420b1c32551ec1bee9ba2cf145887fea6362fd80.pdf

AGRADECIMENTOS

Gostaríamos de agradecer à CAPES, ao CNPq e à Fapemig pelo suporte financeiro.

REFERÊNCIAS

  • 1
    Gamow, G.; Z. Phys. 1928, 51, 204.
  • 2
    Gurney, R. W.; Condon, E. U.; Phys. Rev. 1929, 33, 127.
  • 3
    Fowler, R. W.; Nordheim, L.; Proc. R. Soc. London, Ser. A 1928, 119, 173.
  • 4
    Bell, R. P.; The Tunnel Effect in Chemistry, Chapman and Hall: London, 1980.
  • 5
    Johnston, H. S.; Rapp, D.; J. Am. Chem. Soc. 1961, 83, 1.
  • 6
    Levich, V. G.; Dogonadze, R. R.; Dokl. Akad. Nauk SSSR 1960, 133, 158.
  • 7
    Marcus, R. A.; J. Chem. Phys. 1966, 45, 4493.
  • 8
    Khairutdinov, R. F.; Zamaraev, K. I.; Zhdanov, V. P. Em Comprehensive Chemical Kinetics; Compton, R. G., eds.; Elsevier: Amsterdam, 1989, vol. 30.
  • 9
    Tserkis, S. T.; Moustakidis, Ch. C.; Massen, S. E.; Panos, C. P.; Phys. Lett. A 2014, 378, 497.
  • 10
    Johnson, E. A. Em Low-Dimensional semiconductor structures; Barnham, K., Vvedensky, D., eds.; Cambridge Univ. Press: Cambridge, 2001, cap. 2.
  • 11
    Powell, J. L.; Crasemann, B.; Quantum Mechanics, Addison-Wesley: Massachusetts, 1961.
  • 12
    Cohen-Tannoudji, C.; Diu, B.; Laloë, F.; Quantum Mechanics, Wiley: New York, 1977, vol. 1.
  • 13
    Razavy, M.; Quantum Theory of Tunneling, World Scientific: New Jersey, 2003.
  • 14
    Vatannia S.; Gildenblat, G.; IEEE J. Quantum Electron. 1996, 32, 1093.
  • 15
    Sanchez-Soto, L. L.; Monzona, J. J.; Barriuso, A. G.; Carinena, J. F.; Phys. Rep. 2012, 513, 191.
  • 16
    Mukherjee, S.; Karmakar, R.; Deyasi, A.; International Journal of Soft Computing and Engineering 2011, 1, 41.
  • 17
    Das, S.; Am. J. Phys. 2015, 83, 590.
  • 18
    Rozman, M. G.; Reineker, P.; Tehver, R.; Phys. Rev. A: At., Mol., Opt. Phys. 1994, 49, 3310.
  • 19
    Tikochinsky, Y.; Ann. Phys. 1977, 103, 185.
  • 20
    Viterbo, V. D.; Lemes, N. H. T.; Braga, J. P.; Rev. Bras. Ensino Fis. 2014, 36, 1310.
  • 21
    Calogero, F.; Variable phase approach to potential scattering, Academic Press: New York, 1967.
  • 22
    Braga, J. P.; Murrell, J. N.; Mol. Phys. 1984, 53, 295.
  • 23
    Lemes, N. H. T.; Braga, J. P.; Alves, M. O.; Costa, É. D'M.; J. Mol. Model. 2014, 20, 2317.
  • 24
    Kidun, O.; Fominykh, N.; Berakdar, J.; Phys. Rev. A: At., Mol., Opt. Phys. 2005, 71, 022703.
  • 25
    Pliego Jr., J. R.; Braga, J. P.; Bosanac, S. D.; Phys. Rev. A: At., Mol., Opt. Phys. 1995, 53, 342.
  • 26
    Calogero, F.; Degasperis, A. Em Studies in Mathematics and its Applications; Lions, J. L., Papanicolaou, G., Rockafeliar, R. T., Fujita, H., eds.; North-Holland: Amsterdam, 1982, vol. 13.
  • 27
    Ronveaux, A.; Am. J. Phys. 1969, 37, 135.
  • 28
    Zill, D. G.; Cullen, M. S.; Equações Diferenciais, Makron Books: São Paulo, 2007.
  • 29
    Eckart, C.; Phys. Rev. 1930, 35, 1303.
  • 30
    Forsythe, G. E.; Malcolm, M. A.; Moler, C. B.; Computer Methods for Mathematical Computations, Prentice Hall: New Jersey, 1977.

Datas de Publicação

  • Publicação nesta coleção
    Ago 2016

Histórico

  • Recebido
    22 Out 2015
  • Aceito
    28 Jan 2016
Sociedade Brasileira de Química Instituto de Química, Universidade Estadual de Campinas (Unicamp), CP6154, 13083-0970 - Campinas - SP - Brazil
E-mail: quimicanova@sbq.org.br