Collective Dynamics of a Generalized Dicke Model

Transcription

1 Collective Dynamics of a Generalized Dicke Model J. Keeling, J. A. Mayoh, M. J. Bhaseen, B. D. Simons Harvard, January 212 Funding: Jonathan Keeling Collective dynamics Harvard, January / 25

2 Coupling many atoms to light Old question: What happens to radiation when many atoms interact collectively with light. Superradiance dynamical and steady state. Jonathan Keeling Collective dynamics Harvard, January / 25

3 Coupling many atoms to light Old question: What happens to radiation when many atoms interact collectively with light. Superradiance dynamical and steady state. New relevance Superconducting qubits Quantum dots Nitrogen-vacancies in diamond Ultra-cold atoms κ κ Cavity Pump Rydberg atoms Pump Jonathan Keeling Collective dynamics Harvard, January / 25

4 Dicke effect: Enhanced emission H int = k,i ) g k (ψ k S i e ik r i + H.c. Jonathan Keeling Collective dynamics Harvard, January / 25

5 Dicke effect: Enhanced emission H int = k,i ) g k (ψ k S i e ik r i + H.c. If r i r j λ, use i S i S Collective decay: dρ dt = Γ 2 [ S + S ρ S ρs + + ρs + S ] Jonathan Keeling Collective dynamics Harvard, January / 25

6 Dicke effect: Enhanced emission H int = k,i ) g k (ψ k S i e ik r i + H.c. If r i r j λ, use i S i S Collective decay: dρ dt = Γ 2 [ S + S ρ S ρs + + ρs + S ] If S z = S = N/2 initially: I Γ d Sz dt = ΓN2 4 sech2 [ ] ΓN 2 t S z N/2 -N/2 t D t D ΓN 2 /2 I=-Γd S z /dt Jonathan Keeling Collective dynamics Harvard, January / 25

7 Dicke effect: Enhanced emission H int = k,i ) g k (ψ k S i e ik r i + H.c. If r i r j λ, use i S i S Collective decay: dρ dt = Γ 2 [ S + S ρ S ρs + + ρs + S ] If S z = S = N/2 initially: I Γ d Sz dt = ΓN2 4 sech2 [ ] ΓN 2 t S z N/2 -N/2 t D t D ΓN 2 /2 I=-Γd S z /dt Problem: dipole interactions dephase. [Friedberg et al, Phys. Lett. 1972] Jonathan Keeling Collective dynamics Harvard, January / 25

8 Collective radiation with a cavity: Dynamics H int = i ( ψ S i + ψs + i ) Single cavity mode: oscillations [Bonifacio and Preparata PRA 7; JK PRA 9] Jonathan Keeling Collective dynamics Harvard, January / 25

9 Collective radiation with a cavity: Dynamics H int = i ( ψ S i + ψs + i ) Single cavity mode: oscillations If S z = S = N/2 initially: ψ(t) 2 16 T=2ln( N )/ N / N Time [Bonifacio and Preparata PRA 7; JK PRA 9] Jonathan Keeling Collective dynamics Harvard, January / 25

10 Collective radiation with a cavity: Dynamics H int = i ( ψ S i + ψs + i ) Single cavity mode: oscillations If S z = S = N/2 initially: ψ(t) 2 16 T=2ln( N )/ N / N Time [Bonifacio and Preparata PRA 7; JK PRA 9] ψ(t) 2 16 N= Time Jonathan Keeling Collective dynamics Harvard, January / 25

11 Dicke model: Equilibrium superradiance transition ( H = ωψ ψ + ω S z + g ψ S + ψs +). Coherent state: Ψ e λψ +ηs + Ω Energy: E = ω λ 2 + ω N η η gn η λ + λ η 1 + η 2 Small g, min at λ, η = [Hepp, Lieb, Ann. Phys. 73] Jonathan Keeling Collective dynamics Harvard, January / 25

12 Dicke model: Equilibrium superradiance transition ( H = ωψ ψ + ω S z + g ψ S + ψs +). Coherent state: Ψ e λψ +ηs + Ω Energy: E = ω λ 2 + ω N η η gn η λ + λ η 1 + η 2 Small g, min at λ, η = [Hepp, Lieb, Ann. Phys. 73] Jonathan Keeling Collective dynamics Harvard, January / 25

13 Dicke model: Equilibrium superradiance transition ( H = ωψ ψ + ω S z + g ψ S + ψs +). Coherent state: Ψ e λψ +ηs + Ω Energy: E = ω λ 2 + ω N η η gn η λ + λ η 1 + η 2 Small g, min at λ, η = [Hepp, Lieb, Ann. Phys. 73] Jonathan Keeling Collective dynamics Harvard, January / 25

14 Dicke model: Equilibrium superradiance transition ( H = ωψ ψ + ω S z + g ψ S + ψs +). Coherent state: Ψ e λψ +ηs + Ω Energy: E = ω λ 2 + ω N η η gn η λ + λ η 1 + η 2 Small g, min at λ, η = ω ḡ N SR Spontaneous polarisation if: Ng 2 > ωω [Hepp, Lieb, Ann. Phys. 73] Jonathan Keeling Collective dynamics Harvard, January / 25

15 No go theorem and transition Spontaneous polarisation if: Ng 2 > ωω [Rzazewski et al PRL 75] Jonathan Keeling Collective dynamics Harvard, January / 25

16 No go theorem and transition Spontaneous polarisation if: Ng 2 > ωω No go theorem:. Minimal coupling (p ea) 2 /2m i e m A p i g(ψ S + ψs + ), i A 2 2m Nζ(ψ + ψ ) 2 [Rzazewski et al PRL 75] Jonathan Keeling Collective dynamics Harvard, January / 25

17 No go theorem and transition Spontaneous polarisation if: Ng 2 > ωω No go theorem:. Minimal coupling (p ea) 2 /2m i e m A p i g(ψ S + ψs + ), For large N, ω ω + 2Nζ. (RWA) i A 2 2m Nζ(ψ + ψ ) 2 [Rzazewski et al PRL 75] Jonathan Keeling Collective dynamics Harvard, January / 25

18 No go theorem and transition Spontaneous polarisation if: Ng 2 > ωω No go theorem:. Minimal coupling (p ea) 2 /2m i e m A p i g(ψ S + ψs + ), For large N, ω ω + 2Nζ. (RWA) Need Ng 2 > ω (ω + 2Nζ). i A 2 2m Nζ(ψ + ψ ) 2 [Rzazewski et al PRL 75] Jonathan Keeling Collective dynamics Harvard, January / 25

19 No go theorem and transition Spontaneous polarisation if: Ng 2 > ωω No go theorem:. Minimal coupling (p ea) 2 /2m i e m A p i g(ψ S + ψs + ), i A 2 2m Nζ(ψ + ψ ) 2 For large N, ω ω + 2Nζ. (RWA) Need Ng 2 > ω (ω + 2Nζ). But Thomas-Reiche-Kuhn sum rule states: g 2 /ω < 2ζ. No transition [Rzazewski et al PRL 75] Jonathan Keeling Collective dynamics Harvard, January / 25

20 Dicke phase transition: ways out Problem: g 2 /ω < 2ζ for intrinsic parameters. Solutions: Non-solution Ferroelectric transition in D r gauge. [JK JPCM 7 ] See also [Nataf and Ciuti, Nat. Comm. 1; Viehmann et al. PRL 11] Grand canonical ensemble: If H H µ(s z + ψ ψ), need only: g 2 N > (ω µ)(ω µ) Incoherent pumping polariton condensation. Dissociate g, ω, e.g. Raman scheme: ω ω. [Dimer et al. PRA 7; Baumann et al. Nature 1. Also, Black et al. PRL 3 ] Jonathan Keeling Collective dynamics Harvard, January / 25

21 Dicke phase transition: ways out Problem: g 2 /ω < 2ζ for intrinsic parameters. Solutions: Non-solution Ferroelectric transition in D r gauge. [JK JPCM 7 ] See also [Nataf and Ciuti, Nat. Comm. 1; Viehmann et al. PRL 11] Grand canonical ensemble: If H H µ(s z + ψ ψ), need only: g 2 N > (ω µ)(ω µ) Incoherent pumping polariton condensation. Dissociate g, ω, e.g. Raman scheme: ω ω. [Dimer et al. PRA 7; Baumann et al. Nature 1. Also, Black et al. PRL 3 ] Jonathan Keeling Collective dynamics Harvard, January / 25

22 Dicke phase transition: ways out Problem: g 2 /ω < 2ζ for intrinsic parameters. Solutions: Non-solution Ferroelectric transition in D r gauge. [JK JPCM 7 ] See also [Nataf and Ciuti, Nat. Comm. 1; Viehmann et al. PRL 11] Grand canonical ensemble: If H H µ(s z + ψ ψ), need only: g 2 N > (ω µ)(ω µ) Incoherent pumping polariton condensation. Dissociate g, ω, e.g. Raman scheme: ω ω. [Dimer et al. PRA 7; Baumann et al. Nature 1. Also, Black et al. PRL 3 ] Jonathan Keeling Collective dynamics Harvard, January / 25

23 Dicke phase transition: ways out Problem: g 2 /ω < 2ζ for intrinsic parameters. Solutions: Non-solution Ferroelectric transition in D r gauge. [JK JPCM 7 ] See also [Nataf and Ciuti, Nat. Comm. 1; Viehmann et al. PRL 11] Grand canonical ensemble: If H H µ(s z + ψ ψ), need only: g 2 N > (ω µ)(ω µ) Incoherent pumping polariton condensation. Dissociate g, ω, e.g. Raman scheme: ω ω. [Dimer et al. PRA 7; Baumann et al. Nature 1. Also, Black et al. PRL 3 ] Jonathan Keeling Collective dynamics Harvard, January / 25

24 Dicke phase transition: ways out Problem: g 2 /ω < 2ζ for intrinsic parameters. Solutions: Non-solution Ferroelectric transition in D r gauge. [JK JPCM 7 ] See also [Nataf and Ciuti, Nat. Comm. 1; Viehmann et al. PRL 11] Grand canonical ensemble: If H H µ(s z + ψ ψ), need only: g 2 N > (ω µ)(ω µ) Incoherent pumping polariton condensation. Dissociate g, ω, e.g. Raman scheme: ω ω. [Dimer et al. PRA 7; Baumann et al. Nature 1. Also, Black et al. PRL 3 ] κ Pump κ Cavity Pump Jonathan Keeling Collective dynamics Harvard, January / 25

25 Outline 1 Introduction: Dicke model and superradiance Rayleigh scheme: Generalised Dicke model 2 Attractors of dynamics (fixed points) Phase diagram: comparison to experiment 3 Approach to attractors: timescales Why slow timescales emerge 4 Attractors of dynamics (oscillations) Reaching other parameter ranges Jonathan Keeling Collective dynamics Harvard, January / 25

26 Outline 1 Introduction: Dicke model and superradiance Rayleigh scheme: Generalised Dicke model 2 Attractors of dynamics (fixed points) Phase diagram: comparison to experiment 3 Approach to attractors: timescales Why slow timescales emerge 4 Attractors of dynamics (oscillations) Reaching other parameter ranges Jonathan Keeling Collective dynamics Harvard, January / 25

27 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Jonathan Keeling Collective dynamics Harvard, January / 25

28 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Jonathan Keeling Collective dynamics Harvard, January / 25

29 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± Feedback: U g2 ω c ω a H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Jonathan Keeling Collective dynamics Harvard, January / 25

30 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± Feedback: U g2 ω c ω a H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Jonathan Keeling Collective dynamics Harvard, January / 25

31 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± Feedback: U g2 ω c ω a H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Semiclassical EOM ( S = N/2 1) Ṡ = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω+us z )] ψ ig(s + S + ) Jonathan Keeling Collective dynamics Harvard, January / 25

