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, H_N, é construída por meio do produto direto dos autovetores do Hamiltoniano H_n-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, S^z, prática comum em outros métodos, usufruímos da conservação com o quadrado do spin total, S². 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.

In this Thesis we develop a numerical method to diagonalize the Heisenberg model iteratively. In essence, we diagonalize spin chains in steps, each one corresponding to an addition of a spin to a smaller chain. The basis vectors to calculate the Hamiltonian of a N-spin chain, H_N, is built by means of the direct product of the eigenvectors of the (N-1)-spin Hamiltonian, diagonalized on the previous step, by the eigenstates of the N-th added spin. Besides the common use of the conservation of the z-component of the total spin, S^z, we also exploit the conservation of the squared total spin, S². For a specific class of spin systems we also implemented the reflection symmetry. We obtain the entire energy spectrum of spin-1/2 chains up to 20 sites, for which we show the temperature dependence of the magnetic susceptibility and specific heat, for systems with substitutional impurity spins, bond defects, border effects, i.e., for systems without translational invariance. This normally imposes enormous restrictions in many traditional methods. In order to diagonalize chains with a larger number of sites we implemented a procedure that selects lower energy states to be used in the basis vector on the next step. Using this truncation scheme, we are able to obtain low-lying energy states for chains with more than a hundred sites, up to five significant figures of accuracy. Our results reproduce those of the literature for the known cases, in general homogeneous systems. The approaches recently developed in the context of Density Functional Theory to the Heisenberg model, which also apply to inhomogeneous systems, are consistent with our numerical results. We generalize the method to diagonalize spin-1/2 ladders. We calculate the ground-state and the energy gap of this system, for arbitrary ratio of the couplings along the lags or over the rungs of the ladder. We also present the implementation of our iterative method to the Hubbard model, which describes a system of itinerant spins. It is known that in the regime of high Coulomb repulsion between the spins and unitary density (number of spins equal to the number of sites in the chain), this model is mapped onto Heisenberg one, a result which is verified numerically in our procedure by calculating the energy spectrum of both models in na appropriated parametric regime.