Transcript arXiv

《第六届全国冷原子物理和量子信息青年学者学术讨论会》
2012年8月14日-18日, 浙江师范大学
Exact solvability and unified analytical treatments
to qubit-oscillator system
Qing-Hu Chen
(陈 庆 虎 )
Center for Statistical and Theoretical Condensed Matter Physics,
Zhejiang Normal University, Jinhua 321004, China
& Department of Physics, Zhejiang University, Hangzhou 310027, China
arXiv: 1204.3668, Phys. Rev. A, in press
arXiv:
1204.0953
本人
1966
年出生， 早就不属于青年学者
arXiv: 1203.2410向青年朋友请教了
Collaborators
Prof. Ke-Lin Wang
Department of Modern Physics, University of Science and Technology of China, Hefei
230026
Chen Wang (Ph. D student)
Department of Physics, Zhejiang University, Hangzhou 310027
Dr. Yu-Yu Zhang (Former Ph. D student)
Center for Modern Physics, Chongqing University， Congqing 400044
Shu He ( MS student), Prof. Tao Liu
School of Science, Southwest University of Science and Technology, Mianyang 621010
Prof. Shi-Yao Zhu
Beijing Computational Science Research Center, Beijing 100084
Brief introduction to quantum Rabi model (QRM)
Rabi, Phys. Rev. 49, 324 (1936); 51, 652 (1937).
□ The interaction of two-level atom (qubit) with a bosonic mode
H

 z   a  a  g (a   a ) x
2
ω is the resonant frequency of the cavity, Δ is the the transition frequency of the qubit,
and g is the coupling strength, σx,z is usual Pauli matrix, a(a+) is the boson
annihilation (creation) operator. δ=Δ- ω is the detuning.
quantum Rabi model
qubit-oscillator system
(Cavity QED)
(Circuit QED)
In the fully quantum mechanical version
Analytically unsolvable !
□ Jaynes-Cummings (JC) model (1963) under the rotating-wave approximation
(RWA) is analytically solvable. The counter rotating terms (CRTs) is omitted

H   z   a  a  g (a    a  )
2
RWA
 g (a   a  )
CRTs
The Rabi model (RM) describes the simplest interaction
between light and matter.
Although this model has had an impressive impact on many
fields of physics --- many physicists may be surprised to know
that the quantum Rabi model has never been solved exactly. In
other words, it has not been possible to write a closed-form,
analytical solution for it
Outline
1. Exact solution for Qubit-Oscillator Systems
(1) Numerically exact
(2) Analytically exact
(3) Applications
No explicitly expression
E ( n)  E(n,  , , , g )
2. Unified analytical treatments to qubit-oscillator systems
explicitly expression but complicated
3. Concise first-order corrections to the RWA
explicitly expression but very simple
Part I, Exact solution to the Quantum Rabi model (QRM)
□ In RWA, the n-th
1 1
E1,n   (n  )  Rn
2 2
1 1
E2,n   (n  )  Rn
2 2
an | n  

eigenstate is | 

n 
b
|
n

1

 n

at resonance, δ=0
     , Rn   2  4 g 2 (n  1)
  sin  n | n 
|1,n  

cos

|
n

1

n


 cos  n | n 
| 2,n  

sin

|
n

1

n


cos  n 
 Rn   
1
E1,n   (n  )  g n  1
2
1
E2,n   (n  )  g n  1
2
1   | n 
|1,n 


2  | n  1 
| 2,n
2g n 1
2
n  0,1, 2.
 4 g (n  1)
2
1  | n 



2  | n  1 
The ground-state

E0  
2
| 0 | 0, 
Vacuum Rabi splitting in the JC model
The atom is excited by the operator
V  e
GS
e
 g, 0
V GS


g  g
Measured transmission spectrum
showing the vacuum Rabi mode splitting
 e, 0
e, 0  g , 0
e, 0  g , 0

/2
/2
spontaneous emission to GS state
The emission spectrum has two peaks with
equal height (the distance of the two peaks,
2g, is the vacuum Rabi splitting).
2g is the energy difference of the 1st and 2nd
eigenstates
Wallraff et al., Nature 431,
162(2004).
The collapses and revivals in the evolution of the atomic
population inversion
 a(t ) | n 
