skip to main content

Title: Electrically and magnetically switchable nonlinear photocurrent in РТ-symmetric magnetic topological quantum materials

Nonlinear photocurrent in time-reversal invariant noncentrosymmetric systems such as ferroelectric semimetals sparked tremendous interest of utilizing nonlinear optics to characterize condensed matter with exotic phases. Here we provide a microscopic theory of two types of second-order nonlinear direct photocurrents, magnetic shift photocurrent (MSC) and magnetic injection photocurrent (MIC), as the counterparts of normal shift current (NSC) and normal injection current (NIC) in time-reversal symmetry and inversion symmetry broken systems. We show that MSC is mainly governed by shift vector and interband Berry curvature, and MIC is dominated by absorption strength and asymmetry of the group velocity difference at time-reversed ±kpoints. Taking$${\cal{P}}{\cal{T}}$$PT-symmetric magnetic topological quantum material bilayer antiferromagnetic (AFM) MnBi2Te4as an example, we predict the presence of large MIC in the terahertz (THz) frequency regime which can be switched between two AFM states with time-reversed spin orderings upon magnetic transition. In addition, external electric field breaks$${\cal{P}}{\cal{T}}$$PTsymmetry and enables large NSC response in bilayer AFM MnBi2Te4, which can be switched by external electric field. Remarkably, both MIC and NSC are highly tunable under varying electric field due to the field-induced large Rashba and Zeeman splitting, resulting in large nonlinear photocurrent response down to a few THz regime, suggesting bilayer AFM-zMnBi2Te4as a more » tunable platform with rich THz and magneto-optoelectronic applications. Our results reveal that nonlinear photocurrent responses governed by NSC, NIC, MSC, and MIC provide a powerful tool for deciphering magnetic structures and interactions which could be particularly fruitful for probing and understanding magnetic topological quantum materials.

