Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
  • Editors' Suggestion

Influence of local symmetry on lattice dynamics coupled to topological surface states

Jonathan A. Sobota, Samuel W. Teitelbaum, Yijing Huang, José D. Querales-Flores, Robert Power, Meabh Allen, Costel R. Rotundu, Trevor P. Bailey, Ctirad Uher, Tom Henighan, Mason Jiang, Diling Zhu, Matthieu Chollet, Takahiro Sato, Mariano Trigo, Éamonn D. Murray, Ivana Savić, Patrick S. Kirchmann, Stephen Fahy, David A. Reis, and Zhi-Xun Shen
Phys. Rev. B 107, 014305 – Published 12 January 2023
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • METHODS
  • RESULTS
  • DISCUSSION AND CONCLUSIONS
  • ACKNOWLEDGMENTS
  • APPENDICES
    • Supplemental Material
    References

    Abstract

    We investigate coupled electron-lattice dynamics in the topological insulator Bi2Te3 with time-resolved photoemission and time-resolved x-ray diffraction. It is well established that coherent phonons can be launched by optical excitation, but selection rules generally restrict these modes to zone-center wave vectors and Raman-active branches. We find that the topological surface state couples to additional modes, including a continuum of surface-projected bulk modes from both Raman and infrared branches, with possible contributions from surface-localized modes when they exist. Our calculations show that this surface vibrational spectrum occurs naturally as a consequence of the translational and inversion symmetries broken at the surface, without requiring the splitting-off of surface-localized phonon modes. The generality of this result suggests that coherent phonon spectra are useful by providing unique fingerprints for identifying surface states in more controversial materials. These effects may also expand the phase space for tailoring surface state wave functions via ultrafast optical excitation.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    2 More
    • Received 17 September 2021
    • Accepted 22 December 2022

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

    ©2023 American Physical Society

    Physics Subject Headings (PhySH)

    1. Research Areas
    Electron-phonon couplingPhonons
    1. Physical Systems
    Topological materials
    1. Techniques
    Time & angle resolved photoemission spectroscopyUltrafast pump-probe spectroscopyX-ray diffraction
    Condensed Matter, Materials & Applied Physics

    Authors & Affiliations

    Jonathan A. Sobota1,*, Samuel W. Teitelbaum2,3, Yijing Huang2,4, José D. Querales-Flores5, Robert Power6, Meabh Allen6, Costel R. Rotundu1, Trevor P. Bailey7, Ctirad Uher7, Tom Henighan2,4, Mason Jiang2,4, Diling Zhu8, Matthieu Chollet8, Takahiro Sato8, Mariano Trigo1,2, Éamonn D. Murray5, Ivana Savić5, Patrick S. Kirchmann1, Stephen Fahy5,6, David A. Reis1,2,4,9, and Zhi-Xun Shen1,4

    • 1Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
    • 2Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
    • 3Department of Physics, Arizona State University, Tempe, Arizona 85281, USA
    • 4Department of Applied Physics, Stanford University, Stanford, California 94305, USA
    • 5Tyndall National Institute, Lee Maltings, Dyke Parade, Cork T12 R5CP, Ireland
    • 6Department of Physics, University College Cork, College Road, Cork T12 K8AF, Ireland
    • 7Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
    • 8Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
    • 9Department of Photon Science, Stanford University, Stanford, California 94305, USA

    • *sobota@slac.stanford.edu

    Article Text (Subscription Required)

    Click to Expand

    Supplemental Material (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand
    Issue

    Vol. 107, Iss. 1 — 1 January 2023

    Reuse & Permissions
    Access Options
    CHORUS

    Article Available via CHORUS

    Download Accepted Manuscript
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      (a) Schematic of time-resolved ARPES measurements on Bi2Te3 to measure electronic binding energy oscillations. (b) ARPES spectrum along the ΓK¯ direction before and (c) after excitation at the pump-probe delays indicated. The bulk conduction band and most of the surface state are unoccupied in equilibrium, but become partially filled by the excitation. (d) Time-dependent shift in the binding energy of the bulk conduction band and (e) surface state within the energy-momentum windows indicated by dashed lines in (c). (f) Magnitude of the Fourier transform of the bulk conduction band and (g) surface state dynamics after subtracting a slowly varying background (6th order polynomial for the bulk and exponential for the surface). Points are from the data and lines are fits (see Sec. pp2 for fitting methodology and Table 1 for fit parameters). The arrows highlight two modes observed in the surface state but not in the bulk.

      Reuse & Permissions
    • Figure 2
      Figure 2

      (a) Time-resolved x-ray diffraction measurement to measure coherent lattice dynamics. [(b)–(d)] Time-dependent relative change in the intensities of the (hkl) = (4 4 5), (5 5 6), and (3 3 4) Bragg peaks. [(e)–(g)] The corresponding Fourier transforms after exponential background subtraction. Points are from the data, and solid lines are a fit. Two modes are observed, with (h+k+l)-dependent amplitudes consistent with modes of A1g1 and A1g2 symmetries.

      Reuse & Permissions
    • Figure 3
      Figure 3

      Theoretical calculation of surface and bulk lattice dynamics in Bi2Te3 and Bi2Se3 using DFT and a dynamical matrix formalism. (a) Fourier transforms of the Bi atomic motion at the surface and center (“bulk”) of the slab. (b) The calculated response of the surface state compared to the trARPES data for the surface state [reproduced from Fig. 1], showing overall agreement. (c) Γ¯-projected density of states of the full slab and for the top QL only. Arrows denote surface modes outside the bulk continuum. (d) Bulk phonon dispersion curves along Γ-Z. [(e)–(h)] Same set of analysis repeated for Bi2Se3, exhibiting similar phenomenology. Experimental data for Bi2Se3 in (f) reproduced from Ref. [29].

      Reuse & Permissions
    • Figure 4
      Figure 4

      Classical semi-infinite one-dimensional chain of trimers, representing the minimum model containing Raman and infrared phonons. (a) Cartoon of the model. All atoms have identical masses m, connected by spring constants k1/m=(2π)2(1THz)2 and k2/m=(2π)2(0.8THz)2. At t=0 the trimers are uniformly distorted throughout the chain (see Δx arrows) corresponding to excitation of the A1g optical phonon with qz=0. (b) Phonon dispersion and (c) density of states for the bulk. (d) Surface density of states and (e) FFT of resulting motion for a bulk (red) and surface (blue) atom for the case in which the surface k1 are softened by 2%. A damping time of 20 ps is used. (f) The amplitude of the low-frequency and high-frequency peaks relative to the central peak, plotted versus layer number. Corresponding plots for the unperturbed surface are shown in (g)–(i), and for the 2% hardened surface in (j)–(l).

      Reuse & Permissions
    • Figure 5
      Figure 5

      Complementary sensitivities of multimodal probes. Near the surface (within the top quintuple layer), atomic motion involves the surface-projected bulk phonon spectrum with infrared- and Raman-active branches in addition to surface-localized modes. These oscillations are probed by trARPES measurements of the surface state. Deep within the bulk, only Raman-active modes at qz=0 are excited. These are probed by trARPES measurements of the bulk states, as well as trXRD and transient optical reflectivity. The different frequency of the A1g2 mode in XRD is due to its different measurement temperature and fluence.

      Reuse & Permissions
    • Figure 6
      Figure 6

      Overview of the Fourier transform fitting methodology. (a) Real and imaginary parts of the Fourier transform of a simulated data set of two oscillators with equal amplitudes and opposite phases. The frequencies are f1=1.0 THz and f2=1.5 THz with lifetimes τ1=τ2=1 ps. (b) Magnitude of the same Fourier transform. The dashed lines depict the magnitude of each oscillator separately. Note that the total magnitude is not simply the sum of the individual magnitudes due to interference. (c) Simultaneous fit to the real and imaginary parts of the FFT for the bulk band dynamics measured by trARPES. (d) The FFT magnitude resulting from this fit, which is also shown in Fig. 1. [(e) and (f)] Same analysis but for the surface band.

      Reuse & Permissions
    • Figure 7
      Figure 7

      Structure factor analysis of trXRD data. [(a)–(f)] Time-resolved x-ray diffraction measurement for six Bragg peaks excited by 800 nm with an incident fluence of 8.2mJ/cm2. [(g)–(h)] Amplitude of both modes extracted by fitting the time-dependent intensities, plotted as a function of (h+k+l) for the six measured peaks. The solid line is a fit from a structure function model assuming modes of A1g1 and A1g2 symmetries, respectively. This allows for quantitative extraction of the displacements of all five atoms in the unit cell, sketched as arrows in (i) and (j).

      Reuse & Permissions
    • Figure 8
      Figure 8

      Fluence-dependence of mode frequencies measured by time-resolved XRD at room temperature. (a) Fourier transform of the (5 5 6) Bragg peak dynamics as a function of incident fluence. Solid curves are fits and dashed lines are guides to the eye. (b) Fluence-dependent frequencies extracted from the fits. The A1g1 mode is fluence-independent at 1.84 THz, while the A1g2 mode extrapolates to 3.94±0.05 THz at zero fluence with a slope of 0.040THz/(mJ/cm2).

      Reuse & Permissions
    • Figure 9
      Figure 9

      Time-resolved optical reflectivity measurement on Bi2Te3. (a) Relative time-dependent change in the reflectivity at 800 nm. (b) Fourier transform after background subtraction. Points are from the data, and solid lines are a fit. Two peaks corresponding to the A1g1 and A1g2 modes are observed at 1.85 and 3.89 THz.

      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
    ×