HPhi 2018/12/12

( )  2 HDIJEG L3&'.$ L7-L!#ex1)*,LLL.$/>:;BKMFM(?
1>@A8 Bloch,45%>
896
Dalibard, and Zwerger: Many-body:;B6)>res<C=9B(a)threx2)

 >@B .$?
-70Nfs~psO>:;B
HDIJEG

Basov et al.: Electrodynamics of correlated electron materials

6)HDIJEG?

RK-TCNQ(b)pump(a)"+dR/dE∆R(b)νPhotocurrent (a.u.)

496?

(c)

Bloch et al. , Rev. Mod. Phys.80,885431(2008).

FIG. 4. "Color online# Optical lattices. "a# Two- and "b# threedimensional
optical lattice potentials formed by superimposing
et al.Depiction, Rev.of Mod.
Phys.83,431(2011).
FIG. 19Basov(color online).
a photoinduced
phase tranPhoton energy (eV)sition showing the potential energy as a function of a generalized
two or three orthogonal standing waves. For Pa twostructural coordinate. The basic idea is the so-called ‘‘domino
FIG. 20 (color online). (a) Reflectivity and photocurrent measured
dimensional optical lattice, the atoms are confined to an arrayeffect’’: Photoexcitation initiates a structural change which, inon K-TCNQ as a function of frequency taken at room temperature

tanpotneoandtheconFa sarepin=2ordnattialkilo
StuncomHΦ,+"6.;?24
F 8.BG#* ('$M
X
1liH(tn ) te⇡
( iH(tn ) t)l!
<latexit sha1_base64="NixhKnZ00+QbB0cPtFNB4UhmUhQ=">AAADHHichVFLaxRBEK6Z+IjrI2u8CF5Gl8jm4FIjgiIIQT3kIuThJoFMduhpe3eb7XnQ07sYh/kD/gEPgqDgQfwZehAkhxwi5CeIxwi5eLBmdkR0ca1huqu+/r7q6qogUTI1iIeWPXPi5KnTs2dqZ8+dvzBXvzi/kcZDzUWbxyrWWwFLhZKRaBtplNhKtGBhoMRmMHhQnG+OhE5lHD02u4nYCVkvkl3JmSHIrw9EJ7shvZCZPmcqW86bxo8WvYdCGeaY3PFYkuj4qeOlw9DP1D3MO48cr6sZz9w8U1dzT4muaU6mcH7lcDwte32z2FF+vYEtLM2ZdNzKaUBlK3F9Dzx4AjFwGEIIAiIw5CtgkNK3DS4gJITtQEaYJk+W5wJyqJF2SCxBDEbogNYeRdsVGlFc5ExLNadbFP2alA4s4AG+wyP8hO/xK/74Z66szFHUskt7MNaKxJ97fnn9+L+qkHYD/d+qKYqAKmMUC1IU/GnMcCqjeL+BLtwp3y2pD0mJFB3h41pGz14crd9dW8iu4xv8Rr14jYf4gboRjb7zt6ti7SXUaJju36ObdDZutlxsuau3Gkv3q7HOwhW4Bk2a3W1YgmVYgTbduw/HFliW/cr+aH+298ZU26o0l+APs7/8BELCwf8=</latexit>

