재미있는 숫자 상식

해왕성의 궤도 섭동 계산: 르베리에는 천왕성의 공전 오차로 위치를 어떻게 찾았나

천왕성의 관측 위치와 뉴턴 역학의 예측 차이를 이용해 르베리에가 미지 행성의 질량·궤도·하늘 위치를 역산하고 해왕성 발견으로 이어진 과정을 설명합니다.

모두의 계산기2026년 9월 22일7분 읽기

해왕성은 먼저 망원경으로 우연히 발견한 뒤 궤도를 계산한 행성이 아닙니다. 천왕성이 예측된 경로에서 조금씩 벗어나는 현상을 분석해 보이지 않는 행성이 어디에 있어야 하는지 계산했고, 그 방향으로 망원경을 돌려 확인했습니다. 관측된 결과에서 원인을 거꾸로 찾아간 대표적인 역문제입니다.

프랑스 수학자 위르뱅 르베리에(Urbain Le Verrier)는 알려진 행성들의 중력 효과를 계산한 뒤에도 남는 천왕성의 위치 오차를 설명할 미지의 외행성을 추정했습니다. 그가 보낸 예측을 바탕으로 요한 고트프리트 갈레와 하인리히 루이 다레스트는 1846년 9월 23일 베를린 천문대에서 해왕성을 확인했습니다. 발견 위치는 르베리에가 계산한 방향에서 약 1도 이내였습니다.

관측한 천왕성 위치 − 알려진 중력으로 계산한 위치 = 설명해야 할 잔차

궤도 섭동은 무엇인가

태양과 행성 하나만 존재한다면 행성은 두 물체 문제의 타원궤도를 따릅니다. 실제 태양계에서는 목성, 토성, 천왕성처럼 여러 행성이 서로를 끌어당깁니다. 다른 행성의 중력이 기준 타원궤도에 만드는 작고 지속적인 변화를 섭동이라고 합니다.

섭동은 행성이 뉴턴의 법칙을 어긴다는 뜻이 아닙니다. 오히려 모든 행성이 동시에 만유인력을 주고받기 때문에 나타납니다. 케플러 법칙과 구심가속도에서 뉴턴의 거리 역제곱 법칙을 유도하는 과정을 여러 천체에 적용하면 각 행성의 위치와 속도가 조금씩 달라집니다.

천왕성의 실제 위치와 천문표가 어긋나다

윌리엄 허셜이 1781년 천왕성을 행성으로 확인한 뒤 천문학자들은 과거 관측 기록까지 모아 궤도를 계산했습니다. 천왕성은 맨눈으로도 조건이 좋으면 희미하게 보이기 때문에 행성으로 인정되기 전부터 별로 기록된 관측값이 있었습니다. 긴 기간의 자료는 궤도를 정하는 데 도움을 줬지만 서로 다른 시기의 위치를 하나의 궤도표로 맞추는 과정에서 오차가 드러났습니다.

알렉시 부바르는 1821년 목성·토성·천왕성의 천문표를 출판했습니다. 목성과 토성이 천왕성에 주는 섭동을 포함해 앞으로의 위치를 계산했지만, 시간이 지나면서 실제 관측과 표의 예측 사이에 체계적인 차이가 나타났습니다. 단순한 한 번의 측정 실수라면 오차 방향이 들쭉날쭉해야 하지만 잔차에는 시간에 따른 흐름이 있었습니다.

잔차 Δλ(t) = 관측 경도 λ_obs(t) − 계산 경도 λ_calc(t)

여기서 경도는 천구나 궤도면에서 행성의 방향을 나타내는 각도입니다. 르베리에는 관측값 자체보다 관측과 계산의 차이인 잔차를 분석했습니다. 미지 행성의 신호는 천왕성의 전체 공전 운동보다 훨씬 작았기 때문에 큰 기준 운동을 먼저 제거하고 남은 작은 패턴을 살펴야 했습니다.

오차가 생겼다고 바로 새 행성을 뜻하지는 않는다

예측과 관측이 어긋나는 원인은 하나가 아닙니다. 관측 기기의 오차, 별 위치표의 오차, 천왕성 자체의 궤도요소 오류, 알려진 행성의 질량 오차, 계산에서 빠진 중력 효과 또는 중력 법칙의 문제를 검토해야 합니다. 르베리에는 알려진 행성의 섭동과 천왕성 궤도부터 다시 계산해 기존 이론만으로 잔차를 없앨 수 있는지 확인했습니다.

