역제곱 중심력과 유효 퍼텐셜

궤도를 그리기 전에 도달 가능한 거리를 묻는다

타원 그림이 보인다는 사실만으로 그 운동이 물리 법칙에 맞는다고 할 수 없습니다. 초기 위치와 속도를 정했을 때 입자가 중심에서 얼마나 가까워지고 멀어질 수 있는지 먼저 알아봅시다. 이 질문은 평면 운동을 하나의 반지름 에너지 조건으로 줄여 줍니다. 다만 허용되는 상태를 찾는 일과 특정 시각의 위치를 예측하는 일은 다릅니다.

옆으로 던진 물체를 생각해 보세요. 인력은 안쪽으로 당기지만 옆 방향 속도가 각운동량을 만듭니다. 더 작은 반지름으로 가려면 접선 속력이 커져야 하고, 그 운동에너지가 전체 에너지 예산을 넘으면 접근할 수 없습니다. 유효 퍼텐셜은 이 경쟁을 표현합니다.

모델과 단위를 먼저 고정한다

관성계의 원점에 고정된 힘의 원천을 놓습니다. 입자의 질량 m과 매개변수 μ는 양수이며 일정합니다.

$$\ddot q=-\frac{\mu q}{r^3},\qquad r=|q|>0.$$

μ의 단위는 m³/s²이며 입자의 질량을 바꾸어도 고정합니다. 실제 힘은 이 가속도에 m을 곱한 값입니다. 물리적 퍼텐셜 에너지는 무한대에서 0인 −mμ/r입니다. 원천도 움직이는 두 물체 문제, 항력, 다른 천체, 상대론과 유한 크기 접촉은 포함하지 않습니다. r=0에서 식이 특이하므로 그 점을 통과시키는 규칙도 없습니다.

초기 위치를 q₀=(r₀,0), 초기 속도를 (vᵣ₀,vₜ₀)로 둡니다. 첫 성분은 바깥쪽을 양수로 하는 반지름 속도, 둘째는 부호 있는 접선 속도입니다. 접선 속도가 음수이면 회전 방향이 바뀌며 반지름 자체가 음수가 되는 것은 아닙니다.

보존량 두 개를 유도한다

중심력은 원점에 대한 토크가 0입니다. 따라서 단위 질량당 각운동량 h와 전체 각운동량 L이 보존됩니다.

$$h=q_xv_y-q_yv_x=r^2\dot\theta=r_0v_{t0},\qquad L=mh.$$

가속도 식과 속도를 내적하면 시간에 무관한 퍼텐셜에서 역학적 에너지가 보존됨을 얻습니다.

$$\epsilon=\frac{v^2}{2}-\frac\mu r =\frac{v_{r0}^2+v_{t0}^2}{2}-\frac\mu{r_0}.$$

ε는 단위 질량당 에너지로 단위가 m²/s²입니다. 전체 에너지는 mε이며 J 단위입니다. 같은 μ와 초기 조건에서 m만 바꾸면 전체 에너지와 L은 달라져도 가속도와 궤도는 같습니다. 질량과 μ를 동시에 바꾸는 실험으로 해석하면 안 됩니다.

접선 속도를 제거해 유효 퍼텐셜을 만든다

속력 제곱을 ṙ²+r²θ̇²로 분해하고 θ̇=h/r²를 대입합니다.

$$\frac{\dot r^2}{2}+U_{\rm eff}(r)=\epsilon,\qquad U_{\rm eff}(r)=-\frac\mu r+\frac{h^2}{2r^2}.$$

첫 항은 실제 퍼텐셜을 질량으로 나눈 값입니다. 둘째 항은 보존되는 각운동량을 유지하는 데 필요한 접선 운동에너지입니다. 각운동량 장벽이라고 부르지만 원래 관성계의 힘 식에 새로운 바깥 방향 상호작용을 더한 것이 아닙니다.

이 축약에서 각운동 자체가 사라진 것은 아닙니다. θ̇=h/r²라는 관계에 남아 있습니다. 반지름이 작아지면 같은 h에서 각속도가 커지는 이유도 이 식으로 설명됩니다.

제곱근을 취하기 전에 허용 여부를 판정한다

조회한 양의 반지름에서 ε−Ueff≥0이어야 합니다. 허용되는 경우에는

