Thermodynamic sampling of materials using neutral-atom quantum computers

기술 요약: 중성 원자 양자 컴퓨터를 이용한 재료의 열역학적 샘플링

문제 정의
기술적으로 중요한 재료들은 종종 구성적 무질서(예: 고용체, 공공, 치환형 합금)를 나타낸다. 이러한 재료의 특성을 예측하기 위해서는 정확한 에너지 모델뿐만 아니라, 특정 온도와 화학적 포텐셜에서의 평형 상태를 결정하기 위해 방대한 원자 구성 공간을 샘플링할 수 있는 방법론이 필요하다. 몬테카를로 샘플링 및 시뮬레이티드 어닐링과 같은 고전적 접근 방식은 시스템 크기가 증가함에 따라 계산적 병목 현상에 직면하며, 이로 인해 전수 조사는 불가능해진다. 양자 어닐링이 이 목적으로 탐구되어 왔으나, 실제 재료 모델(특히 밀도 범함수 이론(DFT)에서 유도된 모델)을 현재의 양자 하드웨어에 매핑하는 것은 상당한 과제를 안겨준다. 여기에는 상호작용 범위, 에너지 스케일에 대한 하드웨어 제약과 주기적 경계 조건을 네이티브하게 지원하지 못하는 중성 원자 어레이의 특성이 포함된다. 또한, 양자 최적화기의 출력 분포에 대한 엄격한 열역학적 해석을 확립하는 것은 여전히 미해결 과제로 남아 있다. 왜냐하면 많은 장치들이 평형 열 앙상블보다는 결맞는 유니터리 역학(coherent unitary dynamics)으로 더 잘 설명되기 때문이다.

방법론
저자들은 질소가 도핑된 그래핀의 열역학적 특성을 중성 원자 양자 하드웨어(구체적으로 QuEra Aquila 장치)에 적합한 리드베리 원자 해밀토니안(Rydberg-atom Hamiltonian)으로 매핑하는 프레임워크를 제안한다. 방법론은 다음 단계로 진행된다:

1. DFT 데이터 생성: 78개 원자로 구성된 슈퍼셀에 대해 주기적 DFT(CRYSTAL23 코드)를 사용하여 질소 도핑된 그래핀의 형성 에너지를 계산한다. 대칭 독립 구성(SIC)을 생성하고 완화(relaxation)하며, 훈련 및 테스트 세트를 확장하기 위해 대칭 동등 구성(SEC)을 도출한다.
2. 해밀토니안 매핑: DFT 형성 에너지를 리드베리 해밀토니안에 피팅한다. 이는 온사이트 항($V$)과 거리 의존적 쌍 상호작용($C_6/R^6$)으로 구성된 모델과 DFT 에너지 사이의 오차를 최소화하는 과정을 포함한다. 이 피팅을 통해 온사이트 에너지 $V^{DFT}$와 유효 최근접 이웃 거리 $R^{DFT}_{NN}$을 얻는다.
3. 하드웨어 리스케일링 전략: DFT에서 유도된 파라미터와 하드웨어 역량 사이에는 결정적인 불일치가 존재한다. 하드웨어는 최소 원자 간격($R^{min}_{NN} = 4\,\mu\text{m}$)을 요구하는 반면, 피팅된 $R^{DFT}_{NN}$은 이보다 훨씬 작은 약 $1.6\,\mu\text{m}$이다. 결과적으로, 하드웨어의 네이티브 2체 상호작용 강도는 대상 재료의 상호작용보다 수십 배 더 약하다.
   • 이를 해결하기 위해 저자들은 리스케일링 인자 $\alpha_v = (R^{min}_{NN} / R^{DFT}_{NN})^6 \approx 236.69$를 도입한다.
   • 하드웨어의 접근 가능한 에너지 창 내에 맞추기 위해 온사이트 에너지와 쌍 상호작용 모두를 $\alpha_v$만큼 스케일링한다.
   • 이 리스케일링은 하드웨어 샘플링 온도 $T$와 유효 재료 온도 $T' = \alpha_v T$ 사이의 대응 관계를 구축한다.
4. 화학적 포텐셜 제어: 양자 어닐러의 글로벌 레이저 디튜닝($\Delta_g$)은 재료의 그랜드 캐노니컬 화학적 포텐셜 차이($\Delta\mu$)에 매핑된다. 리스케일링 인자를 통해 하드웨어의 제한된 디튜닝 범위를 재료 모델의 화학적으로 유의미한 $\Delta\mu$ 범위까지 커버할 수 있다.
5. 검증:
   • 28-사이트 시스템: 모든 가능한 구성을 전수 조사(exhaustive enumeration)하여 매핑을 검증한다.
   • 78-사이트 시스템: 전수 조사가 불가능한 더 큰 시스템의 경우, 무작위 몬테카를로(MC) 샘플링과 비교하여 벤치마킹한다.
   • 온도 튜닝: 하드웨어에서 최근접 이웃 간격($R_{NN}$)을 변화시켜 유효 온도를 실험적으로 조절함으로써 샘플링된 분포에 대한 제어 능력을 입증한다.

