Modelamento do transporte eletrônico em dispositivos moleculares

No presente trabalho, simulamos as propriedades de transporte e espectro de absorção do composto orgânico Vermelho de Etila. Este é o primeiro estudo teórico de um indicador específico de pH utilizado como nanodispositivo, com base na teoria quântica e no modelo de transporte não-difuso. A distribuição de carga ao longo da molécula é determinada através da técnica, Ab initio, como uma função de um campo elétrico externo. Baseado em um modelo de multiníveis ressonantes também calculamos a corrente como função da tensão de polarização. O acúmulo de carga e a corrente apresentam comportamento semelhante, como a condução do tipo ressonante e curvas carga-tensão e corrente-tensão assimétricas. Os principais resultados sugerem que o sistema presente poderia funcionar como um transistor molecular bi-direcional. Estendemos esta metodologia de análise para outro dispositivo molecular, mas composto de três terminais. Para este sistema, nossa descoberta principal é a resistência diferencial negativa (RDN) na carga Q como uma função do campo elétrico externo. Para explicar este efeito RDN, aplicamos um modelo capacitivo fenomenológico, também baseado em um sistema de multiníveis localizados (que podem ser os LUMOs – Lowest Unoccupied Molecular Orbital – Orbitais moleculares desocupados mais baixos). A capacitância descreve, por efeito de carregamento, a causa do bloqueio de Coulomb (BC) no transporte. Mostramos que o efeito BC dá origem a uma RDN para um conjunto adequado de parâmetros fenomenológicos como: taxa de tunelamento e energia de carregamento. O perfil da RDN obtida nas duas metodologias, ab initio e fenomenológica, estão em comum acordo.

Investigação das características I x V e C x V de NCPS puro, com nitrogênio substitucional carregado (-1 e +1) e com grupos doador (NO2)-aceitador (NH2) através de métodos derivados de hartree-fock

Neste trabalho, fizemos uma investigação sobre o estudo teórico das características I x V e C x V de Nanotubo Carbono de Parede Simples (NCPS) puro, com Nitrogênio substitucional carregado com cargas -1 (caracterizando um indicativo de dopagem tipo n) e +1 (caracterizando um indicativo de dopagem tipo p) e na presença de grupos doador (NO₂)-aceitador (NH₂), através da simulação computacional do estado fundamental de NCPS, bem como de sua estrutura eletrônica e propriedades ópticas, utilizando parametrizações semi-empíricas AM1 (Austin Mudel 1) e ZINDO/S-ClS (Zerner´s lntermediate Neglect of Differential Orbital/Spectroscopic - Cunfiguration lnteraction Single) derivadas da Teoria de Hartree-Fock baseada em técnicas de química quântica. Por meio deste modelo teórico analisamos as propriedades ópticas e eletrônicas, de maior interesse para esses materiais, a fim de se entender a melhor forma de interação desses materiais na fabricação de dispositivos eletrônicos, tais como TECs (Transistores de Efeito de Campo) ou em aplicações em optoeletrônica tais como DEL (Dispositivo Emissor de Luz). Observamos que NCPS com Nitrogênio substitucional apresentam defeitos conformacionais do tipo polarônico. Fizemos as curvas dos espectros UV-visível de Absorção para NCPS armchair e zigzag puro, com Nitrogênio substitucional carregado com cargas (-1 e +1) e na presença de grupos doador (NO₂)-aceitador (NH₂), quando perturbados por intensidades diferentes de campo elétrico. Verificamos que em NCPS zigzag ao aumentarmos a intensidade do campo elétrico, suas curvas sofrem grandes perturbações. Obtivemos as curvas p x E, I x V e C x V para esses NCPS, concluímos que NCPS armchair possui comportamento resistor, pois suas curvas são lineares e zigzag possui comportamento semelhante ao dos dispositivos eletrônicos importantes para o avanço tecnológico. Assim, nossos resultados estão de bom acordo com os resultados experimentais e teóricos de NCPS puro e com Nitrogênio encontrados na literatura.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Ensaios sobre a contribuição do ar nos resultados da modelagem do modo transverso-elétrico do método magnetotelúrico

Apresentamos a solução numérica do modelo da falha infinita para o modo TE utilizando o método dos elementos finitos para investigar o efeito do ar na solução TE. A comparação de nossa solução para o campo elétrico normalizado com a solução de Weaver evidenciou uma discrepância esperada devido ao efeito do ar, que foi intencionalmente negligenciado por Weaver, com o intento de facilitar a solução do problema. O problema analítico, levando em conta a presença do ar, é substancialmente mais difícil. Uma solução do problema real, parcialmente analítica e parcialmente numérica, foi apresentada por Sampaio. Devido a uma discrepância inesperada entre a nossa solução por elementos finitos e a solução exata de Sampaio, recalculamos a sua solução. Iterando a parte numérica da solução de Sampaio até ordens mais altas pudemos verificar que esta apresenta comportamento divergente. Para examinar o efeito do ar na resistividade aparente, comparamos três perfis de resistividade aparente obtidos a partir de nossa solução por elementos finitos com os mesmos perfis obtidos pela solução de Weaver. A discrepância encontrada entre estas soluções foi muito pequena, situando-se na faixa de erro instrumental de campo.