| t  

b
(
t
)
|
n

1



If initially in Photonic Fock state |e, n>
This is the quantum Rabi oscillation.
b(t )  sin 2
2
If initially in photonic coherent state
| 0 (0) |  | g  e
 a  2 / 2

 gt
| 0 | g 
n   | a a |  |  |2
Population inversion under RWA can be evaluated analytically
p(t )  e| |
2

n  1
a(t )  cos 2 gt n  1
2
|  |2 n
n n! cos(2 gt n  1)  F (rwat )
M. O. Scully and M. S. Zubairy, Quantum Optics,
Cambridge University Press, Cambridge, 1997
collapses
<σz(t)>
revivals
Strong coupling Qubit-Oscillator System
Deppe et al., Nature physics 4, 686(2008)
Circuit quantum electrodynamics (QED) system
Experiments: T. Niemczyk et al., Nature Physics 6, 772 (2010)
Title: Beyond the Jaynes-Cummings model: circuit QED in the ultrastrong coupling regime
FIG. 1: Quantum circuit and experimental setup.
P. Forn-Diaz and J. E. Mooij et al., PRL 105, 237001 (2010).
arXiv: 1005.1559.
Spectrum of the flux qubit coupled to the LC resonator.
Theory of Qubit-Oscillator Systems (Biased QRM)
QHC, Tao Liu, and Kelin Wang, arXiv: 1007.1747,
Chin. Phys. Lett. 29, 014208 (2012)
The flux qubit behaves effectively as a two-level system
Hq    z   x  / 2
Δ is the tunnel coupling
The model Hamiltonian can be expressed as
H 
q
2
(cos  z  sin  x )   a  a  g (a   a) z
ωq is the atomic Larmor frequency, ω is the cavity frequency.
q   2  (2 I p   x )2
sin    q  
 2  (2 I p   x ) 2
g is the qubit-resonator coupling strength, enhanced by Josephson junction inductance
Circuit QED
Cavity QED:
g/ω ~
10-6
10-3
0.01
0.1
Nature 431, 162 (2004). Nature Physics 4, 686 (2008) Nature Physics 6, 772 (2010).
Numerically exact solution to QRM for ε≠0
QHC et al, arXiv: 1007.1747
A  a ,
B  a 
Transformation (ω=1)
g
 A A  g 2  
H 
 2




2
B B g  
 2
Ansatz for the wavefunction
  cn n
A
   N n 0 n
tr

  n 0  1 d n n
N tr



B


2
m

g

cm   n Dmn d n  Ecm


2


2
m

g

d m   n Dmn cn  Ed m


2
n
0
(δΦ≠0)
A
A

 A 
n
n!
e
1
 g 2  ga 
2
Dmn 
B
0
mn
A
0

 a   
n!
0
A
a

m
 1 B m n
2
A
n
A
exp  2 g 2 
m! nm
L m (4 g 2 )
n!
Laguerre polynomial
  2g 
nm
◇ For strong coupling or highly excited states, much better than exact diagonalization in a-space
x  6 6
The optimum fitted parameters of the experimental results
I p  515nA
g 2  0.82GHz
q   2  (2 I p   x )2
r 2  8.13GHz
sin    q  
 / h  4.25GHz
 2  (2 I p   x )2
□ The theoretical results are in good agreement with the experimental observations
 x  0

H    x  a  a  g (a   a ) z
2
Hamiltonian
Parity opertor
Numerically exact solution to unbiased QRM: ε=0
QHC et al, PRA 82, 052306(2010)
i ( a  a  z 2 1 2)
e
The system is of even (+) or odd (-) parity.
dn=±cn
The wavefunction is reduced to

  Ntr cn n
A
n0

  Ntr (1) n c n
n
  n0
(m  g  E
2
S-equation
[ H , ]  0
()
)cm

Ntr
n 0
The level transition is only allowed between the even and odd parity
E
()
E(
)



