Tailored coupled cluster method with sample-based quantum diagonalization: Application to titanium-based metallocene catalytic reactions for 1-hexene production

기술 요약: 샘플 기반 양자 대각화법을 결합한 맞춤형 커플드 클러스터 방법

문제 정의
1-헥센 생산을 위한 티타늄 기반 메탈로센 촉매 반응의 선택성을 정확하게 예측하기 위해서는, 경쟁하는 전이 상태(에틸렌 삽입 vs 베타-수소 이동) 사이의 자유 에너지 차이를 1 kcal/mol 이내의 정확도로 계산해야 한다. 현재의 전자 구조 방법론은 이분법적 한계에 직면해 있다: 밀도 범함수 이론(DFT)은 계산 효율성은 높지만 범함수 의존성으로 인해 필요한 정확도가 부족한 경우가 많으며, 고수준 파동 함수 방법(예: CCSD(T))은 실제 전이 금속 시스템에 적용하기에는 계산 비용이 너무 많이 든다. 또한, 샘플 기반 양자 대각화(SQD)와 같은 양자 컴퓨팅 접근 방식은 현재 하드웨어 제약으로 인해 화학적으로 정확한 결과를 얻는 데 필요한 전체 동적 상관관계를 포착하기에는 활성 공간(active space)이 너무 작다는 한계가 있다. 결과적으로, 활성 공간 내의 정적 상관관계(static correlation)와 외부 궤도에서의 동적 상관관계(dynamical correlation)를 동시에 처리하는 것은 여전히 중요한 과제로 남아 있다.

방법론
본 연구는 활성 공간 크기와 상관관계 처리의 한계를 해결하기 위해 샘플 기반 양자 대각화(SQD)와 맞춤형 커플드 클러스터(Tailored Coupled Cluster, TCC) 이론을 통합하는 하이브리드 양자-고전 프레임워크를 제안한다.

1. 대상 시스템: 양이온성 ansa-사이클로펜타다이엔릴-아렌 티타늄 촉매의 최소 분자 모델을 사용하여 두 가지 전이 상태, T3(에틸렌 삽입)와 T4(베타-수소 이동)의 상대 에너지를 조사하였다. 일관된 궤도 프레임워크를 보장하기 위해, T3 상태에 명시적인 에틸렌 분자(T3(+C2H4))를 추가하여 T4의 원자 및 전자 수를 일치시켰다.
2. 정적 상관관계를 위한 SQD: 활성 공간 전자 구조를 처리하기 위해 SQD를 사용하였다. IBM의 Heron 양자 프로세서에서 LUCJ(Local-Unitary Cluster Jastrow) 안사(ansatz)를 사용하여, SQD는 부분 공간 내에서 구성된 해밀토니안을 구축하기 위해 비트스트링을 확률적으로 샘플링한다. 이를 통해 고전적인 완전 활성 공간 구성 상호작용(CASCI)으로는 조합론적으로 불가능한 크기의 활성 공간(최대 30개의 공간 궤도)을 처리할 수 있다. 이 방법은 하드웨어 노이즈를 완화하고 정확한 전자 수와 스핀을 보장하기 위해 자기 일관적 구성 회복(self-consistent configuration recovery)을 활용한다.
3. 동적 상관관계를 위한 TCC: SQD로부터 유도된 파동 함수는 TCC 계산의 참조 상태 역할을 하였다. 이 프레크임워크에서 활성 공간의 클러스터 진폭($\hat{T}_{AS}$)은 SQD 결과에 따라 "맞춤형"(고정)되어 활성 공간 내의 정적 상관관계와 고차 들뜸(triples, quadruples)을 포착한다. 나머지 외부 공간의 진폭($\hat{T}_{Ext}$)은 표준 커플드 클러스터 싱글 및 더블(CCSD) 방정식을 통해 최적화되어 동적 상관관계를 포착한다.
4. 확장: 정확도를 향상시키기 위해 섭동 삼중항 확장인 TCC(T)를 추가로 적용하였다. 활성 공간 에너지가 낮은 수준의 CCSD 근사가 아닌 SQD 수준을 반영하도록 에너지 보정 체계를 활용하였다.
5. 계산 설정: 모든 post-DFT 단일 지점 계산은 전자 상관 효과를 격리하기 위해 기상(gas phase)에서 수행되었으나, 기하 구조 최적화에는 극성 연속체 모델(PCM)을 통한 용매 효과(톨루엔)가 포함되었다.

주요 결과

• SQD 단독 사용 시의 수렴 문제: SQD를 단독으로 적용했을 때, 현재 가능한 활성 공간 내에서 전이 상태의 상대 에너지가 수렴되지 않았다. 결과는 활성 공간 크기에 민데 민감하게 반응하였으며, 이는 정적 상관관계만으로는 이 시스템을 설명하기에 불충분함을 나타낸다.
• SQD-TCC의 성공: 하이브리드 SQD-TCC 접근 방식은 상대 에너지에 대해 합리적으로 수렴된 결과를 도출하였다. TCC 프레임워크를 통해 동적 상관관계를 통합함으로써, 이 방법은 전자 상관관계의 균형 잡힌 처리를 성공적으로 수행하였다.
• 정적 상관관계의 영향: 표준 CCSD(T)와 SQD-TCC(T) 결과 사이에서 상당한 차이(1 kcal/mol 초과)가 관찰되었다. 이는 이 티타늄 기반 시스템에서 정적 상관관계 효과가 매우 중요하며, 이를 무시하는 것(표준 단일 참조 CCSD(T)와 같이)이 예측된 선택도에 상당한 오류를 초과할 수 있음을 강조한다.
• 활성 공간 확장성: 본 연구는 SQD가 가용한 계산 자원(최대 192개 CPU 코어 및 4 TiB 메모리) 내에서 고전적인 CASCI 방법으로는 실질적으로 접근 불가능한 활성 공간(예: 30개 궤도)을 처리할 수 있음을 입증하였다.

의의 및 주장
본 논문은 양자 유도 다중 참조 파동 함수와 고전적 커플드 클러스터 이론의 통합이 복잡한 촉매 시스템에 대한 화학적으로 정확한 시뮬레이션을 향한 실행 가능한 경로를 제공한다고 주장한다. 구체적으로, 본 연구는 다음을 입증하였다:

1. 하이브리드의 필요성: 티타늄 메탈로센 촉매 작용의 신뢰할 수 있는 설명을 위해서는 정적 및 동적 전자 상관관계가 모두 필수적이다.
2. 방법론적 발전: SQD-TCC 프레임워크는 현재의 양자 컴퓨팅 능력(제한된 활성 공간)과 현실적인 촉매 시뮬레이션의 요구사항(큰 활성 공간 + 동적 상관관계) 사이의 간극을 성공적으로 메운다.
3. 선택도 예측: 전자 상관관계의 처리는 반응의 예측된 선택도에 큰 영향을 미치며, 하이브리드 방법은 강한 정적 상관관계를 가진 시스템에 대해 표준적인 고수준 고전 방법보다 더 견고한 접근 방식을 제공한다.

저자들은 본 연구가 현재 기상에서의 퍼텐셜 에너지 표면에 집중하고 있으나, 이 방법론이 상관 파동 함수 수준에서 자유 에너지 기여 및 용매 효과를 포함하는 향후 연구를 위한 토대를 마련했다고 결론지었다.