IMPACT ionization and thermalization in photo-doped Mott insulators

Similar documents
The Hubbard model out of equilibrium - Insights from DMFT -

Dynamical mean field approach to correlated lattice systems in and out of equilibrium

Nonequilibrium Physics of Correlated Electron Materials IV: Nonequilibrium Phase Transitions

Electronic correlations in models and materials. Jan Kuneš

Diagrammatic Monte Carlo methods for Fermions

DMFT for correlated bosons and boson-fermion mixtures

Realistic Materials Simulations Using Dynamical Mean Field Theory

Diagrammatic Monte Carlo simulation of quantum impurity models

Surprising Effects of the Interaction between Electrons in Solids. Dieter Vollhardt

Magnetic Moment Collapse drives Mott transition in MnO

Spin and orbital freezing in unconventional superconductors

From Materials to Models and Back. Dieter Vollhardt

NiO - hole doping and bandstructure of charge transfer insulator

Dual fermion approach to unconventional superconductivity and spin/charge density wave

Carriers Concentration and Current in Semiconductors

Role of Hund Coupling in Two-Orbital Systems

Many body physics issues with the core-hole propagator and resonant inelastic X-ray scattering

Inelastic light scattering and the correlated metal-insulator transition

Organic Conductors and Superconductors: signatures of electronic correlations Martin Dressel 1. Physikalisches Institut der Universität Stuttgart

Cluster Extensions to the Dynamical Mean-Field Theory

Surprising Effects of Electronic Correlations in Solids

arxiv: v4 [cond-mat.str-el] 24 Jan 2018

Lattice modulation experiments with fermions in optical lattices and more

1 G. Kotliar: Lecture 2

Dynamical Mean Field Theory and Numerical Renormalization Group at Finite Temperature: Prospects and Challenges

Cooperative Phenomena

Theory of carbon-based magnetism

?What are the physics questions?

Quantum dynamics in many body systems

Electronic charge reconstruction of doped Mott insulators in multilayered nanostructures

FROM NODAL LIQUID TO NODAL INSULATOR

Lecture 1. OUTLINE Basic Semiconductor Physics. Reading: Chapter 2.1. Semiconductors Intrinsic (undoped) silicon Doping Carrier concentrations

From Gutzwiller Wave Functions to Dynamical Mean-Field Theory

w2dynamics : operation and applications

Engineering of quantum Hamiltonians by high-frequency laser fields Mikhail Katsnelson

Angle-Resolved Two-Photon Photoemission of Mott Insulator

Metal-insulator transition with Gutzwiller-Jastrow wave functions

Electronic Squeezing by Optically Pumped Phonons: Transient Superconductivity in K 3 C 60. With: Eli Wilner Dante Kennes Andrew Millis

Correlatd electrons: the case of high T c cuprates

arxiv:cond-mat/ v1 [cond-mat.str-el] 4 Jul 2000

O. Parcollet CEA-Saclay FRANCE

The Gutzwiller Density Functional Theory

Dynamics of fluctuations in high temperature superconductors far from equilibrium. L. Perfetti, Laboratoire des Solides Irradiés, Ecole Polytechnique

SOLID STATE PHYSICS. Second Edition. John Wiley & Sons. J. R. Hook H. E. Hall. Department of Physics, University of Manchester

Electron Correlation

arxiv:cond-mat/ v1 [cond-mat.str-el] 11 Aug 2003

ESE 372 / Spring 2013 / Lecture 5 Metal Oxide Semiconductor Field Effect Transistor

Bose-Hubbard Model (BHM) at Finite Temperature

3.1 Introduction to Semiconductors. Y. Baghzouz ECE Department UNLV

Dynamical Mean-Field Theory for Correlated Electron Materials Dieter Vollhardt

Green's Function in. Condensed Matter Physics. Wang Huaiyu. Alpha Science International Ltd. SCIENCE PRESS 2 Beijing \S7 Oxford, U.K.

Density matrix renormalization group study of a three- orbital Hubbard model with spin- orbit coupling in one dimension

Ultrafast surface carrier dynamics in topological insulators: Bi 2 Te 3. Marino Marsi

A typical medium approach to Anderson localization in correlated systems.

Lecture 2: Ultracold fermions

Spin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University

Out-of-equilibrium electron dynamics in photoexcited topological insulators studied by TR-ARPES