B
Dmn cn  0
First application to the entanglement dynamics
□ Entanglement : highly nonlocal—
shared among pairs of atoms,
photons, electrons, etc.,
they may be remotely located and
not interacting with each other.
□ Entanglement as a resource in
Yu and Eberly
new approaches to both
Science 2009
computation and communication
CAB in two identical JC atoms without RWA
(a) Bell state 1 for α=π/4
(b) Bell state 2 for α=π/12
□ Initiated from Bell state 1, ESD appear in
non-RWA; disappear in RWA
□ No periodicity of entanglement evolution
for large g
Effect of photonic number
on ESD
□ CAB and Nph
opposite behavior
□ Nph suppress CAB
Bell state 1
Possible origin of ESD
Bell state 2
Our analytically exact solution to the Rabi model
QHC, Tao Liu, Yu-Yu Zhang, and Ke-Lin Wang, EPL 96, 14003 (2011)
arXiv: 1011.3280
The original JC model
H JC
Parity

  z   a  a  g (a   a ) x
2
 JC  e
i ( a  a  z / 21/ 2)
 ( z )e
i a  a
[ H JC ,  JC ]  0
1  1 1


V
2   1 1
To solve it easily, we introduce

H  VH JCV    x   a  a  g (a   a ) z
2

  V JCV   x e

i a  a
[ H , ]  0
One dimensional large Frolich polarons
Qing-Hu Chen, Kelin Wang, and Shaolong Wan, J. Phys.: Condens. Matter 6, 6599(1994)
The wavefunction
Ntr=2
Ntr=3
where
to find the solution for α(q)
Ntr=2
Ntr=3
The wavefunction

  Ntr cn (a  ) n exp( a  ) 0 
n 0


N
tr
 n



c
(

a
)
exp(


a
)
0

n
n 0


Eigenstate of the pairty:
+1-> even parity
- 1-> odd parity
Schroedinger equation
    
H    E ( )  





a
a

g
(
a

a
)


2




 a  a  g (a   a
2

  Ntr cn (a  ) n ea  0 

n 0
 E () 

  Ntr c ( a  ) n e a 0 
 n  0 n


  Ntr c (a  ) n ea  0 

n 0 n


N tr
 n a

c
(

a
) e
0 

)  n0 n

Analytical solution
Comparing with the coefficients

  Ntr cn (a  ) n exp( a  ) 0 
n 0


  Ntr c (a  ) n exp( a  ) 0 
  n 0 n

(a )m exp( a ) 0
 m (1) m (2 ) j
(m  g )cm  (  g )cm1  (m  1) gcm1  
cm j  E (  ) cm
2 j 0
j!
Without normalization, we can set c0=1
The linear term in a+ in the Fock space can be determined by eigenvalue of the
pure coherent state exp(αa+ )|0>, we can set c1=0
energy
E ( )  g   / 2
Recurrence equation
j
m

1 


(2

)
m
cm1 
cm j  m  1...Ntr  1
(m  )cm  (  g )cm1 (1)

(m  1) g 
2
2 j 0 j !

Analytical solution
ci (i  0,1, 2,...Ntr )
can be expressed by α,
cN
tr 1
0
Then we have (g, ∆,ωare given model parameters)
Ntr
j


(2

)
f ( )  ( Ntr  )cNtr  (  g )cNtr 1 (1) Ntr 
cNtr  j  0
2
2 j 0 j !
A polynomial equation for only one variable α,
only real value of αis reasonable.
Zeros of f(α) give α
Eigenenergies
Eeigenfunctions
E (  )  g

/2
  Ntr cn (a  ) n exp( a  ) 0 
n 0


  Ntr c (a  ) n exp( a  ) 0 
  n 0 n

First order approximation (Ntr=2)
Ground-state energy
1  
4g 2
E0  
1  1 
2 
(1  ) 2
 

 2

Irish PRL07: Generalized RWA,
ED: exact diagonalization
in Bosonic Fock space
◇ Much better than GRWA method
Full solutions
Even parity, Ntr=59
y=f(α)
odd parity, Ntr=60
convergence

  Ntr cn (a  ) n exp( a  ) 0 
n 0


N
tr
 n



c
(

a
)
exp(


a
)
0

n
n 0