그 과정을 거친 뒤에도 남는 체계적인 오차를 천왕성 바깥의 미지 행성이 만든다고 가정했습니다. 이 가정 아래에서 문제는 ‘주어진 행성으로 천왕성의 흔들림을 구하라’는 순방향 계산이 아니라 ‘관측된 흔들림을 만드는 행성의 궤도를 찾아라’는 역방향 계산이 됩니다.

순방향 문제와 역문제의 차이

구분	알고 있는 값	구하려는 값
순방향 궤도 계산	행성의 질량·위치·속도	미래 위치와 섭동
해왕성 역산 문제	천왕성의 관측 잔차	미지 행성의 질량·궤도·현재 위치

역문제는 여러 원인 조합이 비슷한 결과를 만들 수 있어 어렵습니다. 멀리 있는 큰 행성과 조금 가까운 작은 행성이 특정 기간에는 비슷한 섭동을 줄 수 있습니다. 관측 기간도 미지 행성의 긴 공전 주기에 비해 짧았으므로 가능한 궤도를 하나로 좁히려면 가정과 근사가 필요했습니다.

미지 행성의 중력이 천왕성에 만드는 가속도

천왕성 위치벡터를 r_U, 미지 행성 위치벡터를 r_N, 미지 행성 질량을 m_N이라고 하겠습니다. 미지 행성이 천왕성에 직접 만드는 가속도는 두 행성 사이 거리의 제곱에 반비례합니다. 태양도 그 행성에 끌리므로 태양 중심 좌표에서는 태양이 받는 가속도를 함께 빼야 합니다. 현대 벡터식으로는 다음처럼 나타낼 수 있습니다.

δa_U = Gm_N[(r_N−r_U)/|r_N−r_U|³ − r_N/|r_N|³]

첫 항은 미지 행성이 천왕성을 끄는 효과이고 두 번째 항은 같은 행성이 태양을 끄는 효과입니다. 중요한 것은 천왕성과 태양이 받는 가속도의 차이입니다. 이 작은 차이가 오랜 시간 누적되면 천왕성의 경도와 태양까지의 거리에 관측 가능한 변화가 생깁니다.

미지 행성을 찾기 위해 조정해야 했던 값

행성 궤도는 현재 하늘 방향 하나로 정해지지 않습니다. 궤도의 크기, 찌그러진 정도, 궤도면의 방향과 특정 시점의 위치가 필요합니다. 여기에 미지 행성의 질량도 천왕성에 주는 섭동의 크기를 바꿉니다.

궤도 긴반지름 또는 평균 거리
궤도 이심률
근일점의 방향
기준 시점의 평균 경도
궤도면과 관련된 요소
미지 행성의 질량
천왕성 기존 궤도요소의 보정값

르베리에는 먼저 가능한 평균 거리에 가정을 두고 섭동식을 전개한 다음 관측 잔차와 가장 잘 맞도록 미지수들을 조정했습니다. 후속 계산에서는 미지 행성의 평균 거리도 풀어야 할 값으로 다뤘습니다. 계산은 오늘날처럼 컴퓨터로 수치 적분한 것이 아니라 긴 급수 전개와 연립 조건방정식을 손으로 처리한 작업이었습니다.

조건방정식으로 관측 잔차를 맞추는 원리

각 관측 시점 t_i에서 계산한 잔차 Δλ_i는 미지 행성의 궤도요소와 질량에 따라 달라집니다. 기준으로 잡은 가상 궤도 주변에서 변화가 작다고 보고 1차 항까지 전개하면 여러 미지수에 관한 선형에 가까운 식을 만들 수 있습니다.

Δλ_i ≈ c_i1δp₁ + c_i2δp₂ + … + c_inδpₙ

δp₁, δp₂ 등은 질량이나 궤도요소의 보정량이고 c_ij는 각 보정이 해당 시점의 천왕성 위치에 얼마나 영향을 주는지 계산한 계수입니다. 관측값이 미지수보다 많으면 모든 식을 정확히 만족시키기보다 잔차 제곱합을 작게 만드는 해를 찾습니다. 르베리에 시대의 실제 전개는 훨씬 복잡했지만 기본 구조는 관측과 모델의 차이를 가장 잘 줄이는 매개변수를 찾는 것입니다.

