Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Spin transport through a spin-12 XXZ chain contacted to fermionic leads

Florian Lange, Satoshi Ejima, Tomonori Shirakawa, Seiji Yunoki, and Holger Fehske
Phys. Rev. B 97, 245124 – Published 14 June 2018

Abstract

We employ matrix-product state techniques to numerically study the zero-temperature spin transport in a finite spin-12 XXZ chain coupled to fermionic leads with a spin bias voltage. Current-voltage characteristics are calculated for parameters corresponding to the gapless XY phase and the gapped Néel phase. In both cases, the low-bias spin current is strongly suppressed unless the parameters of the model are fine tuned. For the XY phase, this corresponds to a conducting fixed point where the conductance agrees with the Luttinger-liquid prediction. In the Néel phase, fine tuning the parameters similarly leads to an unsuppressed spin current with a linear current-voltage characteristic at low bias voltages. However, with increasing the bias voltage, there occurs a sharp crossover to a region where the current-voltage characteristic is no longer linear and a smaller differential conductance is observed. We furthermore show that the parameters maximizing the spin current minimize the Friedel oscillations at the interface, in agreement with the previous analyses of the charge current for inhomogeneous Hubbard and spinless fermion chains.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
1 More
  • Received 23 April 2018

DOI:https://doi.org/10.1103/PhysRevB.97.245124

©2018 American Physical Society

Physics Subject Headings (PhySH)

Condensed Matter, Materials & Applied Physics

Authors & Affiliations

Florian Lange1,2, Satoshi Ejima1,2, Tomonori Shirakawa2,3,4,5, Seiji Yunoki2,4,5, and Holger Fehske1

  • 1Institut für Physik, Ernst-Moritz-Arndt-Universität Greifswald, 17489 Greifswald, Germany
  • 2Computational Condensed Matter Physics Laboratory, RIKEN Cluster for Pioneering Research (CPR), Saitama 351-0198, Japan
  • 3International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
  • 4Computational Quantum Matter Research Team, RIKEN Center for Emergent Matter Science (CEMS), Saitama 351-0198, Japan
  • 5Computational Materials Science Research Team, RIKEN Center for Computational Science (R-CCS), Hyogo 650-0047, Japan

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 97, Iss. 24 — 15 June 2018

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1

    Schematic depiction of the setup defined by the Hamiltonian Ĥ0+ĤV according to Eqs. (1) and (6). Blue (green) circles indicate the spin chain (left and right leads). The red dashed line denotes the spin bias potential, which linearly decreases inside the spin chain.

    Reuse & Permissions
  • Figure 2
    Figure 2

    Time evolution of the spin current jjz(τ) in a junction composed of an isotropic spin chain (Δ=1) of NS=8 sites (shaded region) coupled to tight-binding leads of NL=100 sites for spin bias V/t=1 and several values of J/t. (a) Spin current profile at three different times τt=10, 20, and 30. (b) Time dependence of the spin current j0z(τ) between the spin chain and the left lead (solid lines) and estimated steady-state value (dashed lines). The result for J/t=2.4 with a larger size of the leads NL=500 is indicated by the green line. In the inset, the amplitude a1 of the current oscillations [see Eq. (9)] is shown for several different lead sizes NL. The solid line is a fit to a11/NL.

    Reuse & Permissions
  • Figure 3
    Figure 3

    Steady-state current jz as a function of J/t for four different values of Δ. Other parameters are NS=8, NL=500, and V/t=0.2. The dashed line shows the current V/(4π) expected for a Luttinger liquid with smooth interfaces.

    Reuse & Permissions
  • Figure 4
    Figure 4

    Magnetization profile Ŝjz around the interface for an applied magnetic field h/J=0.05. The dashed line indicates the interface between the tight-binding lead (j0) and the spin chain (j>0). The systems sizes are NL=400 and NS=400 for Δ=1 and NL=400 and NS=800 for Δ=2. The inset in (f) is a magnified view of the region close to the interface, highlighting the Friedel oscillations with wave number π. Solid black lines in (a) and (c) are fits of the data in the form of Eq. (10). The results are obtained by finite temperature calculations at the inverse temperature β.

    Reuse & Permissions
  • Figure 5
    Figure 5

    Strength of the Friedel oscillations OF defined in Eq. (11) around the interface inside the spin chain at the inverse temperature β. The system sizes are the same as in Fig. 4.

    Reuse & Permissions
  • Figure 6
    Figure 6

    Current-voltage curve in the two-lead setup described in Eq. (1) for different parameters in the gapless phase. The dashed line is the conductance G=1/(4π) of a Luttinger liquid smoothly connected to noninteracting leads.

    Reuse & Permissions
  • Figure 7
    Figure 7

    Current-voltage curve for a spin chain with NS=8, anisotropy Δ=2, and different values of J/t. The dashed line corresponds to the ideal conductance G=1/(4π) obtained in the Luttinger-liquid regime.

    Reuse & Permissions
  • Figure 8
    Figure 8

    Same as in Fig. 7 but for a spin chain with J/t=1.7 and various chain lengths.

    Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review B

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×