Coulomb interactions make electrodynam tansport of Dirac semi

Coulomb interactions make electrodynam tansport of Dirac semi – metals nonloca Baruch Rosenstein Collaborators: Theory: Hsien-Chong Kao (Nat. Taiwan Normal University) Meir Lewkowicz (Ariel University, Israel) Experiment: Wen – Bin Jian, Jian-Jhong Lai (Nat. Chiao Tung University) Tours, June, 13, 2016

Outline 1. Free relativistic massless fermions in graphene: finite conductivity without either carriers or impurities. 2. Problems with calculation of the electron – electron interaction effects and how they were resolved. 3. Nonlocal electrodynamics in graphene. 4. Experiment. 5. 3 D Weyl - semi metals.

Carbon atoms in graphene arranged in honeycomb 2 D crystal by the covalent bonds between nearest neighbours. E. Andrei et al, Nature Nano 3, 491(08)

1. Finite conductivity and no screening in noninteracting graphene Suspended samples are clearly undoped and exhibit the minimal resistivity at zero temperature. Andrei et al, Nature Nano 3, 491(08) On optical frequencies conductivity is the same Nair et al, Science 102 10451(08)

Tight binding model Spectrum consists of two bands with Fermi surface pinched right between them: not an obvious band insulator or a metal.

Fermi surface consists of two non - equivalent points of the Brillouin zone with sufficient little group to support a two dimensional representation Around the Hamiltonian is that of the left and right two component Weyl spinors:

The Dirac point in band structure was pointed out very early by Wallace: electrons in graphene should behave like a 2 D analogue of relativistic massless particles. The L, R Weyl spinors can be combined into two 4 Dirac spinor described by that possesses the chiral (sublattice) symmetry

DOS of relativistic massless fermions For the spectrum the number of electrons is Wang et al, APL 101, 183110 (12) At Dirac point the charge of carriers density flips sign: no “free” carriers at the neutrality (Dirac) point.

Explanation: Schwinger’s electron – hole pair creation by electric field The basic picture of the pseudo – diffusive resistivity in pure graphene is the creation of the electron – hole pairs by electric field. The pairs carry current that can be further increased by reorientation of moving particles Linear response for current

Heuristic interpretation of the conductivity The pair can be created with back to back momenta So that the current is proportional to electric field, number of available states and to the life time of the virtual pair Therefore

Dynamical approach to linear response demonstrating the absence of scale separation Let’s try to understand qualitatively how massless fermions react on electric field by just switching a homogeneous electric field and observing the creation of electron - hole pairs and the induced charge motion. Lewkowicz, B. R. , PRL 102, 106802 (09) One indeed observes that the current stabilizes at finite value. In DC this requires the use of tight binding model rather than its Dirac continuum theory.

In addition to the finite term converging to , there is also a linear acceleration term that vanishes due to cancellation of part near the Dirac points and far away. It can be presented as a full derivative Kao, Lewkowicz, B. R. , PRB 81, 041416 (2010) Therefore there is no perfect scale separation when massless fermions are involved. This is also the source of the chiral anomaly and doubling.

Accidental nature of the “pseudo-Ohmic” linear response B. R, Lewkowicz, PRL 77, (09) ; B. R. , Kao, Lewkowicz, PRB 81, R 041416 (10) Stays there till a crossover time to a linearly increasing Schwinger’s regime

At each one solves, using the WKB approximation, a tunneling problem similar to that in the Landau – Zener transition. Gavrilov, Gitman, PRD 53, 7162 (1986) The interband transition probability at large times is Schwinger, PR 82, 664 (1951) Hence Schwinger, PR (1962)

Coulomb repulsion potential remains long range The vacuum polarization Is too weak to exponentially screen the 3 D Coulomb with 2 D Fourier transform The relation between the polarization and conductivity obeys the usual charge conservation relation that is sometimes naively written as.

