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

arXiv:cond-mat/9807266 · Physical review. B, Condensed matter

δ-function Bose-gas picture of S=1 antiferromagnetic quantum spin chains near critical fields

Kouichi Okunishi, Yasuhiro Hieida, Yasuhiro Akutsu

Published Mar 1, 1999• 51 citations
arXivHTMLPDF

Abstract

We study the zero-temperature magnetization curve $(M\ensuremath{-}H$ curve) of the one-dimensional quantum antiferromagnet of spin one. The Hamiltonian H we consider is of the bilinear-biquadratic form: $H={\ensuremath{\sum}}_{i}f({s}_{i}\ensuremath{\cdot}{s}_{i+1})$ (+Zeeman term) where ${s}_{i}$ is the spin operator at site i and $f(X)=X+\ensuremath{\beta}{X}^{2}$ with $0<~\ensuremath{\beta}<1.$ We focus on validity of the $\ensuremath{\delta}$-function Bose-gas picture near the two critical fields: upper-critical field ${H}_{s}$ above which the magnetization saturates and the lower-critical field ${H}_{c}$ associated with the Haldane gap. As for the behavior near ${H}_{s},$ we take ``low-energy effective S matrix'' approach, where the correct effective Bose-gas coupling constant c is extracted from the two down-spin S matrix in its low-energy limit. We find that the resulting value of c differs from the spin-wave value. We draw the $M\ensuremath{-}H$ curve by using the resultant Bose gas, and compare it with numerical calculation where the product-wave-function renormalization-group (PWFRG) method, a variant of White's density-matrix renormalization group method, is employed. Excellent agreement is seen between the PWFRG calculation and the correctly mapped Bose-gas calculation. We also test the validity of the Bose-gas picture near the lower-critical field ${H}_{c}.$ Comparing the PWFRG-calculated $M\ensuremath{-}H$ curves with the Bose-gas prediction, we find that there are two distinct regions, I and II, of $\ensuremath{\beta}$ separated by a critical value ${\ensuremath{\beta}}_{c}(\ensuremath{\approx}0.41).$ In region I, $0<\ensuremath{\beta}<{\ensuremath{\beta}}_{c},$ the effective Bose coupling c is positive but rather small. The small value of c makes the ``critical region'' of the square-root behavior $M\ensuremath{\sim}\sqrt{H\ensuremath{-}{H}_{c}}$ very narrow. Further, we find that in the $\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\beta}}{\ensuremath{\beta}}_{c}\ensuremath{-}0,$ the square-root behavior transmutes to a different one, $M\ensuremath{\sim}(H\ensuremath{-}{H}_{c}{)}^{\ensuremath{\theta}}$ with $\ensuremath{\theta}\ensuremath{\approx}1/4.$ In region II, ${\ensuremath{\beta}}_{c}<\ensuremath{\beta}<1,$ the square-root behavior is more pronounced as compared with region I, but the effective coupling c becomes negative.

DOI
10.1103/physrevb.59.6806

Related papers

Magnetization process of bilinear-biquadratic spin chains at finite temperature

Kouichi Okunishi

We investigate magnetization process $(M\ensuremath{-}H\phantom{\rule{0ex}{0ex}}\mathrm{c}\mathrm{u}\mathrm{r}\mathrm{v}\mathrm{e})$ of $S=1$ bilinear-biquadratic (BLBQ) spin chain near critical fields at finite temperatures. We use the density matrix renormalization group (DMRG) for the two-dimensional classical lattice model mapped by the Trotter decomposition, with a help of the Baxter's variational principle. By comparing the DMRG result of $S=1/2 \mathrm{XY}$ chain with the exact solution, we show that the DMRG is an efficient tool to calculate the $M\ensuremath{-}H$ curve at finite temperatures. Further, we compute the $M\ensuremath{-}H$ curve of the BLBQ chain. We compare the DMRG curve of the magnetization process of the BLBQ chain with those obtained by analytic approaches: correctly mapped $\ensuremath{\delta}$-function Bose-gas approach and ``Bethe-ansatz approximation'' approach. Near the saturation field, we show that the $\ensuremath{\delta}$-function Bose gas and the Bethe-ansatz approximation describe the $M\ensuremath{-}H$ curve well in both at zero temperature and finite temperatures. Near the lower critical field, we find that the $\ensuremath{\delta}$-function Bose gas is a good effective model at the Heisenberg point. We further find that the $\ensuremath{\delta}$-function Bose gas cannot describe the $M\ensuremath{-}H$ curve at finite temperatures, near the special point where the $M\ensuremath{-}H$ curve changes qualitatively.

