<<

Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 1

Collective Effects

in

Equilibrium and Nonequilibrium Physics

Website: http://cncs.bnu.edu.cn/mccross/Course/ Caltech Mirror: http://haides.caltech.edu/BNU/

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 2

Today’s Lecture

Superfluids and superconductors

• What are superfluidity and ? • Review of phase dynamics • Description in terms of a macroscopic phase • Supercurrents that flow for ever • Josephson effect • Four sounds

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 3

The Amazing World of Superfluidity and Superconductivity

• Electric currents in loops that flow for ever (measured for ∼ decade) • Beakers of fluid that empty them- selves • Fluids that flow without resistance through tiny holes • Flow in surface films less than an atomic layer thick • Flow driven by temperature differ- ences (fountain effect)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 4

History of Superfluidity and Superconductivity 1908 Liquefaction of 4He by Kamerlingh Onnes

1911 Discovery of superconductivity by Onnes (resistance drops to zero)

1933 : superconductors expel magnetic field

1937 Discovery of superfluidity in 4He by Allen and Misener

1938 Connection of superfluidity with Bose-Einstein condensation by London

1955 Feynman’s theory of quantized vortices

1956 Onsager and Penrose identify the broken symmetry in superfluidity ODLRO

1957 BCS theory of superconductivity

1962 Josephson effect

1973 Discovery of superfluidity in 3He at 2mK by Osheroff, Lee, and Richardson

1986 Discovery of high-Tc superconductors by Bednorz and Müller 1995- Study of superfluidity in ultracold trapped dilute

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 5

Review of Phase Dynamics with a Conserved Quantity

Rotational symmetry in the XY plane (angle 2)

Sz The XY and Z components of the spin have different prop- erties: S⊥ P −1 Θ Sz =  hsizi is a conserved quantity i inP −1 S⊥ =  hsi⊥i is the XY order parameter i in 

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 5

Review of Phase Dynamics with a Conserved Quantity

Rotational symmetry in the XY plane (angle 2)

Sz The XY and Z components of the spin have different prop- erties: S⊥ P −1 Θ Sz =  hsizi is a conserved quantity i inP −1 S⊥ =  hsi⊥i is the XY order parameter i in 

Sz and 2 are canonically conjugate variables, so that with the free energy Z " # 1 S2 F = dd x K (∇2)2 + z − S b 2 2χ z z

we get δF S˙ =− giving S˙ =−∇·j with j =−K ∇2 z δ2 z Sz Sz δ ˙ F ˙ −1 2 = giving 2 = χ (Sz − χbz) δSz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 6

Spin Current ˙ =−∇· Sz jSz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 6

Spin Current ˙ =−∇· Sz jSz

This is a conservation law with a current jSz of the conserved quantity Sz given by a phase gradient =− ∇2 jSz K

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 6

Spin Current ˙ =−∇· Sz jSz

This is a conservation law with a current jSz of the conserved quantity Sz given by a phase gradient =− ∇2 jSz K

For example

Θ(r)

L

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

• No dynamics in full thermodynamic equilibrium: Sz = χb0z

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

• No dynamics in full thermodynamic equilibrium: Sz = χb0z

• Add an additional external field b1z = γ B1z ˙ 2 =−b1z =−γ B1z

the usual precession of a magnetic moment in an applied field (Larmor precession).

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

• No dynamics in full thermodynamic equilibrium: Sz = χb0z

• Add an additional external field b1z = γ B1z ˙ 2 =−b1z =−γ B1z

the usual precession of a magnetic moment in an applied field (Larmor precession).

• Note that this is an equilibrium state: Sz 6= χ(b0z + b1z) but is a conserved quantity

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

• No dynamics in full thermodynamic equilibrium: Sz = χb0z

• Add an additional external field b1z = γ B1z ˙ 2 =−b1z =−γ B1z

the usual precession of a magnetic moment in an applied field (Larmor precession).

• Note that this is an equilibrium state: Sz 6= χ(b0z + b1z) but is a conserved quantity— no approximations

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 7

Phase Dynamics

˙ −1 2 = χ (Sz − χbz)

• No dynamics in full thermodynamic equilibrium: Sz = χb0z

• Add an additional external field b1z = γ B1z ˙ 2 =−b1z =−γ B1z

the usual precession of a magnetic moment in an applied field (Larmor precession).