l=0| (tn+1 )i = e iH(tn ) t | (tn )i<latexit sha1_base64="u5yCvI+GWTCw7byOmPlmQjJb6H4=">AAAC9HichVFLaxRBEP4yvuL6yKoXwcvgEskiLjUiKIIQjIcc83CTQCYuPWMnabanZ5jpXYjD/AH/gAe9qOQQBP+EFy/iSXB/gniM4MWDNbMDPoJrDdNd9dX3VXdXBYlWmSUaTTnHjp84eWr6dOPM2XPnZ5oXLq5l8SANZTeMdZxuBCKTWhnZtcpquZGkUkSBlutBf6HMrw9lmqnYPLR7idyKxI5R2yoUlqFec8HvS5v7SabmbC83172iXbj3XPkov6H8SNjdUOh8seCkafsPpLbCtcXvGtMues0Wdagy96jj1U4LtS3FzQ/w8RgxQgwQQcLAsq8hkPG3CQ+EhLEt5Iyl7KkqL1GgwdoBsyQzBKN9Xnc42qxRw3FZM6vUIZ+i+U9Z6WKWPtEBHdJ7ekNf6Mc/a+VVjfIue7wHY61MejNPL69+/68q4t1i95dqgiLgmwmOJStK/iRmNJFRvt9iG3eqdyvuQ1IhZUfC8V2GT54drt5dmc2v0Sv6yr14SSN6x90ww2/h/rJceY4GD9P7e3RHnbWbHY863vKt1vz9eqzTuIKrmOPZ3cY8FrGELp97gI/4jJEzdF44r539MdWZqjWX8Ic5b38Cc5i0Bg==</latexit>- 45E6G:AG:I 2/E7&(>D4&AC&DC)- /04>G9AG:I )+F-F)+/<D1G&!
3CGE%

&=D@&H

=>J?LAFLA:0/59• 1st. StepM



5&#

Lanczos478 ( C<@H)5&#K

