Paramagnetism and Diamagnetism
Total Page:16
File Type:pdf, Size:1020Kb
Paramagnetism and Diamagnetism Here we discuss the magnetic response of systems that do not feature (as opposed to, say, ferromagnets or superconductors) any order relevant to the interaction with the magnetic field. We will confine ourselves to the static response to a homogeneous magnetic field. Basic definitions The magnetization M is defined as the magnetic moment per unit vol- ume. In what follows, we will set system volume equal to unity so that the quantities per unit volume and per given volume V will coincide. (It is always straightforward to restore the volume, when necessary). The general statistical mechanical formula for M reads @F (V; N; T; B) M = − ; (1) @B where F (V; N; T; B) is the free energy as a function of volume, V , total number of particles, N, temperature, T , and the magnetic field, B. The magnetic susceptibility χ is defined as the ratio M χ = : (2) B The idea behind this definition is that in the limit B ! 0 the magnetization M is normally a linear function B; in this case, χ is a B-independent coeffi- cient. A system with χ > 0 are called paramagnetic and systems with χ < 0 are called diamagnetic. The response of paramagnetic systems is such that the external magnetic field gets enhanced while the response of diamagnetic systems reduces the external field. Absence of magnetic response in equilibrium classical statistics The equilibrium magnetic response is a purely quantum effect. The free energy of classical systems is fundamentally independent of magnetic field. Indeed, F = −T ln Z; and the partition function for a classical particle1 is given by Z Z = eH[p−(q=c)A(r); r] d3p d3r : 1The generalization to many particles is trivial. 1 Here H is the Hamiltonian function, p is the canonical momentum2 of the particle, r is the coordinate of the particle, A(r is the vector potential of the magnetic field, q is the electric charge of the particle, and c is the speed of light. The dependence of H on p and A is always through the variable p~ = p − (q=c)A(r) : This variable can be used as the new integration variable over canonical momentum, because, at any fixed r, it is nothing but a shift of the origin of p by the constant (q=c)A(r), and we get exactly the same expression as at A = 0: Z Z Z = eH[p~; r] d3p~ d3r ≡ eH[p ;r] d3p d3r : Two-level system: Paramagnetic response. Isentropic demagneti- zation Consider a macroscopic system formed by N independent two-level systems each of which having the Hamiltonian (written as the diagonal matrix in the representation of the energy eigenstates): ! µ B 0 H = ∗ ; (3) 0 −µ∗B where µ∗ is a certain constant characterizing the strength of the interaction with magnetic field. Without loss of generality, we set µ∗ > 0, since changing the sign of µ∗ leads to nothing but re-labeling the eigenvalues, which is irrelevant to the physical properties of the system. The partition function of a single two-level system is −E−=T −E+=T x −x Z0 = e + e = e + e ; where E∓ = ∓µ∗B and µ B x = ∗ : T The probabilities to find the system in the states `-' and `+' are 1 w = e−E∓=T =Z = : ∓ 0 e∓x + 1 2Not to be confused with mass times velocity, mv, which is true only in the absence of magnetic field; the relation between v and p is mv = p − (q=c)A. 2 Calculating @Z0 @F @ @B M = − = TN ln Z0 = TN ; @B @B Z0 we find M = Nµ∗ tanh x : This formula has a very transparent physical meaning: M = µ∗(N− − N+) ;N∓ = Nw∓ : The numbers N− and N+ are the numbers of two-level systems in the states `-' and `+' , respectively. In the limit of B ! 0 we have Nµ2B M ≈ ∗ (B ! 0) ; T and Nµ2 χ ≈ ∗ (B ! 0) : T The system demonstrates a paramagnetic response. Now consider the entropy of the system. By definition, the entropy of a single two-component system is S0 = −w− ln w− − w+ ln w+ ; and the entropy of the whole system is S = NS0. Note that S ≡ S(x), mean- ing that the dependence of S on B and T is only in terms their ratio. It is straightforward to dramatically reduce B without changing the probabilities w∓, because for the Hamiltonian (3) with time-dependent B the states `-' and `+' remain the eigenstates and their probabilities cannot change during the evolution. This means that the state remains an equilibrium state, but with the new temperature that simply follows the evolution of B, preserv- ing the ratio B=T . This process is called isentropic demagnetization. It is one of the standard cryogenic protocols. First one applies a magnetic field to polarize a two-level subsystem (normally, a spin subsystem). Then the magnetic field is quickly reduced to reduce the temperature of the two-level subsystem. The two-level subsystem then equilibrates with the rest of the system. The temperature of the two-level subsystem increases while the temperature of the rest of the system decreases. 3.