• Note that this is an equilibrium state: Sz 6= χ(b0z + b1z) but is a conserved quantity— no approximations • For formal proof see Halperin and Saslow, Phys. Rev. B 16, 2154 (1977), Appendix: “the Larmor precession theorem”

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 8

Hydrodynamic Approach Hydrodynamics: a formal derivation of long wavelength dynamics of conserved quantities and broken symmetry variables in a thermodynamic approach

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 8

Hydrodynamic Approach Hydrodynamics: a formal derivation of long wavelength dynamics of conserved quantities and broken symmetry variables in a thermodynamic approach Starting points

• generalized rigidity: extra contribution to the energy density from gradients of the broken symmetry variable 1 ε = K (∇2)2 2 • thermodynamic identity

dε = Tds+ µzdsz + 8 · d(∇2) with 8 = K ∇2

• equilibrium phase dynamics (Larmor precession theorem) ˙ 2 = µz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 8

Hydrodynamic Approach Hydrodynamics: a formal derivation of long wavelength dynamics of conserved quantities and broken symmetry variables in a thermodynamic approach Starting points

• generalized rigidity: extra contribution to the energy density from gradients of the broken symmetry variable 1 ε = K (∇2)2 2 • thermodynamic identity

dε = Tds+ µzdsz + 8 · d(∇2) with 8 = K ∇2

• equilibrium phase dynamics (Larmor precession theorem) ˙ 2 = µz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 8

Hydrodynamic Approach Hydrodynamics: a formal derivation of long wavelength dynamics of conserved quantities and broken symmetry variables in a thermodynamic approach Starting points

• generalized rigidity: extra contribution to the energy density from gradients of the broken symmetry variable 1 ε = K (∇2)2 2 • thermodynamic identity

dε = Tds+ µzdsz + 8 · d(∇2) with 8 = K ∇2

• equilibrium phase dynamics (Larmor precession theorem) ˙ 2 = µz

Derive • dynamical equations for conserved quantities and broken symmetry variables for slowly varying disturbances

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 9

Rigidity and the Thermodynamic Identity In terms of the energy density

dε = Tds + µzdsz + 8 · d(∇2)

• conjugate fields are     ∂ε ∂ε µ = and 8 = z ∂s ∂∇2 z s,∇2 s,sz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 9

Rigidity and the Thermodynamic Identity In terms of the energy density

dε = Tds + µzdsz + 8 · d(∇2)

• conjugate fields are     ∂ε ∂ε µ = and 8 = z ∂s ∂∇2 z s,∇2 s,sz

Or with the free energy density f = ε − Ts

df =−sdT + µzdsz + 8 · d(∇2)

• conjugate fields are     ∂ f ∂ f µ = and 8 = z ∂s ∂∇2 z T,∇2 T,sz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 9

Rigidity and the Thermodynamic Identity In terms of the energy density

dε = Tds + µzdsz + 8 · d(∇2)

• conjugate fields are     ∂ε ∂ε µ = and 8 = z ∂s ∂∇2 z s,∇2 s,sz

Or with the free energy density f = ε − Ts

df =−sdT + µzdsz + 8 · d(∇2)

• conjugate fields are     ∂ f ∂ f µ = and 8 = z ∂s ∂∇2 z T,∇2 T,sz

These give −1 µz = χ (Sz − χbz) and 8 = K ∇2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 10

Entropy Production

Tds = dε − µzdsz − 8 · d(∇2)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 10

Entropy Production

Tds = dε − µzdsz − 8 · d(∇2)

• Form time derivative of entropy density ε µ 8 (∇2) ds = 1 d − z dsz − · d dt T dt T dt T dt

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 10

Entropy Production

Tds = dε − µzdsz − 8 · d(∇2)

• Form time derivative of entropy density ε µ 8 (∇2) ds = 1 d − z dsz − · d dt T dt T dt T dt ε • Conservation laws and dynamics of broken symmetry variable (j , jsz unknown)

ds 1 ε µ 8 =− ∇·j + z ∇·jsz − ·∇µ dt T T T z

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 10

Entropy Production

Tds = dε − µzdsz − 8 · d(∇2)

• Form time derivative of entropy density ε µ 8 (∇2) ds = 1 d − z dsz − · d dt T dt T dt T dt ε • Conservation laws and dynamics of broken symmetry variable (j , jsz unknown)

