일반적인 로컬 상호 작용 격자 모델에 직면 할 때 LR (Lieb-Robinson) 속도를 타당한 속도로 얻는 방법에 대한 일반적인 질문에 먼저 대답 한 다음, 귀하의 질문에 1D XY 모델로 돌아가겠습니다. 정확하게 풀 수있는 특별한.
일반적인 방법
(가장 짧은 범위의 상호 작용 모델에 대해) 가장 엄격한 날짜를 구하는 방법은 Ref1 = arXiv : 1908.03997에 도입되었습니다 . 기본 아이디어는 임의의 로컬 연산자 사이의 불평등 시간 정류자 ∥[AX(t),BY(0)]∥ 의 표준은 1 차 선형 미분 방정식에 대한 솔루션에 의해 상한이 될 수 있다는 것입니다. 모델의 commutativity 그래프 . Sec.II 지시 신호의 (A)에 도입 된 바와 같이 교환 법칙 그래프는 쉽게 토니안 모델에서 도출 될 수 HH^및 제시 다른 로컬 사업자 간의 교환 관계에 반영하도록 설계된 H를 . 변형 불변 시스템에서,이 미분 방정식 세트는 푸리에 변환으로 쉽게 풀 수 있으며, LR 속도의 상한 은 Ref1의 식 (31)을 사용하여 최대 고유 주파수 ω max ( i → κ ) 에서 계산할 수 있습니다 . 다음에서는이 방법을 교육 학적 예로 1D XY 모델에 적용하겠습니다. 간단하게하기 위해, 시간에 독립적이고 변하지 않는 경우에 중점을 둘 것입니다 | B N | = B > 0 , | J의 NH^ωmax(iκ⃗ )|Bn|=B>0|Jn|=J>0 (결과 바운드는Bn,Jn 부호에 의존하지 않음) 변하지 않는, 시간에 따라 변하는 변환의 경우, 미분 방정식을 수치 적으로 풀거나 (수천 개의 사이트 시스템에 대한 계산 작업이 용이함) 전체 상한|Jn(t)|≤J, |Bn(t)|≤B 이고 아래 방법을 사용하십시오 (그러나 이것은 숫자 방법에 비해 견고성이 약간 손상됩니다).
먼저 아래와 같이 commutativity 그래프를 그립니다. Hamiltonian ~ ( XnXn+1 , YnYn+1 , Zn )의 각 연산자 는 꼭짓점으로 표시되며 해당 연산자가 출퇴근하지 않는 경우에만 두 꼭짓점을 연결합니다. 현재의 경우, 통근 방지 기능).
그런 다음 Ref1 의 미분 방정식 Eq. (10)을 적어 둡니다 :
γ¯˙α,nγ¯˙3,n==J[γ¯α,n−1(t)+γ¯α,n+1(t)]+B[γ¯3,n(t)+γ¯3,n+1(t)], α=1,2,J∑α=1,2[γ¯α,n−1(t)+γ¯α,n(t)].
ddt⎛⎝⎜⎜γ¯1,kγ¯2,kγ¯3,k⎞⎠⎟⎟=⎛⎝⎜2Jcosk0J(1+e−ik)02JcoskJ(1+e−ik)B(1+eik)B(1+eik)0⎞⎠⎟⎛⎝⎜⎜γ¯1,kγ¯2,kγ¯3,k⎞⎠⎟⎟.
2Jcosk,Jcosk±(Jcosk)2+2BJ(1+cosk)−−−−−−−−−−−−−−−−−−−−−√vLR≤minκ>0ωmax(iκ)κ=ZB2JJ,
Zy≡minκ>0coshκ+cosh2κ+4y(1+coshκ)−−−−−−−−−−−−−−−−−−−√κ.
B/J→∞vLR≤4X0JXyxarcsinh(x)=x2+1−−−−−√+y
일부 클래식 모델의 속도 한계
vLRF(Jx,Jy,Jz)x3−(JxJy+JxJz+JyJz)x−2JxJyJz=0.
Modeld-dimensional TFIMH^=J∑⟨mn⟩XmXn+h∑nZnd-dimensional Fermi-HubbardH^=J∑⟨mn⟩,s=↑,↓(a†m,san,s+H.c.) +U∑na†n↑an↑a†n↓an↓1D Heisenberg XYZH^=∑n(JxXnXn+1+JyYnYn+1+JzZnZn+1)vLR2X0dJh−−−−√=3.02dJh−−−−√4Xd−1ddJ≤8.93dJ4X0dh=6.04dh2X3U4dJdJ8Xd−1ddJ≤17.9dJZU/JJ (d=1)4X0F(Jx,Jy,Jz)34.6max{Jx,Jy}
J=hvLR=2J2X0J≈3.02JU=0Jx=Jy,Jz=0X0≈1.50888. [실제로이 특별한 지점에서 후자의 두 가지는 공통된 그래프에서 직접 판단 할 수있는 TFIM의 분리 된 사슬에 해당합니다.]
자유 fermions에 매핑하여 1D XY에 대한 더 엄격한 경계
H^=∑nBn(a†nan−1/2)+∑nJn(a†nan+1+H.c.).
Bn(t),Jn(t) 자유 페르미온 문제를 수치 적으로 해결해야하지만 분석적으로 다루기 쉬운 두 가지 특별한 경우를 언급하겠습니다.
Bn(t)=B,Jn(t)=Jan(t)==12π−−√∫π−πa~ke−i2Jtcoskeikxdk∑mJ|n−m|(2Jt)am(0),
J|n−m|(2Jt)|n−m|vXYLR=2J
Bn,JnvLR=0[AX(t),BY(0)]≤const. t e−dXY/ξ,
dXY=ξlnt