On high-temperature expansions of the Hubbard model

We discuss the application of the high-temperature expansion method to the Hubbard model. We recalculate the expansion series of the susceptibility up to the sixth order in the transfer matrix element,t, in the strong correlation limit, and up to the fourth order int in case that the repulsive potential,U, is finite, butt/U 1. It is seen that the convergence of the series is very poor.     

Summation of perturbation series in the generalized Hubbard model

The cumulant expansion is proposed for the summation of perturbation series in the generalized Hubbard model, which considers strongly correlated d-electrons hybridized with nearly free conduction electrons. It is shown, that summation of the principle series for the Green's functions leads to the Kondo renormalization of d-level hybridized with s-band.     

High-temperature expansions to fifteenth order

High-temperature series expansions of the susceptibility and second moment to 15th order are calculated for zero external field on the linear chain (LC), plane square (PSQ), simple cubic (SC), and body-centered cubic (BCC) lattices. Checks for specific models against pertinent work in the literature are detailed.     

High-temperature series expansion for the relaxation times of the two dimensional Ising model

We derive the high temperature series expansions for the two relaxation times of the single spin-flip kinetic Ising model on the square lattice. The series for the linear relaxation time _l is obtained with 20 non-trivial terms, and the analysis yields 2.183±0.005 as the value of the critical exponent _l , which is equal to the dynamical critical exponentz in the two-dimensional case. For the non-linear relaxation time we obtain 15 non-trivial terms, and the analysis leads to the results _nl = 2.08 ± 0.07. The scaling relation _l – _nl = ( being the exponent of the order parameter) seems to be fulfilled, though the error bars of _nl are still quite substantial. In addition, we obtain the series expansion of the linear relaxation time on the honeycomb lattice with 22 non-trivial terms. The result for the critical exponent is close to the value obtained on the square lattice, which is expected from universality.     

Systematic series expansions for processes on networks

We use series expansions to study dynamics of equilibrium and nonequilibrium systems on networks. This analytical method enables us to include detailed nonuniversal effects of the network structure. We show that even low order calculations produce results which compare accurately to numerical simulation, while the results can be systematically improved. We show that certain commonly accepted analytical results for the critical point on networks with a broad degree distribution need to be modified in certain cases due to disassortativity; the present method is able to take into account the assortativity at sufficiently high order, while previous results correspond to leading and second order approximations in this method. Finally, we apply this method to real-world data.     

Flow equations for the ionic Hubbard model

Taking the site-diagonal terms of the ionic Hubbard model (IHM) in one and two spatial dimensions, as H₀, we employ Continuous Unitary Transformations (CUT) to obtain a “classical” effective Hamiltonian in which hopping term has been renormalized to zero. For this Hamiltonian spin gap and charge gap are calculated at half-filling and subject to periodic boundary conditions. Our calculations indicate two transition points. In fixed Δ, as U increases from zero, there is a region in which both spin gap and charge gap are positive and identical; characteristic of band insulators. Upon further increasing U, first transition occurs at U=Uc₁, where spin and charge gaps both vanish and remain zero up to U=Uc₂. A gap-less state in charge and spin sectors characterizes a metal. For U>Uc₂ spin gap remains zero and charge gap becomes positive. This third region corresponds to a Mott insulator in which charge excitations are gaped, while spin excitations remain gap-less.     

Random dispersion approximation for the Hubbard model

We use the Random Dispersion Approximation (RDA) to study the Mott-Hubbard transition in the Hubbard model at half band filling. The RDA becomes exact for the Hubbard model in infinite dimensions. We implement the RDA on finite chains and employ the Lanczos exact diagonalization method in real space to calculate the ground-state energy, the average double occupancy, the charge gap, the momentum distribution, and the quasi-particle weight. We find a satisfactory agreement with perturbative results in the weak- and strong-coupling limits. A straightforward extrapolation of the RDA data for L ≤ 14 lattice results in a continuous Mott-Hubbard transition at U_c≈W. We discuss the significance of a possible signature of a coexistence region between insulating and metallic ground states in the RDA that would correspond to the scenario of a discontinuous Mott-Hubbard transition as found in numerical investigations of the Dynamical Mean-Field Theory for the Hubbard model.     

Time-dependent Gutzwiller approximation for the Hubbard model

We develop a time-dependent Gutzwiller approximation (GA) for the Hubbard model analogous to the time-dependent Hartree-Fock (HF) method. This new formalism incorporates ground state correlations of the random phase approximation (RPA) type beyond the GA. Static quantities like ground state energy and double occupancy are in excellent agreement with exact results in one dimension up to moderate coupling and in two dimensions for all couplings. We find a substantial improvement over traditional GA and HF+RPA treatments. Dynamical correlation functions can be computed and are also substantially better than HF+RPA ones and obey well behaved sum rules.     

Ladder operator for the one-dimensional Hubbard model

The one-dimensional Hubbard model is integrable in the sense that it has an infinite family of conserved currents. We explicitly construct a ladder operator which can be used to iteratively generate all of the conserved current operators. This construction is different from that used for Lorentz invariant systems such as the Heisenberg model. The Hubbard model is not Lorentz invariant, due to the separation of spin and charge excitations. The ladder operator is obtained by a very general formalism which is applicable to any model that can be derived from a solution of the Yang-Baxter equation.     

High-temperature expansions in statistical theory of atomic ordering

In the present article a high-temperature series for the statistical sum of a binary alloy is considered. The study is carried out within a statistical ensemble, corresponding to fixing atomic concentration in certain chosen groups of lattice sites. Diagrammatic technique for calculating the statistical sum of the ensemble indicated is introduced.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 2, pp. 104–108, February, 1982.     

Low-temperature series expansions for lattices of non-integral dimensions

Low-temperature series expansions for Ising models on lattices of non-integral dimensions are studied. Critical exponents β for various dimensions are extrapolated from series expansions for the generalized equivalent neighbor lattice.     

High-temperature series analysis of the p-state Potts glass model on d-dimensional hypercubic lattices

We analyze recently extended high-temperature series expansions for the “Edwards-Anderson” spin-glass susceptibility of the p-state Potts glass model on d-dimensional hypercubic lattices for the case of a symmetric bimodal distribution of ferro- and antiferromagnetic nearest-neighbor couplings . In these star-graph expansions up to order 22 in the inverse temperature , the number of Potts states p and the dimension d are kept as free parameters which can take any value. By applying several series analysis techniques to the new series expansions, this enabled us to determine the critical coupling K_c and the critical exponent of the spin-glass susceptibility in a large region of the two-dimensional (p,d)-parameter space. We discuss the thus obtained information with emphasis on the lower and upper critical dimensions of the model and present a careful comparison with previous estimates for special values of p and d. Received: 25 May 1998 / Revised and Accepted: 11 August 1998