32 Extended Dicke model [Baumann et al. Nature 1] z x κ Pump κ gψ Ω 2 Level System 2 Level system,, : : Ψ(x, z) = 1 : Ψ(x, z) = ω = 2ω recoil e ik(σx+σ z) σ,σ =± Feedback: U g2 ω c ω a H = ωψ ψ + ω S z + g(ψ + ψ )(S + S + )+US z ψ ψ. t ρ = i[h, ρ] κ(ψ ψρ 2ψρψ + ρψ ψ) Semiclassical EOM ( S = N/2 1) Ṡ = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω+us z )] ψ ig(s + S + ) ω khz ω, κ, g N MHz. Jonathan Keeling Collective dynamics Harvard, January / 25

33 Outline 1 Introduction: Dicke model and superradiance Rayleigh scheme: Generalised Dicke model 2 Attractors of dynamics (fixed points) Phase diagram: comparison to experiment 3 Approach to attractors: timescales Why slow timescales emerge 4 Attractors of dynamics (oscillations) Reaching other parameter ranges Jonathan Keeling Collective dynamics Harvard, January / 25

34 Fixed points (steady states) = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) ψ =, S = (,, ±N/2) always a solution. If g > g c, ψ too A S y = I[S ] = B ψ = R[ψ] = Jonathan Keeling Collective dynamics Harvard, January / 25

