Even though the Hubbard model is one of the most fundamental models of highly correlated electrons, analytical and numerical data describing its thermodynamics at nonzero magnetization are relatively scarce. We present a detailed investigation of the thermodynamic properties for the one-dimensional repulsive Hubbard model in the presence of an arbitrary magnetic field for all values of the filling fraction and temperatures as low as $T\ensuremath{\sim}0.005t.$ Our analysis is based on the system of integral equations derived in the quantum transfer matrix framework. We determine the critical exponents of the quantum phase transitions and also provide analytical derivations for some of the universal functions characterizing the thermodynamics in the vicinities of the quantum critical points. Extensive numerical data for the specific heat, susceptibility, compressibility, and entropy are reported. The experimentally relevant double occupancy presents an interesting doubly nonmonotonic temperature dependence at intermediate values of the interaction strength and also at large repulsion and magnetic fields close to the critical value. The susceptibility in zero magnetic field has a logarithmic singularity at low temperatures for all filling factors similar to the behavior of the same quantity in the spin-1/2 isotropic Heisenberg model. We determine the density profiles for a harmonically trapped system and show that while the total density profile seems to depend mainly on the value of chemical potential at the center of the trap the distribution of phases in the inhomogeneous system changes dramatically as we increase the magnetic field.
Even though the Hubbard model is one of the most fundamental models of highly correlated electrons, analytical and numerical data describing its thermodynamics at nonzero magnetization are relatively scarce. We present a detailed investigation of the thermodynamic properties for the one-dimensional repulsive Hubbard model in the presence of an arbitrary magnetic field for all values of the filling fraction and temperatures as low as T\ensuremath∼0.005t. Our analysis is based on the system of integral equations derived in the quantum transfer matrix framework. We determine the critical exponents of the quantum phase transitions and also provide analytical derivations for some of the universal functions characterizing the thermodynamics in the vicinities of the quantum critical points. Extensive numerical data for the specific heat, susceptibility, compressibility, and entropy are reported. The experimentally relevant double occupancy presents an interesting doubly nonmonotonic temperature dependence at intermediate values of the interaction strength and also at large repulsion and magnetic fields close to the critical value. The susceptibility in zero magnetic field has a logarithmic singularity at low temperatures for all filling factors similar to the behavior of the same quantity in the spin-1/2 isotropic Heisenberg model. We determine the density profiles for a harmonically trapped system and show that while the total density profile seems to depend mainly on the value of chemical potential at the center of the trap the distribution of phases in the inhomogeneous system changes dramatically as we increase the magnetic field.
We present analytical results of fundamental properties of one-dimensional (1D) Hubbard model with a repulsive interaction, ranging from fractional excitations to universal thermodynamics, interaction-driven criticality, correlation functions, Contact susceptibilities and quantum cooling. Using the exact solutions of the Bethe Ansatz equations of the Hubbard model, we first rigorously calculate the gapless spin and charge excitations, exhibiting exotic features of fractionalized spinons and holons. Based on the analysis on the fractional charge and spin excitations, the spin-incoherent Luttinger liquid with only the charge propagation mode is elucidated by the asymptotic of the two-point correlation functions with the help of the conformal field theory. Near quadruple critical point, we then further analytically obtain the thermodynamical properties, dimensionless ratios and scaling functions near quantum phase transitions in terms of chemical potential, magnetic field and interaction. In particular, we determine additivity rules of spin and charge susceptibilities, and derive explicit forms of thermodynamics of spin-incoherent Luttinger liquid. Finally, in order to capture deeper insight into the Mott insulator and interaction driven criticality, we further study the double occupancy and its associated Contact and Contact susceptibilities through which an adiabatic cooling scheme upon the quantum criticality is introduced.
Abstract We present analytical results of the fundamental properties of the one-dimensional (1D) Hubbard model with a repulsive interaction. The new model results with arbitrary external fields include: (I) using the exact solutions of the Bethe ansatz equations of the Hubbard model, we first rigorously calculate the gapless spin and charge excitations, exhibiting exotic features of fractionalized spinons and holons. We then investigate the gapped excitations in terms of the spin string and the k − Λ string bound states at arbitrary driving fields, showing subtle differences in spin magnons and charge η -pair excitations. (II) For a high-density and high spin magnetization region, i.e. near the quadruple critical point, we further analytically obtain the thermodynamic properties, dimensionless ratios and scaling functions near quantum phase transitions. (III) Importantly, we give the general scaling functions at quantum criticality for arbitrary filling and interaction strength. These can directly apply to other integrable models. (IV) Based on the fractional excitations and the scaling laws, the spin-incoherent Luttinger liquid (SILL) with only the charge propagation mode is elucidated by the asymptotic of the two-point correlation functions with the help of conformal field theory. We also, for the first time, obtain the analytical results of the thermodynamics for the SILL. (V) Finally, to capture deeper insights into the Mott insulator and interaction-driven criticality, we further study the double occupancy and propose its associated contact and contact susceptibilities, through which an adiabatic cooling scheme based upon quantum criticality is proposed. In this scenario, we build up general relations among arbitrary external- and internal-potential-driven quantum phase transitions, providing a comprehensive understanding of quantum criticality. Our methods offer rich perspectives of quantum integrability and offer promising guidance for future experiments with interacting electrons and ultracold atoms, both with and without a lattice.
The exact phase diagrams of the one-dimensional Hubbard model, both attractive and repulsive, in the presence of an arbitrary magnetic field h for various electron concentrations n and on-site interaction strengths U <0 or U >0 are calculated and investigated. The exact ground-state properties, namely, the ground-state energy, the average spin (magnetization), the concentration of the doubly occupied sites, the kinetic energy, the chemical potential, the spin (magnetic) susceptibility and the charge compressibility, are calculated and examined over a wide range of interaction strengths U for various h and n . It is found that the spin susceptibility at half-filling is non-analytic and changes discontinuously as U →0. The exact theory shows the absence of a charge energy gap in the U ⩽0 region for all n and provides, for the chemical potential, the rigorous upper and lower bounds for half-filled and empty bands respectively. The analytical results derived for the weak-coupling limit and asymptotic expansions for strong coupling are in full agreement with the numerical calculations.
The thermodynamics of the inhomogeneous one-dimensional repulsive fermionic Hubbard model with parabolic confinement is studied by a density-functional-theory approach, based on Mermin's generalization to finite temperatures. A local-density approximation (LDA), based on exact results for the homogeneous model, is used to approximate the correlation part in the Helmholtz free-energy, comprising the thermodynamic Bethe ansatz LDA (TBALDA). The general presentation of the method is given and some properties of the homogeneous model that are relevant to the DFT approach are analyzed. Extensive comparison between TBALDA and numerical exact diagonalization results for thermodynamic properties of small inhomogeneous chains is discussed. In the remaining, a classical thermodynamic treatment of the confined system is developed with the focus on global properties of large systems. It is observed that, depending on the temperature and specific volume, the system can have its temperature increased under isentropic expansion. Such unusual behavior becomes more pronounced with the increase of the intrasite interaction U and can be understood considering the peculiar behavior of the correlation entropy per site of the Hubbard model around half-filling. A related nonmonotonic variation of the temperature under isentropic change of U is also studied. At the end, the dependence of the double occupancy on temperature and volume per particle is discussed.
The phase diagram, the Bethe-ansatz ground-state properties, including the chemical potential μ, the spin (magnetic) and charge susceptibilities, are calculated within the one-dimensional Hubbard model in entire range of interaction strength (-∞< U/t<+∞), magnetic field (h≥ 0) and all electron concentrations (0≤n≤1). The continuous and smooth variation of μ with n and h in the vicinity of n=1 points on the gapless character of charge excitations at U<0 and provides rigorous upper and lower bounds for μ. The spin (magnetic) susceptability χ at half-filling changes discontinuously as U→0 and is strongly enhanced by electron repulsion, comparing with that of the non-interactig case. The compressibility κ ch increases with n at U<0 and shows non-monotonous behavior with a dramatic increase at U>0. Variations of κ ch -1 in both repulsive and attractive cases qualitatively well reproduces corresponding behavior of charge stiffness.
The magnetic susceptibility of the one-dimensional Hubbard model with open boundary conditions at arbitrary filling is obtained from field theory at low temperatures and small magnetic fields, including leading and next-leading orders. Logarithmic contributions to the bulk part are identified as well as algebraic-logarithmic divergences in the boundary contribution. As a manifestation of spin-charge separation, the result for the boundary part at low energies turns out to be independent of filling and interaction strength and identical to the result for the Heisenberg model. For the bulk part at zero temperature, the scale in the logarithms is determined exactly from the Bethe ansatz. At finite temperature, the susceptibility profile as well as the Friedel oscillations in the magnetisation are obtained numerically from the density-matrix renormalisation group applied to transfer matrices. Agreement is found with an exact asymptotic expansion of the relevant correlation function.