CNDO2와 INDO 모델을 위한 포트란90 기반 하트리 폭 계산 프로그램

초록

전산화학에서 전자구조를 정확히 기술하는 ab initio 방법은 큰 시스템에 적용하기 어려운 한계가 있다. 이러한 문제를 해결하기 위해 제로 차분 겹침(ZDO) 접근법에 기반한 반경험적 모델 해밀토니안, 예를 들어 Pariser‑Parr‑Pople, CNDO, INDO 등이 널리 사용된다. 본 논문에서는 CNDO/2와 INDO 방법을 이용해 분자계의 Hartree‑Fock 방정식을 풀어주는 Fortran 90 프로그램을 소개한다. INDO는 1주기 원소를 포함한 분자에, CNDO/2는 1·2주기 원소를 모두 포함한 분자에 적용 가능하다. 코드의 계산 효율성을 높여 C₆₀ 같은 대형 분자도 일반 PC에서 몇 초 만에 계산할 수 있다. 프로그램은 궤도와 총 에너지뿐 아니라 전기 쌍극자 모멘트, Mulliken 전하분포, 선형 광흡수 스펙트럼 등 다양한 물성을 제공한다. 또한 폴리머 단위셀당 에너지를 구하는 방법도 제시한다. 적용 사례로는 소규모 유기·무기 분자, 풀러렌 C₆₀, 그리고 교대 B‑N 사슬과 C 사슬 형태의 모델 폴리머가 있다.

상세 요약

이 논문이 제시하는 Fortran 90 기반 Hartree‑Fock 구현은 반경험적 ZDO 모델을 현대적인 고성능 컴퓨팅 환경에 맞게 재구성한 점에서 의미가 크다. 전통적으로 CNDO/2와 INDO는 1970‑80년대에 개발된 오래된 방법으로, 파라미터화된 두‑전자 적분을 사용해 계산량을 크게 줄이는 대신 정확도는 제한적이었다. 그러나 대규모 시스템, 특히 풀러렌(C₆₀)이나 폴리머와 같이 수백 개 이상의 원자를 포함하는 경우, 전통적인 ab initio 방법은 메모리와 CPU 시간의 폭발적 증가로 실용성이 떨어진다. 여기서 저자들은 Fortran 90의 모듈화와 배열 연산 최적화를 활용해, 핵심 연산인 Fock 행렬 구축과 SCF 반복을 효율적으로 구현하였다. 특히, 대칭성 이용, 스팸(스파스) 행렬 저장, 그리고 루프 전개 최소화를 통해 메모리 사용량을 크게 절감하고, 연산 속도를 “수 초” 수준으로 끌어올렸다.

프로그램이 제공하는 부가 기능—전기 쌍극자 모멘트, Mulliken 전하분포, 선형 광학 흡수 스펙트럼—은 전자 구조 계산 결과를 물리적·화학적 해석에 바로 연결시켜 준다. 이는 사용자가 별도의 후처리 코드를 작성할 필요 없이 한 번의 실행으로 다양한 물성을 얻을 수 있게 함으로써 연구 생산성을 크게 향상시킨다. 또한 폴리머 단위셀당 에너지를 계산하는 기능은 주기적 경계조건을 적용한 전자 구조 연구에 직접 활용될 수 있어, 전도성 고분자나 나노구조 물질 설계에 유용하다.

하지만 몇 가지 한계도 존재한다. 첫째, CNDO/2와 INDO는 파라미터가 실험 데이터에 맞춰져 있기 때문에, 새로운 원소 조합이나 비표준 결합 환경에서는 정확도가 떨어질 수 있다. 둘째, 현재 구현은 스핀-편극(SUHF)이나 다중참조(MR) 방법을 지원하지 않아, 강한 전자 상관이 중요한 시스템(예: 전이금속 복합체)에는 적용이 어려울 것이다. 셋째, Fortran 90은 현대 고성능 컴퓨팅에서 OpenMP나 MPI와 같은 병렬화 도구와의 연계가 비교적 제한적이므로, 다중 코어·다중 노드 환경에서의 확장성은 추가 개발이 필요하다.

그럼에도 불구하고, 이 프로그램은 교육용 및 초기 탐색 단계에서 매우 유용한 도구가 될 수 있다. 저비용 컴퓨터에서도 대형 분자와 폴리머를 빠르게 모델링할 수 있기 때문에, 실험적 설계와 이론적 검증 사이의 피드백 루프를 단축시킨다. 앞으로 파라미터 재조정, 병렬화, 그리고 전자 상관 보정(예: MP2, CI) 모듈을 추가한다면, 반경험적 모델의 적용 범위를 크게 확대할 수 있을 것으로 기대된다.