Time-dependent properties of one-dimensional spin-systems : a DMRG-study

Friedrich, Andreas; Schollwöck, Ulrich

doi:1674

Time-dependent properties of one-dimensional spin-systems : a DMRG-study = Zeitabhängige Eigenschaften eindimensionaler Spin-Systeme : eine DMRG Anwendung

Friedrich, Andreas (Author)

2006

Verantwortlichkeitsangabevorgelegt von Andreas Friedrich

ImpressumAachen : Publikationsserver der RWTH Aachen University 2006

Umfang102 S. : Ill., graph. Darst.

Aachen, Techn. Hochsch., Diss., 2006

Genehmigende Fakultät
Fak01

Hauptberichter/Gutachter
Schollwöck, Ulrich (Thesis advisor)

Tag der mündlichen Prüfung/Habilitation
2006-10-27

Online
URN: urn:nbn:de:hbz:82-opus-16747
URL: https://publications.rwth-aachen.de/record/61435/files/Friedrich_Andreas.pdf

Einrichtungen

Fakultät für Mathematik, Informatik und Naturwissenschaften (100000)

Inhaltliche Beschreibung (Schlagwörter)
Physik (frei)

Thematische Einordnung (Klassifikation)
DDC: 530

Kurzfassung
In den letzten Jahren sind große Fortschritte bei der experimentellen Untersuchung dynamischer Größen niedrigdimensionaler Systeme erzielt worden. Dies trifft sowohl für verschiedene Realisierungen von Quantenpunkten und -drähten zu als auch für Systeme, die durch Anisotropien niedrigdimensionale Strukturen enthalten, wie beispielsweise organische Strukturen wie NTENP und NDMAP, in denen Spinketten gefunden wurden. Die theoretische Beschreibung dieser Systeme ist nach wie vor eine Herausforderung. In dieser Arbeit wurde die Methode der Dichtematrix-Renormierungsgruppe (DMRG) benutzt, um zeitabhängige Größen eindimensionaler Spinsysteme zu untersuchen. Hierzu wurden zwei verschiedene Arten von DMRG-Algorithmen verbessert beziehungsweise verallgemeinert. Der zeitaufwändigste Schritt bei der Anwendung der sogenannten Korrekturvektor-DMRG zur Bestimmung dynamischer Strukturfaktoren, das Lösen eines komplexen linearen Gleichungssystems A|x> = |b> mit nicht hermitescher Matrix A, konnte durch die Implementierung des sogenannten „GMRES”-Algorithmus deutlich beschleunigt werden. Mit dieser verbesserten Version der Korrekturvektor-DMRG wurde die Spektralfunktion der Spin-1-Heisenberg-Kette im starken äußeren Magnetfeld untersucht (s. unten). Der zweite verwendete DMRG-Algorithmus ist die erst kürzlich entwickelte adaptive zeitabhängige DMRG (adaptive t-DMRG). Sie erlaubt die Simulation der Zeitentwicklung eines gegebenen quantenmechanischen Vielteilchenzustandes. In ihrer ursprünglichen Formulierung ist sie allerdings auf Hamilton-Operatoren beschränkt, die nur sehr kurzreichweitige Wechselwirkungen enthalten. In dieser Arbeit wurden Verallgemeinerungen basierend auf Ideen von Steven R. White implementiert, die eine Verwendung beliebiger Hamilton-Operatoren gestatten. Es wurden verschiedene Zeitentwicklungsalgorithmen getestet und verglichen. Als Anwendung der adaptiven t-DMRG wurde die Zeitentwicklung eines einzelnen Elektronenspins in einem Halbleiter-Quantenpunkt berechnet, der mit den nuklearen Spins des Gitters wechselwirkt (s. unten). Eines der grundlegenden Konzepte der Festkörperphysik ist das der Tomonaga-Luttinger-Flüssigkeit (TLL). Obwohl sie in vielen unterschiedlichen Systemen existiert, ist ihre direkte Beobachtung schwierig. Einer der Gründe für diese Schwierigkeiten ist die Tatsache, dass die Anregungen des Systems, die die TLL bilden, keine Anregungslücke haben. Daher finden sich ihre Beiträge zu dynamischen Strukturfaktoren bei extrem kleinen Energien, die experimentell schlecht zugänglich sind. In der Spin-1-Heisenberg-Kette ist die energetisch niedrigste Anregung durch eine Energielücke vom Grundzustand getrennt. Es handelt sich um ein S=1-Triplett, das von einem von außen angelegten Magnetfeld aufgespalten wird. Die Energie der Sz=1-Magnonen wird mit steigendem Feld abgesenkt und schließt die Anregungslücke, wenn das Feld einen kritischen Wert erreicht. Das System durchläuft dann einen Quantenphasenübergang und es bildet sich eine TLL aus Sz=1-Magnonen. Im dynamischen Strukturfaktor bildet sich oberhalb des TLL-Kontinuums ein weiteres Kontinuum aus, das von den Sz=0-Magnonen aus dem S=1-Triplett gebildet wird, die als mobile Verunreinigungen der TLL verstanden werden können. Ihr Beitrag zur Spektralfunktion kann aufgrund starker Wechselwirkung zwischen ihnen und der darunter liegenden TLL einige TLL-Eigenschaften zeigen. In dieser Arbeit wurde das Hochenergie-Kontinuum im longitudinalen dynamischen Strukturfaktor untersucht. Es konnte gezeigt werden, dass es für Energien in der Größenordnung der Haldane-Anregungslücke Delta durch eine Singularität dominiert wird, ähnlich der im TLL-Kontinuum bei kleinen Energien. Es wurden sowohl der zugehörige Exponent alpha' als auch die Randfrequenz des Kontinuums als Funktion des Magnetfeldes untersucht. Der Exponent alpha' zeigt dabei ein deutlich anderes Verhalten als der des TLL-Kontinuums alpha. Während letzterer mit steigendem Magnetfeld zunimmt, sinkt alpha' ab. Ausserdem konnte die effektive Wechselwirkungsstärke zwischen den Sz=0-Magnonen des Hochenergie-Kontinuums und den $S^z=1$-Magnonen der TLL als Funktion des Magnetfeldes bestimmt werden. Besonderes Interesse gilt auf dem Gebiet der Quanteninformationsverarbeitung möglichen physikalischen Realisierungen von Quantenbits (Qubits). Eine der vielversprechendsten Realisierungen besteht aus einem in einem Quantenpunkt eingeschlossenen Elektron, dessen zwei Spinzustände als natürliche Qubit-Zustände betrachtet werden. Wichtig ist eine lange Dekohärenzzeit des Qubitzustands, so dass möglichst viele Rechenoperationen durchgeführt werden können, bevor der Zustand zerfällt. Daher ist ein tieferes Verständnis der Mechanismen, die zu Dekohärenz führen, unerlässlich. In dieser Arbeit wurde die Zeitentwicklung eines elektronischen Spins untersucht, der mit den nuklearen Spins im Halbleitergitter des Quantenpunktes wechselwirkt. Generell hängt der Verlauf der Zeitentwicklung stark vom verwendeten Anfangszustand ab. Einige dieser Zustände sind in dieser Arbeit untersucht worden. Zunächst wurde der Langzeitlimes der z-Komponente des Spins als Funktion der nuklearen Magnetisierung und eines Magnetfeldes betrachtet. Er durchläuft für beide Parameter ein Minimum wenn der jeweils andere konstant gehalten wird. Dieses Verhalten konnte durch die im System herrschende Energieerhaltung erklärt werden. Als eine Möglichkeit, die Dekohärenzzeit eines Zustandes in einem Quantenpunkt zu verlängern, wurde vor kurzem vorgeschlagen, den Zustand des elektronischen Spins in den nuklearen Spins zu speichern. Das Protokoll funktioniert perfekt, wenn die Kopplungen zwischen dem Elektronenspin und allen Kernspins des Halbleiters gleich stark sind und das nukleare Spinbad anfangs vollständig polarisiert ist. Beide Voraussetzungen sind in realen Systemen nicht gegeben. Hier wurde die Güte F des Algorithmus untersucht. Die Abhängigkeit von F von der nuklearen Polarisation wurde für zwei verschiedene Verteilungen der Kopplungskonstanten untersucht. Wie erwartet zeigen beide eine mit sinkender Polarisation stark abfallende Güte F. Dennoch sind deutliche Unterschiede zu erkennen. Für eine der beiden Verteilungen der Kopplungskonstanten wurde das Verhalten als Funktion der Systemlänge betrachtet. Anders als in der Literatur angenommen findet man eine starke Längenabhängigkeit von F. Die Güte nimmt mit wachsender Systemlänge deutlich ab. Dies ist für reale Systeme mit ca 100000 Nuklei ein ernsthaftes Problem. Daher sind für experimentelle Anwendungen dieses Zwischenspeichers für elektronische Spinzustände enorme nukleare Polarisationen nötig, die mit dem heutigen Stand der Technik nicht erreichbar sind.