35 Fixed points (steady states) = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) ψ =, S = (,, ±N/2) always a solution. If g > g c, ψ too A S y = I[S ] = B ψ = R[ψ] = S z S y S x Small g:, only. (ω = 3MHz, UN = 4MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

36 Fixed points (steady states) = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) ψ =, S = (,, ±N/2) always a solution. If g > g c, ψ too A S y = I[S ] = B ψ = R[ψ] = S z S y S x Small g:, only. Larger g: SR too. (ω = 3MHz, UN = 4MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

37 Steady state phase diagram UN=, κ= = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z ω SR = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) ḡ N See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 Jonathan Keeling Collective dynamics Harvard, January / 25

38 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) 4 UN= ω UN=, κ= ḡ N SR SR(A): S y = 2 SR -2 SR g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 Jonathan Keeling Collective dynamics Harvard, January / 25

39 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) 4 UN=-1 ω UN=, κ= ḡ N SR SR(A): S y = SRB + SR(B): ψ = g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 Jonathan Keeling Collective dynamics Harvard, January / 25

40 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) 4 UN= SRB g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 ω UN=, κ= ḡ N SR SR(A): S y = + SR(B): ψ = Jonathan Keeling Collective dynamics Harvard, January / 25

41 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) 4 UN=-2 ω UN=, κ= ḡ N SR SR(A): S y = SRB + SR(B): ψ = g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 Jonathan Keeling Collective dynamics Harvard, January / 25

