Desenvolvimento de um conjunto de modelagem de alta resolução para os métodos eletromagnéticos indutivos

A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre.

Modelagens direta e inversa de dados mCSEM 3D

Marine Controlled Source Electromagnetic - mCSEM é um método geofísico eletromagnético que nos últimos dez anos vem sendo usado na prospecção de hidrocarbonetos com bastante êxito. Este método consiste em um dipolo elétrico horizontal (DEH) localizado um pouco acima do assoalho marinho, operando em baixa frequência (0,1-1,0 Hz) e receptores regularmente distribuídos no fundo do mar que captam os campos eletromagnéticos provenientes da difusão de energia gerada pelo dipolo transmissor. Neste trabalho vamos apresentar o problema direto do método mCSEM 3D, propondo soluções numéricas, através do método dos elementos finitos tridimensionais, para modelos geoelétricos mCSEM 3D. Para fins de análise de coerência, os resultados obtidos são comparados com soluções disponíveis na literatura. Em seguida, apresentaremos a inversão de um de seus modelos segundo uma proposta de metodologia de inversão juntamente com a proposta de solução direta para o mCSEM 3D, acima mencionada, realizando assim a inversão de um modelo geoelétrico do mCSEM 3D para duas frequências.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Equipamento eletromagnético para prospecção geofísica e modelamento reduzido

Nas últimas décadas, a prospecção por métodos eletromagnéticos vem-se constituindo numa técnica eficiente para prospecção mineral. O objetivo deste trabalho foi desenvolver um equipamento para prospecção eletromagnética quantitativa de corpos condutores, através do método dipolo-dipolo, podendo ainda ser usado em modelos reduzidos. Eletricamente, o sistema mede grandezas relacionadas ao acoplamento indutivo entre duas bobinas: transmissora e receptora. Elas são dispostas na superfície da terra, afastadas entre si, e a terra, desse modo, constitui o núcleo acoplador. Quando existem corpos condutores nas proximidades, estes são denunciados por alterações no comportamento do sinal induzido na bobina receptora. O equipamento compreende dois conjuntos: o transmissor e o receptor, além de acessórios. O transmissor gera um campo eletromagnético nas freqüências de 520 e 3.090 Hz, e um sinal de referência para o receptor, o qual é enviado através de um cabo. O receptor, inicialmente, separa o sinal induzido pelos campos secundários gerados por condutores, do campo normalmente recebido, quando a condutividade da subsuperfície é relativamente uniforme (campo primário). Em seguida, decompõe esse sinal em duas componentes ortogonais, uma em fase, e outra em quadratura com o campo primário. Através de duas escalas de precisão, as amplitudes dessas componentes são mostradas como percentagens do campo primário, com precisão de 1%. A sensibilidade do receptor é de 0,5 μV. O circuito eletrônico foi rigorosamente testado com preciso instrumental de laboratório. Em seguida, testou-se sua aplicação no Laboratório de Modelo Reduzido Eletromagnético do NCGG, refazendo-se experiências clássicas, encontradas na literatura especializada. No campo, foi experimentado próximo da cidade de Araci no Estado da Bahia, em áreas prospectadas pela "Rio Doce Geologia e Mineração S/A-DOCEGEO". Em ambos os casos, verificaram-se bons resultados.

Filtros otimizados para transformadas seno, co-seno e de Hankel j0, j1 e j2

