Impact of ab initio E1 strength on neutron capture rates along the nickel isotopic chain

별들의 심장부 깊은 곳에서는 우주적 연금술이 끊임없이 작용하며, 우리 세계를 구성하는 무거운 원소들을 벼려내고 있습니다. 탄소나 산소와 같은 가벼운 원소들이 항성의 안정적인 연소 과정에서 탄생하는 반면, 금이나 백금과 같은 가장 무거운 원소들은 중성자가 엄청난 속도로 포획되는 격렬하고 극단적인 사건 속에서 만들어집니다. 'r-과정(rapid neutron-capture process)'이라 알려진 이 과정은 이러한 무거운 물질들을 만들어내는 우주의 주요 공장입니다. 각 원소가 정확히 얼마나 생성되는지 이해하기 위해, 과학자들은 원자핵이 중성자를 포획하는 정밀한 비율을 알아야 합니다. 그러나 이 우주적 춤에 참여하는 불안정한 중성자 과잉 핵 중 다수는 단 1초 미만의 짧은 시간 동안만 존재하기 때문에, 실험실에서 직접 연구하는 것이 불가능합니다. 결과적으로, 연구자들은 직접적인 측정이 불가능한 빈틈을 메우기 위해 이 찰나의 원자들이 어떻게 행동할지 예측하는 이론적 계산에 의존해야 합니다.

이 퍼즐의 핵심적인 조각은 이 핵들이 빛과 어떻게 상호작용하는지, 구체적으로 '전기 쌍극자 강도(electric dipole strength)'라고 불리는 일종의 진동을 통해 어떻게 에너지를 흡수하는지에 달려 있습니다. 간단히 말해, 이는 에너지를 받았을 때 원자 내부의 양성자와 중성자가 서로에 대해 어떻게 이동하는지를 설명합니다. 수십 년 동안 과학자들은 이러한 불안정한 원자들의 행동을 추측하기 위해 광범위하고 일반화된 모델을 사용해 왔으나, 이러한 추측은 우주에서의 원소 생성을 정확하게 시뮬레이션하는 데 필요한 정밀도가 부족한 경우가 많았습니다. 프란체스카 보나이티(Francesca Bonaiti)와 조안나 E. 솝직(Joanna E. Sobczyk)의 새로운 연구는 이러한 도전에 정면으로 맞서며, 단순한 추측에 의존하는 대신 입자 간의 근본적인 힘으로부터 계산을 시작하는 방식을 사용하여 이러한 상호작용을 밑바닥부터 계산합니다.

연구진은 니켈 동위원소라는 특정 원자 가계에 집중했습니다. 니켈은 r-과정에서 핵심적인 역할을 수행하며, 질량수가 100 근처인 원소들의 풍부함에 영향을 미치는 디딤돌 역할을 합니다. 연구팀은 '아비니토(ab initio)', 즉 '처음부터'라는 의미의 강력한 계산적 접근법을 사용했습니다. 기존 데이터에 계산을 맞추는 대신, 그들은 '카이랄 유효장 이론(chiral effective field theory)'에서 유도된 양성자와 중성자 사이의 상호작용, 즉 핵 물리학의 근본 법칙으로부터 시작했습니다. 그들은 이를 안정적인 니켈-58부터 극도로 중성자가 많은 니켈-80에 이르는 니켈 원자 사슬에 적용했습니다. 복잡한 수학적 구조를 가진 이 불안정한 열린 껍질(open-shell) 핵들을 처리하기 위해, 연구진은 '로렌츠 적분 변환 결합 클러스터(Lorentz integral transform coupled-cluster)' 기법의 새로운 확장을 개발했습니다. 이를 통해 계산을 어렵게 만드는 산란 상태(scattering states)를 직접 구축하지 않고도 이 원자들이 에너지에 어떻게 반응하는지 시뮬레이션할 수 있었습니다.

결과는 중성자가 핵에서 탈출할 수 있는 지점 바로 위에서 나타나는 독특하고 낮은 수준의 강도라는 중요한 특징을 드러냈습니다. 흔히 '피그미 쌍극자 공명(pygmy dipole resonance)'과 연관되는 이 현상은, 핵이 더 많은 중성자를 얻을수록 더욱 뚜렷해지는 부드럽고 낮은 에너지의 진동처럼 작용합니다. 니켈-80의 경우, 이 저에너지 진동은 매우 강력하고 탈출 임계값에 매우 가깝게 위치하여 약 0.50억 도라는 놀라울 정도로 낮은 온도에서도 원자의 행동에 영향을 미치기 시작합니다. 다른 동위원소들의 경우, 이 효과는 중성자 수가 증가함에 따라 더 강해지며, 30억 도에서 100억 도 사이의 온도에서 주요한 요인이 됩니다.

연구진이 이 새로운 제일원리(first-principles) 계산 결과를 별의 폭발을 시뮬레이션하는 데 사용되는 표준 컴퓨터 코드에 입력했을 때, 그 영향은 명확했습니다. 이 특정 전기 쌍극자 강도를 포함함으로써 니켈 핵이 중성자를 포획하는 비율이 크게 높아졌습니다. 더 무거운 동위원소들의 경우, 온도와 특정 동위원소에 따라 포획률이 이전 추정치보다 16%에서 최대 3배까지 증가했습니다. 이러한 향상은 사소한 조정이 아니라, 이 핵들이 환경과 어떻게 상호작용하는지에 대한 우리의 이해를 근본적으로 변화시키는 것입니다. 또한 연구는 결과의 불확실성을 초래하는 가장 큰 원인이 물리 모델 자체보다는 계산을 최종 답을 얻기 위해 역산하는 수학적 과정에 있음을 강조했습니다.

이러한 전기 쌍극자 강도에 대한 최초의 엄격한 제일원리 기술을 제공함으로써, 이 연구는 천체 물리학 시뮬레이션을 위한 더 신뢰할 수 있는 토대를 마련했습니다. 이는 우주에서 무거운 원소들이 생성되는 과정이 특히 가장 무거운 니켈 동위원소들에 대해 이전에 생각했던 것보다 더 효율적일 수 있음을 시사합니다. 연구진은 불확실성을 줄이기 위해 수학적 역산 과정을 정교화하는 후속 작업이 필요함을 인정하면서도, 본 연구가 기초 핵 이론과 거대한 규모의 우주적 원소 형성 사이의 간극을 성공적으로 메웠음을 보여주었습니다. 이는 가장 불안정한 원자들에 대해서도 이제 우리가 이전에는 도달할 수 없었던 수준의 세밀함으로 그 행동을 예측할 수 있게 되었으며, 우리 우주의 무거운 원소들이 진정 어디에서 왔는지에 대한 미스터리를 풀기에 한 걸음 더 다가갔음을 증명합니다.