보이지 않는 원인을 데이터로 추정한다

이 문제는 직접 보지 못한 대상을 관측 가능한 흔적으로 추정한다는 점에서 독일 탱크 문제에서 일부 일련번호로 전체 생산량을 추정하는 방법과 구조가 비슷합니다. 탱크 문제에서는 표본의 최댓값이 전체 수량에 관한 정보를 주고, 해왕성 문제에서는 천왕성의 시간별 위치 잔차가 미지 행성의 중력에 관한 정보를 줍니다.

두 문제 모두 계산 전에 표본이나 관측값이 어떤 과정으로 만들어졌는지 모델을 세워야 합니다. 모델이 틀리거나 관측에 체계적인 편향이 있으면 정밀한 계산도 잘못된 결론을 낼 수 있습니다. 해왕성 사례에서는 뉴턴 역학이 태양계 행성 운동에 적용된다는 가정이 계산의 기반이었습니다.

르베리에의 세 단계 연구

르베리에는 1845년과 1846년에 걸쳐 천왕성 문제를 단계적으로 발표했습니다. 먼저 알려진 행성들이 천왕성에 미치는 섭동을 다시 계산해 기존 이론의 오류만으로 관측 차이를 설명할 수 있는지 살폈습니다. 다음으로 천왕성 바깥의 행성을 가정해 가능한 궤도와 위치 범위를 구했고, 마지막에는 관측에 필요한 구체적인 하늘 위치와 물리적 특성까지 제시했습니다.

1846년 8월 31일 발표한 결과에는 1847년 1월 1일의 예측 경도 326°32′가 포함됐습니다. 르베리에는 새 행성이 염소자리 델타별에서 동쪽으로 약 5도 떨어진 부근에 있을 것이라고 안내했습니다. 망원경으로 무작정 황도 전체를 훑는 대신 좁은 영역을 집중적으로 비교할 수 있는 좌표였습니다.

계산이 실제 관측으로 이어진 밤

르베리에는 1846년 9월 18일 베를린 천문대의 갈레에게 계산 결과를 보내 관측을 요청했습니다. 편지는 9월 23일 도착했습니다. 그날 밤 갈레와 다레스트는 예측 영역을 망원경으로 보고 새로 준비된 별 지도와 비교했습니다. 지도에 없는 천체를 찾은 뒤 다음 밤 위치가 움직인 것을 확인해 행성임을 검증했습니다.

NASA의 역사 자료는 새 행성이 계산 위치에서 약 1도 떨어진 곳에서 발견됐다고 설명합니다. 같은 기준 시점과 경도 체계로 비교한 역사 연구에서는 예측 경도 326°32′와 관측으로 정한 327°24′의 차이를 52분각으로 제시합니다. 보름달의 겉보기 지름이 약 0.5도라는 점을 생각하면 하늘에서 매우 좁은 범위였습니다.

327°24′ − 326°32′ = 0°52′

애덤스도 독립적으로 계산했다

영국의 수학자 존 쿠치 애덤스도 천왕성의 불규칙한 운동을 미지 행성의 섭동으로 설명하는 계산을 독립적으로 진행했습니다. 그는 1840년대 중반에 예측 결과를 영국 천문학계에 전달했지만 제때 출판하거나 실제 발견으로 연결하지 못했습니다. 케임브리지 천문대에서는 탐색이 이뤄졌으나 후보를 즉시 행성으로 식별하지 못했습니다.

발견 뒤 영국과 프랑스에서는 예측 공로를 둘러싼 논쟁이 벌어졌습니다. 애덤스의 독립 연구는 수학사에서 중요하지만 베를린 천문대가 실제로 사용한 좌표는 르베리에가 공개하고 갈레에게 보낸 예측이었습니다. 따라서 독립 계산의 공헌과 관측 발견으로 직접 이어진 경로를 나누어 설명하는 것이 정확합니다.

위치는 맞았지만 궤도 전체는 정확하지 않았다