$$\dot r=\pm\sqrt{2(\epsilon-U_{\rm eff}(r))},\qquad v_t=\frac h r$$

입니다. 전환점 밖에서는 안쪽과 바깥쪽 두 가지 반지름 속도가 가능합니다. 에너지 조건만으로 미래 시각의 분기를 선택하지 않습니다. 그 질문에는 초기 방향과 시간에 따른 운동 풀이가 필요합니다.

차이가 음수이면 허용되지 않는 상태이며 반지름 속도를 0으로 표시하지 않습니다. 등호이면 반지름 속도만 0이고 접선 속도는 남을 수 있습니다. 근점과 먼점에서 입자가 반드시 정지하는 것은 아닙니다.

실험은 초기 보존량과 조회 에너지 조건을 정확한 유리수로 계산합니다. 허용 여부를 임의 오차 허용치로 바꾸지 않으며 제곱근과 공간 그림은 수치 평가임을 구분합니다.

원궤도 조건과 안정성

h가 0이 아닐 때 유효 퍼텐셜의 정지점은

$$r_c=\frac{h^2}{\mu},\qquad U_{\rm eff}(r_c)=-\frac{\mu^2}{2h^2},\qquad U_{\rm eff}''(r_c)=\frac\mu{r_c^3}>0$$

입니다. 양의 곡률은 h를 고정했을 때 반지름 방향의 국소 최소임을 뜻합니다. 원운동이 되려면 에너지까지 이 최솟값과 같아야 합니다. 단순히 rc를 조회한다고 임의 초기 상태가 원운동으로 바뀌지 않습니다.

초기 조건으로는 vᵣ₀=0과 vₜ₀²=μ/r₀가 필요합니다. 그 에너지에서는 허용 반지름이 하나뿐이지만 입자는 접선 방향으로 움직입니다. 역제곱 인력의 원운동 주기는 반지름에 의존합니다. 앞서 다룬 등방 조화력의 진폭과 무관한 진동수와 구분하세요. 모든 중심력이 같은 주기 법칙을 갖는 것은 아닙니다.

이심률 벡터로 궤도의 방향을 복원한다

역제곱 법칙에서는 이심률 벡터도 보존됩니다.

$$\mathbf e=\frac{v^2q-(q\cdot v)v}{\mu}-\frac q r.$$

현재 초기 조건에서는

$$\mathbf e=\left(\frac{r_0v_{t0}^2}{\mu}-1, -\frac{r_0v_{r0}v_{t0}}{\mu}\right),\qquad e^2=1+\frac{2\epsilon h^2}{\mu^2}.$$

h≠0일 때 이 벡터 방향을 기준으로 한 각 f와 p=h²/μ를 사용하면 r=p/(1+e cos f)입니다. 분모는 양수여야 합니다. 음의 에너지는 타원, 0은 포물선, 양수는 쌍곡선이며 e=0은 원입니다. 원에서는 근점 방향이 정해지지 않으므로 그림의 시작 방향을 물리적 불변량으로 해석하지 않습니다.

전환 반지름과 유한한 그림 창

방사 운동이 아닌 경우 가장 가까운 반지름은 rp=p/(1+e)입니다. 이 식은 작은 각운동량에서 비슷한 두 수를 빼는 오차를 피합니다. 음의 에너지에서는 가장 먼 반지름이 ra=−μ/ε−rp입니다. 포물선과 쌍곡선에는 유한한 바깥 전환점이 없습니다.

공간 그림은 시간 경로가 아니라 궤도 기하입니다. 표시 범위는 반지름 12m이며 큰 타원과 비결합 궤도는 그 경계에서 잘립니다. 잘린 끝을 전환점으로 해석하지 않습니다. 쌍곡선 점근선 너머를 선으로 잇거나 분모가 음수인 부분을 그리지 않습니다. 초기 위치도 표본에 명시적으로 포함합니다.

