크라메르스-크로니히 관계

편집 역사 토론
전자기학 물리 계산물리 마지막 수정: 2026-08-20 04:12:20

1. 개요[편집]

크라메르스-크로니히 관계
Kramers-Kronig relations
약칭KK
분야선형 응답 이론 × 광학 × 신호처리
본질인과성 → 응답함수의 실수부·허수부가 서로를 결정
수학적 정체힐베르트 변환
예언자크로니히(1926) · 크라메르스(1927)

흡수만 재고도 굴절을 안다. 공짜처럼 보이지만, 값은 인과성이 이미 치렀다.

크라메르스-크로니히 관계(Kramers-Kronig relations, KK)는 어떤 선형 응답함수가 인과적(원인보다 결과가 먼저 올 수 없음)이기만 하면 그 함수의 실수부와 허수부가 서로 힐베르트 변환으로 묶여 하나가 나머지를 완전히 결정한다는 정리다. 물질의 흡수 스펙트럼(허수부)만 넓은 주파수 대역에서 재면, 굴절·분산(실수부)은 계산으로 복원된다. 반대도 성립한다. 두 양이 독립이 아니라 같은 인과적 함수의 두 얼굴이라는 것이 이 관계의 전부다.1

응답함수 χ(ω)=χ′(ω)+iχ′′(ω)\chi(\omega) = \chi'(\omega) + i\chi''(\omega) 를 생각하자. 시간영역 응답 χ(t)\chi(t) 가 인과적이라는 것(t<0t<0 에서 χ(t)=0\chi(t)=0)은 주파수영역에서 χ(ω)\chi(\omega) 가 상반평면 전체에서 해석적(analytic)이라는 조건과 같다. 여기에 코시 적분 정리와 유수 정리로 대표되는 복소해석의 표준 논법을 적용하면 실수부와 허수부를 잇는 적분식이 강제로 튀어나온다. 물리 법칙 하나 더 넣지 않고, 오직 “미래는 과거에 영향을 못 준다”만으로.

2. 인과성이 스펙트럼을 묶는다[편집]

선형·시간불변 시스템의 응답은 합성곱이다. 입력 f(t)f(t) 에 대한 출력이

x(t)=∫−∞∞χ(t−t′) f(t′) dt′x(t) = \int_{-\infty}^{\infty} \chi(t-t')\,f(t')\,dt'

인데, 원인이 결과를 앞설 수 없으므로 t′>tt'>t 인 미래 입력은 기여하면 안 된다. 즉 응답 커널이 χ(τ)=0 (τ<0)\chi(\tau)=0\ (\tau<0) 을 만족해야 한다. 이 한 줄이 핵심이다. 이 시간영역 인과성은 푸리에 변환을 거쳐 주파수영역의 해석성으로 번역된다.

이 조건을 푸리에 변환하면, χ(ω)\chi(\omega) 를 복소 ω\omega 로 확장했을 때 상반평면(Im ω>0\mathrm{Im}\,\omega>0)에서 극점이 없고 유계인 해석함수가 된다. e−iωte^{-i\omega t} 규약에서 상반평면은 e+Im(ω) te^{+\mathrm{Im}(\omega)\,t} 가 되어 t>0t>0 성분만 살아남기 때문이다. 상반평면에서 해석적이고 ∣ω∣→∞|\omega|\to\infty 에서 충분히 빨리 죽는 함수는 그 경계값(실축 위의 값)이 자기 자신에 대한 적분 제약을 받는다 — 이것이 KK 관계의 수학적 씨앗이다.2

3. 관계식[편집]

실축 위 한 점 ω\omega 를 상반평면 반원 경로로 우회하는 코시 적분에서, 극점 ω′=ω\omega'=\omega 를 코시 주요값(principal value, P\mathcal{P})으로 처리하면 다음이 나온다.