2. Electron – electron interactions effect on conductivity and screening They break the (pseudo) relativistic invariance and are still 3 D, namely less long range in 2 D than that of the 2+1 dimensional massless QED.

How strong are the interactions? Coulomb interactions naively are nonperturbatively strong due to coupling constant of order 1 Just rescale the field in Even if screened this would lead to a strongly coupled electronic system that for example might exhibit the chiral (exciton) condensate

This is very similar to the scalar relativistic interaction that exhibits chiral phase transition of non – Wislon universality class This would make quasiparticle massive and the phase insulating. Lattice simulation indeed demonstrated formation of the condensate above Drut, Lande, PRB 79, 79167425 (09) Ulyubishev et al PRL 111, 056801 13). , …

However in the previous discussion, quite complete understanding of transport was achieved ignoring the interactions. The interaction effects have been observed, but are small. Siegel et al, PRL. 110, 146802 (13) The first definitive measurement of the coupling (on the BN substrate) gave the value. The actual expansion parameter is , small enough to use perturbation theory. Elias et al, Science 337, 1196 (12)

Conflicting perturbation theory results in continuum limit (Dirac approximation) The leading correction to conductivity would be of a value, and can be directly calculated via Kubo formula by calculating diagrams First attempt in which sharp momentum cutoff was used gave a disappointingly large correction indicating that interactions are not weak Herbut, Juricic, Vafek, PRL 100 , 046403 (08)

However subsequent calculation of the polarization provided a different value. The charge conservation leads to a general relation between the charge and current transport functions: For a regular (means local) isotropic time reversal invariant (no magnetic field) system, one can expand Mishchenko, PRL 98, 216801)07); EPL 100, 046403 (08)

This value was considered to be favorable phenomenologically. However in order to force the two quantities, the polarization and the “direct” conductivities to be the same, the interaction was modified with the long range part modified to depend on UV cutoff. It was claimed however that this way the Ward identities of the charge conservation were violated by the cutoff procedure The polarization was recalculated with Ward identities obeyed by using a variant of dimensional regularization with yet different answer Juricic, Vafek, Herbut, PRB 82 , 235402 (10) The new calculations indicated problems with the dimensional regularization and reaffirmed the polarization value Mac. Donald et al, PRB 84, 045429 (11); Sodemann, Fogler, PRB 86, 045429 (12)