반지름 에너지 그림에는 별도의 양의 반지름 창을 사용합니다. 아래 경계는 초기 반지름·조회 반지름·방사 운동이 아닌 근점 반지름 중 최솟값의 3/4입니다. 원점 가까이의 거대한 값이 에너지 우물을 가리지 않도록 선택한 표시 창이며 물리적 허용 구간 자체가 아닙니다. 위 경계는 12m입니다. 가로축은 log10(r/m)로 한 칸 증가가 반지름 10배를 뜻합니다. 세로축의 질량당 에너지는 변환하지 않습니다. 아주 작은 근점 부근은 조밀하게 표본화하며 표시 창 안의 근점·먼점·원궤도 반지름도 포함합니다. 가로축이 시간인지 반지름인지, 세로축이 전체 에너지인지 질량당 에너지인지 먼저 읽으세요.

각운동량이 0인 경우를 따로 다룬다

h=0이면 장벽이 사라지고 각도를 독립변수로 바꾸는 과정이 성립하지 않습니다. 그런데 불변량 항등식은 모든 방사 운동에서 e=1을 줍니다. 음의 에너지로 떨어지는 입자도 마찬가지입니다. 따라서 e=1이면 포물선이라는 말에는 h≠0 조건이 필요합니다.

음의 에너지인 방사 운동은 바깥 전환 반지름 −μ/ε를 가지며 특이점에 접근할 수 있습니다. 에너지가 0 이상이면 유한한 바깥 전환점이 없습니다. 초기 속도의 방향도 별도로 중요합니다. 실험은 방사 분류와 허용 반지름을 표시하고, 성립하지 않는 원뿔곡선을 만들어 내거나 중심 충돌 이후를 이어 붙이지 않습니다.

예제와 해설이 있는 연습문제

1. 원궤도 등호. μ=r₀=1, vᵣ₀=0, vₜ₀=1로 두세요. 해설: h=1, ε=−1/2, e=0이고 rc=1입니다. 반지름 2에서는 Ueff=−3/8이므로 ε−Ueff=−1/8입니다. 보존량이 같은 상태로 그 반지름에 갈 수 없습니다.

2. 먼점에서 출발하는 타원. μ=r₀=1에서 vₜ₀=1/2, vᵣ₀=0으로 두세요. 해설: h=1/2, ε=−7/8, e=3/4입니다. 시작 반지름 1이 가장 먼 점이고 근점 반지름은 1/7입니다. 반지름 속도가 0이라는 사실만으로 원운동이 되지 않습니다.

3. 포물선 경계. μ=2, r₀=1, vᵣ₀=0, vₜ₀=2로 두세요. 해설: ε=0, h=2, p=2, e=1이며 근점은 1입니다. 유한한 먼점은 없고 분모가 0으로 가는 방향에서 거리가 발산합니다.

4. 두 반지름 분기. μ=r₀=1, vᵣ₀=vₜ₀=1에서 반지름 1을 조회하세요. 해설: ε=0, h=1, Ueff=−1/2이므로 반지름 속도는 −1과 +1 모두 가능합니다. 초기값이 +1이어도 에너지 표면에서 안쪽 분기를 제거하지 않습니다.

5. 방사 퇴화. μ=r₀=1에서 정지 출발하면 어떨까요? 해설: h=0, ε=−1입니다. e=1이어도 비퇴화 포물선이 아니라 묶인 방사 낙하입니다. 가장 먼 반지름은 1이며 조회 반지름 3은 허용되지 않습니다.

6. 질량과 회전 방향. 질량을 네 배로 바꾼 뒤, 별도로 접선 속도 부호를 반전하세요. 해설: 질량 변화는 전체 에너지와 L만 네 배로 만듭니다. 접선 반전은 h의 부호와 이심률 벡터의 y성분을 반전하지만 h²로 정해지는 반지름 허용 조건은 보존합니다.

AI 연구에 연결해 생각하기

이 단원은 보존량, 허용 상태 판정, 시간 예측을 구분합니다. 학습된 모델이 타원처럼 보이는 곡선을 그려도 에너지 예산을 어기거나 금지 반지름을 통과하거나 방사 퇴화를 탈출 궤도로 오인할 수 있습니다. 이들은 서로 다른 실패입니다. 정확한 초기 불변량과 특이한 경계 조건, 방향 반전을 검사하면 그림의 모양만 확인하는 것보다 더 많은 정보를 얻습니다.

다음 단원에서는 결합 타원 운동에 실제 시간을 부여합니다. 유도 과정의 원자료로 MIT OpenCourseWare의 케플러 문제 강의 노트를 참고하세요. 여기의 계산과 예제는 위에서 선언한 고정 원천 모델을 사용합니다.