Low-temperature static thermodynamic properties of the linear half-filled-band Hubbard model are computed via a numerical-perturbative method which applies near the atomic limit. Curves for the internal energy, entropy, specific heat, and magnetic susceptibility are reported.
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.
A new first-principles method for numerically calculating finite-temperature properties of quantum many-body systems is described. It is used to compute the thermodynamic properties and magnetic correlations of the Hubbard model, with negligible error at temperatures near the bandwidth or above. Results for the half-filled band and a simple cubic or one-dimensional lattice confirm the existence of a high-temperature peak in the specific heat which has been associated with a smooth but rapid change in conductivity. Correlations calculated for the non-half-filled band simple cubic Hubbard model suggest condensation, an entirely unexpected phenomenon whose existence is shown (by unrelated methods) in Paper II.
The dimerized-chain Hubbard model is characterized by a Coulomb energy U and two alternating transfer constants, an intradimer constant t, and a smaller interdimer constant ${t}^{\ensuremath{'}}$. We have used thermodynamic perturbation theory to treat the half-filled-band case, taking the interdimer transfer as the perturbation term. The unperturbed problem is a lattice of isolated dimers. We give the energy, specific heat, and magnetic susceptibility as functions of temperature for the unperturbed case and for the first nonvanishing order in the perturbation. The results are consistent with the U=0 case, which is solved exactly. The interdimer transfer ${t}^{\ensuremath{'}}$ has the effect of distorting the peaks in the specific heat and reducing the susceptibility; this is shown graphically.
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.
In order to determine finite-temperature properties of the infinite one-dimensional Hubbard model for the half-filled band, finite systems such as chains with 2 to 5 atoms and rings with 4 and 6 atoms are studied with the application of the grand canonical ensemble. The obtained temperature dependence of the specific heat, the entropy, the magnetic susceptibility and some correlation functions clearly shows a rearrangement of the electronic state of the system with the increase of the strength of the Coulomb interaction relative to the transfer integral. In the light of our new results the functional integral approach recently proposed by Kimball and Schrieffer as well as Hubbard's approximate theory is examined.