분자역학·DFT·파동함수·주기계·들뜬상태·동역학 방법을 질문·크기·정확도·시간척도에 맞게 선택합니다.

CHEMISTRY · PART 8 · 물리화학 III와 계산·분광07 / 07

MOF 한 구조의 최적화, 금속 노드의 스핀상태, 기체 흡착등온선, 반응장벽은 모두 ‘계산’이지만 같은 방법으로 풀 수 없습니다.

계산 비용이 가장 큰 방법이 아니라 연구 질문에 충분하고 검증 가능한 방법을 어떻게 고를까요?

이번 장의 답목표 관측량과 허용오차를 먼저 정하고 필요한 전자·원자 자유도와 샘플링을 결정합니다. 방법오차와 샘플링오차를 분리해 작은 기준계와 실험으로 검증합니다.

01기초부터 계산까지

CHEMISTRY, FROM ZERO

방법 선택 지도: 처음 배우는 사람을 위한 출발점

좋은 계산은 가장 비싼 방법이 아니라 질문에 필요한 자유도와 오차를 정확히 맞춘 계산입니다. 구조·상대에너지·반응속도·스펙트럼·흡착평형·확산 중 목표 관측량을 먼저 정하고, 모형오차·수치오차·샘플링오차·실험 매핑오차를 분리합니다.

START HERE

모형오차

선택한 물리·화학 근사 때문에 생기는 체계적 오차입니다.

START HERE

샘플링오차

유한한 구성·시간 표본으로 통계평균을 추정해 생기는 오차입니다.

START HERE

기준계산

방법의 정확도를 평가하기 위해 더 신뢰도 높은 이론이나 데이터와 비교하는 계산입니다.

CONCEPT ROUTE

이 글에서 생각이 이어지는 순서

용어를 외우는 대신, 관찰에서 모형과 계산을 거쳐 실제 해석으로 가는 흐름을 먼저 잡습니다.

01질문·관측량

출발점

02방법·모형

핵심 변수

03수렴·샘플링

모델 적용

04검증·불확도

검증과 해석

$$ right model, converged, validated $$

READ THE EQUATION

공식을 계산기에 넣기 전에 읽는 법

$$ \text{total error}\approx\text{model}+\text{numerical}+\text{sampling}+\text{mapping} $$
오차 항은 단순 합으로 독립적이지 않을 수 있지만 검증 계획을 세우는 실용적 분해입니다.
01

계 크기 N에 대한 비용 스케일링을 확인해 셀·기저를 늘릴 때 시간과 메모리가 얼마나 증가하는지 예산을 세웁니다.

02

cutoff·k점·기저·SCF 수렴은 수치오차를 줄이지만 잘못된 범함수·힘장·양성자화·스핀 같은 모형오차를 고치지 않습니다.

03

서로 독립적인 오차는 제곱합으로 1차 결합할 수 있지만 상관과 체계적 편향은 별도로 보고합니다.

WORKED FOUNDATIONS

기본 계산과 판단을 세 단계로 완전히 전개합니다

숫자 대입, 단위, 결과의 물리적 뜻을 한 번에 확인하세요.

01

O(N³) 방법에서 계 크기를 두 배

$$ \frac{C(2N)}{C(N)}=\frac{(2N)^3}{N^3}=8 $$

원자·기저함수 수가 두 배면 이상적인 계산량은 8배입니다. 메모리·통신·SCF 반복 때문에 실제 wall time은 다를 수 있습니다.

02

20, 10, 30 meV의 독립 오차 결합

$$ u_{total}=\sqrt{20^2+10^2+30^2}=37.4\ \mathrm{meV} $$

단순합 60 meV보다 작지만, 세 항이 같은 방향의 체계적 편향이면 제곱합은 과소평가합니다.

03

흡착 ΔG의 5.00 kJ mol⁻¹ 오차가 K에 미치는 영향

$$ \frac{K_{high}}{K_{low}}=e^{5000/(8.314\times298)}=7.52 $$

상온에서 몇 kJ mol⁻¹ 오차가 평형상수를 거의 한 자릿수 바꿉니다. 에너지 benchmark와 유한온도 샘플링을 함께 검증해야 합니다.

스스로 확인하기 · 계산 비용이 가장 큰 방법이 아니라 연구 질문에 충분하고 검증 가능한 방법을 어떻게 고를까요?