해왕성 발견을 ‘계산만으로 궤도를 완벽히 알아냈다’고 표현하면 사실과 다릅니다. 르베리에가 예측한 긴반지름은 약 36.15 AU였지만 현대값은 약 30.06 AU입니다. 이에 따라 예측 공전 주기는 약 217년이었고 실제 주기는 약 165년입니다. 예측 이심률도 실제보다 컸습니다.

항목	르베리에 예측	현대에 알려진 값	평가
1846년 무렵 하늘 위치	실제 위치에서 1도 이내	관측으로 확인	탐색에 충분히 정확
궤도 긴반지름	약 36.15 AU	약 30.06 AU	크게 예측
공전 주기	약 217년	약 165년	길게 예측

관측된 천왕성의 섭동은 제한된 기간의 자료였습니다. 서로 다른 미지 행성 궤도가 그 기간에는 비슷한 잔차를 만들 수 있었고, 실제 해왕성이 1846년 당시 예측 궤도와 비슷한 하늘 구간에 우연히 가까이 놓여 있었습니다. 그래서 현재 위치는 잘 맞았지만 장기 궤도요소에는 큰 오차가 남았습니다.

AU와 각도 단위를 함께 이해하기

AU는 천문단위로 지구와 태양 사이 평균 거리를 기준으로 한 길이 단위입니다. 해왕성의 궤도 긴반지름 30.06 AU는 태양에서 평균적으로 지구보다 약 30배 멀다는 뜻입니다. 반면 326°32′ 같은 값은 거리가 아니라 황도면에서 본 방향을 나타내는 경도입니다.

1도는 60분각이고 1분각은 60초각입니다. 예측과 관측의 차이 52분각은 약 0.867도입니다. 거리 오차와 방향 오차는 서로 다른 양이므로 ‘1도 이내에서 발견됐다’는 말이 태양까지의 거리를 1% 이내로 맞혔다는 뜻은 아닙니다.

현대에는 같은 문제를 어떻게 계산할까

현대 궤도결정에서는 여러 천체의 운동방정식을 컴퓨터로 수치 적분합니다. 후보 행성의 초기 위치, 속도와 질량을 넣어 관측 시점의 천왕성 위치를 계산하고 실제 자료와의 잔차를 구합니다. 그런 다음 최소제곱법이나 베이즈 추정 같은 방법으로 매개변수를 반복 조정합니다.

χ²(p) = Σ[(관측값_i−모델값_i(p))/σ_i]²

σ_i는 각 관측의 불확도이고 p는 미지 천체의 궤도와 질량을 모은 매개변수입니다. χ²가 작아지는 방향으로 p를 바꾸되 여러 해가 가능한지, 자료에 민감한지, 관측 오차를 올바르게 반영했는지도 검사합니다. 르베리에의 손계산과 현대 수치해석은 도구가 다르지만 관측 잔차를 줄이는 원인을 찾는 기본 목표는 같습니다.

간접 측정이 성공하려면 필요한 조건

보이지 않는 천체를 잔차로 찾으려면 기준 모델이 충분히 정확해야 합니다. 알려진 행성의 질량이나 궤도가 틀리면 그 오차가 미지 행성의 신호처럼 보일 수 있습니다. 관측 시계, 별 목록과 좌표 기준도 일관돼야 하며 오래된 관측과 새 관측의 정확도 차이도 고려해야 합니다.

직접 측정할 수 없는 대상을 각도 관계로 추정한다는 점에서는 에라토스테네스가 두 지역의 그림자 각도로 지구 둘레를 계산한 방법과 닮았습니다. 다만 지구 둘레 문제는 단순한 기하 비례식에 가깝고 해왕성 문제는 시간에 따라 누적되는 여러 천체의 중력 효과를 동시에 풀어야 하는 역문제입니다.

수성의 오차에서는 왜 새 행성을 찾지 못했을까

해왕성 발견의 성공 뒤 르베리에는 수성의 근일점 이동에 남는 오차도 분석했습니다. 당시에는 수성 안쪽의 미지 행성이나 물질이 원인일 수 있다고 보았고 가상의 행성에 벌컨이라는 이름까지 붙었습니다. 그러나 그런 행성은 발견되지 않았습니다.

