Spectral analysis of a completely positive map and thermal - - PowerPoint PPT Presentation

spectral analysis of a completely positive map and
SMART_READER_LITE
LIVE PREVIEW

Spectral analysis of a completely positive map and thermal - - PowerPoint PPT Presentation

Spectral analysis of a completely positive map and thermal relaxation of a QED cavity Joint work with Laurent Bruneau (Cergy) CIRMJanuary 2009 p. 1 Overview CIRMJanuary 2009 p. 2 Overview Cavity QED and the Jaynes-Cummings


slide-1
SLIDE 1

Spectral analysis of a completely positive map and thermal relaxation of a QED cavity

Joint work with Laurent Bruneau (Cergy)

CIRM–January 2009 – p. 1

slide-2
SLIDE 2

Overview

CIRM–January 2009 – p. 2

slide-3
SLIDE 3

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian

CIRM–January 2009 – p. 2

slide-4
SLIDE 4

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations

CIRM–January 2009 – p. 2

slide-5
SLIDE 5

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser

CIRM–January 2009 – p. 2

slide-6
SLIDE 6

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances

CIRM–January 2009 – p. 2

slide-7
SLIDE 7

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors

CIRM–January 2009 – p. 2

slide-8
SLIDE 8

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors
  • Ergodic properties of CP maps

CIRM–January 2009 – p. 2

slide-9
SLIDE 9

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors
  • Ergodic properties of CP maps
  • Ergodic properties of the one-atom maser

CIRM–January 2009 – p. 2

slide-10
SLIDE 10

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors
  • Ergodic properties of CP maps
  • Ergodic properties of the one-atom maser
  • Few words about the proof

CIRM–January 2009 – p. 2

slide-11
SLIDE 11

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors
  • Ergodic properties of CP maps
  • Ergodic properties of the one-atom maser
  • Few words about the proof
  • Metastable states of the one-atom maser

CIRM–January 2009 – p. 2

slide-12
SLIDE 12

Overview

  • Cavity QED and the Jaynes-Cummings Hamiltonian
  • Rabi oscillations
  • The one-atom maser
  • Rabi resonances
  • Rabi sectors
  • Ergodic properties of CP maps
  • Ergodic properties of the one-atom maser
  • Few words about the proof
  • Metastable states of the one-atom maser
  • Open questions

CIRM–January 2009 – p. 2

slide-13
SLIDE 13
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The cavity:

CIRM–January 2009 – p. 3

slide-14
SLIDE 14
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The cavity:

  • ne mode of the quantized

EM-field Hcavity = ωa∗a

CIRM–January 2009 – p. 3

slide-15
SLIDE 15
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The cavity:

  • ne mode of the quantized

EM-field Hcavity = ωa∗a The atom: a 2 level system Hatom = ω0b∗b

CIRM–January 2009 – p. 3

slide-16
SLIDE 16
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The cavity:

  • ne mode of the quantized

EM-field Hcavity = ωa∗a The atom: a 2 level system Hatom = ω0b∗b The coupling: dipole approximation Hdipol = λ 2 (b + b∗)(a + a∗)

CIRM–January 2009 – p. 3

slide-17
SLIDE 17
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The cavity:

  • ne mode of the quantized

EM-field Hcavity = ωa∗a The atom: a 2 level system Hatom = ω0b∗b The coupling: dipole approximation Hdipol = λ 2 (b + b∗)(a + a∗) The coupling: rotating wave approximation HRWA = λ 2 (b∗a + ba∗)

CIRM–January 2009 – p. 3

slide-18
SLIDE 18
  • 1. Cavity QED and the Jaynes-Cummings Hamiltonian

The Jaynes-Cummings Hamiltonian HJC = ωa∗a + ω0b∗b + λ 2 (b∗a + ba∗) The cavity:

  • ne mode of the quantized

EM-field Hcavity = ωa∗a The atom: a 2 level system Hatom = ω0b∗b The coupling: dipole approximation Hdipol = λ 2 (b + b∗)(a + a∗) The coupling: rotating wave approximation HRWA = λ 2 (b∗a + ba∗)

