BSc Final Year Report - Mat<strong>the</strong>w HolmesPrevious comparisons with exact values, in Figs. 8 & 10, hinted at <strong>the</strong> trend toward linearity. With nofield, Fig. 8, <strong>the</strong> ground state increases haphazardly with N. In a field of 10T (Fig. 10), when <strong>the</strong>variational estimates first agree for all N, <strong>the</strong> rise in energies shows more linearity. <strong>The</strong> field strengthin Fig. 10 is similar to <strong>the</strong> h=0.9 line in Fig. 13, above. For confidence in <strong>the</strong> estimations beyond 8spins, <strong>the</strong> field parameter h is increased for safety. In this way we can be confident that <strong>the</strong> degree ofaccuracy of <strong>the</strong> modified <strong>Lanczos</strong> estimates for small chains, is retained beyond N=8.Limitations of <strong>The</strong> Modified <strong>Lanczos</strong> Method<strong>The</strong> Modified <strong>Lanczos</strong> method is confined to making estimation of <strong>the</strong> ground state energies, for both<strong>the</strong> ferro- <strong>and</strong> antiferromagnetic cases. It is thus useful only in <strong>the</strong> zero temperature limit <strong>and</strong> cannotevaluate excited states of <strong>the</strong> chain, for which <strong>the</strong> full algorithm is required.ConclusionBy applying <strong>the</strong> <strong>Heisenberg</strong> Hamiltonian to antiferromagnetically coupled spins along a lineardimension, a model of a closed 1D spin chain has been developed. <strong>The</strong> justification for such a model ismade with reference to experimental observations of quasi-1D chains <strong>and</strong> a <strong>the</strong>oretical consideration of<strong>the</strong> nature of antiferromagnetic exchange interactions.<strong>The</strong> ground state eigenvalues of <strong>the</strong> Hamiltonian, calculated with <strong>the</strong> <strong>Lanczos</strong> method were found to bestrongly dependent on <strong>the</strong> eigenvalue spacing. Accuracy was found to be significantly increased byusing <strong>the</strong> modified <strong>Lanczos</strong> method, a variational technique. To achieve complete accuracy, however,required that <strong>the</strong> chain be placed in a magnetic field. With a field of ~10T a complete agreementbetween <strong>the</strong> <strong>Lanczos</strong> estimates <strong>and</strong> a direct evaluation with <strong>the</strong> Jacobi method, was found up to chainlength of 8 spins.Direct evaluation with Jacobi, for a chain of 8 spins was <strong>the</strong> limit on <strong>the</strong> author's computer. <strong>The</strong> fieldhad <strong>the</strong> effect of pulling apart <strong>the</strong> eigenvalues <strong>and</strong> allowing complete accuracy up to a chain of 8 spins.It also permitted larger systems to be treated than <strong>the</strong> 8 spin limit using <strong>the</strong> Jacobi method.<strong>The</strong> ability for accurate eigenvalue estimation with <strong>the</strong> <strong>Lanczos</strong> method, <strong>and</strong> its ability to treat largersystems, has thus been demonstrated.It is suggested that <strong>the</strong> excited states of <strong>the</strong> Hamiltonian be investigated using <strong>the</strong> full <strong>Lanczos</strong>algorithm, <strong>the</strong> modified method being confined to ground state estimation. <strong>The</strong> ground stateeigenvector might also be used to construct <strong>the</strong> operators with which to investigate a wider range ofphysical properties than have been considered herein.AcknowledgmentsI would like to thank my project supervisor, Dr Paul Stevenson, for his excellent guidance (<strong>and</strong> infinitepatience!) throughout. Now is also <strong>the</strong> time to acknowledge my debt to <strong>the</strong> sun; for shining throughmy window one day.References[1] Aharoni A., Introduction to <strong>the</strong> <strong>The</strong>ory of Ferromagnetism, Oxford (2000).[2] Brussard P. J. & Glaudemans P. W. M., Shell Model Applications in Nuclear Spectroscopy,North-Holl<strong>and</strong> Publishing Company, Amsterdam, (1977).[3] Chatelin F., Eigenvalues of Matrices, Wiley (1993).[4] Cloizeaux J. & Pearson J. J., Phys. Rev. 128, 2131 (1962).[5] Coldea R. et al., Phys. Rev. Lett. 86, 5377 (2001).[6] Dagotto E. & Moreo A., Phys. Rev. D 31, 865 (1985).[7] Dender D. C. et al., Phys. Rev. B 53, 2583 (1996).[8] Ellis T. M. R. et al., Fortan 90 Programming, Addison-Wesley (1994).[9] Hihilashvili R., Derivation of <strong>the</strong> <strong>Heisenberg</strong> Hamiltonian, http://phjoan5.technion.ac.il/~riki/,Israel Institute of Technology Dept. of Physics, (2003).[10] Hutchings M. T., Ikeda H., Milne J. M., J. Phys. C 12, L739 (1979)[11] Krupi ka S. et al, Elements of <strong>The</strong>oretical Magnetism, Iliffe (1968).[12] Lake B. et al, ISIS Science, www.isis.rl.ac.uk/science/index.htm, Ru<strong>the</strong>rford AppletonLaboratory, (2001).[13] Lieb E. H. et al., Ma<strong>the</strong>matical Physics in One Dimension, Academic Press (1996).16
BSc Final Year Report - Mat<strong>the</strong>w Holmes[14] Malvezzi A. L., Brazilian J. Phys. 33, 55 (2003).[15] Merzbacher E., Quantum Mechanics, Wiley (1998).[16] Nepomechie R. I., A Spin Chain Primer, www.physics.miami.edu/nepomechie/primer.pdf,University of Miami Phys. Dept. (1998).[17] O'Reilly E., Quantum <strong>The</strong>ory of Solids, Taylor & Francis (2002).[18] S<strong>and</strong>vik A. W., 'Computational Studies of Quantum Spin Systems' & 'Numerical Solutions of<strong>the</strong> Schrödinger Equation', PY502 Computational Physics, Dept. of Physics, BostonUniversity (2004).[19] Scaife B. K. P. et al, Studies in Numerical Analysis, Academic Press (1974)[20] Stevenson P. D., University of Surrey Phys. Dept., Private communication (2005).[21] Tyablikov S. V., Methods in <strong>The</strong> Quantum <strong>The</strong>ory of Magnetisation, Plenum Press (1967).17