Physical review. B, Condensed matter•Aug 1, 1999•17 citations

Middle-field cusp singularities in the magnetization process of one-dimensional quantum antiferromagnets

Kouichi Okunishi, Yasuhiro Hieida, Yasuhiro Akutsu

We study the zero-temperature magnetization process $(M\ensuremath{-}H$ curve) of one-dimensional quantum antiferromagnets using a variant of the density-matrix renormalization group method. For both the $S=1/2$ zig-zag spin ladder and the $S=1$ bilinear-biquadratic chain, we find clear cusp-type singularities in the middle-field region of the $M\ensuremath{-}H$ curve. These singularities are successfully explained in terms of the double-minimum shape of the energy dispersion of the low-lying excitations. For the $S=1/2$ zig-zag spin ladder, we find that the cusp formation accompanies the Fermi-liquid to non-Fermi-liquid transition.

arXiv:cond-mat/9904155•Physical review. B, Condensed matter•Sep 1, 1999•58 citations

Antiferromagnetic zigzag spin chain in magnetic fields at finite temperatures

Nobuya Maeshima, Kouichi Okunishi

We study thermodynamic behaviors of the antiferromagnetic zigzag spin chain in magnetic fields, using the density-matrix renormalization-group method for the quantum transfer matrix. We focus on the thermodynamics of the system near the critical fields in the ground-state magnetization process $(M\ensuremath{-}H$ curve): the saturation field, the lower critical field associated with excitation gap, and the field at the middle-field cusp singularity. We calculate magnetization, susceptibility, and specific heat of the zigzag chain in magnetic fields at finite temperatures, and then discuss how the calculated quantities reflect the low-lying excitations of the system related with the critical behaviors in the $M\ensuremath{-}H$ curve.

arXiv:cond-mat/0004159•Physical review. B, Condensed matter•Jul 1, 2000•51 citations

XXZ spin- 12 representation of a finite- U Bose-Hubbard chain at half-integer filling

Domenico Giuliano, Davide Rossini, Pasquale Sodano +1 more

Using a similarity Hamiltonian renormalization procedure, we determine an effective spin-$\frac{1}{2}$ representation of the Bose-Hubbard model at half-integer filling and at a finite onsite interaction energy $U$. By means of bosonization, we are able to recast the effective Hamiltonian as that of a spin-$\frac{1}{2}$ XXZ magnetic chain with pertinently renormalized coupling and anisotropy parameters. We use this mapping to provide analytical estimates of the correlation functions of the Bose-Hubbard model. We then compare such results with those based on density matrix renormalization group numerical simulations of the Bose-Hubbard model for various values of $U$ and for a number $L$ of lattice sites as low as $L\ensuremath{\sim}30$. We find an excellent agreement up to $10%$ between the output of analytical and numerical computations, even for relatively small values of $U$. Our analysis implies that, also at finite $U$, the one-dimensional Bose-Hubbard model with suitably chosen parameters may be seen as a quantum simulator of the XXZ chain.

arXiv:1210.2685•Physical Review B•Jan 3, 2013•21 citations

One-dimensional frustrated Heisenberg antiferromagnet of spin 1/2 in a strong magnetic field

S. Gluzman

