Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Acta Physica Sinica

Chemical potential-functional-theory about the properties of one-dimensional Hubbard model at finite temperature

Zhanpeng Lu, Xingbo Wei, Liu Tian-Shuai, Ahai Chen, Gao Xianlong

Published Jan 1, 2017• 0 citations
Publisher

Abstract

In this paper, we numerically solve the thermodynamic Bethe-ansatz coupled equations for a one-dimensional Hubbard model at finite temperature and obtain the second order thermodynamics properties, such as the specific heat, compressibility, and susceptibility. We find that these three quantities could embody the phase transitions of the system, from the vacuum state to the metallic state, from the metallic state to the Mott-insulating phase, from the Mott-insulating phase to the metallic state, and from the metallic state to the band-insulating phase. With the increase of temperature, the thermal fluctuation overwhelms the quantum fluctuations and the phase transition points disappear due to the destruction of the Mott-insulating phase. But in the case of the strong interaction strength, the Mott-insulating phase is robust, embodying the compressibility. Furthermore, we study the thermodynamic properties of the inhomogeneous Hubbard model with trapping potential. Making use of the Bethe-ansatz results from the homogeneous Hubbard model, we construct the chemical potential-functional theory (-BALDA) for the inhomogeneous Hubbard model instead of the commonly used density-functional theory, in order to solve the in-convergence problem of the Kohn-Sham equation in the case of the divergence appearing in the exchange-correlation potential. We further point out a multi-dimensional bisection method which changes the Kohn-Shan equation into a problem of finding the fixed points. Through -BALDA we numerically solve the one-dimensional homogeneous Hubbard model of trapping potential. The density profile and the local compressibility are obtained. We find that at a given interaction strength, the metallic phase and the Mott-insulating phase are destroyed and the density profile becomes a Guassian distribution with increasing temperature. To the metallic phase, Friedel oscillation caused by quantum fluctuations is still visible at low temperature. With increasing temperature, Friedel oscillation will disappear. This situation reflects the fact that the thermal fluctuation overwhelms the quantum fluctuations. For the Mott-insulating phase, the Mott-insulating plateau is robust at a certain temperature and only the boundary of the Mott-insulating plateau is destroyed. With increasing temperature, the Mott insulating plateau will be destroyed. And the change of the local compressibility provides the information about such a change. So we conclude that the thermal fluctuation destroys the original quantum phase. Through our analysis, we find that the -BALDA can be used to study the finite temperature properties for the system of the exchange-correlation potential divergence with high efficiency.

DOI
10.7498/aps.66.126701

Related papers

Thermodynamic Relations in Correlated Systems

Shinji Watanabe, Masatoshi Imada

Several useful thermodynamic relations are derived for metal-insulator transitions, as generalizations of the Clausius-Clapeyron and Eherenfest theorems. These relations hold in any spatial dimensions and at any temperatures. First, they relate several thermodynamic quantities to the slope of the metal-insulator phase boundary drawn in the plane of the chemical potential and the Coulomb interaction in the phase diagram of the Hubbard model. The relations impose constraints on the critical properties of the Mott transition. These thermodynamic relations are indeed confirmed to be satisfied in the cases of the one- and two-dimensional Hubbard models. One of these relations yields that at the continuous Mott transition with a diverging charge compressibility, the doublon susceptibility also diverges. The constraints on the shapes of the phase boundary containing a first-order metal-insulator transition at finite temperatures are clarified based on the thermodynamic relations. For example, the first-order phase boundary is parallel to the temperature axis asymptotically in the zero temperature limit. The applicability of the thermodynamic relations are not restricted only to the metal-insulator transition of the Hubbard model, but also hold in correlated systems with any types of phases in general. We demonstrate such examples in an extended Hubbard model with intersite Coulomb repulsion containing the charge order phase.

arXiv:cond-mat/0412602•Journal of the Physical Society of Japan•Dec 15, 2004•14 citations