CIRM–January 2009 – p. 3

slide-19
SLIDE 19
  • 2. Rabi oscillations

Detuning parameter ∆ = ω − ω0 The Rotating Wave Approximation is known (at least from numerical investigations) to be accurate as long as |∆| ≪ ω0 + ω

CIRM–January 2009 – p. 4

slide-20
SLIDE 20
  • 2. Rabi oscillations

Detuning parameter ∆ = ω − ω0 The Rotating Wave Approximation is known (at least from numerical investigations) to be accurate as long as |∆| ≪ ω0 + ω

CIRM–January 2009 – p. 4

slide-21
SLIDE 21
  • 2. Rabi oscillations

Detuning parameter ∆ = ω − ω0 The Rotating Wave Approximation is known (at least from numerical investigations) to be accurate as long as |∆| ≪ ω0 + ω P(t) = ˛ ˛ ˛n − 1; ↑| e−itHJC |n; ↓ ˛ ˛ ˛

2

CIRM–January 2009 – p. 4

slide-22
SLIDE 22
  • 2. Rabi oscillations

Detuning parameter ∆ = ω − ω0 The Rotating Wave Approximation is known (at least from numerical investigations) to be accurate as long as |∆| ≪ ω0 + ω P(t) = ˛ ˛ ˛n − 1; ↑| e−itHJC |n; ↓ ˛ ˛ ˛

2

n photons Rabi frequency ΩRabi(n) = p λ2n + ∆2

CIRM–January 2009 – p. 4

slide-23
SLIDE 23
  • 2. Rabi oscillations

Detuning parameter ∆ = ω − ω0 The Rotating Wave Approximation is known (at least from numerical investigations) to be accurate as long as |∆| ≪ ω0 + ω P(t) = ˛ ˛ ˛n − 1; ↑| e−itHJC |n; ↓ ˛ ˛ ˛

2

n photons Rabi frequency ΩRabi(n) = p λ2n + ∆2 P(t) = " 1 − „ ∆ ΩRabi(n) «2# sin2 ΩRabi(n)t/2

CIRM–January 2009 – p. 4

slide-24
SLIDE 24
  • 3. The one-atom maser

CIRM–January 2009 – p. 5

slide-25
SLIDE 25
  • 3. The one-atom maser

Repeated interaction scheme H = Hcavity ⊗ Hbeam, Hbeam = O

n≥1

Hatom n

CIRM–January 2009 – p. 5

slide-26
SLIDE 26
  • 3. The one-atom maser

Repeated interaction scheme H = Hcavity ⊗ Hbeam, Hbeam = O

n≥1

Hatom n Hn = HJC acting on Hcavity ⊗ Hatom n

CIRM–January 2009 – p. 5

slide-27
SLIDE 27
  • 3. The one-atom maser

Repeated interaction scheme H = Hcavity ⊗ Hbeam, Hbeam = O

n≥1

Hatom n Hn = HJC acting on Hcavity ⊗ Hatom n Cavity state after n interactions ρn = TrHbeam " e−iτHn · · · e−iτH1 ρ0 ⊗

n

O

k=1

ρatom k ! eiτH1 · · · eiτHn #

CIRM–January 2009 – p. 5

slide-28
SLIDE 28
  • 3. The one-atom maser

Thermal beam: ρatom k = ρβ = Z−1e−βHatom

CIRM–January 2009 – p. 6

slide-29
SLIDE 29
  • 3. The one-atom maser

Thermal beam: ρatom k = ρβ = Z−1e−βHatom ρn = Lβ(ρn−1)

CIRM–January 2009 – p. 6

slide-30
SLIDE 30
  • 3. The one-atom maser

Thermal beam: ρatom k = ρβ = Z−1e−βHatom ρn = Lβ(ρn−1) Reduced dynamics Lβ(ρ) = TrHatom h e−iτHJC “ ρ ⊗ ρβ” eiτHJC i Completely positive, trace preserving map on the trace ideal J 1(Hcavity)

