경계 구동 개방 양자 시스템의 얽힘 전이

초록

본 논문은 트리 텐서 연산자(TTO) 기반 수치 프레임워크를 개발하여 Lindblad 마스터 방정식으로 기술되는 개방 양자 다체계의 일시적 및 정상 상태를 시뮬레이션한다. 이를 이용해 경계 구동 XXZ 스핀 체인에서 환경 결합 강도와 이방성에 따른 얽힘 전이를 확인하고, 얽힘 스케일링이 스핀 전류와 깊은 연관을 가짐을 밝힌다.

상세 분석

논문은 먼저 양자 마스터 방정식의 시간 전파 연산자를 2차 시즈키-투터 분해로 분리한 뒤, 단위 연산은 시간 의존 변분 원리(TDVP)로, 비단위 연산은 로컬 크라우스 채널을 통해 구현한다. 핵심은 TTO가 양자 상태를 두 개의 브랜치 P와 P† 로 표현하고, 루트 텐서 R 에 시스템 절반 간 얽힘 정보가 집중된다는 점이다. 루트 텐서만으로 엔탱글먼트 오브 포메이션과 로그 네거티비티 같은 얽힘 단조량을 효율적으로 계산할 수 있다. 압축 단계에서는 SVD를 이용해 결합 차원 χ와 크라우스 차원 K를 사전에 정의된 최대값 χ_max, K_max 이하로 유지한다. 이 절차는 양자 TEA 라이브러리에 구현돼 정확 대각화가 불가능한 수백 스핀 규모까지 확장 가능하다.

실험 모델로는 1차원 XXZ 스핀 체인을 선택하고, 양쪽 끝을 마르코프ian 환경에 연결해 최대 구동 μ=1 을 적용한다. 이 설정은 기존 연구에서 전송 특성이 이방성 Δ 에 따라 ballistic(Δ<1), sub‑diffusive(Δ=1), insulating(Δ>1) 세 단계로 구분된다고 알려져 있다. 초기 상태를 완전 편극 |Z−⟩ 로 두고, 급격히 환경 결합 γ 을 켜는 퀀치 quench 을 수행한다. 시뮬레이션 결과, 스핀 전류는 초기에는 빛 원뿔 형태로 전파되며, 전류가 체인 중앙에 도달할 때마다 양쪽 절반의 엔트로피와 상호 정보가 급증한다. 특히 로그 네거티비티 N_L 은 ballistic 구간에서 부피 법칙을 따르는 큰 값을 보이며, sub‑diffusive 구간에서는 상호 정보와는 달리 거의 변하지 않는다. insulating 구간에서는 전류와 얽힘 모두 거의 성장하지 않아 전송이 국소화됨을 확인한다.

또한, 정상 상태에서의 얽힘 스케일링을 조사했을 때, ballistic 경우 N_L 이 시스템 크기에 비례해 선형적으로 증가해 부피 얽힘을 나타낸다. 반면 sub‑diffusive 구간에서는 N_L 과 I_{L:R} 이 모두 로그 스케일로 증가해 면적 법칙에 가까운 양상을 보이며, 이는 전송이 확산적이지만 완전한 열적 혼합은 아니란 것을 의미한다. 크라우스 차원 K 와 결합 차원 χ 의 성장률도 전송 구간에 따라 크게 달라져, 계산 복잡도가 ballistic 구간에서 가장 높고, insulating 구간에서는 상대적으로 낮다. 이러한 결과는 얽힘과 전류 사이의 정량적 연결고리를 제공하며, 경계 구동 시스템을 얽힘 생성 자원으로 활용할 가능성을 시사한다.