42 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) 4 2 UN=-4 ω UN=, κ= ḡ N SR SR(A): S y = -2 + SRB + SR(B): ψ = g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 Jonathan Keeling Collective dynamics Harvard, January / 25

43 Steady state phase diagram = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z = ig(ψ + ψ )(S S + ) = [κ + i(ω+us z )] ψ ig(s + S + ) SRB SRB+ SRB+ UN=-4 SRB g N (MHz) See also Domokos and Ritsch PRL 2, Domokos et al. PRL 1 ω UN=, κ= ḡ N SR SR(A): S y = + SR(B): ψ = Jonathan Keeling Collective dynamics Harvard, January / 25

44 Comparison to experiment (ω p -ω c ) (2π MHz) g 2 N (MHz) 2 UN = 1MHz Adapted from: [Bhaseen et al. PRA 12] [Baumann et al Nature 1 ] ω = ω c ω p UN, UN = g 2 4(ω a ω c ) Jonathan Keeling Collective dynamics Harvard, January / 25

45 Outline 1 Introduction: Dicke model and superradiance Rayleigh scheme: Generalised Dicke model 2 Attractors of dynamics (fixed points) Phase diagram: comparison to experiment 3 Approach to attractors: timescales Why slow timescales emerge 4 Attractors of dynamics (oscillations) Reaching other parameter ranges Jonathan Keeling Collective dynamics Harvard, January / 25