« less
Award ID(s):
Publication Date:
Journal Name:
npj Computational Materials
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Complete theoretical understanding of the most complex superconductors requires a detailed knowledge of the symmetry of the superconducting energy-gap$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα, for all momentakon the Fermi surface of every bandα. While there are a variety of techniques for determining$$|{\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha |$$Δkα, no general method existed to measure the signed values of$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα. Recently, however, a technique based on phase-resolved visualization of superconducting quasiparticle interference (QPI) patterns, centered on a single non-magnetic impurity atom, was introduced. In principle, energy-resolved and phase-resolved Fourier analysis of these images identifies wavevectors connecting allk-space regions where$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαhas the same or opposite sign. But use of a single isolated impuritymore »atom, from whose precise location the spatial phase of the scattering interference pattern must be measured, is technically difficult. Here we introduce a generalization of this approach for use with multiple impurity atoms, and demonstrate its validity by comparing the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαit generates to the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαdetermined from single-atom scattering in FeSe where s±energy-gap symmetry is established. Finally, to exemplify utility, we use the multi-atom technique on LiFeAs and find scattering interference between the hole-like and electron-like pockets as predicted for$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαof opposite sign.

    « less
  2. Abstract

    We present a proof of concept for a spectrally selective thermal mid-IR source based on nanopatterned graphene (NPG) with a typical mobility of CVD-grown graphene (up to 3000$$\hbox {cm}^2\,\hbox {V}^{-1}\,\hbox {s}^{-1}$$cm2V-1s-1), ensuring scalability to large areas. For that, we solve the electrostatic problem of a conducting hyperboloid with an elliptical wormhole in the presence of anin-planeelectric field. The localized surface plasmons (LSPs) on the NPG sheet, partially hybridized with graphene phonons and surface phonons of the neighboring materials, allow for the control and tuning of the thermal emission spectrum in the wavelength regime from$$\lambda =3$$λ=3to 12$$\upmu$$μm by adjusting themore »size of and distance between the circular holes in a hexagonal or square lattice structure. Most importantly, the LSPs along with an optical cavity increase the emittance of graphene from about 2.3% for pristine graphene to 80% for NPG, thereby outperforming state-of-the-art pristine graphene light sources operating in the near-infrared by at least a factor of 100. According to our COMSOL calculations, a maximum emission power per area of$$11\times 10^3$$11×103W/$$\hbox {m}^2$$m2at$$T=2000$$T=2000K for a bias voltage of$$V=23$$V=23V is achieved by controlling the temperature of the hot electrons through the Joule heating. By generalizing Planck’s theory to any grey body and deriving the completely general nonlocal fluctuation-dissipation theorem with nonlocal response of surface plasmons in the random phase approximation, we show that the coherence length of the graphene plasmons and the thermally emitted photons can be as large as 13$$\upmu$$μm and 150$$\upmu$$μm, respectively, providing the opportunity to create phased arrays made of nanoantennas represented by the holes in NPG. The spatial phase variation of the coherence allows for beamsteering of the thermal emission in the range between$$12^\circ$$12and$$80^\circ$$80by tuning the Fermi energy between$$E_F=1.0$$EF=1.0eV and$$E_F=0.25$$EF=0.25eV through the gate voltage. Our analysis of the nonlocal hydrodynamic response leads to the conjecture that the diffusion length and viscosity in graphene are frequency-dependent. Using finite-difference time domain calculations, coupled mode theory, and RPA, we develop the model of a mid-IR light source based on NPG, which will pave the way to graphene-based optical mid-IR communication, mid-IR color displays, mid-IR spectroscopy, and virus detection.

    « less
  3. Abstract

    The proximity of many strongly correlated superconductors to density-wave or nematic order has led to an extensive search for fingerprints of pairing mediated by dynamical quantum-critical (QC) fluctuations of the corresponding order parameter. Here we study anisotropics-wave superconductivity induced by anisotropic QC dynamical nematic fluctuations. We solve the non-linear gap equation for the pairing gap$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm)and show that its angular dependence strongly varies below$${T}_{{\rm{c}}}$$Tc. We show that this variation is a signature of QC pairing and comes about because there are multiples-wave pairing instabilities with closely spaced transition temperatures$${T}_{{\rm{c}},n}$$Tc,n. Taken alone, each instability would produce a gap$$\Deltamore »(\theta ,{\omega }_{m})$$Δ(θ,ωm)that changes sign$$8n$$8ntimes along the Fermi surface. We show that the equilibrium gap$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm)is a superposition of multiple components that are nonlinearly induced below the actual$${T}_{{\rm{c}}}={T}_{{\rm{c}},0}$$Tc=Tc,0, and get resonantly enhanced at$$T={T}_{{\rm{c}},n}\ <\ {T}_{{\rm{c}}}$$T=Tc,n<Tc. This gives rise to strong temperature variation of the angular dependence of$$\Delta (\theta ,{\omega }_{m})$$Δ(θ,ωm). This variation progressively disappears away from a QC point.

    « less
  4. Abstract

    Hemiwicking is the phenomena where a liquid wets a textured surface beyond its intrinsic wetting length due to capillary action and imbibition. In this work, we derive a simple analytical model for hemiwicking in micropillar arrays. The model is based on the combined effects of capillary action dictated by interfacial and intermolecular pressures gradients within the curved liquid meniscus and fluid drag from the pillars at ultra-low Reynolds numbers$${\boldsymbol{(}}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{7}}}{\boldsymbol{\lesssim }}{\bf{Re}}{\boldsymbol{\lesssim }}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{3}}}{\boldsymbol{)}}$$(107Re103). Fluid drag is conceptualized via a critical Reynolds number:$${\bf{Re}}{\boldsymbol{=}}\frac{{{\bf{v}}}_{{\bf{0}}}{{\bf{x}}}_{{\bf{0}}}}{{\boldsymbol{\nu }}}$$Re=v0x0ν, wherev0corresponds to the maximum wetting speed on a flat, dry surface andx0is the extension length of the liquidmore »meniscus that drives the bulk fluid toward the adsorbed thin-film region. The model is validated with wicking experiments on different hemiwicking surfaces in conjunction withv0andx0measurements using Water$${\boldsymbol{(}}{{\bf{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{25}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{28}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v02m/s,25µmx028µm), viscous FC-70$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{0.3}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{18.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\boldsymbol{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{38.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v00.3m/s,18.6µmx038.6µm)and lower viscosity Ethanol$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{1.2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{11.8}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{33.3}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v01.2m/s,11.8µmx033.3µm).

    « less
  5. Abstract

    Strong electronic interactions can drive a system into a state with a symmetry breaking. Lattice frustration or competing interactions tend to prevent symmetry breaking, leading to quantum disordered phases. In spin systems frustration can produce a spin liquid state. Frustration of a charge degree of freedom also can result in various exotic states, however, experimental data on these effects is scarce. In this work we demonstrate how in a Mott insulator on a weakly anisotropic triangular lattice a charge ordered state melts on cooling down to low temperatures. Raman scattering spectroscopy finds that$$\kappa$$κ-(BEDT-TTF)$${}_{2}$$2Hg(SCN)$${}_{2}$$2Cl enters an insulating “dipole solid” statemore »at$$T=30\,{\mathrm{K}}$$T=30K, but below$$T=15\,{\mathrm{K}}$$T=15Kthe order melts, while preserving the insulating energy gap. Based on these observations, we suggest a phase diagram relevant to other quantum paraelectric materials.

    « less