Quantum Simulation of Non-Hermitian Linear Response

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 선형 응답 이론 (Linear Response Theory) 과 그린 함수 (Green's functions) 는 외부 약한 섭동에 대한 거시적 및 강상관 계의 반응을 이해하는 보편적인 틀을 제공합니다. 최근 Pan 등 [3] 은 개방 양자계 (Open Quantum Systems) 를 설명하기 위해 비유니터리 (Non-Hermitian) 선형 응답 이론의 이론적 기반을 확립했습니다.
• 문제점:
  • 기존 양자 알고리즘은 폐쇄계 (Closed Systems) 의 에르미트 (Hermitian) 해밀토니안 하에서만 작동하며, 비유니터리 과정 (소산, 연속 측정 등) 을 다루는 데 근본적인 장벽이 존재합니다.
  • 비유니터리 시간 진화 연산자 $e^{Lt}$는 양자 게이트로 직접 실행할 수 없으며, 기존 접근법 (확장 선형 시스템 솔버, 스펙트럼 매핑 정리 등) 은 회로 깊이와 상태 준비 (State Preparation) 비용에서 지수적 오버헤드를 발생시킵니다.
  • 따라서, 비유니터리 선형 응답 이론을 실제 양자 하드웨어에서 실행 가능한 알고리즘으로 변환하는 것이 주요 난제였습니다.

2. 제안된 방법론 (Methodology)

저자는 슈뢰딩거화 (Schrödingerization) 기법을 활용하여 비유니터리 다중 시간 상관 함수를 양자 하드웨어에서 실행 가능한 유니터리 형태로 변환하는 체계적인 알고리즘을 제안합니다.

• 벡터화와 리우빌 공간 (Vectorization & Liouville Space):

  • 린드블라드 마스터 방정식 (Lindblad master equation) 을 따르는 밀도 행렬 $\rho(t)$를 리우빌 공간의 상태 벡터 $|\rho(t)\rangle$로 벡터화합니다.
  • 이를 통해 비에르미트 슈퍼연산자 $L$에 의한 동역학 $\frac{d}{dt}|\rho(t)\rangle = L|\rho(t)\rangle$로 표현됩니다.
  • 양자 회귀 정리 (Quantum Regression Theorem) 를 통해 두 점 상관 함수 $\langle A(t)B(0)\rangle$를 계산할 수 있습니다.

• 슈뢰딩거화 (Schrödingerization):

  • 비에르미트 생성자 $L$을 에르미트 성분 $H_1$과 반에르미트 성분 $H_2$로 분해합니다 ($L = H_1 - iH_2$).
  • 실수 1 차원 변수 $\xi > 0$을 도입하고 왜곡 위상 변환 (warped phase transformation) $w(t, \xi) = e^{-\xi}|\rho(t)\rangle$을 적용합니다.
  • $\xi$에 대한 푸리에 변환을 수행하여 새로운 변수 $\eta$를 도입하면, 비유니터리 동역학이 유니터리 슈뢰딩거 방정식으로 변환됩니다:
    $$i\partial_t \tilde{w} = (\eta H_1 + H_2) \tilde{w}$$
  • 여기서 확장된 해밀토니안 $H_{sch}(\eta) = \eta H_1 + H_2$는 $\eta$가 실수일 때 항상 에르미트이므로, 양자 컴퓨터에서 유니터리 진화가 가능합니다.

• 비유니터리 선형 응답 및 그린 함수:

  • 섭동 $L \to L + \delta L(t)$에 대한 응답 함수 $\chi(\tau)$를 계산할 때, 초기 상태를 $V|\rho_{eq}\rangle$로 설정하고 슈뢰딩거화된 진화를 수행함으로써 일반화된 지연 그린 함수를 추출합니다. 이는 폐쇄계의 쿠보 (Kubo) 공식이 소산이 있는 비에르미트 영역으로 확장된 형태입니다.

3. 주요 기여 (Key Contributions)

1. 알고리즘적 매핑: 비유니터리 다중 시간 상관 함수를 유니터리 형태로 변환하는 체계적인 알고리즘을 제시하여, 비에르미트 선형 응답 이론과 양자 시뮬레이션 간의 간극을 메웠습니다.
2. 최적의 상태 준비 비용: 기존 블록 인코딩 (Block-encoding) 기법이나 QSVT(Quantum Singular Value Transformation) 와 달리, 행렬 역연산에 따른 불리한 상태 준비 비용을 피하고 **최적 (Optimal)**인 상태 준비 비용을 달성했습니다.
3. 확장성 및 효율성: 정밀도 $\epsilon$에 대해 $O(\text{poly}(\log(1/\epsilon)))$로 스케일링되며, 기존 확장 (Dilation) 방법의 지수적 오버헤드를 극복했습니다.
4. 구체적 검증: 단일 큐비트 진폭 감쇠 (Amplitude Damping) 모델을 분석적으로 유도하여 프레임워크의 유효성을 증명했습니다.

4. 실험 결과 및 성능 (Results)

• 시뮬레이션 환경: 단일 큐비트 진폭 감쇠 채널을 대상으로 Qiskit Sampler 를 사용한 노이즈 없는 양자 회로 시뮬레이션을 수행했습니다.
• 수렴성 분석:
  • 푸리에 모드 수 ($N$) 와 오차: 푸리에 격자 점수 $N$을 증가시킬 때 오차가 지수적으로 감소함을 확인했습니다. 이는 상대적으로 적은 수의 보조 큐비트 ($M \approx \log_2 N$) 로도 높은 정밀도를 달성할 수 있음을 의미합니다.
  • 절단 한계 ($L$) 와 오차: 왜곡 변수의 절단 한계 $L$에 대한 오차는 $e^{-L}$로 스케일링되어 이론적 예측과 일치했습니다.
• 정확도: 고해상도 설정 ($N=200, L=40$) 에서 시뮬레이션된 응답 함수가 정확한 린드블라드 해석해 (Analytic solution) 와 완벽하게 일치함을 보였습니다.
• 회로 구현: 수정된 하마르드 테스트 (Modified Hadamard test) 를 사용하여 초기 상태 준비 및 유니터리 진화 후 관측값을 추출하는 회로 구조를 제시했습니다.

5. 의의 및 결론 (Significance)

• 이론적 확장: 비에르미트 선형 응답 이론을 양자 컴퓨팅 구현으로 성공적으로 연결하여, 구동 - 소산 (Driven-dissipative) 강상관 전자계 및 물리 화학 분야에서 양자 컴퓨팅의 응용 범위를 크게 확장했습니다.
• 실용적 가치: 복잡한 행렬 역연산 알고리즘을 우회하고 비유니터리 제한을 극복함으로써, 고차원 개방계 시뮬레이션을 위한 확장 가능 (Scalable) 하고 효율적인 알고리즘적 경로를 최초로 제시했습니다.
• 미래 전망: 이 연구는 소산이 있는 양자 시스템의 동역학적 반응을 연구하는 데 있어 새로운 표준이 될 수 있으며, 향후 실제 양자 하드웨어에서의 노이즈 내성 및 리소스 최적화 연구의 기초를 제공합니다.

핵심 요약: 이 논문은 슈뢰딩거화 기법을 통해 비유니터리 (소산이 있는) 양자 시스템의 선형 응답을 유니터리 양자 회로로 변환하는 알고리즘을 제안했습니다. 이를 통해 기존 방법들의 지수적 오버헤드를 제거하고, 최적의 상태 준비 비용으로 비에르미트 그린 함수를 효율적으로 계산할 수 있음을 이론적으로 증명하고 수치적으로 검증했습니다.