46 Dynamics: Evolution from normal state 4 2 (ii) (i) UN=-4 + SRB -2 (iii) g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

47 Dynamics: Evolution from normal state Gray: S = ( N, N, N/2) Black: Wigner distribution of S, ψ (ii) (i) UN=-4 SRB (i) SR(A) ψ t (ms) -2 (iii) g N (MHz) Oscillations:.1ms Decay: 2ms Jonathan Keeling Collective dynamics Harvard, January / 25

48 Dynamics: Evolution from normal state Gray: S = ( N, N, N/2) Black: Wigner distribution of S, ψ (ii) g N (MHz) (i) (iii) Oscillations:.1ms Decay: 2ms,.1ms UN=-4 SRB (i) SR(A) (ii) SR(B) ψ 2 ψ t (ms) t (ms) Jonathan Keeling Collective dynamics Harvard, January / 25

49 Dynamics: Evolution from normal state Gray: S = ( N, N, N/2) Black: Wigner distribution of S, ψ (ii) g N (MHz) (i) (iii) Oscillations:.1ms Decay: 2ms,.1ms, 2ms UN=-4 SRB (i) SR(A) (ii) SR(B) (iii) SR(A) ψ 2 ψ 2 ψ t (ms) t (ms) t (ms) Jonathan Keeling Collective dynamics Harvard, January / 25

50 Dynamics: Evolution from normal state Gray: S = ( N, N, N/2) Black: Wigner distribution of S, ψ (ii) g N (MHz) (i) (iii) Oscillations:.1ms Decay: 2ms,.1ms, 2ms UN=-4 SRB (i) SR(A) (ii) SR(B) (iii) SR(A) ψ 2 ψ 2 ψ t (ms) t (ms) t (ms) Jonathan Keeling Collective dynamics Harvard, January / 25

51 Asymptotic state: Evolution from normal state (Near to experimental UN = 13MHz). All stable attractors: 4 UN= SRB g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

52 Asymptotic state: Evolution from normal state (Near to experimental UN = 13MHz). Starting from All stable attractors: Asymptotic state UN= SRB SRB ψ g N (MHz) g N (MHz) 1-1 Jonathan Keeling Collective dynamics Harvard, January / 25

53 Timescales for dynamics: What are they? Asymptotic state SRB g N (MHz) ψ No unstable directions Initial growth time One unstable direction Two unstable directions g N (MHz) 1s 1s 1ms 1ms 1ms 1µs 1µs Growth Most unstable eigenvalues near S = (,, N/2) Decay Slowest stable eigenvalues near final state Jonathan Keeling Collective dynamics Harvard, January / 25

54 Timescales for dynamics: What are they? Asymptotic state SRB g N (MHz) ψ No unstable directions Initial growth time One unstable direction Two unstable directions g N (MHz) 1s 1s 1ms 1ms 1ms 1µs 1µs Growth Most unstable eigenvalues near S = (,, N/2) Decay Slowest stable eigenvalues near final state Expand in ω /κ: Oscillations: ω, Decay: ω or ω 2 /κ Jonathan Keeling Collective dynamics Harvard, January / 25

