Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger...

58
Novos Metodos Teoricos pa~ Estudos de Espalhamento de Eletrons e Posi- trons envolvendo Processos Multica- Ililliiiil 8·2·001109 Tese apresentada ao Instituto de FTsica de sio Carlos- USP; para obten~io do tTtulo de Doutor em Ciencias FTsica Basica. IIFSC SERV;c;O 0::-: BIBUOTECA E INFGF<\ ;c.,C;i-O

Transcript of Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger...

Page 1: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Novos Metodos Teoricos pa~ Estudosde Espalhamento de Eletrons e Posi-trons envolvendo Processos Multica-

Ililliiiil8·2·001109

Tese apresentada ao Instituto de FTsicade sio Carlos- USP; para obten~io dotTtulo de Doutor em Ciencias FTsicaBasica.

IIFSC SERV;c;O 0::-: BIBUOTECA EINFGF<\ ;c.,C;i-O

Page 2: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Av. Dr. Carlos Botelho, 1465CEP 13560-250 - Sao Carlos - SPBrasil

Fone (0162)72-6222Fax (0162) 72-2218

MEMBROS DA COMISSAO JULGADORA DA TESE DE DOUTORADO DE OTIL LARA

APRESENTADA AO INSTITUTO DE FiSICA DE SAD CARLOS, UNIVERSIDADE DE SAOPAULO, EM 01/09/1995.

-~----- -pro1r. Lee Mu-Tao ,

;;~;:~~~----- 0

--------~-~---- "-Prof. Dr. Mila

Page 3: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

AGRADECIMENTOS

ao prol. Dr. Lee Mu- Tao, meus sinceros agradecimentos pela orientação

segura e amiga com que me conduziu durante todo este trabalho, pelo convívio entre

nossas famílias durante todos estes anos, tendo ao seu lado a prof. DI. lone Iga,

sempre cordial, amiga e solidária. É um privilégio tê-Ios como amigos.

aos colegas e amigos Sérgio E. Michelin, Milton Fujimoto e Rogério pelas

proveitosas discussões e também pelas horas descontraídas que passamos juntos,

meus agradecimentos.

aos amigos de todas as horas: Márcio e Karfa, Otávio Aggio, Helena Bragato,

Mário Galeano, Agda e João Carfos, Mariazinha, Ana, Struckel, Ronaldo Barbieri,

Paulo Eduardo, Luíz Eugênio e Regina, e tantos outros.

à minha esposa Maria Cristina e aos nossos Guilherme, Tiago e Maria Luiza,

pelo apoiO carinhoso e paciente durante todo este tempo, que teve seus momentos

de dificuldades e também de muita alegria.

ao Instituto de Física de São Carfos, onde fui muito bem recebido, e tive o

prazer de fazer grandes amizades.

aos colegas e amigos do Departamento de Física da Universidade de Brasília,

em particular ao prof. Dr. Oyanarte Portilho.

ao prof. Dr. Charles A. Dauwe, amigo que só a distancia nos separa.

às meninas da biblioteca do Instituto de Física de São Carfos, sempre

solícitas, competentes e gentís.

à Wladerez e suas companheiras de trabalho na secção de alunos da Pós­

Graduação, pelo carinho com que todos, a que alí se dirigem, são tratados.

à CAPES, FINEP-PADCT, CNPq e FAPESP cujo auxílio financeiro e material

foram fundamentais para o desenvolvimento deste trabalho.

Page 4: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Dedico este trabalhoa minha esposa MariaCristina e aos nossos

Guilherme, Tiago eMaria Luiza.

Page 5: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

íNDICE

RESUMO

ABSTRACT ii

CAPíTULO 1 - INTRODUÇÃO 1

CAPíTULO 2 - TEORIAS 5

2.1 - FORMALlSMO DE CLOSE-COUPLlNG 62.2 - FORMALlSMO DE TWO-POTENTIAL E APROXIMAÇÕES DE

ACOPLAMENTOSFRACOS 72.3 - SÉRIES DE BORN 92.'4 - A APROXIMAÇÃO DE ONDA DISTORCIDA 102.5 - UNITARIZAÇÃO 10

CAPíTULO 3 - MÉTODOS PROPOSTOS 12

3.1 - TWO-POTENTIAL MULTICANAL3.2 - O FUNCIONAL - t3.3 - O MÉTODO DAS FRAÇÕES CONTíNUAS

121417

CAPíTULO 4 - CÁLCULOS E RESULTADOS 24

4.1 - POTENCIAL DE HUCK 244.2 - ESTUDO DE ESPALHAMENTO e+~H AO NíVEL DE DOIS CANAIS

ACOPLADOS (1S A 2S) 304.3 - ESTUDO DE ESPALHAMENTO e+~H AO NíVEL DE CINCO CANAIS

ACOPLADOS (1S A 2P) 36

CAPíTULO 5 - CONCLUSÃO 40

APÊNDICE

A - CÁLCULO DO TERMO (~f IV2GjV2/ ~j ) 42B - POTENCIAIS DE INTERAÇÃO 48

REFERÊNCIAS 49

Page 6: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

RESUMO

Dada a relevância de efeitos multicanais na interação de elétrons de baixa

energia com átomos e moléculas, a consideração destes, nos cálculos de secção de

choque, se toma indispensável. O objetivo do presente trabalho visa encontrar um

método teórico aplicável aos espalhamentos multicanal de elétrons por moléculas,

que deve reunir qualidades, tais como a exatidão de resultados calculados e a efi­

ciência computacional. Foram propostos dois novos formalismos teóricos, que são,

respectivamente, o Método de "Two-Potential" Multicanal e o Método do Funcional-C

Multicanal. Estes dois métodos foram desenvolvidos com base no Método de Ondas

Distorcidas e no Princípio Variacional de Schwinger. Por outro lado, um método ite­

rativo, conhecido como Método das Frações Contínuas, proposto por Horacek e Sa­

sakawa, em 1983, foi também adaptado, pela primeira vez, para espalhamentos

multicanais.

Estes três novos métodos foram testados para resolver problemas de es­

palhamento multicanais com o modelo de potencial de Huck, e também para o espa­

Ihamento de pósitron por átomos de hidrogênio, a nível de dois e de cinco canais.

Nossos resultados mostraram que, entre estes três métodos, o de Frações

Contínuas apresenta o melhor desempenho, tanto na eficiência computacional, quan­

to na exatidão dos resultados obtidos. Nosso estudo indica o grande potencial do

Método de Frações Contínuas para o estudo de espalhamento de elétrons por molé­

culas. A generalização deste método para cálculos realísticos está em andamento.

Page 7: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

ABSTRACT

with Huck model potential as well as to treat the positron scattering by atomic

hydrogen in the two-channel-couplingand five-channel-coupling levels.

within the numerical precision of the method. Thus, the present study indicates the

potential of this method and some realistic studies of electron-molecule scattering

using the Method of Continued Fractions are now in progress.

Page 8: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

1

CAPÍTULO 1-INTRODUÇÃO

Colisões elétron-átomo e elétron-molécula são processos básicos em física

atômica. O conhecimento das diversas secções de choque -provenientes--dos--dife­

rentes processos de espalhamento, tanto elástico quanto inelástico. é de

