Neste trabalho usamos o Método de Monte Carlo para estudar o modelo quântico de Heisenberg ferromagnético. Foi utilizada a fórmula de Suzuki-Trotter para transformar uma cadeia fechada de N spins (S=1/2) em um sistema clássico equivalente bidimensional com interações de 4 spins, tipo Ising. A decomposição de tabuleiro foi utilizada no mapeamento. Um algoritmo de alteração global de spin, usado anteriormente por Cullen e Landau, foi empregado para estudar o sistema clássico resultante. Obtivemos resultados para a energia, susceptibilidade e calor específico que foram comparados com resultados publicados por vários autores. Uma lei de potência foi encontrada para susceptibilidade de » t-g; g = 1.61 ± 0.02 e para o comprimento de correlação e » t- v; V = 0.60 ± 0.01
Neste trabalho usamos o Método de Monte Carlo para estudar o modelo quântico de Heisenberg ferromagnético. Foi utilizada a fórmula de Suzuki-Trotter para transformar uma cadeia fechada de N spins (S=1/2) em um sistema clássico equivalente bidimensional com interações de 4 spins, tipo Ising. A decomposição de tabuleiro foi utilizada no mapeamento. Um algoritmo de alteração global de spin, usado anteriormente por Cullen e Landau, foi empregado para estudar o sistema clássico resultante. Obtivemos resultados para a energia, susceptibilidade e calor específico que foram comparados com resultados publicados por vários autores. Uma lei de potência foi encontrada para susceptibilidade de » t-g; g = 1.61 ± 0.02 e para o comprimento de correlação e » t- v; V = 0.60 ± 0.01
A transicao de fase do modelo de Ising com ligacoes mistas (Mixed-bond) numa rede cubica e estudada tanto numericamente com analiticamente. Analiticamente estudos do modelo de Ising mixed-bond atraves de tecnica de grupo de renormalizacao predizem a existencia de reentrâncias magneticas em certos valores do parâmetro de competicao α (α positivo). Este fenomeno e tipicamente encontrado em sistemas (a exemplo do composto EupSr1-p) que apresentam fase vidro-de-spin (spin-glass). Neste trabalho, usamos o algoritmo cluster de Wolff para simular a dinâmica do sistema e obter as quantidades termodinâmicas magnetizacao e susceptibilidade magnetica. A temperatura critica do sistema foi estimada do maximo de susceptibilidade e com estes valores construimos o diagrama de fase Tc versus p para diferentes valores de α ( 0 > α ). Nossos resultados foram comparados como aqueles obtidos pelo metodo de teoria de campo efetivo que emprega uma distribuicao de probabilidade similar dentro de uma estrutura de cluster com dois spins.
O modelo de Ising dipolar bidimensional inclui, além da interação ferromagnética entre os primeiros vizinhos, interações de longo alcance entre os momentos de dipolo magnético dos spins. A presença da interação dipolar muda completamente o sistema, apresentando um rico diagrama de fase, cujas características têm originado inúmeros estudos na literatura. Além disso, a possibilidade de explicar fenômenos observados em filmes magnéticos ultrafinos, os quais possuem diversas aplicações em àreas tecnológicas, também motiva o estudo deste modelo. O estado fundamental ferromagnético do modelo de Ising puro é alterado para uma série de fases do tipo faixas, as quais consistem em domínios ferromagnéticos de largura h com magnetizações opostas. A largura das faixas depende da razao delta das intensidades dos acoplamentos ferromagnético e dipolar. Através de simulações de Monte Carlo e técnicas de repesagem em histogramas múltiplos identificamos as temperaturas críticas de tamanho finito para as transições de fase quando delta=2, o que corresponde a h=2. Calculamos o calor específico e a susceptibilidade do parâmetro de ordem, no intervalo de temperaturas onde as transições são observadas, para diferentes tamanhos de rede. As técnicas de repesagem permitem-nos explorar e identificar máximos distintos nessas funções da temperatura e, desse modo, estimar as temperaturas críticas de tamanho finito com grande precisão. Apresentamos evidências numéricas da existência de uma fase nemática de Ising para tamanhos grandes de rede. Em nossas simulações, observamos esta fase para tamanhos de rede a partir de L=48. Para verificar o quanto a interação dipolar de longo alcance afeta as estimativas físicas, nós calculamos o tempo de autocorrelação integrado nas séries temporais da energia. Inferimos daí quão severo é o critical slowing down (decaimento lento crítico) para esse sistema próximo às transições de fase termodinâmicas. Os resultados obtidos utilizando um algoritmo de atualização local foram comparados com os resultados obtidos utilizando o algoritmo multicanônico.
Desde a descoberta de Stanley de que o modelo esférico é o limite de dimensionalidade de spin infinita do modelo de Heisenberg clássico, o modelo esférico tem atraído mais atenção do que meramente ser uma alternativa ao modelo de Ising. Como é bem conhecido, o modelo esférico clássico contradiz a terceira lei da Termodinâmica. Este defeito pode ser corrigido através de uma quantização do modelo. Nesta tese, escolhemos a proposta de Obermair que permite o controle das flutuações quânticas através do parâmetro g. Nesta tese, utilizamos o modelo esférico quântico de spin como laboratório da influência das flutuações quânticas sobre o comportamento magnético e termodinâmico da rede hipercúbica d-dimensional e das cadeias AB e AB2 antiferromagnéticas. A cadeia AB2 tem como motivações experimentais polímeros ferrimagnéticos e compostos de oxocupratos, dentre outras. Através da análise de observáveis como magnetização, susceptibilidade, comprimento de correlação e funções de correlação obtivemos no caso da rede hipercúbica que o modelo apresenta quebra espontânea de simetria a temperatura finita e para o estado fundamental (T=0). Através de cálculos detalhados dos observáveis mencionados acima pudemos obter os expoentes críticos associados às transições de fase térmica e quântica (T=0) da rede hipercúbica d-dimensional. Verificamos também que o calor específico em baixas temperaturas não viola a terceira lei da termodinâmica. Para o caso da cadeia AB, o estado fundamental antiferromagnético é obtido com a presença, e posterior anulamento, de campo magnético staggered e nas ausências de flutuações térmicas e quânticas e de campo uniforme. Por outro lado, o modelo esférico na cadeia AB3 apresenta estado fundamental ferrimagnético para g=T=h=0 mesmo na ausência de campo staggered. Entretanto, apesar das cadeias AB e AB2 não apresentarem transição de fase térmica ou quântica, os valores locais de spin da cadeia AB podem apresentar um comportamento antiferromagnético de curto alcance para flutuações térmicas e quânticas pouco intensas e baixos valores do campo magnético uniforme quando comparados com os do campo magnético staggered; no caso da cadeia AB2, destacamos comportamento ferrimagnético de curto alcance para fracas flutuações quânticas e térmicas e campo magnético baixo. Esses comportamentos antiferromagnético e ferrimagnético desaparecem com o aumento das flutuações térmicas e quânticas ou da presença do campo magnético. Tais resultados foram confirmados pelo comportamento assintótico das funções de correlação. Finalmente, obtivemos que os comportamentos em baixas temperaturas do modelo esférico quântico podem ser comparados com aquele da cadeia Heisenberg quântica XXZ, através da análise de observáveis como entropia e calor específico. A termodinâmica de altas temperaturas do modelo esférico quântico se assemelha à de um oscilador harmônico quântico, cuja freqüência depende da temperatura e do parâmetro quântico g
Nesta tese desenvolvemos um método numérico para diagonalizar o Hamiltoniano de Heisenberg iterativamente. O método consiste basicamente em diagonalizar cadeias de spins, cada vez maiores, em que cada passo da diagonalização corresponde à adição de um novo spin à cadeia. A base de vetores para calcular o Hamiltoniano de uma cadeia de N spins, HN, é construída por meio do produto direto dos autovetores do Hamiltoniano Hn-1 da rede diagonalizada no passo anterior, pelos autoestados correspondentes ao N-ésimo spin adicionado. Além de usar a comutação do Hamiltoniano com a componente azimutal do spin total, Sz, prática comum em outros métodos, usufruímos da conservação com o quadrado do spin total, S2. Para uma classe específica de redes também implementamos a simetria de reflexão. Obtemos o espectro completo de energia de cadeias de spins 1/2 com até 20 sítios, para as quais mostramos resultados da dependência com a temperatura da susceptibilidade magnética e do calor específico, para redes com impurezas tipo spin substitucionais, com defeitos nas ligações ou com efeitos de bordas, isto é, para sistemas sem invariância translacional. Usualmente essa restrição impõe enormes dificuldades em métodos tradicionais. Para diagonalizar cadeias com um número maior de sítios, implementamos um procedimento que seleciona os estados de mais baixa energia para serem usados na base de vetores do passo seguinte. Com esse tipo de truncamento de estados, fomos capazes de obter o estado fundamental e alguns estados de baixa energia de cadeias com mais de uma centena de sítios, com precisão de até cinco algarismos significativos. Nossos resultados reproduzem os da literatura para os casos conhecidos, em geral sistemas homogêneos. As aproximações desenvolvidas recentemente no contexto da Teoria do Funcional da Densidade, aplicada ao modelo de Heisenberg, e que também se aplicam a sistemas inomogêneos, estão em conformidade com nossos resultados numericamente exatos. Generalizamos o método para diagonalizar escadas de spins 1/2. Calculamos o estado fundamental e o gap de energia desse sistema, onde variamos a razão entre os acoplamentos ao longo das pernas da escada e ao longo dos degraus da mesma; nossos resultados são comparados com os da literatura. Apresentamos também a implementação do método iterativo no modelo de Hubbard, que descreve um sistema de spins itinerantes. Sabe-se que no regime de alta repulsão Coulombiana entre os spins e densidade um (número de spins igual ao número de sítios da cadeia), esse modelo é mapeado no modelo de Heisenberg, resultado que é verificado numericamente em nosso procedimento por meio do cálculo de energias de ambos os modelos em um regime paramétrico apropriado.
Investigações da dinâmica crítica em modelos de magnetismo, para tempos curtos, têm aparecido com grande freqüência na literatura. Essa técnica foi descoberta por Li, Schülke e Zheng que, inspirados em trabalhos anteriores de Huse e Janssen et al., mostraram que generalizações de grandezas como a magnetização e o cumulante de Binder exibem comportamento universal já no início da simulação. O estudo da criticalidade em tempos curtos proporciona um caminho alternativo para a estimativa do expoente z, além de permitir o cálculo de um novo expoente dinâmico θ, associado ao comportamento anômalo da magnetização. Da mesma forma, simulações dependentes do tempo tornaram-se ferramenta útil para estudar transições de fase em autômatos celulares e modelos de spin. Em particular, as melhores estimativas para o expoente z do Ising bidimensional foram obtidas por meio da técnica de propagação de danos, introduzida por Kauffman no estudo de autômatos e mais tarde generalizada para modelos de spin. Na primeira parte deste trabalho utilizamos o método Monte Carlo em tempos curtos para investigar o modelo de Baxter-Wu, definido em uma rede bidimensional triangular com variáveis do tipo Ising, acopladas por interações de três corpos. Obtivemos os expoentes críticos dinâmicos z e θ além dos índices críticos estáticos ß e Nû. Os resultados não corroboram aqueles recentemente obtidos por Santos e Figueiredo para o expoente z. Na segunda parte do trabalho, investigamos a propagação de danos no modelo de Ising unidimensional submetido a duas dinâmicas propostas por Hinrichsen e Domany (HD). Em particular, nós estudamos o efeito da atualização síncrona (paralela) e assíncrona (dinâmica contínua) sobre o espalhamento do dano. Mostramos que o dano não se propaga quando a segunda dinâmica é implementada de forma assíncrona. Também mostramos que as regras para atualização do dano produzidas por essa dinâmica, quando a temperatura vai a infinito e um certo parâmetro Lambda é igual a zero, são equivalentes àquelas do bem conhecido autômato celular (modelo A) de Grassberger.
Nesta tese estudamos vários tópicos ligados a simulações Monte Carlo de sistemas clássicos de spin em rede. Estamos interessados mais nos métodos do que nos resultados de aplicações com grande precisão numérica, devido a limitações computacionais. A ênfase é dada principalmente ao método de Grupo de Renormalização via Monte Carlo. Em primeiro lugar levamos a cabo um estudo detalhado do truncamento no espaço de Hamiltonianas para poder controlar numericamente os efeitos de operadores marginais. Um estudo detalhado do modelo Ashkin Teller N= 2 em duas dimensões é apresentado. Procuramos, a seguir, entender melhor o método de Ferrenberg e Swendsen de histogramas estudando-o para poder calcular, em uma só simulação, valores de expoentes não universais. Apresentamos resultados do modelo de Ising com defeito em 2d. Este método é aplicado ao problema da determinação da ordem de transições. Exemplos de modelos de Ising com interação de multispin são apresentados. Mostramos a seguir uma nova e poderosa técnica de investigar transições de fase que é obtida da combinação de idéias do Grupo de Renormalização via Monte Carlo e do método de histograma. Estes métodos são finalmente usados para estudar o modelo Ashkin Teller N= 3 anisotrópico em duas dimensões