← 로드맵 · ← 이전 L28 · 다음 L30 →

L29 · 텐서와 미분형식

Arfken 7판 4장 (텐서의 변환 법칙, 축약·직접곱, 의사텐서와 쌍대 텐서, 일반 좌표의 텐서: 계량 텐서·크리스토펠 기호·공변 미분, 미분형식의 미분과 적분)

이 레슨이 끝나면 할 수 있는 계산
  1. 위·아래 첨자로 식을 쓰고, 좌표를 바꿀 때 반변·공변 성분 바꾸기
  2. 회전에서 2계 텐서 성분 바꾸기($T'=RTR^T$), 직접곱·축약·대칭/반대칭 분해
  3. 반대칭 텐서와 축벡터를 쌍대로 오가기, 반전에서 의사텐서 판정하기
  4. 관성 텐서를 계산하고 주축 구하기
  5. 계량 텐서로 지표를 올리고 내리기, 크리스토펠 기호·공변 미분·공변 발산·측지선 방정식 계산하기
  6. 민코프스키 계량으로 4-벡터 내적을 계산하고 부스트에서 불변임을 확인하기
  7. 미분형식의 쐐기곱과 외미분을 계산하고, 형식 적분에서 야코비안 얻기
이 레슨의 관례 (책마다 다른 곳)

0. 외워 둘 공식

항목공식
반변(contravariant) 벡터$V'^i=\dfrac{\partial x'^i}{\partial x^j}V^j$ (좌표 미분 $dx^i$처럼 변환)
공변(covariant) 벡터$V'_i=\dfrac{\partial x^j}{\partial x'^i}V_j$ (기울기 $\partial_i f$처럼 변환)
2계 텐서$T'^{ij}=\dfrac{\partial x'^i}{\partial x^k}\dfrac{\partial x'^j}{\partial x^l}T^{kl}$ (지표마다 하나씩)
회전 (직교좌표)$A'_i=R_{ij}A_j$, $\ T'_{ij}=R_{ik}R_{jl}T_{kl}\iff T'=RTR^T$. $\operatorname{tr}T$, $\det T$ 불변
직접곱·축약$(A\otimes B)_{ij}=A_iB_j$ (계수가 더해짐), $\ T^i{}_i$ (계수가 2 줄어듦)
대칭·반대칭$S_{ij}=\tfrac12(T_{ij}+T_{ji})$, $\ A_{ij}=\tfrac12(T_{ij}-T_{ji})$
쌍대(dual)$C_i=\tfrac12\varepsilon_{ijk}A_{jk}\iff A_{ij}=\varepsilon_{ijk}C_k$
관성 텐서$I_{ij}=\sum m\,(r^2\delta_{ij}-x_ix_j)$
계량 텐서$ds^2=g_{ij}\,du^idu^j$, $\ g_{ij}=\dfrac{\partial x^a}{\partial u^i}\dfrac{\partial x^a}{\partial u^j}$, $\ g^{ik}g_{kj}=\delta^i_j$
올리고 내리기$V_i=g_{ij}V^j$, $\ V^i=g^{ij}V_j$, $\ |\mathbf V|^2=V^iV_i$
크리스토펠 기호$\Gamma^k_{ij}=\tfrac12g^{kl}\big(\partial_ig_{jl}+\partial_jg_{il}-\partial_lg_{ij}\big)$
공변 미분$\nabla_jV^i=\partial_jV^i+\Gamma^i_{jk}V^k$, $\ \nabla_jV_i=\partial_jV_i-\Gamma^k_{ji}V_k$
공변 발산·라플라시안$\nabla_iV^i=\dfrac{1}{\sqrt g}\partial_i\big(\sqrt g\,V^i\big)$, $\ \nabla^2f=\dfrac{1}{\sqrt g}\partial_i\big(\sqrt g\,g^{ij}\partial_jf\big)$, $\ g=\det(g_{ij})$
측지선$\ddot u^k+\Gamma^k_{ij}\dot u^i\dot u^j=0$ (점은 호의 길이 $s$로 미분)
민코프스키$a\cdot b=\eta_{\mu\nu}a^\mu b^\nu=a^0b^0-a^1b^1-a^2b^2-a^3b^3$
미분형식$dx\wedge dy=-dy\wedge dx$, $\ dx\wedge dx=0$, $\ d(f\,dx^I)=df\wedge dx^I$, $\ d^2=0$

1. 지표 표기와 변환 법칙

레시피
  1. 위 첨자 $V^i$는 반변 성분, 아래 첨자 $V_i$는 공변 성분입니다. 좌표 자체는 위 첨자 $x^i$로 씁니다.
  2. 식을 쓰면 먼저 지표를 셉니다. 자유 지표는 양변에서 같은 자리(위 또는 아래)에 한 번씩 있어야 합니다. 합하는 지표(더미 지표)는 위·아래 한 쌍입니다. 한 항에 같은 지표가 세 번 나오면 틀린 식입니다.
  3. 반변에는 "새 좌표를 옛 좌표로 미분한" $\dfrac{\partial x'^i}{\partial x^j}$를 곱하고, 공변에는 "옛 좌표를 새 좌표로 미분한" $\dfrac{\partial x^j}{\partial x'^i}$를 곱합니다. 두 행렬은 서로 역행렬입니다.
  4. 직교좌표 $(x,y)$ → 평면 극좌표 $(r,\theta)$에서 쓰는 편미분: $$\begin{aligned}\frac{\partial r}{\partial x}&=\cos\theta,\quad \frac{\partial r}{\partial y}=\sin\theta,\\ \frac{\partial\theta}{\partial x}&=-\frac{\sin\theta}{r},\quad \frac{\partial\theta}{\partial y}=\frac{\cos\theta}{r},\\ \frac{\partial x}{\partial r}&=\cos\theta,\quad \frac{\partial y}{\partial r}=\sin\theta,\\ \frac{\partial x}{\partial\theta}&=-r\sin\theta,\quad \frac{\partial y}{\partial\theta}=r\cos\theta.\end{aligned}$$
  5. 검산: 스칼라 $V^iW_i$는 어느 좌표에서 계산해도 같아야 합니다.
왜 되나? (증명 스케치)

반변의 원형은 좌표의 미분입니다. 연쇄법칙으로 $dx'^i=\dfrac{\partial x'^i}{\partial x^j}dx^j$. 속도 $dx^i/dt$ 처럼 "$dx$를 닮은 것"은 모두 이렇게 변환합니다.

공변의 원형은 기울기입니다. 연쇄법칙으로 $\dfrac{\partial f}{\partial x'^i}=\dfrac{\partial x^j}{\partial x'^i}\dfrac{\partial f}{\partial x^j}$. 이번에는 옛 좌표를 새 좌표로 미분한 것이 붙습니다.

2계 텐서는 $dx^i\,dx^j$ 나 $V^iW^j$ 처럼 지표가 두 개인 것이라, 지표마다 같은 인자를 하나씩 붙입니다. 두 변환 행렬이 서로 역행렬인 것도 연쇄법칙입니다: $\dfrac{\partial x'^i}{\partial x^j}\dfrac{\partial x^j}{\partial x'^k}=\dfrac{\partial x'^i}{\partial x'^k}=\delta^i_k$.

예제 1-1 회전 흐름 $\mathbf V=(-y,\ x)$의 극좌표 반변 성분과 공변 성분

반변 성분: 직교좌표에서 $V^x=-y=-r\sin\theta$, $V^y=x=r\cos\theta$ 입니다.

$$\begin{aligned}V^r&=\frac{\partial r}{\partial x}V^x+\frac{\partial r}{\partial y}V^y\\&=\cos\theta(-r\sin\theta)+\sin\theta(r\cos\theta)=0,\\ V^\theta&=\frac{\partial \theta}{\partial x}V^x+\frac{\partial \theta}{\partial y}V^y\\&=\Big(-\frac{\sin\theta}{r}\Big)(-r\sin\theta)+\frac{\cos\theta}{r}(r\cos\theta)=1.\end{aligned}$$

공변 성분: 직교좌표에서는 위·아래가 같으므로 $V_x=-y$, $V_y=x$ 입니다.

$$\begin{aligned}V_r&=\frac{\partial x}{\partial r}V_x+\frac{\partial y}{\partial r}V_y\\&=\cos\theta(-r\sin\theta)+\sin\theta(r\cos\theta)=0,\\ V_\theta&=\frac{\partial x}{\partial \theta}V_x+\frac{\partial y}{\partial \theta}V_y\\&=(-r\sin\theta)(-r\sin\theta)+(r\cos\theta)(r\cos\theta)=r^2.\end{aligned}$$

검산: $V^iV_i=0\cdot0+1\cdot r^2=r^2$. 직교좌표에서도 $V^xV_x+V^yV_y=y^2+x^2=r^2$ 입니다.

같은 화살표인데 $\theta$ 성분이 세 가지입니다. $V^\theta=1$, $V_\theta=r^2$, 그리고 L09식 단위벡터 성분은 $\mathbf V=r\,\hat{\boldsymbol\theta}$ 이므로 $F_\theta=r$ 입니다. 관계는 $F_\theta=h_\theta V^\theta=V_\theta/h_\theta$ ($h_\theta=r$) 입니다. 6절과 8절에서 다시 씁니다.

2. 회전과 2계 텐서

