Abstract
A Green's-function (GF) theory which is different from the usual two-time GF theory (canonically averaged products of the time-dependent spin operators) is given for the system composed of cubic crystalline fields, Zeeman energy for an applied field, and the exchange energy. When interactions between excited spin waves are neglected, the GF's can be solved. The 2J poles or energies of the GF correspond to 2J branches of excited spin-wave modes where J is the total angular momentum of the rare-earth ion. The spontaneous magnetization M, the specific heat Cv, the parallel susceptibility χ, the second-order anisotropy constant Ku, the fourth-order cubic anisotropy constant K1 and the sixth-order cubic anisotropy constant K2 are evaluated by operating the corresponding spin operators to the GF. At very low temperatures and in the case of very small crystalline anistropy fields, KuKu(T=0), K1K1(T=0), and K2K2(T=0) are approximately proportional to {MM(T=0)}n2, {MM(T=0)}n4, and {MM(T=0)}n6, respectively, in which n2=3(J02Je2)(3J02J2){1(JeJ0)}, n4=5{7(J04Je4)6J2(J02Je2)}(35J0430J2J02+3J4)×{1(JeJ0)}, and n6=21{11(J06Je6)15J2(J04Je4)+5J4(J02Je2)}(231J06315J2J04+105J4J025J6){1(JeJ0)}, where Jen and J0n mean the magnetic quantum numbers (Jz)n averaged over the first excited state and the ground state, respectively. They reduce to the well-known third, tenth, and twenty-first laws when cubic crystalline fields disappear. All the thermodynamical properties mentioned above change rapidly at temperatures corresponding to excitation spectra of the spin-wave modes. In the case that the rare-earth ion has two levels, the formulas for M, χ, and Cv become identical with the expressions previously reported by the author and others.