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
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
Através da nova técnica de Busca Automática de parâmetros de ordem e propriedades críticas, estudamos vários sistemas físicos conhecidos, com o intuito de testar a eficiência desta técnica, a fim de aplicá-la para. sistemas complexos em geral.\nO Método de Busca Automática (MBA) é desenvolvido em combinação com o Método de Simulação de Monte Carlo e este por sua vez se utiliza de uma determinada dinâmica para. fazer o sistema evoluir com o tempo. Nos problemas que tratamos, utilizamos as dinâmicas de Metropolis e de Gláuber, que produzem resultados essencialmente equiva¬lentes.\nFazemos uma breve revisão das transições de fase e fenômenos críticos. Descrevemos o MBA e comentamos alguns trabalhos onde a. aplicação deste método já obteve sucesso.\nAplicamos o MBA para o modelo de Ising ferromagnético. Calculamos as fronteiras críticas e o expoente fl da magnetização em várias redes geométricas e várias dimensões espaciais.\nAplicamos também o MBA para o modelo de Ising ferromagnético em uma rede quadrada com diluição de sítios. Calculamos as fronteiras críticas, o expoente 0 e o expoente 7 da susceptibilidade. O MBA calcula simultaneamente os expoentes (5 e 7. Verificamos que a. classe de universalidade do sistema não é alterada quando se introduz diluição.\nEstudamos o problema, da. percolação de ligações de longo alcance na cadeia linear diluída. Discutimos questões de extensividade e não extensividade inerentes à sistemas com interações de longo alcance.\nFinalmente, apresentamos nossas conclusões e possíveis extensões deste trabalho
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 ν.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 λ é igual a zero, são equivalentes àquelas do bem conhecido autômato celular ( modelo A) de Grassberger.Agradecimentos Ao Prof. Dr. J. R. Drugowich de Felício pela orientação e amizade.À "tia" Claudinha pela amizade e sempre agradabilíssima companhia.À Renatinha, Andréa Grupioni e ao Leandro Maranghetti pelas constantes brigas, constantes discussões e principalmente pela sincera amizade desde os tempos da graduação em Química.Para todo pessoal da Eletroquímica e especialmente a "dona" Paulinha e a "tia" Marlene.À Tânia Lassali pela amizade e incrível paciência por me agüentar.À Eleade e a Rosângela Tenório pelas tentativas, sempre frustradas, de me endireitarem.Aos meus amigos advogados, sempre muito úteis, Elienai e Ricardo Santana.À bióloga, cantora, escritora, poetiza, motoqueira.... e amiga Pierângeli.Aos meio químicos, meio físicos, Cláudia e Parada.À futura psicóloga Paula Garducci por sempre esquecer de cobrar as minhas consultas.À Carol que com o seu irresistível sotaque baiano sempre alegra um pouco mais o meu dia.
Neste trabalho nós estudamos alguns sistemas de spins e vértices exatamente solúveis em duas dimensões. A solubilidade exata está ligada ao fato de existirem soluções não triviais das equações de fatorização, o que nos permite obter a energia livre no limite termodinâmico. Introduzimos e resolvemos pelo método de espalhamento inverso, um modelo de dez vértices assimétrico com dois e três estados nas ligações. Obtemos o diagrama de fases e mostramos que o sistema exibe uma transição de fase de primeira ordem. Analisamos um modelo de oito vértices de férmions livres e propomos uma nova relação funcional que nos permite calcular a energia livre por vértice. Mostramos que este sistema de vértices corresponde ao modelo de Ising na rede Union Jack. Apresentamos um método de solução de modelos de spin em redes triangulares a partir da solução do mesmo modelo na rede quadrada. O método se aplica sempre que o modelo de spins envolver interação de primeiros vizinhos e satisfizer a relação triângulo-estrela. Estendemos para a rede triangular, as soluções autoduais de Fateev e Zamolodchikov para a rede quadrada, de modelos de spin com simetria Z(N). Analisamos as conjecturas existentes sobre a criticalidade do modelo de Potts definido na rede de Kagomé. Baseados na simetria e nas degenerescências dessa rede conjecturamos uma expressão para a sua linha crítica.
Neste trabalho nós estudamos alguns sistemas de spins e vértices exatamente solúveis em duas dimensões. A solubilidade exata está ligada ao fato de existirem soluções não triviais das equações de fatorização, o que nos permite obter a energia livre no limite termodinâmico. Introduzimos e resolvemos pelo método de espalhamento inverso, um modelo de dez vértices assimétrico com dois e três estados nas ligações. Obtemos o diagrama de fases e mostramos que o sistema exibe uma transição de fase de primeira ordem. Analisamos um modelo de oito vértices de férmions livres e propomos uma nova relação funcional que nos permite calcular a energia livre por vértice. Mostramos que este sistema de vértices corresponde ao modelo de Ising na rede Union Jack. Apresentamos um método de solução de modelos de spin em redes triangulares a partir da solução do mesmo modelo na rede quadrada. O método se aplica sempre que o modelo de spins envolver interação de primeiros vizinhos e satisfizer a relação triângulo-estrela. Estendemos para a rede triangular, as soluções autoduais de Fateev e Zamolodchikov para a rede quadrada, de modelos de spin com simetria Z(N). Analisamos as conjecturas existentes sobre a criticalidade do modelo de Potts definido na rede de Kagomé. Baseados na simetria e nas degenerescências dessa rede conjecturamos uma expressão para a sua linha crítica.
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.
Sistemas desordenados vêm sendo objeto de extenso estudo devido ao fato de apresentarem características peculiares não presentes em sistemas ordenados.Entre estas características se encontra a controvenida questão da dependênciado valor de expoentes críticos com parâmetros não pertinentes no conceito tradicional de universalidade aplicada a sistemas uniformes. No presente trabalho, estudamos o comportamento do modelo de Ising generalizado com interaçôes leatórias entre primeiros vizinhos em duas dimensões(rede quadrada), através de simulações de Monte Carlo. Utilizamos uma rede quadrada 60x60 com condições de contorno periódicas e constantes de acoplamemo entre primeiros vizinhos distribuídas aleatoriamentede acordo com as distribuições de probabilidade bimodal e gaussiana. Variamos também os valores assumidos pelos spins: S1= ±1 e S1=0,±1. Este modelo só apresenta umatransição de fase em T=O,logo observamos uma divergência da susceptibilidade magnética em função do campo externo aplicado de acordo com uma lei de potência no limite T-+0 calculada para pequenos valores de campo.Observa-seque o expoentecritico associado 11 dependeda forma da distribuição de probabilidade dos Ji;porém para a mesma distribuição,independedo valor quântico do spin. Realizamos,em seguida,um estudo das propriedades do groundstate de vidros de spin unidimensionais através de resultados numéricos, obtidos utilizando relações de recorrência para a Hamiltonianado modelode Ising com interaçôes aleatórias em uma dimensão,na presença de z vizinhos. Utilizamos cadeias lineares com 10 spins e com constantes de acoplamentos seguindo distribuições de probabilidade uniforme, bimodale gaussiana. Estimamos o expoente critico Δ para diferentes distribuiçõese diferentes números de vizinhos. o que permitiu concluir que este sistema apresenta um comportamento dependente da distribuição de probabilidade, para certas distribuições, dependente do número de vizinhos