Acta Physica Sinica

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

Tianbao Zhang, Yu Xuan-Ping, Ahai Chen

Published Jan 1, 2015• 0 citations

Abstract

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.

DOI: 10.7498/aps.64.156402

Related papers

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

Zhanpeng Lu, Xingbo Wei, Liu Tian-Shuai +2 more

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.

Acta Physica SinicaJan 1, 20170 citations

Universal low-energy physics in 1D strongly repulsive multi-component Fermi gases

Yuzhu Jiang, Peng He, Xi‐Wen Guan

Abstract It has been shown (Yang and You 2011 Chin. Phys. Lett. 28 020503 ) that at zero temperature the ground state of the one-dimensional (1D) w -component Fermi gas coincides with that of the spinless Bose gas in the limit . This behavior was experimentally evidenced through quasi-1D tightly trapping ultracold 173 Yb atoms in a recent paper (Pagano et al 2014 Nat. Phys. 10 198 ). However, understanding of low-temperature behavior of Fermi gases with a repulsive interaction requires spin-charge separated conformal field theories of an effective Tomonaga–Luttinger liquid and an antiferromagnetic SU ( w ) Heisenberg spin chain. Here we analytically derive universal thermodynamics of 1D strongly repulsive fermionic gases with SU ( w ) symmetry via the Yang–Yang thermodynamic Bethe ansatz method. The analytical free energy and magnetic properties of the systems at low temperature in a weak magnetic field are obtained through the Wiener–Hopf method. In particular, the free energy essentially manifests the spin-charge separated conformal field theories for high-symmetry systems with arbitrary repulsive interaction strength. We also find that the sound velocity of the Fermi gases in the large w limit coincides with that for the spinless Bose gas, whereas the spin velocity vanishes quickly as w becomes large. This indicates strong suppression of the Fermi exclusion statistics by the commutativity feature among the w -component fermions with different spin states in the Tomonaga–Luttinger liquid phase. Moreover, the equations of state and critical behavior of physical quantities at finite temperature are analytically derived in terms of the polylogarithm functions in the quantum critical region.

arXiv:1602.02507Journal of Physics A Mathematical and TheoreticalMar 18, 20167 citations

Influence of thermal fluctuations on quantum phase transitions in one-dimensional disordered systems: Charge density waves and Luttinger liquids

Andreas Glatz, Thomas Nattermann

The low-temperature phase diagram of one-dimensional weakly disordered quantum systems such as charge or spin density waves and Luttinger liquids is studied by a full finite temperature renormalization-group (RG) calculation. For vanishing quantum fluctuations this approach is amended by an exact solution in the case of strong disorder and by a mapping onto the Burgers equation with noise in the case of weak disorder, respectively. At zero temperature we reproduce the quantum phase transition between a pinned (localized) and an unpinned (delocalized) phase for weak and strong quantum fluctuations, respectively, as found previously by Fukuyama or Giamarchi and Schulz. At finite temperatures the localization transition is suppressed: the random potential is wiped out by thermal fluctuations on length scales larger than the thermal de Broglie wavelength of the phason excitations. The existence of a zero-temperature transition is reflected in a rich crossover phase diagram of the correlation functions. In particular we find four different scaling regions: a classical disordered, a quantum disordered, a quantum critical, and a thermal region. The results can be transferred directly to the discussion of the influence of disorder in superfluids. Finally we extend the RG calculation to the treatment of a commensurate lattice potential. Applications to related systems are discussed as well.

arXiv:cond-mat/0310432Physical Review BMar 24, 20049 citations

Equilibrium thermodynamic properties of interacting two-component bosons in one dimension

Antoine Klauser, Jean-Sébastien Caux

The interplay of quantum statistics, interactions, and temperature is studied within the framework of the bosonic two-component theory with repulsive delta-function interaction in one dimension. We numerically solve the thermodynamic Bethe ansatz and obtain the equation of state as a function of temperature and of the interaction strength, the relative chemical potential, and either the total chemical potential or a fixed number of particles, allowing quantification of the full crossover behavior of the system between its low-temperature ferromagnetic and high-temperature unpolarized regime, and from the low coupling decoherent regime to the fermionization regime at high interaction.

arXiv:1106.4971Physical Review ASep 2, 201120 citations

Fermi gases in one dimension: From Bethe ansatz to experiments

Xi‐Wen Guan, Murray T. Batchelor, Chaohong Lee

This article reviews theoretical and experimental developments for one-dimensional Fermi gases. Specifically, the experimentally realized two-component delta-function interacting Fermi gas---the Gaudin-Yang model---and its generalizations to multicomponent Fermi systems with larger spin symmetries is discussed. The exact results obtained for Bethe ansatz integrable models of this kind enable the study of the nature and microscopic origin of a wide range of quantum many-body phenomena driven by spin population imbalance, dynamical interactions, and magnetic fields. This physics includes Bardeen-Cooper-Schrieffer-like pairing, Tomonaga-Luttinger liquids, spin-charge separation, Fulde-Ferrel-Larkin-Ovchinnikov-like pair correlations, quantum criticality and scaling, polarons, and the few-body physics of the trimer state (trions). The fascinating interplay between exactly solved models and experimental developments in one dimension promises to yield further insight into the exciting and fundamental physics of interacting Fermi systems.

arXiv:1301.6446Reviews of Modern PhysicsNov 27, 2013544 citations

Thermodynamics, contact, and density profiles of the repulsive Gaudin-Yang model

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

We address the problem of computing the thermodynamic properties of the repulsive one-dimensional two-component Fermi gas with contact interaction, also known as the Gaudin-Yang model. Using a specific lattice embedding and the quantum transfer matrix we derive an exact system of only two nonlinear integral equations for the thermodynamics of the homogeneous model which is valid for all temperatures and values of the chemical potential, magnetic field, and coupling strength. This system allows for an easy and extremely accurate calculation of thermodynamic properties circumventing the difficulties associated with the truncation of the thermodynamic Bethe ansatz system of equations. We present extensive results for the densities, polarization, magnetic susceptibility, specific heat, interaction energy, Tan contact, and local correlation function of opposite spins. Our results show that at low and intermediate temperatures the experimentally accessible contact is a nonmonotonic function of the coupling strength. As a function of the temperature the contact presents a pronounced local minimum in the Tonks-Girardeau regime which signals an abrupt change of the momentum distribution in a small interval of temperature. The density profiles of the system in the presence of a harmonic trapping potential are computed using the exact solution of the homogeneous model coupled with the local density approximation. We find that at finite temperature the density profile presents a double shell structure (partially polarized center and fully polarized wings) only when the polarization in the center of the trap is above a critical value which is monotonically increasing with temperature.

arXiv:1512.07485Physical Review AMar 9, 201624 citations