레시피 (직교좌표, 회전 행렬 $R$: $R^TR=I$, $\det R=1$)
  1. 벡터: $A'_i=R_{ij}A_j$, 행렬로 $A'=RA$.
  2. 2계 텐서: 지표마다 $R$이 하나씩 붙습니다. $T'_{ij}=R_{ik}R_{jl}T_{kl}$. 행렬로 쓰면 두 번째 $R$은 전치되어 오른쪽에 붙습니다: $T'=RTR^T$.
  3. 계산 순서: $M=RT$를 먼저 구하고, 그다음 $T'=MR^T$를 구합니다.
  4. 검산: $\operatorname{tr}T'=\operatorname{tr}T$ (이유: $T'_{ii}=R_{ik}R_{il}T_{kl}=\delta_{kl}T_{kl}=T_{kk}$), $\det T'=\det T$, 대칭 텐서는 돌려도 대칭입니다.
예제 2-1 $T=\begin{pmatrix}1&2&0\\0&3&0\\0&0&5\end{pmatrix}$ 를 $z$축 둘레 $\alpha=45^\circ$ 회전으로 변환하기

$\cos45^\circ=\sin45^\circ=\tfrac1{\sqrt2}$ 이므로

$$R=\frac1{\sqrt2}\begin{pmatrix}1&1&0\\-1&1&0\\0&0&\sqrt2\end{pmatrix}.$$

① $M=RT$: 1행은 $T$의 (1행+2행), 2행은 (2행−1행)이고, 둘 다 $\tfrac1{\sqrt2}$배입니다. 3행은 그대로입니다.

$$M=\frac1{\sqrt2}\begin{pmatrix}1&5&0\\-1&1&0\\0&0&5\sqrt2\end{pmatrix}.$$

② $T'=MR^T$: $R^T$를 오른쪽에 곱하면 열에 같은 조작을 합니다. 새 1열은 (1열+2열), 새 2열은 (2열−1열)이고 $\tfrac1{\sqrt2}$배입니다.

$$T'=\frac12\begin{pmatrix}1+5&5-1&0\\-1+1&1+1&0\\0&0&10\end{pmatrix}=\begin{pmatrix}3&2&0\\0&1&0\\0&0&5\end{pmatrix}.$$

검산: $\operatorname{tr}T=1+3+5=9=3+1+5=\operatorname{tr}T'$. $\det T=15=\det T'$ (둘 다 위삼각 행렬이라 대각 성분의 곱).

덤: 반대칭 부분 $\tfrac12(T-T^T)$는 회전 전후 모두 $\begin{pmatrix}0&1&0\\-1&0&0\\0&0&0\end{pmatrix}$ 입니다. 4절에서 보듯 이것은 $z$ 방향 축벡터라서 $z$축 둘레 회전에 변하지 않습니다.

그림 표기로 보기: G3 6절 예제 6-3 — 회전 화살촉: 대각합 불변 = 화살촉 쌍소멸, 반대칭 부분 = 축벡터의 회전.

3. 직접곱, 축약, 대칭·반대칭 분해

레시피
  1. 직접곱(direct product): 성분끼리 곱해 새 텐서를 만듭니다. $T_{ij}=a_ib_j$ (행렬로 $\mathbf a\mathbf b^T$). 계수(rank)는 더해집니다: $1+1=2$. 일반 좌표에서는 $a^ib_j$ 처럼 위·아래가 섞인 혼합 텐서도 됩니다.
  2. 축약(contraction): 위 지표 하나와 아래 지표 하나를 같게 놓고 합합니다. 계수가 2 줄어듭니다. $T^i{}_i$ 는 대각합(스칼라), $a^ib_i$ 는 내적입니다. 직교좌표에서는 그냥 $T_{ii}$ 입니다.
  3. 분해: $T_{ij}=S_{ij}+A_{ij}$, $S=\tfrac12(T+T^T)$, $A=\tfrac12(T-T^T)$. 3차원에서 $S$는 독립 성분 6개, $A$는 3개입니다. 대칭·반대칭 성질은 좌표를 바꿔도 유지됩니다.
  4. (더 쪼개기) $S_{ij}=\tfrac13(\operatorname{tr}T)\,\delta_{ij}+\big(S_{ij}-\tfrac13(\operatorname{tr}T)\,\delta_{ij}\big)$. 대각합 부분 1개 + 대각합이 0인 대칭 부분 5개 + 반대칭 3개 = 9개입니다.
  5. 대칭과 반대칭을 완전히 축약하면 0입니다: $S_{ij}A_{ij}=0$.
왜 되나? (증명 스케치)

축약이 스칼라인 이유: 1절의 변환 법칙을 넣으면

$$V'^iW'_i=\frac{\partial x'^i}{\partial x^j}\frac{\partial x^k}{\partial x'^i}\,V^jW_k.$$

가운데에서 $i$로 합한 것은 연쇄법칙으로 $\dfrac{\partial x^k}{\partial x^j}=\delta^k_j$ 이므로 $V'^iW'_i=V^jW_j$. 위 지표의 인자와 아래 지표의 인자가 서로 상쇄됩니다. 혼합 텐서의 대각합 $T^i{}_i$ 도 같은 이유로 스칼라입니다.

그래서 축약은 위 하나·아래 하나의 짝이어야 합니다. 위 지표끼리 $V^iW^i$ 를 합하면 두 인자가 모두 $\partial x'/\partial x$ 라 상쇄되지 않습니다. 직교좌표의 회전에서만 $R^{-1}=R^T$ 라서 위·아래 구별 없이 괜찮습니다.

대칭·반대칭이 좌표와 무관한 것도 같은 이유입니다: $T'^{ij}$ 와 $T'^{ji}$ 는 같은 인자로 변환되므로 $T^{ij}=\pm T^{ji}$ 가 그대로 유지됩니다.

예제 3-1 $\mathbf a=(1,2,0)$, $\mathbf b=(3,-1,2)$의 직접곱 $T_{ij}=a_ib_j$를 축약하고 분해하기 $$T=\begin{pmatrix}1\\2\\0\end{pmatrix}\begin{pmatrix}3&-1&2\end{pmatrix}=\begin{pmatrix}3&-1&2\\6&-2&4\\0&0&0\end{pmatrix}.$$

축약: $T_{ii}=3-2+0=1$. 내적 $\mathbf a\cdot\mathbf b=3-2+0=1$ 과 같습니다.

분해: $T^T$와 더하고 빼서 반으로 나눕니다.

$$S=\begin{pmatrix}3&\tfrac52&1\\\tfrac52&-2&2\\1&2&0\end{pmatrix},$$ $$A=\begin{pmatrix}0&-\tfrac72&1\\\tfrac72&0&2\\-1&-2&0\end{pmatrix}.$$

검산: $S+A=T$. 그리고 $S_{ij}A_{ij}$ (성분끼리 곱해서 전부 더하기)는 $-\tfrac{35}{4}+1+\tfrac{35}{4}+4-1-4=0$ 입니다.

반대칭 부분의 정체 (4절 미리보기): $C_i=\tfrac12\varepsilon_{ijk}A_{jk}$ 는 $(A_{23},A_{31},A_{12})=(2,-1,-\tfrac72)$ 입니다. 한편

$$\begin{aligned}(\mathbf a\times\mathbf b)_x&=2\cdot2-0\cdot(-1)=4,\\(\mathbf a\times\mathbf b)_y&=0\cdot3-1\cdot2=-2,\\(\mathbf a\times\mathbf b)_z&=1\cdot(-1)-2\cdot3=-7.\end{aligned}$$

즉 $\mathbf C=\tfrac12\,\mathbf a\times\mathbf b$ 입니다. 외적은 직접곱의 반대칭 부분입니다.

그림 표기로 보기: G3 2절 예제 2-1 — 1 + 3 + 5 분해를 대칭화·반대칭화 막대로, 성분 제곱합 검산.

4. 의사텐서와 쌍대

레시피
  1. 의사텐서(pseudotensor)는 변환할 때 $\det R$가 하나 더 붙습니다. 의사벡터는 $C'_i=(\det R)\,R_{ij}C_j$. 회전($\det R=1$)에서는 보통 텐서와 구별되지 않고, 반전·거울 반사($\det R=-1$)에서 차이가 납니다.
  2. 반전(inversion) $R=-I$: 보통 벡터(극벡터, polar vector)는 부호가 바뀝니다. 보통 2계 텐서는 $(-1)^2=1$이라 그대로입니다. 의사벡터(축벡터, axial vector)는 $(\det R)(-1)=+1$이라 그대로입니다.
  3. $\varepsilon_{ijk}$는 모든 직교좌표계에서 같은 값을 갖는 3계 의사텐서입니다.
  4. 쌍대: 반대칭 2계 텐서 $A_{ij}$ ↔ 축벡터 $C_i$. $$\begin{aligned}C_i&=\tfrac12\varepsilon_{ijk}A_{jk}\ \Rightarrow\ \mathbf C=(A_{23},\,A_{31},\,A_{12}),\\ A&=\begin{pmatrix}0&C_3&-C_2\\-C_3&0&C_1\\C_2&-C_1&0\end{pmatrix}.\end{aligned}$$
  5. 쓸모: $A\mathbf v=\mathbf v\times\mathbf C$ (성분으로 $A_{ij}v_j=\varepsilon_{ijk}v_jC_k$). 외적 $\mathbf a\times\mathbf b$, 각운동량 $\mathbf L=\mathbf r\times\mathbf p$, 각속도 $\boldsymbol\omega$, 자기장 $\mathbf B$가 축벡터입니다.

