양자역학

1. 개요[편집]

양자역학 Quantum Mechanics
지배 방정식	슈뢰딩거 방정식
상태 표현	파동함수 ψ (복소 진폭), 확률밀도 |ψ|²
적용 스케일	원자·분자·고체 전자 구조 (~0.1 nm, ~eV)
계산 갈래	파동함수 기반 / 밀도범함수이론 / 양자 몬테카를로

방정식은 선형이고 우아한데, 미지수가 전자 수에 지수적으로 늘어난다. 그래서 이 바닥은 근사의 예술이 됐다.

양자역학(quantum mechanics)은 원자 규모 이하의 물질과 빛의 거동을 확률 진폭인 파동함수로 기술하는 물리 이론이다. 고전역학이 입자의 위치와 운동량을 동시에 확정된 값으로 다루는 반면, 양자역학은 상태를 힐베르트 공간의 벡터로 두고 관측 가능량을 연산자의 고윳값으로 얻는다. 관측 결과가 확률적으로만 예측된다는 점, 그리고 위치와 운동량처럼 교환하지 않는 연산자 쌍에 불확정성 관계가 성립한다는 점이 이론의 근본 성격을 결정한다.¹

이 위키에서 양자역학이 중요한 이유는 철학이 아니라 계산이다. 재료의 밴드갭, 분자의 결합 에너지, 촉매 표면의 반응 장벽 — 실험 없이 이런 값을 내놓는 유일한 원리적 경로가 양자역학이고, 그것을 컴퓨터로 수행하는 것이 계산물리와 양자화학의 본업이다.

2. 슈뢰딩거 방정식[편집]

시간 의존 형태는 다음과 같다.

$$i\hbar \frac{\partial}{\partial t} \Psi(\mathbf{r}, t) = \hat{H} \Psi(\mathbf{r}, t)$$

여기서 해밀토니안 $\hat{H}$는 운동 에너지와 퍼텐셜 에너지 연산자의 합이다. 퍼텐셜이 시간에 무관하면 변수분리로 시간 독립 형태

$$\hat{H}\psi = E\psi$$

가 되고, 이 순간 물리 문제는 정확히 **고유값 문제**로 변신한다. 수치해석 입장에서는 대단히 반가운 소식이다. 이미 잘 발달한 선형대수 기계를 그대로 갖다 쓸 수 있다는 뜻이기 때문. 실제로 전자 구조 계산 코드의 시간 대부분은 거대한 에르미트 행렬의 낮은 고유값 몇 개를 구하는 데 쓰이며, 여기에 데이비드슨 알고리즘 같은 반복 대각화 기법이 동원된다.

반가움은 오래가지 않는다. $N$개 전자의 파동함수는 $3N$차원 공간의 함수라서, 격자점을 각 축에 $M$개씩만 잡아도 미지수가 $M^{3N}$개다. 전자 10개짜리 분자도 저장이 불가능하다. 이 지수 폭발이 양자 계산 전체의 설계 제약이며, 월터 콘이 “지수 장벽(exponential wall)“이라 부른 바로 그것이다.

3. 근사의 계층[편집]

그래서 실제 계산은 근사를 층층이 쌓아 올린다.

• 본-오펜하이머 근사 — 핵이 전자보다 수천 배 무거우니 핵을 정지시킨 채 전자만 푼다. 이 근사가 포텐셜 에너지 표면이라는 개념 자체를 만들어내고, 분자동역학이 성립할 토대가 된다.
• 하트리-폭 방법 — 다전자 파동함수를 단일 슬레이터 행렬식으로 두고, 각 전자가 나머지의 평균장을 느낀다고 본다. 교환은 정확히, 상관은 아예 무시한다. 변분법이 여기서 이론적 보증을 준다.
• 포스트-HF — 빠진 상관 에너지를 되찾는 계열. 배치상호작용과 결합 클러스터 방법이 대표이며, CCSD(T)는 오늘날 “황금 표준”으로 불린다. 정확도의 대가는 계 크기에 대한 가파른 비용 증가다.²
• 밀도범함수이론(DFT) — 미지수를 $3N$차원 파동함수에서 3차원 전자밀도로 갈아탄다. 호엔베르크-콘 정리가 이 교체를 정당화하고, 실제 계산은 콘-샴 방정식이라는 유효 단일입자 방정식으로 수행된다. 모든 어려움은 교환-상관 범함수 하나에 몰아넣었는데, 그 정확한 형태는 아무도 모른다.
• 양자 몬테카를로 — 고차원 적분을 확률 표본추출로 정면 돌파한다. 페르미온 부호 문제라는 고유의 벽이 있지만, 정확도 기준점으로서 가치가 크다.

4. 고체와 주기계[편집]

무한히 반복되는 결정에서는 블로흐 정리가 파동함수를 주기 함수와 평면파의 곱으로 제한한다. 덕분에 무한 결정이 단위 격자 하나의 문제로 축소되고, 계산은 역격자 공간의 $\mathbf{k}$점 표본추출로 바뀐다. 이 표본을 효율적으로 고르는 표준이 몽코스트-팩 격자다. 결과로 얻는 것이 밴드 구조 계산과 상태밀도이며, 여기서 금속·반도체·절연체가 갈린다.

