We study the correlation functions of one-dimensional Hubbard model in the presence of external magnetic field through the conformal field method. The long distance behaviour of the correlation functions and their unusual exponents for the model in the presence of a magnetic field are developed by solving the dressed charge matrix equations and setting the number of occupancies
to one, as alternative to the usual zero used by authors in literatures. This work shows that the exponent of the correlation functions is a monotonous function of magnetic field and the correlation functions decay as powers of these unusual exponents. As the magnetic field goes to zero, we obtain the exponents as 8.125, 11.125, 17.125, 26.125 and 38.125 at kF, 3kF, 5kF, 7kF
and 9kF. Our analytical results will provide insights into criticality in condensed matter physics.
Correlation Functions Magnetic Field Unusual Exponents1. Introduction
Almost twenty five years ago, Frahm and Korepin introduced the calculation of critical exponents for the one-dimensional (1D) Hubbard model, using the finite size scaling and the principle of conformal field theory (CFT) [1] . This enabled theorists to explore the physics of 1D correlated electron systems. Notwithstanding significant works, the understanding of the behaviour of correlated electron systems is not yet complete. In one dimension, the Hubbard Hamiltonian provides opportunity to study correlation effects in 1D models and the correlation functions decay as power of the distance [2] - [4] . It is the calculation of the critical exponents characterizing this power-law behaviour that have attracted constant theoretical interest. Outstanding results in this field
(with correlation exponents for at zero magnetic field around the Fermi points
) have been obtained from conformal field techniques, perturbation calculations and renormalization group methods in different models [1] [5] -[7] . For our calculation, we obtain the correlation exponents as
8.125, 11.125, 17.125, 26.125 and 38.125 around and by setting the parameter characterizing particle-hole excitation to one as the magnetic field goes to zero, and the unusual expo-
nent of the correlation function changes monotonically with change in magnetic field. The progress made in the understanding of critical phenomena in quantum systems as a result of conformal invariance have provided great insights to the problem of calculation of these critical exponents [8] . Although, interacting 1D quantum systems might carry countless low-energy excitations, with linear dispersion relations, but with different Fermi velocities, so the systems are not Lorentz invariant [9] . When the motions of these excitations are decoupled, one can now apply the CFT [10] . Usually, in the application of the conformal field techniques, the nonnegative integer characterizing particle-hole excitations is always taken as zero, but in this paper we shall calculate the electron
field correlation function and the density-density correlation function by setting the parameter to one, and
investigate how this affects the conformal dimensions and critical exponents of the correlation functions. This paper is organized as follows. In Section 2, we review the Bethe Ansatz equations of the Hubbard model and the analytic form of the correlation functions predicted by CFT is given. The dressed charge matrix elements are also calculated with the Wiener-Hopf technique and these elements are used to obtain the magnetic field dependence of the conformal dimensions. The long-distance behaviour of the electron field and density-density correlation functions and their unusual exponents for small magnetic field are calculated in Section 3. The electron field correlation function in momentum space and their Tomonaga-Luttinger (TL) liquid behaviour is examined in Section 4. Finally, Section 5 is devoted to discussion of the properties of the critical exponents for and conclusion.
2. The Hubbard Model and the Dressed Charge
The Hubbard model is basically the simplest model describing interacting spin-1/2 fermions in many-body physics. In the presence of magnetic field it is defined by the Hamiltonian [11]
where is the creation (annihilation) operator with electron spin σ at site j and is the
number operator. u is the on-site Coulomb repulsion, μ is the chemical potential and H is the external magnetic field. The hopping integral t = 1. Lieb and Wu [2] has solved Equation (1) exactly and obtained the Bethe Ansatz equations
where the quantum number Ij and Jα are integers or half-odd integer, with N↑ and N↓ being the number of electrons with spin up and down, and Ns = N↓ down spins are characterized by the moment a kj of holons and rapidities λα of spinons.
In the thermodynamic limit, with continuous momentum and rapidity variables, the Lieb-Wu equations become integral equations for the ground state distribution functions of moment a and of rapidities, obeying the equations
The state corresponding to the solution of Equations (2) and (3) has energy and momentum given by
where the conformal dimensions are given by
The positive integers, for holon and spinon describes particle-hole excitations, with being the number of occupancies that a particle at the right (left) Fermi level jumps to, represents the
change in the number of electrons (down-spin) with respect to the ground state, represents the number of particles which transfer from one Fermi level of the holon to the other and represents the number of particles which transfer from one Fermi level of the spinon to the other, and both and are either integer or half-odd integer values. Finally, the dressed charge matrix Z describing anomalous behaviour of critical exponents is given by
and the elements are defined by the solutions of the following coupled integral equations
where the kernel is defined as
The values of are fixed by
For small magnetic field we solve the dressed charge matrix equations by Wiener-Hopf technique [12] [13] for terms up to order in the strong coupling limit. With Equation (16), we write Equation (13) as
Fourier transforming Equation (17), we obtain
where the kernels are given by
We solve Equation (18) by introducing the function
and expanding it as
where are defined as the solutions of the Wiener-Hopf equations
The driving terms and the solutions becomes smaller as n increases because is large. Our procedure follows Fabian et al. [11] . Assuming the function and are known. We define
Where the functions are analytic on the upper and lower planes respectively, with
Also we assume
In terms of these functions we express the Fourier transform of Equation (23) as
where is the Fourier transform of. Now we split Equation (23) into the sum of two parts that are analytical and non-zero in the upper and lower half planes. To obtain this we use the factorization
where are analytic and non-zero in the upper and lower half planes respectively and are normalized as
Useful special function of are
Using Equations (27) and (28), we obtain
Decompose the right hand side of Equation (32) into the sum of two functions
This implies that
To obtain the solution of Equation (22) for, we set the driving term to be
We decompose the first term by using
The second term of Equation (35) is meromorphic function of with simple poles located at
Note, there is no pole at. The decomposition of the factor gives
Using Equation (39) we can express the function, for any function that is analytic and bounded in the lower half-plane as the sum of two functions analytic in the upper/lower half-plane
Applying the formula Equations (41) to (35) and Equation (33), we obtain
Now,
Therefore,
For
The functions are obtained by using Equation (34)
From Equation (23) for, by setting in Equation(48), we obtain
By definition
where is magnetic field, is critical field, strong coupling, magnetic field at zero temperature and corresponds to Fermi points. Combining the result Equation (49) with Equation (50), we obtain the first order contribution to as follows
As, we use Equations (30) and (31) on Equation (52) to obtain
Simplifying further, we obtain
Using Equation (51), we obtain
Next, the second order contribution to is obtained by taking the Fourier transform of Equation (23) for.
From Equation (28)
From Equation (33),
We have decomposed into which is analytic in the upper and lower half-planes. is given by
where is a small positive constant. has a branch cut along the negative imaginary axis and by deforming the contour of integration we rewrite Equation (58) as
From Equation (29), as
Since,
For the integrand rapidly decrease because, and hence the integral is approximated by expanding the terms other than around. Therefore, we obtain
From Equation (34), we obtain
Using
we obtain
Using the value of from Equation (51) in Equation (68), we obtain
Therefore, with Equations (55) and (69), we obtain
Now to evaluate the dressed charge matrix element we take the Fourier transform of Equations (16) and (18) and obtain
Applying the same process in the determination of Equation (70), we obtain
Similarly, with the same process, we obtain the other two elements of the dressed charge matrix as
and
From Equation (16) together with the property that for and , the down-spin density is obtained as
Using Equation (75) on Equations (70) and (72), we obtain the dressed charge matrix equations as
At half-filling, and by neglecting corrections to order, the elements of the dressed charge become
To obtain the conformal dimensions in terms of small magnetic field we use Equations (80) to (83) on Equations (7) and (8). Note that,
Therefore, the magnetic field dependence of the conformal dimensions are given by
According to the principles of CFT, the general expression for correlation function contains factors from holons and spinons, given by [11]
3. Correlation Functions in Magnetic Field
We now use the results obtained in the last section to obtain the magnetic field dependence of the unusual exponents of the electron field correlation function and density-density correlation function by setting the non-
negative integer. First we consider the electron field correlation function with up-spin which originates
from the quantum numbers (Dc, Ds) = (1/2, −1/2), (3/2, −3/2), (5/2, −5/2), (7/2, −7/2), (9/2, −9/2), ΔNc = 1 and ΔNs = 0. Therefore, the corresponding conformal dimensions for are
where the contributions from and terms of order are neglected. Using Equations (89) and (91) on Equation (87), we obtain
The critical exponent is given by
This implies that
and
Next, we obtain the conformal dimensions for as
Using Equations (97) and (99) on Equation (87), we obtain
The critical exponent is given by
and
Next, for, we obtain the conformal dimensions as
Using Equations (104) and (106) on Equation (87), we obtain
The critical exponent is given by
and
For, we obtain the conformal dimensions as
Using Equations (111) and (113) on Equation (87), we obtain
The critical exponent is given by
and
Finally, for, we obtain the conformal dimensions as
Using Equations (118) and (120) on Equation (87), we obtain
The critical exponent is given by
and
Combining Equations (92), (100), (107), (114) and (121), we obtain the long-distance asymptotic form of the electron field correlation function with up-spin as
Lastly, we consider the density-density correlation function which originates from the quantum numbers and. Here the corresponding conformal dimensions for are
Again contributions from are neglected. Using Equations (125) and (126) on Equation (87), we obtain
The critical exponents are given by
and
For the conformal dimensions are
Using Equations (130) and (131) on Equation (87), we obtain
The critical exponents are given by
and
Next, for the conformal dimensions are
Using Equations (135) and (136) on Equation (87), we obtain
and
Finally, for the conformal dimensions are
Using Equations (140) and (141) on Equation (87), we obtain
and
4. Correlation Function in Momentum Space
The electron field correlation function Equation (124) has singularities at the Fermi points and respectively. Therefore, at, the momentum distribution is given by
The critical exponent
and
Here we neglect logarithmic field dependence. Equation (145) represents the momentum distribution function around for the electron field correlator. It exhibits a typical power-law behaviour of the TL liquid, with critical exponent given by Equation (146). This unusual exponent as the magnetic field goes to zero.
Another singularity is at. The momentum distribution here is
with critical exponent
Equation (148) exhibits a typical power-law singularity of the TL liquid around the Fermi point, and as the magnetic field goes to zero.
Next at, the momentum distribution is given by
with the unusual exponent
Also, Equation (150) represents the momentum distribution function around the Fermi point for the electron field correlator, and it exhibits a typical power-law behaviour of the TL liquid with critical exponent given by Equation (151). This unusual exponent as the magnetic field goes to zero.
At, the momentum distribution is
and
Equation (152) exhibits a typical power-law behaviour of the TL liquid around, with critical exponent given by Equation (153). This unusual exponent as the magnetic field goes to zero.
Finally, at the momentum distribution takes the form
with
and Equation (154) also exhibits typical power-law behaviour of the TL liquid around the Fermi point, with critical exponent as the magnetic field goes to zero.
5. Discussions
In this paper, we have calculated the electron field and density-density correlation functions and their unusual exponents by using the nonnegative integer characterizing particle-hole excitations as 1 in the 1D Hubbard model. Based on the principles of CFT, we obtain expressions for the unusual exponents that describe the long-distance behaviour of the correlation functions in coordinate and momentum space. The unusual behaviour of the exponents depend on the magnetic field. The zero magnetic field exponents at and for in the strong-coupling limit has been obtained before [1] [3] [14] -[16] . At zero magnetic field is a monotonous function of the coupling constant u. In our calculation, we have used and obtain the exponents at and respectively, and observe that is a monotonous function of the magnetic field. The part arises
from the excitation of, the part from , the part from, the part from, and part from . This implies both holon and spinon excitations are responsible for and oscillation parts respectively.
In conclusion, the electron field correlation function and the unusual exponents has been obtained around the Fermi points and respectively, and the density-density correlation function around and is also obtained for. These results show that the correlation function clearly exhibits power-law behaviour of TL liquids as the exponent changes monotonically with change in magnetic field. Setting indicates presence of particle-hole excitations in the asymptotics. Therefore, the results obtained here are contributions from both particle-hole excitations and collective modes (holons and spinons).
ReferencesFrahm, H. and Korepin, V.E. (1990) Critical Exponents for the One-Dimensional Hubbard Model. Physical Review B, 42, 10553-10565. http://dx.doi.org/10.1103/PhysRevB.42.10553Lieb, E.H. and Wu, F.Y. (1968) Absence of Mott Transition in an Exact Solution of the Short-Range, One-Band Model in One Dimension. Physical Review Letters, 20, 1445-1448. http://dx.doi.org/10.1103/PhysRevLett.20.1445Woynarovich, F. (1989) Finite-Size Effects in a Non-Half-Filled Hubbard Chain. Journal of Physics A, 22, 4243-4256. http://dx.doi.org/10.1088/0305-4470/22/19/017Nenuwe, O.N. and Akpojotor, F. (2015) Power-Law Dependence of Correlation Functions in the Tomonaga-Luttinger Liquid. International Journal of Theoretical and Mathematical Physics, 5, 8-15. http://www.sapub.org/global/showpaperpdf.aspx?doi=10.5923/j.ijtmp.20150501.02Parola, A. and Sorella, S. (1990) Asymptotic Spin-Spin Correlations of the U→∞, One-Dimensional Hubbard Model. Physical Review Letters, 64, 1831-1834. http://dx.doi.org/10.1103/PhysRevLett.64.1831Finkel’shtein A.M. ,et al. (1977)Correlation Functions in One-Dimensional Hubbard Model 25, 73-76.Luther, A. and Peschel, I. (1975) Calculation of Critical Exponents in Two Dimensions from Quantum Field Theory in one Dimension. Physical Review B, 12, 3906. http://dx.doi.org/10.1103/physrevb.12.3908Belavin, A.A. Polyakov, A.M. and Zamolodchikov, A.B. (1984) Infinite Conformal Symmetry in Two-Dimensional Quantum Field Theory. Nuclear Physics B, 241, 333-380. http://dx.doi.org/10.1016/0550-3213(84)90052-XKawakami, N. and Yang, S.-K. (1991) Luttinger Liquid Properties of Highly Correlated Electron Systems in One Dimension. Journal of Physics: Condensed Matter, 3, 5983-6008. http://dx.doi.org/10.1088/0953-8984/3/32/007Izergin, A.G., Korepin, V.E. and Reshetikhin, Y. (1989) Conformal Dimensions in Bethe Ansatz Solvable Models. Journal of Physics A, 22, 2615-2620. http://dx.doi.org/10.1088/0305-4470/22/13/052Penc, K. and Solyom, J. (1993) One Dimensional Hubbard Model in a Magnetic Field and the Multicomponent Tomonaga-Luttinger Model. Physical Review B, 47, 6273-6292. http://dx.doi.org/10.1103/PhysRevB.47.6273Fabian, H.L.E., Frahm, H., Frank, G.O.H., Andreas, K. and Korepin, V.E. (2005) The One-Dimensional Hubbard Model. Cambridge University Press, New York, 1-674.Yang, C.N. and Yang, C.P. (1966) One-Dimensional Chain of Anisotropic Spin-Spin Interactions. II. Properties of the Ground-State Energy per Lattice Site for an Infinite System. Physical Review Letters, 150, 327-339. http://dx.doi.org/10.1103/PhysRev.150.327Kawakami, N. and Yang, S.-K. (1990) Luttinger Anomaly Exponent of Momentum Distribution in the Hubbard Chain. Physics Letters A, 148, 359-362. http://dx.doi.org/10.1016/0375-9601(90)90818-9Qin, S.J. and Yu, L. (1996) Momentum Distribution Critical Exponents for the One-Dimensional Large-U Hubbard Model in the Thermodynamic Limit. Physical Review B, 54, 1447-1450. http://dx.doi.org/10.1103/PhysRevB.54.1447Qin, S.J., Liang, S.D., Su, Z.B. and Yu, L. (1995) Density-Matrix Renormalization-Group Calculation of Correlation Functions in the One-Dimensional Hubbard Model. Physical Review B, 52, R5475-R5478. http://dx.doi.org/10.1103/PhysRevB.52.R5475