목표 관측량과 허용오차를 먼저 정하고 필요한 전자·원자 자유도와 샘플링을 결정합니다. 방법오차와 샘플링오차를 분리해 작은 기준계와 실험으로 검증합니다.

02왜 이 개념이 필요한가

MOF 한 구조의 최적화, 금속 노드의 스핀상태, 기체 흡착등온선, 반응장벽은 모두 ‘계산’이지만 같은 방법으로 풀 수 없습니다. 이 장면을 풀려면 용어를 외우는 것보다 어떤 물리량을 고정하고 무엇을 변화시키는지 먼저 정해야 합니다. 분자역학·DFT·파동함수·주기계·들뜬상태·동역학 방법을 질문·크기·정확도·시간척도에 맞게 선택합니다.

읽기 전에 구분할 것
  • 전자에너지 정확도와 유한온도 평균의 샘플링을 구분합니다.
  • 큰 기저·촘촘한 k점은 수치오차를 줄이지만 모형 편향을 없애지 않습니다.
  • 계산과 실험은 같은 조성·온도·압력·관측량을 비교해야 합니다.

03핵심 모델과 식

힘장은 큰 계·긴 시간에 유리하지만 결합 변화가 어렵고, DFT는 전자구조와 반응을 다루지만 범함수 의존성이 있습니다. post-HF·다중참조는 작은 분자의 높은 정확도에, GCMC는 흡착 평형, MD는 시간 상관과 수송에 적합합니다.

$$ \text{total error}\approx\text{model}+\text{numerical}+\text{sampling}+\text{mapping} $$
오차 항은 단순 합으로 독립적이지 않을 수 있지만 검증 계획을 세우는 실용적 분해입니다.

04공학 예제로 계산하기

CO₂@MOF는 DFT로 결합자리와 전하밀도를, 고품질 양자 기준으로 힘장 매개변수를, GCMC로 등온선을, MD로 확산을 구할 수 있습니다. 각 단계의 출력과 불확도를 다음 단계에 전달합니다.

풀이와 해석의 순서
  • 연구 질문을 구조·에너지·속도·스펙트럼·통계평균 중 하나의 관측량으로 명확히 씁니다.
  • 최소 필요한 계 크기·경계조건·전자구조·샘플링 방법을 선택합니다.
  • 수렴검사·대안 방법·기준 데이터로 오차 예산을 만들고 재현 정보를 기록합니다.

05오해·한계·다음 연결

방법 이름만으로 품질이 보장되지 않습니다. 잘못된 양성자 위치·산화상태·스핀·용매·결함 모형은 높은 수준 이론보다 더 큰 오차를 냅니다.

재료·DFT로 이어지는 연결
  • 이 파트까지의 화학 개념이 DFT 입력과 출력 해석의 사전 지식입니다.
  • 멀티스케일 모델은 전자→원자→미세구조→공정 사이 매핑을 명시해야 합니다.
  • 기계학습 힘장은 학습영역과 불확도 감시가 핵심입니다.

06이 글의 용어 사전

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

모형오차model error
선택한 물리·화학 근사 때문에 생기는 체계적 오차입니다.
샘플링오차sampling error
유한한 구성·시간 표본으로 통계평균을 추정해 생기는 오차입니다.
기준계산benchmark calculation
방법의 정확도를 평가하기 위해 더 신뢰도 높은 이론이나 데이터와 비교하는 계산입니다.
속도상수rate constant
주어진 온도에서 속도식의 비례를 정하는 반응 고유의 계수입니다.
흡착adsorption
기체나 용액의 성분이 고체 표면 또는 계면에 농축되는 현상입니다.
과전압overpotential
전극 반응을 유한한 속도로 진행시키기 위해 평형전위보다 추가로 필요한 전위입니다.

07참고자료와 더 읽을거리

  1. KAIST 화학과 학부 교과목 개요물리·유기·분석·무기·재료화학 과목의 공식 범위와 연결을 확인하는 기준입니다.
  2. NIST Chemistry WebBook반응속도·분광·열화학 데이터를 계산과 실험의 기준값으로 사용합니다.
  3. IUPAC Gold Book: electrochemistry and kinetics속도상수·과전압·전극전위·흡착 용어를 표준 정의로 확인합니다.
NEXT CHEMISTRY NOTE고급 양자화학: 공리와 근사