1. 개요
Équations de Lagrange / Lagrange 力 學라그랑주 역학은 라그랑주(J. L. Lagrange;1736~1813)가 1788년에 논문 해석 역학(Mécanique Analytique)에서 발표한 이론이며, 라그랑지언이라는 물리량을 통해서 물체의 운동을 설명하는 역학 체계이다.[1] "라그랑주 역학"이라는 거창한 이름만 보면 상대성 이론이라든가 양자역학처럼 뭔가 기존의 물리학으로는 설명할 수 없었던 새로운 패러다임을 만든 것 같지만, 사실은 그렇지 않다. 라그랑주 역학은 고등학교에서도 배우는 기존의 뉴턴 역학과 완전히 일치하는 결과를 준다. 다시 말하면, 라그랑주 역학으로 풀 수 있는 문제는 뉴턴 역학으로도 풀 수 있다. 에너지 장과 힘과의 관계를 통하여 뉴턴의 운동 방정식으로부터 오일러-라그랑주 정리를 유도할 수 있으며 반대로 오일러-라그랑주 정리로부터 뉴턴의 운동 방정식을 유도할 수 있다.
2. 뉴턴 역학과 라그랑주 역학의 비교
2.1. 계산의 편리함
라그랑주 역학은 계산이 편하다. 라그랑주 역학의 핵심 물리량은 라그랑지언으로, 스칼라이다. 기본적으로 공간을 기술하는 벡터는 성분이 3개이고, 스칼라는 하나이기 때문에 계산이 덜 복잡하다. 실제로 다양한 물리 문제에서 힘, 속도, 각운동량, 토크 같은 개념들이 등장하는데, 이들은 벡터이므로 방향을 고려해줘야 하지만, 라그랑지언은 스칼라이니 방향을 고려할 필요가 없다. 그래서 그림을 그려서 힘을 화살표로 표시하는 짓을 할 필요가 없고, 그냥 기계적으로 에너지만 구한 다음 방정식에 대입만 하면 다 풀리게 된다. 특히 라그랑주 역학에서는 일반화 좌표계라는 것을 도입하는데, 이렇게 좌표 변환을 해도 운동 방정식은 변하지 않는다.[2] 이것은 좌표 변환을 하면 복잡해지는 뉴턴 역학을 생각하면 엄청난 장점이라 할 수 있다.2.2. 철학적 관점의 차이
뉴턴 역학과 라그랑주 역학은 똑같은 자연현상을 기술하는 것처럼 보여도 둘은 물리 현상을 분석하는 시각이 근본적으로 다르다. 기본적으로, 뉴턴 역학에서 가장 중요하게 다루는 방정식은 누구나 들어본 적 있는 [math(\mathbf{F}=m \mathbf{a})]이다. 이 식은 물체에 힘 [math(\mathbf{F})]이 작용하면 가속도 [math(\mathbf{a})]를 만든다는 세계관을 표현하고 있다. 즉, 원인인 힘이 있으면 그 결과로 가속도가 생기는 인과론적인 세계관인 것이다. 이에 따르면 현재 상태를 완벽하게 알고 있으면 머나먼 과거부터 시작해 머나먼 미래까지 한치의 오차 없이 정확하게 예측할 수 있다.하지만 라그랑주 역학의 세계관은 그렇지 않다. 라그랑주 역학을 한 마디로 요약하면 다음과 같다.
자연은 항상
액션이 최소화되는 것을 좋아한다.
해밀턴의 원리
위 원리에서 나온 액션(action)이 무엇인지는 나중에 소개할 것이다. 이는 다시 말해 "신은 항상 자연을 최적화시킨다"[3]는 것이다. 이 놀라운 결과는 마치 자연이 항상 어떤 목적(액션의 최소화)을 가지고 굴러가는 것처럼 보인다는 것을 말해준다. 더 놀라운 것은 이 법칙 하나만 가지고 뉴턴 역학이랑 완전히 똑같은 체계를 만들어낼 수 있다는 것이다. 실제로 이것은 과학자들 사이에서 뉴턴 역학보다 더 근본적인 공리로 채택되고 있으며, 고전 역학을 넘어서 현대 물리학에까지 적용되는 아주 중요한 원리이다. 이는
최소작용의 원리(least action principle)라고도 부른다.해밀턴의 원리
라그랑주 역학에서는 잘 규정된 최초의 시간 · 위치 좌표와 최후의 시간 · 위치 좌표, 그리고 해밀턴의 원리를 잘 만족하는 라그랑지언, 이 두 가지를 필요로 한다.[4] 최소작용의 원리를 기술할 때 라그랑지언을 구성하는 속도와 위치를 편미분을 활용해 독립적으로 기술한다. 즉, 라그랑주 역학에서 특정 시간의 상황을 기술한다는 것은 단순히 속도가 얼마이고, 변위가 얼마인지는 기술하지만 이 둘이 어떻게 연결되어있는가는 보지 않는다는 것을 의미한다.[5] 바로 이 점이 뉴턴 역학과 라그랑주 역학을 나누는 관점의 차이다.
매 순간마다 속도와 위치가 어떤 식으로 연결되어 있는 지를 확인하지 않는다는 것, 과정을 기술하지 않는다는 것이다. 최초의 시간-위치 좌표와 최후의 시간-위치 좌표만 지나가기만 하면 되며, 동시에 해밀턴의 원리를 적용하면 좌표변환에 대해서 구조적으로 항상 불변한 어떤 방정식을 얻게 된다. 이 방정식이 다름이 아니라 오일러-라그랑주 방정식이다. 해밀턴의 원리가 적용된 라그랑지언이 나타내는 물리계의 운동현상을 궤적으로 나타내면, 최초의 좌표와 최후의 좌표 사이를 지나갈 때 항상 오일러-라그랑주 방정식을 만족하는 궤적을 따라 지나간다. 따라서 지정된 최초의 좌표와 최후의 좌표 내에서 해밀턴의 원리에 따라 일이 벌어진다는 것은 알 수 있어도 최초의 좌표와 최후의 좌표를 벗어난 머나먼 미래와 머나먼 과거에 대한 예측 및 추측에 대한 답을 라그랑지언이 직접적으로 주지 않는다. 다만 오일러-라그랑주 방정식을 활용하여 뉴턴 역학적 해석을 실시할 수 있을 뿐이다.
2.3. 현대물리학으로의 확장
바로 이러한 관점의 차이로 인하여 라그랑주 역학은 양자역학의 대수론적 해석에 동원될 수 있었다. 양자역학에서 우리가 유의미하게 정할 수 있는 것은 바로 최초의 양자상태[6]와 최후의 양자상태[7] 두 개 뿐이다. 그리고 어떤 상호작용을 거치는지를 나타내는 슈뢰딩거 방정식[8]을 통해 일이 벌어지는 과정(궤적)을 '이 정도 일어날 수 있다'란 확률로 답을 얻는다. 이로 인해 양자역학은 라그랑주 역학으로 해석할 수 있음을 알 수 있다.[9]이러한 이유 때문에, 라그랑주 역학은 현대물리학의 꽃이라고 할 수 있는 양자역학으로 넘어가는 데 아주 중요한 다리 역할을 하게 된다. 실제로, 미시 세계를 설명하는 양자역학에서는 힘 같은 물리량이나 원인과 결과 같은 개념이 고전적 의미를 상실한다. 대신 에너지를 의미하는 스칼라 함수 해밀토니언이 기본적인 물리량으로 쓰이는데, 이 해밀토니언이라는 것은 해밀턴 역학에서 나온 것이고 해밀턴 역학은 라그랑주 역학에서 나온 것이다. 결국 스칼라를 다루는 라그랑주 역학은 현대물리학의 기초가 될 수 있었던 것이다.
따라서 이러한 라그랑주 역학의 관점에서 물리학에 접근한다면 새로운 세계에 눈을 뜰 수 있을 것이다.
3. 라그랑지언
고전역학에서 어떠한 계(system)의 라그랑지언 [math(L)][10]은 그 계의 운동 에너지 [math(T)]에서 퍼텐셜 에너지 [math(U)]를 뺀 값으로 정의된다.[11] 즉, 다음과 같다.
[math(\displaystyle L \equiv T-U )]
라그랑지언을 표기할 때는 일반화 좌표계를 사용한다. 예를 들어 3차원에서 [math( n )]개의 점입자가 있는 상황을 생각해보자. 이 상황에서 입자 하나의 위치를 표기하려면 3개의 함수가 필요하므로, 모든 입자의 위치를 나타내려면 [math( 3n )]개의 시간에 대한 함수가 필요하다. 이를 [math( q_{i}(t) \,(i=1,\,2,\,3) )]로 표현하자.
이때, 해당 시스템의 라그랑지언은 아래와 같이 [math( 3n )]개의 변수와 그 미분값, 그리고 시간 [math( t )]의 함수로 나타낼 수 있다.
[math(\displaystyle L=L(q_{i}(t),\, \dot{q}_{i}(t) , \, t ) )]
그리고 라그랑지언의 시간 적분은 [math( 3n )]개의 함수를 입력으로 하는 범함수가 되고, 이것을 액션이라고 부른다.
[math(\displaystyle S \equiv \int_{t_{1}}^{t_{2}}L(q_{i}(t),\, \dot{q}_{i}(t) , \, t )\,{\rm d}t )]
보다 더 일반화된 액션(action)으로 아래와 같이 쓰기도 하고,
[math(\displaystyle s=\int L \,{\rm d}^nx )]
간혹
[math(\displaystyle s=\int L \,{\rm d}w_n )]
으로 표기되기도 한다.
4. 해밀턴의 원리
해밀턴의 원리(최소 작용의 원리, Hamilton's principle)는 자연이 선택하는 경로는 라그랑지안의 시간 적분인 액션 [math(S)]을 최소화시키는 경로라는 원리다. 수식으로 쓰면 다음과 같다.[math(\displaystyle \delta S=\delta \int_{t_{1}}^{t_{2}} L\,{\rm d}t=0 )]
이때 [math( \delta )]는 변분법에서의 표기법으로, [math( \delta S )]는 [math( S )]의 변분이라고 한다. 자세한 내용은 해밀턴의 원리 참고.
사실 역사적으로는 해밀턴의 원리(1834)보다 라그랑주 역학(1788)이 먼저 나온 법칙이다. 실제로 라그랑주는 처음에 해밀턴의 원리가 아닌 가상 일(virtual work)의 원리를 통해서 운동 방정식을 유도했다. 하지만 현재에는 해밀턴의 원리를 통해 라그랑주 역학을 설명하는 것이 더 명확하기 때문에 이 문서에서는 해밀턴의 원리를 먼저 소개한다.
5. 오일러-라그랑주 방정식
우선, 어떤 물리량 [math(a)]의 시간 미분을 다음과 같이 표기하기로 약속하자.
[math(\displaystyle \frac{{\rm d}a}{{\rm d}t} \equiv \dot{a} )]
즉, 시간 미분 연산자 [math({\rm d}/{\rm d}t)]를 점([math( \dot{\square})])으로 나타낸 것이다. 같은 논법으로, 시간에 대한 이계도함수는 다음과 같다.
[math(\displaystyle \frac{{\rm d}^2a}{{\rm d}t^2} \equiv \ddot{a} )]
물체의 위치에 대해 [math(x_{i})]축의 좌표를 [math(x_{i})][12]로 나타낼 것이며, 이것의 미분은 각 축에 대한 속도의 성분이 된다. 이때, 퍼텐셜 에너지는 위치에 대한 함수이므로 [math(U(x_{i}))]로 놓을 수 있고, 운동 에너지는 곧 이 위치의 시간 미분인 속도의 함수가 되므로 [math(T(\dot{x}_{i}))]가 된다. 따라서 라그랑지언은 이들의 선형 결합이므로 [math(x_{i})], [math(\dot{x}_{i})]의 함수이고, 라그랑지언이 속도와 위치의 함수라는 것을 명확히 해서 해밀턴의 원리를 다시 쓰면,
[math(\displaystyle \delta \int_{ t_1 }^{ t_2 }{ L(x_{i},\,\dot{x}_{i})\, {\rm d}t} = 0 )]
따라서 위의 해밀턴의 원리를 풀면 물체의 운동을 완벽히 설명할 수 있다. 고맙게도 위 방정식을 푸는 방법은 수학자들이 이미 다 연구해 놓았기 때문에 우리는 이걸 풀려고 노력할 필요가 없다. 놀랍게도 이 식은 다변수 함수의 오일러 방정식과 완전히 똑같은 꼴이다. 이 식의 해는 다음 방정식의 해와 같음이 알려져 있다.
[math(\displaystyle \frac{\partial L}{\partial x_{i}}-\frac{\rm d}{{\rm d}t}\frac{\partial L}{\partial \dot{x}_{i}}=0 )]
위 식이 바로, 입자의 오일러-라그랑주 방정식이며, 도출 과정은 해당 문서나, 변분법 문서를 참조하라.
이 방정식은 라그랑주 역학의 핵심이라고 할 수 있으며, 뉴턴 역학으로 따지면 [math( \mathbf{F} = m \mathbf{a} )]의 역할을 하는 아주 중요한 식이다. 만약 계가 [math( n )]개의 입자로 이루어져 있다면, 3차원에서 각 입자마다 좌표 3개씩 총 [math( 3n )]개의 좌표가 필요할 것이다. 입자 [math( n )]개 중 [math( \alpha )]번째 입자의 [math( i )]번째 좌표를 [math( x_{\alpha,\, i} )]라고 하자. 예를 들어 3번째 입자의 [math( y )]좌표는 [math(x_{3,\,2} )]이다. 그러면 계의 오일러-라그랑주 방정식은 다음과 같은 [math( 3n )]개의 연립방정식이 된다.
[math(\displaystyle \frac{\partial L}{\partial x_{\alpha i}}-\frac{\rm d}{{\rm d}t}\frac{\partial L}{\partial \dot{x}_{\alpha i}}=0 \quad (\alpha=1,\,2,\,3,\,\cdots,\,n ; \, i=1,\,2,\,3) )]
이러한 패러다임에 따라 역학 문제를 해석하는 방법을 라그랑주 역학이라고 부른다. 실제로 이 식은 뉴턴 역학과 같은 결과를 준다. 예를 들어 용수철에 매달려서 1차원 운동하는 입자를 생각해 보자. 이 입자는 [math(x)]만을 이용하여, 위치를 기술할 수 있으므로, 입자의 퍼텐셜 에너지와 운동 에너지는 다음과 같다.
[math(\displaystyle \begin{aligned} U&=\frac{1}{2}kx^{2} \\ T&=\frac{1}{2}m \dot{x}^{2} \end{aligned} )]
따라서 이 계의 라그랑지언은
[math(\displaystyle \begin{aligned} L&=T-U \\ &=\frac{1}{2}m \dot{x}^{2}-\frac{1}{2}kx^{2} \end{aligned} )]
이고, 이 라그랑지언을 오일러-라그랑주 방정식
[math(\displaystyle \frac{\partial L}{\partial x}-\frac{\rm d}{{\rm d}t}\frac{\partial L}{\partial \dot{x}}=0 )]
에 대입하면, 운동 방정식은 다음과 같다.
[math(\displaystyle m \ddot{x} =-kx )]
이것은 뉴턴 역학의 결과와 완전히 일치한다.
오일러-라그랑주 방정식을 풀 때 유념할 점이 하나 있다. 기존 편미분에서는 [math( x )]를 [math( \dot{x} )]로 편미분할 때의 값은 자명하게 0이라고 할 수 없다. 그러나 특별히 라그랑지안을 적용하는 오일러-라그랑주 방정식에서는 이것이 0이다. 즉, 시간에 대한 도함수 꼴로 관련이 있는 변수여도 모양이 다르면 아예 다른 변수로 가정하여 문제를 해결한다.
5.1. 일반화 좌표계에서
오일러-라그랑주 방정식을 풀 때의 핵심은 일반화 좌표계(generalized coordinates)를 도입하는 것이다. 기본적으로 라그랑주 역학은 전부 일반화 좌표계에서 풀게 된다.자유도가 [math( s )]인 일반화 좌표계에서 운동 방정식을 풀려면 그냥 다음 절차에 따르기만 하면 쉽게 풀 수 있다.
|
여기서 2~5번 과정은 그냥 컴퓨터에 넣기만 하면 바로 계산이 된다. 따라서 우리가 고민해야 되는건 1번밖에 없는데, 이것조차 뉴턴 역학으로 푸는 것에 비하면 쉽다. 이런 방법으로 풀어도 되는 이유를 알아보자. 오일러-라그랑주 방정식을 처음 소개할 때, 자유도가 [math( s )]인 계의 라그랑지언은 다음과 같이 물체들의 위치와 속도에 대한 함수로 나타낼 수 있었다.
[math(\displaystyle \begin{aligned} L&=T(\dot{x}_{\alpha,\, i}) - U(x_{\alpha,\, i}) \\ &=L(x_{\alpha,\, i},\,\dot{x}_{\alpha,\,i}) \end{aligned} )]
여기서 [math( x_{\alpha,\,i} , \dot{x}_{\alpha,\,i} )]는 [math( s )]개의 일반화 좌표계와 시간에 대한 함수로 나타낼 수 있다.
[math(\displaystyle \begin{aligned} x_{\alpha,\,i}&=x_{\alpha,\,i}(q_{j},\,t) \\ x_{\alpha,\,i}&=\dot{x}_{\alpha,\,i}(q_{j},\, \dot{q}_{j} ,\, t) \quad (j=1,\,2 ,\,3 ,\, \cdots ,\,s) \end{aligned} )]
이를 대입하면 결국 라그랑지언도 일반화 좌표계로 위치, 속도, 시간에 대한 함수가 된다.
[math(\displaystyle L=L(q_{j},\, \dot{q}_{j} ,\, t) \quad (j=1,\,2 ,\,3 ,\, \cdots ,\,s) )]
그러면 해밀턴의 원리도 다음과 같이 일반화 좌표계로 쓸 수 있다.
[math(\displaystyle \delta \int L(q_{j},\, \dot{q}_{j} ,\, t)\,{\rm d}t=0 \quad (j=1,\,2 ,\,3 ,\, \cdots ,\,s) )]
그런데 이것은 좌표 변환을 하기 전과 완전히 똑같은 꼴이다. 즉, 이것도 다변수 함수의 오일러 방정식으로 풀 수 있다. 그러면 다음과 같이 [math( s )]개의 연립 미분방정식을 얻는다.
[math(\displaystyle \frac{\partial L}{\partial q_{j}}-\frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{j}} \quad (j=1,\,2 ,\,3 ,\, \cdots ,\,s) )]
위 식은 그냥 오일러-라그랑주 방정식에 좌표 알파벳만 바꿔 넣은 꼴이다.
5.2. 구속력을 고려하는 경우
지금까지 다뤘던 라그랑주 역학의 특징은 어려운 구속력을 그냥 무시해 버리고 풀 수 있다는 것이다. 하지만 가끔은 구속력을 무시하면 안 될 때가 있다. 예를 들어 엘리베이터를 설계할 때에는 엘리베이터를 매다는 줄이 어느 정도 장력까지 버틸 수 있는지를 고려해야 하고, 대교를 짓는다면 얼마나 많은 차가 다녀도 하중을 견딜 수 있는지를 계산해야 한다. 그런데 라그랑주 역학이 또 사기적인 것은 구속력도 뉴턴 역학보다 훨씬 쉽게 구할 수 있다는 것이다.지금까지는 자유도가 [math( s )]인 계는 [math( s )]개의 연립방정식을 풀어서 구했다. 예를 들어 평면 위에서 원운동하는 입자는 자유도가 1이므로 하나의 방정식만 있으면 된다. 하지만 이렇게 해서는 구속력을 구할 수 없다. 따라서 다른 방법을 시도해 봐야 한다. 이는 계를 자유도 2라고 간주하고, 구속 조건을 고려하는 것이다. 즉, "자유도 2개에 의한 오일러 방정식 2개 + 구속조건 식 1개" 이렇게 해서 총 3개의 방정식이 나온다. 그러면 식이 3개니까 미지수도 3개를 풀 수 있다. 구할 수 있는 것은 "일반화 좌표 2개 + 구속력"이다. 이렇게 하면 구속력까지 모두 구할 수 있게 된다.
보통 구속 조건은 일반화 좌표로 [math( f(q_j, \, t) = 0 )]의 꼴로 쓸 수 있다.[13] 그런데 이걸 어떻게 풀 수 있을까? 답은 구속 조건이 있는 오일러 방정식에 그대로 나와 있다. 완전히 똑같은 꼴이기 때문이다. 따라서 구속 조건을 [math( f(q_j, \,t)=0 )]이라고 하면 어떤 함수 [math( \lambda (t) )]에 대하여 다음 식이 성립한다.
[math(\displaystyle {\partial L \over \partial q_j} - {{\rm d} \over {\rm d}t} {\partial L \over \partial \dot{q}_j} + \lambda (t) {\partial f \over \partial q_j } = 0 )]
구속 조건이 여러 개일 경우 [math( \lambda )]를 여러 개 쓰면 된다.
[math(\displaystyle {\partial L \over \partial q_j} - {{\rm d} \over {\rm d}t} {\partial L \over \partial \dot{q}_j} + \lambda_1 {\partial f \over \partial q_j } + \lambda_2 {\partial g \over \partial q_j } = 0 )]
마찬가지로 구속 조건이 더 많으면 [math( \lambda )] 항을 계속 늘리면 된다.
이는 오일러-라그랑주 방정식에 그냥 [math( \lambda ({\partial f }/{ \partial q_j }) )]같은 항들 여러 개만 더한 꼴이다. [math( m )]개의 구속 조건이 있을 때 이 항들을 보통 [math( Q_j )]로 써서
[math(\displaystyle Q_j = \sum_{k=1}^{m} { \lambda_k {\partial f_k \over \partial q_j }} )]
로 정의하고, [math( Q_j )]를 일반화 구속력 또는 일반화 힘이라고 한다. 이것은 구속력의 그 좌표 방향 성분이며, 단위는 힘이 아닐 수도 있다.
신기하게도, 일반화 구속력은 그 좌표가 어떤 단위인지에 따라 실제 구속력이 나온다. 예를 들어, 극좌표계에서
[math(\displaystyle Q_{r} = \lambda {\partial f \over \partial r } )]
이면 구속력 [math( Q_r )]는 [math( r )] 방향의 구속력이 된다. 그런데,
[math(\displaystyle Q_{\theta} = \lambda {\partial f \over \partial \theta } )]
이면 구속력 [math( Q_{\theta} )]는 [math( \theta )] 방향의 구속 토크가 된다.
즉, 구속력 항을 더한 오일러-라그랑주 방정식들과 구속 조건까지 다 연립해서 풀면 모든 일반화 좌표랑 구속력까지 빠짐없이 구할 수 있다.
이때, 구속력의 부호는 일반화 좌표의 기저 벡터의 방향을 기준으로 정한다. 즉, 양의 값을 가지면 그 기저 벡터와 같은 방향, 음의 값을 가지면 반대 방향인 것이다.
5.3. 일반적인 오일러-라그랑주 방정식
밑에서 뉴턴 역학과의 등가성을 확인하는 경우에서 다음을 다뤘다.[math(\displaystyle \frac{{\rm d}}{{\rm d}t}\frac{\partial T}{\partial \dot{q}_{j}}-\frac{\partial T}{\partial {q}_{j}}=Q_{j} )]
여기서 [math(Q_{j})]를 다음과 같이 분리하는데,
[math(\displaystyle Q_{j}=C_{j}+N_{j} )]
[math(C_{j})]는 퍼텐셜의 음의 그레이디언트를 통해 얻을 수 있는 힘, 즉 보존력인 경우이다. [math(N_{j})]는 그러한 방식으로 얻을 수 없는 비보존력이다. 밑에서 다룬 것을 이용하면 다음과 같은 일반적인 오일러-라그랑주 방정식을 얻는다.
[math(\displaystyle \frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{j}}-\frac{\partial L}{\partial {q}_{j}}=N_{j} )]
다만 [math(N_{j})]가 속도에 의존하는 경우 등은 매우 풀기 어려워진다. 특히 속도에 비례하는 마찰력같은 것이 작용할 경우 라그랑주 역학으로는 분석하기 까다로워진다.
6. 일반화 운동량
일반화 운동량은 아래와 같이 정의되는 물리량이다.[math( \displaystyle p_{j}=\frac{\partial L}{\partial \dot{q}_{j}} )]
이 운동량은 뉴턴 역학에서 말하는 선운동량과 각운동량 모두 포함되어 있으며, 서로 대응되는 관계에 있다.
오일러-라그랑주 방정식에 의하면
[math( \displaystyle \frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{j}}=\frac{\partial L}{\partial {q}_{j}} )]
일반화 운동량의 정의에 의하여
[math( \displaystyle \dot{p}_{j}=\frac{\partial L}{\partial {q}_{j}} )]
따라서 다음의 두 관계식을 얻는다. 이 두 관계는 해밀토니언을 얻을 때 다시 등장하므로 잘 살펴볼 필요가 있다.
[math( \displaystyle \begin{aligned} p_{j}&=\frac{\partial L}{\partial \dot{q}_{j}} \\ \dot{p}_{j}&=\frac{\partial L}{\partial q_{j}} \end{aligned})]
7. 운동 에너지 정리
뉴턴 역학에서 운동 에너지는[math(\begin{aligned} T=\sum_{\alpha i}\frac{1}{2}m_{\alpha} {\dot x_{\alpha i}}^2 \end{aligned})]
으로 나타낸다.
한편, 일반화 좌표
[math(\begin{aligned} x_{\alpha i}=x_{\alpha i}(q_{j},\,t) \end{aligned})]
를 도입한다면
|
[math(\begin{aligned} \dot{x}_{\alpha i} &=\sum_{j} \frac{\partial x_{\alpha i}}{\partial q_{j}}\dot{q}_{j}+\frac{\partial x_{\alpha i}}{\partial t} \\ {\dot x_{\alpha i}}^2&=\sum_{jk} \frac{\partial x_{\alpha i}}{\partial q_{j}}\frac{\partial x_{\alpha i}}{\partial q_{k}}\dot{q}_{j}\dot{q}_{k}+2\sum_{j} \frac{\partial x_{\alpha i}}{\partial q_{j}} \frac{\partial x_{\alpha i}}{\partial t} \dot{q}_{j}+\left(\frac{\partial x_{\alpha i}}{\partial t} \right)^{2} \end{aligned})] |
|
[math(\begin{aligned} T=\sum_{\alpha} \left[\sum_{ijk} \frac{1}{2} m_{\alpha}\frac{\partial x_{\alpha i}}{\partial q_{j}}\frac{\partial x_{\alpha i}}{\partial q_{k}}\dot{q}_{j}\dot{q}_{k}+\sum_{ij} m_{\alpha} \frac{\partial x_{\alpha i}}{\partial q_{j}} \frac{\partial x_{\alpha i}}{\partial t} \dot{q}_{j}+\frac{1}{2}m_{\alpha}\left(\frac{\partial x_{\alpha i}}{\partial t} \right)^{2} \right] \end{aligned})] |
[math(\begin{aligned} T=\sum_{jk} a_{jk} \dot{q}_{j}\dot{q}_{k}+\sum_{j} b_{j} \dot{q}_{j} + c \end{aligned})]
여기서 좌표와 일반화 좌표 사이의 변환 식이 시간에 의존하지 않을 때, 즉 스클로노믹(scleronomic)일 때
[math(\begin{aligned} \frac{\partial x_{\alpha i}}{\partial t}=0 \end{aligned})]
이므로 [math(b_{j}=c=0)]이 된다.
따라서 일반화 속도로 나타내어진 운동 에너지는
[math(\begin{aligned} T=\sum_{jk} a_{jk} \dot{q}_{j}\dot{q}_{k} \end{aligned})]
이 된다.
다음을 조사하자.
|
[math(\begin{aligned} \frac{\partial T}{\partial \dot{q}_{l}}&=\sum_{jk} a_{jk} \frac{\partial}{\partial \dot{q}_{l}}(\dot{q}_{j} \dot{q}_{k}) \\&=\sum_{jk} a_{jk} \left[ \frac{\partial \dot{q}_{j}}{\partial \dot{q}_{l}} \dot{q}_{k}+\dot q_j \frac{\partial \dot{q}_{k}}{\partial \dot{q}_{l}} \right] \\&=\sum_{jk} a_{jk} \left[ \delta_{jl} \dot{q}_{k}+\dot q_j \delta_{kl} \right] \\ &=\sum_{k} a_{lk} \dot{q}_{k}+\sum_{j} a_{jl }\dot{q}_{j} \end{aligned})] |
[math(\begin{aligned} \frac{\partial T}{\partial \dot{q}_{l}}&=2\sum_{k} a_{lk} \dot{q}_{k} \end{aligned})]
다음을 고려하자.
[math(\begin{aligned} \sum_{l}\dot{q}_{l}\frac{\partial T}{\partial \dot{q}_{l}}&=2\sum_{lk} a_{lk} \dot{q}_{l}\dot{q}_{k} \end{aligned})]
그런데 우변은 [math(2T)]와 같으므로 다음을 얻는다.
[math(\begin{aligned} \therefore \sum_{l}\dot{q}_{l}\frac{\partial T}{\partial \dot{q}_{l}}&=2T \end{aligned})]
8. 뉴턴 역학과의 등가성
라그랑주 역학이 뉴턴 역학과 등가성이 있는지 조사해보자.8.1. 뉴턴 역학에서 라그랑주 역학으로 환원하기
먼저해야 할 일은 달랑베르의 원리를 상기하는 것이다.[math(\displaystyle \delta W=(\mathbf{F}-\mathbf{\dot{p}})\boldsymbol{\cdot}\delta \mathbf{r}=0 )]
이 원리는 평형을 이루는 물체의 가상 변위에 대한 가상 일은 없음을 나타낸다.
이 원리를 사용하여 우리의 목적을 이루어보자. 우선 위 것을 성분 별로 쓴다.
[math(\displaystyle \sum_{i} (F_{i}-\dot{p}_{i})\delta x_{i}=0 )]
뉴턴 역학에서 [math(p_{i}=m\dot{x}_{i})]이기에 [math(\dot{p}_{i}=m\ddot{x}_{i})]이다.
라그랑주 역학에서는 일반화 좌표를 도입하기 때문에 [math(x_{i}(q_{j},\,t))]로 쓸 수 있다. 따라서
[math(\displaystyle \delta x_{i}=\sum_{j} \frac{\partial x_{i}}{\partial q_{j}} \delta q_{j} )]
이다. 이는 가상 변위를 고려하기에 시간에 대한 고려를 할 필요가 없기 때문에 시간에 대한 미분 항은 나타나지 않은 것이다. 이것은 변분법 문서의 [math(\delta)]기호에 관한 부분을 보는 것을 권장한다.
이상에서 위의 달랑베르 원리에서 운동량 항만 보게 되면,
[math(\displaystyle \sum_{i} \dot{p}_{i}\delta x_{i}=\sum_{ij} m\ddot{x}_{i} \frac{\partial x_{i}}{\partial q_{j}} \delta q_{j} )]
그런데
[math(\displaystyle \frac{{\rm d}}{{\rm d}t}\left[ \dot{x}\frac{\partial x_{i}}{\partial q_{j}} \right]=\ddot{x}\frac{\partial x_{i}}{\partial q_{j}}+\dot{x}\frac{{\rm d}}{{\rm d}t}\frac{\partial x_{i}}{\partial q_{j}} )]
인데, 여기서 우변의 제 2항의 미분 연산은 편미분 안으로 넣을 수 있어서,
|
[math(\displaystyle \frac{{\rm d}}{{\rm d}t}\left[ \dot{x}\frac{\partial x_{i}}{\partial q_{j}} \right]=\ddot{x}\frac{\partial x_{i}}{\partial q_{j}}+\dot{x}\frac{\partial \dot{x}_{i}}{\partial q_{j}} )]
따라서
[math(\displaystyle \sum_{i} \dot{p}_{i}\delta x_{i}=\sum_{ij} m \left[\frac{{\rm d}}{{\rm d}t}\left[ \dot{x}\frac{\partial x_{i}}{\partial q_{j}} \right]-\dot{x}\frac{\partial \dot{x}_{i}}{\partial q_{j}} \right] \delta q_{j} )]
여기서 다음의 사실을 이용한다.
[math(\displaystyle \frac{\partial x_{i}}{\partial q_{j}}=\frac{\partial \dot{x}_{i}}{\partial \dot{q}_{j}} )]
|
이에
|
[math(\displaystyle \begin{aligned} \sum_{i} \dot{p}_{i}\delta x_{i}&=\sum_{ij} m \left[\frac{{\rm d}}{{\rm d}t}\left[ \dot{x}\frac{\partial \dot{x}_{i}}{\partial \dot{q}_{j}} \right]-\dot{x}\frac{\partial \dot{x}_{i}}{\partial q_{j}} \right] \delta q_{j} \\ &=\sum_{ij} \left[\frac{{\rm d}}{{\rm d}t} \frac{\partial}{\partial \dot{q}_{j}}\left[ \frac{1}{2}m {\dot x_i}^2 \right]-\frac{\partial}{\partial {q}_{j}}\left[ \frac{1}{2}m {\dot x_i}^2 \right] \right] \delta q_{j} \\&=\sum_{j} \left[\frac{{\rm d}}{{\rm d}t} \frac{\partial}{\partial \dot{q}_{j}}\left[ \frac{1}{2}m \left(\sum_{i} {\dot x_i}^2 \right) \right]-\frac{\partial}{\partial {q}_{j}}\left[ \frac{1}{2}m \left(\sum_{i} {\dot x_i}^2 \right) \right] \right] \delta q_{j} \\&=\sum_{j} \left[\frac{{\rm d}}{{\rm d}t} \frac{\partial}{\partial \dot{q}_{j}}\left[ \frac{1}{2}m |\mathbf{\dot{r}}|^{2}\right]-\frac{\partial}{\partial {q}_{j}}\left[ \frac{1}{2}m |\mathbf{\dot{r}}|^{2} \right] \right] \delta q_{j} \end{aligned} )] |
[math(\displaystyle \begin{aligned} \frac{\partial}{\partial u}(\dot{x}_{i}^{2})=2\dot{x}_{i} \frac{\partial \dot{x}_{i}}{\partial u} \end{aligned} )]
이상에서
[math(\displaystyle \begin{aligned} \sum_{i} \dot{p}_{i}\delta x_{i}&=\sum_{j} \left[\frac{{\rm d}}{{\rm d} t} \frac{\partial T}{\partial \dot{q}_{j}}-\frac{\partial T}{\partial {q}_{j}} \right] \end{aligned} )]
여기서 나온 [math(T)]가 운동 에너지이다.
이번엔 힘과 관련된 항을 보자.
[math(\displaystyle \sum_{i} F_{i}\delta x_{i}= \sum_{ij} F_{i}\frac{\partial x_{i}}{\partial q_{j} }\delta q_{j} )]
여기서 나온
[math(\displaystyle \sum_{i}F_{i}\frac{\partial x_{i}}{\partial q_{j} } \equiv Q_{j} )]
를 일반화 힘(generalized force)라 정의한다.
이상에서
[math(\displaystyle \sum_{j}\left[Q_{j}-\left[\frac{{\rm d}}{{\rm d} t} \frac{\partial T}{\partial \dot{q}_{j}}-\frac{\partial T}{\partial {q}_{j}} \right] \right]\delta q_{j}=0 )]
인데, [math(\delta q_{j} \neq 0)]이므로 일반적으로 위 등식이 성립하려면
[math(\displaystyle \frac{{\rm d}}{{\rm d} t} \frac{\partial T}{\partial \dot{q}_{j}}-\frac{\partial T}{\partial {q}_{j}}=Q_{j} \qquad \cdots \, \small{(\ast)})]
이것이 라그랑주 역학 버전의 달랑베르의 원리이다.
가장 간단하게 힘을 퍼텐셜 에너지의 음의 그레이디언트로 얻을 수 있다 가정하면
[math(\displaystyle F_{i}=-\frac{\partial U}{\partial x_{i}})]
이 때 주어지는 일반화 힘은 다음과 같다.
[math(\displaystyle -\frac{\partial U}{\partial q_{j} } = Q_{j} )]
이것을 [math((\ast))]에 대입하면
[math(\displaystyle \frac{{\rm d}}{{\rm d} t} \frac{\partial T}{\partial \dot{q}_{j}}-\frac{\partial (T-U)}{\partial {q}_{j}}=0 )]
퍼텐셜 에너지가 일반화 속도에 무관하다면 이 식을 다음과 같이 써도 무방하다.
[math(\displaystyle \frac{{\rm d}}{{\rm d} t} \frac{\partial (T-U)}{\partial \dot{q}_{j}}-\frac{\partial (T-U)}{\partial {q}_{j}}=0 )]
여기서 나온 [math(L \equiv T-U)]가 바로 라그랑지언이고,
[math(\displaystyle \frac{\partial L}{\partial {q}_{j}}-\frac{{\rm d}}{{\rm d} t} \frac{\partial L}{\partial \dot{q}_{j}}=0 )]
위 식은 곧 오일러-라그랑주 방정식이다.
위 논의는 라그랑주 역학이 뉴턴 역학과 전혀 다른 새로운 이론이 아닌 서로 등가성을 이룸을 보여준다.
8.1.1. 퍼텐셜이 일반화 속도에 의존하는 경우
일반화 속도에 의존하는 퍼텐셜 [math(\tilde{U}(q_{j},\,\dot{q}_{j}))]을 고려하는데, 여기에 특수한 케이스로써[math(\displaystyle Q_{j}=-\frac{\partial \tilde{U}}{\partial q_{j}}+\frac{{\rm d}}{{\rm d}t}\frac{\partial \tilde{U}}{\partial \dot{q}_{j}} )]
을 만족한다면 오일러-라그랑주 방정식은
[math(\displaystyle \frac{\partial (T-\tilde{U})}{\partial q_{i}}-\frac{{\rm d}}{{\rm d}t}\frac{\partial (T-\tilde{U})}{\partial \dot{q}_{i}}=0)]
형식이 되고, 라그랑지언을 마찬가지로 [math(T-\tilde{U} \equiv L)]로 정의하면 된다.
대표적인 예로 로런츠 힘 하에 있는 입자의 라그랑지언이 있다.
하지만 모든 퍼텐셜이 해당 조건을 만족하는 것이 아님에 유의해야 한다.
8.2. 라그랑주 역학에서 뉴턴 역학으로 환원하기
윗 문단에서 뉴턴 역학으로 라그랑주 역학을 유도해봄으로써 두 이론 사이의 등가성을 확인했다.이번엔 라그랑주 역학에서 뉴턴 역학으로 환원을 해보자. 일반화 좌표로 직교 좌표를 쓰면,
[math(\displaystyle T=\frac{1}{2}m |\mathbf{\dot{r}}|^{2} )]
이다. 오일러-라그랑주 방정식
[math(\displaystyle \frac{\partial (T-U)}{\partial x_{i}}-\frac{{\rm d}}{{\rm d}t}\frac{\partial (T-U)}{\partial \dot{x}_{i}}=0 )]
인데, 일반적인 상황에서 퍼텐셜 에너지는 속도에 의존하지 않고, 운동 에너지는 위에서 본 것 처럼 위치에 의존하지 않는다. 따라서 오일러 라그랑주 방정식을
[math(\displaystyle -\frac{\partial U}{\partial x_{i}}-\frac{{\rm d}}{{\rm d}t}\frac{\partial T}{\partial \dot{x}_{i}}=0 )]
형식으로 쓸 수 있고,
[math(\displaystyle -\frac{\partial U}{\partial x_{i}}=\frac{{\rm d}}{{\rm d}t}(m v_{i}) )]
이것은
[math(\displaystyle -\boldsymbol{\nabla}U=\mathbf{\dot{p}} )]
의 한 성분이고, 힘과 퍼텐셜 사이의 관계를 생각해볼 때
[math(\displaystyle \mathbf{F}=\mathbf{\dot{p}} )]
로 뉴턴의 운동 법칙으로 환원됨을 알 수 있다.
9. 보존 정리
아래의 내용을 요약하면 아래와 같다.| 관성계의 성질 | 라그랑지언의 성질 | 보존되는 양 |
| 시간의 균질성 |
시간 이동에 대한 불변 (시간의 직접적인 함수가 아님) |
총 에너지 |
| 공간의 균질성 | 병진에 대한 불변 | 선운동량 |
| 공간의 등방성 | 회전에 대한 불변 | 각운동량 |
아래는 학부 수준의 증명이다. 더욱 더 엄밀한 증명은 뇌터 정리 문서를 참조한다.
9.1. 에너지 보존
라그랑지언의 시간 이동에 대한 불변은 에너지 보존을 이끌어낸다. [math(\delta t)] 동안의 이동에 불변이라고 하면 다음이 성립한다.[math(\displaystyle \begin{aligned} L(q_{i},\,\dot{q}_{i},\,t)&=L(q_{i},\,\dot{q}_{i},\,t+\delta t) \\&=L(q_{i},\,\dot{q}_{i},\,t)+\frac{\partial L}{\partial t}\delta t \end{aligned})]
따라서 다음을 얻는다.
[math(\displaystyle \begin{aligned} \frac{\partial L}{\partial t}\delta t=0 \end{aligned})]
한편, [math(\delta t \neq 0)]이므로 위 등식을 만족하려면
[math(\displaystyle \begin{aligned} \frac{\partial L}{\partial t}=0 \end{aligned})]
이다.
라그랑지언의 시간 미분을 생각해보자.
[math(\displaystyle \begin{aligned} \frac{{\rm d}L}{{\rm d}t}=\sum_{j} \left[ \frac{\partial L}{\partial q_{j}} \dot{q}_{j}+\frac{\partial L}{\partial \dot{q}_{j}} \ddot{q}_{j} \right] \end{aligned})]
오일러-라그랑주 방정식에 의해
[math(\displaystyle \begin{aligned} \frac{\partial L}{\partial q_{j}}=\frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{j}} \end{aligned})]
제 1항을 교체 하면
[math(\displaystyle \begin{aligned} \frac{{\rm d}L}{{\rm d}t}&=\sum_{j} \left[ \frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{j}} \dot{q}_{j}+\frac{\partial L}{\partial \dot{q}_{j}} \ddot{q}_{j} \right] \\ &=\sum_{j} \frac{{\rm d}}{{\rm d}t}\left[\frac{\partial L}{\partial \dot{q}_{j}} \dot{q}_{j} \right] \end{aligned})]
이상에서
[math(\displaystyle \begin{aligned} \frac{{\rm d}}{{\rm d}t}\left[L-\sum_{j}\frac{\partial L}{\partial \dot{q}_{j}} \dot{q}_{j} \right]=0 \end{aligned})]
이에 각괄호 안의 양은 한 상수 임을 알 수 있는데, 그 상수를 [math(-H)][14]라 하자.
[math(\displaystyle \begin{aligned} L-\sum_{j}\frac{\partial L}{\partial \dot{q}_{j}} \dot{q}_{j}=-H \end{aligned})]
그런데 일반적으로 퍼텐셜 에너지는 일반화 속도에 의존하지 않기 때문에
[math(\displaystyle \begin{aligned} L-\sum_{j}\frac{\partial T}{\partial \dot{q}_{j}} \dot{q}_{j}=-H \end{aligned})]
이고, 계가 스클로노믹하다면
[math(\displaystyle \begin{aligned} \sum_{j}\frac{\partial T}{\partial \dot{q}_{j}} \dot{q}_{j}=2T \end{aligned})]
임을 위에서 증명했다.
[math(\displaystyle \begin{aligned} L-2T=-H \end{aligned})]
라그랑지언의 정의 [math(L=T-U)]를 이용하면 다음을 얻는다.
[math(\displaystyle \begin{aligned} T+U=E=H \end{aligned})]
즉, 라그랑지언의 시간 이동에 대한 불변은 에너지의 보존을 이끌어낸다.
그러나 해밀토니언 [math(H)]가 다음의 조건을 만족할 때만 총 에너지와 같음에 유의해야 한다.
- 직교 좌표와 일반화 좌표 간의 변환 공식에 시간이 포함되지 않아야 한다.
- 퍼텐셜 에너지가 속도에 의존하지 않아야 한다.
9.2. 일반화 운동량 보존
라그랑지언의 무한소 변위의 병진에 대한 불변은 일반화 운동량 보존을 이끌어낸다. [math(\delta \mathbf{q})]의 병진에 불변이라고 하면[math(\displaystyle \begin{aligned} L(q_{i},\,\dot{q}_{i},\,t)&=L(q_{i}+\delta q_{i},\,\dot{q}_{i},\,t+\delta t) \\&=L(q_{i},\,\dot{q}_{i},\,t)+\sum_{i} \frac{\partial L}{\partial q_{i}}\delta q_{i}\\&=L(q_{i},\,\dot{q}_{i},\,t)+\boldsymbol{\nabla}L \boldsymbol{\cdot} \delta \mathbf{q} \end{aligned})]
시간 이동과 마찬가지로 위 등식이 만족하려면
[math(\displaystyle \begin{aligned} \boldsymbol{\nabla}L =\mathbf{0} \quad \Rightarrow \quad \frac{\partial L}{\partial q_{i}}=0 \end{aligned})]
오일러-라그랑주 방정식에 의하면
[math(\displaystyle \begin{aligned} \frac{\partial L}{\partial q_{i}}=\frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{i}} \end{aligned})]
이므로
[math(\displaystyle \begin{aligned} \frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{q}_{i}}=0 \end{aligned})]
일반화 운동량은 [math(p_{i}={\partial L}/{\partial \dot{q}_{i}})]로 정의되므로
[math(\displaystyle \begin{aligned} \dot{p}_{i}=0 \end{aligned})]
무한소 변위의 병진은 일반화 운동량의 보존을 이끌어 낸다.
일반화 운동량은 곧 뉴턴 역학에서 선운동량이나 각운동량을 나타내므로 이들의 보존과도 연관이 있다.
덤으로 라그랑지언에 [math(q_{j})]가 직접적으로 표현되지 않을 때도 이 일반화 운동량은 보존되는데, 이때의 [math(q_{j})]를 순환 좌표라 한다.
9.2.1. 선운동량 보존
이는 일반화 좌표로 직교 좌표를 택한다면[math(\displaystyle \begin{aligned} \frac{{\rm d}}{{\rm d}t}\frac{\partial L}{\partial \dot{x}_{i}}=0 \end{aligned})]
뉴턴 역학에서 퍼텐셜 에너지는 일반적으로 속도에 의존하지 않으므로 라그랑지언의 속도 미분은 운동 에너지와 관련이 된다. 운동 에너지
[math(\displaystyle \begin{aligned} T=\frac{1}{2}m|\mathbf{\dot{r}}|^{2} \end{aligned})]
임을 이용하면
[math(\displaystyle \begin{aligned} \frac{{\rm d}}{{\rm d}t}(m\dot{x}_{i})=0 \end{aligned})]
로 선운동량 보존을 이끌어 낸다.
9.2.2. 각운동량 보존
각운동량 보존은 살짝 증명하기 어렵다. 이때는 무한소 회전을 고려하는데,
[math(\displaystyle \begin{aligned} \delta \mathbf{r}= \delta \boldsymbol{\phi} \times \mathbf{r} \end{aligned})]
이 변환은 속도 벡터에 대해서도 동일하게 이루어진다.
[math(\displaystyle \begin{aligned} \delta \mathbf{\dot{r}}= \delta \boldsymbol{\phi} \times \mathbf{\dot{r}} \end{aligned})]
라그랑지언의 무한소 회전을 고려했을 때, 이것이 동등하다면 각운동량은 보존된다. 즉,
|
[math(\displaystyle \begin{aligned} L(x_{i},\,\dot{x}_{i},\,t)&=L(x_{i}+\delta x_{i},\,\dot{x}_{i}+\delta \dot{x}_{i},\,t) \\&=L(x_{i},\,\dot{x}_{i},\,t)+\sum_{i} \left[ \frac{\partial L}{\partial x_{i}}\delta x_{i}+\frac{\partial L}{\partial \dot{x}_{i}}\delta \dot{x}_{i} \right] \end{aligned})] |
[math(\displaystyle \begin{aligned} \sum_{i} \left[ \frac{\partial L}{\partial x_{i}}\delta x_{i}+\frac{\partial L}{\partial \dot{x}_{i}}\delta \dot{x}_{i} \right]=0 \end{aligned})]
한편 다음의 관계를 기억한다면
[math(\displaystyle \begin{aligned} \frac{\partial L}{\partial \dot{x}_{i}}&=p_{i} \\\frac{\partial L}{\partial {x}_{i}}&=\dot{p}_{i} \end{aligned})]
이것을 대입하면
[math(\displaystyle \begin{aligned} \sum_{i} [ \dot{p}_{i} \delta x_{i}+p_{i}\delta \dot{x}_{i} ]= \mathbf{\dot{p}} \boldsymbol{\cdot}(\delta \boldsymbol{\phi} \times \mathbf{r} )+\mathbf{p} \boldsymbol{\cdot}(\delta \boldsymbol{\phi} \times \mathbf{\dot{r}} ) \end{aligned})]
벡터의 삼중곱은 배열 순서를 유지한 채로 교환 가능하다.
|
[math(\displaystyle \begin{aligned} \mathbf{\dot{p}} \boldsymbol{\cdot}(\delta \boldsymbol{\phi} \times \mathbf{r} )+\mathbf{p} \boldsymbol{\cdot}(\delta \boldsymbol{\phi} \times \mathbf{\dot{r}} )&=\delta \boldsymbol{\phi} \boldsymbol{\cdot}( \mathbf{r} \times \mathbf{\dot{p}})+\delta \boldsymbol{\phi}\boldsymbol{\cdot}( \mathbf{\dot{r}} \times \mathbf{p} )\\ &= \delta \boldsymbol{\phi} \boldsymbol{\cdot} \frac{{\rm d}}{{\rm d}t}(\mathbf{r\times p})\\ &= \delta \boldsymbol{\phi} \boldsymbol{\cdot} \mathbf{\dot{L}} \end{aligned})] |
[math(\displaystyle \begin{aligned} \mathbf{\dot{L}}=\mathbf{0} \end{aligned})]
즉, 공간의 등방성은 각운동량 보존을 이끌어 낸다.
10. 특수 상대성 이론에서 라그랑지언
의
상대론적 라그랑지언
부분을
참고하십시오.11. 비고전역학적인 확장
원칙적으로 뉴턴 역학의 범위에서 벗어나는 어떠한 물리 현상도 라그랑주 역학이 예측해 주지는 않는다.[15] 그럼에도 불구하고 라그랑지언, 액션, 라그랑주 역학에 대한 공부로부터 물리학과의 커리큘럼을 시작하는 궁극적인 이유는 고전역학을 벗어나는 레벨의 물리에서 매우 다양하게 유사한 방법론을 도입하여 사용하게 되기 때문이다.상대성 이론에서는 시간과 공간이 동일한 방식으로 다루어져야 하므로 액션을 라그랑지언의 시간 적분으로 표현하는 것이 적절하지 않다. 이를 위해서 라그랑지언 밀도 [math(\boldsymbol{\mathcal{L}})]라는 개념이 흔히 사용되는데, 기본적으로 라그랑지언 밀도의 3차원적인 공간적분이 보통의 라그랑지언이 된다.
[math( \displaystyle \begin{aligned} L(t) &= \int \mathcal{L} (\mathbf{r},\,t ) \, {\rm d} \mathbf{r} \\ S &= \int \mathcal{L} ( x ) \,{\rm d}x^\mu \end{aligned})]
위의 마지막 수식, 그리고 이 아래 수식들에서는 사차원 벡터 (4-vector)를 이탤릭체로 표기하는 방법을 사용한다. 혼동될 여지가 없는 경우엔 이 라그랑지언 밀도를 그냥 라그랑지언이라고 부르는 경우도 많다. 다음은 각 양자장론의 분야별로 출발점을 주는 라그랑지언 밀도와 그것에 오일러-라그랑주 방정식을 적용하여 유도한 운동방정식의 리스트이다.[16]
11.1. 비상대론적인 양자장론
[math(\displaystyle \begin{aligned} L (\psi, \, \boldsymbol{\nabla} \psi, \, \dot{\psi} ) &= i \hbar \psi^{\ast} {\partial \psi \over \partial t} - {\hbar^2 \over 2m} \boldsymbol{\nabla} \psi^{\ast} \boldsymbol{\nabla} \psi -V (\mathbf{r},\,t ) \psi^* \psi \\ i \hbar \frac{\partial}{\partial t} \psi(\mathbf{r},\,t )&= \left[ -\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf{r},\,t ) \right] \psi(\mathbf{r},\,t ) \end{aligned} )]
11.2. 스핀 0인 중성자유입자의 양자장론
(이 문단 아래로는 [math( \epsilon_0 = \mu_0 = c = \hbar =1 )] 로 놓는 자연 단위계를 사용한다.)[math(\displaystyle \begin{aligned} \mathcal{L} &= {1 \over 2} (\partial^\mu \phi ) (\partial_\mu \phi ) - {1 \over 2} m^2 \phi^2 \\ 0 &= (\partial^\mu \partial_\mu + m^2 ) \phi \end{aligned} )]
위의 방정식을 클라인-고든 방정식이라고 부른다.[17]
11.3. 스핀 1/2인 자유입자의 양자장론
[math(\displaystyle \begin{aligned} \mathcal{L} &= \bar{\psi} ( i \gamma^\mu \partial_\mu - m ) \psi \\ 0&= ( i \gamma^\mu \partial_\mu - m ) \psi \end{aligned} )]
여기에서의 [math( \psi )]는 [math( 4 \times 1 )] 크기의 행렬로 표현되는 스피너(spinor)이고, [math( \gamma^\mu )]는 디랙 행렬이다. 위의 방정식을 디랙 방정식이라고 부른다.
11.4. 스핀 1인 자유입자의 양자장론
[math(\displaystyle \begin{aligned} \mathcal{L} &= -{1 \over 4} F^{\mu \nu} F_{\mu \nu}+ {1 \over 2} m^2 A^\mu A_\mu F^{\mu \nu} \\ &= \partial^\mu A^\nu - \partial^\nu A^\mu \\ 0&= \partial^\nu \partial_\nu A^\mu -\partial^\mu \partial_\nu A^\nu +m^2 A^\mu \end{aligned})]
여기에서의 [math( A^\mu )]가 4차원 퍼텐셜, [math( F^{\mu \nu} )]는 전자기 텐서이다. 보통 아래의 운동방정식의 우변에 0 대신 4차원 전류항 [math( j^\mu )]를 넣은 것을 프로카 방정식이라고 부른다. 위의 방정식에서 로런츠 게이지(Lorenz gauge) [math( \partial^\mu A_\mu =0)]를 적용하여 게이지 고정(gauge fixing)을 하면 위에서 스칼라 장의 운동 방정식과 똑같은 모양이된다. 즉 운동 방정식은 로런츠 게이지에서 그 자체로도 맥스웰 방정식이면서 클라인-고든 방정식이기도 하다.
참고로 질량 [math(m)]이 0이 아닐 때 프로카 방정식은 게이지 변환에 대해 불변이 아니다. 따라서 게이지 불변성을 만족시키려면 광자의 질량이 0이어야 하며 동시에 전류 밀도(current density) [math( j^\mu )]는 0이어야 한다.
여기서 조금만 더 나아가자. 전하를 띤 페르미온[18]의 전자기적 상호작용을 기술하기 위해서는 최소한 자유롭게 움직이는 페르미온과 광자에 대한 라그랑지언이 필요하고, 전하를 띈 입자가 움직이고 있으므로 전류 밀도도 필요하다. 아무런 조건 없이 전류 밀도를 고려하여 라그랑지언의 게이지 변환을 가하면 라그랑지언이 게이지 불변을 만족시키지 않는 것으로 보여, 잘못된 이론을 구상한 것이 아닌가 생각해 볼 수 있다. 그러나 Dirac 장만을 사용해 전류밀도를 나타내고 전자기장의 게이지 변환이 전자기장만 바꾸는 것이 아닌 디락장의 phase를 바꾸는 변환으로 작동한다면, 전류 밀도가 있어도 전체 라그랑지언의 꼴은 게이지 변환에 대해서 불변하다.
즉, 라그랑지언을 구성하는 전류 밀도가 게이지 불변이기 위해서는 광자의 게이지 변환외에 전하를 띤 페르미온(디락장)의 게이지 변환도 동시에 작동해야 한다는 것이다. 전류 밀도가 다른 장의 게이지 변환으로 보완되어 라그랑지언의 항으로 표현될 수 있다는 것을 통해, 전류 밀도를 허용하지 말아야 할 것으로 보였던 문제가 해결된다.
그럼에도 질량이 있으면서 입자와 상호작용하는 벡터장(전자기장도 벡터장의 일부이다.)을 어떻게 기술할 것인가란 문제는 여전히 해결되지 않고 남는다. 그러나 전류 밀도가 존재하기 위해서 Dirac 장의 게이지 변환을 고려한 것과 마찬가지 방법을 활용해, 광자의 게이지 변환과 추가 field의 게이지변환을 고려하여 적절하게 질량항이 존재할 수 있게 만들수 있지 않을까 생각해 볼수 있을 것이다. 이 추가 장이 바로 힉스이다.[19]
질량항은 광자와 Dirac 장으로 구축된 양자 전자기학에서 사용하는 게이지(일명 [math( U(1))] 게이지) 변환보다 더 큰 개념의 게이지 변환으로 부터 얻을 수 있다. 대표적으로 표준모형의 [math( SU(2)_L \times U_Y(1))] [20] 게이지 불변을 예로 들수 있다. 물론 그냥 게이지 불변이기만 하면 질량항을 여전히 만들 수는 없으나, [math( SU_L(2) \times U_Y(1))] 게이지군이 자발적 대칭성 붕괴로 인해 전자기학의 게이지 군인 [math( U_e(1))]으로 변화하면[21] 줄어든 자유도만큼 게이지 장들의 질량항이 늘어난다. 이것이 그 유명한 Higgs mechanism이다.[22]
이를 통해서 (전자기장과 같이) 스핀이 1인 벡터장의 질량항은 게이지 장이 만족시키는 게이지 변환보다 더 큰 구조의 게이지 변환으로부터 불변성을 부여할 수 있으므로, 질량항이 존재하면 안 될 것 같다는 문제도 해결된다.[23]
11.5. 양자 전기역학
[math(\displaystyle \mathcal{L} = \bar{\psi} \left( i \gamma^\mu \partial_\mu - m \right) \psi - {1 \over 4} F_{\mu \nu} F^{\mu \nu} - q \bar{\psi} \gamma^\mu \psi A_\mu )]
[math( \psi )] 와 [math( A^\mu )] 각 방향의 오일러-라그랑주 방정식으로부터 전하 [math( q )]를 가진 스핀 1/2입자(전자와 양전자)와 스핀 1입자(광자)의 운동 방정식이 튀어나오고, 자연스럽게 전류의 정체가 무엇인지, 그것을 통해 두 종류의 입자가 어떻게 서로 영향을 주고받는지까지 결정된다.
양자 전기역학의 라그랑지언이 저런 꼴을 갖도록 하는 원인 중 하나로 게이지 불변성이 있다. 페르미온 장의 위상 변화 [math(\psi \to e^{i \phi} \psi)]에서 [math(\phi)]가 상수이면 자동으로 불변성을 만족시키던 것을 [math(\phi)]가 상수가 아닌 경우로 확장하는 과정에서 게이지 불변성을 만족시키는 스핀 1의 장이 튀어나오고 그 장이 다름 아닌 전자기장인 것을 볼 수 있다. 즉, 전자기장이 존재하는 이유가 게이지 불변성 때문이라는 것이다. 그리고 이 논리를 여러 페르미온들의 다중항에 대한 일반적인 위상 변화에 대한 것으로 확장할 수도 있다. 중력을 제외한 모든 상호작용을 다 설명하는 양-밀스 장이 바로 그것. 자세한 내용은 게이지 장 참고.
12. 예제
13. 관련 문서
[1]
Méchanique analitique , par M. de La Grange,(Lagrange, Joseph-Louis) Auteur du texte 1788 P512,
#
[2]
스칼라는 좌표 변환에 불변이기 때문이다.
[3]
신의 지혜로 작용이 최소화된다는 것이 모페르튀이(Maupertuis)의 원리다.
[4]
밑에 설명되어 있지만 먼저 이야기를 하자면 반드시 위치와 시간을 라그랑지언의 변수로 택할 필요는 없다. 다만 위치역할을 수행할수 있는 함수와 속도역할을 수행할 수 있는 함수라면 어떤 변수를 쓰든 관계없다.
[5]
대수적으로 [math( dv/dx )]는 0이 아니지만 [math( \partial v/\partial x )]는 관점에 따라서 0으로 치부할 수 있다는 점을 떠올리면 편하다.
[6]
라그랑주 역학으로 따지면 최초의 시간 · 위치 좌표
[7]
최후의 시간-위치 좌표
[8]
또는
맥스웰 방정식,
클라인-고든 방정식,
디랙 방정식 등:
오일러-라그랑주 방정식과 같은 역할을 한다.
[9]
양자역학은 라그랑주 역학으로 해석된다는 것을 밝힌 사람이 파인만이며, 이 방법은
경로적분이라고 불린다.
[10]
구별을 위해 흘려쓴 서체인 [math(\mathcal{L})]이나 [math(\mathscr{L})]로 쓰기도 한다.
[11]
사실, 현대물리학에서 라그랑지언은 하술할 최소 작용의 원리에 따라, 오일러-라그랑주 방정식을 통해 '계의 방정식'을 이끌어낼 수 있는 물리량으로 정해지는 경우가 많다. 즉, 어떤 물리 현상들의 물리량들의 결과로 정의되는 게 아니라, 라그랑지언이 더 근본적인 물리량이고 거기서
최소 작용의 원리에 따라 각종 방정식들이 도출된다고 보는 편이 타당하다.
[12]
이를테면, 3차원의 직교 좌표계를 고려하자. 이 경우엔 일반적으로 위치의 좌표를 [math(x, \, y,\,z)]로 나타낼 수 있다. 그러나 이것을 [math(x_{1}\equiv x,\,x_{2}\equiv y, x_{3}\equiv z)]로 쓰자는 것이다.
[13]
예를 들어, 반지름이 [math( R )]인 원운동하는 입자는 구속 조건이 [math( x^2 + y^2 = R^2 )]이므로 [math( f = x^2 + y^2 - R^2 = 0 )]로 쓸 수 있다. 또한, 두 물체가 연결된 작은 도르래는 양쪽 줄의 길이를 [math( X, Y )]라 하고 줄의 총 길이를 [math( l )]이라고 하면 [math( X+Y=l )]이므로 [math( f = X + Y - l = 0 )]로 쓸 수 있다.
[14]
사실 각괄호 안은 라그랑지언의
르장드르 변환인
해밀토니언이다.
[15]
현실적으로는 각각의 해석법이 더 유리한 경우들이 있으므로, 그때 그때 필요한 만큼 라그랑주 역학을 사용하게 된다. 그래도 실제 물리학과 학부 2학년 정도로 들어서면 뉴턴 역학으로 풀기 전에 라그랑주 역학을 도입해보는 게 효율적이다.
[16]
실제로 양자장을 만들기 위해서는 이차 양자화 과정을 거쳐야 한다.
[17]
Field의 Lagrangian 해는 유일하지 않다. 다만 사람들이 주로 많이 쓰는, 혹은 물리적 해석에 더욱 용이한 Lagrangian만을 쓸 뿐이며, 지면에 쓰인 Lagrangian은 여러책이나 논문에서 많이 사용하는 Lagrangian이다.
[18]
스핀이 반정수배인 입자. 대표적으로 전자가 있다.
[19]
한가지 예제로 ablian Higgs mechanism에서 계산을 진행할 경우, 광자는 질량을 갖는다. 표준모형은 Non-abilan Higgs model에 따른다. 이 첨삭에서 언급한 Non-abelian Higgs model은 다음 문단에 개략적으로 설명되어 있다.
[20]
[math( SU_L(2) \times U_Y(1) )] 게이지 불변이란 말에서 [math( SU_L(2))]는 표준모형상 짝이 되는 spin 1/2 Dirac field들(예를 들어, 중성미자와 전자, 또는 up quark와 down quark)중에 왼손잡이(Left chiral, Left handed가 아님을 주의하자.)만을 떼어내서 doublet을 구성하였을 때, [math( SU_L(2) )] 군의 요소에 대한 변환에 불변하다는 의미의 게이지 군이다. [math( U_Y(1))]은 전자기학에서 사용하는 게이지 변환과 유사하게 Hypercharge [math( Y )]라는 값에 따라 spin 1/2 Dirac field들의 phase를 바꾸는 변환이다.
[21]
[math(SU_L(2)\times U_Y(1) )] 군의 부분군으로 [math(U_e(1))]를 구성할 수 있다. 그러나 [math(U_e(1))]의 요소에 따른 변환에 대해 라그랑지언이 불변하기를 요구하면 남은 세 요소의 게이지 변환에 대해서 더 이상 불변하지 않게 된다.
[22]
자발적 대칭성 붕괴후에도 남은, 여전히 게이지 불변이 성립하는 게이지 장은 Massless, 질량이 없다.
[23]
다만 앞선 전류 밀도의 문제와는 해결방법이 다른 것이, 자발적 대칭성 붕괴와 진공 기댓값을 요구한다는 점이다. 이를 통해서 자발적 대칭성 붕괴로 인해 줄어드는 게이지 변환의 갯수(또는 변환에 대한 자유도의 개수)가 massive gauge field 종류의 개수를 결정하기 때문에, Higgs mechanism은 퍼텐셜이 어떤 함수인가에 따라 나타나는 결론이 다르다.