55 Timescales for dynamics: What are they? Asymptotic state SRB g N (MHz) Growth Most unstable eigenvalues near S = (,, N/2) Decay Slowest stable eigenvalues near final state Expand in ω /κ: Oscillations: ω, Decay: ω or ω 2 /κ ψ No unstable directions Initial growth time One unstable direction Two unstable directions g N (MHz) Asymptotic decay time SRB g N (MHz) 1s 1s 1ms 1ms 1ms 1µs 1µs 1s 1s 1ms 1ms 1ms 1µs 1µs Jonathan Keeling Collective dynamics Harvard, January / 25

56 Timescales for dynamics: Consequences for experiment 4 Asymptotic state SRB 1 1 ψ g N (MHz) 1-1 Jonathan Keeling Collective dynamics Harvard, January / 25

57 Timescales for dynamics: Consequences for experiment 6 1ms sweep Asymptotic state SRB ψ g 2 N (MHz 2 ) g N (MHz) 1-1 Jonathan Keeling Collective dynamics Harvard, January / 25

58 Timescales for dynamics: Consequences for experiment 6 1ms sweep Asymptotic state SRB ψ g 2 N (MHz 2 ) ms sweep g N (MHz) g 2 N (MHz 2 ) Jonathan Keeling Collective dynamics Harvard, January / 25

59 Timescales for dynamics: Why so slow and varied? Suppose co- and counter-rotating terms differ a b Ω a Ωb g ψ g ψ H =... + g(ψ S + ψs + ) + g (ψ S + + ψs ) Level System SR(A) near phase boundary at small δg Critical slowing down SR(A), SR(B) continuously connect Jonathan Keeling Collective dynamics Harvard, January / 25

60 Timescales for dynamics: Why so slow and varied? Suppose co- and counter-rotating terms differ a b Ω a Ωb g ψ g ψ H =... + g(ψ S + ψs + ) + g (ψ S + + ψs ) Level System 4 UN= SRB g N (MHz) SR(A) near phase boundary at small δg Critical slowing down SR(A), SR(B) continuously connect Jonathan Keeling Collective dynamics Harvard, January / 25

61 Timescales for dynamics: Why so slow and varied? Suppose co- and counter-rotating terms differ a b Ω a Ωb g ψ g ψ H =... + g(ψ S + ψs + ) + g (ψ S + + ψs ) Level System δg = g g, 2ḡ = g + g 4 2 UN= ḡ N=1 + SRB g N (MHz) δg/ḡ SR(A) near phase boundary at small δg Critical slowing down SR(A), SR(B) continuously connect Jonathan Keeling Collective dynamics Harvard, January / 25

62 Timescales for dynamics: Why so slow and varied? Suppose co- and counter-rotating terms differ a b Ω a Ωb g ψ g ψ H =... + g(ψ S + ψs + ) + g (ψ S + + ψs ) Level System δg = g g, 2ḡ = g + g 4 2 UN= ḡ N=1 + SRB g N (MHz) δg/ḡ SR(A) near phase boundary at small δg Critical slowing down SR(A), SR(B) continuously connect Jonathan Keeling Collective dynamics Harvard, January / 25

63 Outline 1 Introduction: Dicke model and superradiance Rayleigh scheme: Generalised Dicke model 2 Attractors of dynamics (fixed points) Phase diagram: comparison to experiment 3 Approach to attractors: timescales Why slow timescales emerge 4 Attractors of dynamics (oscillations) Reaching other parameter ranges Jonathan Keeling Collective dynamics Harvard, January / 25

64 Regions without fixed points Changing U: gψ Ω U g2 ω c ω a 2 Level System 4 2 UN=-4 + SRB g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

65 Regions without fixed points g ψ Changing U: Ω U g2 ω c ω a 2 Level System 4 2 UN=-4 + SRB g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

66 Regions without fixed points g ψ Changing U: Ω U g2 ω c ω a 2 Level System 4 UN= 2 SR -2 SR g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

67 Regions without fixed points g ψ Changing U: Ω U g2 ω c ω a 2 Level System 4 UN= g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

