6. Auxiliary field continuous time quantum Monte Carlo
|
|
- Garey Ferguson
- 5 years ago
- Views:
Transcription
1 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 model H µ n +Un n + k t k a k a k + k V k a k c k +H.c. 6.1 Diagrammatic Monte Carlo In order to evaluate the partition function, we write the Hamiltonian as H H + V so that we have { β } e βh e βh exp T τ dτ ˆVτ 6.2 with ˆVτ in the interaction picture ˆVτ e +τh Ve τh. and the usual imaginary time ordering operatore T τ. In Diagrammatic Monte-Carlo, the Feynman diagrams of the perturbation expansion are sampled by Monte Carlo. Then, the partition function is Tr Z Tr e βh Z with Z Tre βh. Thus Z 1 { β Tr e βh exp T τ Z Z { } e βh β exp T τ dτ ˆVτ Tre βh dτ ˆVτ } { β exp T τ 6.3 dτ ˆVτ } 1 This chapter is based on a course given by Aaram J. Kim. The master thesis Topological Phases of Interacting Fermions in Optical Lattices with Artificial Gauge Fields by Michael Buchhold is used as a reference. 77
2 6.4 This expectation value will be calculated with quantum Monte Carlo. First, we write more explicitly { β exp T τ dτ ˆVτ } 1 β β dτ k... dτ 1 Vτk... Vτ 1 k! k β β β dτ k... dτ 2 dτ 1 Vτk... Vτ 1 τ k 1 τ 1 k because if τ e +τh e τh, then 6.5 τ τ e +τh H He τh e +τh Ve τh τ Vτ τ This is formally solved by integration: β Iteration via τ β τ 6.6 dτvττ 6.7 dτ Vτ τ 6.8 leads to the expression 6.5. Then Z { β exp T τ dτ Z ˆVτ } 1 β β dτ k... dτ 1 ˆVτ k... ˆVτ 1 k! k Wx dxwx 6.9 x Here, we interpret Wx as a weight. In order to be able to interpret Wx as a probability, we assume it is positive. In the continuous-time auxiliary field quantum Monte-Carlo algorithm, an auxilary field decomposition of the interaction part of the Anderson impurity model is used. In CT-AUX, the choice of interaction is V U n n n + n 2 K β
3 where the second term is a Hartree like term, and the third term is a control parameter for the expansion order. This means that the noninteracting Hamiltonian is modified as H H V 6.11 The choice of V yields the following values for empty, singly occupied and doubly ocupied impurity: K β K β U 2 K β U 2 K β 6.12 Now a Hubbard-tratonovich transformation is applied to V: V U n n n + n K 2 β K e γsn n s± The identity can be checked by applying the last expression on the basis states in For example, for, n n, and the expression yields K β. For, n 1, n V K β cosh γ cosh γ 1 + Uβ 2K 6.14 s represents an Ising type boson that was introduced, an auxiliary spin. Eventually there will be k auxiliary spins over which we sum: Z + K k [ e γs kn τ k... e γs 1n τ 1 ] 6.15 Z {s i } Here, is the physical spin, with numerical value +1 for and 1 for. Then, for Monte Carlo sampling we identify Z dτ k... dτ 1 Wx 6.16 Z k {s i } 79
4 s 2 1 s 4 +1 s 1 +1 s 3 +1 s 5 1 β τ 1 τ 2 τ 3 τ 4 τ 5 Figure 6.1: Example of a CT-AUX configuration at perturbation order 5. The configuration space consists of perturbation order k, time indices τ and auxiliary spin indices s i ; there are k time points. The random walker needs to roam the entire configuration space. CT-AUX is a generalization of the discrete time quantum Monte Carlo algorithm Hirsch-Fye. Explicit calculation of the matrix We now prepare to write the action in terms of Grassmann variables by rewriting the interaction V K e γsn n K e γsn 6.17 s±1 s±1 further. Taylor expansion of the exponential yields because n 2 n e γsn 1 + γsn + 1 2! γs2 n n + n + γsn + 1 2! γs2 n n + e γs n 1 1 e γs n 1 1 e γs 1 c c 1 [ ] 1 e γs c c + e γs c c e γs e γs 1c c e γs 1 1 e γs c c 6.18 where Fermion antiperiodicity c c 1 c c was used. Thus, we have found that V K e γsn n K [ ] e γs e γs 1c c 6.19 s±1 s±1, Now we can go from Hamiltonian to action formalism, writing for the average of an operator A with respect to the noninteracting Hamiltonian 8
5 H A Tr[ ] Ae βh Tr [ ] e βh D[c + c]ae D[c+ c]e 6.2 While in the first expression, c, c are operators, in the second, c and c + are Grassmann variables. We now come back to the effective action of the previous chapter eff dτdτ c + τg 1 τ τ cτ + dτvτ 6.21 where the Weiss function G encodes the result of integrating out the bath degrees of freedom. We will use a number of identities for Grassmann variables: D[c + c]e ij c+ i ijc j det[] D[c + c]c x c + y e ij c+ i ijc j D[c+ c]e ij c+ i ijc j 1 xy e ac+c 1 ac + c 6.22 We now introduce the matrix elements ij of the action with i, j indices for time, spin, lattice sites, etc., i, j,..., M where M ; as we have passed from continous imaginary time τ integrals to sums, the τ i have to be infinitely dense. There are infinitely more times τ i than the k times we have in a given configuration of perturbation order k. Thus D[c + c]e D[c + c]e ij c+ i ijc j with ij G 1 ij τ i τ j 6.23 For the interaction, we introduce the matrices A and B with { δ ij e γs i i for τ i {τ 1... τ k } A ij δ ij otherwise {δ ij 1 e γs i i for τ i {τ 1... τ k } B ij otherwise 6.24 The time mesh for the Feynman path integral now has the imaginary times τ 1 to τ k as well as an infinite number of additional time indices τ i. A has, 81
6 for perturbation order k, k exponentials on the diagonal for times τ 1, τ 2,..., τ k, and otherwise 1: 1 e γs 1 A e γs 2... and consequently, B is 1 e γs 1 B 1 e γs Then we can write the average of Eq. 6.9 as V... V deta D[c + c]e c+ i B ij c j e c+ i ij c j det 6.27 The product over spin reflects the fact that the noninteracting part of the Hamiltonian H contains only quadratic operators with index so that H commutes with H so that integration weights Wx are products of weights for the individual spin components,. If we now rewrite the exponentials as e c+ i B ij c j we find for the integral we need to compute I deta det B det 6.28 Due to the M index, matrices A, are infinite dimensional. The idea is now to use the zeros in the B matrix to make things finite-dimensional by reordering. det B detdet1 B
7 because B 1 B 1, and the determinant of a product splits into a product of determinants. We now exploit that at perturbation order k, only k entries in the B matrix are nonzero: b 1 b 1 B... B b 2... b k b k This rearrangement doesn t produce a sign as rows and columns are exchanged. In the same way, 1 is rearranged into 1 only that this doesn t produce empty submatrices. Then and B 1 x B 1 x x 1 B 1 T 1 B 1 T 1 but in the determinant, only 1 B 1 remains: det1 B 1 det1 B In Eq. 6.28, det cancels with the denominator, so that we have I deta det1 B 1 detã det1 B because the rearranged à has only entries of 1 beyond the reduced space of dimension k. Now we arrive at the final expression I det à à B 1 e Γ e Γ 1G detn where the k k matrix G, with k numbering the time points, was introduced, and the entire matrix that needs to be manipulated is now called N 1 matrix N matrix. The notation for à is à e Γ diag e γs1,..., e γs k 83
8 Later, the Weiss function G will be represented by a very large mesh with many more than k points; it will be interpolated linearly or by cubic methods from k points. Insertion and removal Now, we will assume that we know the N 1 matrix at a given perturbation order k. To implement a random walk through configuration space, we need to either insert a time point or remove a time point, based on the Metropolis algorithm. a k4 β τ 1 τ 2 τ 3 τ 4 b k5 insertion c k3 removal β τ 1 τ 2 τ 3 τ 4 τ β τ 1 τ 2 τ 3 Figure 6.2: chematic of insertion and removal updates to a configuration. Insertion.- The decision to insert an auxiliary spin into the configuration is taken based on a uniform random number r in the range [, 1]. If r > min 1, f W k+1 W k 6.36 where f is factor explained later, the update is accepted, otherwise rejected. Thus, to take the decision, the ratio of weights needs to be calculated, specifically: det N k+1 1 det N k 1 det N k 1 det N k N k Q det RN k Q 6.37 R Nk 1 Nk 1 R Q 84
9 This uses the fact that N k and N k 1 are inverses of each other, and matrices are multiplied using the properties of determinants. During the simulation N k rather than N k 1 is stored. For the enlarged matrix N k 1 we have N k+1 1 N k 1 Q 6.38 R with k 1 vector Q, 1 k vector R and defined by 1 l n: Q l e γs l 1G τ l τ R l e γs 1G τ τ l e γs e γs 1G In case the insertion move is accepted, we now need the missing elements of the enlarged matrix N k+1, expressed in terms of the additional elements Q, R, of N k+1 1 which we know: N k+1[ N k+1] 1 P Q P Q 1I R 6.4 R where we write P for N k 1. This yields the equations PP + QR 1I PQ + Q RP + R T RQ which we have to solve for PP, Q, R, and so on. The third equation yields RP R R RP and inserting into the fourth equation RP 1 Q + 1 RP 1 Q 1 1 RP 1 Q 6.43 with P 1 N k. is a number which we don t need to replace so that we have from the third equation R RN k
10 The second equation Q PQ is inserted into the first PP PQR 1 P P QR where Q is a column vector, R is a row vector. Now we use the herman- Morrison formula for an invertible square matrix A, with v T A 1 u 1 A + uv T 1 A 1 A 1 uv T A 1 1 v T A 1 u 6.47 which allows a cheap update of the already known inverse of A if A is modified by the rank 1 matrix uv T. Here, we have u Q, vt R and P P QR 1 [ P P 1 1 Q ] [ RP 1 ] 1 RP 1 Q [ N N k k Q ] [RN k] + RN k Q Nk + [ N k Q ] [RN k] 6.48 This is the most massive calculation in the update. Now we find for Q Q PQ Nk Q + [ N k Q ] [RN k] Q Nk Q 1 + RNk Q N k Q RN k Q RN k Q N k Q 6.49 Removal.- The decision on insertion versus removal is taken with 5% probability. For the deletion of an auxiliary spin s at time τ from the matrix N k, the situation is N k+1 1 N k 1 Q P N Q k+1 R R with RN k Q 1 Q N k Q R RN k P N k + [ N k Q ] [RN k]
11 The idea is now to use to solve the P equation for N k using the expressions for Q and R: P N k + Q R Nk + Q R N k P Q R 6.51 Monte Carlo procedure.- To summarize, we use the expansion of the partition function β τ k K Z dτ k dτ 1 Z k {si, τ i } 6.52 where k {s i } Z k {si, τ i } Z detn 1 {si, τ i } 6.53 For the specific formulation of the update, we need the condition of detailed balance Wxpx x Wx px x 6.54 Then in general, the Metropolis algorithm is to update, with a random number r, if r < min 1, Wx 6.55 Wx The detailed balance condition can be slightly modified by splitting the transition probability into proposal and acceptance probabilities: p p proposal p acceptance 6.56 For an insertion, the proposal probability is the probability p 1 of choosing a time τ from the interval [1, β], p 1 dτ β, times the probability p of choosing the auxiliary spin s from { 1, 1}, i.e. p proposal k k + 1 p 1 p 2 dτ 6.57 For deletion, the tuple τ i, s i must be chosen from the existing k + 1 auxiliary spins, so that the proposal probability is p proposal k + 1 k 1 n
12 Then, the Metropolis update decision is modified to r < min 1, Wx p proposal x x Wxp proposal x x 6.59 Detailed balance then yields for the acceptance probabilities p acceptance k k + 1 p acceptance k + 1 k because W k+1 dτ 1... dτ k dτ k+1 K W k dτ 1... dτ k K W k+1 K detn 1 k + 1dτ W k k + 1 k+1 detn 1 k + 1 k detn 1 k k + 1 detn 1 k 6.6 This means that p acceptance k k + 1 min p acceptance k + 1 k min { { 1, K k + 1 1, k + 1 K detn 1 k + 1 detn 1 k detn 1 k detn 1 k + 1 Measurement of the Greens function The Greens function, written for QMC without a minus sign, is G τ τ T τ c τc τ } } 6.61 D[c + c]c x c + ye ij c+ i B ijc j D[c+ c]e ij c+ i B ijc j [ B 1] xy 6.62 where x, y stand for arbitrary imaginary time indices. For the matrix inverse, we can write B 1 [ 1 B 1 ] B B 1 1[ 1 B 1 + B 1] B 1 1 B
13 Here we recapitulate the method of getting rid of the zero subspace: The matrices, B are infinite dimensional but B has many zeros: B B B B 1 1 T B 1 1 B 1 T 1 B B B B 1 1 B B 6.64 Now we go back to the definitions of the A, B and matrices: A e Γ B 1 e Γ 1 G 6.65 with the noninteracting Weiss function G ; thus is invertible. Now we insert Ã: 1 B B 1 Ã 1Ã 1Ã B B Ã B 1 1Ã B N {si, τ i } e Γ 1 M 6.66 Then, coming back to the matrix inverse 6.63 B 1 G + G N {si, τ i } e Γ 1 This means that the full Greens function can be evaluated as G 6.67 G τ τ G τ τ + G τ τ i M ijg τ j τ 6.68 where τ i, τ j are configuration time points. The k k matrix M is M N e Γ The Weiss function G does not depend on the configuration, but M does. Therefore, we only need to accumulate terms that depend on the configuration, in particular M ij G τ j τ. This is accumulated for many i, and only a vector needs to be stored. This also explains why the N matrix is stored; it is needed for the Greens function. The matrix update cost is of order Ok 2 because only matrix vector multiplications are needed. 89
Quantum Monte Carlo method in details.
Quantum Monte Carlo method in details. I. CONTINUOUS AND DISCRETE HUBBARD-STRATONOVICH TRANSFORMATIONS The HS transformation is based on identity : exp{ 1 2 A2 } = (2π) 1/2 + dx exp{ 1 2 x2 xa} (1) This
More informationDiffusion Monte Carlo
Diffusion Monte Carlo Notes for Boulder Summer School 2010 Bryan Clark July 22, 2010 Diffusion Monte Carlo The big idea: VMC is a useful technique, but often we want to sample observables of the true ground
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 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 information1 Differentiable manifolds and smooth maps. (Solutions)
1 Differentiable manifolds and smooth maps Solutions Last updated: March 17 2011 Problem 1 The state of the planar pendulum is entirely defined by the position of its moving end in the plane R 2 Since
More information1 Differentiable manifolds and smooth maps. (Solutions)
1 Differentiable manifolds and smooth maps Solutions Last updated: February 16 2012 Problem 1 a The projection maps a point P x y S 1 to the point P u 0 R 2 the intersection of the line NP with the x-axis
More information7 Monte Carlo Simulations of Quantum Spin Models
7 Monte Carlo Simulations of Quantum Spin Models Stefan Wessel Institute for Theoretical Solid State Physics RWTH Aachen University Contents 1 Introduction 2 2 World lines and local updates 2 2.1 Suzuki-Trotter
More informationPhysics 127c: Statistical Mechanics. Application of Path Integrals to Superfluidity in He 4
Physics 17c: Statistical Mechanics Application of Path Integrals to Superfluidity in He 4 The path integral method, and its recent implementation using quantum Monte Carlo methods, provides both an intuitive
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 informationA continuous time algorithm for quantum impurity models
ISSP, Aug. 6 p.1 A continuou time algorithm for quantum impurity model Philipp Werner Department of Phyic, Columbia Univerity cond-mat/512727 ISSP, Aug. 6 p.2 Outline Introduction Dynamical mean field
More informationQuantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC
Quantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC Itay Hen Information Sciences Institute, USC Workshop on Theory and Practice of AQC and Quantum Simulation Trieste, Italy August
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 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 informationQuantum Simulation Studies of Charge Patterns in Fermi-Bose Systems
Quantum Simulation Studies of Charge Patterns in Fermi-Bose Systems 1. Hubbard to Holstein 2. Peierls Picture of CDW Transition 3. Phonons with Dispersion 4. Holstein Model on Honeycomb Lattice 5. A New
More informationBy following steps analogous to those that led to (20.177), one may show (exercise 20.30) that in Feynman s gauge, = 1, the photon propagator is
20.2 Fermionic path integrals 74 factor, which cancels. But if before integrating over all gauge transformations, we shift so that 4 changes to 4 A 0, then the exponential factor is exp[ i 2 R ( A 0 4
More informationLecture 1: The Equilibrium Green Function Method
Lecture 1: The Equilibrium Green Function Method Mark Jarrell April 27, 2011 Contents 1 Why Green functions? 2 2 Different types of Green functions 4 2.1 Retarded, advanced, time ordered and Matsubara
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 informationPhysics 239/139 Spring 2018 Assignment 2 Solutions
University of California at San Diego Department of Physics Prof. John McGreevy Physics 39/139 Spring 018 Assignment Solutions Due 1:30pm Monday, April 16, 018 1. Classical circuits brain-warmer. (a) Show
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationMath 215 HW #9 Solutions
Math 5 HW #9 Solutions. Problem 4.4.. If A is a 5 by 5 matrix with all a ij, then det A. Volumes or the big formula or pivots should give some upper bound on the determinant. Answer: Let v i be the ith
More information1 Last time: determinants
1 Last time: determinants Let n be a positive integer If A is an n n matrix, then its determinant is the number det A = Π(X, A)( 1) inv(x) X S n where S n is the set of n n permutation matrices Π(X, A)
More informationGaussian Elimination and Back Substitution
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 4 Notes These notes correspond to Sections 31 and 32 in the text Gaussian Elimination and Back Substitution The basic idea behind methods for solving
More informationM. A. Gusmão IF-UFRGS
M. A. Gusmão IF-UFRGS - 217 1 FIP164-217/2 Text 9 Mean-field approximation - II Heisenberg Hamiltonian in wave-vector space As we saw in Text 8, the uniform susceptibility does not diverge in the case
More informationPhysics 115/242 Monte Carlo simulations in Statistical Physics
Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationQuantum Monte Carlo Methods.
Quantum Monte Carlo Methods. Not yet. I will start from the alternative formulation of quantum mechanics in terms of path integrals. First, for the two level system, then for the many-body lattice system,
More informationAdvanced Monte Carlo Methods Problems
Advanced Monte Carlo Methods Problems September-November, 2012 Contents 1 Integration with the Monte Carlo method 2 1.1 Non-uniform random numbers.......................... 2 1.2 Gaussian RNG..................................
More informationHomework Set #8 Solutions
Exercises.2 (p. 19) Homework Set #8 Solutions Assignment: Do #6, 8, 12, 14, 2, 24, 26, 29, 0, 2, 4, 5, 6, 9, 40, 42 6. Reducing the matrix to echelon form: 1 5 2 1 R2 R2 R1 1 5 0 18 12 2 1 R R 2R1 1 5
More informationKitaev honeycomb lattice model: from A to B and beyond
Kitaev honeycomb lattice model: from A to B and beyond Jiri Vala Department of Mathematical Physics National University of Ireland at Maynooth Postdoc: PhD students: Collaborators: Graham Kells Ahmet Bolukbasi
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
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 informationlattice QCD and the hadron spectrum Jozef Dudek ODU/JLab
lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab the light meson spectrum relatively simple models of hadrons: bound states of constituent quarks and antiquarks the quark model empirical meson
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationPhysics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion
Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to
More informationVacuum Energy and Effective Potentials
Vacuum Energy and Effective Potentials Quantum field theories have badly divergent vacuum energies. In perturbation theory, the leading term is the net zero-point energy E zeropoint = particle species
More informationlattice QCD and the hadron spectrum Jozef Dudek ODU/JLab
lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab a black box? QCD lattice QCD observables (scattering amplitudes?) in these lectures, hope to give you a look inside the box 2 these lectures how
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More information4 4 and perturbation theory
and perturbation theory We now consider interacting theories. A simple case is the interacting scalar field, with a interaction. This corresponds to a -body contact repulsive interaction between scalar
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 24 Jan 2000
The sign problem in Monte Carlo simulations of frustrated quantum spin systems Patrik Henelius and Anders W. Sandvik 2 National High Magnetic Field Laboratory, 800 East Paul Dirac Dr., Tallahassee, Florida
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationLocalization I: General considerations, one-parameter scaling
PHYS598PTD A.J.Leggett 2013 Lecture 4 Localization I: General considerations 1 Localization I: General considerations, one-parameter scaling Traditionally, two mechanisms for localization of electron states
More informationLecture 9: Numerical Linear Algebra Primer (February 11st)
10-725/36-725: Convex Optimization Spring 2015 Lecture 9: Numerical Linear Algebra Primer (February 11st) Lecturer: Ryan Tibshirani Scribes: Avinash Siravuru, Guofan Wu, Maosheng Liu Note: LaTeX template
More informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationMatrix Operations: Determinant
Matrix Operations: Determinant Determinants Determinants are only applicable for square matrices. Determinant of the square matrix A is denoted as: det(a) or A Recall that the absolute value of the determinant
More informationPolynomial Filtered Hybrid Monte Carlo
Polynomial Filtered Hybrid Monte Carlo Waseem Kamleh and Michael J. Peardon CSSM & U NIVERSITY OF A DELAIDE QCDNA VII, July 4th 6th, 2012 Introduction Generating large, light quark dynamical gauge field
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More information33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM
33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM (UPDATED MARCH 17, 2018) The final exam will be cumulative, with a bit more weight on more recent material. This outline covers the what we ve done since the
More informationTime Evolving Block Decimation Algorithm
Time Evolving Block Decimation Algorithm Application to bosons on a lattice Jakub Zakrzewski Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, Jagiellonian University,
More informationDiagrammatic Green s Functions Approach to the Bose-Hubbard Model
Diagrammatic Green s Functions Approach to the Bose-Hubbard Model Matthias Ohliger Institut für Theoretische Physik Freie Universität Berlin 22nd of January 2008 Content OVERVIEW CONSIDERED SYSTEM BASIC
More informationMagnets, 1D quantum system, and quantum Phase transitions
134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider
More informationFinal EXAM Preparation Sheet
Final EXAM Preparation Sheet M369 Fall 217 1 Key concepts The following list contains the main concepts and ideas that we have explored this semester. For each concept, make sure that you remember about
More informationReview Questions REVIEW QUESTIONS 71
REVIEW QUESTIONS 71 MATLAB, is [42]. For a comprehensive treatment of error analysis and perturbation theory for linear systems and many other problems in linear algebra, see [126, 241]. An overview of
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationGreen s functions: calculation methods
M. A. Gusmão IF-UFRGS 1 FIP161 Text 13 Green s functions: calculation methods Equations of motion One of the usual methods for evaluating Green s functions starts by writing an equation of motion for the
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6
CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6 GENE H GOLUB Issues with Floating-point Arithmetic We conclude our discussion of floating-point arithmetic by highlighting two issues that frequently
More information7. Dimension and Structure.
7. Dimension and Structure 7.1. Basis and Dimension Bases for Subspaces Example 2 The standard unit vectors e 1, e 2,, e n are linearly independent, for if we write (2) in component form, then we obtain
More informationEXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)
EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily
More informationDiagonalization. MATH 1502 Calculus II Notes. November 4, 2008
Diagonalization MATH 1502 Calculus II Notes November 4, 2008 We want to understand all linear transformations L : R n R m. The basic case is the one in which n = m. That is to say, the case in which the
More informationChapter 4 & 5: Vector Spaces & Linear Transformations
Chapter 4 & 5: Vector Spaces & Linear Transformations Philip Gressman University of Pennsylvania Philip Gressman Math 240 002 2014C: Chapters 4 & 5 1 / 40 Objective The purpose of Chapter 4 is to think
More informationTopological insulator part II: Berry Phase and Topological index
Phys60.nb 11 3 Topological insulator part II: Berry Phase and Topological index 3.1. Last chapter Topological insulator: an insulator in the bulk and a metal near the boundary (surface or edge) Quantum
More informationMATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.
MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If
More informationMethods of Fermion Monte Carlo
Methods of Fermion Monte Carlo Malvin H. Kalos kalos@llnl.gov Institute for Nuclear Theory University of Washington Seattle, USA July, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of
More informationMath 215B: Solutions 1
Math 15B: Solutions 1 Due Thursday, January 18, 018 (1) Let π : X X be a covering space. Let Φ be a smooth structure on X. Prove that there is a smooth structure Φ on X so that π : ( X, Φ) (X, Φ) is an
More informationThis MUST hold matrix multiplication satisfies the distributive property.
The columns of AB are combinations of the columns of A. The reason is that each column of AB equals A times the corresponding column of B. But that is a linear combination of the columns of A with coefficients
More informationBeta functions in quantum electrodynamics
Beta functions in quantum electrodynamics based on S-66 Let s calculate the beta function in QED: the dictionary: Note! following the usual procedure: we find: or equivalently: For a theory with N Dirac
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationFermi liquid theory Can we apply the free fermion approximation to a real metal? Phys540.nb Strong interaction vs.
Phys540.nb 7 Fermi liquid theory.. Can we apply the free fermion approximation to a real metal? Can we ignore interactions in a real metal? Experiments says yes (free fermion models work very well), but
More informationMath 200 University of Connecticut
RELATIVISTIC ADDITION AND REAL ADDITION KEITH CONRAD Math 200 University of Connecticut Date: Aug. 31, 2005. RELATIVISTIC ADDITION AND REAL ADDITION 1 1. Introduction For three particles P, Q, R travelling
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationFinite temperature QFT: A dual path integral representation
A dual path integral representation I. Roditi Centro Brasileiro de Pesquisas Físicas February 20, 2009 1 I. Roditi (CBPF) Collaborators: 2 I. Roditi (CBPF) Collaborators: C. Cappa Ttira 2 I. Roditi (CBPF)
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
More informationHartree-Fock-Roothan Self-Consistent Field Method
Hartree-Fock-Roothan Self-Consistent Field Method 1. Helium Here is a summary of the derivation of the Hartree-Fock equations presented in class. First consider the ground state of He and start with with
More informationReview for Exam Find all a for which the following linear system has no solutions, one solution, and infinitely many solutions.
Review for Exam. Find all a for which the following linear system has no solutions, one solution, and infinitely many solutions. x + y z = 2 x + 2y + z = 3 x + y + (a 2 5)z = a 2 The augmented matrix for
More informationChapter 1: Systems of Linear Equations and Matrices
: Systems of Linear Equations and Matrices Multiple Choice Questions. Which of the following equations is linear? (A) x + 3x 3 + 4x 4 3 = 5 (B) 3x x + x 3 = 5 (C) 5x + 5 x x 3 = x + cos (x ) + 4x 3 = 7.
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationActivity: Derive a matrix from input-output pairs
Activity: Derive a matrix from input-output pairs [ ] a b The 2 2 matrix A = satisfies the following equations: c d Calculate the entries of the matrix. [ ] [ ] a b 5 c d 10 [ ] [ ] a b 2 c d 1 = = [ ]
More informationLecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel
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 information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationComplexity of the quantum adiabatic algorithm
Complexity of the quantum adiabatic algorithm Peter Young e-mail:peter@physics.ucsc.edu Collaborators: S. Knysh and V. N. Smelyanskiy Colloquium at Princeton, September 24, 2009 p.1 Introduction What is
More informationTopics in linear algebra
Chapter 6 Topics in linear algebra 6.1 Change of basis I want to remind you of one of the basic ideas in linear algebra: change of basis. Let F be a field, V and W be finite dimensional vector spaces over
More informationLecture 1 Systems of Linear Equations and Matrices
Lecture 1 Systems of Linear Equations and Matrices Math 19620 Outline of Course Linear Equations and Matrices Linear Transformations, Inverses Bases, Linear Independence, Subspaces Abstract Vector Spaces
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationExam 2 Solutions. (a) Is W closed under addition? Why or why not? W is not closed under addition. For example,
Exam 2 Solutions. Let V be the set of pairs of real numbers (x, y). Define the following operations on V : (x, y) (x, y ) = (x + x, xx + yy ) r (x, y) = (rx, y) Check if V together with and satisfy properties
More informationVolume in n Dimensions
Volume in n Dimensions MA 305 Kurt Bryan Introduction You ve seen that if we have two vectors v and w in two dimensions then the area spanned by these vectors can be computed as v w = v 1 w 2 v 2 w 1 (where
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationRegularization Physics 230A, Spring 2007, Hitoshi Murayama
Regularization Physics 3A, Spring 7, Hitoshi Murayama Introduction In quantum field theories, we encounter many apparent divergences. Of course all physical quantities are finite, and therefore divergences
More informationMeron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations
Meron-Cluster and Nested Cluster Algorithms: Addressing the Sign Problem in Quantum Monte Carlo Simulations Uwe-Jens Wiese Bern University IPAM Workshop QS2009, January 26, 2009 Collaborators: B. B. Beard
More informationLecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 10 The Dirac equation WS2010/11: Introduction to Nuclear and Particle Physics The Dirac equation The Dirac equation is a relativistic quantum mechanical wave equation formulated by British physicist
More informationBASIC NOTIONS. x + y = 1 3, 3x 5y + z = A + 3B,C + 2D, DC are not defined. A + C =
CHAPTER I BASIC NOTIONS (a) 8666 and 8833 (b) a =6,a =4 will work in the first case, but there are no possible such weightings to produce the second case, since Student and Student 3 have to end up with
More informationIterative real-time path integral approach to nonequilibrium quantum transport
Iterative real-time path integral approach to nonequilibrium quantum transport Michael Thorwart Institut für Theoretische Physik Heinrich-Heine-Universität Düsseldorf funded by the SPP 1243 Quantum Transport
More informationContinuous-time quantum Monte Carlo algorithms for impurity problems
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
More informationA Whirlwind Introduction to Coupled Cluster Response Theory. 1 Time-Independent Coupled Cluster Theory
A Whirlwind Introduction to Coupled Cluster Response Theory T. Daniel Crawford, Virginia Tech, Blacksburg, Virginia, U.S.A. 1 Time-Independent Coupled Cluster Theory If the Hamiltonian is independent of
More informationIs a system of fermions in the crossover BCS-BEC. BEC regime a new type of superfluid?
Is a system of fermions in the crossover BCS-BEC BEC regime a new type of superfluid? Finite temperature properties of a Fermi gas in the unitary regime Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski
More informationLECTURES 14/15: LINEAR INDEPENDENCE AND BASES
LECTURES 14/15: LINEAR INDEPENDENCE AND BASES MA1111: LINEAR ALGEBRA I, MICHAELMAS 2016 1. Linear Independence We have seen in examples of span sets of vectors that sometimes adding additional vectors
More informationApplied Linear Algebra
Applied Linear Algebra Gábor P. Nagy and Viktor Vígh University of Szeged Bolyai Institute Winter 2014 1 / 262 Table of contents I 1 Introduction, review Complex numbers Vectors and matrices Determinants
More information