arXiv:1904.02283v1 [physics.hist-ph] 3 Apr 2019 nlivnin e oimdaeapiain nmany in physics. applications mathematical immediate of areas to research orig- led Their inventions 1D. in inal problems mathe- Bethe physical of transcendental for the initiation equations of the ansatz study is the impor- in papers key analysis of matical The des series their [11] others. of and Walker tance [12] [10], Pearson by and Griffith earlier Cloizeaux [9], obtained Orbach results the the of Hulth´en [8], to magnetization adding of thereby behaviour model, critical dia- the phase and pressure-volume gram the magnetization, the ergy, h ulrneo nstoi aaee n magnetic and ∆ parameter throughout anisotropic chain field of ansatz spin range Bethe Heisenberg full the 1D the of analysis the rigorous presented for published a 2] equations [1, time papers Yang first of and the Yang series which for in a Heisen- mid-60’s, of the the delta- consists in on 1D chain work of spin The berg thermodynamics gas. Bose the Heisen- interacting and (1D) function chain one-dimensional spin the on berg Yang, professor Ning with collaboration Chen in did Yang Ping Chen fessor time. that at ar- physics research mathematical various of opened problems eas [7] physical hypothesis on are Bethe’s works 60’s using His the influential. in and made con- seminal physics early mathematical his to [1–6], did tributions papers he scientific Although many 1998. published in not State Uni- Ohio retirement The Hopkins his at Johns physics until taught the University then at He PhD 1960. he in his Uni- versity Later, Harvard and un- at 1953, his Science 1948. in of versity pursue of Degree to Master summer his University the completed very Brown in a to studies commu- was went dergraduate physics he He mathematical mind, our our nity. To in physicist May. humble this in away passed rfso hnPn agserysgicn contributions significant early Yang’s Ping Chen Professor h oko h hroyaiso h oegswas gas Bose the of thermodynamics the on work The esalfcshr rmrl ntowrs hc Pro- which works, two on primarily here focus shall We Yang Ping Chen Professor that news sad very a was It H oevr hyotie h rudsaeen- state ground the obtained they Moreover, . ua nttt fPyisadMteais hns Academ Chinese Mathematics, and Physics of Institute Wuhan euiu xeietlvrfiain nrcn years. recent Yan in The verifications le mechanics. experimental interaction statistical beautiful delta-function quantum develo proposed in a subsequent Yang problems with of and to bosons range Yang of led which system chain 1D method stat the spin The and Heisenberg physics hypothesis. 1D mathematical Bethe’s of the fields on interac the work th delta-function in 1D transition, influential of liquid-gas and thermodynamics in the proble and behaviour various chain, for critical hypothesis model, Bethe’s of gas study the on papers nte6’ rfso hnPn agwt rfso hnNing Chen Professor with Yang Ping Chen Professor 60’s the In .INTRODUCTION I. 1 2 tt e aoaoyo antcRsnneadAoi n Mo and Atomic and Resonance Magnetic of Laboratory Key State eateto hoeia hsc,Rsac colo Phys of School Research Physics, Theoretical of Department utainNtoa nvriy abraAT00,Austra 0200, ACT Canberra University, National Australian iWnGuan Xi-Wen ,2, 1, ∗ al eh nazitgal oesi ait ffields of variety no- a of in number models integrable A ansatz Bethe models. table integrable quantum of history Yang-Yang equation. name ansatz the Bethe bears thermodynamic equation permanently The obtained mechanics. influential they statistical quasi- and quantum of profound on has impact terms method and in This potential, particles chemical temperature. the strength, interaction of momenta, energy determines that dressed equation Yang-Yangthe min- called ansatz a so Bethe Such the thermodynamic to densities. rise hole gives en- and condition free particle imization of Gibbs terms the in of be ergy minimization can integrable the thermodynamics from of the signifi- determined that properties showed a macroscopic They to treat systems. led thermodynamics to and exact step systems cant first interacting the many-body was delta-function of with This gas the Bose interaction. calculate 1D Yang to the and of ensemble Yang thermodynamics canonical where grand 1969, the in proposed published [24] paper the Chen Deng, Xian Jia 1949. Yang, in Ning Chicago of Chen University left: at Yang From Ping 1: FIG. n egHe Feng and h 0swr rubyb h otectn iei the in time exciting most the be arguably were 60’s The n-iesoa 1)Hiebr spin Heisenberg (1D) one-dimensional e d osgicn plctosi wide a in applications significant to ads so hsc.Tewrso h lattice the on works The physics. of ms igbsn r infiatyimportant significantly are bosons ting n agtemdnmc a found has thermodynamics Yang and g sia ehnc.I atclr the particular, In mechanics. istical fSine,Whn407,China 430071, Wuhan Sciences, of y mnsi aypolm using problems many in pments otettetemdnmc of thermodynamics the treat to 1 agpbihdsvrlseminal several published Yang c n Engineering, and ics omteaia physics mathematical to eua Physics, lecular lia 2 of physics were solved at that time, including the Lieb- case with ∆ = −1. Walker [11] studied the particular Liniger Bose gas [13], the Yang-Gaudin model [14, 15], case of ∆ ≤−1. the Hubbard model [16], the SU(N) interacting Fermi Since Bethe proposed a particular wave function to ob- gases [17], etc. Professor Chen Ning Yang discovered tain the spectrum of the model (1) in 1931 [7], there were the necessary condition for the Bethe ansatz solvability, very few publications on the Bethe’s method for about which is now known as the Yang-Baxter equation, i.e. the 30 years. In Yang and Yang’s series of papers published factorization condition–the scattering matrix of a quan- in the mid-60’s, they [2] carried out an analytical study tum many-body system can be factorized into a product of the Bethe ansatz equations for the Heisenberg spin of many two-body scattering matrices. In the early 70’s chain throughout the full range of anisotropic parameter Professor Rodney Baxter [18] independently showed that ∆ in a presence of magnetic field. Yang and Yang proved such a factorization relation also occurred as the condi- that the ground state energy is an analytical function of tions for commuting transfer matrices in 2D lattice mod- ∆ and magnetization y denoted as f(∆,y) for |∆| < 1, 1 els in statistical mechanics. A short time later, the study where y = L σz, and L is the total number of sites. In of Yang-Baxter integrable models flourished in Canberra, particular, they built up a rigorous analysis of the Bethe St. Petersburg, Stoney Brook, Kyoto schools and other ansatz equationsP of the model so that the function of the places, see reviews [19–23]. The Yang-Baxter equation ground state f(∆,y) was given explicitly. thus became a key theme in many areas of mathematical In this series of papers, they used the inverse tangent physics and statistical mechanics. function to define the phase shift of the exchange of two In particular, besides the above-mentioned two works, spins, namely Professor Chen Ping Yang published other interesting pa- p−q pers [3–6]. Together with Professor Chen Ning Yang [4] −1 ∆ sin( 2 ) he studied the lattice gas model and found a very inter- Θ(p, q) = 2 tan p+q p−q , (2) " cos( ) − ∆ cos( )# esting feature of the specific heat near the phase tran- 2 2 sition, displaying a sharp peak. The aim of this short where p and q are the quasimomenta of two exchanged communication is to elaborate further on the two works spins. Thus quasimomenta pi’s are within the following introduced above. As a matter of fact, we have personally range: benefited a lot from those works of Professor Chen Ping
Yang with Professor Chen Ning Yang. Here we wish to −π
length is given by G/L = E/L−µn−Ts with the relation to the free energy F = G + µN. Here µ is the chemical potential, n is the linear density. It is important to note that the entropy s, the energy E, and the density n are functions of the particle and hole densities subject to the Bethe ansatz equation (5). The minimization condition δG L = 0 with respect to particle density ρ leads to the so called Yang-Yang thermodynamic Bethe ansatz (TBA) equation [24]
∞ 2 Tc dq − ε(q) ε(k)= k − µ − ln 1+ e T , π −∞ c2 + (k − q)2 Z (7) FIG. 2: Magnetization versus magnetic field H for the 1D which determines the thermodynamics of the system in Heisenberg chain with different values of ∆, indicating differ- a whole temperature regime. Using the Bethe ansatz ent magnetic ordering. The figure from [2]. equation (5) again, the Gibbs free energy per length p = ∂G − ∂L T,µ,c gives III. YANG-YANG THERMODYNAMICS T − p = ln 1+ e ε(k)/T dk. (8) 2π In 1969 Professor Chen Ping Yang in collaboration Z with Professor Chen Ning Yang [24] published his semi- Thus other thermodynamic quantities can be further de- nal work on the thermodynamics of the Lieb-Liniger Bose rived through thermodynamic relations, see [3]. gas. They proposed for the first time the grand canon- We arguably say that the TBA equation (7) provides ical ensemble description of the integrable model using a prototype of quantum statistical mechanics. It encodes the Bethe ansatz equations. Technically, the thermody- not only the quantum statistical effect but also the rich namics of the model is determined from the minimization quantum many-body dynamical interaction effect. As is conditions of the Gibbs free energy by using twice the mentioned in the commentary by Professor Chen Ning Bethe ansatz equations of the model. They started their Yang [28],“It shows the subtlety in the definitions of the formalism from the distribution function of the quasi- vacuum, the interaction, and the excitation spectrum”. momenta subject to the Bethe ansatz equation of the Moreover, Professor Chen Ping Yang [3] showed its con- Lieb-Linger model nection to both bosonic and fermionic statistics using the ∞ ′ ′ TBA equation. The particle and hole densities in terms h 1 c ρ(k )dk ρ(k)+ ρ (k)= + 2 ′ 2 , (5) of the quasimomenta reveal such subtle changes under a 2π π −∞ c + (k − k ) Z change of the temperature and interaction strength. where ρ(k) and ρh(k) are respectively the particle and Building on the Yang-Yang’s approach [3], Professor hole density distribution functions at finite temperatures. Minoru Takahashi has made further important contri- A brilliant consideration of the degeneracies of an equilib- butions to the thermodynamics of 1D integrable models rium state at finite temperatures was proposed by Yang [29, 30]. and Yang through the Fermi statistics in an interval of In view of the grand canonical ensemble, there exists the quasimomenta a quantum phase transition at the chemical potential µc = 0 at zero temperature. Yang and Yang’s TBA (L(ρ + ρh)dk)! dW = . (6) equation (7) provides a precise understanding of the uni- (Lρdk)!(Lρhdk)! versal thermodynamics, the quantum criticality and the quantum liquid in 1D Lieb-Liniger gas, see a review [31]. Consequently, Yang and Yang were able to give the ex- Ultracold bosonic atoms trapped in a quasi-1D geom- pression of entropy per unit length etry are ultimately related to the integrable models of ∞ quantum gases. Based on the TBA, particularly striking h ρ ρ s = (ρ + ρ ) ln 1+ − ρ ln dk. examples were the measurements of the thermodynamics −∞ ρh ρh Z and quantum fluctuations [32–40], the dynamic structure Here we would like to emphasize that such a subtle con- factor [41], the quantum criticality and the Tomonaga- nection between the Bethe ansatz microscopic states and Luttinger liquid (TLL) [42]. the macroscopic state of the system play a key role in the In a recent paper [42], the density profiles of quasi- Yang-Yang method. 1D trapped ultracold 87Rb atoms were measured by in Maximizing the entropy is the next key step in the situ absorption imaging. The density scaling law and Yang-Yang approach. The Gibbs free energy per unit the equation of states were obtained by rescaling these 4
Lett. 19, 588 (1967). [4] C. N. Yang and C. P. Yang, Phys. Rev. Lett. 13, 303 (1964). [5] C. P. Yang, Phys. Rev. A 2,154 (1970). [6] C. N. Yang and C. P. Yang, Phys. Lett. 20,9 (1966) C. P. Yang, Phys. Rev. Lett. 19, 586 (1967). Alfred C. T. Wu, C. P. Yang, Kurt Fuchel and Sidney Heller, Phys. Rev. Lett. 12, 57 (1964). Alfred C. T. Wu and C. P. Yang, Phys. Rev. D 1, 3180 (1970); B. D. Metcalf and C. P. Yang, Phys. Rev. B 18, 2304 (1978); C. N. Yang and C. P. Yang, Trans. NY Acad. Sci. 40, 267 (1980). [7] H. A. Bethe, Z. Physik 71, 205 (1931). [8] L. Hulthen´en, Arkiv. Mat. Mstron. Fysik 26 A, 11 (1938). [9] R. Orbach, Phys. Rev. 112, 309 (1958). [10] R. B. Griffiths, Phys. Rev. 113, A768 (1964). [11] L. R. Walker, Phys. Rev. 116, 1089 (1959). FIG. 3: Contour plot of the dimensionless specific heat [12] J. des Cloizeaux and J. J. Pearson, Phys. Rev. 128, 2131 c˜v = cv/(kB c) in T − µ plane. The dotes denote experimen- (1962). tal data of the specific heat peaks that mark two crossover [13] E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). temperatures. Whereas the dashed lines show the crossover [14] C. N. Yang, Phys. Rev. Lett. 19, 1312 (1967). temperatures predicted from the TBA equation (7). Other [15] M. Gaudin, Phys. Lett. A 24, 55 (1967). predictions on universal scaling functions of thermodynami- [16] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 cal properties and critical exponents from the TAB equation (1968). were found to agree very well with experimental observation. [17] B. Sutherland, Phys. Rev. Lett. 20, 98 (1968). From [42]. [18] R. J. Baxter, Ann. Phys. (N. Y.) 70, 193 (1972); R. J. Baxter, Ann. Phys. (N. Y.) 70, 323 (1972). [19] V. E. Korepin, A. G. Izergin, and N. M. Bogoliubov, Quantum Inverse Scattering Method and Correlation measurements at different temperatures and chemical po- Functions (Cambridge: Cambridge University Press), tentials. Based on the obtained equation of states, two 1993. crossover branches that distinguish the quantum critical [20] B. Sutherland, Beautiful Models: 70 years of exactly regime from the classical gas and the TLL were observed solved quantum many-body problems (Singapore: World through the double-peak structure of the specific heat, Scientific), 2004. see Fig. 3. Furthermore, the measured propagations of [21] M. Takahashi Thermodynamics of One-Dimensional density disturbances, the Luttinger parameters and also Solvable Models (Cambridge: Cambridge University Press), 1999. the power-law behavior in the momentum profiles con- [22] M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac and firm the existence of the TLL. The updated observations M. Rigol, Rev. Mod. Phys. 83 1405 (2011). of such many-body phenomena have revealed the beauty [23] X. W. Guan, M. T. Batchelor and C. Lee, Rev. Mod. of Yang and Yang’s grand canonical ensemble approach Phys. 85, 1633 (2013). to interacting many-body systems. [24] C. N. Yang and C. P. Yang, J. Math. Phys. 10, 1115 (1969). In summary, we have presented two significant works [25] I. Nolden, J. Stat. Phys. 67, 155 (1992). which Professor Chen Ping Yang did in the 60’s.Those [26] J. D. Shore and D. J. Buckman, Phys. Rev. Lett. 72, 604 works leads to significant applications in a wide range (1994); of problems in quantum statistical mechanics and math- D. J. Buckman and J. D. Shore, J. Stat. Phys. 78, 1227 ematical physics. His contributions are a remarkable (1995). legacy to physics. [27] H. Y. Huang, F. Y. Wu, H. Kunz and D. Kim, Physica A 228, 1 (1996). This work has been supported by the key NNSF C [28] C. N. Yang Selected papers 1945-1980 W. H. Freeman grant No. 11534014. and Company, 1983. [29] M. Takahashi, Prog. Theor. Phys. 46, 401 (1971); M. Takahashi, Prog. Theor. Phys. 46, 1388 (1971). [30] M. Takahashi, Prog. Theor. Phys. 47, 69 (1972); M. Takahashi, Prog. Theor. Phys. 50, 1519 (1973); ∗ e-mail:[email protected]; [email protected] M. Takahashi, Prog. Theor. Phys. 52, 103 (1973). [1] C. N. Yang and C. P. Yang, Phys. Rev. 147, 303 (1966). [31] Y.-Z. Jiang, Y.-Y. Chen and X.-W. Guan, Chinese Phys. [2] C. N. Yang and C. P. Yang, Phys. Rev. 150, 321 (1966); B 24, 050311 (2015). C. N. Yang and C. P. Yang, Phys. Rev. 150, 327 (1966); [32] A. H. van Amerongen, J. J. P. van Es, P. Wicke, K.V. C. N. Yang and C. P. Yang, Phys. Rev. 151, 258 (1966). Kheruntsyan, and N. J. van Druten, Phys. Rev. Lett. [3] S. Sutherland, C. N. Yang, and C. P. Yang, Phys. Rev. 100, 090402 (2008). 5
[33] J. Armijo, T. Jacqmin, K.V. Kheruntsyan, and I. Bou- (2010). choule, Phys. Rev. Lett. 105, 230402 (2010). [39] Y. Sagi, M. Brook, I. Almog, and N. Davidson, Phys. [34] J. Armijo, T. Jacqmin, K.V. Kheruntsyan, and I. Bou- Rev. Lett. 108, 093002 (2012). choule, Phys. Rev. A 83, 021605(R) (2011). [40] A. Vogler, R. Labouvie, F. Stubenrauch, G. Barontini, V. [35] J. Armijo, Phys. Rev. Lett. 108, 225306 (2012). Guarrera, and H. Ott. Phys. Rev. A 88, 031603 (2013). [36] T. Jacqmin, J. Armijo, T. Berrada, K.V. Kheruntsyan, [41] F. Meinert, M. Panfil, M. J. Mark, K. Lauber, J.-S. Caux, and I. Bouchoule, Phys. Rev. Lett. 106, 230405 (2011). and H.-C. N¨agerl, Phys. Rev. Lett. 115, 085301 (2015). [37] H.-P. Stimming, N. J. Mauser, J. Schmiedmayer, and I. [42] B. Yang, Y.-Y. Chen, Y.-G. Zheng, H. Sun, H.-N. Dai, E. Mazets, Phys. Rev. Lett. 105, 015301 (2010). X.-W. Guan, Z.-S. Yuan, and J.-W. Pan, Phys. Rev. Lett. [38] P. Kr¨uger, P., S. Hofferberth, I. E. Mazets, I. Lesanovsky, 119, 165701 (2017). and J. Schmiedmayer, Phys. Rev. Lett. 105, 265302