fundamental importância para o entendimento da estrutura eletrônica dos átomos e

moléculas, bem como dos efeitos dinâmicos que afetam os processos de colisão.

Secçães de choque de ionização por impacto de elétrons são dados importantes na

modelagem da estrutura e na dinâmica de plasmas a altas temperaturas, que

ocorrem tanto naturalmente nas estrelas e na atmosfera1-3, inclusive de outros

planetas, quanto artificialmente em aparelhos de fusão de plasmas, e também plas­

mas de baixa temperatura (plasmas frios), que se apresentam como tecnologia4,5 de

ponta neste final de século. Esses avanços tecnológicos que têm ocorrido nos úl­

timos anos gerou uma necessidade, cada vez maior, na precisão das secçóes de

choque em muitas áreas de estudos fundamentais e aplicados.

Apesar de sua importância, medidas de secções de choque são muito difíceis

Page 9: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

de serem realizadas. Existem muito poucos resultados pUblicados6,7 na literatura e

obtenyao de sec96es de choque de processos diffceis, ou mesmo impossfveis, de se

medir experimentalmente, como 0 espalhamento de eletrons por especies quimi-

Nos ultimos dez anos varias teorias multicanais foram propostas, como 0 Metodo

Variacional de Schwinger Multicanal (SMC)8-12, Metodo Variacional de Kohn 13-15e

matriz_R16-19.Estes metodos apresentam, em principio, uma base teorica muito

ab-initio ja realizados, que incluem explicitamente todas as intera¢es, fica ram

restritos a poucos canais abertos17,18,20,e tambem restritos a poucas moleculas, as

pequenas. Muito recentemente, Lima et a121,22extenderam 0 SMC para 0 calculo de

senvolvidos para calculos em estado sOlido. Por outro lado, as aproxima¢es de or-

dem menores, tais como 0 metoda de ondas distorcidas23,24que foi bastante aplicado

eletrons com atomos, deu resultados bastante satisfatorios quando a energia do

eletron incidente era alta. No infcio dos anos 80 0 metoda das ondas distorcidas

tambem foi aplicado para 0 calculo de excitagOes eletronicas de moleculaS25-35,para

Page 10: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

3

particular da excitação X1L; 4C3nu do H2, embora as secções de choque calculadas

com o método das ondas distorcidas estejam muito semelhantes às calculadas pelo

Método Multicanal de Schwinge~8, ao nível de dois canais, quantitativamente ambas

as teorias superestimam as secções de choque, por um fator de 3 a 4, quando com­

paradas com os dados experimentais36 e também com os cálculos realizados com o

Método Variacional de Kohn2o ao nível de 5 canais, e com a matriz-R envolvendo 7

canais 17,18,evidenciando assim, os efeitos de acoplamento multicanal. Na tentativa

de incorporar o acoplamento multicanal, Silva Jr37 introduziu o método da teoria de

muitos corpos de primeira ordem unitarizado para o espalhamento de elétrons por

átomos de He, na faixa de energia de 29 a 80 eV, envolvendo todos os estados do

He com número quântico principal n<3. Em alguns casos, por exemplo, para a ex­

citação de 11S a 23S, a introdução do acoplamento por unitarização parece ter me­

lhorado a forma da curva de secção de choque, e já em outros casos não se re­

gistrou o mesmo efeito.

Na tentativa de entendermos os efeitos de acoplamento multicanal, realizamos

um estudo sistemático para a colisão pésitron + hidrogênio, a nível de 2 e de 5 ca­

nais. Neste estudo, foram aplicados diversos métodos teóricos, desde o mais sim­

ples, que é a primeira aproximação de Bom (PAB) , até um método ab-initio recen­

temente adaptado por nós para estudos de colisões multicanais, conhecido como

Método das Frações Contínuas (MFC)38. Comparamos a sistemática entre os resul­

tados obtidos por diversos métodos do presente estudo, assim como o custo com­

putacional de cada método, o que vai nos possibilitar escolher um método teórico

mais conveniente para cálculos realísticos envolvendo elétrons e moléculas.

Como o alvo mais simples e também o mais estudado é o hidrogênio atômico,

elejemos a colisão e+ 4 H como o objetivo deste estudo.

Page 11: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

No Capitulo 2 apresentamos os formalismos teoricos dos diversos metodos.

No Capitulo 3, dois novos metodos, Two-Potencial Multicanal (MTPM) e Metodo das

FrayOes Continuas (MCF), 580 mostrados com detalhes. Os resultados deste estudo

sao discutidos no Capitulo 4, e as conclusOes 580 apresentadas no Capitulo 5.

Page 12: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

CAPiTULO 2 - TEOR/AS

Neste capitulo, apresentarnos os forrnalisrnos te6ricos que forarn utilizados

neste trabalho para terrnos de cornparac;ao. Ernbora 0 nosso objetivo neste estudo

seja 0 espalharnento de p6sitron por atornos de hidrogenio, os forrnalisrnos forarn

derivados para 0 espalharnento de eletrons por urn atorno ou rnolecula, onde 0 efeito

de troca devido a indistinguibilidade entre 0 eletron espalhado e os eletrons do alvo e

considerada. Para espalharnento de p6sitrons por atornos ou rnoleculas esta

indistinguibilidade nao existe, portanto, 0 efeito de troca entre 0 projetil e os eletrons

do alvo nao precisa ser considerado.

Page 13: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

2.1 FORMALISMO DE CLOSE - COUPLING

A equayao de Schroedinger para a intera980 de um eletron com uma molecula

co

(v2 +k;)r:Cfo) = IUnmCfo)FmCfo)m=O

Page 14: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

sendo En a energia correspondente ao estado <I>n' Se k; > 0 0 canal e aberto, caso

contrario 0 canal e fechado. Os elementos da matriz do potencial sac dados por:

Na pratica, 0 somatorio na eq.(2.3) e truncado para poucos canais abertos. Como os

efeitos de acoplamento intercanal e levado explicitamente em conta, 0 metodo e

conhecido como metoda de acoplamento forte (close-coupling). Em geral, 0 estudo

de espalhamento de eletrons por moleculas com N canais retidos segue a resolugao

de N equa¢es integro-diferenciais do tipo da equagao (2.3), que exige urn grande

esforyo computacional. De fato, os metodos existentes tais como 0 Metodo da

matriz-R, SMC e Metodo Variacional de Khon pertencem a c1asse de acoplamento

forte, embora utilizem maneiras diferentes para a resolugao de equa¢es acopladas.

2.2 FORMALISMO DE TWO-POTENTIAL E

APROXIMACOES DE ACOPLAMENTOS FRACOS

A eq. (2.1), que descreve a interagao de urn eletron com uma molecula, pode

ser convertida em uma equagao integral

conhecida como equagao de lippmann-Schwinger, onde <I>e a solugao da equagao

Page 15: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

G 1_0- H.-E

o

e 0 operador de Green da particula livre. Fisicamente, et> corresponde a soluyao

quando nao existe a interayao entre 0 projetil e 0 alvo, ou seja, v=o.o potencial V pode ser dividido arbitrariamente em duas partes, ou seja,

1onde G1 = - -- e 0 novo operador de Green da onda distorcida.

H1-E

A funyao de onda distorcida (pelo potencial V1) e a soluyao da equayao de Schroe-