그림 표기로 보기: G1 연습 D3 — 반대칭 부분의 쌍대 ½εijkAjk = ½a×b를 ε 점과 팔 꼬기(clank)로 읽기

예제 4-1 반대칭 행렬 $A=\begin{pmatrix}0&3&-2\\-3&0&1\\2&-1&0\end{pmatrix}$ 의 쌍대 벡터, 그리고 반전에서의 행동

① $\mathbf C=(A_{23},A_{31},A_{12})=(1,\,2,\,3)$.

② 확인: $\mathbf v=(1,0,0)$이면 $A\mathbf v$는 $A$의 1열 $(0,-3,2)$ 입니다. 한편 $\mathbf v\times\mathbf C=(0\cdot3-0\cdot2,\ 0\cdot1-1\cdot3,\ 1\cdot2-0\cdot1)=(0,-3,2)$. 같습니다.

③ 반전: $\mathbf r=(1,0,0)$, $\mathbf p=(0,1,0)$ 이면 $\mathbf L=\mathbf r\times\mathbf p=(0,0,1)$. 반전 뒤에는 $\mathbf r'=(-1,0,0)$, $\mathbf p'=(0,-1,0)$ 이고

$$\begin{aligned}\mathbf L'&=\mathbf r'\times\mathbf p'=\big(0,\ 0,\ (-1)(-1)-0\cdot0\big)\\&=(0,0,1)=\mathbf L.\end{aligned}$$

보통 벡터라면 $R\mathbf L=-\mathbf L$ 이어야 하는데 그대로입니다. 의사벡터 법칙 $(\det R)R\mathbf L=(-1)(-1)\mathbf L=\mathbf L$ 과 맞습니다. 반대칭 텐서 $L_{ij}=x_ip_j-x_jp_i$ 로 보면 $(-1)^2=1$ 이라 당연합니다.

5. 관성 텐서와 주축

레시피 (L05의 대각화 절차를 그대로 씁니다)
  1. 점질량: $I_{ij}=\sum_n m_n\big(r_n^2\delta_{ij}-x_{n,i}\,x_{n,j}\big)$. 연속체: $I_{ij}=\int\rho\,(r^2\delta_{ij}-x_ix_j)\,dV$.
  2. 대각 성분은 그 축까지 거리의 제곱: $I_{xx}=\sum m(y^2+z^2)$. 비대각 성분은 부호가 $-$: $I_{xy}=-\sum m\,xy$.
  3. $I$는 실대칭 행렬입니다. 고유값이 주관성모멘트(principal moment), 고유벡터가 주축(principal axis)입니다.
  4. 쓰임: $\mathbf L=I\boldsymbol\omega$, $T=\tfrac12\boldsymbol\omega^TI\boldsymbol\omega$, 단위벡터 $\hat{\mathbf n}$ 축에 대한 관성모멘트 $I_n=\hat{\mathbf n}^TI\hat{\mathbf n}$.
  5. 검산: $\operatorname{tr}I=2\sum m r^2$. 평행축 정리: $I_{ij}(\text{점 O})=I_{ij}(\text{질량중심})+M(d^2\delta_{ij}-d_id_j)$, $\mathbf d$는 O에서 본 질량중심 위치.
왜 되나? (증명 스케치)

강체가 각속도 $\boldsymbol\omega$ 로 돌면 각 점의 속도는 $\boldsymbol\omega\times\mathbf r$ 이고, 각운동량은 BAC−CAB(L07)로

$$\begin{aligned}\mathbf L&=\sum m\,\mathbf r\times(\boldsymbol\omega\times\mathbf r)\\&=\sum m\big(r^2\boldsymbol\omega-(\mathbf r\cdot\boldsymbol\omega)\mathbf r\big).\end{aligned}$$

성분으로 쓰면 $L_i=\sum m(r^2\delta_{ij}-x_ix_j)\,\omega_j=I_{ij}\omega_j$. 관성 텐서의 정의(비대각의 $-$ 부호 포함)가 여기서 나옵니다.

주축은 그 축으로 돌릴 때 $\mathbf L$ 이 $\boldsymbol\omega$ 와 나란한 축입니다: $I\boldsymbol\omega=\lambda\boldsymbol\omega$. 그대로 고유값 문제입니다. $I$ 가 실대칭이므로 고유값은 실수이고 주축은 서로 수직으로 고를 수 있습니다(L05). 주축 좌표에서는 $I$ 가 대각이라 회전 에너지가 $\tfrac12(I_1\omega_1^2+I_2\omega_2^2+I_3\omega_3^2)$ 로 교차항 없이 쓰입니다.

예제 5-1 질량 $m$인 점 세 개 $(1,1,0)$, $(-1,-1,0)$, $(0,0,1)$ (길이 단위 생략)의 원점 관성 텐서와 주축

① 점마다 $r^2\delta_{ij}-x_ix_j$ 를 씁니다. 대각은 $r^2-x_i^2$, 비대각은 $-x_ix_j$ 입니다. $(1,1,0)$은 $r^2=2$ 이므로 대각이 $(2-1,\ 2-1,\ 2-0)$, $xy$ 자리가 $-1$ 입니다. $(0,0,1)$은 $r^2=1$ 입니다.

$$(1,1,0):\ \begin{pmatrix}1&-1&0\\-1&1&0\\0&0&2\end{pmatrix},$$ $$(0,0,1):\ \operatorname{diag}(1,1,0).$$

$(-1,-1,0)$은 $xy=(-1)(-1)=1$ 이라 $(1,1,0)$과 같습니다.

② 더하면

$$I=m\begin{pmatrix}3&-2&0\\-2&3&0\\0&0&4\end{pmatrix}.$$

검산: $\operatorname{tr}I=10m=2m(2+2+1)$.

③ $z$는 따로 떨어져 있으므로 $4m$, 축 $\hat{\mathbf z}$. 남은 2×2 블록: $\lambda^2-6\lambda+5=0\Rightarrow\lambda=1,\,5$. $\lambda=1$이면 $(3-1)v_1-2v_2=0\Rightarrow\mathbf v\propto(1,1)$. $\lambda=5$이면 $\mathbf v\propto(1,-1)$.

④ 주관성모멘트 $m,\ 4m,\ 5m$. 주축 $\tfrac{1}{\sqrt2}(1,1,0)$, $(0,0,1)$, $\tfrac1{\sqrt2}(1,-1,0)$.

물리 확인: $(1,1,0)$ 방향 축 위에는 두 점이 놓여 기여가 0이고, $(0,0,1)$의 점만 거리 1이라 $I=m$ 입니다.

예제 5-2 균일한 정육면체 (질량 $M$, 한 변 $a$, $0\le x,y,z\le a$): 꼭짓점(원점)에 대한 관성 텐서와 주축

밀도 $\rho=M/a^3$. 각 변수의 적분 $\int_0^a dx=a$, $\int_0^ax\,dx=\tfrac{a^2}{2}$, $\int_0^ax^2dx=\tfrac{a^3}{3}$ 만 씁니다.

$$\begin{aligned}I_{xx}&=\rho\int(y^2+z^2)\,dV=\rho\Big(a\cdot\frac{a^3}{3}\cdot a+a\cdot a\cdot\frac{a^3}{3}\Big)\\&=\frac23\rho a^5=\frac23Ma^2,\\ I_{xy}&=-\rho\int xy\,dV=-\rho\cdot\frac{a^2}{2}\cdot\frac{a^2}{2}\cdot a=-\frac14Ma^2.\end{aligned}$$

대칭성으로 대각은 모두 $\tfrac23$, 비대각은 모두 $-\tfrac14$ 입니다.

$$I=\frac{Ma^2}{12}\begin{pmatrix}8&-3&-3\\-3&8&-3\\-3&-3&8\end{pmatrix}.$$

주축: 대칭성으로 대각선 $(1,1,1)$을 짐작합니다. 각 행의 합이 $8-3-3=2$ 이므로 $I(1,1,1)^T=\tfrac{2Ma^2}{12}(1,1,1)^T$, 즉 $\lambda=\tfrac{Ma^2}{6}$. 여기에 수직인 $(1,-1,0)$을 넣으면 $(8+3,\,-3-8,\,-3+3)=11(1,-1,0)$, 즉 $\lambda=\tfrac{11Ma^2}{12}$ (수직인 방향은 모두 이 값, 겹친 고유값).

검산 ①: $\tfrac{2+11+11}{12}Ma^2=2Ma^2=\operatorname{tr}I$.

검산 ② (평행축 정리): 중심에서는 $I=\tfrac{Ma^2}{6}\delta_{ij}$ 입니다. $\mathbf d=\tfrac a2(1,1,1)$, $d^2=\tfrac34a^2$ 이므로 대각 $\tfrac16+\tfrac34-\tfrac14=\tfrac23$, 비대각 $0-\tfrac14=-\tfrac14$ (단위 $Ma^2$). 맞습니다. 대각선 축은 중심을 지나므로 그 축의 값 $\tfrac{Ma^2}{6}$도 중심 값과 같습니다.

6. 계량 텐서