Remarks on QHC et al, EPL 96, 14003 (2011)
arXiv: 1011.3280
The JC model without the RWA can be mapped to
a polynomial equation with a single variable.
Its solutions recover exactly all eigenvalues and eigenfunctions
of the model for all coupling strengths and detunings.
In the past 80 years, it is analytically unsolvable.
4 months later, what happen?
See also arXiv: 1103.2461
E. Solano, Physics 4, 68 (2011).
Supplemental Material to Braak PRL paper
use the representation of bosonic creation and anihilation operators in the Bargmann
space of analytical functions in a complex variable z
R. Koc et al., J. Phys. A: Math.
Gen. 35, 9425(2002).
transcendental function
Zeros of the G-function gives
all eigenenergies with parity ±1
arXiv:1204.3856
A recent work by Braak has renewed the interest in the old problem of coupling a
photon field to a single spin 1/2 state, using the Rabi model. The central statement of
this work is that the eigenfunctions in Bargmann representation must be analytic
functions in the entire complex plane. Based on this condition, a procedure is derived
from the series expansion of the eigenstates which provides a recursive evaluation of
the spectrum. …
In the following, it is shown that the use of the extra condition of analyticity of
the eigenfunction in Bargmann representation is not necessary.
Travenec: PRA 85, 043805 (2012)
Nevertheless, there are disputes on whether the term exact solvability should be used,
if the G functions are given only by some Taylor expansions with coefficients coming
from a recurrence scheme.
In my opinion, the word integrability should be left rather for models where a
sufficient number of integrals of motion are known, which is not the case for Rabi
models.
Exact solvability of the quantum Rabi models within Bogoliubov operators
QHC et al., Phys. Rev. A ( in press, 2012), see also arXiv:1204.3668
A proof of our forthcoming article
Qing-Hu Chen et al., Phys. Rev. A ( in press, 2012)
1. Re-derivation of Braak's solution in a physical way!
two Bogoliubov transformations
If both wavefunctions are the true eigen-function for a non-degenerate eigenstate with
eigenvalue E, they should be in principle only different by a complex constant r
where
Final G- function
Braak, PRL 2011
ε=0
□ Within extended coherent states, a recent exact solution to the quantum Rabi model
[Daniel Braak, Phys. Rev. Lett. 107, 100401(2011)] can be recovered in an alternative
simpler and more physical way, without uses of any extra conditions.
Dear Qinghu,
Many thanks for your interesting paper on the derivation of the G-function for the
generalized Rabi model. Your method is certainly simpler than my approach, which was
based on symmetry considerations,......
The question is: how do you justify this condition of proportionality which in turn gives
the spectrum? In the Bargmann space approach it is the condition of analyticity throughout
the whole complex plane which determines whether a state
is an element of the Hilbert space or not - and this leads together with Z_2-symmetry to the
G-function. What corresponds to the Bargmann condition in your approach? .
Best regards,
Daniel
the extra condition is covered in the vacuum state in the space
of the Bogoliubov operators. These vacuum states are well defined
and known as the coherent states, so the present derivation is more
physical and simpler.
The Juddian solutions originate from the properties of degeneracy
□ Both Hilbert spaces in the two Bogoliubov operators are complete, if
truncation is not done, the proportionality is justified naturally for nondegenerate states.
□ Koc et al in J. Phys. A: Math. Gen. 35, 9425(2002) have obtained isolated
exact solutions in the QRM, which are just the Juddian solutions
with doubly degenerate eigenvalues.
The degenerate eigenstates are excluded in principle in the solutions based on
the proportionality. It naturally follows that the Juddian solutions are exceptional
ones
is not analytic in x but has simple poles
at x=0,1,2….
All exceptional eigenvalues have the form En=n-g2
the necessary and sufficient condition for the occurrence of this eigenvalue is
fn(x)=0
Comparisons with our previous work I
□ Qing-Hu Chen, Tao Liu, Yuan Yang, and Kelin Wang, PRA 82, 052306(2010)
Qing-Hu Chen, Lei Li, Tao Liu, and Kelin Wang, arXiv: 1007.1747
  Ntr cn n
   Nn 0
tr

  n0 d n n



B
A
Braak, PRL 107, 100401(2011). His wavefunctions are unfold by us
  Ntr n !en n
   Nn 0
tr


 n0 n ! f n n

A


A
  Ntr  1n n ! f 'n n 
B

   nN0
n
tr



1
n
!
e
'
n



n
B 
 n 0