68 Regions without fixed points g ψ Changing U: Ω U g2 ω c ω a 2 Level System 4 2 UN= g N (MHz) Jonathan Keeling Collective dynamics Harvard, January / 25

69 Regions without fixed points g ψ Changing U: Ω U g2 ω c ω a 2 Level System 4 UN=4 2 Persistent Oscillations g N (MHz) ψ t (ms) Jonathan Keeling Collective dynamics Harvard, January / 25

70 Persistent (optomechanical) oscillations ψ t (ms) Ṡ = i(ω + U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω + US z )] ψ ig(s + S + ) Jonathan Keeling Collective dynamics Harvard, January / 25

71 Persistent (optomechanical) oscillations ψ t (ms) Ṡ = i(ω + U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω + US z )] ψ ig(s + S + ) Fix ω + US z = if ψ =. Jonathan Keeling Collective dynamics Harvard, January / 25

72 Persistent (optomechanical) oscillations ψ t (ms) Ṡ = i(ω + U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω + US z )] ψ ig(s + S + ) Fix ω + US z = if ψ =. Jonathan Keeling Collective dynamics Harvard, January / 25

73 Persistent (optomechanical) oscillations ψ t (ms) Fix ω + US z = if ψ =. N S = re iθ 2 r = 4 (Sz ) 2 Ṡ = i(ω + U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω + US z )] ψ ig(s + S + ) Get: θ = ω + U ψ 2 ψ + κψ = 2igr cos(θ) Jonathan Keeling Collective dynamics Harvard, January / 25

74 Persistent (optomechanical) oscillations ψ t (ms) Fix ω + US z = if ψ =. N S = re iθ 2 r = 4 (Sz ) 2 ψ ψ 2 S x S y S z Ṡ = i(ω + U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω + US z )] ψ ig(s + S + ) Get: θ = ω + U ψ 2 ψ + κψ = 2igr cos(θ) t(ms) Jonathan Keeling Collective dynamics Harvard, January / S x, S y, S z

75 Tuning g, g, U [Dimer et al. Phys. Rev. A. (27)] Ω a a b Ω Separate pump strength/detuning b g ψ g ψ 2 2 Level System g g Ω b, g g Ω a b, U g g a a 2 b Jonathan Keeling Collective dynamics Harvard, January / 25

76 Tuning g, g, U [Dimer et al. Phys. Rev. A. (27)] Ω a a b Ω Separate pump strength/detuning b g ψ g ψ 2 2 Level System g g Ω b, g g Ω a b Possible realization: Hyperfine levels, U g g a a 2 b σ σ + m = 1 F m = F m =+1 F Jonathan Keeling Collective dynamics Harvard, January / 25

77 Tuning g, g, U [Dimer et al. Phys. Rev. A. (27)] Ω a a b Ω Separate pump strength/detuning b g ψ g ψ 2 2 Level System g g Ω b, g g Ω a b Possible realization: Hyperfine levels, U g g a a 2 b σ σ + Atoms m = 1 F m = F m =+1 F B Jonathan Keeling Collective dynamics Harvard, January / 25

78 Summary Wide variety of dynamical phases 4 2 UN=4 4 2 UN= SR 4 2 UN= ḡ N=1 + SRB -2-2 SR g N (MHz) g N (MHz) g N (MHz) δg/ḡ Slow dynamics ψ t (ms) 4 Initial growth time No unstable directions 2 One unstable direction -2 Two unstable directions g N (MHz) 1s 1s 1ms 1ms 1ms 1µs 1µs 4 Asymptotic decay time 2 SRB g N (MHz) 1s 1s 1ms 1ms 1ms 1µs 1µs Persistent oscillations if U > ψ t (ms) ψ 2 12 ψ 2 S.4 1 x S y S 8 z t(ms) JK et al. PRL 1, Bhaseen et al. PRA 12 S x, S y, S z Jonathan Keeling Collective dynamics Harvard, January / 25

79 Jonathan Keeling Collective dynamics Harvard, January / 3

