Equation of motion approach to the Hubbard model in infinite dimensions

We consider the Hubbard model on the infinite-dimensional Bethe lattice and construct a systematic series of self-consistent approximations to the one-particle Green's function, $G^{(n)}(\omega),\ n=2,3,\dots\$ . The first $n-1$ equations of motion are exactly fullfilled by $G^{(n)}(\omega)$ and the $n$'th equation of motion is decoupled following a simple set of decoupling rules. $G^{(2)}(\omega)$ corresponds to the Hubbard-III approximation. We present analytic and numerical results for the Mott-Hubbard transition at half filling for $n=2,3,4$.Comment: 10pager, REVTEX, 8-figures not available in postscript, manuscript may be understood without figure

Cluster persistence in one-dimensional diffusion--limited cluster--cluster aggregation

The persistence probability, $P_C(t)$, of a cluster to remain unaggregated is studied in cluster-cluster aggregation, when the diffusion coefficient of a cluster depends on its size $s$ as $D(s) \sim s^\gamma$. In the mean-field the problem maps to the survival of three annihilating random walkers with time-dependent noise correlations. For $\gamma \ge 0$ the motion of persistent clusters becomes asymptotically irrelevant and the mean-field theory provides a correct description. For $\gamma < 0$ the spatial fluctuations remain relevant and the persistence probability is overestimated by the random walk theory. The decay of persistence determines the small size tail of the cluster size distribution. For $0 < \gamma < 2$ the distribution is flat and, surprisingly, independent of $\gamma$.Comment: 11 pages, 6 figures, RevTeX4, submitted to Phys. Rev.

Charge-density-wave order parameter of the Falicov-Kimball model in infinite dimensions

In the large-U limit, the Falicov-Kimball model maps onto an effective Ising model, with an order parameter described by a BCS-like mean-field theory in infinite dimensions. In the small-U limit, van Dongen and Vollhardt showed that the order parameter assumes a strange non-BCS-like shape with a sharp reduction near T approx T_c/2. Here we numerically investigate the crossover between these two regimes and qualitatively determine the order parameter for a variety of different values of U. We find the overall behavior of the order parameter as a function of temperature to be quite anomalous.Comment: (5 pages, 3 figures, typeset with ReVTeX4

Many Body Correlation Corrections to Superconducting Pairing in Two Dimensions.

We demonstrate that in the strong coupling limit (the superconducting gap $\Delta$ is as large as the chemical potential $\mu$), which is relevant to the high-$T_c$ superconductivity, the correlation corrections to the gap and critical temperature are about 10\% of the corresponding mean field approximation values. For the weak coupling ($\Delta \ll \mu$) the correlation corrections are very large: of the order of 100\% of the corresponding mean field values.Comment: LaTeX 12 page

Symmetry breaking in the Hubbard model at weak coupling

The phase diagram of the Hubbard model is studied at weak coupling in two and three spatial dimensions. It is shown that the Neel temperature and the order parameter in d=3 are smaller than the Hartree-Fock predictions by a factor of q=0.2599. For d=2 we show that the self-consistent (sc) perturbation series bears no relevance to the behavior of the exact solution of the Hubbard model in the symmetry-broken phase. We also investigate an anisotropic model and show that the coupling between planes is essential for the validity of mean-field-type order parameters

Influence of Hybridization on the Properties of the Spinless Falicov-Kimball Model

Without a hybridization between the localized f- and the conduction (c-) electron states the spinless Falicov-Kimball model (FKM) is exactly solvable in the limit of high spatial dimension, as first shown by Brandt and Mielsch. Here I show that at least for sufficiently small c-f-interaction this exact inhomogeneous ground state is also obtained in Hartree-Fock approximation. With hybridization the model is no longer exactly solvable, but the approximation yields that the inhomogeneous charge-density wave (CDW) ground state remains stable also for finite hybridization V smaller than a critical hybridization V_c, above which no inhomogeneous CDW solution but only a homogeneous solution is obtained. The spinless FKM does not allow for a ''ferroelectric'' ground state with a spontaneous polarization, i.e. there is no nonvanishing $$-expectation value in the limit of vanishing hybridization.Comment: 7 pages, 6 figure

Magnetic properties of the three-dimensional Hubbard model at half filling

We study the magnetic properties of the 3d Hubbard model at half-filling in the TPSC formalism, previously developed for the 2d model. We focus on the N\'eel transition approached from the disordered side and on the paramagnetic phase. We find a very good quantitative agreement with Dynamical Mean-Field results for the isotropic 3d model. Calculations on finite size lattices also provide satisfactory comparisons with Monte Carlo results up to the intermediate coupling regime. We point out a qualitative difference between the isotropic 3d case, and the 2d or anisotropic 3d cases for the double occupation factor. Even for this local correlation function, 2d or anisotropic 3d cases are out of reach of DMF: this comes from the inability of DMF to account for antiferromagnetic fluctuations, which are crucial.Comment: RevTex, 9 pages +10 figure

Effects of Boson Dispersion in Fermion-Boson Coupled Systems

We study the nonlinear feedback in a fermion-boson system using an extension of dynamical mean-field theory and the quantum Monte Carlo method. In the perturbative regimes (weak-coupling and atomic limits) the effective interaction among fermions increases as the width of the boson dispersion increases. In the strong coupling regime away from the anti-adiabatic limit, the effective interaction decreases as we increase the width of the boson dispersion. This behavior is closely related with complete softening of the boson field. We elucidate the parameters that control this nonperturbative region where fluctuations of the dispersive bosons enhance the delocalization of fermions.Comment: 14 pages RevTeX including 12 PS figure

Magnetic and Dynamic Properties of the Hubbard Model in Infinite Dimensions

An essentially exact solution of the infinite dimensional Hubbard model is made possible by using a self-consistent mapping of the Hubbard model in this limit to an effective single impurity Anderson model. Solving the latter with quantum Monte Carlo procedures enables us to obtain exact results for the one and two-particle properties of the infinite dimensional Hubbard model. In particular we find antiferromagnetism and a pseudogap in the single-particle density of states for sufficiently large values of the intrasite Coulomb interaction at half filling. Both the antiferromagnetic phase and the insulating phase above the N\'eel temperature are found to be quickly suppressed on doping. The latter is replaced by a heavy electron metal with a quasiparticle mass strongly dependent on doping as soon as $n<1$. At half filling the antiferromagnetic phase boundary agrees surprisingly well in shape and order of magnitude with results for the three dimensional Hubbard model.Comment: 32 page

Asymptotics of self-similar solutions to coagulation equations with product kernel

We consider mass-conserving self-similar solutions for Smoluchowski's coagulation equation with kernel $K(\xi,\eta)= (\xi \eta)^{\lambda}$ with $\lambda \in (0,1/2)$. It is known that such self-similar solutions $g(x)$ satisfy that $x^{-1+2\lambda} g(x)$ is bounded above and below as $x \to 0$. In this paper we describe in detail via formal asymptotics the qualitative behavior of a suitably rescaled function $h(x)=h_{\lambda} x^{-1+2\lambda} g(x)$ in the limit $\lambda \to 0$. It turns out that $h \sim 1+ C x^{\lambda/2} \cos(\sqrt{\lambda} \log x)$ as $x \to 0$. As $x$ becomes larger $h$ develops peaks of height $1/\lambda$ that are separated by large regions where $h$ is small. Finally, $h$ converges to zero exponentially fast as $x \to \infty$. Our analysis is based on different approximations of a nonlocal operator, that reduces the original equation in certain regimes to a system of ODE