link coefficients in two ansatz of the
wavefunction
◇ No essential differences, except that the avenues to obtain the basically same
coefficients are different !
Comparisons with our previous work II
□ Qing-Hu Chen, Tao Liu, Yu-Yu Zhang, and Ke-Lin Wang, EPL 96, 14003 (2011)
arXiv:1011.3280
Ntr

(2 ) j
Ntr 
f ( )  ( Ntr  )cNtr  (  g )cNtr 1 (1)
cNtr  j  0

2
2 j 0 j !

Braak’s G function
j
(2

)
f  ( )  
cNtr  j  0
j!
j 0
Ntr
Reduced to
E ( )  g  / 2
◇ zeros of the both functions defined through different power series can give the
exact eigenvalues. Both are analytically exact solutions
◇ In our practical evaluation, it is not more difficult to locate the zeros for our
function than those proposed by Braak, because the poles at x = n emerging in the
latter are not present in our earlier solution.
Both are not analytical closed-form solutions!
□ Braak’s solution
the eigenvalues are given by the zeros of the above Heun functions, which
can not be obtained without truncation in the power series.
◇ the expansion can not be closed naturally like in the JC model under the RWA
  Ntr n !en n
   Nn 0
tr

  n0 n ! f n n



A
A
□ Enrique Solano called on in viewpoint Physics 4, 68 (2011):
An intense dialogue between mathematics and physics will be needed to describe and
predict unprecedented physical phenomena, since they might be hidden in the quantum
numbers associated with Braak’s integrability criterion and analytical expressions.
Otherwise, the present achievement will remain but a mathematical monologue.
◇ by our two theories, we can also describe any physical phenomena based on
Braak’s solution. Until now, no exceptions!
Wolf, Kollar, and Braak, PRA 85, 053817 (2012), see also arXiv:1203.6039
Braak’s method is absolutely not the unique method to obtain their physical figures 4-7
By our two earlier methods, we can completely reproduce all their figures
Ntr=36
It is funny that we have studied the dynamics much earlier, which is however still only on arXiv
Yu-Yu Zhang, Qing-Hu Chen, Shi-Yao Zhu, arXiv: 1106.2191,
Wolf, Kollar, and Braak, PRA 85, 053817 (2012), see also arXiv: 1203.6039
□ Counter rotating effect on collapses and revivals
(Yu-Yu Zhang, QHC, Shi-Yao Zhu, arXiv: 1106.2191 )
| 0 (0) |  | g  e
Population inversion
i.e. <σz(t)>
Evolution at different
g
 a   2 /2
2
| 0 | g  ， n    10
P(t) is periodic with 2π，
g=2， periodicity τ0= 4gπ=25.132
the mismatched phase in cos function in the summation at other time ,
leading to collapses.
□ Instead of studying the evolution of < σx> ,
we studied the vacuum Rabi splitting
(Yu-Yu Zhang, QHC, Shi-Yao Zhu, arXiv: 1106.2191 )
Emission spectrum for g=0.1 and 0.8 at resonance
Different heights
More than 2 peaks
V|GS>
No splitting
Remarks for, Yu-Yu Zhang,
QHC, Shi-Yao Zhu, arXiv: 1106.2191
□ It is predicted that two quantum phenomena are fundamentally altered.
◇ In the VRS:
different heights of peaks
more than two peaks
the splittings completely vanish in the strong coupling regime.
◇ The collapses and revivals disappear in the intermediate g
but reappear periodically in the strong g
different from aperiodic collapses and revivals in RWA.
□ Braak et. al. discussed the similar dynamics in PRA 85, 053817 (2012) ( arXiv:
1203.6039) .
◇ Their results for the evolution of < σz > are not beyond ours!
◇ Instead of studying the evolution of < σx> as it is related the emission, we
studied the vacuum Rabi splitting where new result is given and which is easy to be
observed in experiments.
Exact solvability of the two-photon quantum Rabi models
within Bogoliubov operators
Qing-Hu Chen et al., Phys. Rev. A ( in press, 2012), see also arXiv:1204.3668
Bogoliubov transformation
Similarly, the G-function for two-photon QRM is obtained
Tutorial for Bogoliubov operators approach
◇ find the corresponding Bogoliubov operators to remove the linear terms.
◇ expand the wavefunctions in terms of each Bogoliubov operator respectively.
◇ Eliminating the ratio constant of these wavefunctions will give the
transcendental functions, which would be defined through power series in model
parameter dependent quantities with coefficients related recursively.
◇ Finally, zeros of these transcendental functions would give the eigenvalues
exactly, where numerical solutions to the one-variable (or finite variables in other
multi-level systems for example) nonlinear equation must be required.
It can be encouraged to extended to various spin-boson systems with
multi-level, even multi-mode.
Part II, Unified analytical treatments to qubit-oscillator
Shu He, Yu-Yu Zhang, and QHC et al, arXiv: 1204.0953
□ Hamiltonian for the whole system reads
 

