Continuous-time quantum Monte Carlo algorithms for impurity problems
|
|
- Simon Pierce
- 5 years ago
- Views:
Transcription
1 Continuous-time quantum Monte Carlo algorithms for impurity problems Michel Ferrero Centre de Physique Théorique Ecole Polytechnique, France Quantum Monte Carlo methods at work for novel phases of matter Trieste, February 2-3, 2012
2 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
3 References Quantum impurity problems: A. C. Hewson, The Kondo Problem to Heavy Fermions, Cambridge University Press CT-QMC solvers: E. Gull et al., RMP (2011) The interaction-expansion algorithm Rubtsov et al., PRB (2005) and Rubtsov and Lichtenstein, JETP Lett. (2004) The hybridization-expansion algorithm Werner and Millis, PRB (2006) and Werner et al., PRL (2006) Dynamical mean-field theory: A. Georges et al., RMP (1996) Open source codes with CT-QMC solvers: TRIQS: ALPS:
4 Quantum impurity problems They generically describe the behavior of a magnetic impurity embedded in an electronic host The impurity is a set of orbitals carrying local many-body interactions and localized spins. It can exchange electrons with an uncorrelated fermionic bath e «Atomic» problem with several orbitals H loc e Electronic bath
5 The Kondo problem Experiments have shown that magnetic impurities in a metallic host (e.g. iron in gold) induce a minimum in the resistivity versus temperature curve The reason is that at low temperatures the conduction electrons screen the magnetic impurities: this is the Kondo effect
6 The Kondo effect: general idea At high temperature, conduction electrons pass by the impurity magnetic moment and only experience a small antiferromagnetic correlation impurity d y ¾ v F c y k¾ conduction electron At low temperature, the impurity magnetic moment and conduction electrons bind to form an overall non-magnetic state impurity d y ¾ c y k¾ conduction electron
7 The Anderson model A very successful model to understand magnetic impurities in a metallic host is the Anderson model Hamiltonian: ½(²) d y ¾ U ² 0 V k¾ Z = Tr exp( H) H = X c y k¾ ¾=";# ² 0 d y ¾d ¾ + Un d" n d# + X k;¾ V k¾ c y k¾ d ¾ + h:c: + X k;¾ ² ² k¾ c y k¾ c k¾ local many-body interaction Hybridization to the bath free bath states
8 Action for the Anderson model After integrating out the fermionic bath: Z Z = D[d y ; d]e S Z S = X d d 0 d y ¾( )G 1 0¾ ( 0 )d ¾ ( 0 ) + d Un d" ( )n d# ( ) ¾ 0 0 G 1 0¾ (i! n) = i! n ² 0 ¾ (i! n ) non-interacting Green s function ¾ (i! n ) = X k Remember: G 0¾ Generic case: jv k¾ j 2 i! n ² k¾ Z hybridization function: describes the transition between the bath and the orbital contains the information about the structure of the bath! S = X a;b Z 0 d d 0 d y a( )G 1 0;ab ( 0 )d b ( 0 ) + Z 0 d H int (fd y a; d a g)( ) G 1 0;ab (i! n) = i! n ± ab + t ab ab (i! n )
9 A difficult problem but a rich toolbox of methods The Anderson model is a many-body (correlated) problem with an infinite number of degrees of freedom! It has attracted a lot of interest and many techniques have been developed: Analytical methods (Bethe Ansatz, BCFT) Controllable algorithms: Exact diagonalization Numerical renormalization group Density matrix renormalization group Continuous-time quantum Monte Carlo algorithms Approximate solvers (e.g. NCA) All have pros and cons!
10 Impurity model in many contexts Quantum impurity problems appear in many different contexts magnetic impurity in metallic host nanostructures, quantum dots within the DMFT: pnictides, cuprates, cold-atomic systems,
11 A renewed interest within DMFT There are very efficient Monte Carlo algorithms for bosons For fermions, the situation is different and approximation methods are needed dynamical mean-field theory (DMFT) DMFT has brought a renewed interest in quantum impurity models in the late 90 s but also emphasized the limitations of existing algorithms and asked for the development of new tools
12 Dynamical mean-field theory The dynamical mean-field theory makes an approximation of a lattice model using an auxiliary quantum impurity problem It has been used a lot to find an approximate solution of the Hubbard model: H = X k ² k c y k c k + U X i n i" n i# free electrons on a lattice many-body interaction on every lattice site The DMFT has been used in many different contexts: realistic electronic structure calculations, models for superconductivity, cold-atomic systems, etc.
13 Dynamical mean-field theory cont d Lattice Hubbard model A. Georges and G. Kotliar, PRB (1992) A. Georges et al., RMP (1996) Anderson impurity model Bath G 0 Lattice Hubbard model is mapped onto a single-site model with an effective action: S = X ¾ Z 0 d d 0 d y ¾( )G 1 0¾ ( 0 )d ¾ ( 0 ) + Z 0 d Un d" ( )n d# ( ) The bath has to be set self-consistently
14 DMFT equations and limitations S e = Z 0 Z c y ¾( )G 1 0 ( 0 )c ¾ ( 0 ) + 0 (i! n ) = G 1 0 (i! n) G 1 (i! n ) lattice (i! n ; k) Ã (i! n ) Un " ( )n # ( ) G(i! n ) = X k G lattice (i! n ; k) G 1 0 (i! n) = ³ X k (i!n ) i! n + ¹ ² k (i! n ) In DMFT the self-energy of the lattice is local: Z, m*, coherence temperature, lifetimes are constant along the Fermi surface
15 Example: the Mott transition The DMFT is good at describing transfer of spectral weight due to correlation effects Example: the destruction of coherent quasiparticles (small QP coherence scale, short lifetime) near Mott transition Metal Mott transition U Insulator
16 What s the challenge? The impurity solver must compute the local Green s function Bath can have a rich structure, be gapped (insulators, superconductors) Structures appear at all scales (transfer of spectral weight in the Mott transition) The impurity solver must be able to treat many orbitals (e.g. pnictides) The interaction Hamiltonian must be generic (pair-hopping, spin flip terms) The model is studied in different temperature regimes One would like to be able to have real-frequency spectra
17 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
18 Continuous-time quantum Monte Carlo methods They have been a small revolution! There are different flavors: CT-INT: Interaction expansion CT-HYB: Hybridization expansion CT-AUX: Auxiliary-field formulation The underlying principle is the same for all these algorithms They stem from earlier work on diagrammatic Monte Carlo Prokof ev and Svistunov, PRL (1998) Prokof ev et al., Sov. Phys. JETP (1998)
19 Some notable results DMFT used together with electronic structure calculations Rotationally invariant multi-orbital problems Kondo lattice Cluster DMFT (e.g. for cuprate superconductors) Fermionic cold-atomic systems 4-point correlators for extensions of DMFT (dual-fermions) First steps in non-equilibrium
20 Diagrammatic Monte Carlo Write the partition function in a perturbation theory Z = Generate stochastically diagrams so that the statistical weight of a diagram is the absolute value of its contribution in the partition function expansion (Monte Carlo) This is possible even if the diagrams depend on continuous variables Use these diagrams to compute the partition function, Green s functions, etc.
21 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
22 Monte Carlo: a quick summary The Monte Carlo is a method to do sums: X X p(x)f(x); where p(x) > 0; p(x) = 1 x x The idea is to generate stochastically a large set of such that the probability to find a given is p(x n ) The original sum is replaced by an average over the set fx i g X p(x)f(x)» 1 N x NX f(x i ) Note that x should be understood as a set of variables x = ( ; ¾; : : :) One can think of x as a configuration in a phase space i=1 x n x i
23 Markov chain How do we generate configurations with the correct distribution? We don t sample independent configurations, but generate a Markov chain: x 1! x 2! x 3! x 4! : : : The transition probability W x;y to go from configuration x to a configuration y must satisfy: Normalization: X W x;y = 1 y Ergodicity: one must be able to reach any configuration Stationary distribution, balance condition: X d dt p(x) = 0 p(y)w y;x X p(x)w x;y = 0 y6=x y6=x Detailed balance (sufficient but not necessary): W x;y = p(y) W y;x p(x)
24 Metropolis algorithm How do we get a transition probability satisfying these criteria? Proposal rejection scheme: Propose a change with a chosen proposal rate Accept this proposal with a probability Otherwise don t change the configuration The total transition rate is Metropolis algorithm: A x;y W x;y = A x;y P x;y h A x;y = min 1; P i y;xp(y) P x;y p(x) P x;y It satisfies the detailed balance and therefore the Markov chain will be distributed according to p(x)
25 The fermionic sign problem Imagine we want to compute this average: P x hfi = P w(x)f(x) x w(x) We would like to use w(x) as a probability. But what if it can be negative? We must use the absolute value instead: hfi = P x jw(x)jf(x)sign(w(x)) Px jw(x)jsign(w(x))» P N i=1 f(x i)sign(w(x i )) P N i=1 sign(w(x i)) If signs alternate the denominator is very small and there is a big variance! Gets worse at low temperatures, big systems Fermionic problems almost always suffer this sign problem!
26 A Monte Carlo algorithm propose update x! y compute P x;y ; P y;x ; p(y) Metropolis h A x;y = min 1; P i y;xp(y) P x;y p(x) NO x! x initialize simulation leave configuration unchanged measures observables YES x! y update configuration
27 A couple of things to be careful about The measurements must start when the Markov chain has reached a stationary distribution. One must make sure to let the system thermalize! Successive points in the Markov chain are clearly correlated. These autocorrelation effects must be taken into account when computing the statistical errors: Binning analysis Jackknife Ideally one should measure only when configurations have become uncorrelated There are tools around to do these analysis!
28 Partition function: Diagrammatic Monte Carlo Z = Propose diagrams in a Markov chain: Z C w(c)» XMC C sign(w(c)) Accept these proposals with a rate (Metropolis) such that the diagrams appear with probability density w(c) From the generated configurations, compute the observables you are interested in hfi = 1 Z w(c)f(c)» 1 XMC f(c) sign(w(c)) Z C Z The different versions of continuous-time Monte Carlo solvers correspond to different choices of writing the partition function C
29 Getting a series expansion for the partition function Write the action in two parts S = S A + S B Express the partition function as Z Z Z = D[d y ; d]e S A e S B = D[d y ; d]e S A Using that hxi = 1 Z D[d y ; d]e S X Z X n ( 1) n n! S n B we get Z = Z A X n ( 1) n n! D E T SB n A This is an average over the states described by the action A. It generally involves sums and integrals over imaginary time and can have a diagrammatical representation.
30 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
31 Interaction-expansion CT-QMC We focus on the simplest Anderson model (can be generalized to density-density multi-orbital problems) We want to derive an expansion around the non-interacting limit (expansion in the interaction) ½(²) d y ¾ U ² 0 V k¾ c y k¾ H = X ¾=";# ² 0 d y ¾d ¾ + Un d" n d# + X k;¾ V k¾ c y k¾ d ¾ + h:c: + X k;¾ ² ² k¾ c y k¾ c k¾
32 Interaction expansion cont d We work in the imaginary-time formalism Z Z = D[d y ; d]e S hai = 1 Z D[d y ; d]e S A Z The action for the Anderson model: S = X ¾ Z 0 d d 0 d y ¾( )G 1 0¾ ( 0 )d ¾ ( 0 ) + Z 0 d Un d" ( )n d# ( ) The partition function can be written as Z = Z D[d y ; d]e S 0 ³ Z exp 0 d ( U)n d" ( )n d# ( ) this will produce a time-ordered average over the noninteracting state
33 Interaction expansion cont d We write a series expansion for the exponential ³ Z = Z 0 DT Z exp = Z 0 D T X = Z 0 X n n Z 0 0 ( U) n n! E d ( U)n d" ( )n d# ( ) 0 ³ Z ne d n d" ( )n d# ( ) 0 d 1 d n ( U) n n! 0 D E T n d" ( 1 )n d# ( 1 ) n d" ( n )n d# ( n ) 0 At this stage we have a perturbation expansion for the partition function of the type: Z XMC Z = w(c)» sign(w(c)) However there is an obvious sign problem! C C
34 Interaction expansion cont d The trick is to rewrite the interaction differently Un " n # = U X ³n " s" ³n # s# + U ³n " + n # 2 2 s=";# + const s¾ = (2± s¾ 1)± we absorb this term in the chemical potential We eventually get (up and down spins decouple) Z = Z 0 X n ( U) n n! Z d i 1 2 n X s i D E T (n d" ( 1 ) s1 ") (n d" ( n ) sn ") 0 D E T (n d# ( 1 ) s1 #) (n d# ( n ) sn #) By tuning ± we can improve the sign problem a lot! 0
35 Z = Z C w(c)» XMC C Interaction expansion cont d In the end we have X Z X X Z = Z 0 d 1 d n n s 1 s n ( U) n D E T n!2 n (n d" ( 1 ) s1 ") (n d" ( n ) sn ") Product of two non-interacting averages sign(w(c)) hfi = 1 Z Now we need to find a way to compute the averages Z Sum over many (some continuous) variables D T (n d# ( 1 ) s1 #) (n d# ( n ) sn #) C w(c)f(c)» 1 Z XMC C 0 E 0 f(c) sign(w(c))
36 Computing the averages The averages are on a non-interacting state. Therefore we can use Wick s theorem With this definition of the Green s function G 0¾ ( ) = T hc ¾ ( )c y ¾(0)i 0 we get D E (n d¾ ( 1 ) s1 ¾) (n d¾ ( n ) sn ¾) = det Dn ¾ = 0 G 0¾ (0 ) s1 ¾ G 0¾ ( 1 2 ) G 0¾ ( 2 1 ) G 0¾ (0 ) s2 ¾ G 0¾ ( 2 3 ) G 0¾ ( n n 1 ) G 0¾ (0 ) sn ¾ and X Z Z = Z 0 n X X d 1 d n s 1 s n ( U)n n!2 n det D" n det D # n
37 MC sum: Monte Carlo elements X Z X X d 1 d n n 1 > > n s 1 s n The configurations are diagrams of the perturbation expansion. They can be seen as a set of interaction vertices at different imaginary times with an auxiliary spin at every vertex. s i n The weight of every diagram is given by ³ U n det D " 2 n det D n #
38 Generating diagrams We need to create a Markov chain of diagrams We can propose any changes to go from one diagram to another. A simple solution is to use two moves : An insertion of a vertex: we pick a random imaginary time and insert a vertex with a spin randomly up or down (A) A removal of a vertex: pick a random vertex and remove it (B) (B) (A) 1 2 i n 0
39 What is the acceptance rate? Insertion of a vertex h A x;y = min 1; P i y;x½(y) P x;y ½(x) 1 2 n x; n 1 2 n n+1 0 y; n + 1 P x;y ½(x) = d n+1 P y;x ½(y) = n + 1 ³ U 2 Accept move with: ³ U 2 n det D " n det D # n ny d i i=1 n+1 n+1 det D " n+1 det Y D# n+1 d i i=1 h A x;y = min 1; U n + 1 det D" n+1 D# i n+1 det D nd " n #
40 What is the acceptance rate? Removal of a vertex h A x;y = min 1; P i y;x½(y) P x;y ½(x) 1 2 n x; n 1 n 1 0 y; n 1 P x;y ½(x) = n ³ U 2 P y;x ½(y) = d n ³ U 2 Accept move with: n det D " n det D # n ny d i i=1 n 1 det D " n 1 det D# n 1 n 1 Y i=1 d i d 0 i h A x;y = min 1; n U det D" n 1 D# i n 1 det D nd " n #
41 Getting the Green s function We know how to generate a distribution corresponding to the terms in the partition function Now we just need to find how to measure the Green s function from this distribution G( ) = T hd( )d y (0)i = 1 Z D[d y ; d]e S d y (0)d( ) Z S = X ¾ Z 0 Z d d 0 d y ¾( )G 1 0¾ ( 0 )d ¾ ( 0 ) + 0 d Un d" ( )n d# ( ) G 1 0¾ ( 0 ) = (@ 0 ² 0 )±( 0 ) ¾ ( 0 ) We finally see that G ¾ ( ) = 1 ± ln Z ± ¾ ( )
42 Getting the Green s function cont d We have Z = Z 0 Z C ³ U 2 n det D " n det D # n The functional derivative eventually gives: G ¾ (i! n ) = G 0¾ (i! n ) 1 G2 0¾(i! n ) Z X ³ [Dn ¾ ] 1 ij ei! n( i j ) U 2 C ij n det D " n det D # n So we see that we need to compute the following Monte Carlo average to get the Green s function G ¾ (i! n )» G 0¾ (i! n ) 1 Z G2 0¾ (i! n) XMC C X ij [D ¾ n ] 1 ij ei! n( i j ) sign(w(c))
43 Computational effort The effort comes from the calculation of the determinants and of the inverse matrix (needed for the Green s function measure) Z ³ U n Z = Z 0 det D " C 2 n det D n # G ¾ (i! n )» G 0¾ (i! n ) 1 Z G2 0¾ (i! n) XMC It would be very slow to calculate them from scratch at every move They can be updated quickly using the Sherman-Morrison formula The computational effort grows in O(n 3 ) C X ij [D ¾ n ] 1 ij ei! n( i j ) sign(w(c))
44 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
45 Hybridization-expansion CT-QMC We focus on the simplest Anderson model (can easily? be generalized to generic multi-orbital problems) We want to derive an expansion around the atomic limit (expansion in the hybridization) ½(²) d y ¾ U ² 0 V k¾ c y k¾ H = X ¾=";# ² 0 d y ¾d ¾ + Un d" n d# + X k;¾ V k¾ c y k¾ d ¾ + h:c: + X k;¾ ² ² k¾ c y k¾ c k¾
46 Hybridization expansion cont d We work in the imaginary-time formalism Z Z = D[d y ; d]e S hai = 1 Z D[d y ; d]e S A Z The standard action for the Anderson model: Z S = X d d 0 d y ¾( )G 1 0¾ ( 0 )d ¾ ( 0 ) + d Un d" ( )n d# ( ) ¾ 0 0 G 1 0¾ (i! n) = i! n ² 0 ¾ (i! n ) ¾ (i! n ) = X jv k¾ j 2 i! n ² k¾ k Rewrite the action as the sum of an atomic term containing all local terms U; ² 0 and a term with the hybridization to the bath: Z S = S loc + X ¾ Z 0 d d 0 d y ¾( ) ¾ ( 0 )d ¾ ( 0 ) = S loc + X ¾ S ¾ hyb
47 Z = Z = Z = Z Hybridization expansion cont d P Z Y h X 1 D[d y ; d]e S loc ¾ S¾ hyb = D[d y ; d]e S ( 1) n i loc S ¾ n n! hyb ¾ n=0 1X Y ( 1) DT n ¾ E n¾ ^S¾ n ¾! hyb loc n " ;n # =0 ¾ 1X Z d " 1 n " ;n # =0 0 Trace involving both spin up and down operators : : : d 0" n " Z Y ¾ 0 ( 1) n ¾ n ¾! d # 1 Yn ¾ i=1 : : : d 0# n # n Y " Tr he H loc T ¾ ( ¾ i 0¾ i ) i=1 n# d y " ( " i )d "( 0" i ) Y Sum over many (continuous) variables i=1 Product of hybridization functions d y # ( # i )d #( 0# i ) i Z = Z C w(c)» XMC C sign(w(c)) hfi = 1 Z Z C w(c)f(c)» 1 Z XMC C f(c) sign(w(c))
48 Hybridization expansion cont d MC sum: 1X Z n " ;n # =0 0 d " 1 : : : d 0" n " Z 0 d # 1 : : : d 0# n # Diagrams: Weight: ¾ =" ¾ =# n Y " Tr he H loc T " 1 0" 1 " 2 0" 2 i=1 n# d y " ( " i )d "( 0" i ) Y Y ¾ # 1 0# 1 ( 1) n ¾ Unfortunately these diagrams have alternating signs problems! n ¾! n ¾ Y i=1 i=1 d y # ( # i )d #( 0# i ) i ¾ ( ¾ i 0¾ i ) 0
49 Trick: resumming diagrams The idea is to resum diagrams into a determinant. We start from a diagram where " 1 > : : : > n " " and sum all the permutations of f " i g ¾ =" " 1 0" 1 " 2 0" 1 + " 2 0" 2 " 1 0" 2 " 1 > " 2 " 1 $ " 2 h i Tr e H loc d y " ( " 1 )d "( 0" 1 )dy " ( " 2 )d "( 0" 2 ) h i ( 1)Tr e H loc d y " ( " 1 )d "( 0" 1 )dy " ( " 2 )d "( 0" 2 ) "( " 1 0" 1 ) "( " 2 0" 2 ) 1 2 "( " 2 0" 1 ) "( " 1 0" 2 ) = h i Tr e H loc d y " ( " 1 )d "( 0" 1 )dy " ( " 2 )d "( 0" 2 ) 1 det 2 "( " 1 k;l 2 k 0" l )
50 MC sum: Hybridization expansion Monte Carlo 1X n " ;n # =0 Z d " " 1 >:::> " 1 n 0" 0" 1 >:::> n : : : d 0" n Z d # # 1 >:::> # 1 n 0# 0# 1 >:::> n : : : d 0# n # Diagrams: ¾ =" " 1 0" 1 " 2 0" 2 ¾ =# # 1 0# 1 0 Weight: n Y " Tr he H loc T i=1 n# d y " ( " i )d "( 0" i ) Y i=1 d y # ( # i )d #( 0# i ) i ( 1) n "+n # det " ( " k 0" l ) det # ( # k 0# l ) 1 k;l n " 1 k;l n #
51 Generating diagrams New diagrams are generated with two moves : Insertion of an (anti)-link: chose a spin flavor and pick two random imaginary times such that there is no operator between them. Either construct a link (A) or an anti-link (B) Removal of a link: chose a spin flavor and remove a random link (C) " 1 0" 1 " 2 0" 2 (A) ¾ =" (C) " 1 0" 1 " 2 0" 2 " 3 0" 3 (B) " 1 0" 1 " 2 0" 2 " 1 0" 1 " 2 0" 2 " 3 0" 3 0
52 Insertion of an (anti)-link What is the acceptance rate for this move? h A x;y = min 1; P i y;x½(y) P x;y ½(x) ¾ =" " 1 0" 1 " 2 0" 2 " 1 0" 1 " 2 0" 2 " 3 0" 3 x; n y; n + 1 P x;y ½(x) = 1 2 d n+1 d 0 n+1 l max TrC x det Cx P y;x ½(y) = n + 1 TrC y det Cy n+1 Y i=1 ny d i d i 0 i=1 d i d 0 i 0 Accept with probability: h A x;y = min 1; l max n + 1 TrC y det i Cy TrC x det Cx
53 Removal of a link What is the acceptance rate for this move? h A x;y = min 1; P i y;x½(y) P x;y ½(x) ¾ =" " 1 0" 1 " 2 0" 2 " 1 0" 1 n Y 0 x; n y; n 1 P x;y ½(x) = n TrC x det Cx P y;x ½(y) = 1 2 d n i=1 d i d 0 i d 0 n l max TrC y det Cy n 1 Y i=1 d i d 0 i Accept with probability: h A x;y = min 1; n TrC y det i Cy l max TrC x det Cx
54 Measuring the Green s function We know how to sample diagrams with weights corresponding to their contribution in the partition function. Z = Z C How do we get the Green s function? G " ( ) = 1 Z G " ( ) = 1 Z ( 1) n "+n # det "C det #C TrC = Z Z C C G ¾ ( ) = 1 ± log Z ± ¾ ( ) ± det "C ± " ( ) ( 1)n "+n # det #C TrC X ±( " k 0" l + )[ 1 "C ] k;l ( 1) n "+n # det "C det #C TrC k;l Z C w(c)» XMC C sign(w(c)) w(c) Measure: G ¾ ( )» 1 Z XMC C X k;l ±( ¾ k 0¾ l + ) [ 1 ¾C ] k;l sign(w(c))
55 Measuring the Green s function cont d Each configuration give contributions for a discrete set of imaginary times: G ¾ ( )» 1 XMC X ±( k ¾ l 0¾ + ) [ 1 ¾C Z ] k;l sign(w(c)) C k;l ¾ =" " 1 0" 1 " 2 0" 2 " 2 0" 2 0 These contribution can be binned on a very fine imaginary-time grid. This induces high frequency noise in Matsubara frequencies
56 Measuring using Legendre polynomials Legendre polynomials are a basis to express function defined over an interval L. Boehnke et al., PRB (2011) We can express the imaginarytime Green s function in this basis p 2l + 1 G( ) = X l 0 P l [x( )]G l The coefficients G l in this basis decay very quickly
57 Legendre basis acting as a noise filter The noise in the Matsubara frequencies can be reduced by truncating the Legendre coefficients that are zero within their error bars A typical outcome of this procedure:
58 Computational effort Can the contribution of a diagram be computed quickly? Determinants can be updated quickly (Sherman-Morrison) For simple Hamiltonians, the trace is very easy ¾ =" " 1 0" 1 " 2 0" 2 > l " > 0 ¾ =# # 1 0# 1 0 > l # > 0 overlap: l 1 l 2 Tr C = e (l "+l # )¹ U(l 1 +l 2 ) Computational effort grows in O(n 3 )
59 What about non density-density Hamiltonians? The hybridization expansion algorithm can be modified for generic Hamiltonians Configurations are a set of creation / destruction operators of different flavor on a single imaginary-time line " 1 0" 1 " 2 " 0" 3 0" 2 3 0" 4 " 4 0 The main drawback is that there is no longer a quick way to compute the trace Operators are matrices that must be multiplied and traced over all atomic states The number of these atomic states quickly becomes huge with several orbitals
60 CT-INT versus CT-HYB CT-INT & CT-AUX: series in the interaction Many orbitals, weak coupling, high temperatures density-density Hamiltonians Average perturbation order» U CT-HYB: series in the hybridization function Good at low temperatures, strong coupling Can treat generic Hamiltonians Hard to treat many orbitals Average perturbation order is the kinetic energy
61 Pros and cons of the CT-QMC algorithms Pros: Faster than before Monte Carlo can easily be parallelized Flexible Hamiltonians (not CT-INT) Good scaling with number of orbitals if density-density Cons: Many orbitals difficult with generic Hamiltonian They are mainly in imaginary time, so one needs to do analytical continuation, and this is a very delicate procedure! Note: real-time algorithms have been developed. However the sign problem is a lot more present there. Sign problem
62 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
63 What did the CT-QMC algorithms change? In simple model studies, CT-QMC algorithms allow to: Reach much lower temperatures Acquire high-quality imaginary-time data and do more reliable analytical continuations In the study of realistic materials, they allow to treat the full SU(2) exchange in multi-orbital problems (e.g. pnictides, nickelates, ruthenates, ) More involved models can be studied Cluster extensions of DMFT Dual-fermion approaches Holstein-Hubbard model, out-of-equilibrium
64 State of the art: NRG versus CT-QMC Direct comparison between NRG and CT-QMC for a half-filled Anderson impurity model (flat Wilson bath) Parameters are: W = 1; U = 0:5; U=¼ = 2:5 Analytical continuation of the CT-QMC selfenergy using Padé approximants NRG: Rok Zitko CT-QMC: Jernej Mravlje
65 Detailed structures of self-energy Even fine details of the self-energy are captured correctly! Single-band Bethe lattice DMFT calculation show very good agreement, where only small temperature effect induce a difference between NRG and CT-QMC Im (!)!
66 Mott transition in the Hubbard model Single-site DMFT studies have shown that the transition is first order at finite temperature, but second order at zero temperature DMFT cluster extensions have been considered (more later) Werner et al., PRB (2007)
67 Not only single-site DMFT for the Hubbard model Holstein-Hubbard model H = X ² k c y k¾ c k¾ + X Un i" n i# k¾ i + X (b y i + b X i)(n i 1) +! 0 b y i b i Kondo lattice i i H = X k¾ ² k c y k¾ c k¾ + J X i ~S i ~¾ i Extensions of DMFT to include spatial correlations Cluster extensions that ask for multi-impurity models Dual-fermion approaches that require for higher correlators
68 Cuprate superconductors A(k i ;!) ARPES A(k i ;!) Kaminski et al., PRB (2005) A(k;! = 0) (0,0) (, ) Spectral properties are isotropic along the Fermi surface! (, )! Fermi arcs Strong variations over the Brillouin zone! (0,0)
69 k-space variations are a problem for single-site DMFT S e = Z 0 Z c y ¾( )G 1 0 ( 0 )c ¾ ( 0 ) + 0 (i! n ) = G 1 0 (i! n) G 1 (i! n ) lattice (i! n ; k) Ã (i! n ) Un " ( )n # ( ) G(i! n ) = X k G lattice (i! n ; k) G 1 0 (i! n) = ³ X k (i!n ) i! n + ¹ ² k (i! n ) In DMFT the self-energy of the lattice is local: Z, m*, coherence temperature, lifetimes are constant along the Fermi surface
70 Including a k-dependence Dynamical cluster approximation (DCA) M. H. Hettler et al., PRB (1998) T. Maier et al., RMP (2005) Bath G 0 Anderson multiimpurity model Cluster self-energies for every cluster momenta: (K) There are ways to improve DMFT, but then you need to solve a multi-impurity Anderson model!
71 k-dependence of the self-energy The DCA approach has been tested against a direct Monte Carlo evaluation of the Hubbard lattice problem The solution scales nicely with the number of sites and has converged for 32 sites Kozik et al., EPL (2010)
72 Back to cuprates with a 2-site extension M. Ferrero et al., PRB (2009) 1 2 (, ) electron exchange (0,0) 0% pseudogap in antinodal orbital metallic, protected nodal orbital Singlet state dominates Ordinary Fermi liquid metal 16% Hole doping Selective Mott transition
73 Can one describe the phenomenology of the cuprates? ARPES intensity maps (, ) 8% 12% 25% (0,0) Tunneling conductance (STM) 8%
74 Optical conductivity Interplane charge dynamics Resistivity underdoped overdoped
75 Confirmation of the results on larger clusters E. Gull et al., PRB (2010) Computations on clusters up to 16 sites have been performed Results show that the main features of the above results are robust The approach to the Mott insulator is characterized by successive Mott transitions in reciprocal space The regions close to the antinode get insulating first, while regions close to the node remain metallic down to the Mott insulating state
76 Dual-fermion calculations for the Hubbard model Dual-fermion approaches are another way to include spatial correlations. They require the knowledge of 4-point vertex functions These methods have benefitted a lot from the new CT-QMC Spectral intensity map of the Hubbard model at 14% doping Dynamical susceptibility at half-filling with a visible magnon spectrum in the paramagnetic state Rubtsov et al., PRB, (2009) Hafermann et al., PRL, (2009)
77 Stabilization of a Sarma phase in cold-atomic systems Dao et al., PRL (2008) Cold-atomic systems allow to study superfluid phases with great precision The situation is rather clear when there is an attraction between two balanced populations (BCS BEC crossover) But what happens as we introduce a population imbalance and the Fermi surface mismatch? Could we obtain a Sarma phase which is both superfluid and polarized? (¼; ¼) Sarma has shown that it is not stable at weak coupling (0; 0) But what about strong coupling?
78 A study of a negative-u Hubbard model shows that: At weak coupling a phase separated state is stabilized at low temperature At larger coupling a polarized superfluid is stabilized down to very low temperatures The presence of two Fermi surfaces could be observed in cold-atomic experiments Sarma phase cont d
79 Electronic structure calculations There has been a large effort to treat fully many-body effects within density-functional-based (DFT) electronic structure calculations One of the approaches is DFT+DMFT where the many-body effects are computed by DMFT electron exchange Unit cell hybridized to a bath Because realistic systems often involve many bands and because it is sometimes important to consider a full SU(2) Hamiltonian, new solvers have been extremely useful
80 As-p, O-p Fe-d La Correlations in pnictides superconductors LaFeAsO M. Aichhorn et al., PRB (2009) - moderate correlations - m* lifetime effects - shifted bands A(k;!) LDA bands Density of states Correlation effects seen by LDA+DMFT
81 Interaction quench in the Hubbard model Eckstein et al., PRB (2010) What happens it we suddenly turn on an interaction in the Hubbard model? Does it reach a steady state? Is this steady state a thermal state? A real-time CT-QMC algorithm has been used to study this problem within non-equilibrium DMFT 0 < U < 3:2V 3:2V < U
82 Summary Continuous-time quantum Monte Carlo algorithms have allowed for progress in computing the properties of strongly-correlated systems Lower temperatures Generic Hamiltonians, new models Larger number of orbitals / sites The idea of the algorithms is to sample stochastically the diagrams of a series expansion of the partition function According to one s need, different expansions can be used There are still limitations (sign problem, speed, ) and more work has to be done!
83 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
84 A need for tools to do many-body calculations Solving a quantum impurity problem is usually just one step of a larger calculation For example, a typical DMFT calculation might require to Manipulate Green s functions Solve a quantum impurity problem Sum quantities over the Brillouin zone of a lattice Analyze, archive, plot data One might also have to use codes from different origins and glue results together It is therefore useful to have a set of flexible tools to work with the objects of a many-body calculations
85 Illustration: a DFT+DMFT calculation Wien2k Converged? ½ LDA+DMFT! V KS New density Nºº k = X G 0 ºº 0(k; i! n )e i! n0 + n ½ DMFT (r) = X hrjª kº inºº k hª 0 kº 0jri k;ºº 0 DFT Step V KS! ² kº ; ª kº (r) Projectors Wien2k dmftproj P mº (k) = h! km jª kº i Gmm 0 1 = ³ 0 mm 0 + G loc 1 mm 0 Green s functions manipulations Lattice Green s function G 1 ºº 0 (k; i! n ) = (i! n + ¹ ² kº )± ºº 0 ºº 0(k; i! n ) G loc mm 0 (i! n ) = X kºº 0 P mº (k)g ºº 0(k; i! n )P º 0 m 0 (k) Anderson Impurity Model G 0 mm 0 (i! n )! mm 0(i! n ) CT-QMC
86 Toolbox for Research on Interacting Quantum Systems Was born from the need to have a versatile environment to study interacting quantum systems in different contexts Started at Ecole Polytechnique and IPhT Saclay (with O. Parcollet, M. Aichhorn, L. Pourovskii, V. Vildosola, C. Martins, L. Boehnke, ) It is an open-source (GPL license) library
87 What is in the package? Python libraries: Lattice tools Data analysis tools Green s functions Plotting interface Quantum impurity solvers Ready-to-use applications: Wien2TRIQS C++ libraries for developers: Monte Carlo Array library Determinant manipulations
88 Example: Bethe lattice DMFT
89 Example: Bethe lattice DMFT
90 Making codes open source Probably the best way to learn about Monte Carlo methods is to implement one s own algorithm In the hands-on we try to code an interaction-expansion continuous-time algorithm Note however that in general it is not always a good idea to reinvent the wheel! There are good open-source tools around and one should never hesitate to use them Making one s codes available for the rest of the scientific community is a very good idea!
91 Outline Quantum impurity problems Continuous-time quantum Monte Carlo (CT-QMC) methods Generic introduction to Monte Carlo The interaction-expansion algorithm (CT-INT) The hybridization-expansion algorithm (CT-HYB) Some applications of CT-QMC solvers TRIQS: a Toolbox for Research on Quantum Interacting Systems Implementation details for the CT-INT
92 Reminder: A Monte Carlo algorithm propose update x! y compute P x;y ; P y;x ; p(y) Metropolis h A x;y = min 1; P i y;xp(y) P x;y p(x) NO x! x initialize simulation leave configuration unchanged measures observables YES x! y update configuration
93 Implementing the CT-INT algorithm in practice Monte Carlo loop Configuration 1 2 n 0 P y;x p(y) P x;y p(x) Moves acting on the configuration (insert/remove) Measures on the configuration Z; G Results
94 Implementation cont d Good practice to try to separate the algorithm in pieces The Monte Carlo loop does not need to know about the system. It just calls moves and expects back a Metropolis ratio to decide if the move is accepted or rejected What we need to do is implement insertions and removals in the configuration and compute the ratio 1 2 n 0 P y;x p(y) P x;y p(x) This ratio involves a determinant! This is the tricky part!
95 Reminder The weight of a Monte Carlo configuration is p(c n ) = ³ U 2 n det D " n det D # n And the matrix in the determinant has this definition det Dn ¾ = G 0¾ (0 ) s1 ¾ G 0¾ ( 1 2 ) G 0¾ ( 2 1 ) G 0¾ (0 ) s2 ¾ G 0¾ ( 2 3 ) G 0¾ ( n n 1 ) G 0¾ (0 ) sn ¾ with s¾ = (2± s¾ 1)±
96 Inserting a new vertex Inserting a vertex adds a line and column to the matrix h A x;y = min 1; U 1 2 i n + 1 det D" n+1 D# i n+1 n det D nd " n # 0 det D ¾ n = G 0¾ (0 ) s1 ¾ G 0¾ ( 1 2 ) G 0¾ ( 1 3 ) G 0¾ ( 2 1 ) G 0¾ (0 ) s2 ¾ G 0¾ ( 2 3 ) G 0¾ ( 3 1 ) G 0¾ ( 3 2 ) G 0¾ (0 ) s3 ¾ det Dn+1 ¾ = G 0¾ (0 ) s1 ¾ G 0¾ ( 1 2 ) G 0¾ ( 1 3 ) G 0¾ ( 1 4 ) G 0¾ ( 2 1 ) G 0¾ (0 ) s2 ¾ G 0¾ ( 2 3 ) G 0¾ ( 2 4 ) G 0¾ ( 3 1 ) G 0¾ ( 3 2 ) G 0¾ (0 ) s3 ¾ G 0¾ ( 3 4 ) G 0¾ ( 4 1 ) G 0¾ ( 4 2 ) G 0¾ ( 4 3 ) G 0¾ (0 ) s4 ¾
97 Removing a vertex Removing a vertex removes a line and column in the matrix n G 0¾ (0 ) s1 ¾ G 0¾ ( 1 2 ) G 0¾ ( 1 3 ) G 0¾ ( 1 4 ) G 0¾ ( 2 1 ) G 0¾ (0 ) s2 ¾ G 0¾ ( 2 3 ) G 0¾ ( 2 4 ) G 0¾ ( 3 1 ) G 0¾ ( 3 2 ) G 0¾ (0 ) s3 ¾ G 0¾ ( 3 4 ) G 0¾ ( 4 1 ) G 0¾ ( 4 2 ) G 0¾ ( 4 3 ) G 0¾ (0 ) s4 ¾ i h A x;y = min 1; n U det D" n 1 D# i n 1 det D nd " n #
98 No need to recompute a full determinant We see that we need to compute ratios of determinant of matrices that have an added / removed line and column This can be done faster that by computing the full determinant of the matrix before and after the move The trick is to keep a copy of both the matrix and its inverse and to update them Indeed: D 1 ij = ( 1) i+j 1 det D det D[j; i] Means the matrix with line j and column i removed
99 Fast updates for determinant and inverse matrix This is a special case of the Sherman-Morrison formula µ A B C d µ U V W x = µ AU + BW AV + xb CU + dw CV + xd = µ Insertion Removal V = xa 1 B x 1 = d CA 1 B µ A B det C d det A = x 1 A 1 = U x 1 V W W = xca 1 U = A 1 + x 1 V W det A µ = x A B det C d
100 Fast updates cont d One can find the determinant ratios and the new inverse matrix in operations of order n 2 Here we have shown insertions and removals happening at the very end of the matrix. If the insertions or removals are somewhere in the middle, one must add a correction ( 1) i+j An important thing to note is that the determinant ratio can be obtained before actually computing the inverse matrix. We can therefore decide to accept or reject the move and update the inverse matrix only if the move is accepted! Another good thing about have the inverse matrix is that we need it anyway to compute the Green s function These are all the ingredients to write an efficient algorithm!
An introduction to the dynamical mean-field theory. L. V. Pourovskii
An introduction to the dynamical mean-field theory L. V. Pourovskii Nordita school on Photon-Matter interaction, Stockholm, 06.10.2016 OUTLINE The standard density-functional-theory (DFT) framework An
More informationDiagrammatic Monte Carlo methods for Fermions
Diagrammatic Monte Carlo methods for Fermions Philipp Werner Department of Physics, Columbia University PRL 97, 7645 (26) PRB 74, 15517 (26) PRB 75, 8518 (27) PRB 76, 235123 (27) PRL 99, 12645 (27) PRL
More informationDiagrammatic Monte Carlo simulation of quantum impurity models
Diagrammatic Monte Carlo simulation of quantum impurity models Philipp Werner ETH Zurich IPAM, UCLA, Jan. 2009 Outline Continuous-time auxiliary field method (CT-AUX) Weak coupling expansion and auxiliary
More informationDynamical mean field approach to correlated lattice systems in and out of equilibrium
Dynamical mean field approach to correlated lattice systems in and out of equilibrium Philipp Werner University of Fribourg, Switzerland Kyoto, December 2013 Overview Dynamical mean field approximation
More informationIntroduction to DMFT
Introduction to DMFT Lecture 2 : DMFT formalism 1 Toulouse, May 25th 2007 O. Parcollet 1. Derivation of the DMFT equations 2. Impurity solvers. 1 Derivation of DMFT equations 2 Cavity method. Large dimension
More informationw2dynamics : operation and applications
w2dynamics : operation and applications Giorgio Sangiovanni ERC Kick-off Meeting, 2.9.2013 Hackers Nico Parragh (Uni Wü) Markus Wallerberger (TU) Patrik Gunacker (TU) Andreas Hausoel (Uni Wü) A solver
More informationDe l atome au. supraconducteur à haute température critique. O. Parcollet Institut de Physique Théorique CEA-Saclay, France
De l atome au 1 supraconducteur à haute température critique O. Parcollet Institut de Physique Théorique CEA-Saclay, France Quantum liquids Quantum many-body systems, fermions (or bosons), with interactions,
More informationElectronic correlations in models and materials. Jan Kuneš
Electronic correlations in models and materials Jan Kuneš Outline Dynamical-mean field theory Implementation (impurity problem) Single-band Hubbard model MnO under pressure moment collapse metal-insulator
More informationAdvanced Computation for Complex Materials
Advanced Computation for Complex Materials Computational Progress is brainpower limited, not machine limited Algorithms Physics Major progress in algorithms Quantum Monte Carlo Density Matrix Renormalization
More information6. Auxiliary field continuous time quantum Monte Carlo
6. Auxiliary field continuous time quantum Monte Carlo The purpose of the auxiliary field continuous time quantum Monte Carlo method 1 is to calculate the full Greens function of the Anderson impurity
More informationDual fermion approach to unconventional superconductivity and spin/charge density wave
June 24, 2014 (ISSP workshop) Dual fermion approach to unconventional superconductivity and spin/charge density wave Junya Otsuki (Tohoku U, Sendai) in collaboration with H. Hafermann (CEA Gif-sur-Yvette,
More informationThe Hubbard model out of equilibrium - Insights from DMFT -
The Hubbard model out of equilibrium - Insights from DMFT - t U Philipp Werner University of Fribourg, Switzerland KITP, October 212 The Hubbard model out of equilibrium - Insights from DMFT - In collaboration
More informationMagnetic Moment Collapse drives Mott transition in MnO
Magnetic Moment Collapse drives Mott transition in MnO J. Kuneš Institute of Physics, Uni. Augsburg in collaboration with: V. I. Anisimov, A. V. Lukoyanov, W. E. Pickett, R. T. Scalettar, D. Vollhardt,
More informationQuasiparticle dynamics and interactions in non uniformly polarizable solids
Quasiparticle dynamics and interactions in non uniformly polarizable solids Mona Berciu University of British Columbia à beautiful physics that George Sawatzky has been pursuing for a long time à today,
More informationBad metallic behavior, entropy & thermopower: DMFT results on ruthenates
Bad metallic behavior, entropy & thermopower: DMFT results on ruthenates Jernej Mravlje SPICE Workshop, Schloss Waldthausen, 28. June 2015 People A. Georges, CdF,EP,UNIGE & X. Deng, Rutgers,NJ, L. de'medici,
More informationNUMERICAL METHODS FOR QUANTUM IMPURITY MODELS
NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Anrew Mitchell Utrecht University Quantum impurity problems Part 1: Quantum impurity problems
More informationO. Parcollet CEA-Saclay FRANCE
Cluster Dynamical Mean Field Analysis of the Mott transition O. Parcollet CEA-Saclay FRANCE Dynamical Breakup of the Fermi Surface in a doped Mott Insulator M. Civelli, M. Capone, S. S. Kancharla, O.P.,
More informationSpin and orbital freezing in unconventional superconductors
Spin and orbital freezing in unconventional superconductors Philipp Werner University of Fribourg Kyoto, November 2017 Spin and orbital freezing in unconventional superconductors In collaboration with:
More informationAccelerated Monte Carlo simulations with restricted Boltzmann machines
Accelerated simulations with restricted Boltzmann machines Li Huang 1 Lei Wang 2 1 China Academy of Engineering Physics 2 Institute of Physics, Chinese Academy of Sciences July 6, 2017 Machine Learning
More informationTopological order in the pseudogap metal
HARVARD Topological order in the pseudogap metal High Temperature Superconductivity Unifying Themes in Diverse Materials 2018 Aspen Winter Conference Aspen Center for Physics Subir Sachdev January 16,
More informationHow to model holes doped into a cuprate layer
How to model holes doped into a cuprate layer Mona Berciu University of British Columbia With: George Sawatzky and Bayo Lau Hadi Ebrahimnejad, Mirko Moller, and Clemens Adolphs Stewart Blusson Institute
More informationCluster Extensions to the Dynamical Mean-Field Theory
Thomas Pruschke Institut für Theoretische Physik Universität Göttingen Cluster Extensions to the Dynamical Mean-Field Theory 1. Why cluster methods? Thomas Pruschke Institut für Theoretische Physik Universität
More informationInteraction-induced Symmetry Protected Topological Phase in Harper-Hofstadter models
Interaction-induced Symmetry Protected Topological Phase in Harper-Hofstadter models arxiv:1609.03760 Lode Pollet Dario Hügel Hugo Strand, Philipp Werner (Uni Fribourg) Algorithmic developments diagrammatic
More informationPseudogap opening and formation of Fermi arcs as an orbital-selective Mott transition in momentum space
PHYSICAL REVIEW B 8, 645 9 Pseudogap opening and formation of Fermi arcs as an orbital-selective Mott transition in momentum space Michel Ferrero,, Pablo S. Cornaglia,,3 Lorenzo De Leo, Olivier Parcollet,
More informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
More informationQuantum Cluster Methods (CPT/CDMFT)
Quantum Cluster Methods (CPT/CDMFT) David Sénéchal Département de physique Université de Sherbrooke Sherbrooke (Québec) Canada Autumn School on Correlated Electrons Forschungszentrum Jülich, Sept. 24,
More informationQuantum impurities in a bosonic bath
Ralf Bulla Institut für Theoretische Physik Universität zu Köln 27.11.2008 contents introduction quantum impurity systems numerical renormalization group bosonic NRG spin-boson model bosonic single-impurity
More informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationFrom Materials to Models and Back. Dieter Vollhardt
From Materials to Models and Back Dieter Vollhardt 28 th Edgar Lüscher Seminar, Klosters; February 8, 2017 From Materials to Models and Back - The Need for Models in Condensed Matter Physics - Outline:
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationNano-DMFT : the electronic structure of small, strongly correlated, systems
Nano-DMFT : the electronic structure of small, strongly correlated, systems Nanoscale Dynamical Mean-Field Theory for Molecules and Mesoscopic Devices in the Strong-Correlation Regime Author: S. Florens,
More informationQuantum impurity models Algorithms and applications
1 Quantum impurity models Algorithms and applications Collège de France, 5 Mai 21 O. Parcollet Institut de Physique Théorique CEA-Saclay, France Motivations : why do we need specific algorithms? A few
More informationIMPACT ionization and thermalization in photo-doped Mott insulators
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:
More informationThe Hubbard model in cold atoms and in the high-tc cuprates
The Hubbard model in cold atoms and in the high-tc cuprates Daniel E. Sheehy Aspen, June 2009 Sheehy@LSU.EDU What are the key outstanding problems from condensed matter physics which ultracold atoms and
More informationLocal moment approach to the multi - orbital single impurity Anderson and Hubbard models
Local moment approach to the multi - orbital single impurity Anderson and Hubbard models Anna Kauch Institute of Theoretical Physics Warsaw University PIPT/Les Houches Summer School on Quantum Magnetism
More informationA DCA Study of the High Energy Kink Structure in the Hubbard Model Spectra
A DCA Study of the High Energy Kink Structure in the Hubbard Model Spectra M. Jarrell, A. Macridin, Th. Maier, D.J. Scalapino Thanks to T. Devereaux, A. Lanzara, W. Meevasana, B. Moritz, G. A. Sawatzky,
More informationSeconde partie: Quelques questions liées au transport dans les matériaux à fortes corrélations électroniques
``Enseigner la recherche en train de se faire Chaire de Physique de la Matière Condensée Seconde partie: Quelques questions liées au transport dans les matériaux à fortes corrélations électroniques Antoine
More informationContinuous time QMC methods
Continuous time QMC methods Matthias Troyer (ETH Zürich) Philipp Werner (Columbia ETHZ) Emanuel Gull (ETHZ Columbia) Andy J. Millis (Columbia) Olivier Parcollet (Paris) Sebastian Fuchs, Thomas Pruschke
More informationAn efficient impurity-solver for the dynamical mean field theory algorithm
Papers in Physics, vol. 9, art. 95 (217) www.papersinphysics.org Received: 31 March 217, Accepted: 6 June 217 Edited by: D. Domínguez Reviewed by: A. Feiguin, Northeastern University, Boston, United States.
More informationDT I JAN S S"= = 11111'11 I HtI IBlIIIt g ~II. Report: ONR Grant N J Unclassified: For General Distribution LECTF
' DT I,, Final Report: ONR Grant N00014-91-J-1143 Unclassified: For General Distribution LECTF JAN 2 8 1993 Steven R. White, UC Irvine c Summary Over the last two years, we have used several different
More informationContinuous Time Monte Carlo methods for fermions
Continuous Time Monte Carlo methods for fermions Alexander Lichtenstein University of Hamburg In collaboration with A. Rubtsov (Moscow University) P. Werner (ETH Zurich) Outline Calculation of Path Integral
More informationFROM NODAL LIQUID TO NODAL INSULATOR
FROM NODAL LIQUID TO NODAL INSULATOR Collaborators: Urs Ledermann and Maurice Rice John Hopkinson (Toronto) GORDON, 2004, Oxford Doped Mott insulator? Mott physics: U Antiferro fluctuations: J SC fluctuations
More informationPG5295 Muitos Corpos 1 Electronic Transport in Quantum dots 2 Kondo effect: Intro/theory. 3 Kondo effect in nanostructures
PG5295 Muitos Corpos 1 Electronic Transport in Quantum dots 2 Kondo effect: Intro/theory. 3 Kondo effect in nanostructures Prof. Luis Gregório Dias DFMT PG5295 Muitos Corpos 1 Electronic Transport in Quantum
More informationARPES studies of cuprates. Inna Vishik Physics 250 (Special topics: spectroscopies of quantum materials) UC Davis, Fall 2016
ARPES studies of cuprates Inna Vishik Physics 250 (Special topics: spectroscopies of quantum materials) UC Davis, Fall 2016 Goals of lecture Understand why gaps are important and various ways that gap
More informationɛ(k) = h2 k 2 2m, k F = (3π 2 n) 1/3
4D-XY Quantum Criticality in Underdoped High-T c cuprates M. Franz University of British Columbia franz@physics.ubc.ca February 22, 2005 In collaboration with: A.P. Iyengar (theory) D.P. Broun, D.A. Bonn
More informationStrongly correlated Cooper pair insulators and superfluids
Strongly correlated Cooper pair insulators and superfluids Predrag Nikolić George Mason University Acknowledgments Collaborators Subir Sachdev Eun-Gook Moon Anton Burkov Arun Paramekanti Affiliations and
More informationTopological Kondo Insulator SmB 6. Tetsuya Takimoto
Topological Kondo Insulator SmB 6 J. Phys. Soc. Jpn. 80 123720, (2011). Tetsuya Takimoto Department of Physics, Hanyang University Collaborator: Ki-Hoon Lee (POSTECH) Content 1. Introduction of SmB 6 in-gap
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 21 Mar 2006
Non-Fermi-liquid phases in the two-band Hubbard model: Finite-temperature exact diagonalization study of Hund s rule coupling A. Liebsch and T. A. Costi Institut für Festkörperforschung, Forschungszentrum
More informationRole of Hund Coupling in Two-Orbital Systems
Role of Hund Coupling in Two-Orbital Systems Gun Sang Jeon Ewha Womans University 2013-08-30 NCTS Workshop on Quantum Condensation (QC13) collaboration with A. J. Kim, M.Y. Choi (SNU) Mott-Hubbard Transition
More information4 Development of the LDA+DMFT Approach
4 Development of the LDA+DMFT Approach Alexander Lichtenstein I. Institut für Theoretische Physik Universität Hamburg Contents 1 Introduction 2 2 Functional approach: from DFT to DMFT 4 3 Local correlations:
More informationDMFT for correlated bosons and boson-fermion mixtures
DMFT for correlated bosons and boson-fermion mixtures Workshop on Recent developments in dynamical mean-field theory ETH ürich, September 29, 2009 Dieter Vollhardt Supported by Deutsche Forschungsgemeinschaft
More informationNUMERICAL METHODS FOR QUANTUM IMPURITY MODELS
NUMERICAL METODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Andrew Mitchell, Utrecht University Quantum impurity problems Part 1: Quantum impurity problems
More informationFrom Gutzwiller Wave Functions to Dynamical Mean-Field Theory
From utzwiller Wave Functions to Dynamical Mean-Field Theory Dieter Vollhardt Autumn School on Correlated Electrons DMFT at 25: Infinite Dimensions Forschungszentrum Jülich, September 15, 2014 Supported
More informationMott physics: from basic concepts to iron superconductors
Mott physics: from basic concepts to iron superconductors E. Bascones Instituto de Ciencia de Materiales de Madrid (ICMM-CSIC) Outline Mott physics: Basic concepts (single orbital & half filling) - Mott
More informationCooperative Phenomena
Cooperative Phenomena Frankfurt am Main Kaiserslautern Mainz B1, B2, B4, B6, B13N A7, A9, A12 A10, B5, B8 Materials Design - Synthesis & Modelling A3, A8, B1, B2, B4, B6, B9, B11, B13N A5, A7, A9, A12,
More informationTuning magnetic anisotropy, Kondo screening and Dzyaloshinskii-Moriya interaction in pairs of Fe adatoms
Tuning magnetic anisotropy, Kondo screening and Dzyaloshinskii-Moriya interaction in pairs of Fe adatoms Department of Physics, Hamburg University, Hamburg, Germany SPICE Workshop, Mainz Outline Tune magnetic
More informationA brief Introduction of Fe-based SC
Part I: Introduction A brief Introduction of Fe-based SC Yunkyu Bang (Chonnam National Univ., Kwangju, Korea) Lecture 1: Introduction 1. Overview 2. What is sign-changing s-wave gap : +/-s-wave gap Lecture
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationNumerical Methods in Quantum Many-body Theory. Gun Sang Jeon Pyeong-chang Summer Institute 2014
Numerical Methods in Quantum Many-body Theory Gun Sang Jeon 2014-08-25 Pyeong-chang Summer Institute 2014 Contents Introduction to Computational Physics Monte Carlo Methods: Basics and Applications Numerical
More informationMOTTNESS AND STRONG COUPLING
MOTTNESS AND STRONG COUPLING ROB LEIGH UNIVERSITY OF ILLINOIS Rutgers University April 2008 based on various papers with Philip Phillips and Ting-Pong Choy PRL 99 (2007) 046404 PRB 77 (2008) 014512 PRB
More informationLattice modulation experiments with fermions in optical lattices and more
Lattice modulation experiments with fermions in optical lattices and more Nonequilibrium dynamics of Hubbard model Ehud Altman Weizmann Institute David Pekker Harvard University Rajdeep Sensarma Harvard
More informationTRIQS : A Toolbox for Research in Interacting Quantum Systems Introduction to Hands-on
1 TRIQS : A Toolbox for Research in Interacting Quantum Systems Introduction to Hands-on Olivier Parcollet Center for Computational Quantum Physics (CCQ) Flatiron Institute, Simons Foundation New York
More informationIdeas on non-fermi liquid metals and quantum criticality. T. Senthil (MIT).
Ideas on non-fermi liquid metals and quantum criticality T. Senthil (MIT). Plan Lecture 1: General discussion of heavy fermi liquids and their magnetism Review of some experiments Concrete `Kondo breakdown
More informationMany-body effects in iron pnictides and chalcogenides
Many-body effects in iron pnictides and chalcogenides separability of non-local and dynamical correlation effects Jan M. Tomczak Vienna University of Technology jan.tomczak@tuwien.ac.at Emergent Quantum
More informationQuantum Impurities In and Out of Equilibrium. Natan Andrei
Quantum Impurities In and Out of Equilibrium Natan Andrei HRI 1- Feb 2008 Quantum Impurity Quantum Impurity - a system with a few degrees of freedom interacting with a large (macroscopic) system. Often
More information2D Bose and Non-Fermi Liquid Metals
2D Bose and Non-Fermi Liquid Metals MPA Fisher, with O. Motrunich, D. Sheng, E. Gull, S. Trebst, A. Feiguin KITP Cold Atoms Workshop 10/5/2010 Interest: A class of exotic gapless 2D Many-Body States a)
More informationSimulations of Quantum Dimer Models
Simulations of Quantum Dimer Models Didier Poilblanc Laboratoire de Physique Théorique CNRS & Université de Toulouse 1 A wide range of applications Disordered frustrated quantum magnets Correlated fermions
More informationWang-Landau sampling for Quantum Monte Carlo. Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart
Wang-Landau sampling for Quantum Monte Carlo Stefan Wessel Institut für Theoretische Physik III Universität Stuttgart Overview Classical Monte Carlo First order phase transitions Classical Wang-Landau
More informationZhiping Yin. Department of Physics, Rutgers University Collaborators: G. Kotliar, K. Haule
DFT+DMFT to Correlated Electronic Structures: Recent Developments and Applications to Iron-based Superconductors Zhiping Yin Department of Physics, Rutgers University Collaborators: G. Kotliar, K. Haule
More informationPHYSICAL REVIEW B 80,
PHYSICAL REVIEW B 8, 6526 29 Finite-temperature exact diagonalization cluster dynamical mean-field study of the two-dimensional Hubbard model: Pseudogap, non-fermi-liquid behavior, and particle-hole asymmetry
More informationTopological Kondo Insulators!
Topological Kondo Insulators! Maxim Dzero, University of Maryland Collaborators: Kai Sun, University of Maryland Victor Galitski, University of Maryland Piers Coleman, Rutgers University Main idea Kondo
More informationFluctuating exchange theory of dynamical electron correlations and magnetism
Fluctuating exchange theory of dynamical electron correlations and magnetism Václav Drchal Institute of Physics ASCR, Praha, Czech Republic Grant Agency of ASCR: project IAA11616 Workshop Frontiers in
More informationNew perspectives in superconductors. E. Bascones Instituto de Ciencia de Materiales de Madrid (ICMM-CSIC)
New perspectives in superconductors E. Bascones Instituto de Ciencia de Materiales de Madrid (ICMM-CSIC) E. Bascones leni@icmm.csic.es Outline Talk I: Correlations in iron superconductors Introduction
More informationExcitonic Condensation in Systems of Strongly Correlated Electrons. Jan Kuneš and Pavel Augustinský DFG FOR1346
Excitonic Condensation in Systems of Strongly Correlated Electrons Jan Kuneš and Pavel Augustinský DFG FOR1346 Motivation - unconventional long-range order incommensurate spin spirals complex order parameters
More informationMarkov Chain Monte Carlo The Metropolis-Hastings Algorithm
Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Anthony Trubiano April 11th, 2018 1 Introduction Markov Chain Monte Carlo (MCMC) methods are a class of algorithms for sampling from a probability
More informationStrongly Correlated Systems of Cold Atoms Detection of many-body quantum phases by measuring correlation functions
Strongly Correlated Systems of Cold Atoms Detection of many-body quantum phases by measuring correlation functions Anatoli Polkovnikov Boston University Ehud Altman Weizmann Vladimir Gritsev Harvard Mikhail
More informationSecond Lecture: Quantum Monte Carlo Techniques
Second Lecture: Quantum Monte Carlo Techniques Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml aml@pks.mpg.de Lecture Notes at http:www.pks.mpg.de/~aml/leshouches
More informationBose-Hubbard Model (BHM) at Finite Temperature
Bose-Hubbard Model (BHM) at Finite Temperature - a Layman s (Sebastian Schmidt) proposal - pick up Diploma work at FU-Berlin with PD Dr. Axel Pelster (Uni Duisburg-Essen) ~ Diagrammatic techniques, high-order,
More informationMott transition : beyond Dynamical Mean Field Theory
Mott transition : beyond Dynamical Mean Field Theory O. Parcollet 1. Cluster methods. 2. CDMFT 3. Mott transition in frustrated systems : hot-cold spots. Coll: G. Biroli (SPhT), G. Kotliar (Rutgers) Ref:
More informationA New look at the Pseudogap Phase in the Cuprates.
A New look at the Pseudogap Phase in the Cuprates. Patrick Lee MIT Common themes: 1. Competing order. 2. superconducting fluctuations. 3. Spin gap: RVB. What is the elephant? My answer: All of the above!
More informationThe 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007
The 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007 Kondo Effect in Metals and Quantum Dots Jan von Delft
More informationDMFT and beyond : IPAM, Los Angeles, Jan. 26th 2009 O. Parcollet Institut de Physique Théorique CEA-Saclay, France
DMFT and beyond : From quantum impurities to high temperature superconductors 1 IPAM, Los Angeles, Jan. 26th 29 O. Parcollet Institut de Physique Théorique CEA-Saclay, France Coll : M. Ferrero (IPhT),
More informationFate of the Kondo impurity in a superconducting medium
Karpacz, 2 8 March 214 Fate of the Kondo impurity in a superconducting medium T. Domański M. Curie Skłodowska University Lublin, Poland http://kft.umcs.lublin.pl/doman/lectures Motivation Physical dilemma
More informationSpin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University
Spin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University Supported by U.S. DoE Basic Energy Sciences, Materials Sciences & Engineering DE-FG02-08ER46544 Overview
More informationMany body physics issues with the core-hole propagator and resonant inelastic X-ray scattering
Many body physics issues with the core-hole propagator and resonant inelastic X-ray scattering Jim Freericks (Georgetown University) Funding: Civilian Research and Development Foundation In collaboration
More informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationNeutron scattering from quantum materials
Neutron scattering from quantum materials Bernhard Keimer Max Planck Institute for Solid State Research Max Planck UBC UTokyo Center for Quantum Materials Detection of bosonic elementary excitations in
More informationSpectral Density Functional Theory
Spectral Density Functional Theory Sergej Savrasov Financial support NSF US DOE LANL Collaborators and Content Constructing New Functionals to Access Energetics and Spectra of Correlated Solids Phonons
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationDynamical phase transition and prethermalization. Mobile magnetic impurity in Fermi superfluids
Dynamical phase transition and prethermalization Pietro Smacchia, Alessandro Silva (SISSA, Trieste) Dima Abanin (Perimeter Institute, Waterloo) Michael Knap, Eugene Demler (Harvard) Mobile magnetic impurity
More informationANTIFERROMAGNETIC EXCHANGE AND SPIN-FLUCTUATION PAIRING IN CUPRATES
ANTIFERROMAGNETIC EXCHANGE AND SPIN-FLUCTUATION PAIRING IN CUPRATES N.M.Plakida Joint Institute for Nuclear Research, Dubna, Russia CORPES, Dresden, 26.05.2005 Publications and collaborators: N.M. Plakida,
More informationLocal moment approach to multi-orbital Anderson and Hubbard models
Local moment approach to multi-orbital Anderson and Hubbard models Anna Kauch 1 and Krzysztof Byczuk,1 1 Institute of Theoretical Physics, Warsaw University, ul. Hoża 69, PL--681 Warszawa, Poland Theoretical
More informationAnderson impurity model at finite Coulomb interaction U: Generalized noncrossing approximation
PHYSICAL REVIEW B, VOLUME 64, 155111 Anderson impurity model at finite Coulomb interaction U: Generalized noncrossing approximation K. Haule, 1,2 S. Kirchner, 2 J. Kroha, 2 and P. Wölfle 2 1 J. Stefan
More informationSimple Explanation of Fermi Arcs in Cuprate Pseudogaps: A Motional Narrowing Phenomenon
Simple Explanation of Fermi Arcs in Cuprate Pseudogaps: A Motional Narrowing Phenomenon ABSTRACT: ARPES measurements on underdoped cuprates above the superconducting transition temperature exhibit the
More informationA Holographic Model of the Kondo Effect (Part 1)
A Holographic Model of the Kondo Effect (Part 1) Andy O Bannon Quantum Field Theory, String Theory, and Condensed Matter Physics Kolymbari, Greece September 1, 2014 Credits Based on 1310.3271 Johanna Erdmenger
More informationGround State Projector QMC in the valence-bond basis
Quantum Monte Carlo Methods at Work for Novel Phases of Matter Trieste, Italy, Jan 23 - Feb 3, 2012 Ground State Projector QMC in the valence-bond basis Anders. Sandvik, Boston University Outline: The
More informationTensor network simulations of strongly correlated quantum systems
CENTRE FOR QUANTUM TECHNOLOGIES NATIONAL UNIVERSITY OF SINGAPORE AND CLARENDON LABORATORY UNIVERSITY OF OXFORD Tensor network simulations of strongly correlated quantum systems Stephen Clark LXXT[[[GSQPEFS\EGYOEGXMZMXMIWUYERXYQGSYVWI
More informationORBITAL SELECTIVITY AND HUND S PHYSICS IN IRON-BASED SC. Laura Fanfarillo
ORBITAL SELECTIVITY AND HUND S PHYSICS IN IRON-BASED SC Laura Fanfarillo FROM FERMI LIQUID TO NON-FERMI LIQUID Strong Correlation Bad Metal High Temperature Fermi Liquid Low Temperature Tuning parameter
More informationElectron Correlation
Series in Modern Condensed Matter Physics Vol. 5 Lecture Notes an Electron Correlation and Magnetism Patrik Fazekas Research Institute for Solid State Physics & Optics, Budapest lb World Scientific h Singapore
More informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More information