CIRM–January 2009 – p. 6

slide-31
SLIDE 31
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ

CIRM–January 2009 – p. 7

slide-32
SLIDE 32
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ

|n

CIRM–January 2009 – p. 7

slide-33
SLIDE 33
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ

|n

CIRM–January 2009 – p. 7

slide-34
SLIDE 34
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ Dimensionless parameters η = „ ∆τ 2π «2 , ξ = „ λτ 2π «2

CIRM–January 2009 – p. 7

slide-35
SLIDE 35
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ Dimensionless parameters η = „ ∆τ 2π «2 , ξ = „ λτ 2π «2 R(η, ξ) = {n ∈ N∗|ξn + η is a perfect square}

CIRM–January 2009 – p. 7

slide-36
SLIDE 36
  • 4. Rabi resonances

A resonance occurs when the interaction time τ is a multiple of the Rabi period ΩRabi(n)τ ∈ 2πZ Dimensionless parameters η = „ ∆τ 2π «2 , ξ = „ λτ 2π «2 R(η, ξ) = {n ∈ N∗|ξn + η is a perfect square}

  • Definition. The system is
  • Non resonant: R(η, ξ) is empty.
  • Simply resonant: R(η, ξ) = {n1}.
  • Fully resonant: R(η, ξ) = {n1, n2, . . .} i.e. has ∞-many resonances.
  • Degenerate: fully resonant and there exist n ∈ R(η, ξ) ∪ {0} and m ∈ R(η, ξ) such

that n + 1, m + 1 ∈ R(η, ξ).

CIRM–January 2009 – p. 7

slide-37
SLIDE 37
  • 4. Rabi resonances
  • Lemma. If η and ξ are

CIRM–January 2009 – p. 8

slide-38
SLIDE 38
  • 4. Rabi resonances
  • Lemma. If η and ξ are
  • not both irrational or not both rational: non-resonant;

CIRM–January 2009 – p. 8

slide-39
SLIDE 39
  • 4. Rabi resonances
  • Lemma. If η and ξ are
  • not both irrational or not both rational: non-resonant;
  • both irrational: either simply resonant or non-resonant (the generic case);

CIRM–January 2009 – p. 8

slide-40
SLIDE 40
  • 4. Rabi resonances
  • Lemma. If η and ξ are
  • not both irrational or not both rational: non-resonant;
  • both irrational: either simply resonant or non-resonant (the generic case);
  • both rational: η = a/b, ξ = c/d (irreducible) and m = LCM(b, d)

X = {x ∈ {0, . . . , ξm − 1} | x2m ≡ ηm (mod ξm)} then non-resonant if X is empty or fully resonant R(η, ξ) = {(k2 − η)/ξ | k = jmξ + x, j ∈ N∗, x ∈ X} ∩ N∗

CIRM–January 2009 – p. 8

slide-41
SLIDE 41
  • 4. Rabi resonances
  • Lemma. If η and ξ are
  • not both irrational or not both rational: non-resonant;
  • both irrational: either simply resonant or non-resonant (the generic case);
  • both rational: η = a/b, ξ = c/d (irreducible) and m = LCM(b, d)

X = {x ∈ {0, . . . , ξm − 1} | x2m ≡ ηm (mod ξm)} then non-resonant if X is empty or fully resonant R(η, ξ) = {(k2 − η)/ξ | k = jmξ + x, j ∈ N∗, x ∈ X} ∩ N∗ Moreover, if degenerate then η and ξ are integers such that η > 0 is a quadratic residue modulo ξ.

CIRM–January 2009 – p. 8

slide-42
SLIDE 42
  • 4. Rabi resonances
  • Lemma. If η and ξ are
  • not both irrational or not both rational: non-resonant;
  • both irrational: either simply resonant or non-resonant (the generic case);
  • both rational: η = a/b, ξ = c/d (irreducible) and m = LCM(b, d)