ds 1 ε µ 8 =− ∇·j + z ∇·jsz − ·∇µ dt T T T z • Entropy production equation ds =−∇·js + R with R ≥ 0 dt

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 10

Entropy Production

Tds = dε − µzdsz − 8 · d(∇2)

• Form time derivative of entropy density ε µ 8 (∇2) ds = 1 d − z dsz − · d dt T dt T dt T dt ε • Conservation laws and dynamics of broken symmetry variable (j , jsz unknown)

ds 1 ε µ 8 =− ∇·j + z ∇·jsz − ·∇µ dt T T T z • Entropy production equation ds =−∇·js + R with R ≥ 0 dt • Identify the entropy current and production

s −1 ε sz j = T (j − µzj )

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8)·∇µz

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) · ∇T − (j + 8) · ∇µz

(strategy: R should be a function of gradients of the conjugate variables)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8 =−K ∇2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8 =−K ∇2

• Energy current ε sz j = µzj

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8 =−K ∇2

• Energy current ε sz j = µzj =−µz K ∇2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8 =−K ∇2

• Energy current ε sz j = µzj =−µz K ∇2

• Entropy current js = 0

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 11

Equilibrium Dynamics Entropy Production

−1 ε sz sz RT =−T (j − µzj ) ·∇T − (j + 8) ·∇µz

(strategy: R should be a function of gradients of the conjugate variables) In the absence of dissipation the rate of entropy production must be zero. • Spin current jsz =−8 =−K ∇2

• Energy current ε sz j = µzj =−µz K ∇2

• Entropy current js = 0

We will consider adding dissipation later.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

• superfluidity occurs due to Bose condensation

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

• superfluidity occurs due to Bose condensation • the order parameter is “the expectation value of the quantum field operator for destroying a particle” 9 =hψi

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

• superfluidity occurs due to Bose condensation • the order parameter is “the expectation value of the quantum field operator for destroying a particle” 9 =hψi • 9 is a complex variable: 9 =|9|ei2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

• superfluidity occurs due to Bose condensation • the order parameter is “the expectation value of the quantum field operator for destroying a particle” 9 =hψi • 9 is a complex variable: 9 =|9|ei2 2  |9| gives the “condensate density” n0: the fraction of particles in the zero momentum state is n0/n  2 is the phase of the condensate

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 12

Superfluidity

• superfluidity occurs due to Bose condensation • the order parameter is “the expectation value of the quantum field operator for destroying a particle” 9 =hψi • 9 is a complex variable: 9 =|9|ei2 2  |9| gives the “condensate density” n0: the fraction of particles in the zero momentum state is n0/n  2 is the phase of the condensate wave function  There are a macroscopic number of particles in a single wave function and so 2 is a macroscopic thermodynamic variable,andis the broken symmetry variable.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

• Any phase gives an equivalent state; the ordered state is characterized by a particular phase

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

• Any phase gives an equivalent state; the ordered state is characterized by a particular phase • There is an energy cost for gradients of the phase Z 1 h¯ 2 E = n (∇2)2 dd x 2 s m

2  Stiffness constant K is written as ns(h¯ /m) and ns is called the superfluid density

 Stiffness constant not the same as the condensate density ns 6= n0

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

• Any phase gives an equivalent state; the ordered state is characterized by a particular phase • There is an energy cost for gradients of the phase Z 1 h¯ 2 E = n (∇2)2 dd x 2 s m

2  Stiffness constant K is written as ns(h¯ /m) and ns is called the superfluid density

 Stiffness constant not the same as the condensate density ns 6= n0 • Conjugate variable to the phase 2 is the number of particles N

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

• Any phase gives an equivalent state; the ordered state is characterized by a particular phase • There is an energy cost for gradients of the phase Z 1 h¯ 2 E = n (∇2)2 dd x 2 s m

2  Stiffness constant K is written as ns(h¯ /m) and ns is called the superfluid density

 Stiffness constant not the same as the condensate density ns 6= n0 • Conjugate variable to the phase 2 is the number of particles N • Currents and dynamics of the phase are coupled to the density, i.e., mass or electric currents

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 13

Broken Phase (Gauge) Symmetry

• Any phase gives an equivalent state; the ordered state is characterized by a particular phase • There is an energy cost for gradients of the phase Z 1 h¯ 2 E = n (∇2)2 dd x 2 s m