H     z   x    a  a  g (a   a) z
2 
2
◇ Goal: without any numerical effort
E ( n)  E(n,  , , , g )
□ Recent most cited work
E. K. Irish, PRL 99, 173601 (2007)
Generalized rotating-wave approximation (GRWA)
ε=0
E. Solano’s group: PRL105, 263603(2010),
an expansion in deep strong coupling regime (DSC) ε=0
M. Grifoni’s group, Phys. Rev. A 82, 062320 (2010).
Van Vleck perturbation theory (VVP)
ε≠0
…..
Recall: Our previous numerically exact solution:
Qing-Hu Chen, Tao Liu, and Kelin Wang, arXiv: 1007.1747,
Chin. Phys. Lett. 29, 014208 (2012)
After transformation
A  a ,
B  a 
 A A  g 2  
H 
 2

g



2
B B g  
 2
Ansatz for the wavefunction
  Ntr cn n
A
   N n 0 n
tr


 n 0  1 d n n



B
n
0
A
A

 A 
e
n
n!
0
1
 g 2  ga 
2
A
0

 a   
n!
n
0
A
a


2
m

g

cm   n Dmn (2 )d n  Ecm


2


2
m

g

d m   n Dmn (2 )cn  Ed m


2
□ We alternatively present some analytical results in the framework of the above
Formalism, and recover many previous analytical results!
Variational study for ε=0
g to relaxed to a variational parameter α, the trial state is chosen as the most simple state
Ntr=0, i.e. the vacuum state
 0 A
  

0
 B
0
A
e
1
 g 2  a
2
0
a
0
B
e
◇ The energy expectation
◇ Minimizing the energy gives
◇ In the weak coupling limit, we can obtain α and the ground state energy
Which are exactly the same as Eqs. (7) and (8) obtained in
Yuanwei Zhang et al., PRA 83, 065802 (2011)
◇ In the strong coupling limit, we can obtain α and the ground state energy
  g, E0   g 2
1
 g 2  a
2
0
a
Recovery of E. Solano’s group: PRL 105, 263603(2010), ε=0
Perturbation theory based on the exact solution in the strong coupling limit
  a  a g (a   a)

0
neglecting the qubit tunneling term H
 a a  g (a  a) z

H0 0 

0

a
a

g
(
a

a
)

z
the exact eigenstates
the eigenvalues for the m state.
the second-order perturbative study. E. g. for even parity
Finally we have
Which are exactly the same as that in Eq. (5) in E. Solano’s group: PRL105, 263603(2010),
Recovery of M. Grifoni’s group, PRA 82, 062320 (2010). ε≠0
Van Vleck perturbation theory (VVP)
For any value of the qubit bias ε≠0 , the Hamiltonian with a vanishing tunneling
element Δ = 0 can be diagonalized in terms of two eigenstates
the corresponding eigenvalues are
For finite Δ, the perturbative matrix elements becomes
the full Hamiltonian can be diagonalized perturbatively to second-order in Δ by using
VVP theory
l is εdependent, ambiguity
Recovery of E. K. Irish, PRL 99, 173601 (2007), GRWA ε=0
(m  g  E
2
Eg. Ntr=4
()
)cm

Ntr
n 0
Dmn cn  0
  g 2 D00  E
D01
D02
D03

D10
1  g 2 D11  E
D12
D13


D20
D21
2  g 2 D22  E
D23

D30
D31
D32
3  g 2 D33  E


D40
D41
D42
D43
4  g2