레시피
  1. 직교좌표 $x^a$를 새 좌표 $u^i$의 식으로 쓰고, 편미분 벡터 $\mathbf e_i=\partial\mathbf x/\partial u^i$를 구합니다. 그러면 $g_{ij}=\mathbf e_i\cdot\mathbf e_j$ 입니다. 또는 $ds^2$을 전개해 $du^idu^j$의 계수를 읽습니다. 이때 교차항 $du^1du^2$의 계수는 $2g_{12}$ 입니다.
  2. $g^{ij}$는 행렬 $(g_{ij})$의 역행렬입니다. 대각이면 $g^{ii}=1/g_{ii}$. 2×2는 $\begin{pmatrix}p&q\\q&s\end{pmatrix}^{-1}=\dfrac1{ps-q^2}\begin{pmatrix}s&-q\\-q&p\end{pmatrix}$.
  3. $g=\det(g_{ij})$. 부피(넓이) 요소는 $\sqrt g\,du^1du^2\cdots$ 입니다. L09의 $h_1h_2h_3$와 같습니다.
  4. 내리기 $V_i=g_{ij}V^j$, 올리기 $V^i=g^{ij}V_j$. 길이 $|\mathbf V|^2=g_{ij}V^iV^j=V^iV_i$.
  5. 대각 계량에서는 $h_i=\sqrt{g_{ii}}$ 이고, 단위벡터 성분은 $F_i=h_iV^i=V_i/h_i$ 입니다.
예제 6-1 평면 극좌표, 구면좌표, 2-구 표면의 계량 텐서와 지표 올리고 내리기

(a) 평면 극좌표 $x=r\cos\theta$, $y=r\sin\theta$. $\mathbf e_r=(\cos\theta,\ \sin\theta)$, $\mathbf e_\theta=(-r\sin\theta,\ r\cos\theta)$. 서로 내적하면 $g_{rr}=1$, $g_{r\theta}=0$, $g_{\theta\theta}=r^2\sin^2\theta+r^2\cos^2\theta=r^2$.

$$\begin{aligned}g_{ij}&=\begin{pmatrix}1&0\\0&r^2\end{pmatrix},\quad g^{ij}=\begin{pmatrix}1&0\\0&1/r^2\end{pmatrix},\\ \sqrt g&=r.\end{aligned}$$

(b) 구면좌표: L09 예제 1-1의 편미분 벡터는 서로 수직이고 길이가 $1,\ r,\ r\sin\theta$ 입니다.

$$\begin{aligned}g_{ij}&=\operatorname{diag}(1,\ r^2,\ r^2\sin^2\theta),\\ g^{ij}&=\operatorname{diag}\Big(1,\ \frac1{r^2},\ \frac1{r^2\sin^2\theta}\Big),\\ \sqrt g&=r^2\sin\theta.\end{aligned}$$

(c) 반지름 $a$인 2-구 표면: $r=a$로 고정하면 $dr=0$ 이고 $ds^2=a^2d\theta^2+a^2\sin^2\theta\,d\varphi^2$. 좌표는 $(\theta,\varphi)$ 두 개입니다.

$$\begin{aligned}g_{ij}&=\operatorname{diag}(a^2,\ a^2\sin^2\theta),\\ g^{ij}&=\operatorname{diag}\Big(\frac1{a^2},\ \frac1{a^2\sin^2\theta}\Big),\quad\sqrt g=a^2\sin\theta.\end{aligned}$$

검산: 넓이 $\int_0^{2\pi}\!\int_0^\pi a^2\sin\theta\,d\theta\,d\varphi=4\pi a^2$.

(d) 올리고 내리기 (2-구, $a=1$, $\theta=\pi/6$ 이라 $\sin^2\theta=\tfrac14$): $V^i=(V^\theta,V^\varphi)=(1,\,2)$.

$$V_\theta=g_{\theta\theta}V^\theta=1,\qquad V_\varphi=g_{\varphi\varphi}V^\varphi=\tfrac14\cdot2=\tfrac12.$$

길이: $|\mathbf V|^2=V^iV_i=1\cdot1+2\cdot\tfrac12=2$. 되올리기: $V^\varphi=g^{\varphi\varphi}V_\varphi=4\cdot\tfrac12=2$. 단위벡터 성분으로도 $F_\theta=1\cdot1=1$, $F_\varphi=\tfrac12\cdot2=1$ 이라 $1^2+1^2=2$ 로 같습니다.

7. 크리스토펠 기호

레시피
  1. 공식: $\Gamma^k_{ij}=\tfrac12g^{kl}\big(\partial_ig_{jl}+\partial_jg_{il}-\partial_lg_{ij}\big)$. $\Gamma^k_{ij}=\Gamma^k_{ji}$ 이므로 2차원에서는 6개, 3차원에서는 18개만 구하면 됩니다.
  2. 대각 계량의 지름길 (합 없음, $i\ne j$): $$\begin{aligned}\Gamma^i_{ii}&=\frac{\partial_ig_{ii}}{2g_{ii}},\\ \Gamma^i_{ij}=\Gamma^i_{ji}&=\frac{\partial_jg_{ii}}{2g_{ii}},\\ \Gamma^j_{ii}&=-\frac{\partial_jg_{ii}}{2g_{jj}},\end{aligned}$$ 세 지표가 모두 다르면 0입니다.
  3. 어떤 $g_{ii}$가 어떤 좌표에 의존하는지 먼저 적습니다. 상수인 성분의 미분은 모두 0이라 대부분의 $\Gamma$가 바로 지워집니다.
  4. 검산: 평면 직교좌표에서는 모두 0입니다. 극좌표는 평면인데도 $\Gamma\ne0$ 입니다. $\Gamma$는 텐서가 아닙니다.
왜 되나? (증명 스케치)

기저벡터 $\mathbf e_j=\partial\mathbf x/\partial u^j$ 의 미분을 다시 기저로 펼친 계수를 $\Gamma$ 로 정의합니다: $\partial_i\mathbf e_j=\Gamma^k_{ij}\mathbf e_k$. $\partial_i\mathbf e_j=\partial_i\partial_j\mathbf x$ 는 $i,j$ 에 대칭이므로 $\Gamma^k_{ij}=\Gamma^k_{ji}$.

$g_{jl}=\mathbf e_j\cdot\mathbf e_l$ 을 곱의 미분으로 미분하고 $\Gamma_{l,ij}\equiv\mathbf e_l\cdot\partial_i\mathbf e_j=g_{lk}\Gamma^k_{ij}$ 로 쓰면

$$\partial_ig_{jl}=\Gamma_{l,ij}+\Gamma_{j,il}.$$

지표를 돌려 쓴 세 식 $\partial_ig_{jl}$, $\partial_jg_{il}$, $\partial_lg_{ij}$ 에서 앞의 둘을 더하고 셋째를 빼면, 뒤 두 지표의 대칭성 때문에 나머지가 짝지어 지워지고 $\Gamma_{l,ij}$ 만 두 번 남습니다.

$$\partial_ig_{jl}+\partial_jg_{il}-\partial_lg_{ij}=2\Gamma_{l,ij}.$$

$g^{kl}$ 을 곱해 지표를 올리면 레시피 1의 공식입니다. 요점은 계량만 알면 $\mathbf x(u)$ 를 몰라도 $\Gamma$ 를 구할 수 있다는 것입니다.

예제 7-1 평면 극좌표와 2-구의 크리스토펠 기호

(a) 평면 극좌표: $g_{rr}=1$, $g_{\theta\theta}=r^2$. $r$에 의존하는 것은 $g_{\theta\theta}$뿐입니다.

$$\begin{aligned}\Gamma^r_{\theta\theta}&=\tfrac12g^{rr}\big(\partial_\theta g_{\theta r}+\partial_\theta g_{\theta r}-\partial_rg_{\theta\theta}\big)\\&=\tfrac12(0+0-2r)=-r,\\ \Gamma^\theta_{r\theta}&=\tfrac12g^{\theta\theta}\big(\partial_rg_{\theta\theta}+\partial_\theta g_{r\theta}-\partial_\theta g_{r\theta}\big)\\&=\tfrac12\cdot\frac1{r^2}\cdot2r=\frac1r=\Gamma^\theta_{\theta r}.\end{aligned}$$

나머지 $\Gamma^r_{rr}$, $\Gamma^r_{r\theta}$, $\Gamma^\theta_{rr}$, $\Gamma^\theta_{\theta\theta}$는 상수인 $g_{rr}$의 미분이나 $\partial_\theta g_{\theta\theta}=0$ 만 들어 있어서 0입니다. (8절 측지선 방정식에서 $\Gamma^r_{\theta\theta}=-r$ 이 구심 가속도 항 $-r\dot\theta^2$ 을 만듭니다.)

(b) 반지름 $a$인 2-구: $g_{\theta\theta}=a^2$ (상수), $g_{\varphi\varphi}=a^2\sin^2\theta$ ($\theta$에 의존).

$$\begin{aligned}\Gamma^\theta_{\varphi\varphi}&=-\frac{\partial_\theta g_{\varphi\varphi}}{2g_{\theta\theta}}=-\frac{2a^2\sin\theta\cos\theta}{2a^2}\\&=-\sin\theta\cos\theta,\\ \Gamma^\varphi_{\theta\varphi}=\Gamma^\varphi_{\varphi\theta}&=\frac{\partial_\theta g_{\varphi\varphi}}{2g_{\varphi\varphi}}=\frac{2a^2\sin\theta\cos\theta}{2a^2\sin^2\theta}\\&=\cot\theta.\end{aligned}$$