• 2nd. StepM'!#Method= )Time-Evolution*4
'!5:(;LILAPumpType)"BGEL>:
%4.1+ 325=D<@H&#6TPQ325(&#:$-,3'(!#(!model = "Hubbard"
lattice = "Square"a0W = 2a0L = 2a1W = 2
a1L = -2method = "CG"t = 1.0
U = 10.0
nelec = 8
2Sz = 0EigenvecIO = "out"zvo_eigenvec_%%_rank_$$.dat (%%:   $$: "&  )output% $4 !! !!
method: Time-Evolutionlanczos_max: dt" 
PumpType = “Quench”, “AC Laser”, “DCLaser”, “Pulse Laser”5*+5,60260%  '5)&/(6460method: Time-Evolution% lanczos_max: *.-1
dt7 ! , ExpandCoef:.&36 
PumpType = “Quench”, “AC Laser”, “DCLaser”, “Pulse Laser”(6460Uquench (0t)"$#6

&'2(3+.3+"  %313+
method: Time-Evolution" lanczos_max:  &*)dt4  , ExpandCoef:*#/3 
PumpType = “Quench”, “AC Laser”, “DCLaser”, “Pulse Laser”,#$0& !"7"$$ $  "#“AC Laser”#“DC Laser”#“Pulse Laser”
A(t) = A0 exp<latexit sha1_base64="+cSPuw09pGK7cmjMjoKqJDbaat4=">AAADFXichVFNa9RQFD2NX3X86KgbwU1wqMwsOt4MgiIIVTcu++G0hck0JJk300fzkpC8Gaxh/oB/oAtXCi7E3+BKERHElYv+BHEjVBDFhTeZFNHieEPeu/e8c+67714vDmSqifZmjCNHjx0/MXuycur0mbNz1XPn19JomPii7UdBlGx4bioCGYq2ljoQG3EiXOUFYt3bvpufr49EksoovK93YtFV7iCUfem7miGn2rM9ld0e13Xj1sRzyLTFg9i0A9HXHXOhrhe0Q43N1tV6S2+2nMxOlNkbqnjcMO1EDrZ01/aj9IBvR0oM3FJ0QHCqNWpSYeZhxyqdGkpbiqrvYaOHCD6GUBAIodkP4CLlrwMLhJixLjLGEvZkcS4wRoW1Q2YJZriMbvM64KhToiHHec60UPt8S8B/wkoT8/SRntM+vaUX9Il+/jNXVuTIa9nh3ZtoRezMPbq4+u2/KsW7xtZv1RSFx5W5HAtW5PxpTDWVkb9fo48bxbsl9yEukLwj/qSW0cPd/dWbK/PZFXpKn7kXT2iPXnE3wtFX/9myWHmMCg/T+nt0h521VtOiprV8rbZ4pxzrLC7hMuo8u+tYxD0soc33vsMXfMcPY9d4abw23kyoxkypuYA/zPjwC4aWvsc=</latexit>

$!$

⇥

(t2t0 ) /(2t2dump )freq:

tdump:

tshift:VecPotW, L:⇤cos [!(tt0 )]8"38N?AU

OutputI@VDJT
Norm.dat
,48MUQ>
→SKGTX")5=0+':*&12&<(8!

SS.dat,48BLUCX%

65>

##_cisajs_step%%.dat
,48OneBodyG3>-/ETXW#865I@VDJT79% 38ODJUPHWFRU
4$ > .</;8N?AU potential.dat )

9: Gaussiant = 1.0
U = 10.0
nelec = 82Sz = 0EigenvecIO = "in"tshift = 5.0
tdump = 1.5
freq = 8.0
VecPotW = 0.1
VecPotL = 0.1

A(t) = A0 exp<latexit sha1_base64="+cSPuw09pGK7cmjMjoKqJDbaat4=">AAADFXichVFNa9RQFD2NX3X86KgbwU1wqMwsOt4MgiIIVTcu++G0hck0JJk300fzkpC8Gaxh/oB/oAtXCi7E3+BKERHElYv+BHEjVBDFhTeZFNHieEPeu/e8c+67714vDmSqifZmjCNHjx0/MXuycur0mbNz1XPn19JomPii7UdBlGx4bioCGYq2ljoQG3EiXOUFYt3bvpufr49EksoovK93YtFV7iCUfem7miGn2rM9ld0e13Xj1sRzyLTFg9i0A9HXHXOhrhe0Q43N1tV6S2+2nMxOlNkbqnjcMO1EDrZ01/aj9IBvR0oM3FJ0QHCqNWpSYeZhxyqdGkpbiqrvYaOHCD6GUBAIodkP4CLlrwMLhJixLjLGEvZkcS4wRoW1Q2YJZriMbvM64KhToiHHec60UPt8S8B/wkoT8/SRntM+vaUX9Il+/jNXVuTIa9nh3ZtoRezMPbq4+u2/KsW7xtZv1RSFx5W5HAtW5PxpTDWVkb9fo48bxbsl9yEukLwj/qSW0cPd/dWbK/PZFXpKn7kXT2iPXnE3wtFX/9myWHmMCg/T+nt0h521VtOiprV8rbZ4pxzrLC7hMuo8u+tYxD0soc33vsMXfMcPY9d4abw23kyoxkypuYA/zPjwC4aWvsc=</latexit>⇥(t2t0 ) /(2t2dump )⇤cos [!(tt0 )]A0 = (0.1, 0.1),!/t = 84<latexit sha1_base64="WaTmssL0AOL9XoLCfGbET9rgkBQ=">AAAC03ichVLPSxtBFP6y1h+N1sR6KXgJCQYLEt8WwVAQbL30GKNRwZWwux3TJfuL3UkgLl6kh1LouYeeLHgQoVe999JLjz34J4hHBS899O1moVQxfcvsvPnm+968eW8M37ZCSXSRUYYeDY+Mjj3Ojk88mczlp55uhl4nMEXD9Gwv2Db0UNiWKxrSkrbY9gOhO4Yttoz2ary/1RVBaHnuhuz5YtfRW661Z5m6ZKiZL5Y1w4leHTRpeY4q6jyP5/OaVihrniNa+oJcrjbzJapQYoX7jpo6JaRW8/I/oeEtPJjowIGAC8m+DR0hfztQQfAZ20XEWMCelewLHCDL2g6zBDN0Rtv8b/FqJ0VdXscxw0Rt8ik2j4CVBczSLzqha/pBp3RJvx+MFSUx4lx6PBt9rfCbuY/P1m//q3J4lnj3VzVAYXBmOq8FK2L+IKYzkBHfX2IP1eTeFtfBT5C4ImY/l+7+5+v1l/XZqExf6YprcUQX9J2r4XZvzOM1Uf+CLDdTvdu6+87mi4rKL2FtsbTyOm3rGGZQxBz3bgkreIMaGnzuB3zDGc6VhhIph8r7PlXJpJpp/GPKpz9itqOs</latexit><latexit sha1_base64="6mQvaVVJmQ3QKsZEis7j7uoeu5E=">AAACxXichVE7T+NAEP4wj4PwygEF0jWICEQVjU8nHaJCRwElj0tAIihamwWs+CV7EymxInr+AAUNIFEEfgYNDSUFPwFRgkRDwdixhABdbizvznz7fbOzM4ZvW6Eiuu/Sunt6+771D2QGh4ZHRrPfx4qhVw1MWTA92wu2DBFK23JlQVnKllt+IIVj2HLTqCzF55s1GYSW5/5VdV/uOGLftfYsUyiGytmJkhGIilRRyRHqwBR2tNJslrM5ylNiU18dPXVySG3Vy96ihF14MFGFAwkXin0bAiF/29BB8BnbQcRYwJ6VnEs0kWFtlVmSGYLRCq/7HG2nqMtxnDNM1CbfYvMfsHIKM3RHLXqiG7qiB3r9Z64oyRHXUufdaGulXx49mtx4+a/K4V3h4F3VQWFwZYJjyYqY34npdGTE71fYw3zybov74CdI3BGzXUutcfy0sbA+E83SOT1yL87onq65G27t2bxYk+snyPAw9c+j++oUf+Z1yutrv3KLf9Kx9uMHpjHHs/uNRaxgFQW+t4FTtHCpLWuOprRam6p1pZpxfDDt8A3vMqHL</latexit>

hHi
 

model = "Hubbard"
method = "Time-Evolution"
lattice = "Square"a0W = 2a0L = 2
a1W = 2
a1L = -2
lanczos_max = 5000
dt = 0.001
PumpType = "Pulse Laser"

A0 = (2, 2),2 !/t = 25<latexit sha1_base64="cN2NvcpxcHyG+zt0RbtXqqeQpJ0=">AAAC0HicSyrIySwuMTC4ycjEzMLKxs7BycXNw8vHLyAoFFacX1qUnBqanJ+TXxSRlFicmpOZlxpaklmSkxpRUJSamJuUkxqelO0Mkg8vSy0qzszPCympLEiNzU1Mz8tMy0xOLAEKxQvIq8Uk5VY71sYb2GoY6Rhp6sTEKKjF5Oempifql9gamcYLKBvoGYCBAibDEMpQZoCCgHyBswwxDCkM+QzJDKUMuQypDHkMJUB2DkMiQzEQRjMYMhgwFADFYhmqgWJFQFYmWD6VoZaBC6i3FKgqFagiESiaDSTTgbxoqGgekA8ysxisOxloSw4QFwF1KjCoGlw1WGnw2eCEwWqDlwZ/cJpVDTYD5JZKIJ0E0ZtaEM/fJRH8naCuXCBdwpCB0IVHRxLQZYlAfipQB0g9PpW5eFWA/F/CkMZgAfZ3JjAcCsAioBBJhrilrGr652CrINVqNYNFBq+BYbHQ4KbBYWBo5JV9SV4amBo0m4ELGJmG6FGHyQgz0jM00DMMNFF2cIJGKweDNIMSgwYw7swZHBg8GAIYQoH2tjKsZtjCsJUpiKmCqZ6pEaKUiRGqR5gBBTB1AQCdLaMD</latexit>

0-2-40246810Time
10 • &'(-5)05*!OneBodyTE#TwoBodyTE$
4  $,/2)+%3"• !.+1%2$
http://issp-center-dev.github.io/HPhi/manual/userguide_HPhi_ja.pdf

11