2  Stiffness constant K is written as ns(h¯ /m) and ns is called the superfluid density

 Stiffness constant not the same as the condensate density ns 6= n0 • Conjugate variable to the phase 2 is the number of particles N • Currents and dynamics of the phase are coupled to the density, i.e., mass or electric currents • Currents are present in equilibrium, and so are supercurrents

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 14

Supercurrents by Analogy

• One-to-one correspondence at the quantum operator level

hN¯ ≡ Sz and 2phase ≡−2spin

(e.g., ↑≡ particle, ↓≡ no particle)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 14

Supercurrents by Analogy

• One-to-one correspondence at the quantum operator level

hN¯ ≡ Sz and 2phase ≡−2spin

(e.g., ↑≡ particle, ↓≡ no particle) • Gradient of the phase gives a flow of particles ∂n =−∇·j with j = n (h¯/m)∇2 ∂t s

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 14

Supercurrents by Analogy

• One-to-one correspondence at the quantum operator level

hN¯ ≡ Sz and 2phase ≡−2spin

(e.g., ↑≡ particle, ↓≡ no particle) • Gradient of the phase gives a flow of particles ∂n =−∇·j with j = n (h¯/m)∇2 ∂t s

• Often associate a flow with a velocity: introduce superfluid velocity

vs = (h¯/m)∇2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 14

Supercurrents by Analogy

• One-to-one correspondence at the quantum operator level

hN¯ ≡ Sz and 2phase ≡−2spin

(e.g., ↑≡ particle, ↓≡ no particle) • Gradient of the phase gives a flow of particles ∂n =−∇·j with j = n (h¯/m)∇2 ∂t s

• Often associate a flow with a velocity: introduce superfluid velocity

vs = (h¯/m)∇2 and then j = nsvs

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 14

Supercurrents by Analogy

• One-to-one correspondence at the quantum operator level

hN¯ ≡ Sz and 2phase ≡−2spin

(e.g., ↑≡ particle, ↓≡ no particle) • Gradient of the phase gives a flow of particles ∂n =−∇·j with j = n (h¯/m)∇2 ∂t s

• Often associate a flow with a velocity: introduce superfluid velocity

vs = (h¯/m)∇2 and then j = nsvs • Or write in terms of flow of mass ∂ρ =−∇·g with g = ρ v ,ρ = mn ∂t s s s s

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 15

Hydrodynamic Derivation • Free energy expression: generalized rigidity and energy in external potential V h¯ 2n 1 f = s (∇2)2 + Kn2 + Vn 2m 2 (K is bulk modulus)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 15

Hydrodynamic Derivation • Free energy expression: generalized rigidity and energy in external potential V h¯ 2n 1 f = s (∇2)2 + Kn2 + Vn 2m 2 (K is bulk modulus) • Equilibrium phase dynamics (Larmor precession theorem) from dynamics with added constant potential δV : − δ 9(V, t) = 9(0, t)e ihN¯ Vt

gives   ∂ f ¯2˙ =−δ ¯2˙ =− =−µ h V or in general h ∂ n T

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 15

Hydrodynamic Derivation • Free energy expression: generalized rigidity and energy in external potential V h¯ 2n 1 f = s (∇2)2 + Kn2 + Vn 2m 2 (K is bulk modulus) • Equilibrium phase dynamics (Larmor precession theorem) from dynamics with added constant potential δV : − δ 9(V, t) = 9(0, t)e ihN¯ Vt

gives   ∂ f ¯2˙ =−δ ¯2˙ =− =−µ h V or in general h ∂ n T • Entropy production argument from the thermodynamic identity

2 dε = Tds+ µdn + 8 · d(∇2) with 8 = (h¯ ns/m)∇2 gives the current of particles

n˙ =−∇·j with j = ns (h¯/m)∇2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 16

Currents that Flow Forever

Θ

Θ(r)

L

= ρ ( / )( π/ ) g = ρs(h¯/m)(π/L) g s h¯ m 2 L H · = h vs dl m quantum of circulation

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 17

Josephson Effect

y

x 12

• Energy depends on phase difference. For weak coupling

E =−Jc cos(22 − 21)

• Change in number of particles: current I = dN2/dt dE I = d22

d.c. Josephson effect I = Jc sin(22 − 21)

• Time dependence of phase is given by the potential ˙ h¯2i =−µi ˙ ˙ a.c. Josephson effect h¯(22 − 21) =−1µ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 18