나머지 4개는 0입니다. 반지름 $a$는 약분되어 사라집니다.

8. 공변 미분, 공변 발산, 측지선

레시피 A — 공변 미분(covariant derivative)
  1. 곡선좌표에서는 기저벡터 자체가 위치마다 바뀝니다: $\partial_j\mathbf e_k=\Gamma^i_{jk}\mathbf e_i$ (7절). 그래서 $\mathbf V=V^k\mathbf e_k$ 를 미분하면 성분의 미분에 기저의 미분이 더해집니다. $$\begin{aligned}\partial_j\mathbf V&=\big(\partial_jV^i+\Gamma^i_{jk}V^k\big)\,\mathbf e_i\\&\equiv(\nabla_jV^i)\,\mathbf e_i.\end{aligned}$$
  2. 반변 성분: $\nabla_jV^i=\partial_jV^i+\Gamma^i_{jk}V^k$ (더하기). 공변 성분: $\nabla_jV_i=\partial_jV_i-\Gamma^k_{ji}V_k$ (빼기). 결과는 2계 텐서입니다. Arfken은 $V^i{}_{;j}$ 로 쓰는 것으로 기억합니다(기억 기반, 미검증).
  3. 계산 순서: 0이 아닌 $\Gamma$ 목록(7절)을 옆에 두고, 각 $(i,j)$ 마다 $\partial_jV^i$ 에 $\Gamma^i_{jk}V^k$ ($k$로 합)를 더합니다.
  4. 검산: 직교좌표에서 상수인 벡터장은 어느 좌표에서든 $\nabla_jV^i=0$ 입니다(성분은 변해도 화살표는 그대로이므로). 계량도 $\nabla_kg_{ij}=0$ 입니다.
  5. $i=j$ 로 축약하면 발산입니다: $\nabla_iV^i=\partial_iV^i+\Gamma^i_{ik}V^k$ (레시피 B).
예제 8-1 평면 극좌표에서 (가) 상수 벡터 $\hat{\mathbf y}$ (나) 회전 흐름 $(-y,\ x)$ 의 공변 미분

0이 아닌 $\Gamma$ 는 $\Gamma^r_{\theta\theta}=-r$, $\Gamma^\theta_{r\theta}=\Gamma^\theta_{\theta r}=\frac1r$ 뿐입니다(예제 7-1).

(가) 1절 방법으로 $V^r=\frac{\partial r}{\partial y}=\sin\theta$, $V^\theta=\frac{\partial\theta}{\partial y}=\frac{\cos\theta}{r}$.

$$\begin{aligned}\nabla_rV^r&=\partial_r(\sin\theta)=0,\\\nabla_\theta V^r&=\partial_\theta(\sin\theta)+\Gamma^r_{\theta\theta}V^\theta\\&=\cos\theta+(-r)\frac{\cos\theta}{r}=0,\\\nabla_rV^\theta&=\partial_r\Big(\frac{\cos\theta}{r}\Big)+\Gamma^\theta_{r\theta}V^\theta\\&=-\frac{\cos\theta}{r^2}+\frac1r\cdot\frac{\cos\theta}{r}=0,\\\nabla_\theta V^\theta&=\partial_\theta\Big(\frac{\cos\theta}{r}\Big)+\Gamma^\theta_{\theta r}V^r\\&=-\frac{\sin\theta}{r}+\frac1r\sin\theta=0.\end{aligned}$$

성분 $V^r$, $V^\theta$ 는 위치마다 변하지만, 그 변화는 기저벡터가 도는 것($\Gamma$ 항)과 정확히 상쇄됩니다. 상수 벡터라는 검산이 맞습니다.

(나) 예제 1-1에서 $V^r=0$, $V^\theta=1$.

$$\begin{aligned}\nabla_rV^r&=0,\\\nabla_\theta V^r&=0+\Gamma^r_{\theta\theta}V^\theta=-r,\\\nabla_rV^\theta&=0+\Gamma^\theta_{r\theta}V^\theta=\frac1r,\\\nabla_\theta V^\theta&=0+\Gamma^\theta_{\theta r}V^r=0.\end{aligned}$$

축약 $\nabla_iV^i=\nabla_rV^r+\nabla_\theta V^\theta=0$. 직교좌표의 $\nabla\cdot(-y,x)=0$ 과 같습니다. ($\partial_\theta V^\theta=0$ 이어도 $\nabla_\theta V^r=-r$ 처럼 $\Gamma$ 항에서 0이 아닌 성분이 생깁니다.)

레시피 B — 공변 발산
  1. 단위벡터 성분 $F_i$를 좌표(반변) 성분으로 바꿉니다: $V^i=F_i/h_i$. 이 단계를 빼먹는 실수가 가장 많습니다.
  2. $\sqrt g$를 구합니다 (직교 좌표계는 $h_1h_2h_3$).
  3. $\nabla\cdot\mathbf F=\dfrac1{\sqrt g}\sum_i\partial_i\big(\sqrt g\,V^i\big)$.
  4. 라플라시안은 $V^i=g^{ij}\partial_jf$ 를 넣은 것: $\nabla^2f=\dfrac1{\sqrt g}\partial_i\big(\sqrt g\,g^{ij}\partial_jf\big)$.
왜 되나? (증명 스케치)

축약 $\nabla_iV^i=\partial_iV^i+\Gamma^i_{ik}V^k$ 에서 $\Gamma^i_{ik}$ 만 알면 됩니다. 크리스토펠 공식을 $i$로 축약하면 $g^{il}\partial_ig_{kl}$ 과 $-g^{il}\partial_lg_{ik}$ 가 (더하는 글자 $i,l$ 을 맞바꾸면) 서로 지워지고 $\Gamma^i_{ik}=\tfrac12g^{il}\partial_kg_{il}$ 만 남습니다.

대각 계량이면 이것은 $\tfrac12\sum_i\dfrac{\partial_kg_{ii}}{g_{ii}}=\tfrac12\partial_k\ln(g_{11}g_{22}g_{33})=\partial_k\ln\sqrt g$ 입니다. (일반 계량에서도 $\partial_k\ln\det g=g^{il}\partial_kg_{il}$ 이라 결과가 같습니다.)

그러면 $\nabla_iV^i=\partial_iV^i+V^i\,\partial_i\ln\sqrt g=\dfrac1{\sqrt g}\partial_i\big(\sqrt g\,V^i\big)$ 입니다(곱의 미분을 거꾸로 읽음). 라플라시안은 $V^i=g^{ij}\partial_jf$ 를 넣은 것입니다.

예제 8-2 공변 발산으로 구면좌표의 발산 공식 다시 얻기

$\sqrt g=r^2\sin\theta$ 이고, 반변 성분은 $V^r=F_r$, $V^\theta=F_\theta/r$, $V^\varphi=F_\varphi/(r\sin\theta)$ 입니다.

$$\begin{aligned}\nabla\cdot\mathbf F&=\frac{1}{r^2\sin\theta}\Big[\partial_r\big(r^2\sin\theta\,F_r\big)\\&\qquad+\partial_\theta\Big(r^2\sin\theta\,\frac{F_\theta}{r}\Big)\\&\qquad+\partial_\varphi\Big(r^2\sin\theta\,\frac{F_\varphi}{r\sin\theta}\Big)\Big]\\&=\frac1{r^2}\partial_r(r^2F_r)+\frac{1}{r\sin\theta}\partial_\theta(\sin\theta\,F_\theta)\\&\quad+\frac1{r\sin\theta}\partial_\varphi F_\varphi.\end{aligned}$$

첫 항에서는 $\sin\theta$가, 둘째 항에서는 $r$이, 셋째 항에서는 $r\sin\theta$가 미분 밖으로 나왔습니다. L09 0절의 구면좌표 발산 공식과 똑같습니다.

빠른 확인: $\mathbf F=r\,\hat{\mathbf r}=(x,y,z)$ 이면 $F_r=r$ 이고 $\frac1{r^2}\partial_r(r^3)=3$. 직교좌표에서도 $\nabla\cdot(x,y,z)=3$ 입니다.

레시피 C — 측지선 방정식
  1. 곡선을 $u^k(s)$로 둡니다 ($s$는 호의 길이, 또는 그에 비례하는 매개변수).
  2. $\ddot u^k+\Gamma^k_{ij}\dot u^i\dot u^j=0$. $\Gamma$가 대칭이라 $i\ne j$ 항은 2배로 셉니다.
  3. 후보 곡선을 넣어 각 식이 0이 되는지 확인합니다.
예제 8-3 2-구에서 적도가 측지선임을 확인하기

예제 7-1(b)의 $\Gamma$ 두 종류를 넣습니다. $\Gamma^\varphi_{\theta\varphi}=\Gamma^\varphi_{\varphi\theta}$ 이라 두 번 셉니다.

$$\begin{aligned}\ddot\theta-\sin\theta\cos\theta\,\dot\varphi^2&=0,\\ \ddot\varphi+2\cot\theta\,\dot\theta\,\dot\varphi&=0.\end{aligned}$$

적도 $\theta=\tfrac\pi2$, $\varphi=\omega s$: $\dot\theta=\ddot\theta=0$, $\dot\varphi=\omega$, $\ddot\varphi=0$.

따라서 적도(대원)는 측지선입니다.