In this thesis we have applied the DMRG method to calculate dynamical and time-dependent properties of one-dimensional spin systems. Two different kinds of DMRG algorithms have been applied, and improved or generalized. In the correction vector DMRG used to calculate dynamical structure factors, the most time-consuming step is to solve a system of complex linear equations A|x>=|b> with A and |b> given where the matrix A is non-hermitian. In this work we have used the so-called "GMRES"-algorithm so solve for real and imaginary part of |x> simultaneously. Since this avoids squaring the matrix it leads to a speed-up of the algorithm. We have applied this improved correction vector DMRG to a spin-$1$ Heisenberg chain in a strong external magnetic field. The results of this study are summarized below. The second kind of DMRG algorithm that has been used in this thesis is the recently developed adaptive time-dependent DMRG (adaptive t-DMRG). It allows us to calculate the time-evolution of a given quantum many-body state under the action of a certain time-independent system hamiltonian. In its original formulation the adaptive t-DMRG is restricted to hamiltonians containg only short range interactions because the time-evolution operator has to be split up into local operators via a Suzuki-Trotter decomposition. In this work we have implemented a generalization of this algorithm to hamiltonians containing arbitrary interactions basd on ideas by Steven R. White. Several different time-evolution algorithms and sets of target states have been used and compared. The optimal choice among the algorithms tested is a Krylov subspace algorithm allowing for extremely long time steps and hence leading to short calculation times. In this form the adaptive t-DMRG shows at least comparable performance on systems where also Suzukmi-Trotter algorithms could be used. We have used the adaptive t-DMRG to calculate the time-evolution of a single electronic spin in a semiconductor quantum dot due to the hyperfine interaction with the nuclear spins of the underlying lattice. The findings are summarized below. The results of the applications of these different DMRG algorithms presented in this work can be summarized as follows: A Tomonaga-Luttinger liquid (TLL) is one of the paradigmatic concepts in solid state physics. Although it is believed to exist in several different systems its observation is difficult and most of the experimental evidence has been obtained indirectly. One reason for this is that the excitations forming a TLL are gapless and thus their contributions to dynamical structure factors are at extremely low energies. Higher energy features might be easier to measure. A spin-1 Heisenberg chain is a well-known example of a gapped antiferromagnet that, when driven by a magnetic field, undergoes a quantum phase transition to a TLL. In this regime the spectral function, apart from the TLL-continuum, exhibits a second band at higher energies. It is formed by excitations that act as mobile impurities and interact strongly with the excitations of the underlying Tomonaga-Luttinger liquid. Due to this strong interaction their contribution to dynamical structure factors are believed to bear TLL-properties like a power-law decay and a non-universal edge-exponent. In this work we have shown that the high-energy continuum emerging in the longitudinal dynamical structure factor contains an edge singularity similar to the one existing in the low-energy response. We have demonstrated that this high energy part of the spectrum for energies of the order of the Haldane gap Delta is dominated by this power law edge singularity. We have studied the behavior of the corresponding exponent alpha' and its edge frequency omega as functions of the external magnetic field. The behavior of the exponent alpha' as a function of the magnetic field is very different from the one of the TLL edge exponent alpha of the low-energy continuum. While alpha grows with increasing field alpha' decreases. We have also been able to calculate the effective interaction strength between the Sz=0 magnons forming the high-energy continuum and the Sz=1 magnons forming the underlying TLL. The field of quantum information processing and quantum computation has been very active in recent years. Many proposals of solid state realizations of qubits have been made one of the most promising of which is to use the spin of a single electron confined in a semiconductor quantum dot as natural qubit. To actually perform quantum computing several pre-requisites have to be fulfilled. One of the most important is a long decohenerce time of the qubit states such that as many gate operations as possible can be performed. Thus a deeper understanding if the processes leading to decoherence is indispensable. We have investigated the time-evolution and the decoherence behavior of a single electronic spin in a semiconductor quantum dot due to the hyperfine interaction with the nuclear spins of the underlying lattice. In general the time-evolution strongly depends on the initial state chosen. Several such initial states have been studied. The quantity we have analyzed is the long-time limit of the electron spin as a function of an external magnetic field. We have found that it runs through a minimum as a function of the polarization of the nuclear spin system when the magnetic field is kept constant due to energy conservation . The same is true when the magnetic field B is changed at fixed nuclear polarization. In this case the value of B where the minimum is reached strongly depends on the polarization of the nuclear spin bath. One way to extend the decoherence time for electronic states in a quantum dot is to store these states in the nuclear spin bath, which has a much larger intrinsic decoherence time. The coupling between the nuclear spins and the electron induces oscillations in S^z(t). If all coupling constants between the electronic and the nuclear spins were equal and the nuclear spin bath was fully polarized the electronic state would be completely encoded in the nuclear spins after half a period. At that point the electron is ejected from the dot. After some time a new electron, which is in down state too, is injected into the dot. After another half oscillation the electronic state is fully recovered. Real baths are not fully polarized and the coupling constants are not equal. This leads to a loss of fidelity F. We have investigated the fidelity of this protocol as a function of the polarization of the nuclear spins. This has been done for two different distributions of coupling constants. Both sets of coupling constants qualitatively show the expected behavior. The corresponding fidelities decrease with decreasing nuclear polarization. Nevertheless strong differences can be realized. They refer to the slope of the mentioned decrease and the maximum reachable fidelity at full nuclear polarization. For one set of couplings we have investigated the dependence of the fidelity on the number of nuclear spins $L$ for fixed polarization of the nuclear spin bath. Contrary to what is generally assumed we find a rather strong L-dependence. The fidelity of the protocol decreases roughly linearly with increasing number of nuclear spins. For real quantum dots with approximately 100000 nuclei this is a serious problem. For an experimental realization a very large nuclear spin polarization is needed which is not reachable with today's technology.

OpenAccess:
PDF
(additional files)