X = {x ∈ {0, . . . , ξm − 1} | x2m ≡ ηm (mod ξm)} then non-resonant if X is empty or fully resonant R(η, ξ) = {(k2 − η)/ξ | k = jmξ + x, j ∈ N∗, x ∈ X} ∩ N∗ Moreover, if degenerate then η and ξ are integers such that η > 0 is a quadratic residue modulo ξ.

  • Remark. This lemma is elementary but characterizing integers η, ξ for which the system is

degenerate is a very hard (open) problem in Diophantine analysis.

CIRM–January 2009 – p. 8

slide-43
SLIDE 43
  • 5. Rabi sectors

Decomposition into Rabi sectors ℓ2(N) = Hcavity =

r

M

k=1

H(k) where H(k) = ℓ2(Ik)

CIRM–January 2009 – p. 9

slide-44
SLIDE 44
  • 5. Rabi sectors

Decomposition into Rabi sectors ℓ2(N) = Hcavity =

r

M

k=1

H(k) where H(k) = ℓ2(Ik) and r = 1 I1 ≡ N if R(η, ξ) is empty, r = 2 I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, n1 + 1, . . .} if R(η, ξ) = {n1}, r = ∞ I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, . . . , n2 − 1}, . . . if R(η, ξ) = {n1, n2, . . .}.

CIRM–January 2009 – p. 9

slide-45
SLIDE 45
  • 5. Rabi sectors

Decomposition into Rabi sectors ℓ2(N) = Hcavity =

r

M

k=1

H(k) where H(k) = ℓ2(Ik) and r = 1 I1 ≡ N if R(η, ξ) is empty, r = 2 I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, n1 + 1, . . .} if R(η, ξ) = {n1}, r = ∞ I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, . . . , n2 − 1}, . . . if R(η, ξ) = {n1, n2, . . .}. Pk denotes the orthogonal projection onto H(k)

CIRM–January 2009 – p. 9

slide-46
SLIDE 46
  • 5. Rabi sectors

Decomposition into Rabi sectors ℓ2(N) = Hcavity =

r

M

k=1

H(k) where H(k) = ℓ2(Ik) and r = 1 I1 ≡ N if R(η, ξ) is empty, r = 2 I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, n1 + 1, . . .} if R(η, ξ) = {n1}, r = ∞ I1 ≡ {0, . . . , n1 − 1}, I2 ≡ {n1, . . . , n2 − 1}, . . . if R(η, ξ) = {n1, n2, . . .}. Pk denotes the orthogonal projection onto H(k) Partial Gibbs state in H(k): ρ(k)β∗

cavity =

e−β∗HcavityPk Tr e−β∗HcavityPk , β∗ = β ω0 ω

CIRM–January 2009 – p. 9

slide-47
SLIDE 47
  • 6. Ergodic properties of CP maps

Support of a density matrix ρ is the orthogonal projection s(ρ) onto the closure of Ran ρ.

CIRM–January 2009 – p. 10

slide-48
SLIDE 48
  • 6. Ergodic properties of CP maps

Support of a density matrix ρ is the orthogonal projection s(ρ) onto the closure of Ran ρ. µ ≪ ρ means s(µ) ≤ s(ρ)

CIRM–January 2009 – p. 10

slide-49
SLIDE 49
  • 6. Ergodic properties of CP maps

Support of a density matrix ρ is the orthogonal projection s(ρ) onto the closure of Ran ρ. µ ≪ ρ means s(µ) ≤ s(ρ) ρ is ergodic for the CP map L iff, for all µ ≪ ρ, A ∈ B(Hcavity) lim

N→∞

1 N

N

X

n=1

(Ln(µ)) (A) = ρ(A)

CIRM–January 2009 – p. 10

slide-50
SLIDE 50
  • 6. Ergodic properties of CP maps

Support of a density matrix ρ is the orthogonal projection s(ρ) onto the closure of Ran ρ. µ ≪ ρ means s(µ) ≤ s(ρ) ρ is ergodic for the CP map L iff, for all µ ≪ ρ, A ∈ B(Hcavity) lim

