The advent of computer era has opened the possibility to perform large scale
|
|
- Alaina Dennis
- 5 years ago
- Views:
Transcription
1 Chapter 2 Density Matrix Renormalization Group 2.1 Introduction The advent of computer era has opened the possibility to perform large scale numerical simulations of the quantum many-body systems and thus revealing completely new perspectives in the field of condensed matter theory [100]. Indeed, together with the analytic approaches, numerical techniques provide lot of information and details, otherwise inaccessible. However, the simulation of a quantum mechanical system is a very difficult task. It is because of the number of parameters required to represent the quantum state is large. This value usually grow exponentially with the number of constituents of the system due to the corresponding exponential growth of its Hilbert space. This exponential scaling drastically reduces the possibility of a direct simulation of many-body quantum systems. In order to overcome this 12
2 13 limitation many numerical tools have been developed, which include Monte Carlo techniques [101], efficient Hamiltonian digitalization methods such as Lanczos [102] and Davidson [103] procedures and real space renormalization group method [104]. Each of these methods has its own limitations, for example sign problem for fermion models in Monte Carlo, small system size in exact digitalization method, large error in real space renormalization group methods. The Density Matrix Renormalization Group (DMRG) was developed by S. R. White in 1992 [68, 69, 105] to remove some of these limitations and since then it has been proved to be a very powerful method for low dimensional interacting systems. Its remarkable accuracy can be seen in the spin-1 Heisenberg chain; where for a system of hundreds of sites, a precision of for the ground state energy has been achieved [106]. Since then it has been applied to a great variety of systems and problems including spin chains and ladders [107], fermionic [108, 109] and bosonic [34, 75, 110, 111] systems, disordered models [32, 67], impurities [112, 113], molecules [114], to list a few. The details of these can be seen in the recent review by Schollwock [70] and by Hallberg [71]. The global idea of DMRG may be stated in a single sentence: When splitting a system into blocks, these should not be isolated [115]. White [69] realized that, in a real space renormalization group calculation, the lowest energy states of the isolated blocks need not be the best building
3 14 bricks to construct the big block. It is necessary for the block to be related to its environment so as it may choose the best states. Density Matrix Renormalization Group method provides the best method of choosing these best states when going from smaller block to bigger block. This method has since developed into one of the most accurate methods for obtaining not only the ground-state properties and but also the low-lying excited states of onedimensional (d = 1) quantum systems [70, 116, 117]. Extensions to higher dimensions [118] and finite temperature are active areas of research [119]. DMRG algorithms can be classified into two categories viz. Infinite system DMRG and Finite Size DMRG (FS-DMRG). Infinite DMRG can be applied for the quantum system when there is no disorder in the system for which FS-DMRG algorithm has been proved to be better. For the sake of completeness we describe these methods below. The following part of this chapter is organized in the following manner. In the Section (2.2) we briefly discuss the real space renormalization group method and its limitations. Section (2.3) discusses about methods of overcoming these limitations. The details of the Infinite Size Density Matrix Renormalization Method and the Finite Size Density Matrix method are given in Sections (2.4) and (2.5) respectively. The Section(2.6) we describe the physical quantities which are obtained from DMRG method described above and apply this method to pure Bose-Hubbard model.
4 2.2 Real-Space Renormalization Group Method 15 Historically Density Matrix Renormalization Group (DMRG) method has its origin in the analysis of the failure of the real-space renormalization group (RSRG) methods to yield quantitatively acceptable results for the low-energy properties of quantum many-body problems [69]. Since most of the DMRG algorithm can be formulated in the standard real space renormalization group language, we briefly discuss the same below. The tensorial structure of the Hilbert space of a composite system leads to an exponential growth of its dimension and so is the resources needed for the simulation with the number of the system constituents. However, if one is interested in the ground state properties of a one-dimensional system, the number of parameters is limited for noncritical systems or grows polynomially for a critical one. This implies that it is possible to rewrite the state of the system in a more efficient way, i.e., it can be described by using a number of coefficients which is much smaller than the dimension of Hilbert space. Equivalently, a strategy to simulate ground state properties of a system is to consider only a relevant subset of states of the full Hilbert space. This idea is reminiscent of the renormalization group (RG) introduced by Wilson [104]. In the real space renormalization group procedure one typically begins
5 with a small part of a quantum system, called a block, (represented by B, of size L, living on an M - dimensional Hilbert space), and a Hamiltonian which describes the interaction between two identical blocks. Then one projects the composite 2-block system (of size 2L) representation (whose dimension M2) onto the subspace spanned by the M lowest energy eigenstates, thus obtaining a new M-dimensional basis. This procedure is then iteratively repeated, until the desired system size is reached. Below we list the iterative steps of RSRG method Standard Real-Space Renormalization Group algorithm for a 1D quantum system 16 The iterative procedure of the standard real-space renormalization group algorithm contains the following steps: Step 1. Consider a block B of length L (usually L = 2 for the first iteration). Let the basis of this block be represented by WO with p 1,2,, M, where M is the dimension of the Hilbert space of the block B. For example, by taking L = 2, M = 4 and 16, respectively, for spin 1/2 and Hubbard Hamiltonian. We denote B I with represents a single site.
6 17 Step 2. Consider two such blocks B+B' as shown below and its Hamiltonian is given by HBB' = HB HB/ + Interaction between B and B'. (2.1) The HBB' matrix has a dimension of M 2 x M2. BB' Step 3. Diagonalize HBB' and obtain the M eigenvectors of the lowest eigenvalues. Let us represent them by la > with a = 1,2,3,, M 2. i.e. la) = E (2.2) /fir Step 4. Change basis of block B + B' from Iµ/, pr) to la), keeping only the lowest M states, and corresponding Hamiltonian become Hgnew =OHBB/0 1- Step 5. Replace B with Bnew. Thus after one iteration I,u/) is replaced by la) and truncated the Hilbert's space from M 2 to M. Thus BB' B Step 6. Go to Step 1 and continue with the iteration. Renormalization group was successfully applied for the Kondo prob- lem [104], but fails in the description of the strongly interacting systems.
7 18 This failure is due to the procedure followed to increase the system size and to the criterion used to select the representative states of the renormalized block [69, 104]: indeed decimation procedure of the Hilbert space is based on the assumption that the ground state of the entire system will essentially be composed of energetically low-lying states living on smaller subsystems (the forming blocks) which is not always true. A simple counter-example is given by a free particle in a box: the ground state of the system with length 2L has no nodes, whereas any combination of states of two boxes with length L will have node in the middle, thus resulting in higher energy. A convenient strategy to solve this RG breakdown is the following: before choosing the states to be retained for a finite-size block, it is first embedded in some environment that mimics the thermodynamic limit of the system. This is the new key ingredient of the DMRG algorithm. As discussed below, system growth is then slowed down with the number of the algorithm's iterations: from the exponentially fast growth Wilson's procedure to the DMRG linear growth [69, 105]. 2.3 Density-Matrix Approach The Density Matrix Renormalization Group algorithm receives its name from one of the tools, i.e., density matrix, incorporated into it. It is arguable whether the density matrix is the key ingredient of the method, but history
8 19 has made to be considered to be so. The reason for which the density matrix appears-is.the necessity of "fitting": it is necessary to find out which are the block states which reproduce most accurately a chosen global state, which shall be termed the target state. According to some authors, density matrices are the most fundamental way to describe quantum mechanical systems. According to Feynman [121] density matrices are necessary because systems must always be separated from their environment in order to be analyzed. This comment was the one which made White [69, 105] to use them as the foundation of his RG formulation. The superblock method forms the basis of density matrix approach suggested by White and Noack [69]. In the superblock method, one diagonalizes a larger system composed of three or more blocks which includes the two blocks, say BB' which are used to form B,,. The wave function for the superblock are projected onto blocks BB', and these projected states of BB' are kept for the next iteration. For a single particle wave function, this projection is single valued and trivial. The superblock method works quite well in a single particle model [69], with the accuracy increasing rapidly with the number of extra blocks used. However for a many-particle wave function, the "projection" of a wave function onto BB' is many valued, and, in fact, a single many-particle state for the entire lattice generally "projects" onto a complete set of block states. However some of these states are most
9 20 important than others; the density matrix tells us which states are the most important. To understand why density matrix is used to choose the states which are to be kept, consider the following argument by analogy. For an isolated block at finite temperature, the probability that the block is in eigenstates a of the block Hamiltonian is 'proportional to its Boltzmann weight exp(- 13Ea ), 13 = 1/k BT, kb being the Boltzmann's constant. The Boltzmann weight is an eigenvalue of the density matrix exp(--ohb ), and an eigenstates of the Hamiltonian is also an eigenstates of the density matrix. Since lowest energy corresponds to highest probability in the Boltzmann's weight, we can view the standard RG approach as choosing the M most probable eigenstates to represent the block given the assumption that the block is isolated. However in reality the block is not isolated, the density matrix is not exp( OHB), and eigenstates of the block Hamiltonian are not eigenstates of the block's density matrix. For a system which is strongly coupled to the outside universe, it is much more appropriate to use the eigenstates of the density matrix to describe the system rather than the eigenstates of the system's Hamiltonian. Thus a natural generalization of the standard approach is to choose to keep the M most probable eigenstates of the block density matrix. In particular, it can be shown that keeping the most probable eigenstates of the density matrix gives the most accurate representation of the state of
10 21 the system as a whole, i.e. the block plus rest of the lattice [68). Assume that a superblock is diagonalized and one particular state, known as target state, kb) is obtained. Let li), with i = 1, 2,, I be a complete set of states of BB' (the system) and 1j), j = 1, 2,, J be the states of the rest of the superblock, i.e.,the "universe". So kb) = E (2.3) For simplicity assume that oij is real. To define a procedure for producing a set of states of the system lua), a = 1,..., m, with lua) = Ei u'li), which are optimal for representing exactly if I > m. An accurate expansion for > can be constructed in the following form, >= > 1./ > (2.4) a,j In other words one has to minimize the following term, s > - 17-T (2.5) by varying over all ct,,,,j and lua >, subject to < ualua' >= Without loss of generality one can write, >= E aale > Iv- > (2.6) a where v7 =< jlva >= Naac,j, with N chosen to set Ei ItT1 2 = 1. In matrix
11 22 notation one can write, s = (02, - E aau7v7) a=1 2 (2.7) and S to be minimized over all eta, va, and a a, given the specified value of m. The solution to this minimization problem is known from linear algebra [68]. By considering as a rectangular matrix, the solution is produced by the singular value decomposition [68, 120] of 0, = UDV T (2.8) where U and D are I x I matrices, V is an I x J matrix (where j = 1, J, and J > 1), U is orthogonal, V is column orthogonal and the diagonal matrix D contains the singular values of 0. ua, va and aa which minimize S are given as follows; the m largest-magnitude diagonal elements of D are the aa, and the corresponding columns of U and V are the u and v. These optimal states ua are also eigenvectors of the reduced density matrix of the system as part of the universe. This reduced density matrix for the system depends on the state of the universe, which in this case is a pure state 10 >. The density matrix for system is given by, pii = EY i lkiej (2.9) These gives, p = UD2UT (2.10)
12 23 i.e. U diagonalizes p. The eigenvalues of p are wa = 4, and the optimal states u are the eigenstates of p with the largest eigenvalues. Each w a represents the probability of the system being in the state ua, with Ea wa = 1. The deviation of P, = Erna_ l wa from unity measures the accuracy of the truncation to m states Density-Matrix Algorithms A density-matrix algorithm is defined mainly by the form of the superblock and the manner in which the blocks are enlarged (such as by doubling the block, B' = BB, or by adding a single site, B' = B + *), and by the choice of superblock eigenstates i.e. target states, used in constructing the density matrix. An eigenstate of the superblock Hamiltonian used for the determination of density matrix is called the target state. The most efficient algorithms use only a single target state in constructing the density matrix. By targeting only one state, the block states are more specialized for representing that state. The important characteristic of a density matrix algorithm is that the accuracy increases with the number of states kept in the density matrix m. The coefficient governing the increase of accuracy with m is largest with a single target state used for the construction of the density matrix. Superblock can be constructed in different ways. Generally it is enlarged by adding two sites, rather than by adding a block as done in standard R,SRG
13 approach (See Section (2.2)). We construct superblock for a system of L sites with open boundary as shown below; r L Here Be represents left block composed of a - 1 sites. Similarly L represents right block composed of 4 1 sites. The represents a single site. In every iteration of the DMRG algorithm, the left (right) block is up- dated with adding a new site to it, thus increasing its size by one. i.e. BQ and L L _ 1 Br L The standard DMRG algorithm can be used in two different ways viz infinite size method and finite size method, which are briefly discussed below. 2.4 Infinite size density-matrix renormalization group method In the first step of the infinite size density-matrix renormalization group method, one starts with a chain of four sites. The superblock configuration will then be Bi. Here 141. and.13;. both represent a single site. We construct the Hamiltonian for this superblock and diagonalize it to obtain the energy and the wave function of the target state. The Davidson algorithm [103] is used for the diagonalization of the hamiltonian matrix
14 which is highly sparse. Using this target state, the density matrix for the new left (right) block is obtained. The new left and right block is, respectively 25.B1 and 41: The density matrix for the left (right) block is then diagonalized to obtain its eigenvalues and the eigenstates. The new left (right) block is represented by the eigenstates corresponds to largest m eigenvalues of the density matrix of the left (right) block. This is the main idea behind the density matrix renormalization group method. The effective hamiltonian is formed for the new left and right blocks, which has now two sites each. With this new left and right blocks, the new superblock, with length L = 6 is B2 13?. This iteration is continued till desired length is achieved. At each iteration step, both blocks size is increased by one site and the total length of the chain is increased by 2 sites. The infinite size method is used to find ground state properties of infinite chain of pure system. Each step of the iteration pushes the ends of the chain farther from the two sites at the center. After many steps, each block approximately represents one-half of an infinite chain.
15 2.4.1 Graphical presentation of Infinite Size DMRG Method 26 I I rol L= 4 L= 6. I L = 8 (continued till the required length is reached) Algorithm For Infinite size DMRG method 1. represent a site. The basis of single site is represented by la) with a = 1, 2,, n where n is the number of states per site. 2. BL 1-1 and B, represents left and right block respectively. The basis of the left and the right blocks are represented by IN) and Iµ,.) 7.) respec- tively. Here {pi, p,r} = 1, 2,, m where m is the number of states in the left and right blocks. In the first iteration, L = 4 and the left and right blocks have one site each. i.e. m = n in the first iteration. 1, Br 1 B, So the 3. Consider the superblock configuration L _i.b.1 L _ BT 1. The
16 27 basis states for this superbiock is written as jilt at O r >, where at and a, represent the basis for the single sites close to left and right block respectively. The dimension of the Hilbert space of this superblock is m 2 x n2. However, considering fixed density of bosons in the lattice, this dimension can be further reduced. Only states with fixed density of bosons needed to be considered because the number operator commutes with the Hamiltonian. 4. The Hamiltonian of the superbiock H = HLB H LS HRS + HRB HLB,LS HLS,RS HRS,RB Here first four terms represent Hamiltonian of left block(lb), left single site(ls), right single site(rs) and right block(rb) respectively. The last three terms represent the interaction between LB-LS, LS-RS and RS-RB respectively. 5. The Hamiltonian matrix Hij is highly sparse. Diagonalize this Hamiltonian to obtain the energy EL(N) and wavefunction RoLN) of the target state. Note that L is the length of the system and N is the number of particle in this system. Thus WOLN) = Ecili >= E al ar (2.11) ar
17 6. Construct the Density Matrix of new left block which is represented by 28 B 2 whose states can be represented by lkitat) The density matrix for the left block is calculated using Eq. (2.11). ni c*c la la r 11 r. Crr (2.12) Similarly the density matrix for the new right block, represented by Br* -1, is given by P(artir)1(Trvr) Cillatarttr CiltatTrur Mat ( 2. 13) 7. Diagonalize p( No., ),(,0.1) to obtain its eigenvalues co c, 1 k and eigenstates laik ) 1,4 E oatati,a/0-0, (2.14) Mat where k = 1, 2,, m x n. Note that the dimension of the new left block is m x n. Similarly we diagonalize the density matrix for the right block p'(.0., A, ),(7.,,,,, ) to obtain its eigenvalues war k and eigenstates 1c4) lark E > (2.15) anur 8. The eigenvalues of density matrix obey the condition 7 a t(r) Wat(r) = 1 and (.4.) l(r) a > 0. k 9. The new left block is represented by the eigenstates of the density matrix lalk ) with k = 1, 2, m corresponding to m largest c.o c, / k. In other
18 29 words we truncate the new left block with states m x n to m states. This truncation procedure is the most important part of the DMRG algorithm. Similarly we represent the new right block by 14) with k = 1, 2, m corresponding to m largest wok. 10. Go to step 3. and continue with the iteration. 2.5 Finite-Size Density Matrix Renormalization Group method The Infinite Size Density-Matrix Renormalization Group (DMRG) method has proven to be very useful in the studies of pure one-dimensional quantum systems [70, 71]. Open boundary conditions are normally preferred for such calculations since the loss of accuracy with increasing system size is much less than system with periodic boundary conditions. However, this approach has many serious draw back applying it to disorder system. In the infinitesystem DMRG method outlined above the left and right-block bases are not optimized in the following sense: The DMRG estimate for the target-state energy, at the step when the length of the system is L, is not as close to the exact value of the target-state energy for this system size as it can be. It has been found that the FS-DMRG method overcomes this problem [69, 68, 75]. In addition if one is interested in the scaling behavior of various quantities, it is necessary to use Finite-Size Density-Matrix Renormalization Group(FS-
19 DMRG) as discussed below [75]. The conventional FS-DMRG method consists of the following two steps: The infinite-system density-matrix renormalization group method (DMRG), where we start with a system with four sites, add two sites at each step of the iteration, and continue till we obtain a system with the desired number L of sites as discussed in the Sec. (2.4). 2. The finite-system method in which the system size L is held fixed, but the energy of a target state is improved iteratively by a sweeping procedure, described below, till convergence is obtained. First we use the infinite-system DMRG iterations to build up the system to a certain desired size, say L. The L-site superblock configuration is now given by 13, -1 B7* -1. In the next step of the FS-DMRG method, the superblock L configuration E17.2, which clearly keeps the system size fixed at L, is used. This step is called sweeping in the right direction since it increases (decreases) the size of the left (right) block by one site. For this superblock the system is L 1, the universe is B 2, the associated density matrix can be found, and from its most significant states the new effective Hamiltonian for the left block, with (2 + 1) sites, is obtained. We sweep again, in this way, to obtain a left block with (1 + 2) sites and so on till the left block has (L 3) sites and the right block has 1 site so that, along
20 31 with the two sites in between these blocks, the system still has size L; or, if a pre-assigned convergence criterion for the target-state energy is satisfied, this sweeping can be terminated earlier. Note that, in these sweeping steps, for the right block we need to BL-3, which we have already obtained in earlier steps of the infinite-system DMRG. Next we sweep leftward: the size of the left (right) block decreases (increases) by one site at each step. Furthermore, in each of the right- and left- sweeping steps, the energy of the target state decreases systematically till it converges (we use a six-figure convergence criterion in our calculations) Graphical presentation of Finite Size DMRG Method Let us assume we want system size L = 12. First using infinite size DMRG, build a system of L = 12 as represented below. Infinite size DMRG I a I I I
21 After reaching the desired system size (in this case L = 12), we begin finite size sweeping as given below. 32 I I I Finite size sweeping I I I I I I I I I I I 2.6 Method of calculation In the last two sections we had described the algorithm of both infinite size DMRG and FS-DMRG methods. This session we apply this method to the pure Bose-Hubbard model given by
22 33 = t E (alai + Eni (ni <i, j> i 1). (2.16) The number of possible states per site n in the case of Bose- Hubbard model is infinite since there can be any number of bosons on given a site. In a practical DMRG calculation, however, we must restrict the number n = nmax. The smaller the interaction parameters U in the model (2.16), the larger must be nmax. We have taken nmax = 4 and is sufficient for the values of U considered here; we have checked in representative cases that our results do not change significantly if nmax = 5. Using the FS-DMRG method described in the previous sections (2.4) and (2.5) the ground state energy EL(N) and its ground state wave functions i'/poln) of model (2.16) on a chain of length L having N bosons are obtained. From them we calculate a number physically relevant quantities as given below. These quantities are used to distinguish different phases possible in the model (2.16). The single particle energy gap GL for a system of size L is defined as = (2.17) where At (pi) is the energy required to add (remove) one boson for a system of length L; i.e., = EL(N + 1) EL (N) (2.18)
23 = EL(N) E L (N 1). (2.19) 34 The single particle gap of model (2.16) for two values of on-site interactions, U = 2 and 8 and for density p = NIL = 1 is given in Fig It is obvious from this figure that the single particle energy gap vanishes in the thermodynamic limit (L > oo) for U = 2 and remain finite for U = 8. This implies a phase transition from gapless superfluid to finite gaped Mott insulator phase as on-site interaction is increased /L 0.25 Figure 2.1: Single particle energy gap CI, as a function of 1/L of model (2.16) for U = 2 and 8 with density p = 1. The critical value of on-site interaction UG, for the SF to. MI transition is best obtained from finite size analysis of single particle energy gap [75] where plots of LG L versus U for different values of length L will coalesce in
24 the superfluid region of the phase diagram. Such a figure is given in Fig. 2.2 and the critical on-site interaction UG, fz L= fr o- ********* I 5 I I 8 I U Figure 2.2: Finite size scaling of single particle energy gap LGL versus U for different lengths showing SF to MI quantum phase transition. The plots of LG L coalesce below U. < Uc 3.4. The correlation function that characterizes the SF phase is rs,f(r) (PoLNI atoarlooln), (2.20) where I b, al,n) is the ground-state wavefunction of the system with size L and N bosons. The associated correlation length can be obtained from the second moment of this correlation function, namely, (sf = Er r2r1f(r)] 1/2 (r) (2.21)
25 36 Note that (f,f is the correlation length for SF ordering in a system of size L; it remains finite so long as L < oo. In the superfluid phase the correlation function II,F(r) shows power law decay when plotted against r which is as follows, risf 1 r -K/2 (2.22) where K is Luttinger Liquid parameter [62, 63]. For onsite interaction U = 2 the power law decay of the SF correlation is given in Fig Plots similar to this is used to obtain the Luttinger Liquid parameter K and the variation of K as a function of onsite interaction is given in Fig The Bosonization study of model (2.16) [62, 63] predict the SF to MI transition at K = 1/2 which has been verified in Fig Another way of showing the SF-MI transition in the model (2.16) is from the study of finite size scaling of correlation length defined by Eq. (2.21). The SF correlation function decay exponentially in the Mott insulator phase due to the existence of finite gap in its energy spectrum and the correlation length 1/GL. Thus finite size scaling of correlation length behaves similar to that of gap and such plots are given in Fig. 2.5 showing the coalescence of L/(15! for different lengths L below U = Uc =, 3.4 which is consistent with finite size scaling of LG L as given in Fig Other physical quantities which we calculate include the expectation value
26 r Figure 2.3: Variation of SF correlation function rsf(r) with respect to r for U = 2. The continuous line is the fit r s-e(r) ti r-1 /2 which gives the LL parameter K = 0.33 ± N E. aims N U Figure 2.4: Variation of Luttinger Liquid parameter K with respect to the onsite interaction U for model (2.16).
27 L= U Figure 2.5: Finite size scaling of correlation length LAE versus U for different lengths showing SF to MI quantum phase transition. The plots of Lgi, coalesce below U < Uc ti 3.4. of number operator (n i) which measures the local density and is defined as, (ni) = (,bolninilooln) (2.23) Using this expectation value, we calculate CDW order parameter OMDW = - E(-1)z[(n%) - (2.24) This order parameter is finite in density wave phase but zero in superfluid and Mott insulator phases. We can also distinguish superfluid phase from Mott insulator or charge density wave phases using compressibility defined as change in the density
28 39 with respect to change in the chemical potential, i.e., 8n (2.25) where the chemical potential is define as the change in the ground state energy with respect to change in the total number of bosons, the is obtained from the relation SE 5N (2.26) U= Pa OOOOOOOOOOOOOOOOOOO OOOOOO OOOOOOOOO Figure 2.6: Boson density p and compressibility is; as a function of chemical potential,u for U = 2. In Fig. 2.6 density p and compressibility lc are plotted as function of chemical potential for onsite interaction U = 2. Finite value of compressibility for all values of 1.1, for this value of onsite interaction suggest that the
29 U= , Figure 2.7: Boson density p and compressibility is as a function of chemical potentialµ for U = 8. system is always in the superfluid phase. However, as we increase the onsite interaction U > Uc, Mott insulator phase exist for some values of chemical potential where the compressibility vanishes. For example in Fig. 2.7 we have plotted both p and n for various values of chemical potential,u. It is obvious that for 1.45 < µ < 4.23, the compressibility is zero and corresponding value of density p = 1. This region has density pinned to one. The value ofµ at upper and lower knees, respectively, are equal to A+ and bi-, the energy require for adding and removing one particle from the ground state (see Eq From the calculation of p+ and tr we obtain the phase diagram of pure Bose-Hubbard model in the µ U plane and is given in Fig. 2.8(a).
30 41 The Mott insulator phase has finite gap and is given in Fig. 2.8(b). 7 6: 5: 3: 2: 1: 0 (a) 4 8 SF MI U 10 Figure 2.8: (a) Phase diagram of Bose-Hubbard model (2.16) for p = 1 and (b) the variation of gap G = 11+,u in the Mott insulator phase.
In this chapter we focus our attention on the model (3.1) in the hard-core. the model (3.1) can be mapped into spin-2 XX model in the presence of
Chapter 4 Hard core Bosons 4.1 Introduction In this chapter we focus our attention on the model (3.1) in the hard-core limit, i.e., taking U = oo. In this limit with density of bosons p = 1/2 the model
More informationBose-Hubbard model with random and Fibonacci potentials
Chapter 6 Bose-Hubbard model with random and Fibonacci potentials 6.1 Introduction Experimentalists are able to produce random and pseudo-random optical potentials in addition to periodic potentials by
More informationThe density matrix renormalization group and tensor network methods
The density matrix renormalization group and tensor network methods Outline Steve White Exploiting the low entanglement of ground states Matrix product states and DMRG 1D 2D Tensor network states Some
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 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 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 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 informationQuantum spin systems - models and computational methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Quantum spin systems - models and computational methods Anders W. Sandvik, Boston University Lecture outline Introduction
More informationR. Citro. In collaboration with: A. Minguzzi (LPMMC, Grenoble, France) E. Orignac (ENS, Lyon, France), X. Deng & L. Santos (MP, Hannover, Germany)
Phase Diagram of interacting Bose gases in one-dimensional disordered optical lattices R. Citro In collaboration with: A. Minguzzi (LPMMC, Grenoble, France) E. Orignac (ENS, Lyon, France), X. Deng & L.
More informationLow-lying excitation spectrum of quantum many-body systems
PHYSICAL REVIEW B VOLUME 54, NUMBER 21 1 DECEMBER 1996-I Low-lying excitation spectrum of quantum many-body systems P. Monthoux Department of Physics and National High Magnetic Field Laboratory, Florida
More informationQuantum Lattice Models & Introduction to Exact Diagonalization
Quantum Lattice Models & Introduction to Exact Diagonalization H! = E! Andreas Läuchli IRRMA EPF Lausanne ALPS User Workshop CSCS Manno, 28/9/2004 Outline of this lecture: Quantum Lattice Models Lattices
More informationIntroduction to the Mathematics of the XY -Spin Chain
Introduction to the Mathematics of the XY -Spin Chain Günter Stolz June 9, 2014 Abstract In the following we present an introduction to the mathematical theory of the XY spin chain. The importance of this
More information21 Density-Matrix Renormalization Group Algorithms
21 Density-Matrix Renormalization Group Algorithms Eric Jeckelmann Institut für Theoretische Physik, Leibniz Universität Hannover, 30167 Hannover, Germany In this chapter I will introduce the basic Density
More informationReal-Space Renormalization Group (RSRG) Approach to Quantum Spin Lattice Systems
WDS'11 Proceedings of Contributed Papers, Part III, 49 54, 011. ISBN 978-80-7378-186-6 MATFYZPRESS Real-Space Renormalization Group (RSRG) Approach to Quantum Spin Lattice Systems A. S. Serov and G. V.
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 informationDensity Matrix Renormalization: A Review of the Method and its Applications arxiv:cond-mat/ v1 26 Mar 2003
Density Matrix Renormalization: A Review of the Method and its Applications arxiv:cond-mat/0303557v1 26 Mar 2003 Karen Hallberg Centro Atómico Bariloche and Instituto Balseiro 8400 Bariloche, Argentina
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 informationMomentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model
Momentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model Örs Legeza Reinhard M. Noack Collaborators Georg Ehlers Jeno Sólyom Gergely Barcza Steven R. White Collaborators Georg Ehlers
More informationSpin liquids on ladders and in 2d
Spin liquids on ladders and in 2d MPA Fisher (with O. Motrunich) Minnesota, FTPI, 5/3/08 Interest: Quantum Spin liquid phases of 2d Mott insulators Background: Three classes of 2d Spin liquids a) Topological
More information4 Matrix product states
Physics 3b Lecture 5 Caltech, 05//7 4 Matrix product states Matrix product state (MPS) is a highly useful tool in the study of interacting quantum systems in one dimension, both analytically and numerically.
More informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More informationNON EQUILIBRIUM DYNAMICS OF QUANTUM ISING CHAINS IN THE PRESENCE OF TRANSVERSE AND LONGITUDINAL MAGNETIC FIELDS
NON EQUILIBRIUM DYNAMICS OF QUANTUM ISING CHAINS IN THE PRESENCE OF TRANSVERSE AND LONGITUDINAL MAGNETIC FIELDS by Zahra Mokhtari THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE
More informationNumerical Renormalization Group Methods
Numerical Renormalization Group Methods Martin Andersson 1 Licenciate thesis Institute of Theoretical Physics Chalmers University of Technology and Göteborg University S-412 96 Göteborg, Sweden Abstract
More informationPhysics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory 1. Introduction Bound state perturbation theory applies to the bound states of perturbed systems,
More informationJournal Club: Brief Introduction to Tensor Network
Journal Club: Brief Introduction to Tensor Network Wei-Han Hsiao a a The University of Chicago E-mail: weihanhsiao@uchicago.edu Abstract: This note summarizes the talk given on March 8th 2016 which was
More informationH ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org Simulations of
More informationMany-Body Localization. Geoffrey Ji
Many-Body Localization Geoffrey Ji Outline Aside: Quantum thermalization; ETH Single-particle (Anderson) localization Many-body localization Some phenomenology (l-bit model) Numerics & Experiments Thermalization
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 informationQuantum many-body systems and tensor networks: simulation methods and applications
Quantum many-body systems and tensor networks: simulation methods and applications Román Orús School of Physical Sciences, University of Queensland, Brisbane (Australia) Department of Physics and Astronomy,
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 4 Sep 2006
Advances in Physics Vol. 00, No. 00, January-February 2005, 1 54 arxiv:cond-mat/0609039v1 [cond-mat.str-el] 4 Sep 2006 New Trends in Density Matrix Renormalization KAREN A. HALLBERG Instituto Balseiro
More informationQuantum Spin-Metals in Weak Mott Insulators
Quantum Spin-Metals in Weak Mott Insulators MPA Fisher (with O. Motrunich, Donna Sheng, Simon Trebst) Quantum Critical Phenomena conference Toronto 9/27/08 Quantum Spin-metals - spin liquids with Bose
More informationLECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS
LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS Path-integral for lattice bosons*: oriented closed loops, of course LECTURE 3 WORM ALGORITHM
More information1 Quantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More informationNew trends in density matrix renormalization
Advances in Physics, Vol. 55, Nos. 5 6, July October 2006, 477 526 New trends in density matrix renormalization KAREN A. HALLBERG Instituto Balseiro and Centro Ato mico Bariloche, Comisio n Nacional de
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 informationLecture notes for QFT I (662)
Preprint typeset in JHEP style - PAPER VERSION Lecture notes for QFT I (66) Martin Kruczenski Department of Physics, Purdue University, 55 Northwestern Avenue, W. Lafayette, IN 47907-036. E-mail: markru@purdue.edu
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 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 informationQuantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More informationImpurity corrections to the thermodynamics in spin chains using a transfer-matrix DMRG method
PHYSICAL REVIEW B VOLUME 59, NUMBER 9 1 MARCH 1999-I Impurity corrections to the thermodynamics in spin chains using a transfer-matrix DMRG method Stefan Rommer and Sebastian Eggert Institute of Theoretical
More informationNumerical diagonalization studies of quantum spin chains
PY 502, Computational Physics, Fall 2016 Anders W. Sandvik, Boston University Numerical diagonalization studies of quantum spin chains Introduction to computational studies of spin chains Using basis states
More informationThe density-matrix renormalization group
REVIEWS OF MODERN PHYSICS, VOLUME 77, JANUARY 2005 The density-matrix renormalization group U. Schollwöck RWTH Aachen University, D-52056 Aachen, Germany Published 26 April 2005 The density-matrix renormalization
More informationExcursion: MPS & DMRG
Excursion: MPS & DMRG Johannes.Schachenmayer@gmail.com Acronyms for: - Matrix product states - Density matrix renormalization group Numerical methods for simulations of time dynamics of large 1D quantum
More informationTensor network methods in condensed matter physics. ISSP, University of Tokyo, Tsuyoshi Okubo
Tensor network methods in condensed matter physics ISSP, University of Tokyo, Tsuyoshi Okubo Contents Possible target of tensor network methods! Tensor network methods! Tensor network states as ground
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationRenormalization of Tensor Network States
Renormalization of Tensor Network States I. Coarse Graining Tensor Renormalization Tao Xiang Institute of Physics Chinese Academy of Sciences txiang@iphy.ac.cn Numerical Renormalization Group brief introduction
More informationConsistent Histories. Chapter Chain Operators and Weights
Chapter 10 Consistent Histories 10.1 Chain Operators and Weights The previous chapter showed how the Born rule can be used to assign probabilities to a sample space of histories based upon an initial state
More informationPHYS852 Quantum Mechanics II, Spring 2010 HOMEWORK ASSIGNMENT 8: Solutions. Topics covered: hydrogen fine structure
PHYS85 Quantum Mechanics II, Spring HOMEWORK ASSIGNMENT 8: Solutions Topics covered: hydrogen fine structure. [ pts] Let the Hamiltonian H depend on the parameter λ, so that H = H(λ). The eigenstates and
More informationMain matrix factorizations
Main matrix factorizations A P L U P permutation matrix, L lower triangular, U upper triangular Key use: Solve square linear system Ax b. A Q R Q unitary, R upper triangular Key use: Solve square or overdetrmined
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
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 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 informationFully symmetric and non-fractionalized Mott insulators at fractional site-filling
Fully symmetric and non-fractionalized Mott insulators at fractional site-filling Itamar Kimchi University of California, Berkeley EQPCM @ ISSP June 19, 2013 PRL 2013 (kagome), 1207.0498...[PNAS] (honeycomb)
More informationQuantum dynamics in ultracold atoms
Rather don t use Power-Points title Page Use my ypage one instead Quantum dynamics in ultracold atoms Corinna Kollath (Ecole Polytechnique Paris, France) T. Giamarchi (University of Geneva) A. Läuchli
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 informationChapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found.
Chapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found. In applying quantum mechanics to 'real' chemical problems, one is usually faced with a Schrödinger
More informationHeisenberg Antiferromagnet on a Triangular Lattice* ABSTRACT
SLAC-PUB-4880 February 1989 (T) Heisenberg Antiferromagnet on a Triangular Lattice* D. HORN School of Physics and Astronomy Tel-Aviv Tel-Aviv University 69978 Israel and H.R.QuINN AND M.WEINSTEIN Stanford
More informationThermodynamics of quantum Heisenberg spin chains
PHYSICAL REVIEW B VOLUME 58, NUMBER 14 Thermodynamics of quantum Heisenberg spin chains 1 OCTOBER 1998-II Tao Xiang Research Centre in Superconductivity, University of Cambridge, Madingley Road, Cambridge
More informationGapless Spin Liquids in Two Dimensions
Gapless Spin Liquids in Two Dimensions MPA Fisher (with O. Motrunich, Donna Sheng, Matt Block) Boulder Summerschool 7/20/10 Interest Quantum Phases of 2d electrons (spins) with emergent rather than broken
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationIntroduction to tensor network state -- concept and algorithm. Z. Y. Xie ( 谢志远 ) ITP, Beijing
Introduction to tensor network state -- concept and algorithm Z. Y. Xie ( 谢志远 ) 2018.10.29 ITP, Beijing Outline Illusion of complexity of Hilbert space Matrix product state (MPS) as lowly-entangled state
More informationParallelization Strategies for Density Matrix Renormalization Group algorithms on Shared-Memory Systems
Parallelization Strategies for Density Matrix Renormalization Group algorithms on Shared-Memory Systems G. Hager HPC Services, Computing Center Erlangen, Germany E. Jeckelmann Theoretical Physics, Univ.
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, Institute for Theoretical Physics, ETH Zurich Bosons vs Fermions P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009) : fermionic 2D MERA P. Corboz, R. Orus,
More informationJ. Phys.: Condens. Matter 10 (1998) L159 L165. Printed in the UK PII: S (98)90604-X
J. Phys.: Condens. Matter 10 (1998) L159 L165. Printed in the UK PII: S0953-8984(98)90604-X LETTER TO THE EDITOR Calculation of the susceptibility of the S = 1 antiferromagnetic Heisenberg chain with single-ion
More informationEfficient time evolution of one-dimensional quantum systems
Efficient time evolution of one-dimensional quantum systems Frank Pollmann Max-Planck-Institut für komplexer Systeme, Dresden, Germany Sep. 5, 2012 Hsinchu Problems we will address... Finding ground states
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 informationThe 1+1-dimensional Ising model
Chapter 4 The 1+1-dimensional Ising model The 1+1-dimensional Ising model is one of the most important models in statistical mechanics. It is an interacting system, and behaves accordingly. Yet for a variety
More informationSusumu YAMADA 1,3 Toshiyuki IMAMURA 2,3, Masahiko MACHIDA 1,3
Dynamical Variation of Eigenvalue Problems in Density-Matrix Renormalization-Group Code PP12, Feb. 15, 2012 1 Center for Computational Science and e-systems, Japan Atomic Energy Agency 2 The University
More informationQuantum Metropolis Sampling. K. Temme, K. Vollbrecht, T. Osborne, D. Poulin, F. Verstraete arxiv:
Quantum Metropolis Sampling K. Temme, K. Vollbrecht, T. Osborne, D. Poulin, F. Verstraete arxiv: 09110365 Quantum Simulation Principal use of quantum computers: quantum simulation Quantum Chemistry Determination
More informationTime-dependent DMRG:
The time-dependent DMRG and its applications Adrian Feiguin Time-dependent DMRG: ^ ^ ih Ψ( t) = 0 t t [ H ( t) E ] Ψ( )... In a truncated basis: t=3 τ t=4 τ t=5τ t=2 τ t= τ t=0 Hilbert space S.R.White
More informationPhysics 221A Fall 1996 Notes 21 Hyperfine Structure in Hydrogen and Alkali Atoms
Physics 221A Fall 1996 Notes 21 Hyperfine Structure in Hydrogen and Alkali Atoms Hyperfine effects in atomic physics are due to the interaction of the atomic electrons with the electric and magnetic multipole
More informationThe Mott Metal-Insulator Transition
Florian Gebhard The Mott Metal-Insulator Transition Models and Methods With 38 Figures Springer 1. Metal Insulator Transitions 1 1.1 Classification of Metals and Insulators 2 1.1.1 Definition of Metal
More informationQuantum Phases in Bose-Hubbard Models with Spin-orbit Interactions
Quantum Phases in Bose-Hubbard Models with Spin-orbit Interactions Shizhong Zhang The University of Hong Kong Institute for Advanced Study, Tsinghua 24 October 2012 The plan 1. Introduction to Bose-Hubbard
More informationDensity Matrices. Chapter Introduction
Chapter 15 Density Matrices 15.1 Introduction Density matrices are employed in quantum mechanics to give a partial description of a quantum system, one from which certain details have been omitted. For
More informationQuantum Many Body Systems and Tensor Networks
Quantum Many Body Systems and Tensor Networks Aditya Jain International Institute of Information Technology, Hyderabad aditya.jain@research.iiit.ac.in July 30, 2015 Aditya Jain (IIIT-H) Quantum Hamiltonian
More informationNumerical Methods in Many-body Physics
Numerical Methods in Many-body Physics Reinhard M. Noack Philipps-Universität Marburg Exchange Lecture BME, Budapest, Spring 2007 International Research Training Group 790 Electron-Electron Interactions
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
More information221B Lecture Notes Spontaneous Symmetry Breaking
B Lecture Notes Spontaneous Symmetry Breaking Spontaneous Symmetry Breaking Spontaneous Symmetry Breaking is an ubiquitous concept in modern physics, especially in condensed matter and particle physics.
More informationQuantum Monte Carlo Simulations of Exciton Condensates
Quantum Monte Carlo Simulations of Exciton Condensates J. Shumway a and D. M. Ceperley b a Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 8583 b Dept. of Physics, University of Illinois,
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationQuantum simulation with string-bond states: Joining PEPS and Monte Carlo
Quantum simulation with string-bond states: Joining PEPS and Monte Carlo N. Schuch 1, A. Sfondrini 1,2, F. Mezzacapo 1, J. Cerrillo 1,3, M. Wolf 1,4, F. Verstraete 5, I. Cirac 1 1 Max-Planck-Institute
More informationThe end is (not) in sight: exact diagonalization, Lanczos, and DMRG
The end is (not) in sight: exact diagonalization, Lanczos, and DMRG Jürgen Schnack, Matthias Exler, Peter Hage, Frank Hesmer Department of Physics - University of Osnabrück http://www.physik.uni-osnabrueck.de/makrosysteme/
More informationSommerfeld (1920) noted energy levels of Li deduced from spectroscopy looked like H, with slight adjustment of principal quantum number:
Spin. Historical Spectroscopy of Alkali atoms First expt. to suggest need for electron spin: observation of splitting of expected spectral lines for alkali atoms: i.e. expect one line based on analogy
More informationFrustration-free Ground States of Quantum Spin Systems 1
1 Davis, January 19, 2011 Frustration-free Ground States of Quantum Spin Systems 1 Bruno Nachtergaele (UC Davis) based on joint work with Sven Bachmann, Spyridon Michalakis, Robert Sims, and Reinhard Werner
More informationarxiv: v2 [cond-mat.str-el] 3 Jan 2011
The density-matrix renormalization group in the age of matrix product states Ulrich Schollwöck arxiv:1008.3477v2 [cond-mat.str-el] 3 Jan 2011 Department of Physics, Arnold Sommerfeld Center for Theoretical
More informationParallel Singular Value Decomposition. Jiaxing Tan
Parallel Singular Value Decomposition Jiaxing Tan Outline What is SVD? How to calculate SVD? How to parallelize SVD? Future Work What is SVD? Matrix Decomposition Eigen Decomposition A (non-zero) vector
More informationTime-Independent Perturbation Theory
4 Phys46.nb Time-Independent Perturbation Theory.. Overview... General question Assuming that we have a Hamiltonian, H = H + λ H (.) where λ is a very small real number. The eigenstates of the Hamiltonian
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 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 informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More information2015 Summer School on Emergent Phenomena in Quantum Materials. Program Overview
Emergent Phenomena in Quantum Materials Program Overview Each talk to be 45min with 15min Q&A. Monday 8/3 8:00AM Registration & Breakfast 9:00-9:10 Welcoming Remarks 9:10-10:10 Eugene Demler Harvard University
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 informationMany Body Quantum Mechanics
Many Body Quantum Mechanics In this section, we set up the many body formalism for quantum systems. This is useful in any problem involving identical particles. For example, it automatically takes care
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 informationA new perspective on long range SU(N) spin models
A new perspective on long range SU(N) spin models Thomas Quella University of Cologne Workshop on Lie Theory and Mathematical Physics Centre de Recherches Mathématiques (CRM), Montreal Based on work with
More informationNewton s Method and Localization
Newton s Method and Localization Workshop on Analytical Aspects of Mathematical Physics John Imbrie May 30, 2013 Overview Diagonalizing the Hamiltonian is a goal in quantum theory. I would like to discuss
More informationRenormalization of Tensor- Network States Tao Xiang
Renormalization of Tensor- Network States Tao Xiang Institute of Physics/Institute of Theoretical Physics Chinese Academy of Sciences txiang@iphy.ac.cn Physical Background: characteristic energy scales
More informationThe state of a quantum ideal gas is uniquely specified by the occupancy of singleparticle
Ideal Bose gas The state of a quantum ideal gas is uniquely specified by the occupancy of singleparticle states, the set {n α } where α denotes the quantum numbers of a singleparticles state such as k
More information