수성의 잔여 근일점 이동은 훗날 아인슈타인의 일반상대성이론으로 설명됐습니다. 관측 잔차가 있다고 해서 언제나 보이지 않는 천체가 존재하는 것은 아니라는 사례입니다. 해왕성 문제에서는 미지 행성 가설이 맞았지만 다른 문제에서는 기본 이론을 수정해야 할 수 있습니다.

해왕성 섭동 계산에서 자주 생기는 오류

천왕성이 뉴턴 법칙을 어겨서 해왕성을 찾았다고 설명합니다. 실제로는 여러 천체에 뉴턴 중력을 적용했을 때 생기는 추가 섭동을 찾았습니다.
르베리에가 해왕성을 망원경으로 직접 발견했다고 말합니다. 실제 관측자는 베를린 천문대의 갈레와 다레스트였습니다.
애덤스의 독립 계산을 생략하거나 반대로 그의 좌표가 베를린 발견에 사용됐다고 설명합니다.
예측 위치가 1도 이내였다는 사실을 궤도 전체와 질량까지 정확했다는 뜻으로 확대합니다.
경도 차이와 태양까지의 거리 차이를 같은 단위의 오차처럼 비교합니다.
한 번의 관측 오차만으로 미지 행성을 추정했다고 말합니다. 실제로는 오랜 기간 누적된 잔차 패턴을 사용했습니다.
모든 관측 잔차는 반드시 새 천체 때문이라고 가정합니다. 측정·모델·상수의 오류도 검토해야 합니다.

해왕성 발견 계산 FAQ

르베리에는 해왕성을 보지 않고 위치를 어떻게 알았나요?

알려진 행성들의 중력을 반영해 천왕성 위치를 계산한 뒤 실제 관측과의 잔차를 구했습니다. 천왕성 바깥의 미지 행성을 가정하고 그 질량과 궤도요소를 바꿔 가며 잔차 패턴을 가장 잘 재현하는 해를 찾았습니다.

해왕성은 계산 위치와 정확히 같은 곳에 있었나요?

완전히 같지는 않았습니다. 관측에 필요한 하늘 위치는 약 1도 이내로 맞았지만 긴반지름과 공전 주기 같은 전체 궤도요소에는 상당한 차이가 있었습니다. 계산은 망원경의 탐색 범위를 좁히기에 충분히 정확했습니다.

애덤스와 르베리에 중 누가 먼저 해왕성을 발견했나요?

두 사람은 미지 행성의 위치를 독립적으로 계산했습니다. 르베리에의 결과가 먼저 공개됐고 그의 좌표가 베를린 관측으로 직접 이어졌습니다. 실제 망원경 관측으로 해왕성을 확인한 사람은 갈레와 다레스트입니다. 계산의 우선권과 관측 발견을 구분해야 합니다.

갈릴레오가 해왕성을 먼저 봤다는 말은 사실인가요?

갈릴레오는 1612년과 1613년 목성 주변을 관측한 기록에 해왕성을 별처럼 남겼습니다. 그러나 움직이는 행성으로 확인하지 않았기 때문에 해왕성의 발견으로 분류되지는 않습니다. 1846년 관측은 계산으로 예측된 천체의 이동을 확인해 행성임을 알아본 사건입니다.

천왕성의 오차에서 해왕성까지 핵심 정리

천문학자들은 천왕성의 관측 위치에서 알려진 행성과 뉴턴 역학으로 계산한 위치를 빼 잔차를 만들었습니다. 르베리에는 알려진 섭동과 천왕성 궤도요소를 다시 점검한 뒤에도 남는 잔차를 바깥쪽 미지 행성의 중력으로 설명했습니다. 후보 행성의 질량과 궤도요소를 조정해 오랜 기간의 잔차를 줄이고 1846년 하늘에서 찾아야 할 방향을 계산했습니다.

갈레와 다레스트는 1846년 9월 23일 그 예측을 이용해 해왕성을 약 1도 이내에서 확인했습니다. 애덤스도 독립 계산을 수행했지만 베를린 발견은 르베리에의 공개 예측에서 직접 이어졌습니다. 예측 위치는 탁월했으나 예측 궤도의 크기와 주기는 실제와 달랐습니다. 해왕성 발견은 정확한 물리 법칙, 긴 관측 기록, 잔차 분석과 망원경 검증이 결합된 과학적 역추론의 사례입니다.