N→∞

1 N

N

X

n=1

(Ln(µ)) (A) = ρ(A) ρ is mixing for L iff, for all µ ≪ ρ, A ∈ B(Hcavity) lim

n→∞ (Ln(µ)) (A) = ρ(A),

CIRM–January 2009 – p. 10

slide-51
SLIDE 51
  • 6. Ergodic properties of CP maps

Support of a density matrix ρ is the orthogonal projection s(ρ) onto the closure of Ran ρ. µ ≪ ρ means s(µ) ≤ s(ρ) ρ is ergodic for the CP map L iff, for all µ ≪ ρ, A ∈ B(Hcavity) lim

N→∞

1 N

N

X

n=1

(Ln(µ)) (A) = ρ(A) ρ is mixing for L iff, for all µ ≪ ρ, A ∈ B(Hcavity) lim

n→∞ (Ln(µ)) (A) = ρ(A),

and exponentially mixing iff |(Ln(µ)) (A) − ρ(A)| ≤ CA,µ e−αn, for some constants CA,µ and α > 0.

CIRM–January 2009 – p. 10

slide-52
SLIDE 52
  • 7. Ergodic properties of the one-atom maser

Main Theorem. 1. If the system is non-resonant then Lβ has no invariant state for β ≤ 0 and a unique ergodic state ρβ∗

cavity =

e−β∗Hcavity Tr e−β∗Hcavity , β∗ = β ω0 ω for β > 0. In the latter case any initial state relaxes in the mean to this thermal equilibrium state lim

N→∞

1 N

N

X

n=1

“ Ln

β(µ)

” (A) = ρβ∗

cavity(A)

for any A ∈ B(Hcavity).

CIRM–January 2009 – p. 11

slide-53
SLIDE 53
  • 7. Ergodic properties of the one-atom maser

Main Theorem. 2. If the system is simply resonant then Lβ has the unique ergodic state ρ(1) β∗

cavity if β ≤ 0 and two ergodic states ρ(1) β∗ cavity, ρ(2) β∗ cavity if β > 0. In the latter case, for any

state µ, one has lim

N→∞

1 N

N

X

n=1

“ Ln

β(µ)

” (A) = µ(P1) ρ(1) β∗

cavity(A) + µ(P2) ρ(2) β∗ cavity(A),

for any A ∈ B(Hcavity).

CIRM–January 2009 – p. 11

slide-54
SLIDE 54
  • 7. Ergodic properties of the one-atom maser

Main Theorem. 3. If the system is fully resonant then for any β ∈ R, Lβ has infinitely many ergodic states ρ(k) β∗

cavity, k = 1, 2, . . . Moreover, if the system is non-degenerate,

lim

N→∞

1 N

N

X

n=1

“ Ln

β(µ)

” (A) =

X

k=1

µ(Pk) ρ(k) β∗

cavity(A),

holds for any state µ and all A ∈ B(Hcavity).

CIRM–January 2009 – p. 11

slide-55
SLIDE 55
  • 7. Ergodic properties of the one-atom maser

Main Theorem. 4. If the system is fully resonant and degenerate there exists a finite set D(η, ξ) ⊂ Z such that the conclusions of 3. still hold provided the non-resonance condition (NR) ei(τω+ξπ)d = 1 is satisfied for all d ∈ D(η, ξ).

  • 5. In all the previous cases any invariant state is diagonal and can be represented as a

convex linear combination of ergodic states, i.e., the set of invariant states is a simplex whose extremal points are ergodic states. In the remaining case, i.e., if condition (NR) fails, there are non-diagonal invariant states.

  • 6. Whenever the state ρ(k) β∗

cavity is ergodic it is also exponentially mixing if the Rabi sector

H(k)

cavity is finite dimensional.

CIRM–January 2009 – p. 11

slide-56
SLIDE 56
  • 7. Ergodic properties of the one-atom maser