χ′(ω)=1π P ⁣∫−∞∞χ′′(ω′)ω′−ω dω′\chi'(\omega) = \frac{1}{\pi}\,\mathcal{P}\!\int_{-\infty}^{\infty}\frac{\chi''(\omega')}{\omega'-\omega}\,d\omega' χ′′(ω)=−1π P ⁣∫−∞∞χ′(ω′)ω′−ω dω′\chi''(\omega) = -\frac{1}{\pi}\,\mathcal{P}\!\int_{-\infty}^{\infty}\frac{\chi'(\omega')}{\omega'-\omega}\,d\omega'

즉 χ′\chi' 와 χ′′\chi'' 는 서로 힐베르트 변환 쌍이다(같은 이름을 쓰는 힐베르트 공간론과는 다른, 별개의 적분 변환이다). 물리적 응답함수는 시간영역에서 실수이므로 χ(−ω)=χ∗(ω)\chi(-\omega)=\chi^{*}(\omega), 곧 실수부는 우함수·허수부는 기함수다. 이 대칭을 넣어 양의 주파수만 남기면 실측에 바로 쓰는 형태가 된다.

χ′(ω)=2π P ⁣∫0∞ω′ χ′′(ω′)ω′2−ω2 dω′\chi'(\omega) = \frac{2}{\pi}\,\mathcal{P}\!\int_{0}^{\infty}\frac{\omega'\,\chi''(\omega')}{\omega'^{2}-\omega^{2}}\,d\omega'

허수부(흡수·소산)를 전 대역에서 알면 실수부(분산)가 적분 하나로 나온다는 뜻이다. 이 비대칭이 KK의 실용적 매력이다 — 흡수는 재기 쉽고, 위상·굴절은 재기 어렵다.

4. 유전율과 굴절률[편집]

전자기에서 가장 많이 쓰는 응답함수는 유전율 ε(ω)\varepsilon(\omega) 다. 진공 응답을 뺀 ε(ω)−1\varepsilon(\omega)-1 이 인과적 감수율 역할을 하므로 KK가 그대로 적용된다.

ε′(ω)−1=2π P ⁣∫0∞ω′ ε′′(ω′)ω′2−ω2 dω′\varepsilon'(\omega) - 1 = \frac{2}{\pi}\,\mathcal{P}\!\int_{0}^{\infty}\frac{\omega'\,\varepsilon''(\omega')}{\omega'^{2}-\omega^{2}}\,d\omega'

ε′′\varepsilon'' 은 흡수(빛이 물질에 먹히는 양)를, ε′\varepsilon' 은 굴절·분산을 지배한다. 복소 굴절률 n~=n+iκ=ε\tilde{n}=n+i\kappa=\sqrt{\varepsilon} 에서 κ\kappa 는 흡수계수, nn 은 위상속도를 정한다. 흡수선(공명) 근처에서 n(ω)n(\omega) 이 정상적으로 증가하다가 흡수 봉우리를 관통하며 거꾸로 꺾이는 비정상 분산(anomalous dispersion)은, KK 적분에서 좁은 ε′′\varepsilon'' 봉우리가 ε′\varepsilon' 에 만드는 S자 파형으로 정확히 설명된다. 굴절과 흡수가 붙어 다니는 것은 우연이 아니라 인과성의 명령이다.3

이 얽힘은 음의 굴절률 문서에서 본 “광대역·저손실 음굴절은 금지”라는 결론의 뿌리이기도 하다. 가파른 ε′\varepsilon' 변화는 KK를 통해 반드시 근처에 큰 ε′′\varepsilon''(흡수)를 동반하기 때문이다.

5. 합 규칙[편집]

KK 적분에 고주파 극한을 넣으면 합 규칙(sum rule)이 떨어진다. 모든 물질은 충분히 높은 주파수에서 전자가 자유입자처럼 굴어 ε(ω)→1−ωp2/ω2\varepsilon(\omega)\to 1-\omega_p^2/\omega^2 로 간다. 이 점근을 KK 관계와 맞추면 f-합 규칙이 나온다.

∫0∞ω ε′′(ω) dω=π2 ωp2,ωp2=ne2ε0m\int_{0}^{\infty}\omega\,\varepsilon''(\omega)\,d\omega = \frac{\pi}{2}\,\omega_p^{2}, \qquad \omega_p^{2}=\frac{n e^{2}}{\varepsilon_0 m}

흡수 스펙트럼을 진동수로 가중해 전 대역 적분하면 그 물질의 전자 밀도 nn 이 나온다는 뜻이다. 좌변은 실험으로 재고 우변은 원자 개수로 아는 값이라, 이 등식은 광학 데이터의 정합성 검사로 쓰인다. 재어 놓은 스펙트럼이 합 규칙을 어기면 어딘가 흡수를 놓쳤거나 눈금이 틀린 것이다. 양자역학의 슈뢰딩거 방정식에서 유도되는 토머스-라이헤-쿤(TRK) 합 규칙이 이것의 미시적 형제다. 플라스마 진동수 ωp\omega_p 와 드루드 모델의 연결 고리가 바로 여기다.

6. 실무: 반사율에서 위상 복원, 그리고 EELS[편집]

KK가 밥값을 하는 대표 무대는 측정 불가능한 위상을 복원하는 일이다.

  • 반사율 KK 분석. 정상 입사 반사율 R(ω)=∣r(ω)∣2R(\omega)=|r(\omega)|^2 는 세기라서 재기 쉽지만, 반사 계수의 위상 θ(ω)\theta(\omega) 는 직접 못 잰다. ln⁡r=12ln⁡R+iθ\ln r = \tfrac{1}{2}\ln R + i\theta 자체가 인과적 함수의 로그라 KK가 성립하므로, 넓은 대역의 RR 하나로 위상을 적분해 얻는다.
θ(ω)=−ωπ P ⁣∫0∞ln⁡R(ω′)ω′2−ω2 dω′\theta(\omega) = -\frac{\omega}{\pi}\,\mathcal{P}\!\int_{0}^{\infty}\frac{\ln R(\omega')}{\omega'^{2}-\omega^{2}}\,d\omega'

위상이 나오면 nn, κ\kappa, ε\varepsilon 이 전부 복원된다. 반사만으로 물질의 전 광학상수를 뽑는 이 절차는 고체 광학분광의 표준 작업이다.

  • 전자에너지손실분광(EELS). 시료를 통과한 전자의 에너지 손실 확률은 손실함수 Im[−1/ε(ω)]\mathrm{Im}[-1/\varepsilon(\omega)] 에 비례한다. 여기에 KK를 걸어 Re[−1/ε]\mathrm{Re}[-1/\varepsilon] 을 복원하면 ε(ω)\varepsilon(\omega) 전체와 유효 전자수·플라스몬 에너지가 나온다. 투과전자현미경 분광의 정량 분석이 이 위에 서 있다.
  • 광대역 물성 외삽·보간. 협대역에서만 잰 nn 과 κ\kappa 를 KK로 상호 검증하고, 물성 테이블의 빈 구간을 메운다.

7. 수치·실무의 함정[편집]

KK 적분은 원리는 우아하지만 실측에 걸면 골치가 여럿이다.

  • 꼬리가 지배한다. 적분은 00 부터 ∞\infty 인데 측정은 유한 대역이다. 잘린 고주파·저주파 꼬리를 물리적 점근형(고주파 ω−2\omega^{-2}, 저주파 상수)으로 외삽하지 않으면 실수부가 통째로 편향된다. 특히 관심 주파수가 대역 끝에 가까울수록 꼬리 오차가 지배적이다.
  • 주요값의 특이점. ω′=ω\omega'=\omega 에서 피적분함수가 발산하므로 코시 주요값을 수치적으로 다뤄야 한다. 특이점을 사이에 둔 대칭 구적이나 마우다치(Maclaurin) 방법을 쓰지 않고 순진하게 적분하면 값이 튄다.
  • 이산 데이터의 합 규칙 위반. 잡음과 눈금 오차로 f-합 규칙이 안 맞으면 그 자체가 데이터 품질 경고다. 무시하고 KK를 돌리면 없는 흡수선을 만들어내기도 한다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 크로니히(Ralph Kronig)가 1926년에 X선 분산·흡수 맥락에서, 크라메르스(Hendrik Kramers)가 1927년에 독립적으로 정식화했다. 둘 다 광학 분산에서 출발했지만, 정리 자체는 유전율·자기감수율·역학 임피던스·전달함수 등 인과적 선형 응답이면 무엇에나 적용된다. 응답함수의 이름만 바꿔 끼우면 그대로 동작한다. ↩

  2. 여기서 “충분히 빨리 죽는다”가 은근히 중요하다. ∣ω∣→∞|\omega|\to\infty 에서 χ(ω)\chi(\omega) 가 상수로 남으면(예: 즉각 응답 성분) 그 상수를 먼저 빼고 KK를 걸어야 한다. 유전율에서 ε−1\varepsilon-1 이나 ε−ε∞\varepsilon-\varepsilon_\infty 를 쓰는 이유가 이것이다. 안 빼면 적분이 안 죽어서 관계식이 발산한다. ↩

  3. “왜 프리즘은 파란빛을 더 꺾을까”라는 초등학교 질문의 정직한 답이 사실 KK다. 유리의 흡수선은 자외선에 있고, 그 흡수의 KK 그림자가 가시광 영역의 n(ω)n(\omega) 를 파장이 짧을수록 크게 만든다. 무지개의 순서는 자외선 흡수띠가 인과성을 통해 남긴 발자국인 셈이다. ↩