Apresentamos dois algoritmos automáticos, os quais se utilizam do método dos mínimos quadrados de Wiener-Hopf, para o cálculo de filtros lineares digitais para as transformadas seno, co-seno e de Hankel J₀, J₁ e J₂. O primeiro, que otimiza os parâmetros: incremento das abscissas, abscissa inicial e o fator de deslocamento utilizados para os cálculos dos coeficientes dos filtros lineares digitais que são aferidos através de transformadas co-seno, seno e o segundo, que otimiza os parâmetros: incremento das abscissas e abscissa inicial utilizados para os cálculos dos coeficientes dos filtros lineares digitais que são aferidos através de transformadas de Hankel J₀, J₁ e J₂. Esses algoritmos levaram às propostas de novos filtros lineares digitais de 19, 30 e 40 pontos para as transformadas co-seno e seno e de novos filtros otimizados de 37 , 27 e 19 pontos para as transformadas J₀, J₁ e J₂, respectivamente. O desempenho dos novos filtros em relação aos filtros existentes na literatura geofísica é avaliado usando-se um modelo geofísico constituído por dois semi-espaços. Como fonte usou-se uma linha infinita de corrente entre os semi-espaços originando, desta forma, transformadas co-seno e seno. Verificou-se melhores desempenhos na maioria das simulações usando o novo filtro co-seno de 19 pontos em relação às simulações usando o filtro co-seno de 19 pontos existente na literatura. Verificou-se também a equivalência de desempenhos nas simulações usando o novo filtro seno de 19 pontos em relação às simulações usando o filtro seno de 20 pontos existente na literatura. Adicionalmente usou-se também como fonte um dipolo magnético vertical entre os semi-espaços originando desta forma, transformadas J₀ e J₁, verificando-se melhores desempenhos na maioria das simulações usando o novo filtro J₁ de 27 pontos em relação ao filtro J₁ de 47 pontos existente na literatura. Verificou-se também a equivalência de desempenhos na maioria das simulações usando o novo filtro J₀ de 37 pontos em relação ao filtro J₀ de 61 pontos existente na literatura. Usou-se também como fonte um dipolo magnético horizontal entre os semi-espaços, verificando-se um desempenho análogo ao que foi descrito anteriormente dos novos filtros de 37 e 27 pontos para as respectivas transformadas J₀ e J₁ em relação aos filtros de 61 e 47 pontos existentes na literatura, destas respectivas transformadas. Finalmente verificou-se a equivalência de desempenhos entre os novos filtros J₀ de 37 pontos e J₁ de 27 pontos em relação aos filtros de 61 e 47 pontos existentes na literatura destas transformadas, respectivamente, quando aplicados em modelos de sondagens elétricas verticais (Wenner e Schlumberger). A maioria dos nossos filtros contêm poucos coeficientes quando comparados àqueles geralmente usados na geofísica. Este aspecto é muito importante porque transformadas utilizando filtros lineares digitais são usadas maciçamente em problemas numéricos geofísicos.

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Investigação por eletroresistividade de ambientes cársticos no município de São João de Pirabas - PA

A presente tese procura investigar ambientes cársticos, usando os métodos eletroresistivos de sondagem elétrica vertical e caminhamento dipolo-dipolo. A área piloto para essa investigação está localizada na porção nordeste do Estado do Pará, no Município de São João de Pirabas. Nessa localidade, durante a perfuração de um poço para abastecimento de água para a cidade, a equipe de prospecção da Fundação Nacional de Saúde encontrou problemas devido a presença de cavidades nos calcários da Formação Pirabas. Para analisar esse problema foram realizadas 15 sondagens elétricas verticais, usando o arranjo Schlumberger, em diversos pontos da cidade. A abertura máxima dos eletrodos de corrente, para a realização das sondagens, foi em torno de 600 metros. As curvas de resistividade aparente das sondagens foram, inicialmente, interpretadas usando programas computacionais que calculam as resistividade do meio investigado com modelos unidimensionais, sem variações laterais de resistividade. Porém, as sondagens realizadas próximas ao poço, onde ocorrem cavernas, foram interpretadas usando modelos bidimensionais para representar as variações laterais de resistividade devido às cavernas. O algoritmo usado para essa interpretação se baseia na técnica dos elementos finitos. Os resultados das interpretações das sondagens elétricas, para os vários modelos de cavernas estudados, mostraram que as anomalias de resistividade devido a presença de cavernas são muito pequenas. Devido a este fato, o método de sondagem elétrica vertical não apresentou boa resolução para identificar com precisão a presença deste tipo de caverna em subsuperfície. Por outro lado, os resultados numéricos obtidos pelo caminhamento dipolo-dipolo mostraram que esse método é mais eficiente e menos ambíguo, para estudar cavernas, que o método de sondagem elétrica vertical. Contudo o método de caminhamento dipolo-dipolo é muito mais dispendioso do que o de sondagem elétrica verticais, tanto do ponto de vista de trabalho de campo como no processamento e interpretação de dados.

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

Modelagem de eletrorresistividade 2-D a partir do potencial elétrico secundário

Um dos métodos clássicos da geofísica de exploração é o Método de Eletrorresistividade, estabelecido há um século pelos irmãos Schlumberger e desde então amplamente empregado em prospecção mineral, estudos ambientais e hidrogeologia e em pesquisa de fontes geotermais. Conceitualmente o método consiste de injeção de corrente elétrica na subsuperfície e de medida de diferença de potencial elétrico, resultante da interação da corrente com o meio. As localizações dessas fontes e receptores são determinadas pelo arranjo escolhido para o levantamento. Após o processamento, obtém-se pseudo-seções de resistividade aparente que indicam a distribuição de condutividade em subsuperfície. Devido à simplicidade dos fundamentos físicos de sua formulação, o método apresenta fácil implementação computacional quando comparado aos métodos eletromagnéticos de fonte controlada. Na literatura há inúmeros trabalhos de modelagem computacional, onde se calcula a resposta para problemas 2-D e 3-D. Nestes trabalhos, as pseudo-seções são obtidas a partir do cálculo do potencial elétrico total. Neste trabalho, apresentaremos a resposta da modelagem de eletroresistividade 2-D com o arranjo dipolo-dipolo, obtida a partir do potencial elétrico secundário. A solução é calculada através do método de elementos finitos usando malhas não estruturadas. Para efeito de validação, os resultados são comparados com a resposta 2-D obtida a partir dos potencias totais.

