Accepted Manuscript The Generalized Coherent State Ansatz: Application to Quantum Electron-Vibrational Dynamics Raffaele Borrelli, Maxim F. Gelin PII: DOI: Reference:
S0301-0104(16)30271-3 http://dx.doi.org/10.1016/j.chemphys.2016.05.013 CHEMPH 9557
To appear in:
Chemical Physics
Please cite this article as: R. Borrelli, M.F. Gelin, The Generalized Coherent State Ansatz: Application to Quantum Electron-Vibrational Dynamics, Chemical Physics (2016), doi: http://dx.doi.org/10.1016/j.chemphys.2016.05.013
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
!" #"$"%&'()"* +,!"%"$- .-&-" !"#$%/ 011'(2&-(,$ -, 34&$-45 6'"2-%,$78(9%&-(,$&' :;$&5(2< !"!#$# %&''#$$( ) *!+(, -. /#$(0! !
!"#$#% &'()*+,(-. /( 01+('1% !234456 7+89:(;,<1% !-;:=> *?;+-*@*'- 1A BC*@(,-+=% 0*
!"#$%
1 0#2 !" #$ 3&' ,&$#45$!' 6(7'&0(4 2!6# 35048(&09 7!9#: &0 ! 95;#';&9(8(&0 &3 8(,#<:#;#0:#08 =#0#'!$(>#: 4&?#'#08 98!8#9 (9 :#6#$&;#: !0: !0!$@9#:. A?# ,#8?&:&$&=@ (9 9;#4(B4!$$@ 8!($&'#: 8& :#94'(7# 8?# 8(,# #6&$58(&0 &3 8?# 2!6# 35048(&0 &3 ! 9@98#, (0 2?(4? 9#6#'!$ (08#'!48(0= #$#48'&0(4 98!8#9 !'# 4&5;$#: 8& ! 7!8? &3 ?!',&0(4 &94($$!8&'9. A?# #C5!8(&09 &3 ,&8(&0 3&' 8?# 2!6# ;!4D#8 ;!'!,#8#'9 !'# &78!(0#: 7@ 59(0= 8?# E('!4<-'#0D#$ 8(,#< :#;#0:#08 6!'(!8(&0!$ ;'(04(;$#. A?# ,#8?&:&$&=@ (9 59#: 8& :#94'(7# 8?# C5!085, :@0!,(4!$ 7#?!6(&5' &3 ! ,&:#$ ;&$!'&0 9@98#, !0: (89 94!$(0= !0: 4&06#'=#04# ;'&;#'8(#9 !'# :(94599#: !0: 4&,;!'#: 2(8? 05,#'(4!$$@ #+!48 '#95$89.
'( )*#$+,-.+*
F0# &3 8?# 9(,;$#98 2!@9 8& :#94'(7# 8?# ,&8(&0 &3 054$#( (0 ,&$#45$!' 9@98#,9 (9 7!9#: &0 8?# 59# &3 4&,;$#+ /!599(!0 2!6# ;!4D#89) %&'& 4&< ?#'#08 98!8#9. G#$$#' 2!9 !,&0= 8?# B'98 8& '#4&=0(># 8?# ;'(,!'@ '&$# &3 954? ! 9(,;$# !;;'&+(,!8(&0 &3 8?# 054$#!' 2!6# 35048(&0 (0 8?# 985:@ &3 4?#,(4&<;?@9(4!$ ;'&7$#,9) ,!(0$@ 3&459(0= ?(9 2&'D &0 9#,(4$!99(4!$ 8?#&'(#9.HI) J) K) LM 18 8?# 9!,# 8(,# E!6@:&6 59#: 4&?#'#08 98!8#9 8& 985:@ #0#'=@ 8'!09;&'8 (0 7(&4?#,(4!$ 9@98#,9) 8?&5=? 2(8? ! :("#'#08 8?#&'#8< (4!$ !;;'&!4?.HNM O(04# 8?#0 ! 05,7#' &3 #+8#09(&09 &3 8?#9# ,#8?&:9 ?!9 7##0 ;'&;&9#: 3&' :#!$(0= 2(8? ,5$8(:(,#09(&0!$ 9@98#,9 !0: 0&0
#$$()*+(,'%(,
./*0*#1# 2!""#11(3
E;= 33% N43O
!"#$ "%&'() %&' *+,))-+# .+)'/0!,1%-2"#3$,(+%-"# %-!'0/'4'#/'#% 5+(%('' 6*0789:5;< %&' /-('2% /=#+!-2) >+(-+%-"#+1 !,1%-02"#3$,(+%-"#+1 *+,))-+# 6::0>7*8;< +#/ %&' +. -#-%-" !,1%-41' )4+?#-#$ 6 @7A;< +(' 'B0 %('!'1= $'#'(+1 +#/ 4"?'(C,1 !'%&"/"1"$-') %&+% ,%-1-D'< -# /-E'('#% ?+=)< 1-#'+( 2"!.-#+%-"#) "C *+,))-+# ?+>' 4+2F'%) %" /'3#' + %-!'0/'4'#/'#% )'% "C "(%&"#"(!+1 .+)-) C,#2%-"#) ?&-2& +(' %&'# ,)'/ %" /')2(-.' %&' ?+>'0 C,#2%-"# "C %&' '#%-(' #,21'+( )=)%'!GHI< J< K< L< MNO 8"&'('#% )%+%')< 2"!0 .-#'/ ?-%& + %-!'0/'4'#/'#% :PQ ('4(')'#%+%-"# +(' '!41"='/ -# !'%&"/0 "1"$-') ),2& +) %&' R"2+1 8"&'('#% A%+%' 44("B-!+%-"# 6R8A ;HMMO< +#/ %&' !,1%-2"#3$,(+%-"#+1 S&('#C')% 678S;GHMT< MUO 9&' )"02+11'/ :+>=/"> !" #$ +#/ -%) 'B%'#)-"#)< ?&-2& &+>' .''# 'B%'#)->'1= ,)'/ -# %&' )%,/= "C 4"1+("# +#/ )4-#0.")"# /=#+!-2) +(' +1)" .+)'/ "# %&' ,)' "C ),4'(4")-%-"# "C 2"&'('#% )%+%') HMV< MW< MI< MJ< MK< ML< TNO "( !,1%-41' 2"&'('#% )%+%')GHTMO 5'11'( ?+) +1)" %&' 3()% %" ),$$')% %&+% + )-$#-32+#% -!4(">'!'#% "C %&' )'!-21+))-2+1 /')2(-4%-"# .+)'/ "# *+,))-+# ?+>'4+2F'%) 2",1/ &+>' .''# ".%+-#'/ .= -#%("/,2-#$ + 4"1=#"!-+1 4('C+2%"( %" %&' ?+>' C,#2%-"#GHM< TO 5+$'/"(# 1+%'( ('2"$#-D'/ %&+% %&' #+%,(+1 'B%'#)-"# "C %&' /')2(-4%-"# "C X,+#%,! /=#+!-2) .= *+,))-+# ?+>'4+2F'% ?+) %&' ,)' "C $'#'(+1-D'/ *+,))05'(!-%' C,#2%-"#)HTT< TU< TV< TWO F#"?# -# %&' X,+#%,! "4%-2) 1-%'(0 +%,(' +) *'#'(+1-D'/ 8"&'('#% A%+%') 6*8A;HTIO< +#/ -# !"1'2,1+( /=#+!0 -2) +) 5+$'/"(#Y) ?+>'4+2F'%)GHTJ< TKO )-!-1+( %'2-X,' ?+) /'>'1"4'/ 1+%'( .= Z-11-#$HTL< UN< UM< UTO -# %&' C(+!'?"(F "C )'!-21+))-2+1 !'2&+#-2) +) ?'11 +) #"#0+/-+.+%-2 %(+#)-%-"#)GHUUO 9&",$& + $'#'(+1-D'/ C"(! "C #"# )4('+/-#$< %&'& 2"&'('#%< ?+>'04+2F'%) ?+) 3()% -#%("/,2'/ .= A'#-%DF= -# %&' 3C%-')HUVO< +) 4"-#%'/ ",% .= "%&'( +,%&"()' #"% ='% ('2'->'/ %&' /,' +%%'#%-"#G 5'(' ?' 4("4")' +#/ /'>'1"4 + #'? %=4' "C !" #$ %" ('4(')'#% %&' %-!' '>"1,%-"# "C !"1'2,1+( >-.("#-2 )%+%') ,)-#$ *8A)G ["11"?-#$ + >'(= ('2'#% +44("+2& C"( %&' )%,/= "C X,+#%,! /=#+!-2) "C 2&'!-2"04&=)-2+1 )=)%'!) HUIO ?' '!41"= %&' %-!'0/'4'#/'#% :-(+20[('#F'1 >+(-+%-"#+1 4(-#2-41' %" /'(->' + #'? )'% "C 'X,+%-"#) -# ?&-2& 4")-%-"# +#/ !"!'#%,! "C %&' 2"0 &'('#% )%+%') +(' /=#+!-2+11= 2",41'/ %" %&' '>"1,%-"# "C %&' 'B4+#)-"# 2"'\2-'#%) "C %&' .+)-) C,#2%-"#)G 9&' !'%&"/"1"$= -) +441-'/ %" %&' )%,/= "C + "#'0/-!'#)-"#+1 5"1)%'-# 4"1+("# !"/'1 ?&'(' N '1'2%("#-2 )-%') +(' 2",41'/ %" 1"2+1 4&"#"#)< +#/ -%) ('),1%) +(' 2"!4+('/ ?-%& #,!'(-2+11= 'B+2% 789:5 2+12,1+%-"#)GHUJO T
! "#$ %$&$'()*+$, -.#$'$&/ 0/(/$ 1$/#.,
!"# $%&'& ()'*+&( ,-' .%& /'#. .*0& *1 /,.* +2 "&1*.3425678 91( :9.&' ;&1&'9:*3&( +2 <*&.-56=8 *1 .%& ,-'0 -, #>?&&3&( 1?0+&' #.9.@ A1 .%&*' +9#*) ,-'0 !"# )91 +& -+.9*1&( +2 9BB:2*1; .%& .'91#:9.*-1 -B&'9.-' .- 912 %9'0-1*) -#)*::9.-' &*;&1#.9.& n 56C8 n, α = D(α) n DEF $%&'& a 91( a+
D(α) = exp(αa+ − α ¯ a),
9'& +-#-1 911*%*:9.*-1 91( )'&9.*-1 -B&'9.-'#G 91( α *# 9 )-0B:&H 1?0+&'@ I%&1 n = 0 #.91(9'( :9?+&' )-%&'&1. #.9. 9'& -+.9*1&(@57J8 K1:*4& .%& :9..&'G !"# 9'& 1-. &*;&1#.9. -, .%& 911*%*:9.*-1 -B&'9.-'@ A. )91 +& &9#*:2 L&'*/&( .%9. √ DMF a n, α = n n − 1, α + α n, α √ ¯ n, α . D6F a† n, α = n + 1 n + 1, α + α !" %9L& .%& )--'(*19.& '&B'&1.9.*-157E8 D7F $%&'& .%& .$- '&9: L9'*9+:G q, p 9'& '&:9.&( .- .%& )-0B:&H L9'*9+:& α +2 ϕn (x, q, p) = 2−n/2 (n!)−1/2 π −1/4 e−(x−q)
2 /2+ip(x−q/2))
Hn [(x − q)]
q + ip α= √ . 2
N*19::2G *. *# B-##*+:& .- ;*L& .%& 1?0+&' #.9.& &HB91#*-1 -, 9 !" *1 .%& ,-'0 ∞ −|α|2 /2 X n+k 1/2 n−k 2 n−k n, α = e (−1) (k!) L (|α |)(¯ α ) k k (n!)1/2 k=0
DOF
n−k Ln (·) +&*1; .%& 9##-)*9.&( P9;?&''& B-:21-0*9:# -, -'(&' n@ N-' n = 0G 0, α ≡ α G 91( $& '&)-L&' .%& ,-'0?:9 ,-' :9?+&' )-%&'&1. #.9.@57J8
I*.%-?. :--#*1; 912 ;&1&'9:*.2 $& )91 $'*.& .%& L*+'-1*) $9L& ,?1).*-1 *1 .%& ,-'0 X Ψ(t) = DQF ClK (t) K, αl (t) l l,K
6
!"#"
l $%"&'("$ ")"&*#+,'& $*-*"$ . l = 1, 2, ..., N /0 K '$ - 12)*'3',4"5 K = (k1 , k2 , ..., kd ) ∈ Nd0 0 -,4 *!" K, αl (t) -#" - $"* +6 !"#$%#'%#' 789 4"(,"4 -$ *!" %#+42&* +6 d +,"34'1",$'+,-) 789 d Y ki , αli (t) . K, αl (t) ≡ k1 · · · kd , αl1 (t), · · · , αld (t) =
.:/
i=1
;* '$ ")) <,+ , *!-* *!" ,21="# +6 &+">&'",*$ +6 "5%-,$'+, ? '$ $2=@"&* *+ -, "5%+,",*'-) A#+ *!0 -,4 &-,,+* =" 2$"4 6+# %#-&*'&-) %2#%+$"$ -$ $2&!B C, *!" +*!"# !-,40 =D ',*#+42&',A -, "5%)'&'* *'1"34"%",4",* =-$'$ 62,&*'+,$ -,4 2$',A *!" E'#-&3F#",<") G-#'-*'+,-) -%%#+-&! *!" H2-,*21 4D,-1'&-) %#+=)"1 '$ 6+#12)-*"4 +, - I')="#* =2,4)"0 -,4 *!" =-$'$ $"* '$ )+&-))D -4-%*"4 *+ %#+G'4" *!" J="$*K -%%#+5'1-*'+, *+ *!" 9&!#L4',A"# "H2-*'+, +, *!" %-#-1"*"#$ 1-,'6+)4B ;* '$ *!"#"6+#" %!D$'&-))D $+2,4 *!-* - *#2,&-*'+, +6 *!" -=+G" "5%-,$'+, *+ - #")-*'G")D $1-)) $2=$"* +6 $*-*"$ -)#"-4D %#+G'4"$ - A++4 4"$'%*'+, +6 *!" $*-*" +6 *!" $D$*"1B M!" -%%#+-&! 4"G")+%"4 ', *!'$ +#< '$ =-$"4 +, *!" '4"- +6 %-#*'*'+,',A *!" ",*'#" I')="#* $%-&" ', - $"* +6 $2=$%-&"$ !'&! 4'N"# ', *!" ,21="# +6 G'=#-*'+,$ !'&! -#" -))+ "4 *+ =" $'12)*-,"+2$)D "5&'*"4B M!2$ *!" ",*'#" I')="#* $%-&" H $%-,,"4 =D *!" I-1')*+,'-, +6 OHB PQ &-, =" %-#*'*'+,"4 -$ [ H =
Sc
c
!"#" Sc '$ *!" $%-&" $%-,,"4 =D *!" $*-*"$ ', !'&! +,)D c =-$'$ 62,&*'+,$ -#" $'12)*-,"+2$)D "5&'*"40 !(#( !-G" - ,+,3R"#+ G-)2" +6 kiB C6 &+2#$" 6+# -,D %#-&*'&-) %2#%+$"$ *!" =-$'$ '$ *#2,&-*"4 *+ - 1-5'121 H2-,*21 ,21="#B ! 9',&" ', - $D$*"1 '*! d G'=#-*'+,-) 4"A#""$ +6 6#""4+1 *!"#" ')) =" dc 4'$*',&* &+1=',-*'+,$ !-G',A c "5&'*"4 1+4"$0 2$',A $2&! - %-#*'*'+, *!" -G"62,&*'+, +6 OHB ? &-, =" #'**", ', *!" 6+#1 Ψ(x, t) =
N X l=1
(dc) d X X
X Clki1 ...kic (t) ki1 , αli1 (t) . . . kic , αlic (t) l
wi1 ...wic
c=0 i1 ...iC ki1 ...kic =1
wi N h d X X X = Cl (t) αl (t) + Clki (t) ki , αli (t) αl0 (t) l=1
() X wj wi X X d 2
+
i,j=1 ki =1 kj =1
i=1 ki =1
i Clki kj (t) ki , αli (t) kj , αlj (t) αl00 (t) + ... l ,
T
.S/
!"#" " !$%" &$#'('()*"+ '!" αl $, αl = (αl0 , αli) (* '!" -#,' ,.//$'()* $*+ αl = (αl00 , αli, αlj ) (* '!" ,"0)*+ ,.//$'()*1 2* $*$3)4)., &$#'('()* !$, 5""* &"#6)#/"+ 6)# '!" 7.$*'./ *./5"#, K 8 $*+ (*+(0", )6 '!" "9&$*,()* 0)":0("*', ('! %$3." ;"#) $#" )/(''"+< wi #"","*' '!" /$9(/./ 7.$*'./ *./5"# )* '!" i='! /)+"1 >,(*4 &$#'('()* ? (' (, &),,(53" ') #"+.0" '!" ,(;" )6 '!" 5$,(, ,"' 5@ 3(/('(*4 (', 0)/5(*$')#($3 "9&3),()*1 A*+""+8 (6 "$0! /)+" (, $33) "+ ') !$%" w 5$,(, 6.*0'()*, '!"* $ 0)/5(*$'()* ('! c ,(/.3'$*").,3@ ! d "90('"+ /)+", (33 (*'#)+.0" wc *" 5$,(, ,'$'",8 $*+ ,(*0" !'!"#" $#" c ,.0! 0)/5(*$'()*, '!" +(/"*,()* )6 '!" Sc ,.5,&$0" (33 5" dc wc1 B!.,8 6)# c = 1 " (33 !$%" &),,(53" 0)/5(*$'()*, $*+ dw 5$,(, ,'$'",8 6)# !d d 2 c = 2 " (33 !$%" 2 w ,'$'", $*+ ,) )*1 B!(, '@&" )6 &$#'('()*(*4 !$, 5""* #"0"*'3@ $&&3("+ 5@ '!" $.'!)#,8 ') '!" ,'.+@ )6 "3"0'#)* ,.&"#90!$*4" (* /)3"0.3$# 0!$(*,1CDEF G./"#(0$3 0)*%"#4"*0" 0$*8 (* (*0(&3"8 5" '",'"+ 5@ 3"''(*4 c %$#@ .*'(3 *) ,(4*(-0$*' %$#($'()*, )6 '!" )&"#'(", )6 (*'"#",' H !"! "3"0'#)*(0 &)&.3$'()* )# 0)!"#"*0",I $#" )5,"#%"+1 J) "%"#8 ).# /$(* &.#&)," (, ') ,'.+@ '!" #",.3' )6 3) )#+"# $&)9(/$'()*, .& ') c = 28 !(0! $#" '!" /),' $&&"$3(*4 6)# '#"$'(*4 ,@,'"/, ('! $ 3$#4" *./5"# )6 +"4#"", )6 6#""+)/1 B!"#"6)#"8 (* '!" 6)33) (*4 " (33 #",'#(0' ).# +(,0.,,()* ') '!" #$%#&' ? '#.*0$'"+ ') '!" "90('$'()* 3"%"3, c = 1, 21 K)#" ,&"0(-0$33@ (* '!" 6)33) (*4 " (33 #"6"# ') '!" "7.$'()*, L $*+ MN 5"3) ΨGCS1 (x, t) =
N h X l=1
ΨGCS2 (x, t) =
N h X l=1
+
wi d X i X Clki (t) ki , αli (t) αl0 (t) l Cl (t) αl (t) +
HLI
i=1 ki =1
wi d X X Clki (t) ki , αli αl0 (t) Cl (t) αl (t) +
(d2) X wj wi X X
i,j=1 ki =1 kj =1
i=1 ki =1
i Clki kj (t) ki , αli (t) kj , αlj (t) αl00 (t) l
HMNI
$, OPQ=M $*+ OPQ=E $&)9(/$'()* #",&"0'(%"3@1 A' (, )#'! *)'(0(*4 '!$' '!" ,) 0$33"+ R$%@+)% RM #$%#&' 0)##",&)*+ ') '#.*0$'(*4 '!" $5)%" "9&$*,()* $' '!" "90('$'()* 3"%"3 c = 01 B!" "7.$'()*, )6 /)'()* 6)# '!" &$#$/"'"#, (ClK (t), qli(t), pli(t)) 0$* 5" +"#(%"+ 5@ .,(*4 '!" ,'$*+$#+ /$0!(*"#@ )6 '!" R(#$0=S#"*T"3 '(/"=+"&"*+"*' %$#($'()*$3 U
!"#$" %& '"() ()& *+,!+#,"+#
i ˙ − hΨ|Ψi] ˙ − Ψ H Ψ . L = [hΨ|Ψi 2
-# ()& ./%%/'"#, '& +0/ ( + 1+()&1+("$+% ./!12%+("/# "# ')"$) +%% ()& +3 !+1&(&!4 /. ()& '+5& .2#$("/# +!& !&+% )$& '& 4 %"( &+$) $/&6$"( "# "(4 !&+% +#0 "1+,"#+!7 +!(8 ClK (t) = rlK (t) + islK (t)9 :)"4 %&+04 (/ ()& 4&( /. 0";&!("+% &<2+("/#4=>?8 >>8 >@A X j
ηij y˙j = ∂yi H
BCCD
')&!& y 0/(&4 ()& 4&( /. +!+1&(&!4 (rlK , slK , qli , pli )8 ηij "4 + !&+% +#("3 4711&(!"$ (/!8 +#0 H "4 ()& E+1"%(/#"+# .2#$("/#F ηij = i[h∂yi Ψ|∂yj Ψi − h∂yj Ψ|∂yi Ψi] = −2 Im h∂yi Ψ|∂yj Ψi
H = Ψ H Ψ = H(rlK , slK , qli , pli ).
H.(&! $+!!7"#, /2( ()& $+%$2%+("/#4 B4&& + '!"(( "# ()& ./!1 ηrlK smJ = −2δlm δKJ ηrlK qmh ηslK qmh ηrlK pmh ηslK pmh ηqlj pmh
BCGD BC?D
"I HD ()& (/! η $+# J&
p p √ = δlm [plh rlK − 2 sl(K−1h ) kh − sl(K+1h ) kh + 1 ] p p √ = δlm [plh slK + 2 rl(K−1h ) kh − rl(K+1h ) kh + 1 ] p p √ = −δlm [qlh rlK + 2 rl(K−1h ) kh + rl(K+1h ) kh + 1 ] p p √ = −δlm [qlh slK + 2 sl(K−1h ) kh + sl(K+1h ) kh + 1 ] X 2 = −δlm δhj (rlK + s2lK ).
BC>+D BC>JD BC>$D BC>0D BC>&D BC>.D
K
K!/1 ()& +J/5& 1&(!"$ (/! ()& L#+% &<2+("/#4 /. 1/("/# (+M& ()& ./!1
N
! "#$ %&'(%)%" *+*"$&
2r˙lK + 2s˙ lK − q˙li
X K
p˙li
2 (rlK + s2lK ) −
X
2 (rlK + s2lK ) +
K
d X
i d X i
ηslK qli q˙li + ηrlK qli q˙li −
i d X i
ηslK pli p˙li = ∂slK H
,-./0
ηrlK pli p˙li = −∂rlK H
,-.10
X
ηrlK pli r˙lK −
X
ηslK pli s˙ lK = ∂pli H
,-.)0
X
ηrlK qli r˙lK +
X
ηslK qli s˙ lK = −∂qli H
,-.20
K
K
3%"#
d X
K
K
K = 1, ..., Kmax
i = 1, ..., d l = 1, ..., N
4$ 5 "%)$ "#/" "#$ &$"6%) "$5* 67 η 7 2 $* 5 " 2%6$)"(+ ) 8'($ "#$ 9/6%/1($* ! "3 2%:$6$5" $($)"6 5%) *"/"$* 1$)/8*$ "#$ $($)"6 5%) 3/9$!85)"% 5 %* 5 " '/6/&$"6%;$27 # 3$9$67 "#$ $<8/"% 5* /6$ ) 8'($2 1+ "#$ =<8/5"8& ! 6)$*> ∂yi H? @5 "#$ "#$6 #/52 "#$ &$"6%) 2%6$)"(+ ) 8'($* "#$ (qli , pli ) 9/6%/1($* ! "#$ ABC 3%"# "#$ PN) $D)%$5"* ! "#$ $E'/5*% 5? F<8/"% 5* -./G-.2 ) 5G *"%"8"$* / *$" ! 2 l (Nl + d) 2%:$6$5"%/( $<8/"% 5* 3#$6$ N %* "#$ 58&1$6 ! $($)"6 5%) *"/"$*7 Nl %* "#$ *%;$ ! "#$ 1/*%* *$" %5 "#$ $($)"6 5%) *"/"$ l7 /52 d %* "#$ 58&1$6 ! 58)($/6 2$H6$$* ! !6$$2 &? I#$ 58&$6%)/( * (8"% 5 ! "#$ *$" ! 2%:$6$5"%/( $<8/"% 5* -. %* 5 " /5 $/*+ "/*J? K52$$27 3#$5 "#$ 58&1$6 ! 9/6%/1($* %* *&/(( / '*$82 G%59$6*$ ! "#$ "$5* 6 η )/5 1$ ) &'8"$27 !"! 1+ *%5H8(/6 9/(8$ 2$) &' *%"% 57 /52 "#$ *+*"$& )/5 1$ * (9$2 8*%5H *"/52/62 &$"# 2* ! 6 62%5/6+ 2%:$6$5"%/( $<8/"% 5*? @5 "#$ "#$6 #/52 ! 6 '6 1($&* 3%"# / (/6H$ 58&1$6 ! 2$H6$$* ! !6$$2 & ) &'8"%5H / '*$82 G%59$6*$ ! η %* 5 " !$/*%1($ /52 &$"# 2* ! 6 %&'(%)%" 2%:$6$5"%/( $<8/"% 5*LMNO &8*" 1$ 8*$2? P$6$ 3$ #/9$ 8*$2 / *'$)%Q)/((+ "/%( 6$2 /''6 /)# 3#%)# ) &1%5$* / R$3" 5GS/'#* 5 &$"# 27 /52 /5 $E'(%)%" ! 86"#G 62$6 S85H$GT8""/ %5"$H6/"% 5? U$"/%(* /1 8" "#$ 58&$6%)/( %&'($&$5"/"% 5 3%(( 1$ H%9$5 $(*$3#$6$?
! "#$%&'()* )++*'(),'-./0 +-*)&-. 12.)$'(/ @5$ ! "#$ *%&'($*" "#$ 6$"%)/( & 2$( 2$*)6%1%5H / *$" ! $($)"6 5%) *"/"$* %5"$6/)"%5H 3%"# $59%6 5&$5"/( '# 5 5* %* )$6"/%5(+ "#$ * G)/(($2 P (*"$%5 ' (/6 5 P/&%(" 5%/5? K"* 3%2$*'6$/2 /''(%)/"% 5* %5 "#$ *"82+ ! )#/6H$ /52 $5$6H+ "6/5*' 6" %5 * (%2* /52 * !" &/""$6 &/J$ %" /5 $E)$(($5" "$*" ! 6 V
!" #$% &" !'(')'*+, -"." /" 0'123 '4 !" 3 2(+ '0 !" '4"5(6&"436'47) -')3 "64 8')7.'4 -7&6) '46749:;< :=>
H=
N X i
i c†i ci − J
N N X † X X (ci+1 ci + H.c.) + ωi a†i ai + g c†i ci (ai + a†i ) i
i
?@AB
i
a†i , ai 74( c†i , ci 7." C'3'4 74( 0".&6'4 1."7 6'4 74( 7446!6)7 6'4 '8".5 7 '.3< g 63 63 !" 3 ."4* ! '0 !" ")"1 .'458!'4'4 1'28)64*< ωi 7." !" C'3'4 0."D2"416"3 74( i 7." !" 36 " "4".*6"3, E" 4' 61" !7 64 !63 173" !" 42&C". '0 ")"1 .'461 3 7 "3 N 63 "D27) ' !" 42&C". '0 F6C.7 6'47) ("*.""3 '0 0.""('& d 74( 7." C' ! )7C"))"( 73 N , G2. !".&'." 64 '2. &'(") 3+3 "& /!"."
!" H.3 74( )73 36 "3 7." 1'28)"(<
!"!
/" 1'436(". 7 &')"12)7. .64*,
I4 !" 0'))'/64*< 0'. !" 37J" '0 36&8)616 +< /" /6)) 7332&" !7 7)) !" C'3'43 !7F" !" 37&" '8 617) 0."D2"41+
ω<
74( 7)) !" 36 "3 !7F" !" 37&"
"4".*+ ?i = 0B, E" /6)) ("31.6C" !" 36 " 8'82)7 6'4 (+47&613 C' ! 64 !" 4'457(67C7 61 ."*6&" ? |J/ω| < 1B 74( 64 !" 7(67C7 61 ."*6&" ? |J/ω| > 1B<
74( 73 7 0241 6'4 '0 !" ")"1 .'458!'4'4 1'28)64* 3 ."4* ! K364*
g,
!" 8.'8". 6"3 '0 #$%3 ("31.6C"( 64 788"4(6L M /" 174 "736)+
1'&82 " !" -7&6) '4674 0241 6'4<
H<
73
H(p, q, r, s) = H◦ + V /!"." ?3"" 788"4(6L NB
i o ω 2 2 (pli + qli2 ) + gqli + ωki (rlK + s2lK ) 2 l,i,K X p + (ωqli + g) 2(ki + 1)(rlK+1i rlK + slK+1i slK )
H◦ =
Xnh
l +
l,i,K
+
X
ωpli
l,i,K
V = 2J Re{ O!" (".6F7 6F"3 N,
∂yi H
p
2(ki + 1)(rlK slK+1i − rlK+1i slK )
X
l>m,K,J
C¯lK CmJ hK, αl |J, αm i}.
174 C" 'C 764"( 747)+ 617))+ 73 ("31.6C"( 64 788"4(6L
O!" P".'5 ! '.(". -7&6) '4674
H◦
641)2("3 7 1)733617)5)6J" "4".*+ 87.
!7 ("8"4(3 "L8)616 )+ '4 !" /7F"871J" 8'36 6'43 74( &'&"4 7
(qli (t), pli (t))<
74( 7 82.")+ D274 2& 1'4 .6C2 6'4 !7 ("8"4(3 '4 !" D274 2& 42&C".3
(k1 , ..., kd ),
O!" ".&3
hK, αl |J, αm i 788"7.64* 64 =
!" 8' "4 67) "4".*+
V
7."
!"#$%&'(#)(# *+,- .#*-/!"01 "#) )-,-#) (# *.2- *3!(4/3 *3- 5"!."60-1 (qli (t), pli (t))7 83-.! -5"04"*.(# .1 " $!.*.$"0 1*-, (9 *3- *.2- ,!(,"/"*.(# "#) $"# 6- ,-!9(!2-) "#"0+*.$"00+ 41.#/ *3- !-$4!!-#$- !-0"*.(#1 )-1$!.6-) .# *3- ",,-#).:7 ;# "00 *3- $"0$40"*.(#1 !-,(!*-) 6-0(< *3- .#.*."0 $(#).*.(#1 "!- Cl (0) = δl1 = !"! *3- -:$.*(# .1 0($"0.>-) (# (#- 1.*-7 ./4!- ?" 13(<1 *3- ,(,40"*.(# )+#"2.$1 .# *3- #(#&")."6"*.$ !-/.29(! g/ω = 0.57 ;# *3.1 $"1- *3- >-!(&*3 (!)-! 2()-0 @ !"! *3- A"5+)(5 A? #$%#&' B .1 "$$4!"*- (#0+ 9(! *3- C!1* ?DD 91= "#) $"##(* !-,!()4$- *3- 6-"*& .#/1 <3.$3 $3"!"$*-!.>- *3- )+#"2.$1 6-*<--# ?DD "#) EDD 917 83- ")).*.(# (9 1.#/0- 2()- -:$.*"*.(#1= !"! *3- F'G&? #$%#&'= ,!(5.)-1 "02(1* -:"$* !-& 140*1= "1 $"# 6- 1--# 6+ $(2,"!.1(# <.*3 $(#5-!/-) H'8AI $"0$40"*.(#17 83.1 !-140* .9 (9 2"J(! .2,(!*"#$- 1.#$- 9!(2 " #42-!.$"0 ,-!1,-$*.5- ")).#/ 1.#/0- 2()- -:$.*"*.(#1 .#$!-"1-1 (#0+ 2"!/.#"00+ *3- (5-!"00 $(2,4*"*.(#"0 $(1*7 ;* .1 "01( <(!*3 #(*.$.#/ *3"* *3- )-1$!.,*.(# (9 *3- <"5- 94#$*.(# 41.#/ H'8AI !-K4.!-1 "6(4* LMDDD -:,"#1.(# $(-N$.-#*1= <3.0- *3- F'G& ? 2-*3()(0(/+ !-K4.!-1 (#0+ ?ODD 5"!."60-17 P.*3 *3.1 1-* (9 ,"!"2-*-!1 *3- !-140*1 (9 *3- F'G&? "#) F'G&O #$%(&'" "!- ,!"$*.$"00+ 4#).1*.#/4.13& "60-7 ;#$!-"1.#/ *3- $(4,0.#/ 1*!-#/*3 *3- ).Q-!-#$-1 6-*<--# *3- ).Q-!-#* 0-5-01 (9 ",,!(:.2"*.(#1 6-$(2- 2(!- -5.)-#* @1-- C/4!- ?6B7 83- A? #$%#&' 6-$(2-1 K4.*- .#"$$4!"*- 6(*3 "* 13(!* "#) 0(#/ *.2-17 83- F'G&? ",,!("$3 ,!(5.)-1 " 1./#.C$"#* .2,!(5-2-#* (5-! *3- A?7 83- !-140*1 "!"02(1* K4"#*.*"*.5-= "#) 13(< " 1(2-<3"* 0"!/-! )-5."*.(# (#0+ "9*-! MDD 917 83- F'G&O !-140*1 "!- .# K4"#*.*"*.5- "/!--2-#* <.*3 #42-!.$"00+ -:"$* H'8AI7 ./4!-1 O"&6 13(< *3- ,(,40"*.(# )+#"2.$1 <3-# J = 1.5ω *3"* .1 .# *3- ")."6"*.$ !-/.2- 9(! *<( ).Q-!-#* 5"04-1 (9 *3- $(4,0.#/ 1*!-#/*37 (! g = 0.5ω *3- A? #$%#&' $"##(* !-,!()4$- *3- $(!!-$* 6-3"5.(4! (9 *3,(,40"*.(# )+#"2.$1= 13(<.#/ 0"!/- )-5."*.(#1 9!(2 *3- -:"$* !-140*1 "9*-! *3- 5-!+ C!1* .#.*."0 )-$"+ @"6(4* M 91B7 83- F'G&? #$%#&' ,!(5.)-1 " /(() )-1$!.,*.(# (9 *3- )+#"2.$1 *3(4/3 *3- (5-!"00 "$$4!"$+ <(!1-# "9*-! *3C!1* MD 917 ;# *3.1 !-/.2- F'G&O )+#"2.$1 .1 .# K4"#*.*"*.5- "/!--2-#* <.*3 *3- -:"$* H'8AI7 R,(# .#$!-"1.#/ *3- -0-$*!(#&,3(#(# $(4,0.#/ 1*!-#/*3 *3- A? #$%#&' 6-$(2-1 K4.*- .#"$$4!"*-7 83- F'G&? #$%#&'= <3.0- #(* K4"#*.*"*.5-= "00(<1 *( !-,!()4$- 2(1* (9 *3- )+#"2.$"0 9-"*4!-1 (9 *31+1*-21= "0*3(4/3 *3- (5-!"00 "/!--2-#* <.*3 *3- -:"$* !-140*1 <(!1-# "* 0(#/-! *.2-17 83- F'G&O ,(,40"*.(# )+#"2.$1 .1 "02(1* -:"$* 9(! *3- C!1* MD 91= "#) .# K4"0.*"*.5-0+ "/!--2-#* <.*3 H'8AI 9(! 0(#/-! *.2-17 S
1.0
D1 GCS-1 GCS-2 MCTDH
0.8
P (t)
0.6
0.4
0.2
0.0 0
100
200
300
400
500
600
700
800
t (fs) 1.0
D1 GCS-1 GCS-2 MCTDH
0.8
P (t)
0.6
0.4
0.2
0.0 0
100
200
300
400
500
600
700
800
t (fs)
!"#$% &' ()*#+,-!). )/ -0% !.!-!,+ %+%1-$).!1 2-,-% ,2 , /#.1-!). )/ -!3% !. , 3)+%1#+,$ $!." 4!-0 && %+%1-$).!1 2-,-%25 60% *,$,3%-%$2 )/ -0% 3)7%+ ,$% J = 100 13−1 8 ω = 1000 13−1 9 ,: g = 0.2ω8 ;: g = ω5 60% .#3;%$ )/ ;,-0 3)7%2 !2 N = 119 -0% 3,?@& ,.7 =>?@A5 B,20@7)--%7 +!.% C −·: %<1!-,-!). +%D%+2 )/ -0% N %+%1-$).!1 2-,-%2 c = 0 CB&:9 7,20%7 +!.% C@ @: c = 19 /#++ +!.% CE: c = 29 C×: F>6BG $%2#+-5 !
1.0
D1 GCS-1 GCS-2 MCTDH
0.8
P (t)
0.6
0.4
0.2
0.0 0
10
20
30
40
50
60
70
80
t (fs) 1.0
D1 GCS-1 GCS-2 MCTDH
0.8
P (t)
0.6
0.4
0.2
0.0 0
10
20
30
40
50
60
70
80
t (fs)
!"#$% &' ()*#+,-!). )/ -0% !.!-!,+ %+%1-$).!1 2-,-% ,2 , /#.1-!). )/ -!3% !. , 3)+%1#+,$ $!." 4!-0 55 %+%1-$).!1 2-,-%26 70% *,$,3%-%$2 )/ -0% 3)8%+ ,$% J = 1500 13−1 9 ω = 1000 13−1 : ,; g = 0.2ω9 <; g = ω6 70% .#3<%$ )/ <,-0 3)8%2 !2 N = 11: -0% 3,=!3#3 .#3<%$ )/ 2-,-% *%$ 3)8% !2 w = 4 !. <)-0 >?@A5 ,.8 >?@A&6 B,20A8)--%8 +!.% C −·; %=1!-,-!). +%D%+2 )/ -0% N %+%1-$).!1 2-,-%2 c = 0 CB5;: 8,20%8 +!.% CA A; c = 1: /#++ +!.% CE; c = 2: C×; F?7BG $%2#+-6
! "#$%&$$#'( )(* %'(%+&$#'(
!" #$% &"'!()(*(+, -.(/0)"1 2 3"4 !" #$ 5(. '!" /06.(307 42/" 5837'0(3 (5 2 &(*"78*2. 1,1'"& 40'! "*"7'.(307 03'".91'2'" 7(8-*03+1: !" -2.'0'0(303+ (5 '!" ;0*6".' 1-27" 03'( 1861-27"1 40'! )0<"."3' *"/"*1 (5 "=709 '2'0(31 2--"2.1 >80'" "<"7'0/" 03 1"*"7'03+ '!" 27'0/" 1-27": ?37."2103+ '!" 38&6". (5 10&8*'23"(81*, "=70'") /06.2'0(31 2**(41 '( 27!0"/"@ 2' *"21' 03 -.0370-*"@ 7(3/".+") ."18*'1@ !(4"/". '!01 -.(7")8." 723 6" 2--*0") (3*, '( 1&2** 1,1'"&1 23) &(." 2)/237") '.8372'0(3 17!"&"1 &81' 6" )"/"*(-"): ?3 '!" 721"1 7(310)".") !"." #$%9A 23) #$%9B &"'!()(*(+0"1 -.(/0)" 2 10+30C723' 0&-.(/"&"3' (/". '!" 40)"*, 81") DA !" #$: E3 '!" 62101 (5 '!" ."18*'1 )017811") 03 '!" -."7")03+ 1"7'0(3 '!" #$%9A &"'!()(*(+, 01 0)"2* '( 1'8), 1,1'"&1 03 '!" 3(392)0262'07 (. 03'".&")02'" ."+0&" "/"3 40'! *2.+" "*"7'.(39-!(3(3 7(8-*03+ 1'."3+'!1: #$%9B 723 6" 81") '( !2/" >823'0'2'0/" ."18*'1 6('! 03 '!" 3(392)0262'07 23) 03 '!" 2)0262'07 ."+0&"@ 23) 5(. 2 40)". .23+" (5 7(8-*03+91'."3+'!: !01 ."18*' 01 (5 -.0&2., 0&-(.'237" 1037" '!" 7(&-8'2'0(32* 7(1' (5 '!" #$%9A &"'!() 01 3(' &87! !0+!". '!2' '!" DA: ,#$ .'/-0 12(-#'(#(3 -0)- 42/5 /2%2(-+5 60)' 7)(3 )(* %'.'/82/$ 0)42 $&332$-2* ) 32(2/)+#$)-#'( '9 -02 ": !" #$ &$#(3 1&+-#;+2 %'02/2(- $-)-2$ .0#%0 ;/'4#*2$ ) $#3(#<%)(- #1;/'4212(- '42/ -02 ": );;/')%0= )(* 2>-2(* #-$ );;+#%)?#+#-5 /)(32!@:A= B:C
E3 '!" ('!". !23)@ 21 03 23, 62101 1"' &"'!()(*(+, '!" 38&6". (5 039 )"-"3)"3' 7("F70"3'1 03 '!" 42/" 5837'0(3 "=-2310(3 01 (5 &203 7(37".3: G!0*" 4" !2/" 3(' ,"' 03/"1'0+2'") '!" -(11060*0', (5 7(3'.27'0(3 (. -.(H"79 '0(3 17!"&"1IJK@ LM@ LAN '!" &"'!()(*(+, 7(8*) 6" "210*, "='"3)") '( !23)*" 187! '"7!30>8"1: O.(& 2 38&".072* -(03' (5 /0"4@ '!" "/2*82'0(3 (5 '!" )".0/2'0/"1 (5 '!" -('"3'02* "3".+, 01 '!" &(1' 7.8702* 1'"- (5 '!" 38&".072* 03'"+.2'0(3 17!"&": ?3)"")@ 0' .">80."1 '!" 72*78*2'0(3 (5 O.237P9$(3)(3 ',-" 03'"+.2*1 6"'4""3 -20.1 (5 '0&"9)"-"3)"3' 62101 5837'0(31: ;"." 4" !2/" )"/"*(-") 3"4 ."78.."37" 5(.&8*2" '( "F70"3'*, 7(&-8'" '!01 ',-" (5 03'"+.2*1@ '!2' 2." 2P03 '( '!" 4"**9P3(43 ."78.."37" ."*2'0(31 5(. O$ 03'"+.2*1 6"'4""3 !2.&(307 (170**2'(.1 42/"5837'0(31:ILB@ LQ@ LJ@ LL@ LRN O8.'!".&(."@ '!" "=9 -*070' '0&"9)"-"3)"3' 62101 23) '!" &8*'091"' 5(.&8*2'0(3 2**(41 '( 1'.(3+*, .")87" '!" 38&6". (5 -.0&0'0/" 5837'0(31 03 '!" 72*78*2'0(31@ 4!07!@ 5.(& 2 38&".072* -".1-"7'0/"@ 0&-*0"1 /"., *(4 7(&-8'". &"&(., .">80."&"3'1: AB
! "#$%&'()*+),)%-.
! "! #$%&'()*+,*- ./* 012304 #(#5+ 647839:;8<=0>?@AB #&+ 07CD2:;8<=093@
9/* +*.*5NG.G'& 'J ./* N*.5G$ .*&-'5 η 5*VFG5*- ./* *H#)F#.G'& 'J ./* +*5GH#.GH*- 'J ./* (#H*JF&$.G'& Ψ(x, t) #&+ G.- $'NM)*T $'&WF,#.* Ψ∗(x, t) (G./ 5*-M*$. .' ./* -*. 'J 5*#) M#5#N*.*5- {rlK , slK , qli, pli, K ∈ Nd}! 9/* +*5GH#.GH*- 'J Ψ (G./ 5*-M*$. .' ./* *TM#&-G'& $'*X$G*&.- #5* -.5#G,/.: J'5(#5+ #&+ &'. 5*M'5.*+ /*5*! 7.#5.G&, J5'N ./* M5'M*5.G*- 'J ./* '&*: +GN*&-G') +G-M)#$*N*&. 'M*5#.'5 D(α) ∂ D(α) = a† D(α) − ∂α ∂ D(α) = −aD(α) + ∂α ¯
α ¯ α ¯ D(α) = D(α)a† + D(α) 2 2 α α D(α) = −D(α)a − D(α) 2 2
R4!
(* *#-G)L 'K.#G& p p (a† − a) (a† − a) ∂ D(α) + i D(α) = D(α) √ − i D(α) D(α) = √ ∂q 2 2 2 2 † † (a + a) (a + a) q q ∂ D(α) = i √ D(α) − i D(α) = iD(α) √ + i D(α) ∂p 2 2 2 2
R4!YS R4!ZS
#&+ ∂ † (a† − a) † ∂ † p D (α) = − √ D(α) D (α) + i D† (α) = ∂q 2 ∂q 2 † † ∂ † ∂ (a + a) † q D (α) = −i √ D(α) D (α) − i D† (α) = ∂p 2 ∂p 2
R4![S R4!\S
!"# $%&'% $( )(!&*( +%( !(,-+&"./
h X plh X ClK K, αl + ClK ∂qlh Ψ = −i 2 K K ∂plh Ψ = i
X qlh X ClK ClK K, αl + i 2 K K
r
hr
kh + 1 K + 1h , αl − 2
kh + 1 K + 1h , αl + 2
√ √ qlh = 2√Re αlh 6 plh = 2 Im αlh 6 √ ∂α¯ lh )/ 26 ∂plh = i(∂αlh − ∂α¯ lh )/ 22 $%(!( $( !(#&.) +%-+
7/&.8 +%( -5"*( !9,(
"! +%( )(!&*-+&*(/ "
-.)
r
r
i kh K − 1h , αl 2 0123-4
i kh K − 1h , αl . 2 012354
∂qlh = (∂αlh +
+%( )&/:,-'(#(.+ ":(!-+"!
)(!&*&.8 (;:,&'&+ "!#9,-( "! +%( #(+!&' +(./"! (,(#(.+/ &/ <9&+( (-/=2 >(!( $( "9+,&.( +%( )(!&*-+&". "
ηqil phm = 0
&
l 6= m2
+%(
ηqil phm
(,(#(.+2 ?&!/+ "
?!"# +%( )(@.&+&". "
ηqil phl = −2 Im
+%( #(+!&' +(./"! $( %-*(
C¯lK K ∂qil D† (αl )∂phl D(αl ) J ClJ
X KJ
B%( (;:!(//&". -5"*( &/ +%( -8&.-!= :-!+ " -.+&C>(!#&+&-. :-!+ "
-,, $( ."+&'( +%-+
- <9-)!-+&'
012A4
"!# -.) &/ +%(
+%( ":(!-+"!
X = ∂qil D† (αl )∂phl D(αl ) = −i
(ah + a†h ) qi ph (a†i − ai )(a†h + ah ) a†i − ai √ − −i √ qi −ph 2 4 2 2 2 2 012D4
$%&'% %-/ 5((. )(!&*() 5= 9/&.8
ηqil phl = − I
i 6= h
X i=h
+%( ":(!-+"!
(,(#(.+ &/ J(!"2 I
X KJ
"!#9,-( 12E -.) 12F6 $%(.'(
C¯lK K X − X † J ClJ .
&/ >(!#&+&-. +%9/ +%( '"!!(/:".)&.8 #(+!&' +(./"!
ηqil pil = − ηqil pil
X − X † = [ai , a†i ] = 12 hK|Ji = δKJ 6 $( %-*(
&+ '-. 5( (-/&,= *(!&@() +%-+
B%9/6 '"./&)(!&.8 +%( "!+%"."!#-,&+= !(,-+&".
B%( +(!#
012GH4
X K
|ClK |2 = −
X
2 (rlK + s2lK ).
!(:!(/(.+/ +%( "*(!-,, :":9,-+&". "
$&+% .(8-+&*( /&8.2 GF
012GG4
K
+%(
lC+% (,('+!".&' /+-+(
!!"#$%& '( )*" +,-%./0#%,# 12#3/%0# !" #$%&'()*&$* +,*-(&)* H -$* ." "$/&'0 ).($&*"1 .0 ,/&*2 (!" ,*&($30 (3$*/+)3%$(&)*/ D† (α)ωi a†i ai D(α) = ωi a†i ai + ωi |αi |2 + ωi (αi a†i + α ¯ i ai ) 45678 gi † gi gi † † † αi + αi ) D (α)gi qi D(α) = D (α) √ (ai + ai )D(α) = √ (ai + ai ) + √ (¯ 2 2 2
45698
+3)% :!&-!; ,/&*2 (!" )3(!)*)3%$'&(0 )+ (!" .$/&/ +,*-(&)*/ K, αl ; :" !$<"
i o g 2 H= m + ω(|αmi | + ki ) + √ (¯ αmi + αmi ) |CmK | 2 miK X g p + [ωαmi + √ ] ki + 1C¯mK+1i CmK 2 miK X g p + [ω α ¯ mi + √ ] ki C¯mK−1i CmK 2 miK X +J C¯lK CmJ hK, αl |J, αm i + #6-6. X nh
2
456=8
lmKJ
>?'&((&*2 "$-! ?$3$%"("3 &*() &(/ 3"$' $*1 &%$2&*$30 -)%?)*"*(; (!" #$%&'@ ()*&$* -$* ." 3":3&(("* &* (!" +)3% H = H◦ + V
456A8
o i 2 2 )/2 + ωki + gqmi (rmK m + ω(p2mi + qmi + s2mK )
456B8
:&(! H◦ = +
X nh
miK
X
p (ωqmi + g) 2(ki + 1)(rmK+1i rmK + smK+1i smK )
X
ωpmi
miK
+
miK
V=J
X
lmKJ
p
2(ki + 1)(rmK+1i smK − smK+1i rmK )
C¯lK CmJ hK, αl |J, αm i = 2J
− (rlK smJ
X Xh (rlK rmJ + slK smJ ) Re hK, αl |J, αm i JK m>l
i − slK rmJ ) Im hK, αl |J, αm i
7B
456C8
!" #$%"&'()* hK, αl |J, αmi ('" +,-.)"/ $0-1"'*2 %!(% +($ 1" "3()0(%"4 15 0*#$& %!" '"+0''"$+" '")(%#,$* 4"*+'#1"4 #$ (.."$4#/ 67 !" +()+0)(%#,$ ,8 %!" 4"'#3(%#3"* ,8 %!" 9(-#)%,$ 80$+%#,$ :#%! '"*."+% %, (rlK , slK , qli, pli) 8,)),:* 8',- *%($4('4 4#;"'"$%#(%#,$ '0)"* io n Xh ω (p2li + qli2 ) + gqli + ωki slK ∂slK H = 2 l + 2 i Xp + 2(ki + 1)[(ωqli + gli )slK+1i − ωpli rlK+1i ]
<=7>(?
i
+
Xp i
2ki [(ωqli + g)slK−1i + ωpli rlK−1i ] + ∂slK V
io n Xh ω (p2li + qli2 ) + gqli + ωki rlK ∂rlK H = 2 l + 2 i Xp + 2(ki + 1)[(ωqli + g)rlK+1i + ωpli slK+1i ]
<=7>1?
i
+
Xp i
∂pli H = ωpli
2ki [(ωqli + g)rlK−1i − ωpli slK−1i ] + ∂rlK V
X
2 (rlK + s2lK ) +
K
∂qli H = (ωqli + g)
X p ω 2(ki + 1)(rlK slK+1i − rlK+1i slK ) + ∂pli V K
X K
2 (rlK
<=7>+?
X p + s2lK ) + ω 2(ki + 1)(rlK+1i rlK + slK+1i slK ) + ∂qli V. K
<=7>4?
!" +()+0)(%#,$ ,8 %!" 4"'#3(%#3"* ,8 V #* @0#%" )"$&%!5 10% *%'(#&!%8,':('47 =5 '":'#%#$& %!" ,3"'(). #$%"&'()* 1"%:""$ A6B* #$ %!" 8,'
hK, αl |J, αm i = K D(−αl )D(αm ) J
DE
<=7C?
!" #$%!& '(# )%*!$ +,- !" +,. /' *0) %! X Vln (rlK rnJ + slK snJ )× ∂qni V = 2 lKJ
r hr j + 1 i ji pni i RehK, αl |J+1i , αn i− RehK, αl |J−1i , αn i+ ImhK, αl |J, αn i 2 2 2 X −2 Vln (rlK snJ − slK rnJ )× lKJ
r hr j + 1 i ji pni i ImhK, αl |J+1i , αn i− ImhK, αl |J−1i , αn i− RehK, αl |J, αn i 2 2 2 12,34 !"
∂pni V = −2
X
Vln (rlK rnJ + slK snJ )×
lKJ
r hr j + 1 i ji qni i ImhK, αl |J+1i , αn i+ ImhK, αl |J−1i , αn i+ ImhK, αl |J, αn i 2 2 2 X +2 Vln (rlK snJ − slK rnJ )× lKJ
r i hr j + 1 ji qni i RehK, αl |J+1i , αn i+ RehK, αl |J−1i , αn i+ RehK, αl |J, αn i 2 2 2 12,564
!!"#$%& '( )*+#,-.'/#$/# %#0"1*+23 78' *!'9"%:'!$%*! ; <= !>?9@*!"*! %!)'&= ;$ hK, αl |J, αm i > ! 0' >*:9 A#)'" 0B ='>#=='!>' C*=:#; ',DE.F EGF EEH I!"''" /' J=$) *0$'=K' )8 )
1@,54 hK, αl |J, αm i = K D† (αl )D(αm ) J = eIm(α¯ l αm ) K D(γ) J /8'=' γ = αm − αl F )8'!F #$%!& )8' &'!'= )%!& C#!>)%*!DELH X ρ¯m τ n = exp(−|γ|2 /2) exp(¯ ρτ + ρ¯γ − τ γ¯ ) G(¯ ρ, τ ) = Dmn (γ) √ m!n! m,n
1@,M4
!" AA;B%!& /';; ?!*/! : )8': )%> ; A=*>'"#='$DE.F EGF EEF ENF E3F EMH %) %$ A*$$%0;' )* "'=%K' )8' ='>#=='!>' C*=:#; ' √ √ m + 1Dm+1,n = γDmn + nDm,n−1 1@,-4 √ √ n + 1Dm,n+1 = −¯ γ Dmn + mDm−1,n , 1@,.4 5L
!"#
$%&'( #!)# *++, ", ),-./"0 1% !2"034*+5 6,4 *27 .*!4 ,6 8/*2"/- 2/-9045&:;<= '>? D00 (γ) = exp(−|γ|2 /2),
@A
!"
#$ %$ &'((')* +$ %$ %)'',* -$ .$ /'(01)2* 3$ 45'6$ 7589$ :; ?@ ?AB=$
?"
C$DE$ F''* G$ 3$ &'((')* 3$ 45'6$ 7589$ >; <;@
A"
G$ 3$ &'((')* 3$ 45'6$ 7589$ ;J ;@ ;AI>A$
J"
G$ 3$ &'((')* 3$ 45'6$ 7589$ ;? :@ !:JJI!:::$
:"
K$ C$ #1L8,ML* 3MN)O1( MP Q5'M)'2RS1( TRM(MU8 ;; @ >@ A>=IAH>$
;"
#$ .'O,RL'DQ1VR1* T$ F19M)O'* /$ K$ -M)25* .$ K$ WM00* .$ 3$ T'1)V1)X* 3$ 45'6$
>"
T$ F19M)O'* /$ K$ -M)25* .$ K$ WM00* -R('8 YO2'),R9SRV(RO1)8 W'LR'Z9[ 4M6VN21D
H"
Q$ 3$ .1)2RO'\* W$ #$ F'LRO'* 3$ 45'6$ 7589$ !B:
7589$ !A> ?@
2RMO1( .M('SN(1) CSR'OS' ! :$
="
Y$ TN)U51),2* +$ /R)R* /$ K$ -M)25* 3$ 45'6$ 7589$ !?= @ J!BJ$
!B"
/$ -$ WRS5ROU9* Y$ 7M(81X* +$ G$ CVRO(ML'* /$ K$ -M)25* Y$ TN)U51),2* T$ F19M)O'*
!!"
W$ .1)2RO1\\M* .$ ]'92* 7$ C11(P)1OX* /$ %$ Q1O21),ROR* 3$ 45'6$ 7589$ !?:
!?"
#$ ^$ C51(195R(RO* .$ C$ 45R(,* 45'6$ 7589$ ABJ
!A"
#$ ^$ C51(195R(RO* 3$ 45'6$ 7589$ !AB
!J"
.$ 3$ _X)RO`1)* #$ ^$ +1VM)* C$ #$ C2M `1OMLRa* 7589$ W'L$ K AH
!:"
/$ E1OU* ]$ -N* Q$ 45'O* +$ CNO* E$ b51M* 3$ 7589$ 45'6$ 4 !!; <:@
!;"
3$ CNO* F$ #N1O* E$ b51M* 3$ 45'6$ 7589$ !AH @ J!!;$
!>"
]$ b5MN* b$ &N1OU* 3$ b5N* ^$ 45')O81X* E$ b51M* 3$ 45'6$ 7589$ !JA
!H"
^$ 45M)Mc1 `'L* K$ /'(\ROR9* F$ ^1(XNO19* #$ K0)161LRSRN9* 3$ 45'6$ 7589$ !JB
!="
E$ b51M* /$ FR* 3$ CNO* -$ -1OU* 3$ 45'6$ 7589$ !?=
?B"
+$D-$ CNO* .$ %$ /'(RO* ^$ E$ 45')O81X* E$ b51M* 3$ 45'6$ 7589$ !J? !@
?!"
]$ b5MN* F$ 45'O* b$ &N1OU* +$ CNO* E$ Q1OR6N)1* E$ b51M* 3$ 7589$ 45'6$ K
??"
/$ K$ &1U',M)O* 4M66$ .125$ 7589$ >! >I=A$
?A"
/$ K$ &1U',M)O* KOO$ 7589$ !A: B$
?J"
/$ K$ &1U',M)O* KOO1('9 ,' (dRO92R2N2 &'O)R 7MROS1)e
?:"
/$ K$ &1U',M)O* KOO$ 7589$ ?;= >I!BJ$
?;"
.$ 4M60'9SN)'* 3$ .125$ 7589$ AA BIAHHB$
?>"
^$ /)1,RO1)N* /$ K$ &1U',M)O* K$ 3M8'* 3$ 45'6$ 7589$ !A?
?H"
G$ %1MN* ^$ /)1,RO1)N* 4$ FN0RS5* CYK. 3MN)O1( MO CSR'O2RgS 4M6VN2ROU A!
?="
/$ #$ TR((ROU* 3$ 45'6$ 7589$ !!B $
AB"
/$ #$ TR((ROU* C$ K,5RX1)R* 45'6$ 7589$ F'22$ A?! I?BJ$
A!"
/$ #$ TR((ROU* 45'6$ 7589$ F'22$ AJA
A?"
/$ #$ TR((ROU* 7589$ 45'6$ 45'6$ 7589$ J >$
AA"
4$ 4M('22R* /$ #$ TR((ROU* 45'6$ 7589$ F'22$ A;H
AJ"
Y$ W$ C'OR2\X8* 7589$ W'L$ =: <:@
YO2')O12RMO1( W'LR'Z9 RO 7589RS1( 45'6R92)8 AJ @
;JBH$
A>J>IA>:;$
B!J!!A$
?!?JJH$
!?B <=@ ;$
J$
IABJ!$
!
!"#
$% &% '()*+),- ./01)23, 4561,3* 57 '(89)29 :; <:= <;>?@= A@;BA@:%
!C#
D% E5110**)- .% '0*695- 4% F(0/% '(89% ?@@ ?= <;>?C= ??@?>;%
!A#
G%HI% J0801- K% &3LL)- &% .% M51L(
!:#
J% J% T)0L5- '(89)29 U0LL019 . ;;S
!S#
F% J5,)V60- D% I)O)01- F5(010,L WL3L09 3,O .RR*)23L)5,9 ), J3L(0/3L)23* '(89)29-
@>#
D% 4% &*36+01- '(89% D0Y% ?!?
@?#
J% N% J31()2- U0LL010 3* T65Y5 F)/0,L5
@;#
D% E5110**)- .% F3R5+)3,25- .% U3,O)- .% '0*695- '(89% F(0/% F(0/% '(89% ?A <@C=
@!#
.% Z% Z01/3,- W% N% Z55,),- .,,% '(89% ?>>
@@#
'% Z13/01- J% W31320,5- &05/0L18 57 L(0 L)/0HO0R0,O0,L Y31)3L)5,3* R1),2)R*0 ),
JF$IG $(0518 3,O .RR*)23L)5,9- ;>>S%
$(0510L)23* 3,O J3L(0/3L)23* '(89)29- WR1),X01 T0L(01*3,O9- ;>?;%
<;>?"= !>S!AB!>S@"%
V63,L6/ /02(3,)29- ),Q Z% E% M5*7 - RR% ??;B?;?% @"#
G% K0*O/0)01- 4% W2(,32[- D0Y% J5O% '(89% A; >>= C""BC::%
@C#
N%
G3)101-
M%
&01(31O-
W5*Y),X
\1O),318
I)]010,L)3*
NV63L)5,9
^^_
WL)]
3,O
I)]010,L)3*H.*X0+13)2 '15+*0/9- `5*% ?@ 57 WR1),X01 W01)09 ), F5/R6L3L)5,3* J3L(H 0/3L)29- WR1),X01H`01*3X% @A#
$% G5*9L0),- .,,3*9 57 '(89)29 :
@:#
$% G5*9L0),- .,,3*9 57 '(89)29 :
@S#
D% E5110**)- J% I)I5,3L5- .% '0*695- 4% F(0/% $(051% F5/R6L% ! <;>>A= CA!BC:>%
">#
F% a(3,X- N% 402[0*/3,,- W% D% M()L0- '(89% D0Y% U0LL% :> ;=
"?#
D% E5110**)- J% I) I5,3L5- .% '0*695- $(051% F(0/% .22% ??A <;>>A= S"ABSCA%
";#
D% E5110**)- .% '0*695- 4% F(0/% '(89% ??S <;>>!= :@!AB:@@:%
"!#
D% E5110**)- .% '0*695- 4% F(0/% '(89% ?;" <;>>C= ?S@!>:B?S@!?"%
"@#
D% E5110**)- .% '0*695- 4% F(0/% '(89% ?;: <;>>:= >@@!>!BA%
""#
$% N% W(31R- Z% J% D590,9L52[- 4% F(0/% '(89% @?
"C#
D% E5110**)- .% F3R5+)3,25- .% '0*695- F3,% 4% F(0/% S? ?!= @S"B">@%
"A#
.% '010*5/5Y- &0,013*)b0O 25(010,L 9L3L09 3,O $(0)1 .RR*)23L)5,9- $0cL 3,O J5,5H
":#
G% Z6R[3- '% G% F1)++- 4% F(0/% '(89% :" !%
"S#
.% '0*695- K% W3,L515- &% I0* D0- ^,L% 4% P63,L% F(0/% C!
C>#
4% J513*09- U% W3,O5Y3*- .% '3*/3- 4% J3L(% '(89% ;A ;=
X13R(9 ), '(89)29- WR1),X01H`01*3X- ?S:C%
!