MOTTNESS AND STRONG COUPLING

Ultrashort Lifetime Expansion for Resonant Inelastic X-ray Scattering. Luuk Ament

Spectral Density Functional Theory

Nano-DMFT : the electronic structure of small, strongly correlated, systems

Classification of Solids

Field-Driven Quantum Systems, From Transient to Steady State

Interaction-induced Symmetry Protected Topological Phase in Harper-Hofstadter models

LEC E T C U T R U E R E 17 -Photodetectors

Unit III Free Electron Theory Engineering Physics

An efficient impurity-solver for the dynamical mean field theory algorithm

Strongly correlated multilayered nanostructures near the Mott transition

Mott transition : beyond Dynamical Mean Field Theory

Lecture 15: Optoelectronic devices: Introduction

Quantum Phase Transition

SEMICONDUCTOR PHYSICS REVIEW BONDS,

Introduction to DMFT

The Oxford Solid State Basics

SUPPLEMENTARY INFORMATION

The Peierls distortion seen in 1D chains: The simplest model for a gap.

Numerical Methods in Quantum Many-body Theory. Gun Sang Jeon Pyeong-chang Summer Institute 2014

Transfer of spectral weight in spectroscopies of correlated electron systems

EE495/695 Introduction to Semiconductors I. Y. Baghzouz ECE Department UNLV

Supraconductivité à haute température dans les cuprates et les organiques: Où en est-on?

Local moment approach to the multi - orbital single impurity Anderson and Hubbard models

Electronic and Optoelectronic Properties of Semiconductor Structures

Strongly Correlated Materials: Insights From Dynamical Mean-Field Theory

Relativistic magnetotransport in graphene

arxiv:cond-mat/ v1 26 Jan 1994

Linear temperature dependence of electron spin resonance linewidths in La 0.7 Ca 0.3 MnO 3 and YBaMn 2 O 6

An introduction to the dynamical mean-field theory. L. V. Pourovskii

Current mechanisms Exam January 27, 2012

First-order Mott transition at zero temperature in two dimensions: Variational plaquette study

Semiconductor Detectors

The Mott Metal-Insulator Transition

Semiconductor Physics. Lecture 3

Probing dynamical properties of Fermi-Hubbard systems with a quantum gas microscope Peter Brown Bakr Lab Solvay workshop, February 19 th 2019

Metal-Insulator Transitions at Surfaces

PN Junction

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS

Calculating Band Structure

Mott physics: from basic concepts to iron superconductors

Mean field theories of quantum spin glasses

Optical and transport properties of small polarons from Dynamical Mean-Field Theory

Transcription:

IMPACT ionization and thermalization in photo-doped Mott insulators Philipp Werner (Fribourg) in collaboration with Martin Eckstein (Hamburg) Karsten Held (Vienna) Cargese, September 16

Motivation Photo-doping: nonequilibrium phase transition from a correlation induced insulator to a non-thermal conducting state S. Iwai et al. (3), H. Okamoto et al. (7),... Thermalization of large-gap insulators Impact ionization in small-gap insulators Cooling by magnon scattering U Mott insulating solar cells

Model and method Hubbard model: simplified model for a correlated electron material Gutzwiller, Kanamori, Hubbard (1963) t U Sign problem / exponential scaling: lattice model not solvable use approximate description

Model and method Dynamical mean field theory DMFT: mapping to an impurity problem Metzner & Vollhardt (1989); Georges & Kotliar (199) lattice model impurity model t G latt G imp t k latt imp Formalism can be extended to nonequilibrium systems Schmidt & Monien (); Freericks et al. (6) Impurity solver: computes the dynamics on the correlated site Strong-coupling perturbation theory: Eckstein & Werner (9)

Model and method Equilibrium DMFT phase diagram (half-filling) Paramagnetic calculation: Metal - Mott insulator transition at low T Smooth crossover at high T T metal Mott insulator U

Pulse excited Mott insulator Photo-excitation of carriers across the Mott gap Eckstein & Werner (11) Question: How quickly does the electronic system thermalize? T metal Mott insulator U

Pulse excited Mott insulator Photo-excitation of carriers across the Mott gap Eckstein & Werner (11) Question: How quickly does the electronic system thermalize? pulse form total energy T e 4.1 E(t) 3 1-1 - -3-4 4 6 gap x gap 8 1 1 14 energy(t) -.1 -. -.3 -.4 -.5 -.6 4 6 8 1 1 14 t t