Calculation on the lattice Note that unlike in relativistic QFT, in graphene (or numerous other “Dirac” systems in condensed matter the natural regularization exists: a microscopic theory. To resolve the ambiguity we have performed the calculation within the tight binding model (lattice 3 D with Coulomb interactions between orbitals). B. R. , Maniv, Lewkowicz, PRL 110, 066602 (13); B. R. , H. C. Kao, Lewkowicz, PRB 90, 045137 (2014) The result for the conductivity is plus conservation of charge gives , while the polarization How to reconcile this?

The way out is to reconsider the basic assumptions involved in the derivation of the customary relations Undoped graphene is not a “regular” system due to lack of separation of scales and an additional nonanalytic term appears

3. Local vs. nonlocal electrodynamics When and why the Ohm’s low (locality) is obeyed If a conducting system has a space scale like the mean free path in metals with significant disorder, one expects smooth limit leading for an isotropic time reversal invariant (no magnetic field) system to For a very clean metal the situation might in principle become complex since on the free electron level there is no obvious scale. The strongly screened interactions might be sometimes neglected (or “renormalized away” in the Fermi liquid theory), but the electric response generally becomes nonlocal.

Knows it all… However it is possible to show that in any free lattice model conductivity is local since the difference can be again written as a full derivative and thus vanishes B. R. , H. C. Kao, in preparation (2016)

Long range interactions are required for a conductor to ne nonlocal Short range interactions metals or semi-metals do not change the situation, the screening length providing the necessary length scale. However in Weyl semi – metals (at neutrality point), despite the pseudo-Ohmic conductivity, Coulomb interaction remains long range. This creates a unique situation of the interaction correction creating the nonlocal electric response characterized in the rotation, time reversal invariant case by two scalar constants

For actual computation the trace of the conductivity tensor generally is simpler with an exact relation to polarization being The naïve locality assumption is thus violated

and , From the perturbative values one obtains

, and Extention to other Dirac semi-metals and universality For a Weyl semi-metal the contributions to nonlocality due to interactions is sum of all the nodes: graphene – like, double – graphene like… The integrals converge well in continuum after the appropriate surface term is identified. For the double layer graphene node one obtains: No known “topological explanation for these numbers How this additional term affects physics?

Potentials for the nonlocal electrodynamics Generally a vector field can be decomposed using two potentials One deduces the voltage between two points via sources

4. Experimental observation Anomalous dependence on geometry of electrical resistance in a rectangular flake The correction to the electrical resistance becomes dependent on geometry of the source, drain and points between which the voltage is measured V 21 Source Drain r 2 r 1 Without nonlocal electrodynamics due to interactionsinteraction

The correction to the electrical resistance becomes site dependent y C X=0 D y=W/2 A B X=L y=- W/2 x

An anomalous Hall effect y C X=0 D y=W/2 A B X=L y=- W/2 x

The asymmetric flake set up VTS Vd graphene on Si. O 2 substrate VBS Vg IDS Constant Current Source & 300 K Voltage Meter Mobility: 4000 ~10000 cm 2/V-s

Current (m. A) S D B T S D B Voltage (m. V) Total (m. V) (VDS=9. 04) 10 VTS=4. 87 VDT=4. 09 8. 96 10 T Volage (m. V) VBS=4. 53 VDB=4. 47 9

Nonlocality voltage near the Dirac Point T=300 K T=80 K

Temperature Dependence low. T temperature destroys the nonlocality

Disorder also destroys the nonlocality strong nonlocality weak nonlocality wide carrier range narrow carrier range

Anomalous reflection and transmission The additional term changes in a different way the two polarizations of the electromagnetic wave passing (bouncing off) a graphene layer The and the polarization components have the following reflection coefficients

5. 3 D Dirac semi-metals Topological transition from TI to nontopological insulator was shown to result in 3 D Dirac semi – metal:

Na 3 Bi Liu ZK et al, Science (2014)

Dirac or Weyl points in 3 D appeared occasionally in band structure calculations away from Fermi surface, but not always. Properties of a (poor) metal with strong spin – orbit, Bi, led to its (approximately) Dirac Hamiltonian with pseudospin replaced by spin with all three Pauli matrices now presentpoint theory

For free quasi-particles in addition to the real part B. R. , Lewkowicz PRB 88, 045108 (13) , In many cases of interest the imaginary part dominates over small dissipation.

Electrodynamics of such a material becomes rather peculiar

The interaction correction The correction as in graphene renormalizes coefficients of both the real and the imaginary parts More importantly the difference between the longitudinal and the transverse conductivities is. B. R. , Lewkowicz PRB 88, 045108 (13) ,

The charging phenomenon Fluxon induces charges proportional to

Incident light splits into the transversal and the longitudinal ones, both very weakly attenuated. L T Moreover at certain angle and frequency the monochromatic wave appears that is both transverse and longitudinal.

For certain angle the reflection disappears: total absorption phenomenon

Summary 1. Dirac semi-metal like graphene at Dirac point is a neutral plasma that has finite conductivity, but does not makes the interactions between quasiparticles short range. 2. There is no complete scale separation of the energy scales also in the interacting case and some quantities like the conductivity tensor are nonanalytic. Low energy continuum Weyl model requires regularization that is sufficiently precise to avoid effects of chiral symmetry violation by the cutoff. 3. Transport and electromagnetic properties of WSM are unusual including interaction effects leading to nonlocal electrodynamics resulting in several optical and electric peculiarities.

A radial and a rotational current in a graphene layer with Corbino geometry exhibits different electrical resistance: S D B(t)