비교: 적도가 아닌 위도선 $\theta=\theta_0$ ($0<\theta_0<\tfrac\pi2$)이면 첫 식이 $-\sin\theta_0\cos\theta_0\,\omega^2\ne0$ 이라 측지선이 아닙니다.

9. 민코프스키 계량과 4-벡터

레시피
  1. $x^\mu=(x^0,x^1,x^2,x^3)=(ct,x,y,z)$, $\eta_{\mu\nu}=\operatorname{diag}(1,-1,-1,-1)$.
  2. 내리기: $a_\mu=\eta_{\mu\nu}a^\nu=(a^0,-a^1,-a^2,-a^3)$. 공간 성분만 부호가 바뀝니다.
  3. 내적: $a\cdot b=a^\mu b_\mu=a^0b^0-a^1b^1-a^2b^2-a^3b^3$.
  4. $x$방향 부스트 ($\beta=v/c$, $\gamma=1/\sqrt{1-\beta^2}$): $a'^\mu=\Lambda^\mu{}_\nu a^\nu$, $$\Lambda=\begin{pmatrix}\gamma&-\gamma\beta&0&0\\-\gamma\beta&\gamma&0&0\\0&0&1&0\\0&0&0&1\end{pmatrix}.$$
  5. 불변 조건은 $\Lambda^T\eta\Lambda=\eta$ 입니다. 이 식의 $(0,0)$ 성분이 $\gamma^2-\gamma^2\beta^2=1$, $(0,1)$ 성분이 $-\gamma^2\beta+\gamma^2\beta=0$ 입니다.
예제 9-1 $\beta=\tfrac35$ 부스트에서 $a=(2,1,0,0)$, $b=(3,0,1,0)$의 내적이 그대로인지 확인하기

$\gamma=\dfrac1{\sqrt{1-9/25}}=\dfrac54$, $\gamma\beta=\dfrac34$.

부스트 전: $a\cdot b=2\cdot3-1\cdot0-0\cdot1-0=6$, $a\cdot a=4-1=3$.

변환: $a'^0=\gamma a^0-\gamma\beta a^1$, $a'^1=-\gamma\beta a^0+\gamma a^1$, 나머지는 그대로.

$$\begin{aligned}a'&=\Big(\tfrac54\cdot2-\tfrac34\cdot1,\ -\tfrac34\cdot2+\tfrac54\cdot1,\ 0,\ 0\Big)\\&=\Big(\tfrac74,\ -\tfrac14,\ 0,\ 0\Big),\\ b'&=\Big(\tfrac54\cdot3,\ -\tfrac34\cdot3,\ 1,\ 0\Big)=\Big(\tfrac{15}4,\ -\tfrac94,\ 1,\ 0\Big).\end{aligned}$$

부스트 후:

$$\begin{aligned}a'\cdot b'&=\tfrac74\cdot\tfrac{15}4-\big(-\tfrac14\big)\big(-\tfrac94\big)-0\cdot1=\tfrac{105-9}{16}=6,\\ a'\cdot a'&=\tfrac{49}{16}-\tfrac1{16}=3.\end{aligned}$$

둘 다 그대로입니다. 지표를 내려서 계산해도 같습니다: $b_\mu=(3,0,-1,0)$ 이므로 $a^\mu b_\mu=2\cdot3+1\cdot0+0\cdot(-1)+0=6$.

10. 미분형식

레시피
  1. 쐐기곱(wedge product): $dx\wedge dy=-dy\wedge dx$, $dx\wedge dx=0$. 함수 계수는 앞으로 뺍니다: $(f\,dx)\wedge(g\,dy)=fg\,dx\wedge dy$. 이웃한 두 개를 맞바꿀 때마다 $-1$.
  2. 2차원: $(p\,dx+q\,dy)\wedge(s\,dx+t\,dy)=(pt-qs)\,dx\wedge dy$. 계수가 행렬식입니다.
  3. 외미분(exterior derivative) $d$: 0-형식은 $df=\partial_xf\,dx+\partial_yf\,dy+\partial_zf\,dz$. 그 위로는 $d(f\,dx^I)=df\wedge dx^I$ ($dx^I$는 $dx\wedge dy$ 같은 묶음).
  4. $d^2=0$: $d(df)$의 $dx^i\wedge dx^j$ 계수는 $\partial_i\partial_jf-\partial_j\partial_if=0$ 입니다. 편미분 순서는 바꿔도 같고 쐐기곱은 반대칭이기 때문입니다.
  5. 3차원 대응표:
    형식$d$를 걸면벡터 연산
    0-형식 $f$$\partial_xf\,dx+\partial_yf\,dy+\partial_zf\,dz$$\nabla f$
    1-형식 $A_x\,dx+A_y\,dy+A_z\,dz$$(\nabla\times\mathbf A)_x\,dy\wedge dz+(\nabla\times\mathbf A)_y\,dz\wedge dx+(\nabla\times\mathbf A)_z\,dx\wedge dy$$\nabla\times\mathbf A$
    2-형식 $F_x\,dy\wedge dz+F_y\,dz\wedge dx+F_z\,dx\wedge dy$$(\nabla\cdot\mathbf F)\,dx\wedge dy\wedge dz$$\nabla\cdot\mathbf F$
    $d^2=0$ 이 바로 $\nabla\times\nabla f=0$ 과 $\nabla\cdot(\nabla\times\mathbf A)=0$ 입니다.
  6. 적분에서 변수 바꾸기: $dx$, $dy$를 새 변수의 전미분으로 쓰고 쐐기곱을 전개합니다. 야코비안이 저절로 나옵니다. 순서를 바꾸면 부호가 바뀌는데, 이것이 방향(orientation)입니다.

그림 표기로 보기: G2 7절 — ∫Mdω = ∮∂Mω를 벡터 그림의 "풍선 뒤집기"로 보기

왜 되나? (증명 스케치)

쐐기곱 → 야코비안: 반대칭이므로 $(a\,dx+b\,dy)\wedge(c\,dx+d\,dy)=ad\,dx\wedge dy+bc\,dy\wedge dx=(ad-bc)\,dx\wedge dy$. "선형 + 두 개를 바꾸면 부호 반대"라는 성질은 행렬식을 정하는 성질과 같아서, 계수가 저절로 행렬식이 됩니다. $dx\wedge dx=0$ 은 "같은 행이 두 개면 행렬식이 0"과 같은 말입니다. 그래서 $dx\wedge dy$ 를 새 변수로 펼치면 야코비안이 나옵니다.

$d^2=0$: $d(df)=\sum_{i,j}\partial_j\partial_if\,dx^j\wedge dx^i$. $\partial_j\partial_if$ 는 $i,j$ 에 대칭, $dx^j\wedge dx^i$ 는 반대칭이라 $(i,j)$ 항과 $(j,i)$ 항이 서로 지웁니다(L07의 "대칭×반대칭=0"). 더 높은 형식도 계수마다 같은 논법입니다.

3차원 대응표로 옮기면 $d^2=0$ 은 $\nabla\times\nabla f=0$, $\nabla\cdot(\nabla\times\mathbf A)=0$ 이고, $\int_Md\omega=\oint_{\partial M}\omega$ 는 L08의 기울기 정리·스토크스 정리·발산정리를 한 식으로 쓴 것입니다(연습문제 G4).

예제 10-1 $\omega=x^2y\,dz$ 의 $d\omega$, 그리고 $d(d\omega)=0$

① $d\omega=d(x^2y)\wedge dz=(2xy\,dx+x^2\,dy)\wedge dz=2xy\,dx\wedge dz+x^2\,dy\wedge dz$.

순환 순서로 정리합니다: $dx\wedge dz=-dz\wedge dx$ 이므로

$$d\omega=x^2\,dy\wedge dz-2xy\,dz\wedge dx+0\,dx\wedge dy.$$

대응표 확인: $\mathbf A=(0,0,x^2y)$ 이면 $\nabla\times\mathbf A=(\partial_y(x^2y),\ -\partial_x(x^2y),\ 0)=(x^2,\ -2xy,\ 0)$. 같습니다.

② $d(d\omega)$: 각 계수의 $d$를 앞에 붙입니다.

$$\begin{aligned}d(d\omega)&=d(x^2)\wedge dy\wedge dz\\&\quad+d(-2xy)\wedge dz\wedge dx\\&=2x\,dx\wedge dy\wedge dz\\&\quad+(-2y\,dx-2x\,dy)\wedge dz\wedge dx.\end{aligned}$$

$dx\wedge dz\wedge dx=0$ ($dx$가 두 번). $dy\wedge dz\wedge dx$는 $dx\wedge dy\wedge dz$의 순환이라 같은 부호입니다. 따라서

$$\begin{aligned}d(d\omega)&=2x\,dx\wedge dy\wedge dz-2x\,dx\wedge dy\wedge dz\\&=0.\end{aligned}$$

벡터로는 $\nabla\cdot(x^2,-2xy,0)=2x-2x=0$ 입니다.

예제 10-2 $dx\wedge dy=r\,dr\wedge d\theta$ 와 원판 적분

$x=r\cos\theta$, $y=r\sin\theta$ 의 전미분:

$$\begin{aligned}dx&=\cos\theta\,dr-r\sin\theta\,d\theta,\\ dy&=\sin\theta\,dr+r\cos\theta\,d\theta.\end{aligned}$$

쐐기곱을 전개합니다. $dr\wedge dr=d\theta\wedge d\theta=0$ 이라 엇갈린 두 항만 남습니다.