Pulse excited Mott insulator Photo-excitation of carriers across the Mott gap Eckstein & Werner (11) Question: How quickly does the electronic system thermalize?.18.17.16 thermal value.15 d(t).14.13.1.11.1 4 gap x gap 6 8 1 1 14 t

Pulse excited Mott insulator Photo-excitation of carriers across the Mott gap Eckstein & Werner (11) Question: How quickly does the electronic system thermalize? -1 log 1 d(t)-d(t eff ) - -3-4 U=5 U=3 U=.5 T -5 U=1.5 U= 5 1 15 t U

Pulse excited Mott insulator Photo-excitation of carriers across the Mott gap Eckstein & Werner (11) Strong correlation regime: Relaxation time depends exponentially on U -1 3 log 1 d(t)-d(t eff ) - -3-4 U=5 U=3 U=.5 log 1 relax 1-5 U=1.5 U= 5 1 15 t 3 4 5 U

Small-gap Mott insulator Pulse energy dependence of the relaxation rate Werner, Held & Eckstein (14) thermal value extrapolated value A( ).3.5..15.1 U= U=.5 U=3 U=3.5 U=4 U=4.5 D(t)/D(t=1) 1 U=.5.5-6 -4-4 6 =3 / =.5 / = / =1.5 / 1 3 4 5 6 t

Small-gap Mott insulator Pulse energy dependence of the relaxation rate Werner, Held & Eckstein (14) thermal value extrapolated value A( ).3.5..15.1 U= U=.5 U=3 U=3.5 U=4 U=4.5 D(t)/D(t=1) 1 U=3.5.5-6 -4-4 6 =3.5 / =3 / =.5 / = / 1 3 4 5 6 t Evidence for fast and slow relaxation time

Small-gap Mott insulator Pulse energy dependence of the relaxation rate Werner, Held & Eckstein (14) thermal value extrapolated value 5 U=4 U=3.5 U=3 U=.5 U= relaxation time 15 1 D(t)/D(t=1) 1 U=3.5 5.5 3 3.5 4 /( /) =3.5 / =3 / =.5 / = / 1 3 4 5 6 t Evidence for fast and slow relaxation time

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Fast doublon-hole production by the scattering process doublon high! doublon low + doublon low + hole low

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Fast doublon-hole production by the scattering process doublon high! doublon low + doublon low + hole low hole high! hole low + doublon low + hole low

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Fast doublon-hole production by the scattering process doublon high! doublon low + doublon low + hole low hole high! hole low + doublon low + hole low Consider only upper Hubbard band: doublon high! 3 doublon low

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Fast doublon-hole production by the scattering process doublon high! doublon low + doublon low + hole low hole high! hole low + doublon low + hole low Consider only upper Hubbard band: doublon high! 3 doublon low fast time-scale associated with these processes? Slow time scale related to multi-particle scattering processes

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Time evolution of the spectral function.1 U=3.5 t=18 t=4 t=3 t=36 =4 / t=4 I(, t) -4-3 - -1 1 3 4

Small-gap Mott insulator Impact ionization Werner, Held & Eckstein (14) Time evolution of the spectral function.1 t=3 t=36 t=4 I(, t)-i(, t=4) gain in low-energy weight =.3 x loss in high-energy weight.5 1 1.5.5 3 3.5 4 4.5

Simple model Impact ionization Werner, Held & Eckstein (14) High (low) energy population dd1 dt dd dt d dt D imp imp D 1 (D ) [D = D 1 + D ] = therm = 1 1 D1 dd1 = 3 dt imp D th D Two exponential relaxations D th D(t) = + D 1 (t s )e t ts + D th D(t s ) + D 1 (t s ) e t t s Obtain,, D 1 (t s ) by fitting

Simple model Impact ionization Werner, Held & Eckstein (14) High (low) energy population D 1 (D ) [D = D 1 + D ] D 1 (t s ) D(t s ) U 3.5.88 7. 18.8.5.5.67 7.75 19..5.44 9.35 19.6 3.5 3.46 13.4 6.3 3 3.4 15. 61.4.5 3.6 16.5 64.9 3.5 3.5.15 44. 376 3 3.5.83 48.4 57 4 (4.19 86.9 599 ) initial high-energy populations impact ionization thermalization very small D1: single-exponential relaxation

