Bound states in a double delta potential via Fourier sine and cosine transforms

rbef Revista Brasileira de Ensino de Física Rev. Bras. Ensino Fís. 1806-1117 1806-9126 Sociedade Brasileira de Física São Paulo, SP, Brazil The problem of bound states in a double delta potential is revisited by means of Fourier sine and cosine transforms. ARTIGOS GERAIS Estados ligados em um potencial delta duplo via transformadas seno e cosseno de Fourier Bound states in a double delta potential via Fourier sine and cosine transforms A.S. de Castro1 Departamento de Física e Química, Universidade Estadual Paulista "Júlio de Mesquita Filho", Guaratinguetá, SP, Brasil RESUMO O problema de estados ligados em um potencial delta duplo é revisto com o uso do método das transformadas seno e cosseno de Fourier. Palavras-chave: duplo delta, estado ligado, transformadas seno e cosseno de Fourier. ABSTRACT The problem of bound states in a double delta potential is revisited by means of Fourier sine and cosine transforms. Keywords: double delta, bound state, Fourier sine and cosine transforms. 1. Introdução Em 1927, a quantidade δ(x) definida por foi utilizada por P.A.M. Dirac para mostrar a equivalência entre três diferentes formalismos da mecânica quântica [1]. Embora seja atualmente conhecida como função delta de Dirac, tal quantidade apareceu pela primeira vez em trabalhos de O. Heaviside já em 1895 (veja, e.g., [2]). Conquanto δ(x) não seja uma função no sentido usual da palavra, ela pode ser entendida como o limite de uma sequência de funções, e a teoria das distribuições permite que ela possa ser considerada efetivamente como uma função ordinária [3]. A utilidade da função delta estende-se muito além de seu uso original em mecânica quântica e na definição de função impulsiva. É proveitosa no cálculo de soma de séries de Fourier [4] e de integrais [5], além de ser fundamental na busca de soluções particulares de equações diferenciais não-homogêneas com o método da função de Green (veja, e.g., [6]). Outrossim, a função delta tem serventia para simular funções com valores extremamente grandes em um intervalo extremamente pequeno e para generalizar o conceito de densidade associado com quantidades contínuas para quantidades discretas, sendo a densidade de massa de um objeto pontual um dos exemplos mais elementares. Em muitas circunstâncias, a razão da substituição de funções ordinárias por funções delta é porque, via de regra, tais modelos tornam-se mais simples, e às vezes exatamente solúveis [7, 8]. De fato, as funções delta têm sido utilizadas em uma pletora de trabalhos em situações variadas em física estatística [9], eletromagnetismo [10], ótica [mar], modelos nucleares [11], física do estado sólido [12, 13], física de partículas elementares [15], teoria do espalhamento [16], simulação do comportamento de átomos e moléculas [17, 18], fotoionização [19], efeito Hall quantizado [20], física multidimensional [21], efeito Casimir [22], regularização e renormalização em teoria quântica de campos [23]. A equação de Schrödinger com um potencial constituído de uma soma de duas funções delta de Dirac, doravante denominado potencial delta duplo, tem sido usada na descrição da transferência de um núcleon de valência durante uma colisão nuclear [12] tanto quanto para para modelar as forças de troca entre os dois núcleos no íon de hidrogênio molecular [18]. A bem da verdade, os estados estacionários de uma partícula em um potencial delta duplo ocupa as páginas de muitos livros-texto [24]-[30]. Costumeiramente, os possíveis estados ligados são encontrados pela localização dos polos complexos da amplitude de espalhamento ou por meio de uma solução direta da equação de Schrödinger baseada na descontinuidade da derivada primeira da autofunção, mais a continuidade da autofunção e seu bom comportamento assintótico. Os estados ligados da equação de Schrödinger com potenciais delta também tem sido alvo de investigação tanto com a transformada de Laplace [31] quanto com a transformada exponencial de Fourier [32], e no caso de potenciais constituídos de funções delta de Dirac a equação de Schrödinger independente do tempo transmuta-se numa equação algébrica de primeira ordem para a transformada da autofunção. O uso de transformadas integrais para resolver uma equação diferencial é louvável se a equação transformada puder ser resolvida com maior simplicidade. No entanto, deve-se executar a inversão da transformada para obter a função original do problema. Esta última tarefa pode ser penosa mas é extremamente facilitada pelo uso de tabelas de integrais. Neste trabalho usamos as transformadas seno e cosseno de Fourier na busca de soluções de estados ligados da equação de Schrödinger com o potencial delta duplo. A abordagem da equação de Schrödinger com o potencial delta duplo via transformadas seno e cosseno de Fourier fornece uma aplicação adicional do método de transformadas integrais em um problema físico simples que pode ser do interesse de professores e estudantes de física matemática e mecânica quântica dos cursos de graduação em física. O procedimento adotado envolve contato com equações diferenciais e o comportamento assintótico de suas soluções, função delta de Dirac, condições de contorno, paridade e extensões simétrica e antissimétrica de autofunções, degenerescência em sistemas físicos unidimensionais, et cetera. 2. As transformadas seno e cosseno de Fourier As transformadas seno e cosseno de Fourier de ϕ(x), denotadas por ΦS(k) e ΦC(k) respectivamente, são definidas como [6]-[33] e as transformadas inversas são dadas por A existência das transformadas seno e cosseno de Fourier, e de suas inversas, é assegurada se as funções envolvidas nos integrandos forem absolutamente integráveis no intervalo [0,∞), e para isto acontecer ɸ(x)→0 quando x →∞, e também ΦS(k) e ΦC(k) devem tender a zero quando k →∞. Neste caso, o teorema de Parseval garante que Além disto, se dɸ(x)/dx for também absolutamente integrável valerão as seguintes propriedades diferenciais onde Somente quatro integrais relacionadas com as transformadas seno e cosseno de Fourier de sen kc/(k2±d2) e cos kc/(k2±d2), serão necessários na análise dos estados ligados de um potencial delta duplo. Apesar de que apenas duas integrais sejam utilizados na expressão final das autofunções. 3. Os estados ligados de um potencial delta duplo O potencial delta duplo simétrico pode ser visto como o caso limite de um potencial de poço (ou barreira) quadrado duplo como está ilustrado na Fig. 1, onde L e θ são quantidades positivas. O limite apropriado deve ser realizado impondo θ→0 e |V0|→∞ de tal forma que o produto θV0 permaneça constante. Usando as definições e a equação de Schrödinger independente do tempo para uma partícula de massa m sujeita a um potencial delta duplo simétrico pode ser escrita na forma Para estados ligados, devemos procurar autofunções que se anulam à medida que |x|→∞. Tendo em vista que a Eq. (9) é invariante sob reflexão através da origem (x→-x), se ϕ(x) satisfaz à equação de Schrödinger para um dado E, assim acontece com ϕ(-x), e portanto também satisfazem as combinações lineares ϕ(x)± ϕ(-x). Daí, assegura-se que autofunções com paridades bem definidas podem ser construídas, sendo suficiente concentrar a atenção sobre o semieixo e impor condições de contorno adicionais sobre ϕ na origem. Autofunções em todo o eixo podem ser construídas tomando combinações lineares simétricas e antissimétricas de ϕ definida no lado positivo do eixo x onde θ(x) é a função degrau de Heaviside (1 para x>0, e 0 para x<0). Estas duas autofunções linearmente independentes possuem a mesma energia, então, em princípio, existe uma dupla degenerescência. Por causa da continuidade da autofunção e sua derivada em x=0, as condições de contorno sobre ϕ na origem podem ser combinadas de duas formas distintas: a função par obedece à condição de Neumann homogênea ϕ'(0)=0, enquanto a função ímpar obedece à condição de Dirichlet homogênea ϕ(0)=0. Haja vista que ϕ(x) e dϕ(x)/dx são nulas no infinito, temos a garantia da existência de FS{d2ϕ(x)dx2} e FC{d2ϕ(x)dx2} caso existam as transformadas de ϕ(x) e dϕ(x)/dx. As transformadas seno e cosseno de Fourier da Eq. (10) resultam em equações algébricas de primeiro grau A condição de contorno na origem dita qual o tipo de transformada é mais apropriado: transformada seno de Fourier para uma autofunção ímpar e transformada cosseno de Fourier para uma autofunção par. A transformada seno (cosseno) de Fourier atende à conveniência da autofunção ímpar (par) porque a transformada inversa requer ϕ(0)=0 (ϕ'(0)=0) como pode ser depreendido da Eq. (3). Destarte, obtemos as soluções Das integrais tabeladas na Ref. [33] e transcritas no Apêndice em (A3) e (A4) concluímos sobre a inexistência de transformadas inversas no caso em que . Isto sucede porque a presença de cosκx ou sen κx faz com que ϕ(x) não seja absolutamente integrável. Entretanto, diante das integrais tabeladas expressas por (A1) e (A2), concluímos que a existência das transformadas inversas é assegurada para κ=iξ/L com ξ>0 de tal forma que Neste caso, a transformada inversa da Eq. (13) pode ser obtida usando os resultados tabelados expressos por (A1) e (A2). Assim, para ϕ(-x)=+ϕ(x) temos onde Por outro lado, para ϕ (-x)=-ϕ (x) temos onde ε(x) é a função sinal de x (+1 para x>0, e -1 para x<0), e A consistência das Eqs. (16), (18) e (19) requer que ξ seja solução das equações transcendentais A condição de quantização surgiu, além da exigência de normalizabilidade das autofunções, como uma necessidade de consistência no ajuste dos valores de ϕ (0), ϕ′(0) e ϕ(L). Há de se notar que a condição de quantização também poderia ter sido obtida pela fórmula de conexão entre dϕ/dx à direita e à esquerda de x=L. Tal fórmula pode ser avaliada pela integração da Eq. (10) numa pequena região em redor x=L, e pode ser sumarizada por Já que as funções 1+e-2ξ e 1-e-2ξ são limitadas a valores positivos ao passo que 2aξ é limitada a valores negativos quando a<0, podemos inferir que não há possibilidade de solução para estados ligados se a<0 (potencial repulsivo). Para um potencial atrativo (a>0), a natureza do espectro resultante das soluções das equações transcendentais (20) podem ser visualizadas na Fig. 2, onde constam esboços do membro direito da Eq. (20), e do lado esquerdo da Eq. (20) para a=3/2 e a=1/4. As curvas representadas por 1+e-2ξ e 2aξ sempre se interceptam. Note que 1-e-2ξ≈2ξ para |ξ|≪1 e assim a curva representada por 2aξ oscula a curva representada por 1-e-2ξ em ξ=0 quando a=1, e estas duas curvas se interceptam em algum ponto com abscissa ξ>0 se e somente se a<1. As abscissas das interseções de 1±e-2ξ e 2aξ correspondem aos valores permitidos de ξ. Pode-se depreender da Fig. 2 que sempre há uma e somente uma solução para o caso de uma autofunção simétrica, correspondendo à interseção das curvas 1+e-2ξ e 2aξ. A existência de uma adicional solução, necessariamente para o caso de uma autofunção antissimétrica, sucede tão somente quando a<1, correspondendo à interseção das curvas 1-e-2ξ e 2aξ. Em todo caso, os valores permitidos de ξ aumentam à medida que a diminui. Além disto, o valor permitido de ξ para a solução ímpar é menor que o valor permitido de ξ para a solução par. A relação entre E e ξ dada pela Eq. (14) permite concluir que a energia dos estados ligados diminui com o aumento de ξ e consequentemente a solução com autofunção par corresponde ao ubíquo estado fundamental. Relembrando a relação entre a, L e α dada pela Eq. (8), podemos concluir que a existência de uma solução ímpar, correspondendo ao estado excitado, ocorre somente se o potencial for suficientemente forte (α>ħ2/2mL). É também instrutivo observar que, para a<1, o espectro torna-se degenerado no limite ξ→∞ (E→-∞) com a→0 (α→∞). 4. Comentários finais Este trabalho apresentou uma abordagem alternativa para a busca de estados ligados do potencial delta duplo baseada nas transformações seno e cosseno de Fourier. Com essa abordagem a equação de Schrödinger independente do tempo transmutou-se em equações algébicas de primeira ordem para as transformadas da autofunção. O processo de inversão das transformadas tornou-se amigável porque as integrais envolvidas no procedimento estão presentes em tabelas, do contrário a abordagem, ainda que lícita, perderia o seu valor pedagógico. Agradecimentos O autor é grato ao CNPq pelo apoio financeiro. Um árbitro atencioso contribuiu com comentários e sugestões pertinentes para proscrever incorreções na primeira versão e facetar a presente versão deste trabalho. Apêndice - Integrais úteis As fórmulas 3.742.1, 3.742.3, 3.742.6 e 3.742.8 da Ref. [33] são as integrais úteis para nossos propósitos: com x > 0, c > 0 e d - . Referências [1] P.A.M. Dirac, Poc. R. Soc. London A 113, 621 (1927). [2] J.D. Jackson, Am. J. Phys. 76, 704 (2008). [3] L. Schwartz, Theorie des distributions (Hermann, Paris, 1950); J.P. Marchand, Distributions (North-Holland, Amsterdam, 1962); I.M. Gel'fand and G.E. Shilov, Generalized Functions (Academic, New York, 1964). [4] M.R. Spiegel, J. Appl. Phys. 23, 906 (1952). [5] M.R. Spiegel, J. Appl. Phys. 25, 1302 (1954). [6] E. Butkov, Física Matemática (LTC, Rio de Janeiro, 1988). [7] Y.N. Demkov and V.N. Ostrovskii, Zero-Range Potentials and Their Application in Atomic Physics (Plenum, New York, 1988). [8] S. Albeverio et al., Solvable Models in Quantum Mechanics (Springer-Verlag, New York, 1988). [9] A.E. Siegman, Am. J. Phys. 47, 545 (1979). [10] V. Namias, Am. J. Phys. 45, 624 (1977); D.T. Gillespie, Am. J. Phys. 51, 520 (1983); D.G. Hall, Am. J. Phys. 63, 508 (1995). [11] D. Marcuse, Theory of Dielectric Optical Waveguides (Academic, New York, 1974). [12] G. Breit, Proc. Second Conf. Reactions between Complex Nuclei, Gatlinburg, p.1 (Wiley, New York, 1960); G. Breit, Ann. Phys. (N.Y.) 34, 377 (1965). [13] R. de L. Kronig and W.G. Penney, Proc. Roy. Soc. London A 130, 499 (1931); S. Weinberg, Phys. Lett. B 251, 288 (1990); S. Weinberg, Nucl. Phys. B 363, 3 (1991). [14] S. Albeverio and R. Hoegh-Kronh, J. Oper. Theor. 6, 313 (1981); J.E. Avron et al., Phys. Rev. Lett. 72, 896 (1994). [15] E.U. Condon, Phys. Rev. 49, 459 (1936); C. Thorn, Phys. Rev. D 19, 639 (1979). [16] I.R. Lapidus, Am. J. Phys. 37, 930 (1969); I.R. Lapidus, Am. J. Phys. 37, 1064 (1969); E.Kujawski, Am. J. Phys. 39, 1248 (1971); P. Senn, Am. J. Phys. 56, 916 (1988). [17] C.M. Rosenthal, J. Chem. Phys. 55, 2474 (1971); I.R. Lapidus, Am. J. Phys. 43, 790 (1975). [18] A.A. Frost, J. Chem. Phys. 25, 1150 (1956). [19] S. Geltman, J. Phys. B 10, 831 (1977); E.J. Austin, J. Phys. B 12, 4045 (1979); H. Fearn and W.E. Lamb Jr., Phys. Rev. A 43, 2124 (1991). [20] R.E. Prange, Phys. Rev. B 23, 4802 (1981). [21] E. Demiralp and H. Beker, J. Phys. A 36, 7449 (2003); S. Giardino, Revista Brasileira de Ensino de Física 35, 3307 (2013). [22] L.A.S. Carrilho and J.A. Nogueira, Revista Brasileira de Ensino de Física 31, 2311 (2009). [23] I. Mitra et al., Am J. Phys. 66, 1101 (1998); R.M. Cavalcanti, Revista Brasileira de Ensino de Física 21, 336 (1999). [24] S. Gasiorowicz, Física Quântica (Guanabara Dois, Rio de Janeiro 1974). [25] C. Cohen-Tannoudji et al., Quantum Mechanics, V. 1 (Hermann, Paris 1977). [26] A. Galindo and R. Pascual, Quantum Mechanics I (Springer-Verlag, Berlin 1990). [27] K. Gottfried and T.-M. Yan, Quantum Mechanics: Fundamentals, 2nd. ed. (Springer, New York, 2003). [28] K. Tamvakis, Problems & Solutions in Quantum Mechanics (Cambridge University Press, Cambridge, 2005). [29] R.W. Robinett, Quantum Mechanics, 2nd. ed. (Oxford University Press, Oxford, 2006). [30] D.J. Griffiths, Mecânica Quântica, 2ª ed. (Pearson Prentice Hall, São Paulo, 2011). [31] A.S. de Castro, Revista Brasileira de Ensino de Física 34, 4301 (2012). [32] A.S. de Castro, Revista Brasileira de Ensino de Física 34, 4304 (2012). [33] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, A. Jeffrey and D. Zwillinger (eds.) (Academic Press, New York, 2007), 7th ed. Recebido em 27/8/2013 Aceito em 16/1/2014 Publicado em 11/5/2014 1 E-mail: [e-mail].b [1] P.A.M. Dirac, Poc. R. Soc. London A 113, 621 (1927). Poc. R. Soc. London A 1927 621 113 Dirac P.A.M. [2] J.D. Jackson, Am. J. Phys. 76, 704 (2008). Am. J. Phys. 2008 704 76 Jackson J.D. [3] L. Schwartz, Theorie des distributions (Hermann, Paris, 1950); Theorie des distributions Schwartz L. J.P. Marchand, Distributions (North-Holland, Amsterdam, 1962); Distributions 1962 Marchand J.P. I.M. Gel'fand and G.E. Shilov, Generalized Functions (Academic, New York, 1964). Generalized Functions 1964 Gel'fand I.M. Shilov G.E. [4] M.R. Spiegel, J. Appl. Phys. 23, 906 (1952). J. Appl. Phys. 1952 906 23 Spiegel M.R. [5] M.R. Spiegel, J. Appl. Phys. 25, 1302 (1954). J. Appl. Phys. 1954 1302 25 Spiegel M.R. [6] E. Butkov, Física Matemática (LTC, Rio de Janeiro, 1988). Física Matemática 1988 E. Butkov [7] Y.N. Demkov and V.N. Ostrovskii, Zero-Range Potentials and Their Application in Atomic Physics (Plenum, New York, 1988). Zero-Range Potentials and Their Application in Atomic Physics 1988 Demkov Y.N. Ostrovskii V.N. [8] S. Albeverio et al., Solvable Models in Quantum Mechanics (Springer-Verlag, New York, 1988). Solvable Models in Quantum Mechanics 1988 Albeverio S. [9] A.E. Siegman, Am. J. Phys. 47, 545 (1979). Am. J. Phys. 1979 545 47 Siegman A.E. [10] V. Namias, Am. J. Phys. 45, 624 (1977); Am. J. Phys. 1977 624 45 Namias V. D.T. Gillespie, Am. J. Phys. 51, 520 (1983); Am. J. Phys. 1983 520 51 Gillespie D.T. D.G. Hall, Am. J. Phys. 63, 508 (1995). Am. J. Phys. 1995 508 63 Hall D.G. [11] D. Marcuse, Theory of Dielectric Optical Waveguides (Academic, New York, 1974). Theory of Dielectric Optical Waveguides 1974 Marcuse D. [12] G. Breit, Proc. Second Conf. Reactions between Complex Nuclei, Gatlinburg, p.1 (Wiley, New York, 1960); Proc. Second Conf. Reactions between Complex Nuclei, Gatlinburg 1960 1 Breit G. G. Breit, Ann. Phys. (N.Y.) 34, 377 (1965). Ann. Phys. (N.Y.) 1965 377 34 Breit G. [13] R. de L. Kronig and W.G. Penney, Proc. Roy. Soc. London A 130, 499 (1931); Proc. Roy. Soc. London A 1931 499 130 Kronig R. de L. Penney W.G. S. Weinberg, Phys. Lett. B 251, 288 (1990); Phys. Lett. B 1990 288 251 Weinberg S. S. Weinberg, Nucl. Phys. B 363, 3 (1991). Nucl. Phys. B 1991 3 363 Weinberg S. [14] S. Albeverio and R. Hoegh-Kronh, J. Oper. Theor. 6, 313 (1981); J. Oper. Theor. 1981 313 6 Albeverio S. Hoegh-Kronh R. J.E. Avron et al., Phys. Rev. Lett. 72, 896 (1994). Phys. Rev. Lett. 1994 896 72 Avron J.E. [15] E.U. Condon, Phys. Rev. 49, 459 (1936); Phys. Rev. 1936 459 49 Condon E.U. C. Thorn, Phys. Rev. D 19, 639 (1979). Phys. Rev. D 1979 639 19 Thorn C. [16] I.R. Lapidus, Am. J. Phys. 37, 930 (1969); Am. J. Phys. 1969 930 37 Lapidus I.R. I.R. Lapidus, Am. J. Phys. 37, 1064 (1969); Am. J. Phys. 1969 1064 37 Lapidus I.R. E.Kujawski, Am. J. Phys. 39, 1248 (1971); Am. J. Phys. 1971 1248 39 Kujawski E. P. Senn, Am. J. Phys. 56, 916 (1988). Am. J. Phys. 1988 916 56 Senn P. [17] C.M. Rosenthal, J. Chem. Phys. 55, 2474 (1971); J. Chem. Phys. 1971 2474 55 Rosenthal C.M. I.R. Lapidus, Am. J. Phys. 43, 790 (1975). Am. J. Phys. 1975 790 43 Lapidus I.R. [18] A.A. Frost, J. Chem. Phys. 25, 1150 (1956). J. Chem. Phys. 1956 1150 25 Frost A.A. [19] S. Geltman, J. Phys. B 10, 831 (1977); J. Phys. B 1977 831 10 Geltman S. E.J. Austin, J. Phys. B 12, 4045 (1979); J. Phys. B 1979 4045 12 Austin E.J. H. Fearn and W.E. Lamb Jr., Phys. Rev. A 43, 2124 (1991). Phys. Rev. A 1991 2124 43 Fearn H. Lamb Jr. W.E. [20] R.E. Prange, Phys. Rev. B 23, 4802 (1981). Phys. Rev. B 1981 4802 23 Prange R.E. [21] E. Demiralp and H. Beker, J. Phys. A 36, 7449 (2003); J. Phys. A 2003 7449 36 Demiralp E. Beker H. S. Giardino, Revista Brasileira de Ensino de Física 35, 3307 (2013). Revista Brasileira de Ensino de Física 2013 3307 35 Giardino S. [22] L.A.S. Carrilho and J.A. Nogueira, Revista Brasileira de Ensino de Física 31, 2311 (2009). Revista Brasileira de Ensino de Física 2009 2311 31 Carrilho L.A.S. Nogueira J.A. [23] I. Mitra et al., Am J. Phys. 66, 1101 (1998); Am J. Phys. 1998 1101 66 Mitra I. R.M. Cavalcanti, Revista Brasileira de Ensino de Física 21, 336 (1999). Revista Brasileira de Ensino de Física 1999 336 21 Cavalcanti R.M. [24] S. Gasiorowicz, Física Quântica (Guanabara Dois, Rio de Janeiro 1974). Física Quântica 1974 Gasiorowicz S. [25] C. Cohen-Tannoudji et al., Quantum Mechanics, V. 1 (Hermann, Paris 1977). Quantum Mechanics 1977 1 Cohen-Tannoudji C. [26] A. Galindo and R. Pascual, Quantum Mechanics I (Springer-Verlag, Berlin 1990). Quantum Mechanics I 1990 Galindo A. Pascual R. [27] K. Gottfried and T.-M. Yan, Quantum Mechanics: Fundamentals, 2nd. ed. (Springer, New York, 2003). Quantum Mechanics: Fundamentals, 2nd. ed 2003 Gottfried K. Yan T.-M. [28] K. Tamvakis, Problems & Solutions in Quantum Mechanics (Cambridge University Press, Cambridge, 2005). Problems & Solutions in Quantum Mechanics 2005 Tamvakis K. [29] R.W. Robinett, Quantum Mechanics, 2nd. ed. (Oxford University Press, Oxford, 2006). Quantum Mechanics, 2nd. ed. 2006 Robinett R.W. [30] D.J. Griffiths, Mecânica Quântica, 2ª ed. (Pearson Prentice Hall, São Paulo, 2011). Mecânica Quântica 2011 Griffiths D.J. [31] A.S. de Castro, Revista Brasileira de Ensino de Física 34, 4301 (2012). Revista Brasileira de Ensino de Física 2012 4301 34 Castro A.S. de [32] A.S. de Castro, Revista Brasileira de Ensino de Física 34, 4304 (2012). Revista Brasileira de Ensino de Física 2012 4304 34 Castro A.S. de [33] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, A. Jeffrey and D. Zwillinger (eds.) (Academic Press, New York, 2007), 7th ed. 2007 Gradshteyn I.S. Ryzhik I.M. Jeffrey A. Zwillinger D.