Modelagem e filtragem de distorções estáticas em levantamentos de dipolos contínuos

As distorções estáticas são um grave problema que afeta o método magnetotelúrico e muitas tentativas têm sido feitas para eliminar ou minimizar os seus efeitos. Este trabalho trata de uma técnica nova que trata deste problema, o método EMAP, que é uma adaptação do método magnetotelúrico na qual as medidas de campo elétrico são feitas em uma linha contínua de dipolos conectados entre si, e os dados coletados desta maneira são tratados com um filtro espacial passa-baixa dependente da frequência. Este trabalho é composto de duas partes principais, a simulação numérica dos dados contaminados com as distorções estáticas e a filtragem destes dados com o filtro espacial passa-baixa. Na primeira parte, aplicamos o método dos elementos finitos para simular a resposta dos dipolos elétricos, os quais fazem as medidas do campo elétrico. Na segunda parte aplicamos a janela de Hanning como filtro passa-baixa apropriado para tratar os dados.

Avaliação quantitativa da dinâmica espaço-temporal da precipitação na região hidrográfica Tocantins-Araguaia

A análise do comportamento da precipitação em uma bacia hidrográfica é fundamental para a engenharia e gerenciamento dos recursos hídricos. A Região Hidrográfica Tocantins-Araguaia (RHTA) pela sua ocupação recente e potencialidades econômicas, ganha destaque no cenário nacional. Este trabalho avalia quantitativamente a dinâmica espaço-temporal da precipitação anual nesta região durante um período de 30 anos de dados. A dinâmica da precipitação pode ser analisada pelo cálculo da precipitação média em uma dada área, compondo assim mapas de isoietas de precipitação anual. No entanto, a confecção destes mapas requer um método de interpolação que melhor represente as características pluviométricas em locais não amostrados para posterior análise quantitativa do comportamento da precipitação. Para tanto, foram realizados análises exploratórias descritivas amostral e espacial como requisito de estacionaridade do método de interpolação geoestatístico, ajuste e validação do modelo teórico que se adéque ao variograma de precipitação anual. Após a confecção do mapa de isoietas pelo método de Krigagem Ordinária (sem tendência) e Krigagem Universal (com tendência) foi realizado o cálculo do volume precipitado na região hidrográfica pelo método dos contornos. A dinâmica espacial da precipitação foi realizada com base na análise da estatística descritiva, mapa de isoietas, mapa hipsométrico, Índice de Irregularidade Meteorológica (IMM) e Coeficiente de Variação. A dinâmica temporal foi analisada pela distribuição dos totais anuais de precipitação volumétrica para cada sub-bacia da RHTA, pelo Índice de Anomalia Padronizada, na variação interanual de precipitação e teste de tendência e magnitude representados respectivamente pelo Teste de Mann Kendall e Sen’s. Os resultados correlacionados com as anomalias meteorológicas do Oceano Atlântico (Dipolo) e Pacífico (ENOS) indicam o comportamento da precipitação bastante heterogêneo e com grande variabilidade temporal principalmente na sub-bacia Tocantins-Alto (TOA) (14%). Diminuição da amplitude pluviométrica, em anos de anomalia meteorológica intensa ocasionando um incremento de precipitação ao sul das sub-bacias TOA e ARA e diminuição da precipitação na sub-bacia TOB, em eventos de El Niño. Não se pode comprovar pelo teste de Mann Kendall que há uma tendência estatisticamente significativa no volume precipitado na RHTA, mas o estimador Sen’s dá indícios de queda na precipitação na sub-bacia TOA (-1,24 Km^³/ano) e Araguaia (ARA) (-1,13 Km^³/ano) e aumento da precipitação na sub-bacia do Tocantins Baixo (TOB) (0,53 Km^³/ano) e para a RHTA (-1,5 Km^³/ano). Assim a variabilidade espacial e temporal nas sub-bacias está intimamente relacionada aos eventos de anomalia meteorológica, na qual, a sua ação ocorre de maneira irregular ao longo da área de estudo e pode influenciar as diversas atividades sócio-econômicas na RHTA de acordo com sua magnitude e área de ocorrência.

Índices de extremos climáticos e tendências da precipitação anual e sazonal na Amazônia Oriental