Thermodynamic properties of one-dimensional Gaudin-Yang model at finite temperature

Tianbao Zhang, Yu Xuan-Ping, Ahai Chen

The one-dimensional system interacting via a delta-function interparticle interaction is a very important one in cold atomic systems and has fundamental importance in many-body physics. In one dimension, due to the geometric confinement induced quantum correlations and quantum fluctuations, there may exist a number of unusual phenomena, such as spin-charge separation, effective fermionization and quantum criticality. This paper studies the basic properties of a uniform one-dimensional Gaudin-Yang model for fermions by solving the thermodynamic Bethe-ansatz equations by a numerical method. Numerically, we use the many-variable Newton’s method to solve the coupled equations. We analyze the physical properties, including density, interaction, temperature and entropy at a given temperature and a given interaction, separately. We know that a lot of researches are limited to zero temperature. However, we cannot reach the absolute zero temperature in the real cold atomic experiment. So it is important to deal with the finite temperature problems. We study the density and entropy as a function of the chemical potential, temperature and interaction and, then give the phase diagrams, respectively. We found that there is a quantum critical zone in the phase diagram of entropy, including the high temperature zone with thermal fluctuations and the Luttinger liquid zone with quantum fluctuations. For a given temperature and low chemical potential, the thermal fluctuations are the main factor in the entropy. With the increase of chemical potential, the system enters the quantum critical zone where the competitive effect between the thermal fluctuations and the quantum fluctuations exists. When the chemical potential is large enough, the quantum fluctuations become the main factor in the system’s entropy, and we get the Luttinger liquid phase. Our results can be further used in the finite temperature density-functional theory and to analyze the collective phenomena at a finite temperature.

Acta Physica Sinica•Jan 1, 2015•0 citations

Quantum critical behavior and thermodynamics of the repulsive one-dimensional Hubbard model in a magnetic field

Ovidiu I. Pâţu, Andreas Klümper, Angela Foerster

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.

arXiv:1910.09535•Physical review. B./Physical review. B•Jan 29, 2020•10 citations

Thermodynamic Properties of One-Dimensional Hubbard Model at Finite Temperatures

Chi Yang, ARMEN N. KOCHARIAN, You-Ling Chiang

The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.

International Journal of Modern Physics B•Aug 10, 2003•0 citations

Thermodynamic properties of one-dimensional Hubbard model at finite temperatures

Chia‐Wei Yang

The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.

•Aug 10, 2003•1 citation

Variational cluster approach to thermodynamic properties of interacting fermions at finite temperatures: A case study of the two-dimensional single-band Hubbard model at half filling

Kazuhiro Seki, Tomonori Shirakawa, Seiji Yunoki

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.

