행렬곱 기반 투사 파동함수 앙상블을 이용한 양자 다체계 바닥상태

초록

본 논문은 행렬곱 연산자(MPO) 형태의 투사 연산자를 이용해 기존 시도 파동함수의 단거리 얽힘을 변분적으로 강화하는 새로운 방법을 제안한다. 변분 몬테카를로(VMC)로 MPO의 행렬 원소를 최적화하고, 장거리 얽힘은 초기 파동함수에 맡긴다. 1차원 스핀리스 페르미온 모델에 적용해 효율성을 확인했으며, 텐서곱 기반 투사 연산자를 사용해 고차원 일반화 가능성을 제시한다.

상세 분석

이 연구는 기존 변분 파동함수 접근법의 한계를 보완하기 위해 행렬곱 연산자(MPO)를 투사 연산자로 활용한다는 점에서 혁신적이다. 전통적인 변분 몬테카를로(VMC)에서는 파동함수의 형태를 고정하고 파라미터만 조정하는데, 여기서는 초기 파동함수에 MPO 형태의 연산자를 곱함으로써 파동함수 자체를 동적으로 변형한다. MPO는 1차원 양자 시스템에서 효율적인 얽힘 표현을 가능하게 하는 행렬곱 상태(MPS)의 연산자 버전으로, 제한된 차원(보통 작은 bond dimension) 내에서 복잡한 상관관계를 포착한다. 논문은 두 단계의 전략을 제시한다. 첫 번째는 장거리 얽힘을 포함한 베이스 파동함수(예: 자유 페르미온 Slater determinant 혹은 Jastrow‑보강 파동함수)를 선택하는 것이며, 두 번째는 해당 파동함수에 MPO 투사를 적용해 단거리 얽힘을 보강한다. MPO의 행렬 원소는 VMC 샘플링을 통해 에너지 기대값을 최소화하도록 변분 최적화된다. 이 과정은 기존 MPS 기반 변분 방법과 달리, 초기 파동함수가 이미 장거리 상관을 내포하고 있기 때문에 작은 bond dimension의 MPO만으로도 높은 정확도를 얻을 수 있다. 실험적으로는 1차원 스핀리스 페르미온 모델(특히 Luttinger liquid 영역)에서 기존 Jastrow 파동함수 대비 에너지 오차가 크게 감소했으며, MPO 차원을 늘릴수록 수렴 속도가 급격히 빨라지는 것을 확인했다. 또한, MPO 투사는 텐서 네트워크 형태로 일반화될 수 있어 2차원 이상에서도 비슷한 원리로 적용 가능하다는 점을 제시한다. 이 방법은 기존 DMRG이나 PEPS와 같은 고차원 텐서 네트워크 방법이 갖는 계산 복잡도 문제를 완화하면서도, 변분 몬테카를로의 샘플링 효율성을 유지한다는 장점을 가진다.