Configurações magnéticas colineares e não-colineares em currais de Fe, Cr e Mn sobre a Pt (111)

Neste trabalho investigamos as propriedades magnéticas de currais de Fe, Cr e Mn adsorvidos sobre a superfície de Pt(111) utilizando o método RS-LMTO-ASA (Real Space Linear Muffin Tin Orbital - Atomic Sphere Approximation), o qual é um método de primeiros princípios baseado na Teoria do Funcional da Densidade (DFT-Density Functional Theory), que permite o cálculo de estruturas magnéticas não-colineares. Obtivemos que os átomos de Fe apresentam momentos magnéticos elevados, da ordem de 3.5µB /átomo, e têm uma interação de troca entre primeiros vizinhos forte e ferro-magnética. Isto leva a um arranjo magnético colinear no curral. Para os currais de Mn e Cr encontramos que estes possuem elevado momento magnético, da ordem de 4.51µB /átomo e 4.15µB /átomo, respectivamente, e interações de troca entre primeiros vizinhos antiferro-magnéticas. Isto conduz a arranjos magnéticos colineares em currais simples, assim como interessantes ordenamentos não-colineares, tais como estruturas tipo vértice (skyrmions), para os currais com uma geometria particular onde o antiferromagnetismo se apresenta frustado.

Propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas

Neste trabalho, utilizamos o método de primeiros princípios RS-LMTO-ASA (Real Space – Linear Muffin-Tin Orbital - Atomic Sphere Approximation) baseado na Teoria do Funcional da Densidade (DFT - Density Functional Theory) e implementado para o cálculo de estruturas magnéticas não-colineares, para investigar as propriedades magnéticas de nanoestruturas adsorvidas em superfícies metálicas. Consideramos aglomerados com diferentes geometrias e tamanhos como adátomos, dímeros, trímeros, nanofios e nanoestruturas de geometria triangular de Fe, Fe-Co e Fe-Pt adsorvidos sobre a superfície de Pt(111) e tratamos também nanoestruturas de Mn sobre a superfície de Ag(111). Mostramos que os nanofios de Fe-Co sobre a superfície de Pt(111) apresentam um ordenamento ferromagnético. Devido à redução do número de coordenação presente na superfície, os momentos de spin e orbital nos sítios de Fe e Co mostram-se elevados comparados com os respectivos valores dos momentos destes metais como bulk. Analisamos também como estes momentos variam em função da concentração destes elementos nos nanofios. Para os sistemas compostos por nanofios Fe-Pt adsorvidos em Pt(111), mostramos que é possível sintonizar as interações de troca entre os adátomos magnéticos Fe através da introdução de um diferente número de átomos Pt para ligá-los. Por exemplo, a interação de troca entre os adátomos de Fe pode ser consideravelmente aumentada pela introdução de cadeias de Pt que os conectem e tanto configurações ferromagnéticas, antiferromagnéticas ou não-colineares entre os adátomos de Fe podem ser estabilizadas, dependendo da espessura do espaçador Pt. Para os aglomerados Mn sobre a Ag(111) mostramos que a interação de troca entre os sítios de Mn depende não somente da distância entre os átomos, mas também do número de coordenação de cada sítio. Desta forma, verificamos um magnetismo não-colinear nestas nanoestruturas causado tanto por frustração geométrica, quanto pela competição de interações de curto e longo alcance. Nossos resultados estão em boa concordância com os resultados experimentais da literatura e com os resultados teóricos obtidos por outros métodos, quando existentes.

Propriedades magnéticas de nanoestruturas de metais de transição 3d adsorvidas na superfície de Pt(111)

