질량이 비슷한 두 천체가 둘의 질량 중심 주위를 타원궤도로 돌고 있다. 고전역학 에서 이체 문제 (二體問題, two-body problem )는 서로 상호작용하는 두 물체의 운동을 다루는 문제이다. 보통 이 상호작용은 만유인력 과 같은 역제곱 법칙 이다. 행성 을 공전하는 위성 , 항성 을 공전하는 행성, 쌍성 계 등이 이체 문제에 해당한다.
이체 문제는 식의 변형을 통해 두 개의 독립적인 일체문제 로 변형될 수 있다. 일체문제는 퍼텐셜 안에서의 한 물체의 운동에 관한 식을 푸는 것으로 기술되며, 언제나 해석적인 해를 구할 수 있다. 그러므로 이체 문제도 언제나 해를 구할 수 있다. 그러나 삼체문제 이상의 다체문제 에서는 특수한 경우를 제외하고는 해석적인 해를 구할 수 없고, 수치적인 방법을 쓴다.
질량이 약간 차이나는 두 천체가 둘의 질량 중심 주위를 타원궤도로 돌고 있다. 이 형태는 명왕성 -카론 계나 지구-달 계와 비슷하다. 단지 지구-달 계는 질량 중심이 지구 내부에 있다는 차이가 있다. 이체 문제에 관한 연구는 천체 역학으로부터 비롯되었다. 이체 문제의 역사는 독일의 천문학자 요하네스 케플러 로부터 시작된다. 케플러는 덴마크의 천문학자 튀코 브라헤 의 정밀 천체 관측 자료를 분석하여 행성의 운동에 관한 케플러의 법칙 을 발견하였다. 이후 영국의 물리학자 아이작 뉴턴 , 에드먼드 핼리 등은 천체 간의 거리의 제곱에 반비례하는 힘 으로 케플러의 법칙 을 설명할 수 있다는 사실을 발견하고, 뉴턴은 이를 이용해 유명한 저서 프린키피아 에서 수학적으로 행성들이 타원 궤도를 돈다는 것을 보였다. 또한, 뉴턴은 행성의 질량, 위치, 속도 정보를 이용하여 궤도를 계산하는 법을 개발했다. 뉴턴 이후 레온하르트 오일러 는 '행성과 혜성의 운동'에서 오일러 방법 을 이용하여 천체의 포물선 운동을 수학적으로 보였다. 핼리는 뉴턴의 법칙을 통해 핼리 혜성 등 혜성 의 주기를 알아내는 데 성공했으며 이후 요한 람베르트 는 오일러 방법을 일반화하여 행성의 타원, 쌍곡선 궤도를 수학적으로 계산해내었다. 이후 조제프루이 라그랑주 는 행성 간의 각도들을 이용하여 궤도를 계산하는 새로운 방법을 개발했다. 이 방법은 카를 프리드리히 가우스 에 의해 훨씬 개선되었고, 현재도 계속 사용되고 있다.
이체문제 해결을 위한 야코비 좌표계. 야코비 좌표계는
M
=
m
1
+
m
2
{\displaystyle M=m_{1}+m_{2}}
R
=
m
1
M
x
1
+
m
2
M
x
2
{\displaystyle {\boldsymbol {R}}={\frac {m_{1}}{M}}{\boldsymbol {x}}_{1}+{\frac {m_{2}}{M}}{\boldsymbol {x}}_{2}}
r
=
x
1
−
x
2
{\displaystyle {\boldsymbol {r}}={\boldsymbol {x}}_{1}-{\boldsymbol {x}}_{2}}
[ 1] 두 물체 각각의 위치를
x
1
{\displaystyle \mathbf {x} _{1}}
x
2
{\displaystyle \mathbf {x} _{2}}
m
1
{\displaystyle \mathbf {m} _{1}}
m
2
{\displaystyle \mathbf {m} _{2}}
x
1
(
t
=
0
)
{\displaystyle \mathbf {x} _{1}(\mathbf {t} =0)}
x
2
(
t
=
0
)
{\displaystyle \mathbf {x} _{2}(\mathbf {t} =0)}
v
1
(
t
=
0
)
{\displaystyle \mathbf {v} _{1}(\mathbf {t} =0)}
v
2
(
t
=
0
)
{\displaystyle \mathbf {v} _{2}(\mathbf {t} =0)}
x
1
(
t
)
{\displaystyle \mathbf {x} _{1}(\mathbf {t} )}
x
2
(
t
)
{\displaystyle \mathbf {x} _{2}(\mathbf {t} )}
t
{\displaystyle \mathbf {t} }
각각의 물체에 대하여 뉴턴의 운동 제2 법칙 을 적용해보면
F
12
(
x
1
,
x
2
)
=
m
1
x
¨
1
(
1
)
{\displaystyle \mathbf {F} _{12}(\mathbf {x} _{1},\mathbf {x} _{2})=m_{1}{\ddot {\mathbf {x} }}_{1}\quad \quad \quad (1)}
F
21
(
x
1
,
x
2
)
=
m
2
x
¨
2
(
2
)
{\displaystyle \mathbf {F} _{21}(\mathbf {x} _{1},\mathbf {x} _{2})=m_{2}{\ddot {\mathbf {x} }}_{2}\quad \quad \quad (2)}
F
12
{\displaystyle \mathbf {F} _{12}}
F
21
{\displaystyle \mathbf {F} _{21}}
위의 두 식을 더하고 빼는 과정을 통해 이체 문제를 두 개의 독립적인 일체 문제로 분리할 수 있다. (1)과 (2)를 더하면 이 계의 질량 중심 의 운동을 시간에 대해 나타내는 식이 나오고, 두 식을 빼면 두 물체를 잇는 벡터
r
=
x
1
−
x
2
{\displaystyle \mathbf {r} =\mathbf {x} _{1}-\mathbf {x} _{2}}
x
1
(
t
)
{\displaystyle \mathbf {x} _{1}(\mathbf {t} )}
x
2
(
t
)
{\displaystyle \mathbf {x} _{2}(\mathbf {t} )}
(1)과 (2)를 더하여 아래의 식을 얻을 수 있다.
m
1
x
¨
1
+
m
2
x
¨
2
=
(
m
1
+
m
2
)
R
¨
=
F
12
+
F
21
=
0
{\displaystyle m_{1}{\ddot {\mathbf {x} }}_{1}+m_{2}{\ddot {\mathbf {x} }}_{2}=(m_{1}+m_{2}){\ddot {\mathbf {R} }}=\mathbf {F} _{12}+\mathbf {F} _{21}=0}
여기서 뉴턴의 운동 제3 법칙 에 의하여
F
12
=
−
F
21
{\displaystyle \mathbf {F} _{12}=-\mathbf {F} _{21}}
R
¨
≡
m
1
x
¨
1
+
m
2
x
¨
2
m
1
+
m
2
{\displaystyle {\ddot {\mathbf {R} }}\equiv {\frac {m_{1}{\ddot {\mathbf {x} }}_{1}+m_{2}{\ddot {\mathbf {x} }}_{2}}{m_{1}+m_{2}}}}
R
{\displaystyle \mathbf {R} }
결과적으로는
R
¨
=
0
{\displaystyle {\ddot {\mathbf {R} }}=0}
이다.
이 식을 통해 질량 중심 의 속도
V
=
d
R
/
d
t
{\displaystyle \mathbf {V} =\mathbf {d} \mathbf {R} /\mathbf {d} \mathbf {t} }
운동량
m
1
v
1
+
m
2
v
2
{\displaystyle \mathbf {m} _{1}\mathbf {v} _{1}+\mathbf {m} _{2}\mathbf {v} _{2}}
운동량 보존 법칙 ). 그러므로 질량 중심 의 변위
R
(
t
)
{\displaystyle \mathbf {R} (\mathbf {t} )}
(1)과 (2)를 각각 물체의 질량으로 나누어주고, (1)에서 (2)를 뺀 후에 정리하면
r
¨
=
x
¨
1
−
x
¨
2
=
(
F
12
m
1
−
F
21
m
2
)
=
(
1
m
1
+
1
m
2
)
F
12
{\displaystyle {\ddot {\mathbf {r} }}={\ddot {\mathbf {x} }}_{1}-{\ddot {\mathbf {x} }}_{2}=\left({\frac {\mathbf {F} _{12}}{m_{1}}}-{\frac {\mathbf {F} _{21}}{m_{2}}}\right)=\left({\frac {1}{m_{1}}}+{\frac {1}{m_{2}}}\right)\mathbf {F} _{12}}
여기서 뉴턴의 운동 제3 법칙
F
12
=
−
F
21
{\displaystyle \mathbf {F} _{12}=-\mathbf {F} _{21}}
r
{\displaystyle \mathbf {r} }
두 물체 간에 작용하는 힘은 두 물체 간의 상호작용에 의해서 만들어진 것이기 때문에
x
1
{\displaystyle \mathbf {x} _{1}}
x
2
{\displaystyle \mathbf {x} _{2}}
r
{\displaystyle \mathbf {r} }
μ
r
¨
=
F
12
(
x
1
,
x
2
)
=
F
(
r
)
{\displaystyle \mu {\ddot {\mathbf {r} }}=\mathbf {F} _{12}(\mathbf {x} _{1},\mathbf {x} _{2})=\mathbf {F} (\mathbf {r} )}
여기서
μ
{\displaystyle \mu }
환산 질량
μ
=
1
1
m
1
+
1
m
2
=
m
1
m
2
m
1
+
m
2
{\displaystyle \mu ={\frac {1}{{\frac {1}{m_{1}}}+{\frac {1}{m_{2}}}}}={\frac {m_{1}m_{2}}{m_{1}+m_{2}}}}
R
(
t
)
{\displaystyle \mathbf {R} (\mathbf {t} )}
r
(
t
)
{\displaystyle \mathbf {r} (\mathbf {t} )}
x
1
(
t
)
=
R
(
t
)
+
m
2
m
1
+
m
2
r
(
t
)
{\displaystyle \mathbf {x} _{1}(t)=\mathbf {R} (t)+{\frac {m_{2}}{m_{1}+m_{2}}}\mathbf {r} (t)}
x
2
(
t
)
=
R
(
t
)
−
m
1
m
1
+
m
2
r
(
t
)
{\displaystyle \mathbf {x} _{2}(t)=\mathbf {R} (t)-{\frac {m_{1}}{m_{1}+m_{2}}}\mathbf {r} (t)}
이 식들은
R
{\displaystyle \mathbf {R} }
r
{\displaystyle \mathbf {r} }
이체 문제에서 두 물체는 언제나 한 평면 위에서 움직이게 된다. 운동량
p
{\displaystyle \mathbf {p} }
각운동량
L
{\displaystyle \mathbf {L} }
L
=
r
×
p
=
r
×
μ
d
r
d
t
{\displaystyle \mathbf {L} =\mathbf {r} \times \mathbf {p} =\mathbf {r} \times \mu {\frac {d\mathbf {r} }{dt}}}
L
{\displaystyle \mathbf {L} }
돌림힘
N
{\displaystyle \mathbf {N} }
N
=
d
L
d
t
=
r
˙
×
μ
r
˙
+
r
×
μ
r
¨
,
{\displaystyle \mathbf {N} ={\frac {d\mathbf {L} }{dt}}={\dot {\mathbf {r} }}\times \mu {\dot {\mathbf {r} }}+\mathbf {r} \times \mu {\ddot {\mathbf {r} }}\ ,}
방향이 같은 두 벡터의 외적은 0이라는 외적의 성질을 이용하면
N
=
d
L
d
t
=
r
×
F
,
{\displaystyle \mathbf {N} \ =\ {\frac {d\mathbf {L} }{dt}}=\mathbf {r} \times \mathbf {F} \ ,}
이고,
F
=
μ
d
2
r
/
d
t
2
{\displaystyle \mathbf {F} =\mu \mathbf {d} ^{2}\mathbf {r} /\mathbf {dt} ^{2}}
두 물체 사이의 힘이 두 물체를 잇는 선과 같은 방향이라고 가정하면,
r
{\displaystyle \mathbf {r} }
F
=
0
{\displaystyle \mathbf {F} =0}
각운동량 벡터
L
{\displaystyle \mathbf {L} }
각운동량 보존 ) 그러므로 변위 벡터
r
{\displaystyle \mathbf {r} }
v
{\displaystyle \mathbf {v} }
L
{\displaystyle \mathbf {L} }
대부분의 물리 문제에서 중심력
F
(
r
)
{\displaystyle \mathbf {F} (\mathbf {r} )}
F
(
r
)
=
F
(
r
)
r
^
{\displaystyle \mathbf {F} (\mathbf {r} )=F(r){\hat {\mathbf {r} }}}
여기서 r = |
r
{\displaystyle \mathbf {r} }
r
^
{\displaystyle {\hat {\mathbf {r} }}}
r
{\displaystyle \mathbf {r} }
단위벡터 이다. 그러므로
μ
r
¨
=
F
(
r
)
r
^
,
{\displaystyle \mu {\ddot {\mathbf {r} }}={F}(r){\hat {\mathbf {r} }}\ ,}
이고,
F
(
r
)
{\displaystyle \mathbf {F} (\mathbf {r} )}
이체문제에서 두 물체가 서로에게 작용한 힘에 의한 일은 한 힘이 물체 사이의 상대변위 에 곱해진 것과 같다.
이체 문제에서 한 물체가 다른 물체를 관측할 때, 그 다른 물체는 그 물체의 속도와 위치에 따라 원 , 타원 , 포물선 , 쌍곡선 의 네 가지 이차곡선 궤도 중 한 궤도를 그린다.