80 Extra slides 5 Ferroelectric phase transition 6 Lipkin-Meshkov-Glick 7 Tricritical point Jonathan Keeling Collective dynamics Harvard, January / 3

81 Ferroelectric transition Atoms in Coulomb gauge H = ω k a k a k + i [p i ea(r i )] 2 + V coul Jonathan Keeling Collective dynamics Harvard, January / 3

82 Ferroelectric transition Atoms in Coulomb gauge H = ω k a k a k + i [p i ea(r i )] 2 + V coul Two-level systems dipole-dipole coupling H = ω S z + ωψ ψ + g(s + + S )(ψ + ψ ) + Nζ(ψ + ψ ) 2 η(s + S ) 2 (nb g 2, ζ, η 1/V ). Jonathan Keeling Collective dynamics Harvard, January / 3

83 Ferroelectric transition Atoms in Coulomb gauge H = ω k a k a k + i [p i ea(r i )] 2 + V coul Two-level systems dipole-dipole coupling H = ω S z + ωψ ψ + g(s + + S )(ψ + ψ ) + Nζ(ψ + ψ ) 2 η(s + S ) 2 (nb g 2, ζ, η 1/V ). Ferroelectric polarisation if ω < 2ηN Jonathan Keeling Collective dynamics Harvard, January / 3

84 Ferroelectric transition Atoms in Coulomb gauge H = ω k a k a k + i [p i ea(r i )] 2 + V coul Two-level systems dipole-dipole coupling H = ω S z + ωψ ψ + g(s + + S )(ψ + ψ ) + Nζ(ψ + ψ ) 2 η(s + S ) 2 (nb g 2 Ferroelectric polarisation if ω, ζ, η 1/V ). < 2ηN Gauge transform to dipole gauge D r H = ω S z + ωψ ψ + ḡ(s + S )(ψ ψ ) Dicke transition at ω < Nḡ 2 /ω 2ηN But, ψ describes electric displacement Jonathan Keeling Collective dynamics Harvard, January / 3

85 Timescales for dynamics: U = and Lipkin-Meshkov-Glick Since κ ω, can consider eliminating ψ If U =, simple result: Ṡ = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω+us z )] ψ ig(s + S + ) t S = {S, H} ΓS (S ẑ), H = ω S z Λ + S 2 x Λ S 2 y with: Λ ± ω (g ± g ) 2, Γ 2κ (g 2 g 2 ) κ 2 +ω 2 κ 2 +ω 2 NB, g = g, no dissipation. Dissipation restored by finite κ. Jonathan Keeling Collective dynamics Harvard, January / 3

86 Timescales for dynamics: U = and Lipkin-Meshkov-Glick Since κ ω, can consider eliminating ψ If U =, simple result: Ṡ = i(ω +U ψ 2 )S + 2ig(ψ + ψ )S z Ṡ z = ig(ψ + ψ )(S S + ) ψ = [κ + i(ω+us z )] ψ ig(s + S + ) t S = {S, H} ΓS (S ẑ), H = ω S z Λ + S 2 x Λ S 2 y with: Λ ± ω (g ± g ) 2, Γ 2κ (g 2 g 2 ) κ 2 +ω 2 κ 2 +ω 2 NB, g = g, no dissipation. Dissipation restored by finite κ. Jonathan Keeling Collective dynamics Harvard, January / 3

87 Tricritical point UN=-4 SRB g N (MHz) SRB SRB SRB SRB g N (MHz) If ω u = UN/2 > κ, crossing of boundaries at: ω = ωu 2 κ 2 g ω UN N = 4 For g > g, 2 region exists at edge of SRB. ψ 2 Arg(ψ) S z /N 15 (a) 1 5 π (c) π/2 -π/2 -π.5 (e) g N =1. (MHz) SRB (b) 2 SRB (d) (f) π/16 π/2 7π/ Jonathan Keeling Collective dynamics Harvard, January / 3