Simple model Impact ionization Werner, Held & Eckstein (14) High (low) energy population D 1 (D ) [D = D 1 + D ] D 1 (t s ) D(t s ) U 3.5.88 7. 18.8.5.5.67 7.75 19..5.44 9.35 19.6 3.5 3.46 13.4 6.3 3 3.4 15. 61.4.5 3.6 16.5 64.9 3.5 3.5.15 44. 376 3 3.5.83 48.4 57 4 (4.19 86.9 599 ) initial high-energy populations impact ionization thermalization

Simple model Impact ionization Werner, Held & Eckstein (14) Two-step relaxation predicted by the model. thermal normalized doublon population 1.5 1..5 small high energy population contributes significantly to doublon production U=3.5, =3.5 / DMFT data D 1 (t-t s )/D(t s ) D (t-t s )/D(t s ) D(t-t s )/D(t s ) D th /D(t s ) 1 3 4 5 t-t s

Simple model Impact ionization Werner, Held & Eckstein (14) Fluence dependence: impact ionization timescale shows little fluence dependence thermalization timescale shows larger fluence dependence amplitude D(t s ) D th th.5.18.36 4.884 19.9 14.5.49.917 4.593 19.5 194 1.167.334 3.879 18.3 147.593.15.854 15.6 85. 6.165.5 1.996 11. 46.1 amplitude > : doublon-conserving scattering processes start to deplete the high-energy population

Competing effects Impact ionization Werner, Held & Eckstein (14) Fluence dependence: increasing role of doublon-conserving scattering processes doublon high + doublon low! doublon intermediate 1 amplitude=5 amplitude= amplitude=.5 1 amplitude=5 amplitude= amplitude=.5 I(, t=36)-i(, t=4) [a. u.].5 -.5 U=3.5 I(, t=36)-i(, t=4) [a. u.].5 -.5 U=4-1 1 3 4 5-1 1 3 4 5

Competing effects Impact ionization Werner, Held & Eckstein (14) Scattering with external degrees of freedom (phonons, spins...): (dd 1 /dt) imp+ph/mag = ( 1/ 1/apple)D 1 (dd /dt) imp+ph/mag = (3/ +1/apple)D 1 Reduction in high-energy population decreases effect of impact ionization How effective is the cooling of photo-doped carriers by scattering with phonons/spins?

Cooling of carriers Effect of short-ranged antiferromagnetic correlations Eckstein & Werner (14) 4-site cluster calculations for D Hubbard give cooling rate =3 S NN NN spin correlations

Mobility of photo-doped carriers Mott insulating solar cells LaVO 3 on top of SrTiO 3 has suitable gap size Strong internal fields (carrier separation) Assmann, Held, Sangiovanni... (13) Strong correlations (impact ionization)

Mobility of photo-doped carriers Mott insulating solar cells LaVO 3 on top of SrTiO 3 has suitable gap size Nonequilibrium DMFT simulations show Assmann, Held, Sangiovanni... (13) Localization by strong internal fields a A(,z) [offset].7.6.5.4.3 z=6 z=5 z=4 b c 1.4 1.996 1.8 density ( E: magnetic for each hopping: energy gain Ea, but kinetic energy is bounded..1 E F z=3 z= z=1 - -1 1.6.4. T=1/8 1/5, 1 1 3

Mobility of photo-doped carriers Mott insulating solar cells LaVO 3 on top of SrTiO 3 has suitable gap size Nonequilibrium DMFT simulations show Assmann, Held, Sangiovanni... (13) Localization by strong internal fields Efficient separation of carriers in the presence of AFM order =3 vdrift a A(,z) [offset].7.6.5.4.3 z=6 z=5 z=4 b c 1.4 1.996 1.8 density ( E: magnetic..1 E F z=3 z= z=1 - -1 1.6.4. T=1/8 1/5, 1 1 3

Summary Relaxation of photo-doped carriers - some insights from DMFT Exponential scaling of thermalization time with gap size If gap < width of Hubbard bands: pulse-energy dependent initial relaxation due to impact ionization Decay of high-energy population due to scattering with spins 1/ m Drift velocity in polar heterostructures limited by scattering with spins v drift m Contribution of impact ionization to the efficiency of Mott solar cells requires more detailed analysis