Nonlocal magic of fermions is bounded by entanglement fluctuations and is logarithmic for free fermions

기술 요약: 페르미온의 비국소적 매직은 얽힘 변동량에 의해 제한된다

문제 정의
하트리-포크(Hartree–Fock) 및 BCS 평균장 상태를 포함한 자유 페르미온 상태(가우시안 상태)는 클래식하게 시뮬레이션 가능하지만, 큐비트로 표현될 경우 일반적으로 상당한 "매직"(non-stabilizerness)을 가집니다. 매직은 스테빌라이저 형식(stabilizer formalism)을 통한 클래식 시뮬레이션의 난이도와 결함 허용 양자 컴퓨팅(특히 비-클리포드 $T$ 게이트의 필요 개수)을 위한 자원 비용을 정량화합니다. 국소 유니터리(local unitaries)는 국소적 매직을 제거할 수 있지만, 두 서브시스템 사이의 얽힘에 저장된 "비국소적 매직(nonlocal magic)"은 남아 있습니다.

본 연구가 다루는 핵심 질문은 다음과 같습니다: 자유 페르미온 상태는 실제로 얼마나 많은 비국소적 매직을 보유하고 있는가? 기존 연구에서는 만약 국소 연산이 가우시안 유니터리(자유 페르미온을 자유 페르미온으로 매핑하는 연산)로 제한될 경우, 비국소적 매직이 부분적으로 얽힌 페르미온 쌍의 개수 $m$에 따라 선형적으로 증가한다는 것을 밝혀냈습니다. 그러나 일반적인 페르로의 국소 유니터리는 오직 페르미온 패리티 초선택 규칙(각 영역의 패리티와 교환 법칙을 만족해야 함)에 의해서만 제약을 받으며, 반드시 가우시안일 필요는 없습니다. 본 논문은 일반적인 패리티 보존(non-Gaussian) 국소 유니터리를 허용할 경우 비국소적 매직이 감소하는지, 그렇다면 어떤 메커니즘과 범위 내에서 감소하는지를 조사합니다.

방법론
저자들은 이분할 $A|B$에 대한 순수 페르미온 상태를 분석합니다. 이들은 두 가지 주요 비국소적 매직 척도를 사용합니다:

1. 비국소 스테빌라이저 렌니 엔트로피($M_{\alpha}^{NL}$): 국소 유니터리에 대한 최적화를 거친 후 달성 가능한 최소 SRE.
2. 비국소 스테빌라이저 충실도($\Phi$): (유니터리 적용 후의) 상태와 임의의 스테빌라이저 상태 사이의 최대 중첩. 이는 비국소 최소 상대적 매직 엔트로피($D_{min}^{NL} = -\log_2 \Phi$)와 관련이 있습니다.

방법론은 다음과 같습니다:

• 패리티를 포함한 슈미트 분해(Schmidt Decomposition with Parity): 상태를 명확한 패리티(짝수/홀수 섹터)를 가진 슈미트 벡터로 분해합니다. 저자들은 "평탄한(flat)" 타겟 스펙트럼(벨 쌍, Bell pairs)에 맞추기 위해 짝수 및 홀수 섹터에서 가장 큰 슈미트 값을 선택하는 과정을 통해, 패리티-짝수 국소 유니터리에 의해 달성 가능한 최대 충실도 $\Phi_S^{NL}$에 대한 정확한 공식을 도출합니다.
• 얽힘 용량($C_E$): 얽힘 해밀토니안의 분산(또는 얽힘 에너지 $X_i = -\log_2 \lambda_i$의 분산)을 $C_E$로 정의합니다. 저자들은 얽힘 에너지를 확률 변수로 취급하고, 통계적 논증(중심한계정리, 베리-에센 정리)을 사용하여 이 분포의 퍼짐(spread)을 바탕으로 충실도를 제한합니다.
• 연산의 비교: 자유 페르미온 구조를 유지하는 국소 가우시안 유니터리와 일반적인 패리티 보존 유니터리(포크 상태 치환을 포함하는)의 성능을 대조합니다.

주요 기여 및 결과

1. 비국소적 매직의 선형 감소에서 로그 감소로의 전환:
   본 논문은 가우시안 제한 하에서 관찰된 외연적(extensive, $m$에 비례하는) 비국소적 매직이 해당 제한의 산물임을 입증합니다. 일반적인 패리티 보존 국소 유니터리가 허용될 때, 비국소적 매직은 크게 감소합니다.

   • 충실도: 최적의 스테빌라이저 상태와의 충실도는 부분적으로 얽힌 쌍들에 대해 지수적 감소($e^{-\Theta(m)}$)가 아닌 $O(1/\sqrt{m})$로 감소합니다.
   • 매직 척도: 결과적으로, 매직 척도는 $\Theta(m)$에서 $\Theta(\log m)$으로 떨어집니다. 구체적으로, $D_{min}^{f,NL} \approx \frac{1}{2} \log_2 m$ 이고 $M_2^{f,NL} \approx \frac{3}{2} \log_2 C_E$ 입니다.