The properties of the one-dimensional antiferromagnet of spin 1/2 in a strong magnetic field h with nearest- and next-nearest-neighbor interactions are studied in terms of the equivalent Bose-gas problem. We find that below the first critical field ${\mathit{h}}_{\mathit{c}\mathit{u}}$ the phase with a gap in the elementary excitations spectrum and corresponding plateau in the magnetization curve should exist due to the two-particle attraction in the vicinity of the commensurate wave number ${\mathit{p}}_{0\mathit{u}}$=\ensuremath{\pi}. In the interval between the first and second critical fields ${\mathit{h}}_{\mathit{c}\mathit{u}}$h${\mathit{h}}_{\mathit{c}\mathit{f}}$, the system behaves as the one-dimensional Bose gas with particle repulsion and the magnetization has a typical square-root dependence on field.

Physical review. B, Condensed matter•May 1, 1994•6 citations

Phase boundary of Ising antiferromagnets near H=Hc and T=0 : Results from hard-core lattice gas calculations

Zoltàn Ràcz

Ising antiferromagnets in a near-critical magnetic field at low temperatures are equivalent to hard-core lattice gases. Using this connection and the existing series-expansion results for hard-core lattice gases, we determine the slope of the phase boundary at $T=0$ for the square (sq), plane-triangular (pt), simple cubic (sc), and body-centered cubic (bcc) antiferromagnets. The slope is negative for the sq and pt lattices and nearly zero for the sc case. For the bcc lattice a positive slope is obtained, indicating that the phase boundary bulges above the zero-temperature critical field. We also test M\"uller-Hartmann and Zittarz's postulate for the critical curve of the sq Ising antiferromagnet. A renormalization-group treatment of the hard-square lattice gas yields a critical activity ${z}^{*}=3.7959\ifmmode\pm\else\textpm\fi{}0.0001$, which is in agreement with series-expansion and finite-lattice estimates but at variance with the postulated ${z}^{*}=4$. The same calculation gives $\ensuremath{\nu}=0.999\ifmmode\pm\else\textpm\fi{}0.001$ for the correlation-length exponent, thus supporting the conjecture that the transition of the hard-square lattice gas belongs to the Ising universality class.