$$\begin{aligned}dx\wedge dy&=\cos\theta\cdot r\cos\theta\,dr\wedge d\theta\\&\quad+(-r\sin\theta)\sin\theta\,d\theta\wedge dr\\&=r\cos^2\theta\,dr\wedge d\theta+r\sin^2\theta\,dr\wedge d\theta\\&=r\,dr\wedge d\theta.\end{aligned}$$

야코비안 $r$이 계산 중에 저절로 나왔습니다. 이것으로 반지름 $R$인 원판에서

$$\begin{aligned}\iint(x^2+y^2)\,dx\wedge dy&=\int_0^{2\pi}\!\!\int_0^Rr^2\cdot r\,dr\,d\theta\\&=2\pi\cdot\frac{R^4}{4}=\frac{\pi R^4}{2}.\end{aligned}$$

주의: $dx\wedge dy=-r\,d\theta\wedge dr$ 입니다. 적분 변수의 순서를 바꾸면 부호(방향)도 바뀝니다.

그림 표기로 보기: G7 3절 — 2차원 쐐기곱의 비가 곧 고정 변수 편미분입니다: $(\partial u/\partial x)_v=(du\wedge dv)/(dx\wedge dv)$ (L09 6절과 같은 식).

11. 자주 하는 실수

12. 연습문제

종이에 먼저 풀고 답을 여세요. 모든 답은 SymPy로 검산했습니다.

A. 변환 법칙

A1. 벡터 $\mathbf V=(x,\ y)$ (반지름 방향)의 극좌표 반변 성분과 공변 성분

답$V^r=r$, $V^\theta=0$, $V_r=r$, $V_\theta=0$ — $V^r=\cos\theta\cdot x+\sin\theta\cdot y=r$, $V^\theta=-\frac{\sin\theta}{r}x+\frac{\cos\theta}{r}y=0$.

A2. 상수 벡터 $\mathbf V=(1,\ 0)$의 극좌표 반변·공변 성분과 $V^iV_i$

답$V^r=\cos\theta$, $V^\theta=-\dfrac{\sin\theta}{r}$, $V_r=\cos\theta$, $V_\theta=-r\sin\theta$, $V^iV_i=\cos^2\theta+\sin^2\theta=1$.

A3. $f=x^2-y^2$의 기울기(공변 벡터) $(\partial_rf,\ \partial_\theta f)$를 변환 법칙 $V'_i=\frac{\partial x^j}{\partial x'^i}V_j$로 구하고, $f=r^2\cos2\theta$를 직접 미분한 것과 비교

답$(2r\cos2\theta,\ -2r^2\sin2\theta)$ — $\partial_rf=\cos\theta\cdot2x+\sin\theta\cdot(-2y)=2r(\cos^2\theta-\sin^2\theta)$, $\partial_\theta f=(-r\sin\theta)(2x)+(r\cos\theta)(-2y)=-4r^2\sin\theta\cos\theta$.

B. 회전, 축약, 분해

B1. $z$축 $90^\circ$ 회전 $R=\begin{pmatrix}0&1&0\\-1&0&0\\0&0&1\end{pmatrix}$ 로 $T=\begin{pmatrix}1&2&0\\3&4&0\\0&0&5\end{pmatrix}$ 를 변환하고, $\operatorname{tr}T$와 $T_{ij}T_{ij}$ (성분 제곱합)가 그대로인지 확인

답$T'=\begin{pmatrix}4&-3&0\\-2&1&0\\0&0&5\end{pmatrix}$, 대각합 $10$, 제곱합 $55$로 그대로 — $RT=\begin{pmatrix}3&4&0\\-1&-2&0\\0&0&5\end{pmatrix}$. 제곱합이 불변인 이유: $T'_{ij}T'_{ij}=R_{ik}R_{jl}R_{im}R_{jn}T_{kl}T_{mn}=T_{kl}T_{kl}$.

B2. $\mathbf a=(1,0,1)$, $\mathbf b=(0,2,1)$: $T_{ij}=a_ib_j$, 축약 $T_{ii}$, 대칭 부분 $S$, 반대칭 부분 $A$, 그리고 $A$의 쌍대 벡터를 $\tfrac12\mathbf a\times\mathbf b$와 비교

답$T=\begin{pmatrix}0&2&1\\0&0&0\\0&2&1\end{pmatrix}$, $T_{ii}=1=\mathbf a\cdot\mathbf b$, $S=\begin{pmatrix}0&1&\frac12\\1&0&1\\\frac12&1&1\end{pmatrix}$, $A=\begin{pmatrix}0&1&\frac12\\-1&0&-1\\-\frac12&1&0\end{pmatrix}$, 쌍대 $(A_{23},A_{31},A_{12})=(-1,-\frac12,1)$, $\mathbf a\times\mathbf b=(-2,-1,2)$ 이므로 정확히 절반.

C. 의사텐서와 쌍대

C1. 축벡터 $\boldsymbol\omega=(2,-1,3)$의 쌍대 반대칭 행렬 $\Omega_{ij}=\varepsilon_{ijk}\omega_k$ 를 쓰고, $\mathbf r=(1,0,0)$일 때 $\Omega\mathbf r$과 $\mathbf r\times\boldsymbol\omega$를 비교

답$\Omega=\begin{pmatrix}0&3&1\\-3&0&2\\-1&-2&0\end{pmatrix}$, $\Omega\mathbf r=(0,-3,-1)=\mathbf r\times\boldsymbol\omega$ — $\Omega_{13}=\varepsilon_{132}\omega_2=-(-1)=1$ 에서 부호에 주의.

C2. 반전 $\mathbf r\to-\mathbf r$에서 부호가 바뀌는 것은? (가) $\mathbf r$ (나) $\mathbf r\times\mathbf p$ (다) $\mathbf r\cdot(\mathbf p\times\mathbf q)$ ($\mathbf r,\mathbf p,\mathbf q$ 모두 극벡터) (라) 2계 텐서 $x_ix_j$

답(가) 바뀜, (나) 안 바뀜(축벡터), (다) 바뀜(의사스칼라, $(-1)^3$), (라) 안 바뀜($(-1)^2$).

C3. $xy$평면 거울 반사 $R=\operatorname{diag}(1,1,-1)$. $\mathbf a=(1,0,1)$, $\mathbf b=(0,1,1)$을 반사한 뒤의 $\mathbf a'\times\mathbf b'$을 직접 계산하고, 보통 벡터 법칙 $R(\mathbf a\times\mathbf b)$와 의사벡터 법칙 $(\det R)R(\mathbf a\times\mathbf b)$ 중 어느 것과 맞는지 보기

답$\mathbf a\times\mathbf b=(-1,-1,1)$, $\mathbf a'\times\mathbf b'=(1,0,-1)\times(0,1,-1)=(1,1,1)$. 보통 법칙은 $(-1,-1,-1)$로 틀리고, 의사벡터 법칙 $-(-1,-1,-1)=(1,1,1)$이 맞습니다.

D. 관성 텐서

D1. 질량 $m$인 점 두 개 $(a,0,0)$, $(0,a,a)$의 원점 관성 텐서와 주축

답$I=ma^2\begin{pmatrix}2&0&0\\0&2&-1\\0&-1&2\end{pmatrix}$, 주관성모멘트 $ma^2,\ 2ma^2,\ 3ma^2$, 주축 차례로 $\frac1{\sqrt2}(0,1,1)$, $(1,0,0)$, $\frac1{\sqrt2}(0,1,-1)$ — 2×2 블록 $\begin{pmatrix}2&-1\\-1&2\end{pmatrix}$의 고유값 $1,3$.

D2. 균일한 얇은 정사각형 판 (질량 $M$, $0\le x,y\le a$, $z=0$)의 꼭짓점(원점) 관성 텐서와 주축

답$I=Ma^2\begin{pmatrix}\frac13&-\frac14&0\\-\frac14&\frac13&0\\0&0&\frac23\end{pmatrix}$, 주관성모멘트 $\frac1{12}Ma^2$ (축 $\frac1{\sqrt2}(1,1,0)$), $\frac7{12}Ma^2$ (축 $\frac1{\sqrt2}(1,-1,0)$), $\frac23Ma^2$ (축 $\hat{\mathbf z}$). 검산: $\frac1{12}+\frac7{12}=\frac23$ (수직축 정리).

E. 계량, 크리스토펠 기호, 공변 미분, 측지선

E1. 원통좌표 $(\rho,\varphi,z)$의 $g_{ij}$, $g^{ij}$, $\sqrt g$, 그리고 0이 아닌 크리스토펠 기호

답$g_{ij}=\operatorname{diag}(1,\rho^2,1)$, $g^{ij}=\operatorname{diag}(1,\rho^{-2},1)$, $\sqrt g=\rho$. $\Gamma^\rho_{\varphi\varphi}=-\rho$, $\Gamma^\varphi_{\rho\varphi}=\Gamma^\varphi_{\varphi\rho}=\frac1\rho$ (평면 극좌표와 같고 $z$는 상관없음).

E2. 축 사이 각이 $60^\circ$인 비스듬한 좌표 $x=u+\tfrac12v$, $y=\tfrac{\sqrt3}{2}v$: $g_{ij}$, $g^{ij}$, 그리고 $V^i=(1,1)$의 $V_i$와 $|\mathbf V|^2$