dinger

Agora, a matriz Tfi de transiyao completa podera ser escrita como

Page 16: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

e no formalismo de two-potential, Tfi sera

Tfi = (4)f I~/~)+ (~f /V;/lP;)

2.3 SERIES DE BORN

Agora, se tomarmos V1=O, teremos em consequencia que V=V2, G1=GO,

H1=Ho e ~= 4> e a equayao (2.10) reduzir-se-a a equayao (2.6). Substituindo a

equayao (2.6) na equayao (2.13), sucessivamente, entao a matriz-T sera

representada por uma sarie

que a conhecida como sarie de Born. Quando esta serie a truncada no primeiro

termo, teremos

que a a matriz-T na primeira aproximayao de Born (PAB).

Oa mesma forma, a matriz-T na segunda aproxima980 de Born (SAB) a dada

Page 17: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

2.4 A APROXIMACAO DE ONDA DISTORCIDA

Escolhemos ~ '* 0 de tal forma que a eq. (2.11) apresenta a solugao exata ~.

A substituigao sucessiva da equagao (2.10) na equagao (2.14) nos fornece uma se-

rie semelhante a serie de Born, conhecida como serie de ondas distorcidas

A matriz Tfi, na aproximagao de ondas distorcidas de varias ordens sao dadas por

2.5 UNITARIZACAO

Nas aproxima¢es de Born e de ondas distorcidas, anteriormente discutidas,

os elementos de matriz- T sac calculados para cada transigao, separadamente. Desta

forma, a matriz-S de espalhamento global (incluindo todos os canais abertos) nao e

unitaria. Entretanto, nas express5es anteriores, se tomarmos apenas a parte prin-

cipal da fungao de onda, como tambem a parte principal da fun~o de Green,

Page 18: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

s= l+iKl-iK

T=- 2iKl-iK

Page 19: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

CAPITUL03-METODOSPROPOSTOS

Neste capitulo, apresentamos os dois novos formalismos te6ricos que foram

utilizados neste trabalho, que sao, respectivamente, Two-Potential Multicanal

(MTPM) e 0 Metodo das Fray5es Contfnuas (MCF), tratados com detalhes, bem

como 0 Metodo do Funcional-C Multicanal que e uma extensao imediata do MTPM.

Nesta secy8o, apresentaremos urn novo formalismo te6rico, desenvolvido por

nosso grupo, denominado de Metodo "Two-Potential" Multicanal (MTPM). 0 MTPM euma combinayao da aproximayao de ondas distorcidas e do metoda variacional de

Schwinger. Os acoplamentos de canais sao introduzidos naturalmente e a matriz-Sobtida e unitaria. Reescrevendo a equayao (2.14) utilizando somente os valores

principais das funyOes, obteremos a matriz-K para uma determinada transi~o, quepode ser escrita de tres formas diferentes

Page 20: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

onde '¥: ('¥}) e a funyao de onda tentativa

de onda, isto e, 4:T] = 0

Page 21: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Nestes termos, a matriz Tfi de transiyao toma a forma

