흑체복사의 사례인 빛을 내뿜는 용암 |
[clearfix]
1. 개요
흑체복사(黑體輻射, black body radiation)는 흑체에서 온도에 따라 빛이 나오는 현상이다. 방출 에너지량은 절대온도의 네제곱에 비례하며, 평균 에너지는 절대 온도에 비례한다. 세기가 극대가 되는 파장은 절대 온도에 반비례하고 [math(c = \lambda\nu)]에서 [math(\nu = \dfrac c\lambda)]이므로 진동수의 경우엔 절대온도에 비례한다는 결론이 얻어지지만 구체적인 값을 구할 때에는 저 공식을 쓰면 안 된다. 흑체 복사의 특성은 플랑크 법칙을 따르는 것으로 알려져 있다.2. 열 복사
일반적으로 (열)복사(熱輻射)는 물리학에서 물체로부터 열이나 전자기파가 사방팔방으로 완전히 (흡수 또는) 방출되는 것을 가리키는 용어이다. 이를 전제로 어떤 임의의 가상의 물체를 상정하고 빛을 방출하지 않는 표면을 가지면서 한 줄기 빛만을 통과시킬 수 있고 그 내부는 반사작용이 여전히 있는 공동(cavity) 터널(hole)을 가정하는 이상적인 흑체를 모델로 구현해볼 수 있다.
3. 배경
19세기 말 독일 제국은 알자스-로렌 지방에서 제철 공업을 일으켰다. 강철을 만드는 과정에서 철광석을 코크스 등으로 연소시켜서 고온으로 만들 때 그 빛깔로 온도를 파악하려 했다.붉은 색보다 파란색일 때 온도가 더 높다는 것은 경험적으로는 알았지만, 빛깔과 온도와의 관계에 대해 이론 설명이 필요해졌다. 쇳물 만드는데 그냥 녹을 때까지 퍼부으면 안되나 싶겠지만, 온도가 낮으면 철광석의 산화철에서 산소가 빠져나오지 못해 저질의 철이 되며, 너무 높으면 고온에서 폭발적으로 반응하는 산소의 특성상 다시 산화철이 되기 때문에 철 원자에서 산소만 뽑아내는 매우 정밀한 온도 조절이 필요하다. 그리고 강철이 굳은 뒤에 다시 가열하여 원하는 물성을 얻는 열처리 과정에서도 정밀한 온도 조절이 필요하다.[1] 즉, 조금이라도 더 세밀한 온도 조절을 위해 흑체 복사 이론이 태동했다는 것이다.
이러한 관계에 대한 연구를 하던 독일 하이델베르크 대학의 물리학자 구스타프 키르히호프는 흑체라는 이상적인 물체를 가정하고, 흑체가 방출하는 복사의 세기가 오로지 흑체의 온도와 빛의 파장, 두 요인에 의해서만 결정된다는 것을 밝혀냈다. 이로부터 주어진 온도에서 모든 파장 범위에 해당하는 흑체복사 스펙트럼 분포를 실험을 통해 알아내고, 또한 그 결과를 공식으로부터 알아내는 것이 한동안 당대 물리학자들을 괴롭힌 ‘흑체복사 문제’로 남게 되었으며, 아래 나열되는 복사법칙들이 그 문제에 대한 일련의 해답이라고 볼 수 있다.
3.1. 용어 정리
'세기'라고 하더라도 에너지 그 자체를 가리키는 경우가 있는가 하면, 단위 공간에 가해지는 에너지(에너지 밀도)를 쓰는 경우도 있고, 단위 시간당 가해지는 에너지, 즉 일률을 다루면 \'복사속'이라고 하며, 복사속은 다시 단위 면적에 가하는 양을 가리키는 것인지(복사조도), 점광원으로 환산된 물리량인지(복사강도), 점광원으로 환산된 단위 면적 당 복사속인지(복사휘도), 분석하고자 하는 대상이 복사원이냐 아니면 피복사체가 반사하거나 재복사한 빛이냐에 따라서(복사발산도 / 복사조도)[2] 또 다른 용어를 쓴다. 과거에는 이러한 개념이 세분화되지 않은 상태로 관습적으로 ' 강도'(영어: intensity, 독어: Intensität)라는 말로 통일되다시피 해서, 후술할 역사적인 인물들이 출판한 논문에서 쓰이는 용어를 읽을 때 주의를 요하며 문맥에 따라 구분해서 파악해야한다.물리학적으로는 복사로 방출되는 것도 결국 진동수와 파장을 갖는 전자기파이기 때문에, 위의 물리량을 진동수, 혹은 파장으로 미분해서 얻어지는 분포 함수를 다루기도 하는데 이 경우 앞에 '분광-'(spectral)이라는 용어를 덧붙인다. 이를테면 '분광 복사휘도'는 수학적으로 '복사휘도 분포 함수'와 의미가 같다. 또한 칸델라가 엮인 광도와 관련된 물리량과 구분하기 위해 측광 분야에서 쓰는 물리량과 같은 기호를 쓰는 경우 [math(\rm e)]를 아래첨자로 덧붙여서 구분하기도 한다. 물론 열역학, 통계역학, 양자역학 등에서는 독립적으로 전혀 다른 기호를 쓰는 경우가 일반적이긴 하다.
리스트로 나열하면 아래와 같이 엄청나게 많긴 한데, 본 문서에서는 에너지, 에너지 밀도, 복사휘도를 중심으로 다룰 것이다.
물리량 | 기호 | 차원 | 단위 | 의미 |
(복사) 에너지 Radiant energy |
[math(Q_{\rm e} \\ W \\ E)] | [math(\sf ML^2T^{-2})] | [math(\rm J)] | 전자기파 복사의 에너지 |
(복사) 에너지 밀도 Radiant energy density |
[math(w_{\rm e} \\ u)] | [math(\sf ML^{-1}T^{-2})] | [math(\rm J/m^3)] | 단위 부피당 가해지는 전자기파 복사의 에너지 |
분광 (복사) 에너지 밀도 에너지 밀도 분포 함수 spectral radiant energy density |
[math(u_\nu(\nu,\,T) = \dfrac{\partial w_{\rm e}}{\partial\nu})] | [math(\sf ML^{-1}T^{-1})] | [math(\rm J/(Hz{\cdot}m^3))] | 전자기파의 스펙트럼 중 특정 진동수가 단위 부피에 가하는 복사 에너지 |
[math(u_\lambda(\lambda,\,T) = \dfrac{\partial w_{\rm e}}{\partial\lambda})] | [math(\sf ML^{-2}T^{-2})] | [math(\rm J/(nm{\cdot}m^3))] | 전자기파의 스펙트럼 중 특정 파장이 단위 부피에 가하는 복사 에너지 | |
복사속 Radiant flux |
[math(\varPhi_{\rm e} \\ P \\ F)] | [math(\sf ML^2T^{-3})] | [math({\rm W} = {\rm J/s})] | 단위 시간당 방출, 반사, 투과, 흡수된 복사 에너지 |
분광 복사속 spectral radiant flux |
[math(\varPhi_{{\rm e},\,\nu} = \dfrac{\partial\varPhi_{\rm e}}{\partial\nu})] | [math(\sf ML^2T^{-2})] | [math(\rm W/Hz)] | 전자기파의 스펙트럼 중 특정 진동수의 복사속 |
[math(\varPhi_{{\rm e},\,\lambda} = \dfrac{\partial\varPhi_{\rm e}}{\partial\lambda})] | [math(\sf MLT^{-3})] | [math(\rm W/nm)] | 전자기파의 스펙트럼 중 특정 파장의 복사속 | |
복사강도 Radiant intensity |
[math(I_{\rm e,\,\Omega})] | [math(\sf ML^2T^{-3})] | [math(\rm W/sr)] | 점광원 단위로 환산된[3] 복사속 |
분광 복사강도 spectral radiance specific intensity |
[math(I_{{\rm e,\,\Omega},\,\nu} = \dfrac{\partial I_{\rm e,\,\Omega}}{\partial\nu})] | [math(\sf ML^2T^{-2})] | [math(\rm W/(Hz{\cdot}sr))] | 점광원 단위로 환산된 전자기파 스펙트럼 중 특정 진동수의 복사 에너지 |
[math(I_{{\rm e,\,\Omega},\,\lambda} = \dfrac{\partial I_{\rm e,\,\Omega}}{\partial\lambda})] | [math(\sf MLT^{-3})] | [math(\rm W/(nm{\cdot}sr))] | 점광원 단위로 환산된 전자기파 스펙트럼 중 특정 파장의 복사 에너지 | |
복사휘도 radiance |
[math(L_{\rm e,\,\Omega}\\ B)] | [math(\sf MT^{-3})] | [math(\rm W/(sr{\cdot}m^2))] | 단위 면적에 가해지는 복사강도 |
분광 복사휘도 복사휘도 분포 함수 spectral radiance |
[math(L_{{\rm e,\,\Omega},\,\nu} = B_\nu(\nu,\,T) = \dfrac{\partial B}{\partial\nu})] | [math(\sf MT^{-2})] | [math(\rm W/(Hz{\cdot}sr{\cdot}m^2))] | 단위 면적에 가해지는 전자기파 스펙트럼 중 특정 진동수의 복사강도 |
[math(L_{{\rm e,\,\Omega},\,\lambda} = B_\lambda(\lambda,\,T) = \dfrac{\partial B}{\partial\lambda})] | [math(\sf ML^{-1}T^{-3})] | [math(\rm W/(nm{\cdot}sr{\cdot}m^2))] | 단위 면적에 가해지는 전자기파 스펙트럼 중 특정 파장의 복사강도 | |
복사조도 irradiance |
[math(E_{\rm e})] | [math(\sf MT^{-3})] | [math(\rm W/m^2)] | 어떤 물체가 받은 단위 면적당의 복사속 |
분광 복사조도 spectral irradiance |
[math(E_{{\rm e},\,\nu} = \dfrac{\partial E_{\rm e}}{\partial\nu})] | [math(\sf MT^{-2})] | [math(\rm W/(Hz{\cdot}m^2))] | 어떤 물체가 받은 전자기파 스펙트럼 중 특정 진동수의 단위 면적당 복사속 |
[math(E_{{\rm e},\,\lambda} = \dfrac{\partial E_{\rm e}}{\partial\lambda})] | [math(\sf ML^{-1}T^{-3})] | [math(\rm W/(nm{\cdot}m^2))] | 어떤 물체가 받은 전자기파 스펙트럼 중 특정 파장의 단위 면적당 복사속 | |
복사발산도 radiant exitance |
[math(M_{\rm e})] | [math(\sf MT^{-3})] | [math(\rm W/m^2)] | 어떤 복사원에서 방출되는 단위 면적당의 복사속 |
분광 복사발산도 spectral exitance |
[math(M_{{\rm e},\,\nu} = \dfrac{\partial M_{\rm e}}{\partial\nu})] | [math(\sf MT^{-2})] | [math(\rm W/(Hz{\cdot}m^2))] | 어떤 복사원에서 방출되는 전자기파 스펙트럼 중 특정 진동수의 단위 면적당 복사속 |
[math(M_{{\rm e},\,\lambda} = \dfrac{\partial M_{\rm e}}{\partial\lambda})] | [math(\sf ML^{-1}T^{-3})] | [math(\rm W/(nm{\cdot}m^2))] | 어떤 복사원에서 방출되는 전자기파 스펙트럼 중 특정 파장의 단위 면적당 복사속 | |
복사도 radiosity |
[math(J_{\rm e})] | [math(\sf MT^{-3})] | [math(\rm W/m^2)] | 어떤 복사원에서 방출, 반사, 투과되는 단위 면적당의 복사속[4] |
분광 복사도 spectral radiosity |
[math(J_{{\rm e},\,\nu} = \dfrac{\partial J_{\rm e}}{\partial\nu})] | [math(\sf MT^{-2})] | [math(\rm W/(Hz{\cdot}m^2))] | 어떤 복사원에서 방출, 반사 투과되는 전자기파 스펙트럼 중 특정 진동수의 단위 면적당 복사속 |
[math(J_{{\rm e},\,\lambda} = \dfrac{\partial J_{\rm e}}{\partial\lambda})] | [math(\sf ML^{-1}T^{-3})] | [math(\rm W/(nm{\cdot}m^2))] | 어떤 복사원에서 방출, 반사, 투과되는 전자기파 스펙트럼 중 특정 파장의 단위 면적당 복사속 | |
복사 노출도 radiant exposure |
[math(H_{\rm e})] | [math(\sf MT^{-2})] | [math(\rm J/m^2)] | 어떤 물체가 받은 단위 면적당의 복사 에너지. 복사조도를 시간으로 적분한 물리량. |
분광 복사 노출도 spectral exposure |
[math(H_{{\rm e},\,\nu} = \dfrac{\partial H_{\rm e}}{\partial\nu})] | [math(\sf MT^{-1})] | [math(\rm J/(Hz{\cdot}m^2))] | 어떤 물체가 받은 전자기파 스펙트럼 중 특정 진동수의 단위 면적당 에너지 |
[math(H_{{\rm e},\,\lambda} = \dfrac{\partial H_{\rm e}}{\partial\lambda})] | [math(\sf ML^{-1}T^{-2})] | [math(\rm J/(nm{\cdot}m^2))] | 어떤 물체가 받은 전자기파 스펙트럼 중 특정 파장의 단위 면적당 복사 에너지 |
순수하게 광원에서 나오는 빛의 세기에 대해서만 다루고 인간의 시각 수용체가 느끼는 정도를 반영하지 않기 때문에 물리량의 차원에 [math(\sf J)]가 없고, 따라서 단위에 [math(rm cd)]도 들어가지 않는 것이 특기사항이다. 조도, 휘도 등 측광 분야와 용어가 겹치는 부분이 있으나 전혀 다른 물리량임에 주의.
4. 슈테판-볼츠만 법칙
1864년에 존 틴들(John Tyndall, 1820~1893) 발표한 《보이는 복사와 보이지 않는 복사에 대하여》[5]의 실험 데이터를 바탕으로, 1879년에 오스트리아령 슬로베니아 출신의 요제프 슈테판(1835~1893)이 열복사의 세기가 절대온도의 네제곱에 비례한다는 것을 귀납적으로 발견해서 보고하였고 이후 오스트리아 제국의 루트비히 볼츠만이 1884년에 발표한 《빛의 전자기학 이론을 바탕으로 한 슈테판 법칙, 열복사의 온도 의존성 유도》를 통하여 이것이 참임을 증명하였다. 이 법칙이 슈테판-볼츠만 법칙(Stefan-Boltzmann law)이다.볼츠만은 복사압 [math(P)]와 내부 에너지 밀도 [math(u)]가 [math(P = \dfrac13u)]를 만족한다는 관계를 이용해서 유도해냈는데, 간단한 열역학적인 관계식을 이용하면 어렵지 않게 유도할 수 있다.
먼저 내부 에너지 [math(U)]는 절대온도 [math(T)], 엔트로피 [math(S)], 압력 [math(P)], 부피 [math(V)]에 대하여
[math({\rm d}U = T{\rm\,d}S - P{\rm\,d}V)] |
[math({\left(\dfrac{\partial U}{\partial V}\right)}_T = T{\left(\dfrac{\partial S}{\partial V}\right)}_T-P)] |
[math(\begin{aligned} u &= T{\left(\dfrac{\partial P}{\partial T}\right)}_V - P \\ &= T{\left(\dfrac13\dfrac{\partial u}{\partial T}\right)}_V - \dfrac13u \\ \therefore 4u &= T{\left(\dfrac{\partial u}{\partial T}\right)}_V\end{aligned})] |
[math(\dfrac{4{\rm\,d}T}T = \dfrac{{\rm d}u}u)] |
5. 빈 법칙
빈 법칙(Wien's law)은 1897년에 독일의 물리학자 빌헬름 빈(Wilhelm Wien, 1864~1928)[6]이 제안한 법칙으로 빈 분포 법칙(Wien's distribution law) 혹은 빈 근사(Wien's approximation)라고도 불린다.빈은 흑체의 분광 복사휘도(spectral radiance; 또는 복사휘도 분포 함수)[7] [math(B_\lambda(\lambda,\,T))]가 두 상수 [math(c_1)], [math(c_2)]를 매개로 하여 다음과 같은 꼴로 주어진다고 주장하였다.
[math(B_\lambda(\lambda,\,T) = \dfrac{c_1}{\lambda^5}e^{-\frac{c_2}{\lambda T}})] |
[math(\displaystyle \int_{0{\rm\,m}}^\infty B_\lambda(\lambda,\,T){\rm\,d}\lambda = \int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,d}\nu)] |
[math(\begin{aligned} \int_{0{\rm\,m}}^\infty B_\lambda(\lambda,\,T){\rm\,d}\lambda &= \int_\infty^{0{\rm\,Hz}}B_\lambda(\lambda(\nu),\,T){\left(-\frac c{\nu^2}\right)}{\rm\,d}\nu \\ &= \int_{0{\rm\,Hz}}^\infty \frac c{\nu^2}B_\lambda(\lambda(\nu),\,T){\rm\,d}\nu \end{aligned})] |
[math(B_\nu(\nu,\,T) = \dfrac{c_1\nu^3}{c^4}e^{-\frac{c_2\nu}{cT}})] |
실제로 적분해보면
[math(\begin{aligned} \int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,d}\nu &= \int_{0{\rm\,Hz}}^\infty \frac{c_1\nu^3}{c^4}e^{-\frac{c_2\nu}{cT}}{\rm\,d}\nu \\ &= {\left[-\frac{c_1T\nu^3}{c_2c^3}e^{-\frac{c_2\nu}{cT}}-\frac{3c_1T^2\nu^2}{{c_2}^2c^2}e^{-\frac{c_2\nu}{cT}}-\frac{6c_1T^3\nu}{{c_2}^3c}e^{-\frac{c_2\nu}{cT}}-\frac{6c_1T^4}{{c_2}^4}e^{-\frac{c_2\nu}{cT}}\right]}_{0{\rm\,Hz}}^\infty \\ &= \frac{6c_1}{{c_2}^4}T^4\end{aligned})] |
5.1. 빈 변위 법칙
빈 법칙은 일정 온도 조건에서 흑체로부터 방사되는 빛의 복사휘도가 가장 강한[8] 파장 혹은 진동수을 구하는 데에도 유용하게 쓰였다. 이 값은 곧 분광 복사휘도 [math(B_\nu(\nu,\,T))]의 극댓값이 나타나는 조건을 찾는 것과 같으므로, 예를 들어 [math(B_\lambda(\lambda,\,T))]를 [math(\lambda)]로 편미분해서 그 식이 [math(0{\rm\,W/(m^4sr)})]이 되는 조건을 찾아보면[math(\begin{aligned}\frac{\partial}{\partial\lambda}B_\lambda(\lambda,\,T) &= \frac{c_1}{\lambda^6}e^{-\frac{c^2}{\lambda T}}{\left(\frac{c_2}{\lambda T}-5\right)} = 0{\rm\,W/(m^4sr)} \\ \therefore\frac{c_2}{\lambda_{\sf max} T} &= 5 \\ \lambda_{\sf max} &= \frac{c_2}{5T}\end{aligned})] |
[math(\lambda_{\sf max} \approx \dfrac{2.877\,554{\rm\,mm{\cdot}K}}T)] |
파장에 따른 에너지 밀도를 나타낸 흑체복사곡선(플랑크 곡선) |
이 도표를 보면 온도가 낮으면 파장이 긴 적외선이나 붉은 빛이 많이 나오고 온도가 높으면 파장이 짧은 파란색이나 자외선이 많이 나오는 걸 알 수 있다. 더욱 고온이 되면 자외선 영역도 넘어 X선이나 감마선을 내게되고 반대로 온도가 낮으면 마이크로파, 밀리미터파 등 전파의 영역의 최대 파장이 나온다. 우주배경복사는 저온에서 이런 전파가 나오는 것. 태양의 표면온도는 6000도 정도이고 파장으론 455~500 nm 정도가 가장 강한 피크를 이룬다. 피크의 절반의 세기가 되는 파장은 각각 350 nm 와 840 nm 정도이다.
이때, 특정 파장을 내기 위한 필요한 온도를 색온도(色溫度)라고 한다.
이 법칙은 고온 물체의 색으로부터 그 온도를 추정하는 실용적인 목적으로 매우 많이 쓰인다. 달아오른 금속이나 불꽃의 온도 측정, 태양이나 먼 별의 표면 온도추정 같은 고온 온도 측정에도 널리 쓰이고 열화상카메라 등을 이용해 체온 정도의 낮은 온도도 측정가능하다.
온도에 따라 쇠가 방출하는 색상. | |
달궈진 쇠. 온도가 높을수록 노란(파장이 짧은) 빛을 방출한다. | 800 K부터 12200 K까지 흑체가 방출하는 색상. 표면온도에 따른 항성의 빛깔과 유사하다. |
5.2. 문제점
위와 같은 장점에도 불구하고 빈 법칙은 [math(\nu)]가 작을 때, 즉 [math(\lambda)]가 큰 장파장일 경우에 실험결과에 잘 맞지 않는 문제점이 있었다.6. 레일리-진스 법칙
1900년 6월, 레일리는 흑체복사를 고전적인 관점에서 연구하여, 흑체 속의 전자가 무질서한 운동을 하고 있지만 결국 진공 속의 전자기파와 마찬가지로 취급한다. 이 법칙을 레일리-진스 법칙(Rayleigh-Jeans law)이라 한다. 즉 흑체 양끝을 마디로 하는 정상파라고 하면 이때 어떤 속이 빈 물체(공동, cavity) 안에서 수없이 여러 번 반사되어 기벽과 열적 평형을 이룬 복사선이 바늘 구멍(hole)으로 나온다고 이상적인 물체를 가정했는데, 이것을 흑체라고 한다.6.1. 상자 속 파동
해석을 용이하게 하기 위해 각 변의 길이가 [math(A)], [math(B)], [math(C)]인 직육면체 흑체가 공간 직교좌표계의 각 축 방향에 나란하게 배치되어있고, 흑체 안에 빛이 갇혀서 내부에서 완벽하게 반사되어 존재하는 이상적인 상황을 가정하자.또한, 고전역학적인 관점에 따라 빛을 파동으로 간주하여 변위 [math({\bf r} = (x,\,y,\,z))]와 시간 [math(t)]에 의존하는 함수로서 다루도록 한다. 파동은 [math(t)]축의 양의 방향으로 진행하며 파수 벡터를 [math({\bf k} = (k_x,\,k_y,\,k_z))], 각진동수를 [math(\omega)], 진동수를 [math(\nu)]라고 하면 다음과 같은 사인파로 표현할 수 있다.
[math(\begin{aligned} A({\bf r},\,t) &= A_0\sin{\left(\frac{\bf k}{\rm rad}\bm\cdot{\bf r} - \frac\omega{\rm rad}t\right)} \\ &= A_0\sin{\left(\frac{k_x}{\rm rad}x + \frac{k_y}{\rm rad}y + \frac{k_z}{\rm rad}z - 2\pi\nu t\right)}\end{aligned})] |
[math(\begin{aligned} 2A &= l\lambda_x \\ 2B &= m\lambda_y \\ 2C &= n\lambda_z\end{aligned})] |
[math(\begin{aligned} \frac{k_x}{\rm rad} &= \frac{\pi l}A \\ \frac{k_y}{\rm rad} &= \frac{\pi m}B \\ \frac{k_z}{\rm rad} &= \frac{\pi n}C \end{aligned})] |
이와 더불어 파동방정식을 만족하므로
[math(\begin{aligned} \frac1{c^2}\frac{\partial^2}{\partial t^2}A({\bf r},\,t) &= \nabla^2A({\bf r},\,t) \\ \Rightarrow \frac{(2\pi\nu)^2}{c^2} &= \frac{{k_x}^2 + {k_y}^2 + {k_z}^2}{\rm rad^2} \\ &= {\left(\frac{\pi l}A\right)}^2 + {\left(\frac{\pi m}B\right)}^2 + {\left(\frac{\pi n}C\right)}^2 \\ \therefore \frac{4\nu^2}{c^2} &= {\left(\frac lA\right)}^2 + {\left(\frac mB\right)}^2 + {\left(\frac nC\right)}^2\end{aligned})] |
[math(\dfrac{l^2}{{\left(\dfrac{2\nu A}c\right)}^2} + \dfrac{m^2}{{\left(\dfrac{2\nu B}c\right)}^2} + \dfrac{n^2}{{\left(\dfrac{2\nu C}c\right)}^2} = 1)] |
아울러, 똑같은 진동수라 하더라도 세부적으론 서로 다른 [math(l,\,m,\,n)]의 상태를 가질 수 있다는 것을 의미하며 이는 곧 진동수마다 상태의 축퇴수(degeneracy)가 각기 다르다는 것을 암시한다.
진동수가 [math(\nu)]에서 [math(\nu+{\rm d}\nu)]로 증가할 때, 각 상태 역시 서로 독립이므로 [math({\rm d}l)], [math({\rm d}m)], [math({\rm d}n)]만큼 증가한다고 볼 수 있고 이 증가분의 총합 [math({\rm d}N)]은 [math({\rm d}N = {\rm d}l{\rm\,d}m{\rm\,d}n)]이다.
구의 경우와 비슷하게 [math(\dfrac{2\nu}c = \kappa)]라 놓고 다음과 같이 구면 좌표계 변환을 하자.
[math(\begin{matrix} \begin{cases}l = \kappa A\sin\theta\cos\phi \\ m = \kappa B\sin\theta\sin\phi \\ n = \kappa C\cos\theta \end{cases} & {\left(\begin{aligned} &\kappa>0 \\ &0\le\theta\le\dfrac\pi2 \\ &0\le\theta\le\dfrac\pi2 \end{aligned}\right)}\end{matrix})] |
[math(\begin{aligned} |J| &= \begin{Vmatrix}A\sin\theta\cos\phi & \kappa A\cos\theta\cos\phi & -\kappa A\sin\theta\sin\phi \\ B\sin\theta\sin\phi & \kappa B\cos\theta\sin\phi & \kappa B\sin\theta\cos\phi \\ C\cos\theta & -\kappa C\sin\theta & 0\end{Vmatrix} \\ &= |\kappa^2ABC\sin\theta| \\ &= \kappa^2ABC\sin\theta\end{aligned})] |
[math(\begin{aligned}{\rm d}N &= {\rm d}l{\rm\,d}m{\rm\,d}n \\ &= \kappa^2ABC\sin\theta{\rm\,d}\kappa{\rm\,d}\theta{\rm\,d}\phi \\ &= {\left(\frac{2\nu}c\right)}^2ABC{\rm\,d}{\left(\frac{2\nu}c\right)}\sin\theta{\rm\,d}\theta{\rm\,d}\phi \\ &= \frac{8\nu^2}{c^3}ABC{\rm\,d}\nu\sin\theta{\rm\,d}\theta{\rm\,d}\phi\end{aligned})] |
[math({\rm d}N(\nu) = \dfrac{4\pi\nu^2}{c^3}V{\rm\,d}\nu)] |
[math(g(\nu){\rm\,d}\nu = \dfrac{8\pi\nu^2}{c^3}{\rm\,d}\nu)] |
6.2. 자외선 파탄
문제는 이 다음 단계에서 레일리와 진스가 했던 가정으로, 앞서 흑체 내부는 열적 평형에 도달한 상태라고 했는데 이러한 계는 자유도 하나당 [math(\dfrac12k_{\rm B}T)]의 평균 에너지를 갖는 것이 알려져 있다. [math(g(\nu){\rm\,d}\nu)]가 갖는 의미는 진동수가 [math(\nu)]인 빛의 상태 밀도이므로 평균 에너지를 상태 밀도에 곱하면 흑체의 에너지 밀도가 얻어질 것이라고 추정한 것이다. 이에 따라 빛은 원편광 성질에 의해 자유도가 2이므로 평균 에너지 [math(\langle E\rangle)]를 [math(\langle E\rangle = k_{\rm B}T)]로 계산하여 에너지 밀도 [math(u_\nu(\nu,\,T){\rm\,d}\nu)]를[math(u_\nu(\nu,\,T){\rm\,d}\nu = \dfrac{8\pi\nu^2k_{\rm B}T}{c^3}{\rm\,d}\nu)] |
이 설명은 빈 법칙과는 정반대로 낮은 진동수(장파장)에서는 맞지만 높은 진동수(단파장)에서는 전혀 맞지 않는다. 에너지 밀도가 진동수의 제곱에 비례하게 되므로 높은 진동수인 빛의 에너지가 무한대로 커지는 것이 문제이다. 즉, 아무리 온도가 낮은 물체라도 높은 진동수에서는 무한대의 에너지를 뿜어낸다는 것이다. 쉽게 말하면 X선, 감마선 등이 무한대로 뿜어져 나온다는 것. 이 터무니 없는 결과를 자외선 파탄이라고 부른다.[10] 이 식이 잘못되었다는 것은 해당 분광 에너지 밀도 [math(u_\nu(\nu,\,T))]를 [math(\nu)]의 모든 범위에서 적분하여 총 에너지 밀도만 구해봐도 알 수 있다. 총 에너지 밀도가 수렴하지 않고 발산한다는 것은 곧, 흑체 내부의 에너지가 무한대라는 것이고 이는 고립계의 에너지가 보존된다는 법칙과 모순된다.
[math(\displaystyle E=\int_0^\infty u_\nu(\nu,\,T){\rm\,d}\nu \rightarrow \infty)] |
7. 플랑크 법칙
플랑크는 각각 1900년 10월과 12월에 출판한 두 논문 《빈의 스펙트럼 방정식의 개선에 대하여》[11]와 《정상 스펙트럼에서 에너지 분포 법칙의 이론에 대하여》[12]에서 기존 빈 법칙이 통계역학적인 관점에서 엔트로피 [math(S)]와 에너지 [math(U)]가 다음 특성을 따른다는 것을 발견하였는데[math(\dfrac{{\rm d}^2S}{{\rm d}U^2}=\dfrac{\sf const.}U)] |
[math(\dfrac{{\rm d}^2 S}{{\rm d}U^2} = \dfrac\alpha{U(\beta +U)})] |
[math(E = \dfrac{C\lambda^{-5}}{e^{\frac{c_2}{\lambda T}} - 1})] |
7.1. 양자가설
플랑크는 공식의 물리적인 의미를 설명하기 위하여 진동하는 입자인 '진동자' 모형을 도입하고, 흑체 내벽은 진동자들로 이루어져 있으며 진동자가 진동하면서 복사를 흡수하고 방출한다고 하였다. 이 가정은 당시엔 ad hoc이었으나 공교롭게도 실제 흑체 스펙트럼의 실험 결과에 매우 잘 부합했으며, 근사를 통해 빈 법칙, 레일리-진스 법칙 등을 모두 설명할 수 있다는 점에서 획기적인 발상으로 평가받고 있다.플랑크는 이 진동자가 주고 받을 수 있는 복사에너지가 진동수에 비례하는 에너지 덩어리로 이루어져 있을 수 밖에 없음을 발견하게 되는데, 이것이 바로 그 유명한 양자가설이며, 이 에너지 덩어리를 양자(quantum)라고 부르게 된다.
1900년 12월 14일 베를린, 유수의 물리학자들이 모인 앞에서 플랑크는 자신의 공식과 양자가설을 발표하여 키르히호프가 남긴 흑체문제에 종지부를 찍었지만, '에너지의 양자화'라는 새로운 개념은 당시 물리학자들에게도, 플랑크 본인에게도 받아들이기 힘든 것이었다. '양자'의 도입은 전자기파의 에너지에 불연속성이 존재함을 인정하는 것이고, '전자기파는 진동수에 관계없이 동일한 온도에서는 동일한 에너지를 갖는다'는 고전 열역학의 법칙(에너지 등분배 법칙)과 배치되는 것이기 때문이었다.
어찌됐든 양자의 도입으로 고전 물리학으로는 해결할 수 없는 난제였던 흑체문제를 해결하고, 새로운 이론과 모델의 필요성을 촉발시켰다는 점에서 오늘날 많은 사람들이 플랑크가 양자가설을 발표한 이 시기를 양자역학이 탄생한 시점으로 보고 있다.
그리고 플랑크가 양자가설을 발표하고 나서 5년 후인 1905년, 양자역학의 서막을 완전히 열어제낀 논문[14]이 발표되는데 바로 알베르트 아인슈타인의 등장이었다.
다시 플랑크 법칙으로 돌아와서, 볼츠만 분포에 따르면 분석 대상이 에너지가 [math(E)]인 상태에 있을 확률은 [math(e^{-\frac E{k_{\rm B}T}})]에 비례한다.
여기서 에너지가 양자화되어있다고 가정, 즉 [math(E = E_n = nh\nu)]라 놓으면 확률밀도함수 [math(P(E))]는 다음과 같이 쓸 수 있다.
[math(\begin{aligned} P(E) &= P(nh\nu) \\ &= Ae^{-\frac{nh\nu}{k_{\rm B}T}}\end{aligned})] |
[math(\displaystyle \sum_{n=0}^\infty P(nh\nu) = \sum_{n=0}^\infty Ae^{-\frac{nh\nu}{k_{\rm B}T}} = 1)] |
[math(\begin{aligned} \sum_{n=0}^\infty P(nh\nu) &= \frac A{1-e^{-\frac{h\nu}{k_{\rm B}T}}} = 1 \\ \therefore A &= 1-e^{-\frac{h\nu}{k_{\rm B}T}} \\ P(nh\nu) &= {\left(1-e^{-\frac{h\nu}{k_{\rm B}T}}\right)}e^{-\frac{nh\nu}{k_{\rm B}T}}\end{aligned})] |
[math(\displaystyle \left\langle E\right\rangle = \sum_{n=0}^\infty E_nP(E_n)= \sum_{n=0}^\infty nh\nu{\left(1-e^{-\frac{h\nu}{k_{\rm B}T}}\right)}e^{-\frac{nh\nu}{k_{\rm B}T}})] |
[math(\begin{aligned} \left\langle E\right\rangle &= \sum_{n=0}^\infty nh\nu(1-\epsilon)\epsilon^n \\ &= h\nu(1-\epsilon)\epsilon\sum_{n=0}^\infty n\epsilon^{n-1} \\ &= h\nu(1-\epsilon)\epsilon\sum_{n=0}^\infty\frac{\rm d}{{\rm d}\epsilon}(\epsilon^n) \\ &= h\nu(1-\epsilon)\epsilon\frac{\rm d}{{\rm d}\epsilon}\sum_{n=0}^\infty\epsilon^n \\ &= h\nu(1-\epsilon)\epsilon\frac{\rm d}{{\rm d}\epsilon}{\left(\frac1{1-\epsilon}\right)} \\ &= h\nu\cancel{(1-\epsilon)}\epsilon\frac1{(1-\epsilon)^{\cancel2}} \\ &= \frac{h\nu\epsilon}{1-\epsilon} \\ &= \frac{h\nu}{\epsilon^{-1} - 1} \\ &= \frac{h\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1}\end{aligned})] |
[math(u_\nu(\nu,\,T){\rm\,d}\nu = \dfrac{8\pi h\nu^3}{c^3}\dfrac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}{\rm\,d}\nu)] |
[math(u_\lambda(\lambda,\,T){\rm\,d}\lambda = \dfrac{8\pi hc}{\lambda^5}\dfrac1{e^{\frac{hc}{\lambda k_{\rm B}T}}-1}{\rm\,d}\lambda)] |
플랑크 법칙은 분광 에너지 밀도 대신 분광 복사휘도(spectral radiance) [math(B_\nu(\nu,\,T))]로 다루는 경우가 많은데, 복사휘도는 점광원이 단위 시간 당 단위 면적에 가하는 에너지로 정의된다. 점광원에서 나온 광선속 다발은 공간에 입체각 [math(4\pi{\rm\,sr})]으로 퍼지므로 분광 에너지 밀도의 증분 [math(u_\nu(\nu,\,T){\rm\,d}\nu)]을 [math(4\pi{\rm\,sr})]으로 나누면 점광원의 분광 에너지 밀도 증분이 된다. 여기에 부피의 증분 [math({\rm d}V)]를 곱하면 점광원의 분광 에너지의 증분 [math({\rm d}E = \dfrac{u_\nu(\nu,\,T){\rm\,d}\nu}{4\pi{\rm\,sr}}{\rm\,d}V)]이 되며 [math({\rm d}V)]는 다시 면적의 증분 [math({\rm d}A)]과 변위의 증분 [math({\rm d}r)]을 이용해서 [math({\rm d}V = {\rm d}A{\rm\,d}r)]로, [math({\rm d}r = c{\rm\,d}t)]로 나타낼 수 있으므로
[math(\begin{aligned} {\rm d}E &= \frac{8\pi h\nu^3}{c^3}\frac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}\frac1{4\pi{\rm\,sr}}{\rm\,d}\nu{\rm\,d}A{\cdot}c{\rm\,d}t \\ &= \dfrac{2h\nu^3}{c^2}\dfrac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}{\rm\,sr^{-1}}{\rm\,d}\nu{\rm\,d}A{\rm\,d}t\end{aligned})] |
[math(B_\lambda(\lambda,\,T) = \dfrac{c_{1L}}{\lambda^5}\dfrac1{e^{\frac{c_2}{\lambda T}}-1})] |
위 식들을 이용해서 흑체복사의 평균 에너지 [math(\overline E)]를 구해보자. 분광 에너지 밀도 [math(u_\nu(\nu,\,T))]는 진동수가 [math(\nu)]인 빛의 에너지 밀도를 의미하므로 [math(u_\nu(\nu,\,T))]를 빛의 에너지 [math(h\nu)]로 나눈 물리량은 (상태수가 포함된) 광자의 밀도로 볼 수 있다. 따라서 [math(\dfrac{u_\nu(\nu,\,T)}{h\nu})]를 모든 진동수 범위에서 적분한 값으로 에너지 밀도 [math(u)]를 나누면 평균 에너지가 될 것이다.
[math(\begin{aligned} \overline E &= \frac{\displaystyle\int_0^\infty u_\nu(\nu,\,T){\rm\,d}\nu}{\displaystyle\int_0^\infty \frac{u_\nu(\nu,\,T)}{h\nu}{\rm\,d}\nu} \\ &= \frac{\displaystyle\int_{0{\rm\,Hz}}^\infty\frac{\cancel{8\pi}h\nu^3}{\cancel{c^3}}\frac{{\rm d}\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1}}{\displaystyle\int_{0{\rm\,Hz}}^\infty\frac{\cancel{8\pi}\nu^2}{\cancel{c^3}}\frac{{\rm\,d}\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1}} \\ &= \frac{\displaystyle h\int_{0{\rm\,Hz}}^\infty\frac{\nu^3{\rm\,d}\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1}}{\displaystyle\int_{0{\rm\,Hz}}^\infty\frac{\nu^2{\rm\,d}\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1}}\end{aligned})] |
[math(\begin{aligned} \overline E &= \frac{\displaystyle \cancel h\int_0^\infty{\left(\frac{k_{\rm B}T}{\cancel h}\right)}^{\cancel4}\frac{x^3{\rm\,d}x}{e^x-1}}{\displaystyle\int_0^\infty{\cancel{\left(\frac{k_{\rm B}T}h\right)}^3}\frac{x^2{\rm\,d}x}{e^x-1}} \\ &= \frac{\displaystyle k_{\rm B}T\int_0^\infty\frac{x^3e^{-x}{\rm\,d}x}{1-e^{-x}}}{\displaystyle\int_0^\infty\frac{x^2e^{-x}{\rm\,d}x}{1-e^{-x}}} \\ &= \frac{\displaystyle k_{\rm B}T\int_0^\infty x^3\sum_{n=1}^\infty{\left(e^{-x}\right)}^n{\rm\,d}x}{\displaystyle\int_0^\infty x^2\sum_{n=1}^\infty{\left(e^{-x}\right)}^n{\rm\,d}x} \\ &= \frac{\displaystyle k_{\rm B}T\sum_{n=1}^\infty\int_0^\infty x^3e^{-nx}{\rm\,d}x}{\displaystyle\sum_{n=1}^\infty\int_0^\infty x^2e^{-nx}{\rm\,d}x}\end{aligned})] |
[math(\begin{aligned} \overline E &= \frac{\displaystyle k_{\rm B}T\sum_{n=1}^\infty\int_0^\infty\frac1{n^4}t^3e^{-t}{\rm\,d}t}{\displaystyle\sum_{n=1}^\infty\int_0^\infty\frac1{n^3}t^2e^{-t}{\rm\,d}t} \\ &= \frac{\displaystyle k_{\rm B}T\sum_{n=1}^\infty\frac1{n^4}\int_0^\infty t^3e^{-t}{\rm\,d}t}{\displaystyle\sum_{n=1}^\infty\frac1{n^3}\int_0^\infty t^2e^{-t}{\rm\,d}t} \\ &= \frac{\zeta(4)\Gamma(4)}{\zeta(3)\Gamma(3)}k_{\rm B}T \\ &= \frac{\dfrac{\pi^4}{90}{\cdot}6}{\zeta(3){\cdot}2}k_{\rm B}T \\ &= \frac{\pi^4}{30\zeta(3)}k_{\rm B}T \\ &= 2.701\,178\,032\cdots k_{\rm B}T \\ &= 3.729\,379\times10^{-23}\,({\rm J/K})\,T\end{aligned})] |
나중에 밝혀진 바에 따르면 결과는 올바르지만, 공식 유도에 결함이 내포되어 있다. 공동의 진동자 에너지는 사실 [math(E_n = nh\nu)]가 아닌 영점 에너지 [math(\dfrac12h\nu)]가 포함된 [math({E_n}' = {\left(n+\dfrac12\right)}h\nu = E_n + \dfrac12h\nu)]로 주어지므로, 위 과정에 대입해서 똑같은 방식으로 유도하면 [math(\langle E'\rangle)]은 [math(\dfrac12h\nu)]가 더해진 값이 나온다. 우선 [math(P({E_n}'))]에 대해
[math(\begin{aligned} P({E_n}') &= P{\left(E_n + \dfrac12h\nu\right)} \\ &= e^{-\frac{E_n}{k_{\rm B}T}}e^{-\frac{h\nu}{2k_{\rm B}T}} = Ae^{-\frac{E_n}{k_{\rm B}T}}\end{aligned})] |
다음으로 [math(\langle E'\rangle)]을 계산해보면
[math(\begin{aligned}\sum_{n=0}^\infty {E_n}'P({E_n}') &= \sum_{n=0}^\infty{\left(E_nP({E_n}) + \frac12h\nu P({E_n})\right)} \\ &= \sum_{n=0}^\infty E_nP({E_n}) + \frac12h\nu\sum_{n=0}^\infty P(E_n) \\ &= \langle E\rangle + \dfrac12h\nu\end{aligned})] |
결국 플랑크의 가정에 따라 유도할 경우 진동자의 영점에너지 [math(\dfrac12h\nu)]만큼을 빼야 맞는다는 것인데, 이에 대해서는 광전효과에서 금속의 일함수만큼이 빠진 에너지 [math(h\nu-W)]를 갖는 전자가 방출되는 현상이나, 형광의 원리에서 들뜬 상태와 바닥 상태의 차이만큼의 에너지에 해당하는 빛이 나오는 것과 같은 맥락으로 해석한다. 요컨대, 공동의 바닥 상태와의 차이만큼이 흑체복사로 나타나는 것이기에 영점 에너지를 고려하지 않아도 된다는 해석이다.
7.2. 빈 법칙과의 관계
앞서 빈 법칙은 단파장(고진동수)에서 설명이 잘 되고 장파장(저진동수)에서는 실제와 맞지 않는다고 하였다. 플랑크 법칙의 분광 복사휘도를 다음과 같이 변형하고[math(\begin{aligned} B_\lambda(\lambda,\,T) &= \frac{2hc^2}{\lambda^5}\frac1{e^{\frac{hc}{\lambda k_{\rm B}T}}-1}{\rm\,sr^{-1}} \\ &= \frac{2hc^2}{\lambda^5}\frac{e^{-\frac{hc}{\lambda k_{\rm B}T}}}{1-e^{-\frac{hc}{\lambda k_{\rm B}T}}}{\rm\,sr^{-1}} \\ &= \frac{2hc^2}{\lambda^5}\sum_{n=1}^\infty{\left(e^{-\frac{hc}{\lambda k_{\rm B}T}}\right)}^n{\rm\,sr^{-1}}\end{aligned})] |
[math(B_\lambda(\lambda,\,T) \approx \dfrac{2hc^2}{\lambda^5}e^{-\frac{hc}{\lambda k_{\rm B}T}}{\rm\,sr^{-1}})] |
7.3. 빈 변위 법칙과의 관계
한편 빈 변위 법칙에 대해서, 똑같이 [math(B_\lambda(\lambda,\,T))]를 [math(\lambda)]에 대해 편미분하여 그 값이 [math(0{\rm\,W/(m^4\,sr)})]이 되는 조건을 찾는다.[math(\begin{aligned} \frac{\partial}{\partial\lambda}B_\lambda(\lambda,\,T) &= \frac{\partial}{\partial\lambda}{\left(\frac{2hc^2}{\lambda^5}\frac1{e^{\frac{hc}{\lambda k_{\rm B}T}}-1}{\rm\,sr^{-1}}\right)} \\ &= \frac{2hc^2}{\lambda^6}\frac1{e^{\frac{hc}{\lambda k_{\rm B}T}}-1}{\left\{\frac{hc}{\lambda k_{\rm B}T}\frac{e^{\frac{hc}{\lambda k_{\rm B}T}}}{e^{\frac{hc}{\lambda k_{\rm B}T}}-1} - 5\right\}}{\rm\,sr^{-1}} = 0{\rm\,W/(m^4\,sr)} \end{aligned})] |
[math(\dfrac{\chi e^\chi}{e^\chi - 1} = 5 \Leftrightarrow \chi e^\chi = 5e^\chi - 5)] |
[math(\begin{aligned}(x+5)e^{x+5} &= xe^{x+5} + \cancel{5e^{x+5}} = \cancel{5e^{x+5}} - 5 \\ &\Rightarrow xe^x = -5e^{-5} \\ \therefore x &= W_0(-5e^{-5}) \\ \chi &= 5 + W_0(-5e^{-5})\end{aligned})] |
[math(\begin{aligned} \lambda_{\sf max} &= \frac{6.626\,070\,15\times10^{-34}{\rm\,\cancel{J{\cdot}s}}\cdot299\,792\,458{\rm\,m/\cancel s}}{\chi\cdot1.380\,649\times10^{-23}{\rm\,\cancel J/K}\cdot T} \\ &\approx \frac{2.897\,772{\rm\,mm{\cdot}K}}T\end{aligned})] |
참고로 복사휘도가 극대인 진동수를 계산할 때에는 일반적인 [math(c = \lambda\nu)] 공식을 쓰면 안 된다. [math(B_\nu(\nu,\,T) = \dfrac c{\nu^2}B_\lambda(\lambda(\nu),\,T) = \dfrac{2h\nu^3}{c^2}\dfrac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}{\rm\,sr^{-1}})]로 두 식의 형태가 다르기 때문에 미분해서 [math(\nu_{\sf max})]를 계산해보면
[math(\nu_{\sf max} = \dfrac{\{3+W_0(-3e^{-3})\}k_{\rm B}T}h = 5.878\,926\times10^{10}\,({\rm Hz/K}){\cdot}T)] |
7.4. 레일리-진스 법칙과의 관계
[math(\displaystyle e^x = 1 + \dfrac x{1!} + \dfrac{x^2}{2!} + \cdots)]에서 [math(x\ll1)]이면 [math(e^x \approx 1 + x)]로 근사할 수 있다는 점을 적용하면[math(\begin{aligned} u_\nu(\nu,\,T){\rm\,d}\nu &= \frac{8\pi h\nu^3}{c^3}\frac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}{\rm\,d}\nu \\ &\approx \frac{8\pi h\nu^3}{c^3}\frac1{\dfrac{h\nu}{k_{\rm B}T}}{\rm\,d}\nu \\ &= \frac{8\pi\nu^2k_{\rm B}T}{c^3}{\rm\,d}\nu\end{aligned})] |
7.5. 슈테판-볼츠만 상수 구하기
오늘날의 슈테판-볼츠만 법칙은 엄밀하게는 흑체의 복사발산도 [math(M_{\rm e})]가 절대온도의 네제곱에 비례한다는 법칙으로 정의된다.[15][math(M_{\rm e} = \sigma T^4)] |
이때, 주의할 점이 있는데 흑체에서 나오는 빛의 복사휘도는 람베르트의 코사인 법칙(Lambert's cosine law)[16]을 따르기 때문에 피적분 함수가 [math(B_\nu(\nu,\,T)\cos(\theta/{\rm rad}){\rm\,d}\nu)][17]가 되며, 복사휘도는 관측기에 검출되는 물리량이기 때문에 관측자 방향으로 방사되지 않는 빛은 계산에서 제외되어 [math(0\le(\theta/{\rm\,rad})\le\dfrac\pi2)] 범위로 적분 범위가 제한된다. 모든 공간으로 적분을 했던 에너지 밀도와 비교하면 반구에 대한 적분으로 볼 수 있으며 정리하면 복사발산도에 관한 식은 다음과 같이 주어진다.
[math(\begin{aligned}M_{\rm e} &= \iint_\Omega\int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T)\cos\frac\theta{\rm rad}{\rm\,d}\nu{\rm\,d}\Omega \\ &= \iint_\Omega\int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,sr}\cos\frac\theta{\rm rad}{\rm\,d}\nu{\rm\,d}\frac\Omega{\rm\,sr} \\ &= \int_{0{\rm\,rad}}^{2\pi{\rm\,rad}}\int_{0{\rm\,rad}}^{\frac\pi2{\rm\,rad}}\int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,sr}\cos\frac\theta{\rm rad}{\rm\,d}\nu\sin\frac\theta{\rm rad}{\rm\,d}\frac\theta{\rm rad}{\rm\,d}\frac\phi{\rm rad} \\ &= \int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,sr}{\rm\,d}\nu \int_{0{\rm\,rad}}^{\frac\pi2{\rm\,rad}}\cos\frac\theta{\rm rad}\sin\frac\theta{\rm rad}{\rm\,d}\frac\theta{\rm rad}\int_{0{\rm\,rad}}^{2\pi{\rm\,rad}}{\rm\,d}\frac\phi{\rm rad} \\ &= \int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,sr}{\rm\,d}\nu \frac1{\cancel2}{\cdot}\cancel2\pi \\ &= \pi\int_{0{\rm\,Hz}}^\infty{\left(\frac{2h\nu^3}{c^2}\frac1{e^{\frac{h\nu}{k_{\rm B}T}}-1}{\rm\,sr^{-1}}\right)}{\rm\,sr}{\rm\,d}\nu \\ &= \pi\int_{0{\rm\,Hz}}^\infty\frac{2h\nu^3}{c^2}\frac{{\rm d}\nu}{e^{\frac{h\nu}{k_{\rm B}T}}-1} \\ &= \pi\int_0^\infty \frac{2h}{c^2}{\left(\frac{k_{\rm B}T}h\right)}^4\frac{x^3{\rm\,d}x}{e^x-1} \quad{\left(\frac{h\nu}{k_{\rm B}T} = x \Rightarrow \nu = \frac{k_{\rm B}T}hx,\,{\rm d}\nu = \frac{k_{\rm B}T}h{\rm\,d}x\right)} \\ &= \frac{2\pi{k_{\rm B}}^4T^4}{c^2h^3}\int_0^\infty\frac{x^3e^{-x}{\rm\,d}x}{1-e^{-x}} \\ &= \frac{2\pi{k_{\rm B}}^4T^4}{c^2h^3}\int_0^\infty x^3\sum_{n=1}^\infty e^{-nx}{\rm\,d}x \\ &= \frac{2\pi{k_{\rm B}}^4T^4}{c^2h^3}\sum_{n=1}^\infty\int_0^\infty x^3e^{-nx}{\rm\,d}x \\ &= \frac{2\pi{k_{\rm B}}^4T^4}{c^2h^3}\zeta(4)\Gamma(4) \\ &= \frac{2\pi{k_{\rm B}}^4T^4}{c^2h^3}\frac{\pi^4}{90}{\cdot}6 \\ &= \frac{2\pi^5{k_{\rm B}}^4}{15c^2h^3}T^4\end{aligned})] |
한편, 앞서 [math(B_\nu(\nu,\,T) = \dfrac c{4\pi{\rm\,sr}}u_\nu(\nu,\,T) \Leftrightarrow u_\nu(\nu,\,T) = \dfrac{4\pi}cB_\nu(\nu,\,T){\rm\,sr})]으로 에너지 밀도가 복사휘도와 상수배 관계에 있으므로 에너지 밀도 역시 절대온도의 네제곱에 비례한다. 주의해야할 점은 에너지 밀도는 입체각 [math(4\pi{\rm\,sr})]에 작용하는 복사휘도로 환산된 값과도 같은데 복사휘도처럼 [math(\cos(\theta/{\rm rad}))]를 곱하거나 적분 범위를 제한하는 등의 조건이 들어가지 않기 때문에 비례관계식의 상수까지 [math(\dfrac{4\pi}c)]배가 되지 않는다는 점이다. 식을 전개하면 다음과 같이 된다.
[math(\begin{aligned} u &= \int_{0{\rm\,Hz}}^\infty u_\nu(\nu,\,T){\rm\,d}\nu \\ &= \frac{4\pi}c\int_{0{\rm\,Hz}}^\infty B_\nu(\nu,\,T){\rm\,sr}{\rm\,d}\nu \\ &= \frac{4\cancel\pi}c\frac{M_{\rm e}}{\cancel\pi} \\ &= \frac{4\sigma}cT^4 \\ \end{aligned})] |
8. 기타
색온도, 화이트밸런스는 이 문서로 연결된다. 디스플레이나 조명, 사진 등의 분야로 인해 해당 검색어를 입력했다가 생뚱맞은 물리학 이론으로 리다이렉트가 되어서 당황할 수 있으나, 사실 연관성이 매우 큰 내용이다. 색온도는 빈 변위 법칙에서 유도되며, 화이트밸런스는 여기서 나온 수치 모델로 색감을 보정하는 작업이기 때문이다.
[1]
현대의 강재는 등급마다 열처리 단계에서 필요한 온도와 시간을 설정하고, 이를 지키지 못하면 해당 등급을 충족하지 못하는 것으로 간주한다.
[2]
아래 표를 보면 알 수 있지만 복사발산도와 복사조도는 그 미분인 분광 복사발산도와 분광 복사조도를 포함해서 차원도 단위도 완전히 동일한데, 둘을 구분하는 이유는 흑체가 본질적으로 복사율 [math(\varepsilon)]이 [math(\varepsilon=1)]인 이상적인 존재인 반면, 실제로 존재하는 복사원은 [math(\varepsilon<1)]인 데서 유래한다. 따라서 흑체는 복사발산도와 복사조도가 같은 물질이라고 할 수 있다.
[3]
즉, 공간으로 복사되기 전 한 점에 모여있다고 가정할 때의
[4]
반사되는 복사속을 [math(J_{\rm e,\,r})], 투과되는 복사속을 [math(J_{\rm e,\,tr})]로 나타낼 때, [math(J_{\rm e} = M_{\rm e} + J_{\rm e,\,r} + J_{\rm e,\,tr})]로 정의된다. 즉 상술한 복사발산도는 복사도의 일부에 해당하는 셈. 흑체가 아닌 실제 복사원을 설명하는 물리량이라고 볼 수 있다. 불투명한 물체일 경우 투과되는 복사속이 없어 [math(J_{\rm e,\,tr} = 0{\rm\,W/m^2})]가 된다.
[5]
가시광선 복사와
적외선 복사를 가리키는 말이다.
[6]
풀네임은 빌헬름 카를 베르너 오토 프리츠 프란츠 빈(Wilhelm Carl Werner Otto Fritz Franz Wien)
[7]
여담이지만 빈이 발표한 논문
《흑체복사 스펙트럼의 에너지 분포에 대하여》에서는 복사강도(Intensiät der Strahlung)라고 표현했는데, 내용상으로는 분광 복사휘도이며 오늘날 기준으로 복사강도는 점광원이 단위 시간당 가하는 에너지(
일률), 즉 단위가 [math(\rm W/sr)]으로 다른 물리량이다. 분광 복사휘도는 점광원에서 나오는 전자기파 스펙트럼 중 특정 파장 혹은 진동수의 빛이 단위 시간당 단위 면적에 가하는 에너지(단위는 파장의 경우 [math(\rm W/(nm{\cdot}sr{\cdot}m^2))], 진동수의 경우 [math(\rm W/(Hz{\cdot}sr{\cdot}m^2))])이다.
[8]
후술하겠지만 분광 복사휘도 [math(B_\nu(\nu,\,T))]는 분광 에너지 밀도(혹은 에너지 밀도 분포 함수) [math(u_\nu(\nu,\,T))]와 [math(B_\nu(\nu,\,T) = \dfrac c{4\pi{\rm\,sr}}u_\nu(\nu,\,T))]로 서로 상수배 관계에 있다. 따라서 복사휘도가 가장 강하다는 것은 곧 에너지 밀도가 가장 높다는 말과도 같다.
[9]
위키피디아의 원편광 문서 참고. 빛의 진행 방향에 대하여 왼손방향과 오른손방향의 상태가 가능하다.
[10]
여담으로 '자외선 파탄'의 원문은 ultraviolet catastrophe으로서, 우리 우주가 자외선의 빛으로 뒤덮이게 된다는 해석을 단적으로 드러낸 표현, 즉 '자외선 재앙'이 좀 더 적절한 번역이다. 일본에서 '어떤 이론이 특정 경우에 한하여 실제 현상을 잘 설명하지 못하는(실제 현상에 어긋나는) 것'을 '파탄'(破綻)이라고 표현한다는 점을 감안하면, '자외선 파탄'이라는 용어는 일본에서 쓰이는 '자외파탄(紫外破綻)'이 직역된 것으로 보이며, '자외파탄'은 원문이 의미하는 '자외선 재앙'보단 레일리-진스 법칙이 자외선 영역에서 잘 맞지 않는다는 것을 표현한 것에 가깝다.
[11]
Deutschen Physikalischen Gesellschaft(2), 202-204,
영어 버전
[12]
Deutschen Physikalischen Gesellschaft(2), 237-245,
영어 버전
[13]
정확히는 (정수)[math(\times)](
플랑크 상수) 배
[14]
광양자설, 아인슈타인에게 노벨상의 영예를 안겨준 것이 바로 이 논문이다.
[15]
좀 더 정확하게는 복사율 [math(\varepsilon\,(0<\varepsilon<1))]까지 고려해서 [math(M_{\rm e} = \varepsilon\sigma T^4)]으로 쓰는 게 맞긴 하다. 다만 흑체는 [math(\varepsilon=1)]이기 때문에 고려하지 않아도 된다.
[16]
광원을 똑바로 바라봤을 때의 입사광을 기준으로 [math(\theta)]만큼 비스듬한 각도에서 바라봤을 때 그 세기가 [math(\cos(\theta/{\rm rad}))]배가 된다는 법칙.
[17]
[math(\cos(\theta/{\rm rad}))] 표기에 대해서는
삼각함수 문서의 고찰 항목 참고. 아래 수식에서 [math(\dfrac\Omega{\rm sr})]와 같이 나타내어진 표기 역시 분광 복사휘도의 식에 [math(\rm sr^{-1})]이 포함되어있기 때문에 단위 관계를 맞추기 위해서는 이와 같이 엄밀한 표기가 필요하다.