[논문 리뷰] Quantum simulation of general spin-1/2 Hamiltonians with parity-violating fermionic Gaussian states

Michael Kaicher, Joseph Vovrosh외 2명|arXiv (Cornell University)|2026. 01. 20.

Quantum many-body systems인용 수 0

한 줄 요약

이 논문은 페러티-위반 페르미온 평균장 이론 PV-FMFT를 페러티-위반 페르미온 가우시안 상태 PV-FGS를 기반으로 개발하여 Colpa 매핑과 확장된 힐베르트 공간을 이용해 일반 스핀-1/2 해밀토니언의 실수시- 및 허수시간 다이나믹스를 효율적으로 시뮬레이션한다.

ABSTRACT

We introduce equations of motion for a parity-violating fermionic mean-field theory (PV-FMFT): a numerically efficient fermionic mean-field theory based on parity-violating fermionic Gaussian states (PV-FGS). This work provides explicit equations of motion for studying the real- and imaginary-time evolution of spin-1/2 Hamiltonians with arbitrary geometries and interactions. We extend previous formulations of parity-preserving fermionic mean-field theory (PP-FMFT) by including fermionic displacement operators in the variational Ansatz. Unlike PP-FMFT, PV-FMFT can be applied to general spin-1/2 Hamiltonians, describe quenches from arbitrary initial spin-1/2 product states, and compute local and non-local observables in a straight-forward manner at the same modest computational cost as PP-FMFT -- scaling as $O(N^3)$ in the worst case for a system of $N$ spins or fermionic modes. We demonstrate that PV-FMFT can exactly capture the imaginary- and real-time dynamics of non-interacting spin-1/2 Hamiltonians. We then study the post quench-dynamics of the one- and two-dimensional Ising model in presence of longitudinal and transversal fields with PV-FMFT and compute the single site magnetization and correlation functions, and compare them against results from other state-of-the-art numerical approaches. In two-dimensional spin systems, we show that the employed spin-to-fermion mapping can break rotational symmetry within the PV-FMFT description, and we discuss the resulting consequences for the calculated correlation functions. Our work introduces PV-FMFT as a benchmark for other numerical techniques and quantum simulators, and it outlines both its capabilities and its limitations.

연구 동기 및 목표

• 스핀-1/2 해밀토니언을 페르미로 매핑할 때 나타나는 페러티-위반 페르미온 시스템에 대해 수치적으로 효율적인 평균장 프레임워크를 동기 부여하고 형식화한다.
• 기존의 패러티 보존 FMFT를 확장하여 페러티-위반 가우시안 상태를 통해 임의의 스핀-1/2 해밀토니언을 다룰 수 있도록 한다.
• PV-FMFT 내에서 허용된 실수시(기저 상태) 및 실시간(동역학) 진화에 대한 명시적이고 안정적인 운동방정식을 제공한다.

제안 방법

• 패러티를 깨뜨리기 위해 페르미 이 displacement 항을 포함하여 패러티를 깨뜨리는 변분 Ansatz로 PV-FGS를 도입한다.
• 확장된 PV-FGS 공간에서 공분산 행렬 Γ에 대한 폐쇄형 운동방정식을 도출한다(허수시간: dΓ/dτ = -Hm - Γ Hm Γ; 실시간: dΓ/dt = [Hm, Γ]).
• 확장된 힐베르트 공간(N+1 모드)에서 PV 문제를 PP 프레임워크로 변환하기 위해 Colpa 매핑을 구현하고 선형 케노니컬 변환을 가능하게 한다.
• 에너지 기대값 및 변분 행렬(Hm, M 등)에 대한 명시적 표현을 제공하고 수치적 안정성과 Pfaffian 기반 평가를 논의한다.
• 비상호작용 PV 스핀 해밀타니언에 대해 정확함을 보이고, 1D 및 2D 이징 모델의 종방향 및 횡방향 필드를 가진 퀜치를 연구한다.

Figure 1: Schematic relation between spin mean-field theory (SMFT) and fermionic mean-field theory (FMFT). Each circle represents a family of Hamiltonians and states which are exactly described by a mean-field theory based on spin-1/2 product states (purple), parity-preserving fermionic Gaussian sta

실험 결과

연구 질문

• RQ1 PV-FMFT와 PV-FGS가 비상호작용 스핀-1/2 해밀토니언의 실수시 및 허수시간 다이나믹스를 정확하게 포착할 수 있는가?
• RQ2 스핀-페르미 매핑에서 발생하는 PV 항을 Gaussian 평균장 프레임워크 내에서 어떻게 효율적으로 처리할 수 있는가?
• RQ3 Colpa 매핑이 PV 문제를 PP-FMFT 방법으로 gauge 제약이나 모드 중복 없이 한 대수의 보조 모드 이상으로 확장 없이 다룰 수 있게 해주는가?
• RQ4 1D 및 2D 이징 모델의 종방향 필드를 가진 차퀸에 PV-FMFT를 적용했을 때의 가능성 및 한계는 무엇인가?

주요 결과

• PV-FMFT와 PV-FGS는 비상호작용 스핀-1/2 해밀토니언의 실수시 및 허수시간 다이나믹스를 정확하게 기술할 수 있다.
• 프레임워크는 계산적으로 다룰 수 있는 범위를 유지하며 최악의 경우 N개의 스핀 또는 페르미온 모드에서 O(N^3) 규모로 확장된다.
• PV-FMFT는 종방향 및 횡방향 필드를 갖는 이징 모델에서의 퀜치를 시뮬레이션하고 국소 자기화 및 상관 함수들을 계산할 수 있다.
• Colpa 매핑은 PV-FGS를 확장된 N+1 모드 공간에서 두 개의 PP-FGS의 선형 결합으로 매핑 가능하게 하여 PP-FMFT와 비슷한 계산 비용을 유지한다.
• 2D에서 스핀-페르미 매핑은 PV-FMFT 내에서 회전 대칭을 깨뜨릴 수 있어 계산된 상관 함수에 영향을 주고 매핑으로 인한 artefacts를 강조한다.
• PV-FMFT 프레임워크는 다른 수치 기법 및 양자 시뮬레이터의 벤치마크로 작용하며 역량과 한계를 제시한다.

Figure 2: Workflow diagram describing the steps to represent an arbitrary spin-1/2 Hamiltonian (here a two-dimensional square lattice) as a PP fermionic Hamiltonian. We consider the TFIM model where the spin-1/2 Hamiltonian is mapped to its fermionic description via the Jordan-Wigner transformation.