Quantum Computing for Electronic Circular Dichroism Spectrum Prediction of Chiral Molecules

기술 요약: 키랄 분자의 전자 원편광 이색성(ECD) 스펙트럼 예측을 위한 양자 컴퓨팅

문제 정의
전자 원편광 이색성(ECD) 분광법은 약물 설계 및 거울상 선택적 합성의 안전성을 위해 필수적인 키랄 분자의 절대 구조를 결정하는 데 중요한 도구이다. 그러나 이론적 계산의 높은 비용은 ECD의 실질적인 예측 모델링 적용을 저해한다. 시간 의존 밀도 범함수 이론(TD-DFT) 및 고수준 파동 함수 방법(예: Coupled Cluster)과 같은 기존의 제1원리 접근 방식은 분자의 크기와 복잡성이 커짐에 따라 기하학적 구조 샘플링 및 들뜬 상태 계산에 과도한 비용이 소요된다. 머신러닝은 속도 면에서 이점이 있지만, 복잡하거나 강하게 상관된 시스템에 대해 신뢰할 수 있는 절대 구조 할당에 필요한 물리적 엄밀성과 전이 가능성이 부족한 경우가 많다. 따라서 고전적 알고리즘의 지수적 스케일링 없이도 제1원리의 정확도를 유지할 수 있는 확장 가능한 물리 기반 계산 방법론이 절실히 필요하다.

방법론
저자들은 ECD 스펙트럼을 예측하기 위해 양자 방정식-운동(qEOM) 형식론과 통합된 변분 양자 프레임워크를 소개한다. 이 워크플로우는 멀티 GPU 및 양자 처리 장치(QPU) 가속 기반의 하이브리드 양자-고전 환경 내에서 작동한다. 주요 방법론적 단계는 다음과 같다:

1. 활성 공간 근사(Active-Space Approximation): 큐비트 요구 사항을 관리하기 위해 저자들은 활성 공간 근사법을 사용한다. 고전적 상관 계산(MP2/cc-pVDZ)에서 도출된 자연 궤도 점유수(NOON)를 사용하여, 정적 및 동적 상관 관계를 포착하는 데 필수적인 부분 점유 궤도들의 하위 집합을 선택한다.
2. 해밀토니안 매핑(Hamiltonian Mapping): 활성 공간 전자 해밀토니안은 2차 양자 연산자를 사용하여 구축되며, 조던-위그너(Jordan–Wigner) 변환을 통해 큐비트 표현으로 매핑된다. 이는 약 20~24개의 큐비트(CAS(10,10) 및 CAS(12,12) 활성 공간에 해당)에 작용하는 큐비트 해밀토니안을 생성한다.
3. 바닥 상태 최적화(VQE): 변분 양자 고유치 계산기(VQE)를 사용하여 상관된 바닥 상태 파동 함수($|\Psi_0\rangle$)와 에너지를 결정한다. 유니터리 결합 클러스터(UCC)에서 영감을 얻은 안사츠(ansatz)가 사용되며, 해밀토니안의 기대값을 최소화하기 위해 고전적 최적화 도구인 COBYLA에 의해 최적화된다.
4. 들뜬 상태 계산(qEOM): 양자 방정식-운동 형식을 사용하여 들뜬 상태를 생성한다. 들뜸 연산자(단일 및 이중 들뜸)가 VQE로 준비된 바닥 상태에 작용한다. 들뜸 에너지와 파동 함수 계수는 양자 장치에서 평가된 교환자(commutator)의 행렬 요소를 포함하는 일반화된 고유값 문제를 푸는 과정을 통해 얻어진다.
5. 카이로프티컬(Chiroptical) 성질 계산: 전기적($\mu_{0\mu}$) 및 자기적($m_{0\mu}$) 전이 쌍극자 모멘트는 qEOM 들뜬 상태를 사용하여 계산된다. 이들은 회전 강도($R_{0\mu}$)를 결정하는 데 결합되며, 이 값은 들뜸 에너지와 함께 ECD 스펙트럼을 정의한다.
6. 구현: 시뮬레이션은 CUDA-Q 프레임워크 내에서 NVIDIA의 멀티 GPU 상태 벡터 시뮬레이터를 활용한다. 파울리 문자열(Pauli strings)의 평가를 병렬화하여 실행 시간의 병목 현상을 완화하기 위해, 해밀토니안 기대값은 네 개의 가상 QPU에 분산되어 계산된다.