◇ First-order approximation, the determinant takes the 2-by-2 block form
m  g 2 Dmm  E
Dm 1,m


D14


D24

D34

D44  E 
D04
Dm,m 1
0
m  1  g 2 Dm1,m1  E
give the GRWA results exactly
◇ Moreover, further analytical treatment is easily implemented in this scheme!
e.g. Second-order approximation, the determinant takes the 3-by-3 block form
Beyond the RWA, which can not be implemented within any renormalized RWA
Shu He, Yu-Yu Zhang, and QHC et al, arXiv: 1204.0953
Comparsions for energy levels for zero bias
   
The energy levels as a function of coupling constant for different qubit bias ε
The energy levels as a function of coupling constant for different qubit bias ε
The energy levels as a function of coupling constant for different qubit bias ε
Dynamics for qubit population difference for ε=0
Dynamics for qubit population difference for ε≠0
Summary for, Shu He, Yu-Yu Zhang, and QHC et al, arXiv: 1204.0953
◇ Many previous analytical treatments ( 2 PRL, 2 PRA) can be recovered in the
present scheme unifiedly.
◇ For ε=0, the GRWA is further improved to BRWA, which is more close to the exact
ones at large detuning while the GRWA deviates strongly.
◇ We extend the GRWA to the finite-bias case. The results is much better than VVP in
the weak and intermediate coupling regime, which is more experimentally interesting.
◇ At the experimentally accessible coupling regime, the dynamics of the qubit by
GRWA can always work well, but VVP theory is not valid in the present-day
experimentally interesting coupling regime, except for very large static biases.
Part. III, Concise first-order corrections to the RWA for ε=0
Shu He & Qing-Hu Chen et al, arXiv: 1203.2410
Many approaches are on the basis of various polaron-like transformations or shifted
operators, which are basically photonic coherent states approaches .
Eg. The wavefunction in GRWA
Although vevry accurate solutions
can be obtained, the infinite ‘bare’
states should be involved!

 m
 cm m A
 
 cm 1 m  1
A 


a


 m
m
0
A
A
m!
e
1
 g 2  ga 
2
0
A
0


A

 
m!
m
0
A
a
The original RWA results [from only two bare states] are missed or lost in these
approaches, so the effects of counter rotating-wave terms are not clear!
 cn | n  
|  

d
|
n

1

 n

Only two ‘bare’ states under RWA
We are not interested in the more accurate solutions with infinite ‘bare’ states now
What emerges first beyond the RWA?
The wavefunction in the photonic Fock states
S. Swain, J. Phys. A 6, 1919 (1973).
Transformed Hamiltonian

H    x   a  a  g (a   a ) z
2
Wavefunction with even (odd) parity

  M cn n

n 0


M
n

(

1)
cn n 

n 0

The parity is not considered by Swain, so he developed a continued fraction approach
very complicated!
First approximation
They are just the RWA energy levels
Transforming back to the original frame gives
which are just the eigenstates under the RWA
The second order approximation
the resonant case δ=0
◇ m = 2k+1 with even parity
and m = 2k with odd parity
◇ m = 2k with even parity
and m = 2k + 1 with odd parity
◇ Unified expressions
RWA results
Applications I: vacuum Rabi splitting
◇ Our analytical results for the main peak
height agree excellently with the exact ones in
a wide coupling regime (0 < g < 0.2).
Peak distance
Applications I: Berry phase
The GS Berry phase: Plastinaet al, EPL 76, 182 (2006).
◇ the analytical Berry
phase in the second-order
approximation agrees quite
well with the exact ones for
a wide coupling regime
(0 < g < 0.5).
Summary for Shu He & Qing-Hu Chen et al, arXiv: 1203.2410
□ by the bosonic Fock space and parity symmetrythe analytical results are presented at
different stages:
◇ The first approximation in the present formalism reproduces exactly the RWA results.
◇ The effect of the CRTs emerges clearly just in the second order approximation.
they play dominant role g <= 0.2,
CQED: gmax = 0.12
□ Effects of CRTs on various phenomena in the ultra-strong coupling regime.
Analytical results to the vacuum Rabi splitting and the Berry phase agree well with the
exact ones in a wide coupling regime.
Thank you very much for your attention