Remarks.

CIRM–January 2009 – p. 11

slide-57
SLIDE 57
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.

CIRM–January 2009 – p. 11

slide-58
SLIDE 58
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.
  • Mixing is very slow in infinite dimensional Rabi sectors due to the presence of

∞-many metastable states (1 ∈ spess(Lβ)).

CIRM–January 2009 – p. 11

slide-59
SLIDE 59
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.
  • Mixing is very slow in infinite dimensional Rabi sectors due to the presence of

∞-many metastable states (1 ∈ spess(Lβ)).

  • For given η, ξ it is easy to compute the set D(η, ξ). However it is extremely hard

(and an open problem) to characterize those integers η and ξ for which D(η, ξ) is non-empty.

CIRM–January 2009 – p. 11

slide-60
SLIDE 60
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.
  • Mixing is very slow in infinite dimensional Rabi sectors due to the presence of

∞-many metastable states (1 ∈ spess(Lβ)).

  • For given η, ξ it is easy to compute the set D(η, ξ). However it is extremely hard

(and an open problem) to characterize those integers η and ξ for which D(η, ξ) is non-empty.

  • Degenerate fully resonant systems exist. If η = 241 and ξ = 720 then

720 + 241 = 292, 2 · 720 + 241 = 412, 3 · 720 + 241 = 492 so that 1 and 2 are successive Rabi resonances. In this case D(241, 720) = {1}.

CIRM–January 2009 – p. 11

slide-61
SLIDE 61
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.
  • Mixing is very slow in infinite dimensional Rabi sectors due to the presence of

∞-many metastable states (1 ∈ spess(Lβ)).

  • For given η, ξ it is easy to compute the set D(η, ξ). However it is extremely hard

(and an open problem) to characterize those integers η and ξ for which D(η, ξ) is non-empty.

  • Degenerate fully resonant systems exist. If η = 241 and ξ = 720 then

720 + 241 = 292, 2 · 720 + 241 = 412, 3 · 720 + 241 = 492 so that 1 and 2 are successive Rabi resonances. In this case D(241, 720) = {1}.

  • Another example is η = 1 and ξ = 840 for which 1, 2, 52 and 53 are Rabi

resonances 840 + 1 = 292, 2 · 840 + 1 = 412, 52 · 840 + 1 = 2092, 53 · 840 + 1 = 2112 and D(1, 840) = {51}.

CIRM–January 2009 – p. 11

slide-62
SLIDE 62
  • 7. Ergodic properties of the one-atom maser

Remarks.

  • Numerical experiments support the conjecture that all ergodic states are mixing.
  • Mixing is very slow in infinite dimensional Rabi sectors due to the presence of

∞-many metastable states (1 ∈ spess(Lβ)).

  • For given η, ξ it is easy to compute the set D(η, ξ). However it is extremely hard

(and an open problem) to characterize those integers η and ξ for which D(η, ξ) is non-empty.

  • Degenerate fully resonant systems exist. If η = 241 and ξ = 720 then

720 + 241 = 292, 2 · 720 + 241 = 412, 3 · 720 + 241 = 492 so that 1 and 2 are successive Rabi resonances. In this case D(241, 720) = {1}.

  • Another example is η = 1 and ξ = 840 for which 1, 2, 52 and 53 are Rabi

resonances 840 + 1 = 292, 2 · 840 + 1 = 412, 52 · 840 + 1 = 2092, 53 · 840 + 1 = 2112 and D(1, 840) = {51}.

  • We do not know of any example where D(η, ξ) contains more than one element.

CIRM–January 2009 – p. 11

slide-63
SLIDE 63
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ.

CIRM–January 2009 – p. 12

slide-64
SLIDE 64
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ. Main difficulty: perturbation theory does’nt work! At zero coupling (λ = 0) one has Lβ(ρ) = e−iτHcavityρ eiτHcavity so that sp(Lβ) = sppp(Lβ) = {eiτωd | d ∈ Z} is either finite (if ωτ ∈ 2πQ) or dense on the unit circle, but always infinitely degenerate.

