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

Quantum oscillations in a π-striped superconductor

M. Zelli, Catherine Kallin, and A. John Berlinsky
Phys. Rev. B 86, 104507 – Published 7 September 2012

Abstract

Within Bogoliubov–de Gennes theory, a semiclassical approximation is used to study quantum oscillations and to determine the Fermi surface area associated with these oscillations in a model of a π-striped superconductor, where the d-wave superconducting order parameter oscillates spatially with period 8 and zero average value. This system has a nonzero density of particle-hole states at the Fermi energy, which form Landau-like levels in the presence of a magnetic field B. The Fermi surface is reconstructed via Andreev-Bragg scattering, and the semiclassical motion is along these Fermi surface sections as well as between them via magnetic breakdown. Within the approximation, oscillations periodic in 1/B are found in both the positions and widths of the lowest Landau levels. The area corresponding to these quantum oscillations for intermediate pairing interaction strength is similar to that reported for experimental measurements in the cuprates. A comparison is made of this theory to data for quantum oscillations in the specific heat measured by Riggs et al.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
27 More
  • Received 12 January 2012

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

©2012 American Physical Society

Authors & Affiliations

M. Zelli* and Catherine Kallin

  • Department of Physics and Astronomy, McMaster University, Hamilton, Ontario, Canada

A. John Berlinsky

  • Department of Physics and Astronomy, McMaster University, Hamilton, Ontario, Canada and Perimeter Institute for Theoretical Physics, Waterloo, Ontario, Canada

  • *zellim@mcmaster.ca

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 86, Iss. 10 — 1 September 2012

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

Authorization Required


×