Neste trabalho, utilizamos o método de primeiros princípios, RS-LMTO-ASA (“Real Space - Linear Muffin-Tin Orbital - Atomic Sphere Approximation”), baseado na Teoria do Funcional da Densidade (DFT) e implementado para o cálculo de estruturas magnéticas não-colineares, para investigar as propriedades magnéticas de nanoestruturas de metais de transição 3d (Cr, Mn, Fe, Co e Ni) adsorvidas na superfície de Pt(111). Diferentes geometrias como adátomos, dímeros, trímeros, fios lineares e zig-zag foram consideradas e, o tamanho dos aglomerados foi variado de 2 a 7 átomos. Mostramos que os aglomerados de Fe, Co e Ni sobre a superfície de Pt(111), para todas as geometrias simuladas, apresentam um ordenamento ferromagnético. Devido à redução do número de coordenação presente na superfície, os momentos de spin e orbital nos sítios de Fe, Co e Ni, para as diferentes geometrias, mostram-se elevados comparados com os respectivos valores dos momentos destes metais como bulk. Para os glomerados de Cr e Mn mostramos que a interação de troca antiferromagnética entre primeiros vizinhos leva a um ordenamento antiferromagnético colinear no caso de geometrias lineares. No entanto, se o antiferromagnetismo é frustrado por restrição geométrica imposta aos aglomerados pela superfície triangular do substrato, obtém-se um comportamento magnético não-colinear para aglomerados de Cr e Mn sobre a Pt(111). Nossos resultados estão em boa concordância com os resultados experimentais da literatura e com os resultados teóricos obtidos por outros métodos, quando existentes.

Método de otimização não-linear para impor PSDs arbitrárias em modems DSL: análise estatística e prática

Este trabalho tem como objetivo apresentar um método para solucionar o problema de mapeamento entre as soluções teóricas de gerenciamento dinâmico de espectro (DSM) e os parâmetros de controle das densidades espectrais de potência (PSDs) de transmissão dos modems comerciais. O método utiliza algoritmos genéticos (AG) codificado em inteiros para solucionar o problema. O AG é responsável por achar os melhores parâmetros para representar uma PSD arbitrária, considerando as restrições impostas pelos equipamentos e padrões atuais DSL. O trabalho apresenta ainda um estudo comparativo do método proposto com um concorrente, além de estudo estatístico do método proposto, considerando média, desvio padrão e intervalo de confiança. Adicionalmente, são apresentados dois setups para uso em laboratório, sendo um para medição de PSDs e o outro para funções de transferência, os quais podem ser reaproveitados em outros trabalhos.

Wave equation depth migration using complex Padé approximation

Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domínio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido às ondas evanescentes e produz imagens com menos ruídos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos íngremes. Os resultados foram comparados àqueles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD).

Evaluation via multivariate techniques of scale factor variability in the rietveld method applied to quantitative phase analysis with X ray powder diffraction

ABSTRACT: The present work uses multivariate statistical analysis as a form of establishing the main sources of error in the Quantitative Phase Analysis (QPA) using the Rietveld method. The quantitative determination of crystalline phases using x ray powder diffraction is a complex measurement process whose results are influenced by several factors. Ternary mixtures of Al2O3, MgO and NiO were prepared under controlled conditions and the diffractions were obtained using the Bragg-Brentano geometric arrangement. It was possible to establish four sources of critical variations: the experimental absorption and the scale factor of NiO, which is the phase with the greatest linear absorption coefficient of the ternary mixture; the instrumental characteristics represented by mechanical errors of the goniometer and sample displacement; the other two phases (Al2O3 and MgO); and the temperature and relative humidity of the air in the laboratory. The error sources excessively impair the QPA with the Rietveld method. Therefore it becomes necessary to control them during the measurement procedure.

Modos quasinormais e pólos de Regge para os buracos acústicos canônicos