Breakdown of Superfluidity n∆µ = −s∆T + ∆P

vs

˙ h¯12 =−1µ, vs = (h¯/m)ns 12/L

• pressure or temperature difference accelerates superflow • constant superflow does not require pressure drop

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 18

Breakdown of Superfluidity n∆µ = −s∆T + ∆P

vs

˙ h¯12 =−1µ, vs = (h¯/m)ns 12/L

• pressure or temperature difference accelerates superflow • constant superflow does not require pressure drop n∆µ = −s∆T + ∆P

vs ÿ

• pressure drop (dissipation) requires passage of vortex topological defects (“quantized vortex lines”) across flow channel • presence of dissipation depends on whether vortices can be produced by thermal activation or other mechanism

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 19

Josephson Effect for a Superconductor • 2 is phase of pair wave function • expressions must be gauge invariant in presence of vector potential

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 19

Josephson Effect for a Superconductor • 2 is phase of pair wave function • expressions must be gauge invariant in presence of vector potential For bulk material   h¯ 2e Supercurrent: j = n ∇2(x)+ A s 2m hc¯ Josephson equation: h¯2˙ = 2eV

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 19

Josephson Effect for a Superconductor • 2 is phase of pair wave function • expressions must be gauge invariant in presence of vector potential For bulk material   h¯ 2e Supercurrent: j = n ∇2(x)+ A s 2m hc¯ Josephson equation: h¯2˙ = 2eV R For Josephson junction, current is I = j(y, z) dydz with  Z  2e 2 j(y, z) = jc sin 22 − 21 + Ax dx hc¯ 1 and ˙ ˙ V = (h¯/2e)(22 − 21)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 20

Josephson Junction in a B y x 12

Z ! 2e 2 Josephson current density: j(y, z) = jc sin 22 − 21 + Ax dx hc¯ 1

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 20

Josephson Junction in a Magnetic Field B y x 12

Z ! 2e 2 Josephson current density: j(y, z) = jc sin 22 − 21 + Ax dx hc¯ 1 For field Bzˆ in junction the vector potential is A =−Byxˆ, so that Z

I ∝ dy jc sin[22 − 21 − (2e/hc¯ )Byd]

giving sin(πφ/φ0) I = Ic(B) sin(22 − 21) with Ic(B) = πφ/φ0 where φ = Bld is the flux through the junction and φ0 = hc/2e is the flux quantum (2.1 × 10−7gauss cm2)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 21

Experimental Discovery of the dc Josephson Effect

Rowell, Phys. Rev. Lett. 11, 200 (1963)

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 22

SQUID Superconducting Quantum Interference Device

δΘ 1 B

δΘ2

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 22

SQUID Superconducting Quantum Interference Device

δΘ 1 B

δΘ2

  = h¯ ∇2( )+ 2e Integrate j ns 2m x hc¯ A around whole loop using fact that current j is small I 2e δ2 − δ2 = A · dl = 2πφ/φ 1 2 hc¯ 0 with φ = B × area, the flux through the loop.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 22

SQUID Superconducting Quantum Interference Device

δΘ 1 B

δΘ2

  = h¯ ∇2( )+ 2e Integrate j ns 2m x hc¯ A around whole loop using fact that current j is small I 2e δ2 − δ2 = A · dl = 2πφ/φ 1 2 hc¯ 0 with φ = B × area, the flux through the loop. Total current

I = Jc [sin δ21 + sin δ22] = (πφ/φ ) 1 (δ2 + δ2 ) 2Jc sin 0 sin[ 2 1 2 ] Maximum current varies periodically with applied field — very sensitive .

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 23

Four Sounds in a Superfluid Equations of motion for conserved quantities

ρ˙ =−∇·g g˙ =−∇P s˙ = 0

and the dynamics of the broken symmetry variable

h¯2˙ =−µ

which can be written as

ρv˙s = s∇T −∇P

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 23

Four Sounds in a Superfluid Equations of motion for conserved quantities

ρ˙ =−∇·g g˙ =−∇P s˙ = 0

and the dynamics of the broken symmetry variable

h¯2˙ =−µ

which can be written as

ρv˙s = s∇T −∇P

Need to connect the momentum density to the superfluid velocity.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn:

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density

g = ρsvs

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density (0) g = ρsvs

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density (0) g = ρsvs + ρvn

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density (0) g = ρsvs + ρvn

