온도와 상태 방정식
온도계를 대면 숫자가 나온다. 그 숫자는 무엇을 재는 것일까? 압력계는 무엇을 재고, 기체 1몰이 22.4 L를 차지한다는 말은 언제 맞고 언제 틀릴까? 이 책은 이런 거시적인 양을 분자들의 운동에서 끌어내는 방법을 배운다. 첫 장에서는 그 출발점인 온도와 상태 방정식을 다룬다. 아래 상자에는 원판 분자 120개가 있다. 분자가 벽에 부딪칠 때마다 벽은 운동량을 받고, 이 충격량을 시간과 벽 길이로 나누면 압력이 된다. 오른쪽 피스톤을 끌어 넓이를 줄이거나, 왼쪽 열벽의 온도 막대를 끌어 분자를 데워 보라. 벽에서 잰 압력은 늘 \(NkT/A\)를 따라간다. 이것이 2차원 이상 기체 법칙이다. 둘의 비율이 1보다 조금 큰 까닭은 분자에 크기가 있기 때문이고, 4절에서 다룬다.
- 열평형과 열역학 제0법칙으로 온도를 정의하고, 단열벽·투열벽, 상태 변수와 상태 방정식의 뜻을 설명한다.
- 정적 기체 온도계의 압력을 \(P_\mathrm{tp}\to0\)으로 외삽해 이상 기체 온도를 얻고, 켈빈·섭씨·화씨 눈금과 2019년 SI 재정의를 설명한다.
- 분자 충돌에서 \(PV = NkT\)를 이끌어 내고, 돌턴 분압 법칙과 몰분율로 혼합 기체의 분압을 구하며, MFC 표준 유량(sccm)을 몰 유량과 챔버 분압으로 바꾼다.
- 압축 인자 \(Z\)와 비리얼 전개로 실제 기체의 편차를 기술하고, 제2 비리얼 계수 \(B(T)\)와 보일 온도로 인력과 반발의 균형을 읽으며, 고압 실린더의 실제 기체량을 계산한다.
- 판데르발스 식의 두 보정항을 해석하고, 등온선의 불안정 고리와 임계점 조건에서 \(T_\mathrm c, P_\mathrm c, V_\mathrm c\)를 \(a, b\)로 나타낸다.
- 환산 변수와 대응 상태 원리로 일반화 압축 인자 차트를 읽고, 열팽창 계수 \(\alpha\)와 등온 압축률 \(\kappa_T\)를 계산해 \((\partial P/\partial T)_V = \alpha/\kappa_T\)를 응용한다.
열평형과 열역학 제0법칙
열역학은 관심 있는 부분인 계(system)와 나머지인 주위(surroundings)를 나누는 데서 시작한다. 계는 챔버 속 기체일 수도, 웨이퍼 한 장일 수도, 버블러 속 액체 전구체일 수도 있다. 계와 주위를 가르는 벽의 성질이 중요하다. 물질도 에너지도 통과시키지 않고 열도 막는 벽을 단열벽(adiabatic wall), 움직이지 않고 물질도 막지만 열은 통과시키는 벽을 투열벽(diathermal wall)이라 한다. 보온병의 진공 이중벽은 단열벽에 가깝고, 얇은 금속판은 투열벽에 가깝다.
두 계를 투열벽으로 맞대 두면 처음에는 압력, 부피, 전기 저항 같은 거시적 성질이 변하다가 결국 더 변하지 않는 상태에 이른다. 이 상태를 열평형(thermal equilibrium)이라 한다. 평형 상태는 몇 개의 거시 변수만으로 완전히 정해진다. 기체 한 덩어리라면 \(P, V, n\)과 아직 정의하지 않은 한 가지 변수로 충분하다. 이런 변수를 상태 변수(state variable)라 부른다. 분자 \(10^{23}\)개의 위치와 속도를 몰라도 평형 상태는 몇 개의 숫자로 기술된다는 것, 이것이 열역학이 가능한 이유다.
분자 수준에서는 무슨 일이 일어날까? 아래 시뮬레이터에서 왼쪽 칸에는 빠른(뜨거운) 원판이, 오른쪽 칸에는 느린(차가운) 원판이 있다. 칸막이를 투열벽으로 바꾸면 칸막이가 양쪽 원판과 부딪치며 에너지를 나른다. 칸막이는 열용량이 작은 얇은 막이다. 부딪친 분자를 막 자신의 온도에 맞는 속도로 되튕겨 내면서 에너지를 받거나 내주고, 막의 온도는 그 에너지 장부에 따라 오르내린다(전체 에너지는 보존된다).
움직이지 않는 투열벽으로 나뉜 두 기체가 있다. 원판 수는 양쪽이 같지만 왼쪽 칸이 더 좁고(넓이 35 %), 처음에는 왼쪽이 더 뜨겁다. 오래 기다리면 양쪽에서 같아지는 것은 무엇일까?
투열벽을 사이에 둔 두 기체의 열평형
평형에 이른 두 계 A와 B를 각각 제3의 계 C(온도계)와 맞대 보자. 실험은 다음 사실을 알려 준다.
계 A와 B가 각각 계 C와 열평형에 있으면, A와 B도 서로 열평형에 있다.
열역학 제0법칙(zeroth law of thermodynamics)은 너무 당연해 보여서 제1·2법칙이 나온 뒤에야 이름이 붙었다. 하지만 이 법칙이 있어야 온도를 정의할 수 있다. "열평형에 있다"는 관계가 전이적이므로, 서로 열평형인 계들을 한 묶음으로 모을 수 있다. 각 묶음에 붙이는 이름표가 온도(temperature)다. 온도계는 C 역할을 하는 작은 계이고, 그 이름표를 숫자로 읽을 수 있게 만든 장치다. 위 시뮬레이터에서 보았듯 분자 수준에서 "같은 온도"란 분자당 평균 운동 에너지가 같다는 뜻이다. 5장에서는 이것을 더 일반적인 통계적 정의 \(1/T = \partial S/\partial U\)로 바꾼다.
온도계가 쓸모 있으려면 재는 대상보다 훨씬 작아야 한다. 온도계도 계이므로 맞대는 순간 대상과 열을 주고받아 대상의 온도를 바꾼다. 열용량이 \(C_s\)인 시료(온도 \(T_s\))에 열용량 \(C_p\)인 탐침(온도 \(T_p\))을 대면, 열이 바깥으로 새지 않는 한 둘은 \(T_f = (C_sT_s + C_pT_p)/(C_s + C_p)\)에서 평형을 이룬다(열용량은 2장에서 정식으로 다룬다).
열용량 20인 작은 시료가 150이다. 여기에 열용량 2.0인 25 °C 열전대 탐침을 대면, 탐침이 읽는 평형 온도는 이고 오차는 다. 오차를 1 °C 아래로 줄이려면 탐침의 열용량이 이하여야 한다.
평형 상태가 몇 개의 변수로 정해진다면, 그 변수들 사이에는 관계식이 있어야 한다. 순수한 기체나 액체 한 덩어리에서는 \(P, V, T, n\) 가운데 셋을 정하면 나머지 하나가 정해진다. 이 관계 \(f(P, V, T, n) = 0\)을 상태 방정식(equation of state)이라 한다. 열역학 법칙은 상태 방정식을 알려 주지 않는다. 상태 방정식은 실험으로 재거나, 분자 모형과 통계역학(6장)으로 계산해야 한다. 이 장의 나머지는 가장 중요한 상태 방정식들인 이상 기체, 비리얼 전개, 판데르발스 식과 그 일반화를 다룬다.
온도 눈금과 기체 온도계
제0법칙은 온도가 있다는 것만 알려 주고 숫자는 정해 주지 않는다. 숫자를 붙이려면 온도에 따라 단조롭게 변하는 성질(온도 속성)을 하나 골라 기준점에 값을 정해야 한다. 수은 기둥의 길이, 백금선의 저항, 열전대의 기전력이 모두 쓰인다. 문제는 이 눈금들이 기준점(예: 물의 어는점 0, 끓는점 100)에서는 일치하지만 그 사이와 바깥에서 서로 조금씩 어긋난다는 것이다. 수은의 팽창과 백금의 저항이 온도에 똑같이 선형일 이유가 없기 때문이다.
기체는 예외적이다. 부피를 일정하게 두고 압력을 재는 정적 기체 온도계(constant-volume gas thermometer)에서는 기체를 묽게 할수록 어떤 기체를 쓰든 같은 눈금이 나온다. 그 공통 눈금이 이상 기체 온도 눈금이고, 3장에서 보듯 카르노 기관으로 정의하는 열역학적 온도와 정확히 같다.
물의 삼중점(triple point)은 얼음·물·수증기가 함께 있는 유일한 온도와 압력(611.657 Pa)이라 재현성이 아주 좋다. 온도계를 삼중점 칸에 넣었을 때의 압력을 \(P_\mathrm{tp}\), 측정할 계에 넣었을 때의 압력을 \(P\)라 하면 이상 기체 온도는
왜 외삽이 필요할까? 기체 분자 사이에는 인력과 반발이 있어서 압력이 온도에 정확히 비례하지 않는다. 그 편차는 밀도에 비례하므로(4절의 비리얼 계수), 밀도를 0으로 보내면 사라진다. 아래 시뮬레이터에서 직접 외삽해 보자.
- 유리구에 기체를 넣고 0–100 °C에서 압력을 재면 직선이 나온다(샤를·게이뤼삭).
- 기체 양을 줄이면 기울기는 작아지지만 여전히 직선이다.
- 세 직선을 압력 0까지 연장하면 모두 −273.15 °C 한 점에서 만난다. 기체 종류를 바꿔도 같다.
- 그 점을 원점으로 삼은 눈금이 켈빈 눈금이다. 압력은 절대 온도에 비례한다.
기체 온도계 외삽 실험
질소를 채운 정적 기체 온도계로 끓는 물(1 atm)의 온도를 쟀다. 삼중점 압력이 \(P_\mathrm{tp} = 80.00\) kPa일 때 끓는 물 속 압력은 109.343 kPa, 기체를 덜어 \(P_\mathrm{tp} = 40.00\) kPa로 만들었을 때는 54.6559 kPa였다. 이상 기체 온도를 구하라.
풀이 보기
켈빈과 2019년 SI. 2019년부터 켈빈은 볼츠만 상수를 정확히 \(k = 1.380649\times10^{-23}\) J/K로 고정해 정의한다. 온도는 이제 "분자 한 개의 에너지 척도 \(kT\)"를 재는 양이 되었고, 물의 삼중점 273.16 K는 측정값(불확도 약 0.1 mK)으로 바뀌었다. 실용 눈금은 다음과 같다.
실험실과 팹의 온도 측정은 ITS-90의 고정점(물의 삼중점, 갈륨 녹는점 302.9146 K, 인듐 어는점 429.7485 K, 주석 505.078 K, 아연 692.677 K, 알루미늄 933.473 K, 은 1234.93 K 등)으로 교정한 백금 저항 온도계(Pt100)와 열전대를 쓴다. 1000 °C가 넘는 RTP 챔버처럼 접촉할 수 없는 곳에서는 복사로 재는 파이로미터를 쓰며, 그 원리(플랑크 법칙)는 7장에서 다룬다.
이상 기체와 돌턴 분압 법칙
머리 그림에서 본 것을 식으로 옮겨 보자. 넓이 \(A\)인 2차원 상자에서 질량 \(m\), 속도 \(v_x\)인 원판이 오른쪽 벽에 부딪치면 벽은 운동량 \(2mv_x\)를 받는다. 이 원판은 시간 \(2L_x/|v_x|\)마다 같은 벽에 돌아오므로 벽이 받는 평균 힘은 \(mv_x^2/L_x\)다. 모든 원판에 대해 더하고 벽 길이 \(L_y\)로 나누면 (2차원) 압력은
이렇게 얻은 이상 기체(ideal gas) 상태 방정식은 분자의 크기와 분자 사이 힘을 무시한 결과다.
피스톤 상자로 목표 상태 만들기
압력의 SI 단위는 Pa(= N/m²)이지만, 현장에서는 bar(10⁵ Pa), atm(101 325 Pa), Torr(= 1/760 atm = 133.322 Pa), mTorr, psi를 섞어 쓴다. 이상 기체 식에 넣기 전에 반드시 Pa로 바꾸는 습관을 들이자. 반도체 공정 압력은 대기압의 10⁻⁹에서 실린더 압력 150 bar까지 열네 자릿수에 걸친다.
건식 식각은 1–100 mTorr, 스퍼터링은 수 mTorr, LPCVD는 0.1–1 Torr, PECVD는 1–10 Torr, 상압 공정은 760 Torr 근처에서 일한다. 장비 화면은 Torr, 진공 게이지 사양은 Pa나 mbar, 배관 압력계는 psi나 kPa(게이지 압력)로 적는 경우가 많아 단위 혼동이 흔한 사고 원인이다. 압력은 곧 분자 수밀도 \(n_V = P/kT\)이고, 평균 자유 행로와 플라스마 밀도, 이온의 충돌 횟수가 모두 여기서 정해진다(8장). → EtchBook · 플라스마: 식각의 재료
부피 50 L인 식각 챔버의 압력이 5.0 mTorr, 기체 온도가 25 °C다. 챔버 속 분자 수와 수밀도를 구하라.
풀이 보기
혼합 기체. 서로 반응하지 않는 기체를 섞으면 각 성분은 다른 성분이 없는 것처럼 행동한다. 성분 \(i\)가 혼자 전체 부피를 차지할 때 낼 압력을 분압(partial pressure) \(P_i\)라 하면
같은 온도에서 가벼운 원판 30개가 든 상자와, 질량이 16배인 무거운 원판 30개가 든 상자가 있다. 무거운 원판 상자의 압력은 가벼운 원판 상자의 몇 배일까?
돌턴 분압 법칙 실험: A만, B만, A+B
MFC와 표준 유량. 공정 장비는 기체 유량을 질량 유량 제어기(MFC)로 맞춘다. MFC는 열식 센서로 질량(몰) 유량을 재지만 표시는 sccm(standard cubic centimeter per minute)으로 한다. 1 sccm은 "표준 상태로 환산했을 때 1분에 1 cm³"인 몰 유량이다. 반도체 업계의 표준 상태는 보통 0 °C, 1 atm이다(일부 장비·가스 회사는 20 °C나 25 °C를 쓰므로 사양서를 확인해야 한다).
MFC가 100을 흘리면 몰 유량은 , 분자로는 초당 개다. 이 기체가 400, 1.00인 챔버에 들어가면 실제 부피 유량은 로, 표준 유량의 배다.
여러 기체를 MFC로 섞어 챔버에 흘리고 펌프로 일정한 압력을 유지하면, 잘 섞인 정상 상태에서 각 성분의 몰분율은 유량 비와 같다(모든 성분이 같은 비율로 배기된다고 가정). 따라서 성분 \(i\)의 분압은 \(P_i = (\dot n_i/\sum\dot n_j)\,P\)다. 막의 성장 속도나 도핑 농도는 전체 압력이 아니라 반응 기체의 분압에 따라 정해지는 경우가 많다.
SiH₄ 50 sccm과 N₂ 950 sccm을 섞어 전체 압력 2.0 Torr로 유지하는 CVD 챔버가 있다. (a) SiH₄의 몰 유량과 초당 분자 수, (b) SiH₄의 분압(Torr, Pa)을 구하라.
풀이 보기
CVD에서 막의 성장 속도는 표면에 도달하는 반응 기체의 플럭스 \(\propto P_i/\sqrt{2\pi mkT}\)(8장)와 표면 반응 속도로 정해진다. 그래서 레시피는 "SiH₄ 50 sccm / N₂ 950 sccm / 2 Torr"처럼 유량 비와 전체 압력을 함께 적는다. 도펀트 기체는 보통 희석해서 공급하므로(예: H₂ 속 1 % B₂H₆), 실제 도펀트 분압은 희석률 × 유량 비 × 전체 압력으로 계산해야 한다. 액체 전구체(TEOS, TMA)는 버블러 온도로 증기압을 정해 운반 기체에 실어 보내며, 증기압은 9장에서 다룬다. → ProcessBook · 박막 증착
실제 기체: 압축 인자와 비리얼 전개
이상 기체 법칙에서 얼마나 벗어났는지 재는 가장 간단한 양은 압축 인자(compressibility factor)다.
상온에서 대부분의 기체는 수십 bar까지 \(Z\)가 1보다 약간 작다가, 수백 bar에서 1을 넘는다. 분자 사이 인력은 수 Å 거리에서 작용하는 약한 힘이라 저압에서 먼저 드러나고, 반발은 분자가 서로 닿을 만큼 밀도가 높을 때 우세해진다. H₂와 He은 인력이 아주 약해 상온에서 처음부터 \(Z > 1\)이다.
밀도가 낮을 때의 편차는 밀도에 대한 거듭제곱 급수로 정리할 수 있다. 이것이 비리얼 전개(virial expansion)다.
\(B\)의 부호는 인력과 반발의 줄다리기를 그대로 보여 준다. 반지름 \(\sigma/2\)인 단단한 공(강체 구)이면 \(u\)가 \(r<\sigma\)에서 무한대이고 바깥에서 0이므로 \(B = \tfrac{2\pi}{3}N_\mathrm A\sigma^3\), 분자 부피의 4배인 양수다. 인력이 있으면 \(e^{-u/kT} - 1 > 0\)인 영역이 생겨 \(B\)를 음수 쪽으로 끌어내리고, 이 효과는 저온에서 클수록 커진다. 따라서 \(B(T)\)는 저온에서 음수, 고온에서 양수이고, 그 사이에 \(B = 0\)인 보일 온도(Boyle temperature) \(T_\mathrm B\)가 있다. 보일 온도에서는 저압 영역의 \(Z\)가 1에 가장 가깝게 오래 머물러 기체가 넓은 압력 범위에서 이상 기체처럼 보인다.
| 기체 | B(273 K) [cm³/mol] | B(600 K) [cm³/mol] | 보일 온도 TB [K] |
|---|---|---|---|
| He | +12.0 | +10.4 | 22.6 |
| N₂ | −10.5 | +21.7 | 327.2 |
| Ar | −21.7 | +11.9 | 411.5 |
| CO₂ | −149.7 | −12.4 | 714.8 |
위 표는 실측 대표값이다. 상온의 N₂는 보일 온도 바로 아래라 \(B\)가 작아서 150 bar에서도 \(Z\approx1\)인 반면, CO₂는 \(B\)가 크게 음수라 수십 bar만 되어도 이상 기체와 크게 다르다. 분자 크기 효과가 압력을 높인다는 것을 2차원 원판 기체로 직접 확인해 보자.
상온(300 K)의 CO₂를 1 bar에서 50 bar로 압축했다. 압축 인자 \(Z\)는 어떻게 될까?
압축 인자 Z와 비리얼 계수 B(T)
고압 실린더 속 기체의 양. 47 L 실린더에 150 bar로 충전된 기체가 몇 mol인지, 압력이 30 bar로 떨어졌을 때 얼마가 남았는지는 매일 쓰는 계산이다. \(n = PV/ZRT\)이므로 \(Z\)를 모르면 많게는 수십 %씩 틀린다. 아래 \(Z\)는 실측에 잘 맞는 펭–로빈슨 상태 방정식(판데르발스 식을 개량한 3차 상태 방정식)으로 계산한 대표값이다.
25 °C, 절대 압력 150인 47 L 실린더 속 Ar은 \(Z\) = , 양은 로 이상 기체 계산보다 많다. 같은 조건의 N₂는 \(Z\) = 이고, SiH₄라면 \(Z\) = 로 이상 기체 값의 배가 들어간다.
47 L 실린더에 25 °C, 150 bar(절대 압력)의 Ar이 들어 있다. 이 조건의 압축 인자는 \(Z = 0.925\), 30 bar에서는 \(Z = 0.975\)다. (a) 실린더 속 Ar의 몰수와 질량, (b) 압력이 30 bar로 떨어졌을 때 처음 양의 몇 %가 남았는가?
풀이 보기
특수가스 실린더는 가스 캐비닛에서 압력 변환기로 잔량을 감시하고, 정해진 압력 아래로 떨어지면 교체한다. 압력이 몰수에 정비례한다는 이상 기체 가정은 N₂·He·H₂ 같은 영구 기체에는 대체로 맞지만, CO₂·N₂O처럼 상온에서 액화하는 기체나 임계 온도가 상온 근처인 SiH₄(\(T_\mathrm c\approx270\) K)에서는 크게 틀린다. 그래서 이런 기체는 충전량을 압력이 아니라 무게로 관리한다. 압력이 떨어지는 마지막 구간에서는 불순물(수분·잔류 기체)의 비율이 커지므로 실린더를 완전히 비우지 않고 교체 압력을 둔다. → InfraBook · 가스 공급
판데르발스 기체와 임계점
1873년 판데르발스는 분자 그림에서 바로 나오는 두 보정으로 이상 기체 식을 고쳤다. 첫째, 분자에는 크기가 있어 한 분자가 움직일 수 있는 부피는 \(V\)가 아니라 \(V - nb\)다. 둘째, 기체 안쪽 분자는 사방에서 같은 힘으로 끌리지만 벽 근처 분자는 안쪽으로만 끌려 벽에 덜 세게 부딪친다. 이 압력 감소는 벽에 오는 분자 수(\(\propto n/V\))와 끌어당기는 이웃 수(\(\propto n/V\))의 곱에 비례한다.
| 기체 | a [Pa·m⁶/mol²] | b [cm³/mol] | 실측 Tc [K] | 실측 Pc [bar] | 실측 Zc |
|---|---|---|---|---|---|
| He | 0.00346 | 23.8 | 5.2 | 2.27 | 0.302 |
| H₂ | 0.0245 | 26.5 | 33.1 | 13.0 | 0.303 |
| N₂ | 0.137 | 38.7 | 126.2 | 34.0 | 0.289 |
| Ar | 0.1355 | 32.0 | 150.7 | 48.6 | 0.291 |
| CH₄ | 0.2303 | 43.1 | 190.6 | 46.0 | 0.286 |
| CO₂ | 0.364 | 42.7 | 304.1 | 73.8 | 0.274 |
| H₂O | 0.5536 | 30.5 | 647.1 | 220.6 | 0.229 |
두 항이 압력을 어느 쪽으로 움직이는지 손으로 확인해 보자.
CO₂ 1.00 mol을 0.500 L 용기에 넣고 300 K로 유지한다. 이상 기체 식과 판데르발스 식(\(a = 0.364\) Pa·m⁶/mol², \(b = 4.267\times10^{-5}\) m³/mol)으로 압력과 \(Z\)를 구하라.
풀이 보기
판데르발스 등온선을 여러 온도에서 그려 보면 놀라운 일이 생긴다. 높은 온도에서는 이상 기체처럼 부드러운 쌍곡선이지만, 온도를 낮추면 곡선에 S자 고리가 생긴다.
판데르발스 식에서 온도를 임계 온도보다 낮게 고정하고 큰 부피에서 출발해 부피를 줄여 간다. 식이 주는 압력은?
판데르발스 등온선과 임계점
고리가 처음 생기는 온도에서 극대와 극소가 한 점으로 합쳐진다. 그 점에서 등온선은 기울기도 0, 곡률도 0인 수평 변곡점을 가진다. 이 점이 임계점(critical point)이고, 그 온도가 임계 온도(critical temperature) \(T_\mathrm c\)다.
- 임계 온도 위에서 등온선은 단조 감소한다. 압축하면 기체는 점점 진해질 뿐 액체와 기체의 구분이 없다.
- 임계 온도에서는 극대와 극소가 한 점에서 만나 수평 변곡점이 된다. 두 조건을 연립하면 \(V_\mathrm c, T_\mathrm c, P_\mathrm c\)가 나온다.
- 임계 온도 아래에서는 고리가 생긴다. 극대·극소(스피노달) 사이의 빨간 구간은 \(\partial P/\partial V > 0\)인 불안정 상태다.
- 실제 계는 고리 대신 액체와 기체가 공존하는 수평선을 따른다. 고리의 위·아래 넓이가 같도록 높이를 정하는 것이 맥스웰 작도다.
임계점 조건 \((\partial P/\partial V_\mathrm m)_T = 0\), \((\partial^2P/\partial V_\mathrm m^2)_T = 0\)을 판데르발스 식에 적용하면
판데르발스 상수가 \(a\) = 0.137 Pa·m⁶/mol², \(b\) = 38.7인 기체는 \(T_\mathrm c\) = , \(P_\mathrm c\) = , \(V_\mathrm c\) = , 보일 온도 를 가진다. 25 °C에서 압축만으로 액화할 수 .
N₂의 판데르발스 상수 \(a = 0.137\) Pa·m⁶/mol², \(b = 3.87\times10^{-5}\) m³/mol로 \(T_\mathrm c, P_\mathrm c, V_\mathrm c\)와 \(Z_\mathrm c\)를 구하고 실측값 126.2 K, 34.0 bar, 90 cm³/mol, \(Z_\mathrm c = 0.289\)와 비교하라.
풀이 보기
임계 온도 위의 유체는 아무리 압축해도 액체와 기체로 나뉘지 않는다. 밀도는 액체만큼 높아질 수 있지만 경계면이 생기지 않는 이 상태를 초임계 유체(supercritical fluid)라 한다. 상온의 N₂, Ar, He, SiH₄(270 K)는 초임계이므로 실린더 안에서 액체가 없고, CO₂(304 K)와 NH₃(405 K)는 상온에서 압축하면 액화한다. 초임계 CO₂는 표면 장력이 없어 미세 패턴을 무너뜨리지 않고 말리는 건조 공정에 쓰인다(9장).
도전: 측정 자료에 판데르발스 등온선 맞추기
대응 상태 원리
판데르발스 식에 \(P = P_\mathrm rP_\mathrm c\), \(V_\mathrm m = V_\mathrm rV_\mathrm c\), \(T = T_\mathrm rT_\mathrm c\)를 넣고 임계 상수 공식을 쓰면 \(a\)와 \(b\)가 사라진다.
판데르발스 식이 정확하지 않은데도, 실제 기체들은 놀랄 만큼 이 성질을 따른다. 같은 \(T_\mathrm r\)와 \(P_\mathrm r\)에 있는 기체들은 거의 같은 \(Z\)를 가진다. 이것이 대응 상태 원리(principle of corresponding states)다. 분자 사이 퍼텐셜이 모양은 같고 에너지 척도 \(\varepsilon\)와 길이 척도 \(\sigma\)만 다르다면(예: 레너드–존스 퍼텐셜), 통계역학에서 이 원리가 엄밀히 나온다. \(kT_\mathrm c\propto\varepsilon\), \(V_\mathrm c\propto N_\mathrm A\sigma^3\)이기 때문이다. Ar, Kr, Xe, CH₄, N₂, O₂처럼 구에 가까운 무극성 분자는 아주 잘 맞고, 극성이 강하거나 길쭉한 분자(H₂O, NH₃)와 양자 효과가 큰 He, H₂는 덜 맞는다.
공학에서는 이 원리를 일반화 압축 인자 차트로 쓴다. \(Z\)를 \(T_\mathrm r\)별 곡선으로 \(P_\mathrm r\)에 대해 그려 두면, 어떤 기체든 임계 상수만 알면 \(Z\)를 읽을 수 있다. 정밀도를 높이려면 분자의 비구형성을 나타내는 세 번째 변수인 피처의 이심 인자 \(\omega\)를 더한다(\(Z = Z^{(0)} + \omega Z^{(1)}\)).
일반화 압축 인자 차트
메테인(\(T_\mathrm c = 190.6\) K, \(P_\mathrm c = 46.0\) bar)이 25 °C, 150 bar에 있다. 환산 변수를 구하고 일반화 차트로 \(Z\)와 몰부피를 추정하라.
풀이 보기
열팽창 계수와 등온 압축률
상태 방정식 \(V(T, P)\)를 알면 온도나 압력을 조금 바꿨을 때 부피가 얼마나 변하는지 알 수 있다. 부피의 전미분은
이상 기체 \(V = nRT/P\)에서는 \(\alpha = 1/T\), \(\kappa_T = 1/P\)다. 상온·1 bar 기체의 \(\alpha\approx3.4\times10^{-3}\) K⁻¹은 액체 물(\(2.6\times10^{-4}\) K⁻¹)의 13배, 실리콘(\(3\times2.6\times10^{-6}\) K⁻¹)의 430배다. 물은 특이하게 0–4 °C에서 \(\alpha < 0\)이다. 온도가 오르는데 부피가 줄어든다.
\(V\)가 \(T\)와 \(P\)의 함수이므로 세 편미분 사이에는 순환 관계 \((\partial P/\partial T)_V(\partial T/\partial V)_P(\partial V/\partial P)_T = -1\)이 성립한다. 이것을 \(\alpha\)와 \(\kappa_T\)로 쓰면
밸브 두 개 사이에 물이 꽉 찬 채 갇힌 배관이 5 데워지면(25 °C 물: \(\alpha = 2.57\times10^{-4}\) K⁻¹, \(\kappa_T = 4.52\times10^{-10}\) Pa⁻¹) 압력은 오른다. 기체(N₂, 1 bar)였다면 올랐을 것이다.
양쪽 밸브가 닫힌 초순수 배관 구간이 물로 완전히 차 있다. 실내 온도가 올라 물의 온도가 5.0 K 높아지면 압력은 얼마나 오를까? 배관은 늘어나지 않는다고 본다. 25 °C 물의 \(\alpha = 2.57\times10^{-4}\) K⁻¹, \(\kappa_T = 4.52\times10^{-10}\) Pa⁻¹이다.
풀이 보기
| 물질(25 °C 부근) | 선팽창 계수 αL [10⁻⁶ K⁻¹] | 등온 압축률 κT [10⁻¹¹ Pa⁻¹] |
|---|---|---|
| 실리콘 Si | 2.6 | 1.0 |
| 용융 실리카 SiO₂ | 0.5 | 2.7 |
| 구리 Cu | 16.5 | 0.73 |
| 알루미늄 Al | 23.1 | 1.3 |
| 텅스텐 W | 4.5 | 0.32 |
| 물(액체) | 86 (= α/3) | 45.2 |
| 이상 기체(1 bar) | 1120 (= 1/3T) | 10⁶ (= 1/P) |
고체의 열팽창은 원자 사이 퍼텐셜이 비대칭(비조화)이기 때문에 생긴다. 그 미시적 그림과 그뤼나이젠 관계는 고체물리 · 열적 성질에서 다룬다. 4장에서는 \(\alpha\)와 \(\kappa_T\)가 열용량의 차이 \(C_P - C_V = TV\alpha^2/\kappa_T\)를 정한다는 것을 증명한다. 막과 기판의 열팽창 계수가 다르면 증착 온도에서 식을 때 막에 응력이 생기고 웨이퍼가 휜다.
노광 장비는 웨이퍼와 레티클, 스테이지의 위치를 nm 단위로 맞춰야 한다. 실리콘의 선팽창 계수 \(2.6\times10^{-6}\) K⁻¹을 300 mm에 곱하면 1 K에 780 nm, 0.1 K에도 78 nm가 늘어난다. 그래서 노광 구역의 공조는 ±0.1 °C 이내로 온도를 묶고, 웨이퍼는 노광 전에 척 위에서 온도를 맞춘다(열평형 대기). 같은 계산이 CMP 패드, 식각 척, 프로브 카드의 정렬에도 쓰인다. → InfraBook · 공조: 온도·습도·압력
이 장에서 상태 방정식은 "평형 상태 사이의 관계"였다. 다음 장에서는 한 평형 상태에서 다른 평형 상태로 갈 때 주고받는 일과 열을 다룬다. → 2장 · 제1법칙과 엔탈피
핵심 정리
- 열평형과 제0법칙: 투열벽으로 맞댄 계는 거시 변수가 변하지 않는 열평형에 이르고, "열평형"은 전이적 관계다. 서로 열평형인 계들의 공통 이름표가 온도다. 분자 수준에서는 분자당 평균 운동 에너지가 같다는 뜻이다. 평형 상태는 몇 개의 상태 변수로 정해지고, 그 관계 \(f(P,V,T,n)=0\)이 상태 방정식이다.
- 이상 기체 온도: \(T = 273.16\ \mathrm K\times\lim_{P_\mathrm{tp}\to0}(P/P_\mathrm{tp})\)(부피 일정). 묽은 극한에서는 기체 종류와 무관하다. 2019년부터 켈빈은 \(k = 1.380649\times10^{-23}\) J/K로 정의한다. \(t[{}^\circ\mathrm C] = T[\mathrm K] - 273.15\).
- 이상 기체와 돌턴: 벽 충격량에서 \(PV = NkT = nRT\), \(n_V = P/kT\). 혼합 기체는 \(P = \sum P_i\), \(P_i = x_iP\). 같은 온도에서 압력은 분자 개수만 센다. 1 sccm = \(7.436\times10^{-7}\) mol/s(0 °C, 1 atm 기준). 혼합 유량의 분압은 \(P_i = (\dot n_i/\sum\dot n_j)P\).
- 실제 기체: \(Z = PV_\mathrm m/RT = 1 + B/V_\mathrm m + C/V_\mathrm m^2 + \cdots\approx1 + BP/RT\). \(B<0\)이면 인력, \(B>0\)이면 반발이 우세하고 \(B(T_\mathrm B) = 0\)인 보일 온도에서 바뀐다. 고압 실린더의 몰수는 \(n = PV/ZRT\)이다(Ar 150 bar에서 \(Z\approx0.93\)).
- 판데르발스 식: \(P = RT/(V_\mathrm m-b) - a/V_\mathrm m^2\), \(B = b - a/RT\), \(T_\mathrm B = a/Rb\). 임계점은 \(\partial P/\partial V = \partial^2P/\partial V^2 = 0\): \(V_\mathrm c = 3b\), \(T_\mathrm c = 8a/27Rb\), \(P_\mathrm c = a/27b^2\), \(Z_\mathrm c = 3/8\)(실측 0.23–0.30). \(T<T_\mathrm c\)의 고리 중 \(\partial P/\partial V>0\) 구간은 불안정하고, 실제로는 액체–기체 공존선으로 바뀐다(맥스웰 작도는 9장). \(T>T_\mathrm c\)이면 초임계 유체다.
- 대응 상태 원리: 환산 변수 \(T_\mathrm r = T/T_\mathrm c\), \(P_\mathrm r = P/P_\mathrm c\)가 같으면 기체가 달라도 \(Z\)가 거의 같다. 판데르발스 환산식 \((P_\mathrm r + 3/V_\mathrm r^2)(3V_\mathrm r - 1) = 8T_\mathrm r\)에는 물질 상수가 없다. 일반화 차트로 임계 상수만 알면 \(Z\)를 추정한다.
- 응답 계수: \(dV/V = \alpha\,dT - \kappa_T\,dP\), 이상 기체 \(\alpha = 1/T\), \(\kappa_T = 1/P\). \((\partial P/\partial T)_V = \alpha/\kappa_T\) — 갇힌 물은 1 K에 약 5.7 bar 오른다. 임계점에서 \(\kappa_T\to\infty\). 실리콘 \(\alpha_L = 2.6\times10^{-6}\) K⁻¹.
연습문제
\(R = 8.314\) J/(mol·K), \(k = 1.381\times10^{-23}\) J/K, \(N_\mathrm A = 6.022\times10^{23}\) mol⁻¹, 1 atm = 101 325 Pa, 1 Torr = 133.32 Pa, 0 °C = 273.15 K를 쓴다. 따로 말이 없으면 답이 ±1 % 이내면 정답이다.
RTP 장비가 웨이퍼를 1050 °C까지 가열한다. 이 온도를 켈빈으로 나타내라.
\(T = t + 273.15 = 1323.15\) K ≈ 1323 K.
열전대 하나를 웨이퍼 척에 대니 87.3 °C, 같은 열전대를 챔버 벽에 대니 87.3 °C를 읽었다. 제0법칙에서 곧바로 결론지을 수 있는 것은?
척과 열전대, 벽과 열전대가 각각 열평형이므로 제0법칙에 따라 척과 벽도 열평형, 곧 같은 온도다. 맞대도 알짜 열 흐름이 없다. 압력이나 열용량은 온도와 무관하다.
100 mTorr, 25 °C인 챔버 속 기체의 수밀도를 cm⁻³ 단위로 구하라.
\(P = 13.33\) Pa, \(n_V = 13.33/(1.381\times10^{-23}\times298.15) = 3.24\times10^{21}\) m⁻³ \(= 3.24\times10^{15}\) cm⁻³.
정적 기체 온도계의 삼중점 압력이 50.00 kPa였다. 이 온도계를 어떤 계에 넣으니 압력이 68.30 kPa가 되었다. 이 기체 양에서 읽은 온도 \(T^* = 273.16\,P/P_\mathrm{tp}\)는? (소수 둘째 자리까지)
\(T^* = 273.16\times68.30/50.00 = 373.14\) K. 이 값은 기체 양에 따라 조금씩 달라지므로, 정확한 온도를 얻으려면 기체 양을 바꿔 가며 잰 \(T^*\)를 \(P_\mathrm{tp}\to0\)으로 외삽한다.
N₂ 3.0 mol과 O₂ 1.0 mol을 섞은 기체의 전체 압력이 2.0 atm이다. O₂의 분압은?
\(x_{\mathrm{O_2}} = 1.0/4.0 = 0.25\), \(P_{\mathrm{O_2}} = 0.25\times2.0 = 0.50\) atm.
MFC가 N₂를 200 sccm(0 °C, 1 atm 기준) 흘린다. (a) 몰 유량(mol/s)과 (b) N₂ 1 mol을 흘리는 데 걸리는 시간(분)을 구하라.
(a) \(\dot n = (101325\times200\times10^{-6})/(8.314\times273.15)/60 = 1.49\times10^{-4}\) mol/s. (b) \(1/\dot n = 6724\) s = 112 분. 같은 뜻으로, 표준 상태 몰부피 22 414 cm³를 200 cm³/min으로 나눠도 112 분이다.
H₂에 1.0 % 희석한 B₂H₆ 혼합 기체 20 sccm과 순수 H₂ 2000 sccm을 챔버에 흘려 전체 압력을 10 Torr로 유지한다. B₂H₆의 분압을 mTorr로 구하라. 모든 기체가 같은 비율로 배기된다고 본다.
B₂H₆ 유량 \(= 0.010\times20 = 0.20\) sccm, 전체 2020 sccm. \(x = 0.20/2020 = 9.90\times10^{-5}\), \(P = 9.90\times10^{-5}\times10\) Torr \(= 9.90\times10^{-4}\) Torr = 0.990 mTorr.
CO₂ 1.00 mol을 300 K에서 0.250 L 용기에 넣었더니 압력이 69.8 bar였다. 압축 인자 \(Z\)는?
\(Z = PV/nRT = (69.8\times10^5)(2.50\times10^{-4})/(8.314\times300) = 1745/2494 = 0.700\). 인력 때문에 이상 기체보다 압력이 30 % 낮다.
300 K의 CO₂는 제2 비리얼 계수가 \(B = -122\) cm³/mol이다. 5.0 bar에서 저압 근사 \(Z\approx1+BP/RT\)로 \(Z\)와 몰부피 \(V_\mathrm m\)를 구하라.
\(BP/RT = (-122\times10^{-6})(5.0\times10^5)/(8.314\times300) = -0.0245\), \(Z = 0.976\). \(V_\mathrm m = ZRT/P = 0.976\times2494/5.0\times10^5 = 4.87\times10^{-3}\) m³/mol = 4.87 L/mol(이상 기체 4.99 L/mol).
N₂의 판데르발스 상수 \(a = 0.137\) Pa·m⁶/mol², \(b = 3.87\times10^{-5}\) m³/mol로 보일 온도를 구하라. (실측은 327 K다.)
\(B = b - a/RT = 0\)에서 \(T_\mathrm B = a/Rb = 0.137/(8.314\times3.87\times10^{-5}) = 426\) K. 판데르발스 식은 \(T_\mathrm B = 3.375\,T_\mathrm c\)를 주지만, 실제 기체는 \(2.4\)–\(2.7\,T_\mathrm c\)(N₂ 2.6, Ar 2.7, CO₂ 2.4)라 모형이 30 % 안팎 높게 준다.
N₂ 1.00 mol을 0.100 L 용기에 넣고 300 K로 유지한다. 판데르발스 식(\(a = 0.137\) Pa·m⁶/mol², \(b = 3.87\times10^{-5}\) m³/mol)과 이상 기체 식으로 압력을 각각 구하라.
\(V_\mathrm m = 1.00\times10^{-4}\) m³/mol. 반발 항 \(RT/(V_\mathrm m-b) = 2494.3/6.13\times10^{-5} = 4.069\times10^7\) Pa, 인력 항 \(a/V_\mathrm m^2 = 1.37\times10^7\) Pa. \(P = 2.70\times10^7\) Pa = 270 bar. 이상 기체는 \(2494.3/10^{-4} = 249\) bar. 이 밀도에서는 반발이 우세해 \(Z\approx1.08\)이다.
25 °C의 47 L Ar 실린더 압력이 150 bar에서 30 bar로 떨어졌다(둘 다 절대 압력). 두 압력에서의 압축 인자가 각각 0.925와 0.975라면, 처음 들어 있던 Ar의 몇 %가 남았는가?
\(\dfrac{n_2}{n_1} = \dfrac{P_2/Z_2}{P_1/Z_1} = \dfrac{30/0.975}{150/0.925} = \dfrac{30.77}{162.16} = 0.190\), 곧 19.0 %. 압력 비(20 %)로 계산하면 남은 양을 5 % 과대 평가한다.
CO₂의 실측 임계 상수 \(T_\mathrm c = 304.13\) K, \(P_\mathrm c = 73.77\) bar로 판데르발스 상수 \(a\)와 \(b\)를 구하라.
\(a = 27R^2T_\mathrm c^2/(64P_\mathrm c) = 27(8.314)^2(304.13)^2/(64\times7.377\times10^6) = 0.366\) Pa·m⁶/mol². \(b = RT_\mathrm c/(8P_\mathrm c) = 8.314\times304.13/(8\times7.377\times10^6) = 4.28\times10^{-5}\) m³/mol = 42.8 cm³/mol.
대응 상태 원리에 따르면 다음 중 압축 인자가 가장 비슷한 두 상태는? (임계 상수: Ar 150.7 K, 48.6 bar / CH₄ 190.6 K, 46.0 bar / CO₂ 304.1 K, 73.8 bar)
Ar 226 K, 97 bar는 \(T_\mathrm r = 1.50\), \(P_\mathrm r = 2.00\)이고 CH₄ 286 K, 92 bar도 \(T_\mathrm r = 1.50\), \(P_\mathrm r = 2.00\)이다. 같은 \(T, P\)라도 임계 상수가 다르면 환산 상태가 달라 \(Z\)가 다르다(CO₂ 300 K는 \(T_\mathrm r\approx0.99\)라 액화 근처다).
300 K, 2.0 bar의 이상 기체의 열팽창 계수 \(\alpha\)와 등온 압축률 \(\kappa_T\)를 구하라.
\(V = nRT/P\)에서 \(\alpha = (1/V)(nR/P) = 1/T = 3.33\times10^{-3}\) K⁻¹, \(\kappa_T = -(1/V)(-nRT/P^2) = 1/P = 1/(2.0\times10^5) = 5.0\times10^{-6}\) Pa⁻¹.
양쪽 밸브가 닫혀 물이 꽉 찬 배관이 25 °C에서 30 °C로 데워졌다. 배관이 늘어나지 않는다면 압력은 몇 bar 오르는가? \(\alpha = 2.57\times10^{-4}\) K⁻¹, \(\kappa_T = 4.52\times10^{-10}\) Pa⁻¹을 일정하다고 본다.
\(\Delta P = (\alpha/\kappa_T)\Delta T = (2.57\times10^{-4}/4.52\times10^{-10})\times5 = 2.84\times10^6\) Pa = 28.4 bar.
개념 확인 퀴즈
움직이지 않는 투열벽으로 나뉜 두 기체가 열평형에 이르렀다. 반드시 같은 것은?
정적 기체 온도계에서 \(P_\mathrm{tp}\to0\)으로 외삽하는 이유는?
상온 50 bar에서 어떤 기체의 \(Z = 0.85\)이다. 옳은 설명은?
25 °C에서 압축만으로 액화할 수 없는 기체는? (임계 온도: SiH₄ 270 K, CO₂ 304 K, NH₃ 405 K, H₂O 647 K)
판데르발스 등온선에서 \(\partial P/\partial V > 0\)인 구간에 대한 설명으로 옳은 것은?