How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

우주를 모든 방향으로 뻗어 있는, 주로 수소로 이루어진 거대하고 보이지 않는 가스의 바다라고 상상해 보세요. 빅뱅 이후 처음 몇 억 년 동안, 이 바다는 어둡고 차가우며 고요했습니다. 그러다 첫 번째 별과 은하들이 안개 속의 등대처럼 불을 밝히기 시작했습니다. 그들의 강렬한 복사 에너지는 가스를 가열하고 원자로부터 전자를 떼어내기 시작했는데, 이 과정을 '재이온화'라고 부릅니다. 천문학자들은 이 시대가 어떻게 펼쳐지는지 관찰하기 위해 필사적입니다. 왜냐하면 이 시대는 우리 우주가 어떻게 성장했는지에 대한 비밀을 간직하고 있기 때문입니다. 이를 위해 그들은 거대한 귀 역할을 하며 희미한 속삭임인 '21cm 신호'를 찾아낼 거대한 라디오 망원경인 스퀘어 킬로미터 어레이(SKA)를 구축하고 있습니다. 이 신호는 중성 수소가 방출하는 특정한 라디오파이며, 그 강도는 가스가 얼마나 뜨거운지에 따라 달라집니다. 가스가 차가우면 신호가 따뜻할 때와는 다르게 나타납니다. 이 신호를 측정함으로써 과학자들은 우주의 가열 역사를 지도화하기를 희망합니다.

하지만 문제가 있습니다. 우주는 너무 광대하고 복잡해서 실시간으로 직접 관찰할 수 없기에, 우리는 수십억 년 전에 일어난 일을 시뮬레이션하기 위해 컴퓨터 모델을 만들어야 합니다. 이 시뮬레이션들은 디지털 타임머신과 같습니다. 하지만 어떤 타임머신과 마찬가지로, 이 모델들도 따르는 규칙이 얼마나 정교한가에 따라 성능이 결정됩니다. 어떤 시뮬레이션은 매우 상세하여 모든 가스 입자와 빛의 줄기를 추적하지만, 이를 실행하는 데 슈퍼컴퓨터로 몇 주가 걸립니다. 반면 다른 것들은 더 빠르고 단순하며, 지름길을 사용하여 사용하기는 쉽지만 잠재적으로 정확도가 떨어질 수 있습니다. 큰 질문은 이것입니다. 과연 이 지름길들이 결과를 신뢰할 수 있을 만큼 충분히 잘 작동하는가? 만약 두 개의 서로 다른 시뮬레이션 코드가 가스가 얼마나 뜨거워졌는지에 대해 서로 다른 이야기를 한다면, 우리는 어느 것을 믿어야 할까요? 로맹 메리오(Romain Meriot)와 그의 팀이 수행한 새로운 연구가 해결하고자 하는 퍼즐이 바로 이것입니다.

연구팀은 이 두 가지 매우 다른 '타임머신'을 맞붙여, 이들이 X선에 의한 우주 가스의 가열을 어떻게 처리하는지 확인하기로 했습니다. 한쪽에는 고해상도 3D 시뮬레이션인 **리코리스(Licorice)**가 있었습니다. 이는 마치 풀 모션 비디오 게임처럼 가스가 어떻게 움직이고 빛이 어떻게 이동하는지를 3차원으로 계산합니다. 다른 한쪽에는 더 빠르고 단순한, '1차원적' 접근 방식을 사용하는 **베오른(Beorn)**이라는 코드가 있었습니다. 리코리스를 모든 재료를 개별적으로 저으며 복잡한 스튜를 요리하는 마스터 셰프로 본다면, 베오른은 평균적인 레시피를 바탕으로 전체 음식을 데우는 전자레인지와 같습니다. 연구진은 전자레인지가 셰프의 스튜만큼 맛있는 음식을 만들어낼 수 있는지, 특히 X선이 수소 가스를 데우는 방식에 있어서 말입니다.

공정한 대결을 위해, 과학자들은 단순히 두 코드를 각자의 무작위 설정으로 실행하게 두지 않았습니다. 대신, 그들은 상세한 리코리스 시뮬레이션에서 가져온 '재료'(별과 그 빛의 출력)를 베오른에 직접 입력했습니다. 이를 통해 두 코드 모두 서로 다른 방법을 사용하되, 정확히 똑같은 우주를 요리하도록 만들었습니다. 그런 다음 그들은 가스의 온도, 21cm 신호의 밝기, 그리고 그 신호가 하늘에서 어떻게 변동하는지에 대한 통계적 패턴을 비교했습니다.

