열적 구슬에서 Q볼까지 첫 번째 차수 전이 중 수치 연구

초록

본 논문은 Friedberg‑Lee‑Sirlin 모델을 이용해 우주 초기의 1차 강상전이 동안 Q볼이 어떻게 형성되는지를 3차원 격자 시뮬레이션으로 직접 확인한다. 붕괴하는 거짓 진공 영역이 먼저 ‘열적 구슬’(thermal ball)을 만들고, 마찰과 입자 소멸을 통해 냉각되면서 전하를 보존하는 회전 Q볼로 전이한다. 시뮬레이션 결과는 저질량 Q볼이 풍부히 생성되고, 드물게 분석적 예측보다 수배 무거운 Q볼이 나타나는 넓은 질량 스펙트럼을 보여준다. 전체 Q볼 밀도는 기존 분석보다 약 50 % 높아, 다크 물질 후보로서의 가능성을 재평가한다.

상세 분석

본 연구는 Friedberg‑Lee‑Sirlin(FLS) 모델의 복소 스칼라 χ와 실수 스칼라 ϕ 두 장을 기반으로 한다. χ는 전역 U(1) 대칭을 가지고 있어 전하 Q가 보존되며, ϕ의 진공 기대값에 따라 χ의 유효 질량 mχ(ϕ)=gϕ가 결정된다. ϕ가 0에 가까운 거짓 진공 영역에서는 χ가 거의 무질량이 되므로, χ 입자들은 ϕ가 큰 진짜 진공으로 이동하려는 경향이 있다. 이 메커니즘이 바로 ‘전하 포획(trapping)’ 효과이며, 거짓 진공 구역이 축소될 때 χ 입자들이 중심에 모여 고밀도 전하 구역을 만든다.

시뮬레이션은 CosmoLattice를 이용해 512³ 격자에 Δx와 Δt를 적절히 선택하고, 초기에는 N_b개의 임계 버블을 무작위 위치에 삽입한다. 버블 사이의 겹침을 방지하기 위해 중심 간 거리를 2R_c 이상으로 두고, 초기 전하 Q는 χ의 운동량 πχ에 비대칭을 부여해 구현한다. 마찰 계수 Γ는 0으로 두어 플라즈마 마찰을 무시했지만, ϕ 방정식에 포함된 ‘γ·∂tϕ’ 항은 실제 물리적 마찰을 모사한다.

시뮬레이션 결과는 크게 두 단계로 나뉜다. 첫 번째는 거짓 진공 구역이 급격히 수축하면서 ‘열적 구슬’이 형성되는 단계이다. 이때 χ 입자들은 양·음 전하가 혼재된 상태로, 열적 압력에 의해 구슬이 일시적으로 팽창한다. 두 번째는 마찰에 의해 입자들의 운동 에너지가 소멸되고, 양전하와 음전하가 서로 소멸하면서 전하가 순수하게 남는다. 이 과정에서 구슬은 점차 냉각되고, 최종적으로 전하를 보존하는 안정적인 Q볼이 된다. 특히, 구형이 아닌 ‘덤벨’ 형태를 보이는 Q볼은 각운동량이 크게 축적된 회전 Q볼이며, 이는 버블 충돌과 전하 흐름이 비대칭적으로 진행된 결과로 해석된다.

질량 스펙트럼 분석에서는 Q볼 질량 M을 무차원 변수 ˜M=M ω* η²로 정규화하고, t=1000 ω*⁻¹ 시점에 3,352개의 Q볼을 식별했다. 질량 분포는 낮은 질량 영역에서 급격히 상승하고, ˜M≈200 이상에서는 P(˜M)∝exp(−B ˜M) 형태의 지수 꼬리를 보인다. 여기서 A≈−8.21, B≈3.85×10⁻³이 최적화된 파라미터이다. 분석적 예측값 ˜M*≈173.8과 비교했을 때, 평균 질량은 비슷하지만, 저질량 Q볼이 전체 수의 54 %를 차지하고, 고질량 꼬리 역시 전체 질량의 약 10 %를 기여한다는 점에서 단색(단일 질량) 가정이 크게 틀렸음을 보여준다.

우주론적 함의는 두드러진다. 시뮬레이션에서 측정한 Q볼 에너지 밀도 ρ_num은 분석식 ρ_ana에 비해 r=ρ_num/ρ_ana≈1.5±0.24의 보정 인자를 요구한다. 이는 전하 비대칭 Δχ≈7×10⁻¹¹(전기 약 100 GeV 전이) 혹은 Δχ≈1.5×10⁻¹⁴(10 GeV 전이) 정도면 관측된 다크 물질 밀도를 설명할 수 있음을 의미한다. 또한, 질량이 충분히 큰 Q볼은 Q∝M^{3/2} 관계에 따라 중력 붕괴가 가능해, 원시 블랙홀(PBH) 형성 메커니즘과 연결될 여지가 있다.

기술적 한계로는 Γ=0 설정, 유한한 격자 크기, 그리고 통계적으로 희귀한 초대질량 Q볼의 샘플 부족이 있다. 향후 연구에서는 플라즈마 마찰을 포함한 전산 유체 역학, 더 큰 격자와 높은 N_b, 그리고 다양한 파라미터 스캔을 통해 Q볼 형성 효율과 스펙트럼의 민감도 분석이 필요하다.