Baseado no conjunto de dados diários de precipitação e temperatura do ar de oito estações meteorológicas sobre o Pará, pertencentes ao INMET, disponíveis no período 1961-2007, foram calculados índices de extremos climáticos através da metodologia estatística do software RClimdex. Utilizando-se ainda um conjunto de dados mensais de precipitação de 134 estações pluviométricas da Amazônia oriental, provenientes da ANA e INMET, foram analisadas as tendências anual e sazonal num período de 25 a 37 anos (1970-2009), através da aplicação da estatística não-paramétrica (teste de Mann-Kendall e método de Sen). E por fim, utilizaram-se dados de desmatamento do PRODES/INPE no período 2000-2007, para analisar a possível influência do desmatamento nas tendências de precipitação anual e sazonal do Pará. Considerando os índices de extremos climáticos referentes aos dados diários de precipitação, constatou-se que a variabilidade interanual dos índices possui relações com os episódios ENOS e as fases do Dipolo do Atlântico. Para os extremos climáticos de temperatura evidenciaram-se reduções sistemáticas dos dias frios (TX10p) e aumentos sistemáticos da mínima da máxima temperatura (TXN), da percentagem de noites quentes (TX90p) e da temperatura máxima da máxima (TXX). Os resultados das tendências da precipitação anual usando os testes não-paramétricos evidenciaram que do total de 134 estações, em torno de 51% apresentam tendências positivas e 41% tendências negativas sobre a Amazônia oriental. Contudo, essas tendências anuais “mascararam” de certa forma as tendências sazonais, cujos resultados mostraram-se mais diversificados, explicando melhor os aspectos da variabilidade climática regional. Para o trimestre MAM, observaram-se 99 estações (74%) com tendências positivas, indicando que a quantidade de precipitação durante o período chuvoso principal vem aumentando sistematicamente durante a última década. Em DJF, notaram-se 45 estações (34%) com tendências positivas e 78 estações (58%) com tendências negativas. Por outro lado, nos trimestres que englobam o período seco ou menos chuvoso verificaram-se tendências sazonais predominantemente negativas, com 84 estações (62%) em JJA e 89 estações (66%) em SON. Portanto, a evidência observacional de que a precipitação do período seco encontra-se em diminuição gradativa nas últimas décadas corrobora com a hipótese de que o desmatamento associa-se com a redução da precipitação em escala regional, porém sugere-se que isso ocorra em escala de tempo sazonal.

Aspectos regionais da variabilidade de precipitação no estado do Pará: estudo observacional e modelagem climática em alta resolução

O objetivo do presente trabalho foi agregar diferentes redes de estações meteorológicas de superfície para a criação de um novo banco de dados integrado, a partir do qual foi gerada uma climatologia recente (1978-2007) para a precipitação do estado do Pará em alta resolução espacial – 30 km, permitindo melhor identificar a variabilidade climática regional, sobretudo influenciada pelos aspectos da fisiografia e em função de mecanismos climáticos de grande escala dos oceanos Pacífico e Atlântico. Buscou-se, ainda, estabelecer uma configuração otimizada do modelo climático RegCM3 utilizando duas diferentes parametrizações de cumulus: RegCM3/Grell e RegCM3/MIT. Foram realizadas 26 simulações (1982/83 a 2007/08) durante a estação chuvosa na Amazônia oriental (dezembro a maio) para cada esquema de parametrização convectiva, utilizando 30 km de resolução espacial. O modelo mostrou-se capaz de capturar os sinais de anomalia na presença de forçantes climáticas extremas, como o El Niño-Oscilação Sul e o dipolo do Atlântico. O RegCM3/MIT obteve ótimo desempenho na região de Altamira/PA e performance razoável nos setores Nordeste (região de Belém), Leste ( região de Marabá), Sudeste (região de Conceição do Araguaia), e Noroeste (região de Tiriós). O RegCM3/Grell destacou-se nas regiões Nordeste, Leste, Sudeste e Noroeste, com desempenho razoável. O setor Norte (região de Macapá) foi o mais problemático, com pouca ou nenhuma sensibilidade apresentada pelo modelo. Embora o RegCM3 tenha obtido resultados razoáveis na maior parte do domínio, foram detectados erros sistemáticos nas simulações, com viés seco para o RegCM3/Grell e viés úmido para o RegCM3/MIT na porção Sul e viés seco na porção Norte. Estas características denotam a necessidade de ajustes às condições regionais dos esquemas de convecção.

Analysis of modified bowtie nanoantennas in the excitation and emission regimes

In this work, we analyze modified bowtie nanoantennas with polynomial sides in the excitation and emission regimes. In the excitation regime, the antennas are illuminated by an incident plane wave, and in the emission regime, the excitation is fulfilled by infinitesimal electric dipole positioned in the gap of the nanoantennas. Several antennas with different sizes and polynomial order were numerically analyzed by method of moments. The results show that these novel antennas possess a controllable resonance by the polynomial order and good characteristics of near field enhancement and confinement for applications in enhancement of spontaneous emission of a single molecule.

Incorporação de tiofosfatos em vidros calcogenetos: estrutura e fenômenos fotoinduzidos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)