Tfi = (<1>f \v;\ ~;)+2:(~f \V2\~m)(~mIV2- V2G1V21~nt (~nIV2\~;) (3.6)m,n

sendo ;m a fun~o de onda distorcida no canal de espalhamento m.

Finalmente, fazendo V1=O, 0 que implica em que V2=V1, ;=<1>, G1=Go, teremos

Tfi = 2:{<1>f/V/m)(m/V -VGoV/nfl(n/VI<1>;)m,n

que nos da a aproxima~o de Born-Schwinger. Na equa~o (3.6), 0 termo

( ;m /V;G1V;/ ';n) e 0 que exige maior esfor~ computacional. 0 calcule deste termo e os

-3.2 0 FUNCIONAL- C

(MFCM). 0 funcional-C foi introduzido inicialmente por Takatsuka e McKoy, para

espalhamento de eletrons per atomos e moleculas, onde as fun¢es de base utili-

zadas sac do tipo quadraticamente integraveis (L2-functions). Mais tarde, Lee et al

introduziram a versao do funcional-C iterativo e demonstraram que 0 funcional-C

converge em 12 ordem mais rapido que 0 Metodo Variacional de Schwinger normal,

isto e, a solu~o da 12 itera~o do funcional-C corresponde aproximadamente a solu-

~o da 22 itera~o do metoda variacional de Schwinger, e assim por diante. Uma ou-

tra vantagem do funcional-C e que este nao aumenta em muito os esforyos com-

Page 22: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

putacionais em comparayao com 0 MTPM, desde que as integrais necessarias sejam

as mesmas. A seguir, derivaremos resumidamente 0 metodo do funcional-C multi-

canal, baseado na combina980 do metoda de ondas distorcidas e do funcional-C

normal.

Come~mos pela equa~o integral

A parte da funyao de onda que corresponde a parte espalhada da funyao-C I edefinida como sendo

e utilizando as equacoes (3.8) e (3.1a), a matriz-K do espalhamento podera ser re-

escrita, em termos de C como

E tamoom possivel demonstrar as seguintes identidades

(;f I~Gl~I~)=(Cfl~Gl~l;i)= (Cf I~-~Gl~ICi) (3.11)

para a funyao C. Desta forma podemos definir urn funcional variacional associado

com a equayao (2.33), como

[x] =(c; I~Gl~l<I>i)+(<I>f I~Gl~lc/}-(c;I~-~Gl~IC/) (3.12)

Page 23: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

onde C: (C~) sao as func;5es-C tentativas, associadas com os canais de espalha-

qX] =0ain

(qnjV;G1V;/q;)- 2:(qnjV; -V;G1V;jqm)dm =0m

Page 24: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

17

onde à é uma matriz de dimensão 1xN, com elementos A = (~nlV;G2V;I~t), B é uma

matriz quadrada de dimensão NxN com elemento B"" = (~n Iv; - V2G1V; I~n)' e a é a

matriz da solução da equação (3.16). Assim, podemos escrever

(3.18)

Substituindo a equação (3.17) na equação (3.10), obteremos a matriz-Kfj de

espalhamento na forma

onde

x = L(Çf I~Gl~1 ~)(çn I~-~Gl~1çmfl(çml~Gl~I~)11m

(3.19)

(3.20)

Novamente, tomando V,=O, V2=V e ç=<1>,obteremos o MFCM-Bom, cuja

matriz-K de espalhamento é dada por

K~-B = (<I>IVIcI»+ (cI>IVGoVIcI»+ L (cI> f IVGov/<I>n)(cI>n Iv - VGov/cI>mfl (cI>mIVGoVIcI>;)11m

(3.21)

3.3 MÉTODO DAS FRAÇÕES CONTíNUAS

A resolução da equação de Lippmann-Schwinger para o espalhamento de

elétrons por átomos, utilizando o Método das Frações Contínuas (MCF), foi proposto

inicialmente por Horacek e Sasakawa38, em 1.983. Naquele trabalho, estes autores

Page 25: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

estatico e de estatico-troca, e verificaram que 0 MCF e urn metoda muito eficiente

para tratar problemas de eletron-atomos. Desde entao, varias versoes extendidas41-44

e varios testes numericos38,41.45,46, utilizando potenciais modelo, tanto local quanto

MCF e rapida, mesmo para potencial de interayao forte, 0 que sugere a aplicayao do

MCF para estudos de espalhamento eletron-atomo e eletron-molecula. Entretanto,

Neste trabalho, 0 MCF e utilizado pela primeira vez para estudos de

espalhamento de p6sitron por atomo de hidrogenio, ao nivel de 2 canais e de 5

canais e, os resultados obtidos sao extremamente animadores. A seguir mostramos,

resumidamente, 0 metodo para 0 espalhamento par urn potencial.

onde Go denota a funyao de Green da onda plana, dada em termos da funyao de

Ricatti-Bessel u,(kr)[v,(kr)]por:

p/r'>rp/r'<r

Nosso prop6sito e calcular a matriz de espalhamento K definida por:

Page 26: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

1K=--tgo

k

,

Suponharnos que a prirneira aproxirna~o de Born para K, ou seja, (et>lv\et»,

nao seja igual a zero. Pode ser definido urn potencial ja Uenfraquecido" \11) por:

lJI = et>+G vi et» (et> IVI cI» + G V(l) IlJI)o (et>IVllJl) 0

lJI = et>+(1- G V(I))-I G VIet» (et>/V/lJI)o 0 (et>IVIet»

onde fizernos uso da propriedade ortogonal de \11) em rela~o a et>,

V(I)\et»= 0,

(et>/V(I) = 0,

et>1= GoVet>

lJI1= (I - GOV(l)t et>1

(3.30a)

(3.30b)

Page 27: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Multiplicando (<I>\V pela esquerda em ambos os lados da equa<;ao (3.31),

Com estas equag6es terminamos 0 primeiro passo.

No proximo passo, seguimos urn procedimento semelhante. Da equa980

que e do mesmo tipo que aquelas equay5es com que iniciamos 0 primeiro passo, ou

seja, a equa980 (3.23). Portanto, podemos seguir 0 mesmo procedimento do pri-

meiro passo. Analogamente as equay5es (3.26), (3.30a) e (3.30b), definimos y.,2) , <1>2

e '¥2 par:

V(2) = V(l) _ V(I)/<I>I)(<I>I/V(I)(<I>l/V(l)/<I>I)

<1>2 = GOV(I)<I>1

'¥2 = (1- GoV(2)t <1>1

Page 28: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

o potencial \/2) e ortogonal a ambos <I> e <I>I'

V(2)\<I»= V(2)\<I>I) = 0

(<I>lv(2) = (<I>I lv(2) = 0

(<I> /V(N)/<I> )'I' -<1>+'1' N N

N - N N+l (<I>N /V(N)/<I>N )-(<I>N /V(N)/'PN+1)

<I>N+I = GOV(N)<I>N+I

'PN+1 = <I>N+I + GOV(N+I)'PN+1

(3.38b)

(3.38c)

o operador \/N) e ortogonal a <I>, <I>1' <I>2' ... ,<I>N_I

v(N)I<I>i) = 0

(<I>i IV(N) = 0

Portanto, esperamos que 0 operador \/N) venha a ficar cada vez mais frace

com 0 aumento de N. Como urn resultado, seremos capazes de parar 0 proce-

dimento ap6s certo numero de passes M. Desprezando 0 termo GoV(M)'P M' podemos

colocar:

Page 29: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Na pratica, 0 calculo da amplitude de espalhamento e mais simples que a

funyao de onda. Em virtude das equar;5es (3.38b) e (3.39), obtemos as seguintes

(<t>N !V(N)!<t>N+2) = (<t>N+llv(N+I)!<t>N+I)

(<t>N /v(N)I<t>N+M) = 0

K. = (<t>. Iv(l-l) I'P. )I I-I I

Se usamos as equar;5es (3.38) e (3.41), podemos expressar K; na forma de uma

( I (i)1 )2K. = (<t>. \V(i-l)\<t>.)+ <t>i V <t>i

, ,-1 I (<t>./v{t)I<t>.)- K.I , 1+1

A amplitude K1 e dada pela equayao (3.44), enquanto que ela esta relacionada aamplitude de espalhamento por:

Page 30: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Como 0 nDcleodo metodo proposto e a equayao (3.44), devemos chamar este

metodo de Metodo das Frayees Contfnuas (MCF). Uma vez obtida K1, a funyao de

onda'P e calculada com a ajuda de 'Pi pela equayao (3.32)

Em suma, comec;amoscom i=1. A amplitude parcial K1 e calculada colocando-

Armazenamos as duas integrais (<I>/V/<I>I) e (<I>I/V(I)/<I>I) e a funyao <I>1 na memoria.

Para i=N(>2), calculamos a fun~o <I>N' as duas integrais (<I>N_Ilv(N-O/<I>N) e

(<I> N /V(N) / <I> N ). Armazenamos estas na memoria e apagamos <I> N -I. Colocamos

procedimento iterativo. Levamos este valor de K1 na equayao (3.45) para obter a

amplitude de espalhamento. Tendo os valores convergentes de KN, ... , K1, e fazendo

uso da equa980 (3.43) nas equay5es (3.37a), (3.38a) e (3.46) podemos calcular a

fun980 de onda.

Page 31: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

CAP;TUL04 - CALCULOS E RESULTADOS

(MTPM), do Funcional-C Multicanal (MFCM) e das Fra¢es Continuas (MCF) para 0

estudo do espalhamento de eletrons por urn potencial modelo de dois canais, co-

nhecido como potencial de Huck47.

o potencial de Huck e urn potencial multicanal com a forma:

_ {c( 1- oij) para r 5. a " =V;. - I,J 1, 2lJ 0 parar>a

Page 32: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

intensidades mudando-se a constante C. Pela sua simplicidade, este modele tern

sido frequentemente utilizado para testar metodos aproximados42,48 e para comparar

sendo que I'P;) sac as fun~s ortonormais dos estados do alvo,

No nosso estudo, estabelecemos V1 = 0, entao V2 = V, e isto implica em que

Page 33: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

ou seja, a fun980 de Green distorcida pelo potencial V1 e substitufda pela fun980 de

Green da partfcula livre, Go, e a onda distorcida ~m e substitufda pela onda plana 8m

2

Go(~,Fz;~' ,r;) = L:Sm(rl ,r<:t'Jr; ,r»m=!

Sm(lj ,rJ= k;,,~sen(kmr)\ fPm)Cm(r1,rJ = k;"~ cos(kmr)\ fPm)

Com rela980 ao metoda das frag5es contfnuas multicanal (MCFM), a matriz-K

Para podermos comparar os nossos resultados com os dados disponfveis na

Iiteratura, escolhemos os mesmos parametros utilizados por outros autores, i.e.,

Page 34: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Tabela 1 - Evolu980 dos resultados obtidos para os elementos da matriz-K, a dois

canais, utilizando-se do potencial de Huck, e pelo metodo MCF, ap6s cada itera980.

itera9800 itera9801 iterayao 2 itera9803

k11 0.000000000+00 0.243721310+00 0.243933070+00 0.243933100+00

C2 = 2.0 k12 -0.293854420+00 -0.361050800+00 -0.361165730+00 -0.361165750+00

k220ooooo0+00ס0.0 0.976806600-01 0.977888810-01 0.977888950-01

k11 0.000000000+00 0.636256960+00 0.638064250+00 0.638064760+00

C2 = 4.0 k12 -0.415572910+00 -0.663633240+00 -0.664757700+00 -0.664758020+00

k120.000000000+00 0.255004380+00 0.255851360+00 0.255851590+00

k11 OO0000+00סס0.0 0.137368600+01 0.138385270+01 0.138385710+01

C2 = 6.0 k12 -0.508970780+00 -0.116483700+01 -0.117137010+01 -0.117137300+01

k220.000000000+00 0.550557340+00 0.555032720+00 0.555034620+00

k11 0.000000000+00 0.326634440+01 0.333443020+01 0.333447120+01

C2 = 8.0 k12 -0.587708840+00 -0.238830560+01 -0.243209050+01 -0.243211680+01

k220ooooo0+00ס0.0 0.130911280+01 0.133769350+01 0.133771050+01

k11 0ooooo0+00ס0.0 0.188323800+02 0.217681700+02 0.217707240+02

C2 = 10.0 k12 -0.657078460+00 -0.122627980+02 -0.141292750+02 -0.141308970+02

k220.000000000+00 0.754779870+01 0.873502260+01 0.873605320+01

k11 0.000000000+00 -0.865268870+01 -0.808738020+01 -0.808689500+01

C2 = 12.0 k12 -0.719793390+00 0.512092010+01 0.476739510+01 0.476709050+01

k220ooooo0+00ס0.0 -0.346789680+01 -0.324607090+01 -0.324587920+01

Na Tabela 2, apresentamos os elementos de matrizes- K, para 4 processos de

espalhamentos, calculadas com MTPM e MCF, utilizando-se do potencial de Huck,

de dois canais. Para C2 = 10, os resultados exatos SaDtambem mostrados.

Page 35: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Tabela 2 - Elementos de matrizes- K, calculadas utilizando-se do potencial de Huck,

de dois canais. Unicodado disponfvel na literatura e aquele correspondente a C2= 10.C.: 2 4 6 8 10 12

koo= 0.24118598 0.61930358 1.2971913 2.8654397 10.433999 -13.71298

two-potential kol = - 0.35882732 - 0.65151137 -1.1142349 -2.1315463 -6.9422345 8.3289394

k11 = 0.09666422 0.24820887 0.51989746 1.1484312 4.1818123 -5.495986

koo= 0.24393307 0.63806476 1.3838571 3.3344712 21.770724 -8.0868950

MCF kol = - 0.36116575 - 0.66475802 -1.171373 - 2.4321168 -14.130897 4.7670905

k11 = 0.097788895 0.25585159 0.55503462 1.3377105 8.7360532 -3.2458792

21.7653

Exato43 - 14.1274

8.733

Observa-se que para C2 = 10, os nossos resultados, obtidos com MCF, estao

de matriz- K obtidos com MTPM sao, praticamente, metade dos val ores exatos. Por

outro lado, como e esperado para as interac;6es fracas, e.g., C2 = 2.0, os resultados

Page 36: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Tabela 3 - Secc;oes de choque integral para 0 espalhamento elastico 1~1 (Q11) em

unidades de nao2.

C" 2 4 6 8 10 12

two-potential 0.207429 0.771895 1.39246 1.84652 2.12103 2.27435

MCF 0.211328 0.796658 1.44020 1.90054 2.16792 2.30905

Exato~ 0.211324 0.796646 1.44018 1.90053 2.16791 2.30904

Tabela 4 - Secc;oes de choque integral para 0 espalhamento elastico 1~2 (Q12) em

unidades de nao2.

C" 2 4 6 8 10 12

two-potential 0.385458 0.717192 0.862518 0.857831 0.788288 0.704392

MCF 0.388996 0.726205 0.866754 0.849441 0.767458 0.674327

Exato"'" 0.388993 0.726203 0.866754 0.849444 0.767462 0.674333

Tabela 5 - Secc;oes de choque integral para 0 espalhamento elastico 2~1 (Q21) em

unidades de nao2.

C" 2 4 6 8 10 12

two-potential 1.54183 2.86877 3.45007 3.43132 3.15315 2.81757

MCF 1.55598 2.90482 3.46702 3.39776 3.06983 2.69731

Exato"o 1.55597 2.90481 3.46702 3.39777 3.06985 2.69733

Tabela 6 - Secc;oes de choque integral para 0 espalhamento elastico 2~2 (Q22) em

unidades de nao2.

C" 2 4 6 8 10 12

two-potential 0.245029 0.911814 1.64487 2.18124 2.50551 2.68662

MCF 0.249588 0.940709 1.70029 2.24333 2.55844 2.72445

Exato"'" 0.249583 0.940700 1.70028 2.24333 2.55844 2.72445

Page 37: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Novamente, para toda a faixa de intensidade de potencial estudada, a con-

cordancia entre as sec95es de choque calculadas por MCF estao em excelente

concordancia com os valores exatos. Entretanto, e interessante notar que os resul-

tados de sec95es de choque obtidos por MTPM estao, tambem, em boa concor-

dancia com os resultados de MCF e do calculo exato, mesmo para as interagees

mais fortes.

Atraves deste estudo podemos concluir que 0 MCF e urn metoda muito su-

perior ao MTPM, tanto na exatidao dos resultados quanta no esforyo computacio-

nal. Esta disparidade se torna muito mais evidente nos estudos de espalhamento

multicanais de e+ ~ H, conforme veremos mais adiante.

4.2 ESTUDO DE ESPALHAMENTO e+ ~ H AO NivEL DE

DOIS CANAlS ACOPLADOS (is -2S)

Nesta e na proxima se~o apresentamos os estudos de espalhamento de

positrons de baixa energia por atomos de hidrog€mio,a nlvel de dois canais (15 - 25)

e de cinco canais (15 - 25 - 2Po - 2P1 - 2P_1). 0 principal objetivo de nosso estudo

visa uma investiga~o sistematica sobre os diversos metodos aproximados que sao

aplicaveis aos espalhamentos multicanais, tanto no aspecto da exatidao do metoda

quanta a eficiencia do mesmo. Esperamos que, atraves deste estudo, possamos

selecionar dentre os diversos metodos propostos neste trabalho, aquele que une

ambas as qualidades para nossos estudos futuros de espalhamento de eletrons por

moleculas. 0 processo de espalhamento de positron por atomos de hidrogenio foi

Page 38: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

(desconsiderando tanto a forma~o de positronium quanta os processos de

aniquila~o), pelo fate tambem do positron ser distinto do eletron do alvo e, em

consequ€mcia, nao existir 0 efeito de troca. Por outro lado, alguns resultados de

se~o de choque de espalhamento de e+ ~ H ja foram publicados na lite ratura, 0

onde T~L=O) eo elemento parcial da matriz - T de transi~o, que esta relacionada com

T=- 2iKl-iK

onde, por sua vez, os elementos da matriz - K, Kfj sac dados por:

onde \{It e 0 valor principal da func;ao do espalhamento do canal i, e t/Jf e a func;ao

de onda da partfcula livre, e Ufj = 2Vfj eo potencial de interac;ao reduzido, sendo Vfj

dado por:

Page 39: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

\

32

em que o operador de energia potencial V é dado por:

(4.16)

As expressões dos potenciais de interação para os canais considerados neste

trabalho estão Iistados no Apêndice B.

Na Tabela 7, a seguir, apresentamos as secçóes de choque parcial (L = O) pa-

ra o espalhamento de pósitron por H, ao nível de dois canais, obtidos por diversos

métodos, na energia incidente de 1 Ryd (13.606 eV, k = 1.0).

Tabela 7 - Secçóes de choque total parcial (L= O), em unidades de 1tao2,para o espa­

Ihamento de e+ por hidrogênio, ao nível de dois canais, e para ~ = 1.0.

1S ~181S ~2828 ~2S

PAB

1.4240.19626.4

PABU

0.9970.019113.7

SAB

0.3550.04702.61

SABU

0.3110.02586.21

MODPO

0.6670.092912.0

MODPOU

0.6310.0212.6

MODSO

0.5600.0009549.98

MTPM-OP

0.5540.003511.2

MTPM-OD

0.5300.01526.25

Funcional - C OP

0.5540.003511.2

Funcional - C 00

0.5350.00966.28

MCF

0.5980.009315.01

S.C.:n

0.6080.009315.04

Page 40: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Nossos resultados sac comparados com os de Smith e Burke49 aplicando a

aproximayao de acoplamentos fortes (strong-coupling, SC). Como e esperado, os

resultados obtidos pela primeira aproximayao de Born (PAB) esta muito discrepante

daquelas do SC. Entretanto, aplicando 0 processo de unitarizayao, os resultados da

primeira aproxima~o de Born unitarizada (PABU) apresenta concordancia

sensivelmente melhor com SC. Introduzindo os efeitos de segunda ordem,

observamos que os resultados obtidos pela segunda aproxima~o de Born (SAB) e

tambem a unitarizada (SABU), apresentam resultados aproximadamente da mesma

qualidade; ambos subestimam as secc;5esde choque dos canais de espalhamento

elastico e superestimam as do canal de excitayao, embora apresentem resultados

em melhor concordancia com SC do que PAB. Os metodos da famflia de ondas

distorcidas, tais como ondas distorcidas de primeira ordem (MOOPO), de segunda

ordem (MOOSO) e as unitarizadas (MOOPOU), apresentam resultados em

concordancia sensivelmente melhores com as do SC do que as aproximayOes de

Born, especificamente para os canais de espalhamento elasticos. 0 processo de

unitarizayao aplicado no MOOPO levam a uma melhoria tambem na secy80 de

choque do canal de excitayao (1S-.+2S).

A terceira classe de aproxima~o, os metodos de two-potential multicanais e

de funcional-C multicanais, foram aplicados utilizando tanto as funyOesde ondas pla-

nas (OP) como ondas distorcidas como a base da expansao das funyOesde ondas

do espalhamento. Em termos gerais, os resultados obtidos nao representam uma

melhoria significativa em rela~o as aproximayOesde ondas distorcidas.

Finalmente, as seC¢es de choque obtidas pelo metodo de fra~o contfnuas

(MCF) apresenta uma excelente concordancia com os resultados do SC, mostrando

que e 0 metoda muito superior, em termos de exatidao do metodo, a todos os outros.

Page 41: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Na Tabela 9, apresentamos uma comparac;aodas secy5es de choque parciais

(L=O) calculadas por MCF com os valores obtidos pela SC na faixa de energia

compreendida entre k2 = 1.0 Ryd a k2 = 4.0 Ryd. Todos os resultados do MCF sao

convergidos em 3 itera¢es. Observe-se que a concordancia entre estes dois con-

juntos de secy5es de choque a excelente nas energias em estudo. Esta boa concor-

dancia a esperada, pois como ja foi mostrado na secc;aoII, uma vez que 0 processo

iterativo do MCF tenha convergido, os resultados correspondem de fato aos valores

exatos dentro da precisao numarica do calculo.

Em termos de esforyo computacional, sem dlivida nenhuma que a PAB se

constitue 0 mais economico, a PABU a tambam muito economico, entretanto, os

resultados obtidos por estes matodos deixam muito a desejar. Ainda referindo-se as

aproxima¢es de primeira ordem, os MODPO e MODPOU sac muito superiores em

relayao a PAB e PABU. Em particular, os resultados obtidos por MODPOU

apresentam as mesmas qualidades que MTPM e 0 MCFM, entretanto, computa-

cionalmente a muito eficiente e nao requer muito equipamento. Por outro lado, 0

MODPOU nao requer mudan~s grandes no c6digo computacional do MODPO

normal e a, portanto, urn matodo que merece ser investigado com mais atenyao e

profundidade.

Por outro lado, os matodos de segunda ordem (SAB, SABU, MODSO) e os

MTPM e MFCM requerem esforyos computacionais consideravelmente maiores. A

grande dificuldade computacional envolvida nestes matodos a 0 calculo do termo

(Xllv;G1v;lx:) que foi tratado com detalhes no Apendice A. A parte mais demorada

se refere a integrayao sobre 0 momenta k, onde, para obter a precisao desejada foi

necessario computar os integrandos numa malha suficientemente fina em k, e tam-

barn, estender a integrac;aoata k = 4.0 u.a. ::::::218 eV, para garantir a convergencia.

Page 42: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

Os calculos efetuados com base em ondas planas sac substancialmente mais

rapidos porque a fun<;ao de Green da partfcula livre G;(r,r'), que pode ser expan-

Go+(ii,r; ') = -ik LJ/kr'<h(krJy,:"(Pr)Y"m(f;)1,m

e os termos (¢f !VG ovl ¢i ), podem ser integrados diretamente, utilizando metodos

rmente para este estudo de dois canais, 0 tempo de computa<;ao do MCF, com tres

iterayi5es, e comparavel ao MODPO, entretanto, os resultados sao muito melhores.

Page 43: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

4.3 ESTUDO DE ESPALHAMENTO e+~H AO NivEL DE

CINCO CANAlS ACOPLAD05 (15 - 2S - 2Po- 2P1 - 2P-1)

Agora. reportamos nosso estudo comparativo de resultados de secy6es de

choque do espalhamento de positrons por H, ao nrvel de cinco canais (18 - 28 - 2Po-

2P1 - 2P-1). acoplados. Estes resultados sac mostrados na Tabela 8. para uma

energia incidente k2 = 1.0 Ryd. Os resultados de Burke et al51, obtidos pelo metodo

de "close-coupling" (CC). com tras canais, sac tambem mostrados.

Tabela 8 - SeC96eS de choque parcial (L=O), em unidades de 1r.al, para cinco ca-

nais. e para ki = 1.0.

1S 41S 1S 42S 1S 42P

PAS 1.42 0.196 0.00163

PABU 0.997 0.182 0.0065

SAB 0.290 0.0249 0.0869

SABU 0.283 0.0140 0.00436

MOOPO 0.667 0.093 0.00000894

MOOPOU 0.633 0.0186 0.0054

MODSO 0.464 0.00433 0.00486

MOOSOU 0.468 0.0060 0.0011

MTPM OP 0.510 0.000142 0.00198

MTPM 00 0.697 0.0078 0.0077

Funcional- C OP 0.510 0.000143 0.00194

Funcional - C 00 0.695 0.0131 0.0069

MCF 0.576 0.0032 0.0051

CC51 0.541 0.0070 0.011

Page 44: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

para estudos de dois canais, isto e, as aproximayaes do tipo Born nao apresentam

apresentam concordancia bem melhor que os das aproximayaes de Born. Por outro

lado, os resultados obtidos pela classe de MTPM, MFCM e funcional-C sao,

Nas Tabela 9 e 10, e nas Figuras 4.1 a 4.3, sao mostrados as comparayaes

entre nossos resultados obtidos com 0 MCF, a dois e cinco canais acoplados, com

aqueles do CC, a tras canais acoplados, na faixa de energia entre 0.81 e 4.0 Ryd.

Tabela 9 - Secy5es de choque de onda parcial e+- hidroganio para 2 canais.

L=O S~o de Choque Parcial

k(u.a.) 1S~1S 1S~2S 2S~2S

1.0 a 0.598 0.0093 15.01b 0.608 0.0093 15.04

1.2 a 0.487 0.0192 5.71b 0.487 0.0192 5.71

1.5 a 0.356 0.0209 2.58b 0.356 0.0209 2.58

2.0 a 0.221 0.0125 1.11b 0.221 0.0125 1.11

a: resultados deste trabalhob: resultado da aproxima<;8o de strong-coupling de Smith et al (1960)52.

palhamento 1S~1S e 1S~2S. Para 0 canal de excita<;8o 1S~2P, a discrepancia ebem maior, aproximadamente urn fator de 2. Provavelmente, esta discrepania se de-

ve a diferentes maneiras em que os calculos foram conduzidos. No calculo de Burke

Page 45: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

38

et ai, apenas um canal de espalhamento 2P foi explicitamente incluído, enquanto que

em nossos cálculos, todos os canais de 2P foram levados em conta. Para energias

maiores (que 4.0 Ryd), a contribuição da excitação 1S~2Pn (n = O, 1 e -1), são

praticamente iguais, e nosso resultado de MCF coincide com os resultados de CC.

Tabela 10 - Secçóes de choque de onda parcial e+-hidrogênio para 5 canais.

L=O Secção de Choque Parcial

k(u.a.)

15~ 1515~ 2515~ 2P

0.9

a0.6420.001670.00222b

0.5860.00700.0011.0

a0.5760.00320.00509b

0.5410.00700.0111.2

a0.4690.01360.00472b

0.4530.01240.012'1.5

a0.3480.01780.00263b

0.3430.01710.0052.0

a0.2170.01120.00105b

0.2180.01160.001

a: resultados deste trabalhob: resultado da aproximação de close-coupling de Burke e Smith51 (1962).

0,7

1ij

0,6

'i~

0,5

-8(J)

0,4~T""

.§"~o~

0,3

aS- ~0,2

fi) I.I,I.I.I.I

0,51,01,52,02,53,03,54,04,5

k

Figura 4.1 - Curva a: nosso cálculo pelo MCF e, curva b: pelo CC.

- j-.' .••• < •• ".,...-._.~_ •• -_._.,..~

Page 46: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

iti 0,015E~~r::0 0,010(1)000d.;~~8"r::0

~

0,005

0

0,0000,5

0,014

0,012

itiE 0,010il§ 0,0lJ5

~a.~d.; 0,005

~~0

0,0048

~0,002

0

0,0000,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5

k

Page 47: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

CAPiTULO 5 - CONCLUSOES

Neste trabalho, nosso principal objetivo foi desenvolver urn formalismo teorico

para 0 estudo da interayao de eletrons com moleculas. Propomos os formalismos

teoricos MTPM, Funcional- C e MCF, com 0 intuito de, ao final, escolhermos 0 me-

Ihor metodo, ja que fizemos urn estudo comparativo dos mesmos para 0 espa-

Ihamento multicanal de positron por atomos de hidrogenio, a nivel de 2 e 5 canais,

testando, desta maneira, suas eficiencias computacionais bem como a exatidao dos

resultados obidos. Em todos os testes realizados, 0 MCF se mostrou como 0 mais

eficiente dos tres, apresentando resultados "exatos" no nivel de aproximayao impos-

to nos calculos, tanto que ja esta sendo aplicado no estudo do espalhamento de

eletrons por moleculas, ao nivel da aproximayao estatico-troca. A generalizayao do

MCF, para 0 espalhamento multicanal de eletrons por moleculas, tambem ja esta em

andamento.

Por outro lado, este nosso trabalho revelou que algumas aproxima¢es, de

baixa ordem, especificamente 0 MODPOU, tern apresentado resultados bastante

Page 48: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

razo8veis em comparayao com outros metodos. Devido ao baixo custo computacio-

nal, 0 MODPOU serviria como metoda para primeiros estudos de problemas mais

complexos como, por exemplo, 0 calculo de parametros de correlayao e coerencia de

excitayao molecular por impacto de eletrons.

Page 49: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

APENDICE

G1 = 2:J dE/qn){qn/GllqIR)(qIR/nm

Page 50: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

eTomando En = -2!...+ En e tambem ~n = Xn'Pn' onde 'Pn e a fun9Bo de onda do2

Tfn = (~)2:/1-

1j (\}Ik/l/m/ (;)IVfnl\}lklll~ (;))rr;mj {kf )rrl~ (kJII"".Ijmj

matriz T~2) que:

Page 51: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

T; = 2(.;0)' I ;'.-I,.J,-IfY"" (k,)r".Jf f) f k~~:2('I',flf.f IV•. I'P"~"1X'I'~I,..,IV", 1'1',,,,,, )4~~~ n12"'llfmf

T~2)= 2(.~:YL:i4-'f[L:f k;~k('Pkf/fmt IV/h{'Pk111n;)]('Pkflfmt Iv,,; { 'Pkl~n;)r;i~(i; )r;tnV (if) (A.11)1,,,,, 'In; nlfmf

ou, usando ondas parciais:

Page 52: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

1=J kdkA(k) = .!.[J'"A(k)kdk J'" A(k)kdk] =.!.(I +1 )k'l-e 2 k -k k +k 2 I 2

nOn 0 n

Page 53: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

0,014 0,014

0.012

~

0,012

0,010 0,010

0,003 0,003lI) I \ lI)0 '0~ 0,00;

\ ~ 0,00;lI) .•....•..

0.004 0.004

o,lXI! O,lXI!

o,em o,em

-o,lXI! -o,lXI!0 2 3 4 0 2 3 4

K K

Figura A.1 FiguraA.2

0,04

0,00

I 0,03

-0,01

:<l ca. 0,02

lI) -0,02 ~N e-ll) ~.•.. .•..0,01

-0,03

Page 54: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

O,OXl4

~

O,lXm

\o,am O,lXXl6

\O,lXXl2 O,OXl4

LQ. 0,tm1 ~N

~ 0,lXXl20 ~ 11Q. c.N N~O,CXJX) ~O,CXJX)

-o,tm1 -o,am

-o,lXXl2 -O,OXl4 V-o,am -o,lXXl6

0 2 3 4 0 2 3 4

K K

Figura AS Figura A6

0,001

0,20 {\ a,aD r-0,001

0,15

I\

\-0,002

~ 0,10

Q.N -o,CXl3~

N C.'II NN ca-O,lXl4

N

0,00 \\ .Q,005

~- -0,003

0,00-o,OOT

0 2 3 4 0 2 3 4K K

Figura A7 Figura A8

Page 55: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

B - POTENCIAIS DE INTERACAO

A seguir, Iistamos os potenciais de interar;ao utilizados nos calculos de secr;ao

Page 56: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

REFERENCIAS

01. Czyzk S.L., "Atomic Process with Special Aplication of Gaseous Nebulae", in

Nebulae and Interstellar Matter, Eds. Middlehurst, B.M. and Allen, L.H.

(University of Chicago Press, 1969) p. 403.

02. "Atmospheric Physics and Chemistry", vol. 1 in Applied Atomic Collisions, Eds.

Massey, H.S.W. and Bates, D.R. (Academic Press, New York, 1982).

03. "Gas Lasers", vol. 11\ in Applied Atomic Collisions Physics, Eds. McDaniel, E.W.

(Academis Press, New York, 1982).

04. Capitelli M., Dilonardo M. and Gorse C., Chern. Phys. 56, 29 (1981).

05. Moisan M., Barbeau C., Claude R., Ferreira C.M., Margot J., Paraszczak J., Sa

AB., Sauve G. and Wertheimer M.R., J. Vac. ScL Technol. B 9(1), § (1991).

06. Trajmar S., Register D.F. and Chutjian A, Phys. Rep. 97, 219 (1983).

07. Collins L.A and Schneider B.I., in "Electronic and Atomic Collisions", Edited by

Gilbody H.B., Newell W.R., Read F.H. and Smith AC. (Amsterdam, Elsevier).

08. Lima M.AP., Gibson T.L., Huo W.M. and McKoy V., J. Phys. B: At. Mol. Opt.

Phys. 18, 865 (1985).

09. Lima M.AP., Gibson TL., McKoy V. and Huo W.M., Phys. Rev. A 38,4527

(1988).

10. Sun Q., Winsteade C., McKoy V. and Lima M.AP., J. Chern. Phys. 96(5), 3531

(1991).

11. Pritchard H.P., McKoy V. and Lima M.AP., Phys. Rev.A41(1), 546 (1990).

12. Sun Q., Winstead C. and McKoy, Phys. Rev. A 46,6987 (1992).

13. Rescigno TN. and Schneider B.I., J. Phys. B 21, L691 (1988).

14. Lengsfield B.H. and Rescigno TN., Phys. Rev. A 44,2913 (1991).

15. Rescigno TN. and Schneider B.I., Phys. Rev. A 45(5),2894 (1992).

16. Baluja K.L., Noble C.J. and Tennynson J., J. Phys. B: At. Mol. Opt. Phys. 18,

L851 (1985).

Page 57: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

18. Branchet S.E and Tennynson J. and Morgan L.A, J. Phys. B: At. Mol. Opt.

Phys. 23, 4625 (1990).19 Noble C.J. and Burke P.G., Phys. Rev.Lett. 68, 2001 (1992).

20. Parker S.D., McCurdy C.W., Rescigno TN. and Lensfield IIIB.H., Phys. Rev.

43,3514 (1991).

21. Germano J.S.E and Lima M.AP., Phys. Rev. A 47,3976 (1993)

22. Bettega M.H.F., Policastro AP., Ferreira L.G. and Lima M.AP., Abstracts of

Contributed Papers Vol.1, 260, XVII1ICPEAC,AARHUS, Denmark (1993)

23. Andersen N., Gallager J.W. and Hertel I.V., Phys. Rep. 65, 1 (1988).

24. Batschat K. and Madison D.H., J. Phys. B: At. Mol. Opt. Phys. 20,5839 (1987)

25. FlifietAW. and McKoy V., Phys. Rev. A21, 1863 (1980)

26. Lee M.T. and McKoy v., J. Phys. B: At. Mol. Opt. Phys. 15,3971 (1982).

27. Lee M.T and McKoy V., Phys. Rev. A 28,697 (1983).

28. Lee M.T, Brescansin L.M. and Lima M.AP., J. Phys. B: At. Mol. Opt. Phys.

23,3895 (1990).

29. Lee M.T., Machado L.E, Leal EP., Brescansin L.M., Lima M.AP. and Macha-do F.B.C., J. Phys. B: At. Mol. Opt. Phys. 23, L233 (1990).

30 Lee M.T, Machado L.E, Brescansin L.M., and Menezes G.D., J. Phys. B: At.Mol. Opt. Phys. 24,509 (1991).

31. Lee M.T, Machado L.E, Brescansin L.M., and Menezes G.D., J. Phys. B: At.Mol. Opt. Phys. 24, 509 (1991).

32. Menezes G.D., Brescansin L.M., Lee !\II.T.,Michelin S.L., Machado L.E, and

Csanak G., Phys. Rev. Lett. 71,1701 (1993).

34. Lee M.T., Michelin S.L., Machado L.E, and Brescansin L.M., J. Phys. B: At.Mol. Opt. Phys. 26, L203 (1993).

35. Lee M.T, Michelin S.L., Kroin J., Machado L.E, and Brescansin L.M., J. Phys.B: At. Mol. Opt. Phys. 28, 1859 (1995).

36. Khakoo M.A and Trajmar S., Phys. Rev. A 34, 146 (1986).

37. O.S. Silva Jr, Tese de Doutorado, IFGW- UNICAMP, 1993.

Page 58: Ililliiiil - Biblioteca Digital de Teses e Dissertações da USP · Variacional de Schwinger Multicanal (SMC)8-12,Metodo Variacional de Kohn13-15e matriz_R16-19.Estes metodos apresentam,

51

38. Horacek J. and Sasakawa T., Phys. Rev. A 28.2151 (1983).

39. Takatsuka K. and MeKoy V., Phys. Rev. A 23,2352 (1981).

40. Lee M.T.,Takatsuka K. and MeKoy V., J. Phys. B: AI. MoI. OpI. Phys. 14,4115

(1981).

41. Horacek J. and Sasakawa T., Phys. Rev. A 30, 2274 (1984).

42. Znojil M., Phys. Rev. A 30,2080 (1984).

43. Makowski AJ., Raczynski A and Staszewska G., Chem. Phys. Leti. 114,325

(1985).

44. Meyer H.O., Horacek J. and Cederbaum L.S., Phys. Rev. A 43,3587 (1991).

45. Makowski A.J., Raczynski A and Staszewska G., Phys. Rev. A33, 733 (1986).

46. Makowski AJ., Raczynski A and Staszewska G., J. Phys.B:AI. Mol.OpI. Phys.

19,336 (1986).

47. Huek R.J., Proe. Phys. Soe. (London), 70, 369 (1957).

48. Nesbet R.K., Phys. Rev. 179,60 (1969)

49 Smith K. and Burke P.G., Phys. Rev. 123, 174 (1961).

50. Burke P.G., Sehey H.M. and Smith K., Phys. Rev. 129, 1258 (1963).

51. Burke P.G., and Smith K., Rev. Mod. Phys., 34, ª-, 475 (1962).

52. Smith K., Miller W.F. and Mumford AJ.P., Proe. Phys. Soe. (London) 76, 532

(1960).