2. 얽힘 용량($C_E$)의 역할:
   저자들은 다음과 같은 근본적인 경계 조건을 설정합니다: 명확한 패리티를 가진 모든 순수 페르미온 상태의 비국소적 매직은 얽힘의 용량에 의해 제한됩니다.
   $$\Phi_S^{NL} \geq \frac{0.478}{\sqrt{C_E + 5}}$$
   이는 상태가 큰 비국소적 매직을 가지려면 슈미트 가중치(Schmidt weights)가 여러 차수의 크기로 넓게 퍼져 있어야 함을 의미합니다(즉, 얽힘 에너지의 분산이 커야 함). 자유 페르미온 쌍은 독립적으로 $C_E$에 기여하며, 각 쌍은 최대 상수 $c^* \approx 0.44$를 더하므로 $C_E$는 $m$에 따라 선형적으로 성장합니다. 이러한 $C_E$의 선형 성장은 비국소적 매직의 로그 스케일링으로 이어집니다.

3. 가우시안 vs 비-가우시안 유니터리의 최적성:

   • 작은 시스템 ($m \leq 3$): 최대 3개의 얽힌 쌍에 대해서는 국소 가우시안 유니터리가 모든 패리티-짝수 유니터리 중에서 최적임이 증명되었습니다. 여기서 가우시안 연산과 비-가우시안 연산의 구분은 필수적이지 않습니다.
   • 큰 시스템 ($m \geq 4$): 4개 이상의 쌍에 대해서는 일반적인 패리티 보존 유니터리(특히 포크-상태 치환을 포함하는 것)가 가우시안 유니터리보다 엄격하게 우수한 성능을 보이며, 가우시안 제한 하의 최소값보다 낮은 SRE를 달성합니다.

4. 경계의 통계적 기원:
   충실도의 $1/\sqrt{m}$ 감소는 통계적 현상임이 밝혀졌습니다. 벨 쌍(Bell pairs)은 평탄한 얽힘 스펙트럼(제로 용량)을 가집니다. 벨 쌍으로부터 멀어지기 위해서는(즉, 높은 비국소적 매직을 가지려면), 상태의 슈미트 가중치가 매우 넓게 분포되어야 합니다. $C_E$는 이 퍼짐을 측정합니다. 자유 페르미온의 경우, 얽힘 에너지는 독립적인 기여들의 합이므로 표준 편차가 $\sqrt{m}$의 스케일로 변화하며, 이것이 $1/\sqrt{m}$ 충실도 경계를 결정합니다.

5. 임계 계(Critical Systems):
   임계 자유 페르미온 체인(예: XX 체인)의 경우, $C_E$는 서브시스템 크기 $\ell$에 대해 로그 함수적으로 성장합니다. 따라서 비국소적 매직은 이중 로그($\log \log \ell$)로 성장하며, 이는 큐비트 체인에 대한 연구 결과와 일치합니다.

의의 및 주장
본 논문은 자유 페르미온 상태(하트리-포크 및 BCS 상태 포함)의 거의 모든 매직은 국소적이다라고 주장합니다. 스테빌라이저 회로를 사용하여 상태의 상관관계를 준비하는 데 드는 비용인 비국소적 부분은 놀라울 정도로 작으며, 얽힌 쌍의 개수에 대해 로그 함수적으로만 증가합니다.

• 양자 시뮬레이션에 대한 시사점: 자유 페르미온은 매직을 가지고 있지만, 그 비국소적 부분은 외연적이지 않습니다. 이는 자유 페르미온을 스테빌라이저 방법으로 시뮬레이션하는 "어려움"이 장거리 얽힘의 복잡성 때문이 아니라 국소적 구조 때문임을 시사합니다.
• 결함 허용 컴퓨팅에 대한 시사점: 일반적인 패리티 보존 연산이 가능하다면, 자유 페르미온 상태를 준비하는 데 필요한 $T$ 게이트의 수는 선형적이 아니라 로그 단위로 제한됩니다.
• 보편성: $D_{min}^{NL} \leq \frac{1}{2} \log_2 C_E + O(1)$ 경계는 (자유 또는 상호작용하는) 모든 페르미온 상태에 대해 얽힘 용량을 알고 있다면 성립합니다.

저자들은 자신의 결과가 자유 페르미온이 클래식하게 시뮬레이션 가능하다는 사실(매직 척도는 스테빌라이저 시뮬레이션을 위한 비용을 나타내는 것이지, 일반적인 클래식 복잡도를 나타내는 것이 아님)을 바꾸지 않음을 명시하며, 비-가우시안 국소 유니터리 자체를 구현하는 데 드는 비용(이는 각 영역 내에서 비-클리포드 자원을 요구할 수 있음)을 추정하지는 않는다고 밝혔습니다.