arXiv:1808.03807•Physical review. B./Physical review. B•Nov 7, 2018•26 citations
In this paper, we numerically solve the thermodynamic Bethe-ansatz coupled equations for a one-dimensional Hubbard model at finite temperature and obtain the second order thermodynamics properties, such as the specific heat, compressibility, and susceptibility. We find that these three quantities could embody the phase transitions of the system, from the vacuum state to the metallic state, from the metallic state to the Mott-insulating phase, from the Mott-insulating phase to the metallic state, and from the metallic state to the band-insulating phase. With the increase of temperature, the thermal fluctuation overwhelms the quantum fluctuations and the phase transition points disappear due to the destruction of the Mott-insulating phase. But in the case of the strong interaction strength, the Mott-insulating phase is robust, embodying the compressibility. Furthermore, we study the thermodynamic properties of the inhomogeneous Hubbard model with trapping potential. Making use of the Bethe-ansatz results from the homogeneous Hubbard model, we construct the chemical potential-functional theory (-BALDA) for the inhomogeneous Hubbard model instead of the commonly used density-functional theory, in order to solve the in-convergence problem of the Kohn-Sham equation in the case of the divergence appearing in the exchange-correlation potential. We further point out a multi-dimensional bisection method which changes the Kohn-Shan equation into a problem of finding the fixed points. Through -BALDA we numerically solve the one-dimensional homogeneous Hubbard model of trapping potential. The density profile and the local compressibility are obtained. We find that at a given interaction strength, the metallic phase and the Mott-insulating phase are destroyed and the density profile becomes a Guassian distribution with increasing temperature. To the metallic phase, Friedel oscillation caused by quantum fluctuations is still visible at low temperature. With increasing temperature, Friedel oscillation will disappear. This situation reflects the fact that the thermal fluctuation overwhelms the quantum fluctuations. For the Mott-insulating phase, the Mott-insulating plateau is robust at a certain temperature and only the boundary of the Mott-insulating plateau is destroyed. With increasing temperature, the Mott insulating plateau will be destroyed. And the change of the local compressibility provides the information about such a change. So we conclude that the thermal fluctuation destroys the original quantum phase. Through our analysis, we find that the -BALDA can be used to study the finite temperature properties for the system of the exchange-correlation potential divergence with high efficiency.
Several useful thermodynamic relations are derived for metal-insulator transitions, as generalizations of the Clausius-Clapeyron and Eherenfest theorems. These relations hold in any spatial dimensions and at any temperatures. First, they relate several thermodynamic quantities to the slope of the metal-insulator phase boundary drawn in the plane of the chemical potential and the Coulomb interaction in the phase diagram of the Hubbard model. The relations impose constraints on the critical properties of the Mott transition. These thermodynamic relations are indeed confirmed to be satisfied in the cases of the one- and two-dimensional Hubbard models. One of these relations yields that at the continuous Mott transition with a diverging charge compressibility, the doublon susceptibility also diverges. The constraints on the shapes of the phase boundary containing a first-order metal-insulator transition at finite temperatures are clarified based on the thermodynamic relations. For example, the first-order phase boundary is parallel to the temperature axis asymptotically in the zero temperature limit. The applicability of the thermodynamic relations are not restricted only to the metal-insulator transition of the Hubbard model, but also hold in correlated systems with any types of phases in general. We demonstrate such examples in an extended Hubbard model with intersite Coulomb repulsion containing the charge order phase.
The one-dimensional system interacting via a delta-function interparticle interaction is a very important one in cold atomic systems and has fundamental importance in many-body physics. In one dimension, due to the geometric confinement induced quantum correlations and quantum fluctuations, there may exist a number of unusual phenomena, such as spin-charge separation, effective fermionization and quantum criticality. This paper studies the basic properties of a uniform one-dimensional Gaudin-Yang model for fermions by solving the thermodynamic Bethe-ansatz equations by a numerical method. Numerically, we use the many-variable Newton’s method to solve the coupled equations. We analyze the physical properties, including density, interaction, temperature and entropy at a given temperature and a given interaction, separately. We know that a lot of researches are limited to zero temperature. However, we cannot reach the absolute zero temperature in the real cold atomic experiment. So it is important to deal with the finite temperature problems. We study the density and entropy as a function of the chemical potential, temperature and interaction and, then give the phase diagrams, respectively. We found that there is a quantum critical zone in the phase diagram of entropy, including the high temperature zone with thermal fluctuations and the Luttinger liquid zone with quantum fluctuations. For a given temperature and low chemical potential, the thermal fluctuations are the main factor in the entropy. With the increase of chemical potential, the system enters the quantum critical zone where the competitive effect between the thermal fluctuations and the quantum fluctuations exists. When the chemical potential is large enough, the quantum fluctuations become the main factor in the system’s entropy, and we get the Luttinger liquid phase. Our results can be further used in the finite temperature density-functional theory and to analyze the collective phenomena at a finite temperature.
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.
The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.
The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.
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 UUU and temperature TTT. We also calculate the TTT 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 TTT and/or UUU derivative of the grand potential, we obtain a crossover diagram in the (U,T)(U,T)(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)(U,T)(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.