답$g_{ij}=\begin{pmatrix}1&\frac12\\\frac12&1\end{pmatrix}$, $g^{ij}=\begin{pmatrix}\frac43&-\frac23\\-\frac23&\frac43\end{pmatrix}$, $V_i=(\frac32,\frac32)$, $|\mathbf V|^2=3$ — 직접 확인: $\mathbf V=\mathbf e_u+\mathbf e_v=(\frac32,\frac{\sqrt3}2)$, 길이 제곱 $\frac94+\frac34=3$. 축의 길이가 1이어도 축이 직교하지 않으면 $V_i\ne V^i$ 입니다.

E3. 3차원 구면좌표 $(r,\theta,\varphi)$의 0이 아닌 크리스토펠 기호 전부 (대각 계량 지름길 사용)

답$\Gamma^r_{\theta\theta}=-r$, $\Gamma^r_{\varphi\varphi}=-r\sin^2\theta$, $\Gamma^\theta_{r\theta}=\frac1r$, $\Gamma^\theta_{\varphi\varphi}=-\sin\theta\cos\theta$, $\Gamma^\varphi_{r\varphi}=\frac1r$, $\Gamma^\varphi_{\theta\varphi}=\cot\theta$, 그리고 아래 지표를 바꾼 짝들.

E4. 공변 발산으로 원통좌표의 발산 공식을 얻고, $\mathbf F=\rho^2\hat{\boldsymbol\rho}+z\hat{\mathbf z}$의 발산을 구하기

답$\nabla\cdot\mathbf F=\frac1\rho\partial_\rho(\rho F_\rho)+\frac1\rho\partial_\varphi F_\varphi+\partial_zF_z$, 값은 $3\rho+1$ — $V^\varphi=F_\varphi/\rho$, $\sqrt g=\rho$. 직교좌표 확인: $\mathbf F=(x\rho,\ y\rho,\ z)$, $\partial_x(x\rho)+\partial_y(y\rho)+1=2\rho+\frac{x^2+y^2}{\rho}+1=3\rho+1$.

E5. 공변 라플라시안으로 단위 구면($a=1$) 위의 라플라시안을 쓰고 $f=\cos\theta$에 걸기

답$\nabla^2f=\frac1{\sin\theta}\partial_\theta(\sin\theta\,\partial_\theta f)+\frac1{\sin^2\theta}\partial_\varphi^2f$, $\nabla^2\cos\theta=-2\cos\theta$ — $\frac1{\sin\theta}\partial_\theta(-\sin^2\theta)$. L22의 $-l(l+1)$, $l=1$과 같습니다.

E6. 평면 극좌표의 측지선 방정식을 쓰고, 직선 $x=1$ (즉 $r=\sqrt{1+t^2}$, $\theta=\arctan t$, $t$는 호의 길이)이 만족하는지 확인

답$\ddot r-r\dot\theta^2=0$, $\ddot\theta+\frac2r\dot r\dot\theta=0$. $\dot r=\frac tr$, $\ddot r=\frac1{r^3}$, $\dot\theta=\frac1{r^2}$, $\ddot\theta=-\frac{2t}{r^4}$ 를 넣으면 $\frac1{r^3}-\frac{r}{r^4}=0$, $-\frac{2t}{r^4}+\frac2r\cdot\frac tr\cdot\frac1{r^2}=0$. 만족합니다.

E7. 2-구(반지름 $a$)에서 (가) 경선 $\varphi=\varphi_0$, $\theta=s/a$ (나) 위도선 $\theta=\frac\pi6$, $\varphi=\omega s$ 가 측지선인가?

답(가) 예 — $\ddot\theta=0$, $\dot\varphi=\ddot\varphi=0$ 이라 두 식 모두 0. (나) 아니요 — 첫 식이 $-\sin\frac\pi6\cos\frac\pi6\,\omega^2=-\frac{\sqrt3}{4}\omega^2\ne0$.

E8. 평면 극좌표에서 $\mathbf V=(x,\ y)$ ($V^r=r$, $V^\theta=0$, 연습문제 A1)의 $\nabla_jV^i$ 네 성분과 축약 $\nabla_iV^i$

답$\nabla_rV^r=1$, $\nabla_\theta V^r=0+\Gamma^r_{\theta\theta}V^\theta=0$, $\nabla_rV^\theta=0+\Gamma^\theta_{r\theta}V^\theta=0$, $\nabla_\theta V^\theta=0+\Gamma^\theta_{\theta r}V^r=1$. 즉 $\nabla_jV^i=\delta^i_j$, 축약 $2$ = 직교좌표 $\nabla\cdot(x,y)=2$. ($\partial_\theta V^\theta=0$ 인데도 $\Gamma$ 항 때문에 1입니다.)

E9. 구면좌표에서 $\Gamma^i_{ir}$ 과 $\Gamma^i_{i\theta}$ ($i$로 합)를 E3의 목록으로 구하고 $\partial_k\ln\sqrt g$ 와 비교하세요.

답$\Gamma^i_{ir}=\Gamma^r_{rr}+\Gamma^\theta_{\theta r}+\Gamma^\varphi_{\varphi r}=0+\frac1r+\frac1r=\frac2r$, $\Gamma^i_{i\theta}=\Gamma^\varphi_{\varphi\theta}=\cot\theta$. $\ln\sqrt g=\ln(r^2\sin\theta)$ 를 $r$, $\theta$ 로 미분해도 $\frac2r$, $\cot\theta$. 같습니다.

E10. 평면 극좌표에서 $\nabla_rg_{\theta\theta}$ 와 $\nabla_\theta g_{r\theta}$ 를 계산해 0임을 확인하세요. (공변 2계 텐서: $\nabla_kg_{ij}=\partial_kg_{ij}-\Gamma^m_{ki}g_{mj}-\Gamma^m_{kj}g_{im}$)

답$\nabla_rg_{\theta\theta}=2r-\Gamma^\theta_{r\theta}g_{\theta\theta}-\Gamma^\theta_{r\theta}g_{\theta\theta}=2r-2\cdot\frac1r\cdot r^2=0$. $\nabla_\theta g_{r\theta}=0-\Gamma^\theta_{\theta r}g_{\theta\theta}-\Gamma^r_{\theta\theta}g_{rr}=-\frac1r\cdot r^2-(-r)\cdot1=0$.

F. 민코프스키

F1. $\beta=\tfrac45$ 부스트로 $a=(5,3,0,0)$을 변환하고 $a\cdot a$가 그대로인지 확인

답$\gamma=\frac53$, $a'=(\frac{13}3,-\frac53,0,0)$, $a\cdot a=25-9=16=\frac{169-25}{9}$.

F2. 빛의 4-벡터 $k=(1,1,0,0)$에 $\beta=\tfrac35$ 부스트를 걸고 $k\cdot k$를 확인

답$k'=(\frac12,\frac12,0,0)$, $k\cdot k=k'\cdot k'=0$ — 성분이 $\gamma(1-\beta)=\sqrt{\frac{1-\beta}{1+\beta}}=\frac12$ 배가 됩니다(도플러 인자).

G. 미분형식

G1. $(2\,dx+3\,dy)\wedge(dx-dy)$

답$-5\,dx\wedge dy$ — $-2\,dx\wedge dy+3\,dy\wedge dx$, 또는 행렬식 $2\cdot(-1)-3\cdot1$.

G2. $f=x^2yz$의 $df$를 구하고, $d(df)$의 세 계수가 모두 0임을 확인

답$df=2xyz\,dx+x^2z\,dy+x^2y\,dz$. $dx\wedge dy$ 계수 $\partial_x(x^2z)-\partial_y(2xyz)=0$, $dy\wedge dz$ 계수 $\partial_y(x^2y)-\partial_z(x^2z)=0$, $dz\wedge dx$ 계수 $\partial_z(2xyz)-\partial_x(x^2y)=0$.

G3. $\omega=yz\,dx+xz\,dy+xy\,dz$의 $d\omega$

답$0$ — $\omega=d(xyz)$ 이므로 $d\omega=d^2(xyz)=0$. 벡터로는 $\nabla\times\nabla(xyz)=0$.

G4. $\omega=-y\,dx+x\,dy$: $d\omega$를 구하고, 단위원판에서 $\iint d\omega$와 경계 원에서 $\oint\omega$를 각각 계산해 비교 (스토크스 정리 $\int_Dd\omega=\oint_{\partial D}\omega$)

답$d\omega=2\,dx\wedge dy$, $\iint d\omega=2\pi$. 경계 $x=\cos t$, $y=\sin t$에서 $\omega=(\sin^2t+\cos^2t)\,dt$ 이라 $\oint\omega=2\pi$. 같습니다.

G5. 구면좌표에서 $dx\wedge dy\wedge dz=r^2\sin\theta\,dr\wedge d\theta\wedge d\varphi$ 임을 쐐기곱 전개로 확인 (힌트: 계수는 $3\times3$ 행렬식)

답$dx=\sin\theta\cos\varphi\,dr+r\cos\theta\cos\varphi\,d\theta-r\sin\theta\sin\varphi\,d\varphi$ 등 세 1-형식을 쐐기곱하면 계수가 $\det\frac{\partial(x,y,z)}{\partial(r,\theta,\varphi)}=r^2\sin\theta$ 가 됩니다(L09의 야코비안과 같음).

13. 다음

← 로드맵 · ← 이전 L28 · 다음 L30 →