CIRM–January 2009 – p. 12

slide-65
SLIDE 65
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ. Main difficulty: perturbation theory does’nt work! At zero coupling (λ = 0) one has Lβ(ρ) = e−iτHcavityρ eiτHcavity so that sp(Lβ) = sppp(Lβ) = {eiτωd | d ∈ Z} is either finite (if ωτ ∈ 2πQ) or dense on the unit circle, but always infinitely degenerate. Main tools:

CIRM–January 2009 – p. 12

slide-66
SLIDE 66
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ. Main difficulty: perturbation theory does’nt work! At zero coupling (λ = 0) one has Lβ(ρ) = e−iτHcavityρ eiτHcavity so that sp(Lβ) = sppp(Lβ) = {eiτωd | d ∈ Z} is either finite (if ωτ ∈ 2πQ) or dense on the unit circle, but always infinitely degenerate. Main tools:

  • Use gauge symmetry! It follows from [HJC, a∗a + b∗b] = [Hatom, ρβ

atom] = 0 that

Lβ(e−iθa∗aXeiθa∗a) = e−iθa∗aLβ(X)eiθa∗a

CIRM–January 2009 – p. 12

slide-67
SLIDE 67
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ. Main difficulty: perturbation theory does’nt work! At zero coupling (λ = 0) one has Lβ(ρ) = e−iτHcavityρ eiτHcavity so that sp(Lβ) = sppp(Lβ) = {eiτωd | d ∈ Z} is either finite (if ωτ ∈ 2πQ) or dense on the unit circle, but always infinitely degenerate. Main tools:

  • Use gauge symmetry! It follows from [HJC, a∗a + b∗b] = [Hatom, ρβ

atom] = 0 that

Lβ(e−iθa∗aXeiθa∗a) = e−iθa∗aLβ(X)eiθa∗a

  • Use the block structure induced by Rabi sectors.

CIRM–January 2009 – p. 12

slide-68
SLIDE 68
  • 8. Few words about the proof

Main idea: control the peripheral spectrum of Lβ. Main difficulty: perturbation theory does’nt work! At zero coupling (λ = 0) one has Lβ(ρ) = e−iτHcavityρ eiτHcavity so that sp(Lβ) = sppp(Lβ) = {eiτωd | d ∈ Z} is either finite (if ωτ ∈ 2πQ) or dense on the unit circle, but always infinitely degenerate. Main tools:

  • Use gauge symmetry! It follows from [HJC, a∗a + b∗b] = [Hatom, ρβ

atom] = 0 that

Lβ(e−iθa∗aXeiθa∗a) = e−iθa∗aLβ(X)eiθa∗a

  • Use the block structure induced by Rabi sectors.
  • Use Schrader’s version of Perron-Frobenius theory for trace preserving CP maps
  • n trace ideals [Fields Inst. Commun. 30 (2001)].

CIRM–January 2009 – p. 12

slide-69
SLIDE 69
  • 9. Metastable states of the one-atom maser

By gauge symmetry, the subspace of diagonal states is invariant. The action of Lβ on this subspace is conjugated to that of L = I − ∇∗D(N)e−βω0N∇eβω0N

  • n ℓ1(N) where