Images

  • Figure 1
    Figure 1
    One-particle spectral functions at zero energy for the period 8 pair density wave system. The left panels show the spectral weight for adding or removing an electron (or hole) at E=0 in the extended BZ while the right panels show the same spectral weight folded into the reduced BZ and then repeated, which illustrates the complex shape of the Fermi surface. All panels are at 1/8 doping for (a) Δ=0.02 and (b) Δ=0.25.Reuse & Permissions
  • Figure 2
    Figure 2
    Low-lying energy levels for, from top to bottom, modulated d wave ignoring the effects of vortices, modulated d wave including the effects of vortices, and uniform d wave including the effects of vortices. Here Δ=0.25, μ=0.3, and L=256.Reuse & Permissions
  • Figure 3
    Figure 3
    The occupation function in the reduced BZ for 3 cases from left to right: (a) Δ=0, V=0, and μ=0.23; (b) Δ=0.25, V=0, and μ=0.3; (c) Δ=0, V=0.7, and μ=0.05. V is the magnitude of the interactions in CDW with periodicity of 4 sites. The reduced BZ for the CDW spans from 0 to π/2. Here, however, we have folded its FS once for easier comparison to the π-striped superconductor.Reuse & Permissions
  • Figure 4
    Figure 4
    FS for 3 cases from left to right: (a) Δ=0, V=0, and μ=0.23; (b) Δ=0.25, V=0, and μ=0.3; (c) Δ=0, V=0.7, and μ=0.05. V is the magnitude of the interactions in CDW with periodicity of 4 sites. The reduced BZ for the CDW spans from 0 to π/2. Here, however, we have folded its FS once for easier comparison to the π-striped superconductor.Reuse & Permissions
  • Figure 5
    Figure 5
    Low-energy DOS in the presence of a magnetic field of L=1024 for 3 cases: (a) Δ=0, V=0, and μ=0.23; (b) Δ=0.25, V=0, and μ=0.3; (c) Δ=0, V=0.7, and μ=0.05.Reuse & Permissions
  • Figure 6
    Figure 6
    Semiclassical motion of a nearly free particle system in the presence of a weak periodic potential (a) and a weak periodic superconducting pairing potential [(b) and (c)]. The direction of the semiclassical motion for particles is shown by black arrows. Holes (shown by red arrows) precess in the opposite direction. The gray area in the center figure is (c/eH)2(ATAb) where Ab is the area of the small electron pocket in panel (a) and AT is the area of the original circular FS. Starting from the blue cross in panel (b), the particle can either go over the whole unperturbed circular orbit by tunneling at points G, I, B, and C, or tunnel only at points G and I and Andreev scatter twice at points B and C covering the gray area. Another possible path is to Andreev scatter at points G and I and tunnel at points B and C. However, this path covers the same gray area. The change in the phase of the wave function is cAT/eH when the particle goes over the whole circular circuit and c(ATAb)/eH+β when it travels around the shaded area, where β is the phase shift due to two consecutive Andreev scatterings and is assumed to be relatively field independent. This behavior should be contrasted to that of the linked orbit of Pippard, shown on the left, where the particles orbit around the areas AT and Ab. Thus, as discussed in the text, the areas associated with quantum oscillations in the width of the first LL are different for the periodic potential and the periodic pairing models. Panel (c) shows the closed orbit corresponding to four successive Andreev-Bragg scatterings.Reuse & Permissions
  • Figure 7
    Figure 7
    Comparison of the low-energy DOS of a π-striped superconductor with Δ=0.02 and μ=0.23 in the presence of a magnetic field of L=256 with and without vortices, as described in the text.Reuse & Permissions
  • Figure 8
    Figure 8
    Comparison of the low-energy bands for the BdG Hamiltonian of a π-striped superconductor with Δ=0.02 and μ=0.23 in the presence of a magnetic field of L=256. The solid curves are the bands for the full BdG Hamiltonian, including vortices. The dashed curves are the semiclassical results for no vortices, and the flat lines (dash-dotted lines) are the Landau levels in the limit Δ=0.Reuse & Permissions
  • Figure 9
    Figure 9
    The width of the LL closest to E=0 as a function of 1/B or L for Δ=0.02 and μ=0.23, corresponding to 1/8 doping. 1/B is written in terms of the lattice constant a and flux quantum φ0. The solid line is a spline fit to the data that shows the oscillatory behavior more clearly. The inset shows the first LL for L=256.Reuse & Permissions
  • Figure 10
    Figure 10
    Power spectrum associated with the oscillations of the width and position of the first LL for Δ=0.02 at 1/8 doping. The x axis is rescaled so that it corresponds to area in units of the area of BZ.Reuse & Permissions
  • Figure 11
    Figure 11
    (a) Boomerang-shaped FS orbit involving two Andreev-Bragg scatterings and two tunnelings, as shown schematically in Fig. 6b, but for a period 8 modulation. The area of this orbit is denoted ATAb in the text. (b) The corresponding area Ab. (c) The area 2AbAT, corresponding to the difference of figures (a) and (b).Reuse & Permissions
  • Figure 12
    Figure 12
    Comparison of the geometrical area Ab and the area associated with quantum oscillations in the width of the first LL for Δ=0.02 vs μ in the region around 1/8 doping.Reuse & Permissions
  • Figure 13
    Figure 13
    Comparison of the geometrical area, ATAb, the boomerang-shaped area in Fig. 11a, and the area associated with oscillations in the position of the first LL, corresponding to the highest frequency peak in Fig. 10, shown as a function of μ.Reuse & Permissions
  • Figure 14
    Figure 14
    Power spectrum for oscillations of the position of the lowest LL for small values of the pairing potential amplitude, Δ. As discussed in the text, the peak at 0.0625 corresponds to the area AT, the original FS. The peak at 0.103 corresponds to ATAb, the boomerang-shaped area shown in Fig. 11a, while the feature at 0.1065 corresponds to the orbit with area 2AbAT, shown in Fig. 11c.Reuse & Permissions
  • Figure 15
    Figure 15
    Comparison of the low-energy DOS of a π-striped superconductor in the presence of a magnetic field of L=1024 with Δ=0.25 and μ=0.3 corresponding to 1/8 doping with and without vortices.Reuse & Permissions
  • Figure 16
    Figure 16
    The low-energy DOS for Δ=0.4 and L=800 at half filling. Each (double) peak has twice the degeneracy of a LL.Reuse & Permissions
  • Figure 17
    Figure 17
    Half width of the peak closest to E=0 for different values of Δ at half filling. The Fermi surfaces for two of the Δ values in this figure are shown in Fig. 18.Reuse & Permissions
  • Figure 18
    Figure 18
    Areas consistent with the quantum oscillations seen in the width of the first peak in the low-energy DOS are shown in red (dark-shaded) for two values of Δ at half filling. Note that for μ=0 the gray (light-shaded) areas have the same area as the red areas.Reuse & Permissions
  • Figure 19
    Figure 19
    Power spectrum associated with the oscillations in the width for Δ=0.25 and Δ=0.4 at half filling. The x axis is rescaled so that it corresponds to area in units of the area of BZ.Reuse & Permissions
  • Figure 20
    Figure 20
    Comparison of the geometrical area (red or gray area in Fig. 18) and the area associated with quantum oscillations in the width of the lowest energy peak for different values of Δ at half filling.Reuse & Permissions
  • Figure 21
    Figure 21
    Semilog plot of the width of the first LL for Δ=0.4 at half filling as a function of 1/B showing a fairly linear average behavior for not very large fields. This is expected if the broadening is caused by magnetic breakdown. The dashed line is a linear fit to the data.Reuse & Permissions
  • Figure 22
    Figure 22
    Power spectrum associated with the position of the first LL for Δ=0.6 and μ=0.5. The inset shows the position of the first LL for the same parameters.Reuse & Permissions
  • Figure 23
    Figure 23
    FS for Δ=0.6 and μ=0.5. The difference in the area of the gray (light-shaded) and red (dark-shaded) areas gives rise to the strongest peak in the power spectrum of the position of the first LL.Reuse & Permissions
  • Figure 24
    Figure 24
    Position of the first LL for Δ=0.5 and μ=0.4.Reuse & Permissions
  • Figure 25
    Figure 25
    FS for Δ=0.5 and μ=0.4. The red area is associated with short-period oscillations in Fig. 24 for larger magnetic fields and the gray (light-shaded) area is associated with the oscillations in the width of the first LL when magnetic breakdown occurs. The difference in the areas of the two lobes of the figure-eight shape corresponds to long-period oscillations in Fig. 24 at smaller fields. Black (thin) and red (thin) arrows show the two possible semiclassical paths.Reuse & Permissions
  • Figure 26
    Figure 26
    Position and width of the first peak for Δ=0.25 and μ=0.3, corresponding to 18 doping, plotted versus 1/B.Reuse & Permissions
  • Figure 27
    Figure 27
    Power spectrum for Δ=0.25 and μ=0.3.Reuse & Permissions
  • Figure 28
    Figure 28
    FS for Δ=0.25 and μ=0.3 in the quadrant of the first BZ.Reuse & Permissions
  • Figure 29
    Figure 29
    Comparison of the geometrical area and the area associated with quantum oscillations in the width and position of the first LL as a function of μ for Δ=0.2. The geometrical area is the area corresponding to the red (dark-shaded) region in Fig. 28 in the case of Δ=0.2.Reuse & Permissions
  • Figure 30
    Figure 30
    The spectra associated with oscillations in the width and position for Δ=0.25 and t2=0.15 at 18 doping. The peaks correspond to the gray (light-shaded) and red areas shown in Fig. 31. The results are consistent with those for t2=0.Reuse & Permissions
  • Figure 31
    Figure 31
    The areas associated with the first peaks of the position and width spectra in Fig. 30 for Δ=0.25 and t2=0.15 at 18 doping.Reuse & Permissions
  • Figure 32
    Figure 32
    Specific heat versus 1/B for Δ=0.25 and μ=0.3 and t2=0 for different temperatures. Temperatures in units of the hopping term t are shown on the right. Note the π phase shift in the oscillatory behavior of specific heat as T increases through T*0.003t.Reuse & Permissions
  • Figure 33
    Figure 33
    The oscillatory part of the calculated specific heat for Δ=0.25 and μ=0.34 with a zero second-nearest-neighbor hopping shown as a function of the magnetic field and temperature. To plot the data, t=0.16 eV is chosen.Reuse & Permissions
  • Figure 34
    Figure 34
    The oscillatory parts of the specific-heat data by Riggs et al. and calculations for a π-striped superconductor at T=1 K. The left y-axis scale is for the experimental data and the right one is for the model.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
×