Define the “normal fluid density” ρn = ρ − ρs and write the transformed superfluid (0) velocity vs = vs + vn g = ρsvs + ρnvn

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density (0) g = ρsvs + ρvn

Define the “normal fluid density” ρn = ρ − ρs and write the transformed superfluid (0) velocity vs = vs + vn g = ρsvs + ρnvn

• Entropy current s j = svn

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 24

Galilean Invariance Transform to frame with a velocity −vn: • Momentum density (0) g = ρsvs + ρvn

Define the “normal fluid density” ρn = ρ − ρs and write the transformed superfluid (0) velocity vs = vs + vn g = ρsvs + ρnvn

• Entropy current s j = svn

• Momentum equation can be transformed to

ρsv˙s + ρnv˙n =−∇P

and using the equation for v˙s in the form

(ρs + ρn)v˙s = s∇T −∇P

gives

ρn(v˙s − v˙n) = s∇T

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 25

First Sound

Usual coupled density and momentum equations

ρ˙ =−∇·g g˙ =−∇P

and the pressure-density relationship (K is the bulk modulus)

δP = K δρ/ρ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 25

First Sound

Usual coupled density and momentum equations

ρ˙ =−∇·g g˙ =−∇P

and the pressure-density relationship (K is the bulk modulus)

δP = K δρ/ρ

These give first sound waves ∝ ei(q·r−ωt) propagating with the usual ω = sound speed c1q with s = K c1 ρ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 26

Second Sound

Coupled counterflow and entropy wave. Use c2  c1 ⇒ density constant, g = 0

ρsvs + ρnvn = 0 ⇒ vs − vn =−(ρ/ρs)vn

(remember ρs + ρn = ρ).

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 26

Second Sound

Coupled counterflow and entropy wave. Use c2  c1 ⇒ density constant, g = 0

ρsvs + ρnvn = 0 ⇒ vs − vn =−(ρ/ρs)vn

(remember ρs + ρn = ρ).

Entropy equation: s˙ =−s∇·vn

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 26

Second Sound

Coupled counterflow and entropy wave. Use c2  c1 ⇒ density constant, g = 0

ρsvs + ρnvn = 0 ⇒ vs − vn =−(ρ/ρs)vn

(remember ρs + ρn = ρ).

Entropy equation: s˙ =−s∇·vn

Entropy-temperature relationship (C is the specific heat): δs = CδT/T ˙ CT = sT(ρs/ρ)∇·(vs − vn)

Counterflow equation

ρn(v˙s − v˙n) = s∇T

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 26

Second Sound

Coupled counterflow and entropy wave. Use c2  c1 ⇒ density constant, g = 0

ρsvs + ρnvn = 0 ⇒ vs − vn =−(ρ/ρs)vn

(remember ρs + ρn = ρ).

Entropy equation: s˙ =−s∇·vn

Entropy-temperature relationship (C is the specific heat): δs = CδT/T ˙ CT = sT(ρs/ρ)∇·(vs − vn)

Counterflow equation

ρn(v˙s − v˙n) = s∇T These give propagating second sound waves with the speed s 2 ρs s T c2 = ρn ρC

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 27

Fourth Sound Fluid confined in porous media: no conserved momentum, no Galilean

invariance (no vn), temperature constant

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 27

Fourth Sound Fluid confined in porous media: no conserved momentum, no Galilean

invariance (no vn), temperature constant

ρ˙ =−∇·g

g = ρsvs

ρv˙s =−∇P δP =−K δρ/ρ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 27

Fourth Sound Fluid confined in porous media: no conserved momentum, no Galilean

invariance (no vn), temperature constant

ρ˙ =−∇·g

g = ρsvs

ρv˙s =−∇P δP =−K δρ/ρ

These gives a fourth sound wave propagating with the speed s ρ = s K c4 ρ ρ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 28

Third Sound

Wave propagating in thin films down to atomic layer thickness.

Like fourth sound, but involve changes of thickness rather than density, and effective compressibility depends on strength of interaction with surface.

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 29

Sounds in Helium-4 first

fourth Speed c second

Temperature T Tc s s s 2 K ρs s T ρs K c1 = , c2 = , c4 = ρ ρn ρC ρ ρ

Back Forward Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 4, April 7, 2006 30

Next Lecture

Onsager theory and the fluctuation-dissipation theorem

• Derivation and discussion • Application to nanomechanics and biodetectors

Back Forward