Can photoevaporation open gaps in protoplanetary discs?

문제 제기
원시 행성계 원반의 분산은 행성 형성 시기와 조건을 결정하는 결정적인 요인이다. 고에너지 항성 복사에 의해 유발되는 열적 바람을 통해 가스가 제거되는 현상인 광증발(photoevaporation)이 핵심적인 분산 기제로서 널리 인정받고 있으나, 이것이 어떻게 간극(gap)을 열고 유지하는지에 대한 이론적 근거는 여전히 불확실하다. 표준 모델들은 대개 광증발 질량 손실률 $\dot{\Sigma}(r)$을 정적인 2차원(2D) 시뮬레이션에서 유도된 고정된 싱크(sink) 항으로 취급하는 1차원(1D) 점성 진화 계산에 의존한다. 이러한 접근 방식은 원반 구조가 진화함에 따라 질량 손실 프로파일이 일정하게 유지된다고 암묵적으로 가정한다. 저자들은 이것이 중요한 피드백 루프를 간과하고 있다고 주장한다. 즉, 간극이 형성되기 시작하면 국지적인 원반 구조가 변하여 잠재적으로 바람의 방출 조건을 변화시키고, 국지적 질량 손실률을 감소시키며, 점성 유입(viscous inflow)이 결핍된 영역을 다시 채우도록 허용한다는 것이다. 만약 광증발이 이와 같이 자기 제한적(self-limiting)이라면, 거대 행성이 형성되는 반경에서 가스를 제거하는 효율이 기존 가정보다 낮을 수 있으며, 이는 핵 강착(core accretion) 및 이주(migration) 모델에 중요한 시사점을 준다.

방법론
흐름이 원반 구조에 동적으로 반응할 때 광증발만으로 간극을 열고 유지할 수 있는지 조사하기 위해, 저자들은 두 단계의 계산적 접근 방식을 채택하였다:

1. 1단계 (1D 진화): 저자들은 spock 점성 진화 코드를 사용하여 간극 형성이 시작될 때까지 X선 광증발을 겪는 원반을 시뮬레이션하였다. 시뮬레이션은 Picogna 등(2019)의 질량 손실 프로파일을 포함하여 2.1 Myr 동안 가스 표면 밀도 프로파일($\Sigma$)을 진화시켰다. 이 단계는 2D 시뮬레이션을 위한 초기 조건, 구체적으로는 표면 밀도의 부분적 결핍이 발달한 상태를 설정하였다.
2. 2단계 (2D 복사-유체역학): 2.1 Myr 시점의 1D 시뮬레이션으로부터 얻은 표면 밀도 프로파일을 pluto 코드를 사용하여 축대칭 2D ($r-\theta$) 구형 좌표계로 매핑하였다. 이 시뮬레이션들은 광증발 유출과 원반 간의 상호작용을 모델링하였다.
   • Model p-g: 이미 존재하는 밀도 결핍을 포함하여 1단계에서 진화된 표면 밀도를 활용하였다.
   • Model p-ng: 간극 형성에 미치는 초기 형태의 영향을 평가하기 위해, 결핍을 억제하고(거듭제곱 법칙으로 대체) 수정된 초기 표면 밀도 프로파일을 활용하였다.
   • 시뮬레이션은 시선 방향 가열(line-of-sight heating) 근사를 채택하였으며, 계산 비용을 줄이기 위해 복사 전달을 유체역학으로부터 분리하여 복사 평형을 가정하였다.

주요 결과

• 자기 제한적 간극 형성: 2D 시뮬레이션에서 밀도 저하가 형성되면, 간극 영역 내의 국지적 질량 손실률이 급격히 감소하였다. 이는 결핍된 영역이 내부 원반의 그림자에 머물게 되어 항성 복사에 대한 노출이 줄어들기 때문에 발생하였다. 결과적으로, 간극은 완전히 비어 있는 공동(cavity)으로 깊어지지 않았다.
• 부분적 재충전: 원반 표면(및 바람을 통한)을 따르는 점성 유입과 반경 방향 질량 수송이 결핍된 영역을 부분적으로 채우는 역할을 하였다. 시스템은 깨끗한 내부 구멍(inner hole) 대신 지속적인 부분 결핍 구역을 향해 진화하였다.
  এটি।
• 초기 조건에 대한 무감각성: 초기 결핍이 있는 모델(p-g)과 없는 모델(p-ng)을 비교한 결과, 장기적인 진화 결과는 초기 원반 형태와 크게 무관함을 보여주었다. 두 경우 모두 시스템은 부분 결핍 상태로 자기 조절되었다.
• 1D 처방의 한계: 저자들은 국지적 표면 밀도 결핍에 기초하여 원시 질량 손실 프로파일을 스케일링함으로써 2D 거동을 근사하는 1차 처방(first-order prescription)을 테스트하였다. 이는 정적 처방보다는 일치도가 높았으나, 1D 모델은 여전히 2D 시뮬레이션보다 너무 빠르게 진화하고 더 깊은 간극을 생성하였다. 이러한 불일치는 원반 표면을 따른 반경 방향 질량 수송과 복잡한 차폐(shadowing) 작용을 포착하지 못하는 1D 모델의 한계를 부각시킨다.
• 관측적 특징: 가스 간극이 부분적으로 채워져 있음에도 불구하고, 저하된 영역의 외곽 가장자리에 있는 압력 극대(pressure maximum)는 여 still 먼지 입자의 효율적인 트랩 역할을 할 수 있다. 이는 완전한 가스 제거 없이도 전이 원반(transition discs, 연속파 방출에서의 내부 공동)과 유사한 관측적 특징을 생성할 수 있다.

의의 및 주장
본 논문은 광증발이 효율적으로 깨끗한 내부 공동을 깎아내고 직접적으로 전이 원반을 만든다는 표준 패러다임에 도전한다. 대신, 저자들은 광증발이 지속적인 부분 결핍 구역을 특징으로 하는 준정상 상태(quasi-steady state)로 이어질 수 있다고 제안한다. 이는 완전한 내부 제거 없이 중간 정도의 결핍을 보이는 원반이 광증발의 영향을 받았으나 완전히 분산되지는 않은 시스템을 나타낼 수 있음을 시사한다.

저자들은 1D 원반 진화 모델에서 이러한 자기 조절 거동을 근사하기 위한 적절한 수준의 1차 처방을 제공하며, 이는 행성 형성 및 인구 통계 합성 연구에 사용하기에 적합하다. 그러나 저자들은 이 처방이 여전히 근사적이며, 특히 반경 방향의 질량 재분배를 포함한 흐름의 다차원적 특성을 완전히 포착할 수 없음을 명시적으로 밝히고 있다. 본 논문은 전역적 원반 및 행성 인구 모델을 위한 견고한 레시피를 도출하기 위해 추가적인 다차원 시뮬레이션과 매개변수 공간 탐색이 필요하다고 결론짓는다. 이 연구는 내부의 바람과 행성 형성 경로 사이의 상호작용을 이해하기 위해 정적인 질량 손실 처리를 넘어설 필요가 있음을 강조한다.