Usando o formalismo relativístico no estudo da propagação de perturbações lineares em fluidos ideais, obtêm-se fortes analogias com os resultados encontrados na Teoria da Relatividade Geral. Neste contexto, de acordo com Unruh [W. Unruh, Phys. Rev. Letters 46, 1351 (1981)], é possível simular um espaço-tempo dotado de uma métrica efetiva em um fluído ideal barotrópico, irrotacional e perturbado por ondas acústicas. Esse espaço-tempo efetivo é chamado de espaço-tempo acústico e satisfaz as propriedades geométricas e cinemáticas de um espaço-tempo curvo. Neste trabalho estudamos os modos quasinormais (QNs) e os pólos de Regge (PRs) para um espaço-tempo acústico conhecido como buraco acústico canônico (BAC). No nosso estudo, usamos o método de expansão assintótica proposto por Dolan e Ottewill [S. R. Dolan e A. C. Ottewill, Class. Quantum Gravity 26, 225003 (2009)] para calcularmos, em termos arbitrários do número de overtone n, as frequências QNs e os momentos angulares para os PRs, bem como suas respectivas funções de onda. As frequências e as funções de onda dos modos QNs são expandidas em termos de potências inversas de L = l + 1/2 , onde l é o momento angular, enquanto que os momentos angulares e funções de onda dos PRs são expandidos em termos do inverso das frequências de oscilação do buraco acústico canônico. Comparamos os nossos resultados com os já existentes na literatura, que usam a aproximação de Wentzel-Kramers-Brillouin (WKB) como método de determinação dos modos QNs e dos PRs, e obtemos uma excelente concordância dentro do limite da aproximação eikonal (l ≥ 2 e l > n).

Propriedades magnéticas de nanoestruturas de metais de transição 3d em superfícies de Pd

Motivados por estudos experimentais acerca de monocamadas de metais de transição 3d sobre superfícies de Pd, nesta dissertação investigamos o complexo magnetismo de nanoestruturas, embebidas ou adsorvidas, em superfícies metálicas através de cálculos de primeiros princípios. Utilizamos o método RS-LMTO-ASA (Real Space - Linear MuffinTin Orbital - Atomic Sphere Approximation), o qual é baseado na teoria do funcional da densidade (DFT - Density Functional Theory) e implementado para o cálculo de estruturas magnéticas não colineares. Com este propósito, investigamos nanoestruturas embebidas e ligas (2 x 2) de metais 3d (Cr, Mn, Fe, Co e Ni) na superfície Pd (110), além de nanoestruturas de Cr adsorvidas sobre a superfície de Pd (111). Primeiro, para as nanoestruturas embebidas na superfície Pd (110), analisamos a variação do momento magnético de spin orbital com relação ao número de vizinhos e de valência dos metais 3d. Também mostramos que estas estruturas têm ordenamento magnético colinear, exceto as de Cr e Mn, que apresentam magnetismo não colinear associado à frustração geométrica. Para o caso de nanofios de Cr adsorvidos sobre a superfície de Pd (111), verificamos uma configuração colinear antiferromagnética para cadeias com até 9 átomos. Para o nanofio com 10 átomos obtivemos uma configuração tipo antiferromagnética inclinada (canted). No caso de nanoestruturas de Cr bidimensionais, verificamos complexas configurações magnéticas não colineares com diferentes quiralidades.

Inversão de momentos de fonte em métodos potenciais

A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino.

Análise de um esquema linear para estimativa de anisotropia local a partir de dados de onda P em experimentos de VSP multiazimutal

Neste trabalho é apresentada uma análise do esquema de inversão linear para a estimativa de anisotropia na vizinhança de um receptor situado em um poço a partir de da componente vertical do vetor de vagarosidade e do vetor de polarização de ondas P medidops em experimentos de VSP walkaway multiazimutal. Independente do meio acima do geofone (homogêneo ou heterogêneo) e da forma do poço (pode ser direcional ou curvado, vertical e inclinado), a inversão é feita a partir de uma aproximação de primeira ordem em torno de um meio isotrópico de referência. O esquma da inversão é analisado considerando fatores como: o nível de ruído nos dados, o tipo de onda P, o grau de anisotropia do meio, a escolha dos parâmetros no meio isotrópico de referência e grau de heterogeneidade do meio. Os resultados são apresentados.

Detecção do contato entre camadas através do perfil de indução