(Nx)n = nxn, (∇x)n = ( x0 for n = 0; xn − xn−1 for n ≥ 1; (∇∗x)n = xn − xn+1 and D(N) = 1 1 + e−βω0 ξN ξN + η sin2(π p ξN + η)

CIRM–January 2009 – p. 13

slide-70
SLIDE 70
  • 9. Metastable states of the one-atom maser

By gauge symmetry, the subspace of diagonal states is invariant. The action of Lβ on this subspace is conjugated to that of L = I − ∇∗D(N)e−βω0N∇eβω0N

  • n ℓ1(N) where

(Nx)n = nxn, (∇x)n = ( x0 for n = 0; xn − xn−1 for n ≥ 1; (∇∗x)n = xn − xn+1 and D(N) = 1 1 + e−βω0 ξN ξN + η sin2(π p ξN + η) Rabi resonances are integers n such that D(n) = 0. They decouple L.

CIRM–January 2009 – p. 13

slide-71
SLIDE 71
  • 9. Metastable states of the one-atom maser

There is an increasing sequence mk such that D(mk) = O(k−2). They almost decouple L: Rabi quasi-resonances.

CIRM–January 2009 – p. 14

slide-72
SLIDE 72
  • 9. Metastable states of the one-atom maser

There is an increasing sequence mk such that D(mk) = O(k−2). They almost decouple L: Rabi quasi-resonances. Setting L0 = I − ∇∗D0(N)e−βω0N∇eβω0N with D0(n) = ( if n ∈ {m1, m2, . . .} D(n)

  • therwise,

we get L = L0 + trace class.

CIRM–January 2009 – p. 14

slide-73
SLIDE 73
  • 9. Metastable states of the one-atom maser

There is an increasing sequence mk such that D(mk) = O(k−2). They almost decouple L: Rabi quasi-resonances. Setting L0 = I − ∇∗D0(N)e−βω0N∇eβω0N with D0(n) = ( if n ∈ {m1, m2, . . .} D(n)

  • therwise,

we get L = L0 + trace class. L0 has infinitely degenerate eigenvalue 1: eigenvectors are metastable states of L

CIRM–January 2009 – p. 14

slide-74
SLIDE 74
  • 9. Metastable states of the one-atom maser

The metastable cascade

10 10

1

10

2

10

3

10

4

10

5

10

−1

10 10

1

10

2

CIRM–January 2009 – p. 15

slide-75
SLIDE 75
  • 9. Metastable states of the one-atom maser

Local equilibrium after 5000 interactions

3 5 8 11 15 20 25 31 38 45 53 10

−4

10

−3

10

−2

10

−1

10

CIRM–January 2009 – p. 15

slide-76
SLIDE 76
  • 9. Metastable states of the one-atom maser

Local equilibrium after 50000 interactions

3 5 8 11 15 20 25 31 38 45 53 10

−4

10

−3

10

−2

10

−1

10

CIRM–January 2009 – p. 15

slide-77
SLIDE 77
  • 9. Metastable states of the one-atom maser

Mean photon number

1 2 3 4 x 10

5

7 9 11 13 15 17 19 21

CIRM–January 2009 – p. 15

slide-78
SLIDE 78
  • 10. Open questions

CIRM–January 2009 – p. 16

slide-79
SLIDE 79
  • 10. Open questions
  • Beyond the rotating wave approximation.

CIRM–January 2009 – p. 16

slide-80
SLIDE 80
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?

CIRM–January 2009 – p. 16

slide-81
SLIDE 81
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?
  • Mixing in infinite Rabi sectors.

CIRM–January 2009 – p. 16

slide-82
SLIDE 82
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?
  • Mixing in infinite Rabi sectors.
  • Leaky cavities: coupling the full EM field.

CIRM–January 2009 – p. 16

slide-83
SLIDE 83
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?
  • Mixing in infinite Rabi sectors.
  • Leaky cavities: coupling the full EM field.
  • For number theorists and algebraic geometers: The Diophantine problem.

CIRM–January 2009 – p. 16

slide-84
SLIDE 84
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?
  • Mixing in infinite Rabi sectors.
  • Leaky cavities: coupling the full EM field.
  • For number theorists and algebraic geometers: The Diophantine problem.
  • For probabilitsts: Random interaction times.

CIRM–January 2009 – p. 16

slide-85
SLIDE 85
  • 10. Open questions
  • Beyond the rotating wave approximation.
  • What is the mathematical status of RWA ?
  • Mixing in infinite Rabi sectors.
  • Leaky cavities: coupling the full EM field.
  • For number theorists and algebraic geometers: The Diophantine problem.
  • For probabilitsts: Random interaction times.
  • ...

CIRM–January 2009 – p. 16