Physical review. B, Condensed matter•May 1, 1980•90 citations
We study the zero-temperature magnetization curve (M\ensuremath−H(M\ensuremath{-}H(M\ensuremath−H curve) of the one-dimensional quantum antiferromagnet of spin one. The Hamiltonian H we consider is of the bilinear-biquadratic form: H=\ensuremath∑if(si\ensuremath⋅si+1)H={\ensuremath{\sum}}_{i}f({s}_{i}\ensuremath{\cdot}{s}_{i+1})H=\ensuremath∑i​f(si​\ensuremath⋅si+1​) (+Zeeman term) where si{s}_{i}si​ is the spin operator at site i and f(X)=X+\ensuremathβX2f(X)=X+\ensuremath{\beta}{X}^{2}f(X)=X+\ensuremathβX2 with 0< \ensuremathβ<1.0<~\ensuremath{\beta}<1.0< \ensuremathβ<1. We focus on validity of the \ensuremathδ\ensuremath{\delta}\ensuremathδ-function Bose-gas picture near the two critical fields: upper-critical field Hs{H}_{s}Hs​ above which the magnetization saturates and the lower-critical field Hc{H}_{c}Hc​ associated with the Haldane gap. As for the behavior near Hs,{H}_{s},Hs​, we take low-energy effective S matrix'' approach, where the correct effective Bose-gas coupling constant c is extracted from the two down-spin S matrix in its low-energy limit. We find that the resulting value of c differs from the spin-wave value. We draw the $M\ensuremath{-}H$ curve by using the resultant Bose gas, and compare it with numerical calculation where the product-wave-function renormalization-group (PWFRG) method, a variant of White's density-matrix renormalization group method, is employed. Excellent agreement is seen between the PWFRG calculation and the correctly mapped Bose-gas calculation. We also test the validity of the Bose-gas picture near the lower-critical field ${H}_{c}.$ Comparing the PWFRG-calculated $M\ensuremath{-}H$ curves with the Bose-gas prediction, we find that there are two distinct regions, I and II, of $\ensuremath{\beta}$ separated by a critical value ${\ensuremath{\beta}}_{c}(\ensuremath{\approx}0.41).$ In region I, $0<\ensuremath{\beta}<{\ensuremath{\beta}}_{c},$ the effective Bose coupling c is positive but rather small. The small value of c makes the critical region'' of the square-root behavior M\ensuremath∼H\ensuremath−HcM\ensuremath{\sim}\sqrt{H\ensuremath{-}{H}_{c}}M\ensuremath∼H\ensuremath−Hc​​ very narrow. Further, we find that in the \ensuremathβ\ensuremath→\ensuremathβc\ensuremath−0,\stackrel{\ensuremath{\rightarrow}}{\ensuremath{\beta}}{\ensuremath{\beta}}_{c}\ensuremath{-}0,\ensuremathβ\ensuremath→​\ensuremathβc​\ensuremath−0, the square-root behavior transmutes to a different one, M\ensuremath∼(H\ensuremath−Hc)\ensuremathθM\ensuremath{\sim}(H\ensuremath{-}{H}_{c}{)}^{\ensuremath{\theta}}M\ensuremath∼(H\ensuremath−Hc​)\ensuremathθ with \ensuremathθ\ensuremath≈1/4.\ensuremath{\theta}\ensuremath{\approx}1/4.\ensuremathθ\ensuremath≈1/4. In region II, \ensuremathβc<\ensuremathβ<1,{\ensuremath{\beta}}_{c}<\ensuremath{\beta}<1,\ensuremathβc​<\ensuremathβ<1, the square-root behavior is more pronounced as compared with region I, but the effective coupling c becomes negative.
We investigate magnetization process $(M\ensuremath{-}H\phantom{\rule{0ex}{0ex}}\mathrm{c}\mathrm{u}\mathrm{r}\mathrm{v}\mathrm{e})$ of S=1S=1S=1 bilinear-biquadratic (BLBQ) spin chain near critical fields at finite temperatures. We use the density matrix renormalization group (DMRG) for the two-dimensional classical lattice model mapped by the Trotter decomposition, with a help of the Baxter's variational principle. By comparing the DMRG result of S=1/2XYS=1/2 \mathrm{XY}S=1/2XY chain with the exact solution, we show that the DMRG is an efficient tool to calculate the $M\ensuremath{-}H$ curve at finite temperatures. Further, we compute the $M\ensuremath{-}H$ curve of the BLBQ chain. We compare the DMRG curve of the magnetization process of the BLBQ chain with those obtained by analytic approaches: correctly mapped $\ensuremath{\delta}$-function Bose-gas approach and ``Bethe-ansatz approximation'' approach. Near the saturation field, we show that the $\ensuremath{\delta}$-function Bose gas and the Bethe-ansatz approximation describe the $M\ensuremath{-}H$ curve well in both at zero temperature and finite temperatures. Near the lower critical field, we find that the $\ensuremath{\delta}$-function Bose gas is a good effective model at the Heisenberg point. We further find that the $\ensuremath{\delta}$-function Bose gas cannot describe the $M\ensuremath{-}H$ curve at finite temperatures, near the special point where the $M\ensuremath{-}H$ curve changes qualitatively.
We study the zero-temperature magnetization process $(M\ensuremath{-}H$ curve) of one-dimensional quantum antiferromagnets using a variant of the density-matrix renormalization group method. For both the S=1/2S=1/2S=1/2 zig-zag spin ladder and the S=1S=1S=1 bilinear-biquadratic chain, we find clear cusp-type singularities in the middle-field region of the $M\ensuremath{-}H$ curve. These singularities are successfully explained in terms of the double-minimum shape of the energy dispersion of the low-lying excitations. For the S=1/2S=1/2S=1/2 zig-zag spin ladder, we find that the cusp formation accompanies the Fermi-liquid to non-Fermi-liquid transition.
We study thermodynamic behaviors of the antiferromagnetic zigzag spin chain in magnetic fields, using the density-matrix renormalization-group method for the quantum transfer matrix. We focus on the thermodynamics of the system near the critical fields in the ground-state magnetization process $(M\ensuremath{-}H$ curve): the saturation field, the lower critical field associated with excitation gap, and the field at the middle-field cusp singularity. We calculate magnetization, susceptibility, and specific heat of the zigzag chain in magnetic fields at finite temperatures, and then discuss how the calculated quantities reflect the low-lying excitations of the system related with the critical behaviors in the $M\ensuremath{-}H$ curve.
Using a similarity Hamiltonian renormalization procedure, we determine an effective spin-12\frac{1}{2}21​ representation of the Bose-Hubbard model at half-integer filling and at a finite onsite interaction energy UUU. By means of bosonization, we are able to recast the effective Hamiltonian as that of a spin-12\frac{1}{2}21​ XXZ magnetic chain with pertinently renormalized coupling and anisotropy parameters. We use this mapping to provide analytical estimates of the correlation functions of the Bose-Hubbard model. We then compare such results with those based on density matrix renormalization group numerical simulations of the Bose-Hubbard model for various values of UUU and for a number LLL of lattice sites as low as $L\ensuremath{\sim}30$. We find an excellent agreement up to 1010%10 between the output of analytical and numerical computations, even for relatively small values of UUU. Our analysis implies that, also at finite UUU, the one-dimensional Bose-Hubbard model with suitably chosen parameters may be seen as a quantum simulator of the XXZ chain.
The properties of the one-dimensional antiferromagnet of spin 1/2 in a strong magnetic field h with nearest- and next-nearest-neighbor interactions are studied in terms of the equivalent Bose-gas problem. We find that below the first critical field hcu{\mathit{h}}_{\mathit{c}\mathit{u}}hcu​ the phase with a gap in the elementary excitations spectrum and corresponding plateau in the magnetization curve should exist due to the two-particle attraction in the vicinity of the commensurate wave number p0u{\mathit{p}}_{0\mathit{u}}p0u​=\ensuremath{\pi}. In the interval between the first and second critical fields hcu{\mathit{h}}_{\mathit{c}\mathit{u}}hcu​hhcf{\mathit{h}}_{\mathit{c}\mathit{f}}hcf​, the system behaves as the one-dimensional Bose gas with particle repulsion and the magnetization has a typical square-root dependence on field.
Ising antiferromagnets in a near-critical magnetic field at low temperatures are equivalent to hard-core lattice gases. Using this connection and the existing series-expansion results for hard-core lattice gases, we determine the slope of the phase boundary at T=0T=0T=0 for the square (sq), plane-triangular (pt), simple cubic (sc), and body-centered cubic (bcc) antiferromagnets. The slope is negative for the sq and pt lattices and nearly zero for the sc case. For the bcc lattice a positive slope is obtained, indicating that the phase boundary bulges above the zero-temperature critical field. We also test M\"uller-Hartmann and Zittarz's postulate for the critical curve of the sq Ising antiferromagnet. A renormalization-group treatment of the hard-square lattice gas yields a critical activity ${z}^{*}=3.7959\ifmmode\pm\else\textpm\fi{}0.0001$, which is in agreement with series-expansion and finite-lattice estimates but at variance with the postulated z∗=4{z}^{*}=4z∗=4. The same calculation gives $\ensuremath{\nu}=0.999\ifmmode\pm\else\textpm\fi{}0.001$ for the correlation-length exponent, thus supporting the conjecture that the transition of the hard-square lattice gas belongs to the Ising universality class.