주요 기여

• 프레임워크 개발: 저자들은 바닥 상태 및 들뜬 상태 특성에 대한 완전한 양자 일관적 처리를 사용하여 전자 구조 모델링과 ECD 스펙트럼 예측을 연결하는 엔드 투 엔드 양자-고전 워크플로우를 제시한다.
• 확장성 입증: 이 방법은 20~24개의 큐비트에 매핑된 활성 공간을 사용하여 단일 및 다중 키랄 중심을 가진 12종의 임상적으로 유의미한 키랄 약물을 대상으로 적용되었다.
• 검증: 양자 결과는 확립된 고전적 방법들과 체계적으로 검증되었다. 바닥 상태 에너지는 결합 클러스터 단일 및 이중 들뜸(CCSD) 및 완전 활성 공간 구성 상호작용(CASCI)과 비교 벤치마킹되었다. 들뜬 상태 에너지 및 ECD 스펙트럼은 동일한 활성 궤도 공간 내의 CASCI와 비교되었다.

결과
본 연구는 VQE–qEOM 프레임워크가 고전적 CASCI 참조 계산과 거의 정량적인 일치를 달야함을 보여준다:

• 스펙트럼 정확도: 예측된 ECD 스펙트럼은 단일 및 다중 키랄 중심 분자에 대해 스펙트럼 선형, 코튼 효과(Cotton-effect)의 부호 및 상대적 피크 강도를 정확하게 재현한다.
• 거울상 이성질체 식별: 프레임워크는 R- 및 S-거울상 이성질체 사이의 거울 대칭성을 올바르게 포착하며, 카이로프티컬 응답의 부호 반전을 정확하게 재현한다.
• 성능: D-리보스(CAS(12,12), 24 큐비트) 및 오메프라졸(Omeprazole)과 같은 분자의 경우, 양자 결과는 고전적 참조값과 정량적인 일치를 보인다. 트레오닌(Threonine)과 같이 유연한 분자에서는 상관 효과에 대한 민감도가 높아져 미세한 피크 강도 편차가 관찰되었으나, 전반적인 스펙트럼 경향과 거울상 이성질체 할당은 일관되게 유지되었다.
• 효율성: 멀티 GPU 병렬화 전략은 단일 QPU 실행에 비해 실행 시간을 크게 단축하여 더 큰 활성 공간에서의 시뮬레이션을 가능하게 했다.

의의 및 주장
본 논문은 이 연구가 제1원리 기반의 카이로프티컬 스펙트럼 예측을 위한 실용적이고 확장 가능한 경로로서 양자 알고리즘을 확립한다고 주장한다. VQE의 변분적 견고함과 방정식-운동 기법의 개념적 구조를 결합함으로써, 이 프레임워크는 데이터 중심 모델 및 계산적으로 매우 까다로운 고전적 방법론에 대한 물리적으로 근거 있는 대안을 제공한다. 저자들은 본 접근 방식이 현재의 노이즈가 있는 중간 규모 양자(NISQ) 하드웨어 역량과 제약 조건 사이의 간극을 성공적으로 메운다고 주장한다. 결과는 양자 기반 분자 분광법이 화학 및 제약 분야의 관련 시스템에 적용 가능하다는 것을 시사하며, 향후 양자 하드웨어가 발전함에 따라 더 큰 분자와 더 복잡한 입체 화학 환경으로 확장할 수 있는 토대를 제공한다.