The thermodynamic and correlation functions for one- and three- dimensional systems described by the half-filled band Hubbard Hamiltonians are evaluated in the pair-approximation of cluster variation method. The energy, specific heat, magnetization, susceptibility and various equilibrium correlation functions are computed as functions of temperature as well as the external magnetic field for different values of coulomb energies. Our solutions correspond to Bethe approximation in which the relevant pair-correlation functions are included in our formulation. It is observed that the chemical potential for the half-filled band is exactly half of the coulomb energy at all temperatures, irrespective of the value of the applied magnetic field.
We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion
U and temperature
T. We also calculate the
T dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the
T and/or
U derivative of the grand potential, we obtain a crossover diagram in the
(U,T) plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the
(U,T) plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.
The one-dimensional Hubbard model at half-filling is studied in the framework of the Composite Operator Method using a static approximation. A solution characterized by strong antiferromagnetic correlations and a gap for any non-zero on-site interaction U is found. The corresponding ground-state energy, double occupancy and specific heat are in excellent agreement with those obtained within the Bethe Ansatz. These results show that the Composite Operator Method is an appropriate framework for the half-filled Hubbard chain and can be applied to evaluate properties, like the correlation functions, which cannot be obtained by means of the Bethe Ansatz, except for some limiting cases.
A detailed quantitative study of the physical properties of the infinite-dimensional Hubbard model at half filling is presented. The method makes use of an exact mapping onto a single-impurity model supplemented by a self-consistency condition. This coupled problem is solved numerically. Results for thermodynamic quantities (specific heat, entropy, . . .), one-particle spectral properties, and magnetic properties (response to a uniform magnetic field) are presented and discussed. The nature of the Mott-Hubbard metal-insulator transition found in this model is investigated. A numerical solution of the mean-field equations inside the antiferromagnetic phase is also reported.
In this article we review the properties of the 2D Hubbard model by considering at the same time the cases of repulsive and attractive interaction. The paramagnetic solution is studied by means of the composite operator method in the static approximation for the case of half-filling. Some properties of the two models, as the double occupancy and the spin magnetic susceptibility, are calculated for various values of interaction and temperature and compared. In particular, the different role played by thermal fluctuations is analyzed. Analytical and numerical calculations show that there is a critical value of the interaction, Uc, where the system exhibits a metal-insulator transition. At zero temperature it is found that Uc=−W for the negative-U model and Uc≈1.68W for the positive-U model, where W is the band width. At zero temperature, when the strength of the attractive interaction equals the band width, the system exhibits a phase transition to a pair state, where all the electrons are locally paired. The temperature Tp which controls the crossover to the pair state is calculated as a function of U. For strong attractive interaction χ0 is strongly depressed; increases by increasing T and tends to zero as T→TP.
In order to investigate the thermal properties (specific heat, magnetic susceptibility, entropy, internal energy, and some correlation functions) of the one-dimensional half-filled-band Hubbard model, we have studied linear chains and rings containing two to six atoms, by performing machine calculations. Supplementing the low-temperature behavior obtained from the exact solution for the infinite chain by Lieb and Wu, our results should be suggestive of the properties of the infinite system throughout the entire temperature domain. It is shown that when the ratio of the correlation energy
U to the total width $\ensuremath{\Delta}$ of the band of single-particle excitations is larger than 1, the specific heat has two peaks. The high-temperature peak arises from the gradual metal-insulator transition (or the gradual formation of local moments), while the low-temperature peak is associated with the antiferromagnetic short-range ordering. When $\frac{U}{\ensuremath{\Delta}}$ becomes small, the two peaks merge into one. This picture is consistent with all other thermal properties, including correlation functions. The high-temperature properties are compared with the results predicted by Hubbard's approximate theory based on the truncation of the equations of motion of the Green's functions.