Dois dos principais objetivos da interpretação petrofísica de perfis são a determinação dos limites entre as camadas geológicas e o contato entre fluidos. Para isto, o perfil de indução possui algumas importantes propriedades: É sensível ao tipo de fluido e a distribuição do mesmo no espaço poroso; e o seu registro pode ser modelado com precisão satisfatória como sendo uma convolução entre a condutividade da formação e a função resposta da ferramenta. A primeira propriedade assegura uma boa caracterização dos reservatórios e, ao mesmo tempo, evidencia os contatos entre fluidos, o que permite um zoneamento básico do perfil de poço. A segunda propriedade decorre da relação quasi-linear entre o perfil de indução e a condutividade da formação, o que torna possível o uso da teoria dos sistemas lineares e, particularmente, o desenho de filtros digitais adaptados à deconvolução do sinal original. A idéia neste trabalho é produzir um algoritmo capaz de identificar os contatos entre as camadas atravessadas pelo poço, a partir da condutividade aparente lida pelo perfil de indução. Para simplificar o problema, o modelo de formação assume uma distribuição plano-paralela de camadas homogêneas. Este modelo corresponde a um perfil retangular para condutividade da formação. Usando o perfil de entrada digitalizado, os pontos de inflexão são obtidos numericamente a partir dos extremos da primeira derivada. Isto gera uma primeira aproximação do perfil real da formação. Este perfil estimado é então convolvido com a função resposta da ferramenta gerando um perfil de condutividade aparente. Uma função custo de mínimos quadrados condicionada é definida em termos da diferença entre a condutividade aparente medida e a estimada. A minimização da função custo fornece a condutividade das camadas. O problema de otimização para encontrar o melhor perfil retangular para os dados de indução é linear nas amplitudes (condutividades das camadas), mas uma estimativa não linear para os contatos entre as camadas. Neste caso as amplitudes são estimadas de forma linear pelos mínimos quadrados mantendo-se fixos os contatos. Em um segundo passo mantem-se fixas as amplitudes e são calculadas pequenas mudanças nos limites entre as camadas usando uma aproximação linearizada. Este processo é interativo obtendo sucessivos refinamentos até que um critério de convergência seja satisfeito. O algoritmo é aplicado em dados sintéticos e reais demonstrando a robustez do método.

Migração FFD 3D em profundidade usando aproximação de Padé complexa

Implementações dos métodos de migração diferença finita e Fourier (FFD) usam fatoração direcional para acelerar a performance e economizar custo computacional. Entretanto essa técnica introduz anisotropia numérica que podem erroneamente posicionar os refletores em mergulho ao longo das direções em que o não foi aplicado a fatoração no operador de migração. Implementamos a migração FFD 3D, sem usar a técnica do fatoração direcional, no domínio da frequência usando aproximação de Padé complexa. Essa aproximação elimina a anisotropia numérica ao preço de maior custo computacional buscando a solução do campo de onda para um sistema linear de banda larga. Experimentos numéricos, tanto no modelo homogêneo e heterogêneo, mostram que a técnica da fatoração direcional produz notáveis erros de posicionamento dos refletores em meios com forte variação lateral de velocidade. Comparamos a performance de resolução do algoritmo de FFD usando o método iterativo gradiente biconjugado estabilizado (BICGSTAB) e o multifrontal massively parallel direct solver (MUMPS). Mostrando que a aproximação de Padé complexa é um eficiente precondicionador para o BICGSTAB, reduzindo o número de iterações em relação a aproximação de Padé real. O método iterativo BICGSTAB é mais eficiente que o método direto MUMPS, quando usamos apenas um termo da expansão de Padé complexa. Para maior ângulo de abertura do operador, mais termos da série são requeridos no operador de migração, e neste caso, a performance do método direto é mais eficiente. A validação do algoritmo e as propriedades da evolução computacional foram avaliadas para a resposta ao impulso do modelo de sal SEG/EAGE.