기저함수로는 평면파 기저가 널리 쓰이는데, 핵 근처에서 파동함수가 격렬하게 진동하는 문제 때문에 그대로는 절단 에너지가 감당이 안 된다. 그래서 내각 전자를 유사퍼텐셜로 갈음하거나 PAW 방법으로 재구성한다. 국소화된 관점이 필요하면 블로흐 상태를 와니에 함수로 변환한다. 원자 궤도의 각 부분이 구면조화 함수라는 사실도 기저함수 설계의 출발점이다.

5. 시간 발전과 터널링[편집]

위 화면이 보여주는 것이 슈뢰딩거 시간 발전의 성격이다. 발전 자체는 결정론적이고 유니터리하다 — 규격화 $\int|\psi|^2$는 처음부터 끝까지 1이고, 흡수층이 먹은 몫까지 장부에 넣으면 한 자리도 새지 않는다. 확률이 등장하는 지점은 측정이며, 확률밀도 $|\psi|^2$는 진폭의 제곱이라 경로들이 서로 상쇄될 수 있다. 반면 고전 확률 모형에서는 경우의 수가 더해지기만 하고 상쇄되지 않는다. 이 “간섭이 있느냐”의 차이가 양자와 고전을 가르는 실질적 경계이며, 이중 슬릿 실험이 바로 그 경계를 보여주는 장치다.³

터널링은 이 그림에서 파생되는 가장 실용적인 결론이기도 하다. 에너지가 장벽보다 낮아도 파동함수는 장벽 내부에서 $e^{-\kappa x}$로 감쇠할 뿐 0이 되지 않으므로, 반대편으로 새어 나오는 진폭이 남는다. 감쇠 상수 $\kappa = \sqrt{2m(V_0-E)}/\hbar$가 지수에 들어가 있어서 장벽 폭이 조금만 늘어도 투과확률이 몇 자릿수씩 떨어진다. 주사터널링현미경의 원자 분해능, 알파 붕괴의 반감기 편차, 플래시 메모리의 게이트 누설이 모두 이 지수 의존성의 결과다.

6. 계산 도구와 확장[편집]

실무 코드로는 VASP, Quantum ESPRESSO가 고체 쪽에서, 양자화학 패키지들이 분자 쪽에서 쓰인다. 핵을 고전적으로 움직이면서 전자를 매 스텝 푸는 것이 제일원리 분자동역학이고, 파동함수를 가상 자유도로 함께 굴리는 카-파리넬로 분자동역학이 그 고전적 구현이다. 격자 진동을 다루면 포논 분산이 나오고, 여기서 열용량과 열전도가 따라 나온다.

7. 여담[편집]

양자 계산의 논문에는 “chemical accuracy(1 kcal/mol)“라는 목표치가 자주 등장한다. 전체 에너지의 절댓값은 수천 kcal/mol 규모인데 그 안에서 1 kcal/mol을 맞추자는 것이니, 상대 오차로 따지면 $10^{-4}$ 수준을 요구하는 셈이다. 그것도 정확한 범함수를 모르는 채로. 그래서 이 분야의 실전 감각은 “절대값은 못 믿어도 차이값은 믿는다”에 가깝고, 오차 상쇄를 전략적으로 설계하는 것이 실력의 상당 부분을 차지한다.⁴

8. 관련 문서[편집]

• 슈뢰딩거 방정식 · 본-오펜하이머 근사
• 하트리-폭 방법 · 배치상호작용 · 결합 클러스터 방법
• 밀도범함수이론 · 콘-샴 방정식 · 교환-상관 범함수
• 양자 몬테카를로 · 고유값 문제
• 평면파 기저 · 유사퍼텐셜 · PAW 방법
• 밴드 구조 계산 · 상태밀도 · 포논 분산
• 구면조화 함수 · 와니에 함수
• 계산물리 · 해밀토니안 역학 · 분자동역학

9. Footnotes[편집]

1. 해석 논쟁(코펜하겐·다세계·파일럿 파동 등)은 지금도 진행 중이지만, 계산 결과는 해석과 무관하게 동일하다. 그래서 실무자들 사이에는 “닥치고 계산하라(shut up and calculate)“는 유명한 태도가 있다. 철학은 저녁 자리에서 하고 낮에는 대각화를 돌린다는 뜻. ↩

2. CCSD(T)의 명목 비용은 계 크기의 7제곱에 비례한다. 원자 수를 2배로 늘리면 계산량이 128배가 된다는 소리다. 황금 표준이 황금인 이유는 정확해서이기도 하지만 비싸서이기도 하다. ↩

3. 파인만이 이중 슬릿을 두고 “양자역학의 유일한 미스터리”라고 부른 것은 이 상쇄 때문이다. 확률이 아니라 진폭이 더해진다는 한 줄이 고전과 양자의 전부라 해도 크게 틀리지 않는다. ↩

4. 그래서 논문 심사에서 “다른 범함수로도 해보셨나요?”는 사실상 정형화된 질문이다. 여러 범함수의 결과가 흩어지면 그 물리량은 범함수 의존성이 크다는 뜻이고, 그때부터는 숫자보다 경향을 논하게 된다. ↩