입자수 보존을 활용한 행렬곱 연산자 시뮬레이션 최적화

초록

본 논문은 입자수 보존이라는 전역 대칭을 행렬곱 연산자(MPO) 시뮬레이션에 적용했을 때, 기대와 달리 연산 속도가 반드시 향상되지 않으며, 오히려 고정된 입자수 제약이 엔트랄피를 증가시켜 성능을 제한한다는 점을 밝힌다.

상세 분석

양자 다체 시스템을 텐서 네트워크 형태로 기술할 때, 행렬곱 상태(MPS)는 전역 대칭을 명시적으로 구현하면 차원 축소와 수치 안정성 향상을 얻는다. 이는 각 사이트의 물리적 차원을 대칭 섹터별 블록으로 분해하고, 블록 내부에서만 연산을 수행함으로써 가능한데, 특히 입자수 보존(U(1) 대칭)은 가장 흔히 활용되는 사례이다. 그러나 본 연구는 동적 연산자 자체가 MPO 형태인 경우, 즉 시간 진화 연산자나 열역학적 리듬을 직접 다루는 상황에서 동일한 전략을 적용했을 때 발생하는 구조적 차이를 집중적으로 분석한다.

첫 번째 핵심은 MPO가 MPS와 달리 두 개의 물리 인덱스(입력·출력)를 동시에 보유한다는 점이다. 입자수 보존을 강제하면 각 텐서의 양쪽 물리 인덱스에 대해 동일한 입자수 흐름을 유지해야 하므로, 전체 네트워크에 전역적인 흐름 제약이 부과된다. 이 제약은 “전하 흐름 보존”이라는 형태로 나타나며, 각 텐서의 차원은 단순히 블록 크기의 곱이 아니라, 입자수 분포에 따라 가능한 연결 조합의 수로 제한된다. 결과적으로, 입자수 블록을 하나만 남겨두고 나머지를 버리는 “단일 섹터” 접근법은 MPS에서는 차원 감소 효과가 크지만, MPO에서는 블록 간 연결이 급격히 복잡해져 엔트랄피(양자 얽힘) 증가를 초래한다.

두 번째로, 논문은 수치 실험을 통해 두 가지 대표 모델—1차원 하버드 모델과 장거리 상호작용을 갖는 Bose-Hubbard 체계—에 대해 비교하였다. 입자수 보존을 적용한 MPO는 초기 상태에서는 차원 감소가 관찰되지만, 시간 진화가 진행될수록 전역 입자수 제약 때문에 얽힘 스펙트럼이 급격히 넓어지고, 결국 차원 절감 효과가 사라진다. 특히, 고전적인 TEBD(시간-진화 블록 디코플레이션)와 비교했을 때, 단일 섹터 MPO는 동일 정확도에서 더 많은 보조 차원을 필요로 하며, 메모리 사용량과 계산 시간 모두 증가한다.

세 번째로, 저자들은 “가상 입자수” 기법을 도입해 보조 블록을 추가함으로써 얽힘을 완화하려 시도했지만, 이는 오히려 블록 구조를 복잡하게 만들고, 구현상의 오버헤드가 커져 실용적 이득을 제공하지 못한다는 결론에 도달한다. 따라서 입자수 보존을 MPO에 적용할 때는 단순히 “섹터 차단”이 아니라, 전역 흐름을 고려한 새로운 텐서 분해 전략이 필요함을 강조한다.

결론적으로, 입자수 보존이 MPS 기반 시뮬레이션에서는 강력한 최적화 도구이지만, MPO에서는 전역 제약이 얽힘을 증가시켜 성능을 저해할 수 있다. 이는 텐서 네트워크 설계 시 대칭 활용 방법을 상황에 맞게 선택해야 함을 시사한다.