결과는 희소식과 경고가 섞여 있었습니다. 거시적인 관점, 즉 가스의 평균 온도와 신호의 전체적인 강도를 보았을 때, 두 코드는 매우 잘 일치했습니다. 그것은 마치 전자레인지와 셰프가 모두 스튜를 적절한 평균 온도로 맞춘 것과 같았습니다. 그러나 과학자들이 세부 사항을 자세히 들여다보자 차이점이 나타나기 시작했습니다. 열의 분포가 같지 않았습니다. 상세한 3D 시뮬레이션에서는 열이 얼룩덜룩하여 어떤 곳은 훨씬 뜨겁고 어떤 곳은 더 차가운 복잡한 질감을 형성했습니다. 반면 더 빠른 1D 시뮬레이션에서는 열이 더 매끄럽고 균일했습니다.

열의 '질감'에 나타난 이러한 작은 차이는 놀랍게도 중요한 것으로 드러났습니다. 팀이 신호가 다양한 규모에서 얼마나 변동하는지를 측정하는 방식인 '파워 스펙트럼(power spectrum)'을 계산했을 때, 두 코드 사이에 약 **30%**의 차이가 있음을 발견했습니다. 이를 설명하자면, 케이크 한 입을 먹고 레시피를 추측하려 할 때, 맛의 30% 차이는 재료를 완전히 잘못 추측하게 된다는 것을 의미합니다.

연구진은 왜 코드가 서로 다른지 알아내기 위해 더 깊이 파고들었습니다. 그들은 빠른 코드에서 사용된 지름길이 범인이라는 것을 발견했습니다. 구체적으로, 1D 코드는 가스의 밀도가 어디서나 동일하다고 가정했으며(마치 바다가 완벽하게 평평하다고 가정하는 것처럼), X선이 실제로 얼마나 많은 열을 전달하는지에 대해 단순화된 공식을 사용했습니다. 실제로는 가스가 덩어리져 있으며, 가열은 가스가 정확히 어디에 있는지에 따라 달라집니다. 연구팀이 상세한 3D 코드가 이와 동일한 단순화된 지름길을 사용하도록 강제하자, 결과가 1D 코드와 더 가까워졌으며, 이는 이러한 근사치들이 오류의 원인이었음을 확인시켜 주었습니다.

마지막으로, 팀은 궁극적인 질문을 던졌습니다. "만약 우리가 미래의 SKA 망원경 데이터를 분석하기 위해 이 서로 다른 시뮬레이션들을 사용한다면, 틀린 답을 얻게 될 것인가?" 그들은 망원경이 100시간 동안 관측하는 시나리오를 시뮬레이션했습니다. 그들은 두 코드 사이의 차이가 첫 번째 별들의 특성을 추정하는 데 있어 상당한 편향(bias)을 초래할 것이라는 점을 발견했습니다. 통계적으로 이 오차는 1 시그마(sigma) 이상의 신뢰 수준을 가졌으며, 일부 매개변수의 경우 그보다 더 높았습니다. 이는 만약 천문학자들이 실제 데이터를 해석하는 데 더 빠르고 단순한 코드를 사용한다면, 첫 번째 별들이 실제와는 달랐다고 확신하며 잘못된 결론을 내릴 수도 있음을 의미합니다.

이 연구는 빠른 코드들이 유용하고 '평균적인' 이야기는 잘 전달하지만, SKA 시대에 쓰이기에는 아직 정밀도가 부족하다고 결론짓습니다. 파워 스펙트럼의 30% 차이와 그로 인한 데이터의 편향은, 우리가 극도로 정밀하게 우주를 측정하려 할 때 단순히 지름길을 신뢰할 수 없음을 시사합니다. 저자들은 현재의 덜 민감한 망원경들에게는 이 코드들이 충분히 유사할 수 있지만, SKA가 가동되기 전에 빠른 코드들을 개선하거나 오류를 수정할 새로운 방법을 개발하는 데 많은 노력이 필요하다고 제언합니다. 이는 우주의 여명을 이해하려는 탐구에서, 디테일이 정말 중요하다는 사실을 다시 한번 일깨워줍니다.