주요 기여

• 리스케일링 프레임워크: 본 논문은 DFT 유도 에너지 스케일과 현재 중성 원자 하드웨어의 제한된 동적 범위 사이의 격차를 메우는 체계적인 리스케일링 절차($\alpha_v$)를 도입한다. 이를 통해 하드웨어가 유효 온도 $T'$에서 대상 재료에 대응하는 분포를 샘플링할 수 있게 한다.
• 열역학적 해석: 하드웨어 제어 파라미터(레이저 디튜닝 및 원자 간 간격)를 열역학적 변수(화학적 포텐셜 및 온도)에 연결하는 엄격한 유도를 제공하여, 양자 어닐링 출력을 볼츠만 유사 앙상블로 해석할 수 있게 한다.
• 하드웨어 제약 처리: 기하학적 제약(최소 원자 간격, 2D 배열)과 중성 원자 장치의 주기적 경계 조건 부재를 명시적으로 다루며, 유한한 나노플레이크를 벌크 거동을 대표하는 모델로 취급한다.
• 개념 증명: 프로그래머블 양자 어닐러를 사용하여 도핑된 재료 시스템의 열역학적 특성(평균 조성 및 에너지 분포)을 성공적으로 추출함을 입상적으로 보여준다.

결과

• 28-사이트 벤치마크: 28-사이트 그래핀 나노플레이트에 대해, 양자 어닐러의 출력(평균 질소 농도 대 화학적 포텐셜)은 전수 조사로부터 얻은 정확한 결과와 밀접하게 일치한다. 최적의 적합은 유효 샘플링 온도 $T = 41\,\mu\text{K}$ (이에 따른 재료 유효 온도 $T' \approx 9.7\,\text{mK}$)에 해당한다. 양자 분포와 고전적 분포 사이의 총 변동 거리(TVD)는 0.294로, 저에너지 상태가 온도에 민감하다는 점에 기인한 편차를 포함하여 양호한 일치도를 보인다.
• 78-사이트 벤치마크: 78-사이트 시스템의 경우, 양자 어닐러는 저에너지 구성을 우선적으로 샘플링하는 반면, 무작위 몬테카를로 샘플링은 전체 에너지 지형을 균일하게 탐색한다. 양자 장치는 무작위 샘플링으로는 포착하기 어렵고 조합론적으로 희귀한 에너지 0(비상호작용 질소 배열)의 구성을 성공적으로 식별한다.
• 온도 제어: 최근접 이웃 간격을 $4\,\mu\text{m}$에서 $5\,\mu\text{m}$로 증가시킴으로써, 유효 샘플링 온도를 $\approx 9.7\,\text{mK}$에서 $\approx 36\,\text{mK}$로 튜닝할 수 있음을 성공적으로 입증하였으며, 이는 리스케일링 전략의 열역학적 일관성을 확인시켜 준다.

의의 및 주장
저자들은 본 연구를 현실적인 재료 모델(DFT로부터 유도된)을 리드베리 해밀토니안으로 매핑하는 **개념 증명(proof-of-concept)**으로 규정한다.

• 주요 의의는 심각한 에너지 스케일 불일치에도 불구하고 가용한 양자 하드웨어에 구현 가능한 현실적인 재료 모델을 구현할 수 있는 워크플로우를 확립했다는 점에 있다.
• 본 논문은 제안된 리스케일링 전략을 통해 하드웨어가 열역학적으로 해석 가능한 분포, 즉 유효 온도에서의 볼츠만 유사 앙상블을 샘플링할 수 있게 한다고 주장한다.
• 저자들은 현재의 한계를 인정하며, 본 방법론이 2체 상호작용 및 2D 시스템으로 제한된다는 점을 명시한다. 특히 강한 국소적 상호작용이나 다체 효과(예: 인-질소 공동 도핑)가 있는 시스템은 현재의 2체 해밀토니안으로 정확히 포착할 수 없으며, 더 높은 차수의 결합이나 더 유연한 상호작용 모델이 필요함을 밝히고 있다.
• 본 연구는 특정 시스템에 대해 고전적 방법을 능가한다고 주장하는 것이 아니라, 미래의 재료 발견을 위한 하이브리드 워크플로우의 구성 요소로서 중성 원자 양자 컴퓨터를 사용하는 것의 타당성을 입증하는 데 목적이 있다.