Corresponding author. E-mail: yinxiangl@hotmail.com
Project supported by the National Natural Science Foundation of China (Grant No. 10774035) and the Qianjiang RenCai Program of Zhejiang Province, China (Grant No. 2007R0010).
We use the Schwinger-boson approach to study the anisotropy ferrimagnetic spin-(1/2,1) chain with bond alternation. Based on the effect of bond alternation δ, we obtain energy gap, free energy, and specific heat, respectively. The specific heat with larger bond alternation ( δ > 0.7) displays a peak at low temperature. Based on the effect of XXZ anisotropy parameter Δ, we present excited spectrums, free energy, and specific heat, respectively.
In the past decades, the one-dimensional (1D) magnetic system, for example, the bimetallic compound AMn(S2C2O2)2(H2O)3· 4.5H2O (A = Cu, Ni, Pd, Pt), has been one of the active fields in condensed matter. Many theoretical and experimental researchers have concentrated on the Haldane[1] conjecture that the ground state of a half-odd integer spin chain is gapless and spin correlations decay in the power law, while the ground state of an integer spin chain is a gap and spin correlations decay exponentially. Then, a number of studies were concerned with the one-dimensional mixed spin chain with antiferromagnetic interaction.[2– 6] The ground state of an isotropic mixed spin chain exhibits both ferromagnetic and antiferromagnetic elementary excitations. The excitation of the ferromagnetic aspect is gapless whereas the excitation of the antiferromagnetic aspect has a gap. Various methods have been developed to study the ground states and thermodynamic properties of the Heisenberg alternated spin chain. Pati[7] studied the ground and excited states of uniform and dimerized alternating spin chain by using a density-matrix renormalization-group (DMRG) method and spin wave theory. Abolfath[8] used a numerical exact diagonalization technique to study the phase diagram of one-dimensional quantum ferrimagnets. Sakai[9] investigated the quantum magnetization plateaux of an anisotropic ferrimagnetic spin chain by using size-scaling analyses. Yamamoto[10, 11] used the modified spin wave theory to investigate the thermodynamic properties of a Heisenberg ferrimagnetic spin chain and also studied the magnetic properties of quantum ferrimagnetic spin chains by the quantum Monte Carlo method. Langari[12] used DMRG and exact diagonalization Lanczos methods to study the phase diagram of the XXZ ferrimagnetic chain in the presence of a transverse magnetic field. Chen[13] applied self-consistent mean-field theory and the Green function method to research the ground states of a ferrimagnetic spin chain. In this paper, we study the excitation states and thermodynamic properties of an XXZ ferrimagnetic spin-(1/2, 1) chain with bond alternation by using the Schwinger-boson approach. This paper is organized as follows. In the second section, we present the Schwinger-boson approach to the anisotropy ferrimagnetic spin-(1/2, 1) chain with bond alternation. In the third section, we obtain the energy gaps, zero temperature free energy, and specific heat on the effect of bond alternation δ at Heisenberg point Δ = 1. In the fourth section, we discuss properties such as the energy spectrum, free energy, and specific heat based on the effect of the anisotropic parameter Δ . The conclusion and discussion are made in the final section.
The Hamiltonian of an XXZ ferrimagnetic spin-(1/2, 1) chain with bond alternation can be expressed as
where s = 1/2 and S = 1 are the spin of sublattice A and sublattice B respectively, Δ is the XXZ anisotropic parameter, δ is the bonding alternation and (s · S)Δ = sxSx + sySy + Δ szSz. The spin operator si can be represented by Schwinger bosons ai, ↑ and ai, ↓ :
The spin operator Sj can be represented by Schwinger bosons bj, ↑ and bj, ↓ :
The Hamiltonian can be written as
where
and
To implement the constraint
where
R2i is the coordinate of S2i. According to the diagonalization theory of Boson systems, the Heisenberg equation of motion can be written as
where
At Heisenberg point Δ = 1, the parameters are e = 0 and f = 0. The dynamic matrix can be rewritten as
The energy spectrums can be obtained as
where
The mean-field self-consistent equations can be calculated by minimizing the free energy (i.e., δ F/δ B* = 0, δ F/λ A = 0, δ F/λ B = 0). An infinitesimal external stagger magnetic field was supposed to be upward at site B and downward at site A at zero temperature. Due to the Bose– Einstein condensation, the mean-field self-consistent equations can be obtained as
By solving the above equations numerically, the energy gap with different bond alternation δ at zero temperature are plotted as shown in Fig. 1. The energy gap is monotonically increased while δ is added from 0 to 1. The zero temperature free energy per unit cell is
Due to the degree of freedom being overcounted by a factor of 2, as pointed out by Arovas and Auerbach, [14, 15] the free energy formula is calculated by dividing the part of fluctuation per unit cell F/N + 2sASB by 2 and adding back the classic ground energy per cell − 2sASB. The zero temperature free energy with different bond alternation δ as shown in Fig. 2. The free energy is monotonically decreased with the increase of δ . At finite temperature, the first equation of formula (12) can be modified as
Wu[16] pointed out that the SBMFT is not useful to describe the thermodynamic properties of the ferrimagnetic spin chain at intermediate and high temperatures. Thus, we discuss the thermodynamic properties of this model with T < 0.5. The specific heat with bond alternation versus temperature is calculated and plotted in Fig. 3. The effect of bond alternation δ on a specific heat is shown clearly at low temperature. The specific heat with a larger bond alternation parameter increases more rapidly and exhibits a peak when δ > 0.7 at lower temperature. The peak moves left while increasing the bond alternation δ . This means that the effect of bond alternation causes this abnormal phenomenon of specific heat at low temperature. When the ferrimagnetic spin chain with uniform bond alternation δ = 0, all the results mentioned above equal the results of Wu.[16]
By using the eigenvalue equation Dψ = Eψ , we obtain the energy spectrum
where
and
where Z is the number of nearest neighbor sites. The mean-field self-consistent equations can be obtained by minimizing the free energy (i.e., δ F/δ A* = 0, δ F/δ B* = 0, δ F/λ A = 0, δ F/λ B = 0). At zero temperature, the self-consistent equations are written as follows:
where
By solving the above self-consistent equations numerically, the excited energy spectrums with different anisotropic parameter Δ can be plotted respectively as shown in Fig. 4. The ferromagnetic branch increases more rapidly with smaller XXZ anisotropy from k = 0. The antiferromagnetic branch
where
By solving the new self-consistent equations numerically, the energy gaps with different anisotropic parameter Δ versus temperature are plotted as Fig. 5. It can be seen from the figure that the free energy with smaller anisotropy is bigger at low temperature. The free energy curve with Δ = 1 decreases relative rapidly while that with smaller anisotropy, Δ < 1, descends little and slowly. The specific heat at low temperature is plotted in Fig. 6. It is found that the effect of XXZ anisotropy Δ is not obviously relative, compared with the effect of bond alternation δ . With increasing temperature, the curve with different anisotropy also increases with temperature. The specific heat with smaller anisotropy is reduced at the given temperature.
By using the Schwinger-boson mean-field theory, the XXZ ferrimagnetic spin-(1/2, 1) chain with bond alternation is explored. We present the energy gap, zero temperature free energy, and specific heat of an isotropic ferrimagnetic spin chain with bond alternation δ . Based on the effect of XXZ anisotropy Δ , the energy spectrums and thermodynamic properties at T < 0.4 are calculated respectively. The relation between the properties of a ferrimagnetic spin chain and XXZ anisotropy or bond alternation are shown clearly in the former figures.
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|