DOI: https://doi.org/10.1051/0004-6361/202557579
تاريخ النشر: 2026-01-22
المؤلف: Ayumu Kuwahara وآخرون
الموضوع الرئيسي: علوم الفضاء والكواكب
نظرة عامة
في هذا البحث، يستكشف المؤلفون ديناميات الأغلفة الكوكبية التي تتشكل من خلال تراكم النواة، مع التركيز على التفاعل بين التسخين الناتج عن تراكم المواد الصلبة وعمليات التبريد. لقد عالجت النماذج التقليدية أحادية البعد هذه الأغلفة على أنها ثابتة، ولكن المحاكاة الهيدروديناميكية ثلاثية الأبعاد الأخيرة تكشف أن تبادل الغاز بين الغلاف والقرص المحيط ديناميكي ويتأثر بمعدلات التبريد، التي تتميز بالمعامل غير البعدي $\beta$. تحدد الدراسة ثلاثة أنظمة تبريد متميزة: التبريد السريع (مع $\beta \lesssim 1$)، حيث يبقى الغلاف تقريبًا متساوي الحرارة؛ التبريد المتوسط ($1 \lesssim \beta \lesssim 300$)، الذي يتميز بهيكل ثلاثي الطبقات يتكون من طبقات الحمل، الإشعاع، وإعادة التدوير؛ والتبريد البطيء ($\beta \gtrsim 1000$)، حيث يصبح الغلاف كامل الحمل.
تشير النتائج إلى أن الأغلفة في المناطق الداخلية من الأقراص الكوكبية الأولية (داخل $\sim 1$ وحدة فلكية) من المحتمل أن تكون كاملة الحمل، مما يؤدي إلى تبادل فعال للمواد ولكن قد يؤدي إلى نقص في المواد المتطايرة. في المقابل، فإن الأغلفة في القرص الخارجي (ما بعد $\sim 10$ وحدة فلكية) أكثر عرضة لتطوير طبقات إشعاعية، مما يسمح بتراكم المواد المتطايرة المتبخرة. تؤكد البحث على أهمية النظر في الجوانب الديناميكية والديناميكية لتطور الغلاف لفهم نمو الكواكب وتركيبها، مما يشير إلى أن النماذج المستقبلية يجب أن تدمج تطور الغبار، والنقل الإشعاعي، والديناميات متعددة السوائل لتعزيز فهمنا لعمليات تشكيل الكواكب.
مقدمة
في سياق نموذج تراكم النواة، تحقق هذه الدراسة في تشكيل وديناميات الأغلفة الكوكبية داخل الأقراص الكوكبية الأولية. في البداية، تقوم الأجنة الكوكبية بتراكم الكويكبات والحصى، وعندما تصل إلى كتلة تتجاوز حوالي 10 كتل أرضية ($M_\oplus$)، تبدأ في التقاط الغاز، مما يشكل غلافًا. تعالج النماذج الكلاسيكية أحادية البعد هذه الأغلفة على أنها ثابتة ومعزولة؛ ومع ذلك، تشير المحاكاة الهيدروديناميكية ثلاثية الأبعاد الأخيرة إلى أنها تتفاعل ديناميكيًا مع القرص المحيط من خلال عملية تعرف بإعادة التدوير. تسهل هذه التفاعل تبادل الغاز المستمر، مما يمكن أن يؤثر على التبريد والانكماش للغلاف، وبالتالي يؤثر على توقيت تراكم الغاز السريع.
تخضع نمو الغلاف للتوازن بين التسخين الناتج عن المواد الصلبة المتراكمة والتبريد من خلال الإشعاع أو الحمل، حيث يتأثر كفاءة إعادة التدوير بشكل كبير بكفاءة التبريد الإشعاعي للغاز. تبرز الدراسة أن ديناميات الغاز يمكن أن تؤثر أيضًا على معدلات تراكم المواد الصلبة، وخاصة الحصى، من خلال تقليل التقاطها بسبب سحب الغاز وقد تؤدي إلى تبخر الأنواع المتطايرة. لمعالجة هذه التعقيدات، يقدم المؤلفون تحقيقًا شاملاً لديناميات الغلاف عبر معدلات تسخين وتبريد مختلفة، باستخدام نموذج مبسط لتحديد كفاءة التبريد الإشعاعي. يحدد البحث هيكله، موضحًا الإعداد العددي للمحاكاة، وتصنيف أنظمة التبريد، ونمذجة كفاءة إعادة التدوير للجزيئات والبخار داخل الغلاف.
طرق
في هذا القسم، يصف المؤلفون الطرق العددية المستخدمة لمحاكاة عمليات التسخين والتبريد لكوكب يتراكم الحصى داخل قرص غير جاذب ذاتيًا باستخدام كود Athena++ (Stone et al. 2020). تُجرى المحاكاة في إحداثيات قطبية كروية مركزة على الكوكب، مع تعريف المتغيرات كـ $r$ (المسافة من الكوكب)، $\theta$ (زاوية قطبية)، و$\phi$ (زاوية أفقية). يحل الكود معادلات الاستمرارية، والزخم، والطاقة لسائل قابل للضغط، غير لزج، مع تضمين مصطلحات للاسترخاء الحراري وتسخين التراكم.
يتم نمذجة الجاذبية للكوكب باستخدام دالة ملساء لضمان التوازن الهيدروستاتيكي عند الحدود الداخلية، مما يتجنب التسخين العددي الذي قد ينشأ من صياغة الجاذبية التقليدية. يتم تعريف دالة التنعيم كـ $f_{sm} = \frac{(r – r_{in})^2}{(r – r_{in})^2 + r_{sm}^2}$، حيث $r_{sm} = 0.1 r_{in}$. بالإضافة إلى ذلك، يتم إدخال التأثير الجاذبي للكوكب تدريجيًا في القرص لتخفيف تكوين الصدمات، باستخدام دالة حقن $f_{inj} = 1 – \exp\left(-\frac{1}{2} \frac{t}{t_{inj}}^2\right)$، مع $t_{inj} = 0.5 \Omega^{-1}$، حيث $\Omega$ هو التردد المداري.
نتائج
تصنف نتائج الدراسة الحالات (الحرارية) الديناميكية للأغلفة الكوكبية إلى ثلاثة أنظمة متميزة بناءً على معدلات تبريدها: نظام التبريد السريع ($\beta < 1$)، نظام التبريد المتوسط ($1 \lesssim \beta < 1000$)، ونظام التبريد البطيء ($\beta \gtrsim 1000$). يظهر كل نظام خصائص فريدة، يتم توضيحها في الأقسام اللاحقة (3.1-3.3) لحالة معيارية بكتلة كوكبية قدرها $m = 0.1$ وإضاءة تراكم قدرها $L_{\text{acc}} = 1.1 \times 10^{-5} (m/0.1)^{5/3} (t_{\text{acc}} / 10^5 \text{ سنة})^{-1}$. سيتم مناقشة المزيد من التحليل بشأن تأثير الكتلة الكوكبية وإضاءة التراكم على هذه الأنظمة في القسم 3.4، مما يوفر فهمًا شاملاً لكيفية تأثير هذه المعلمات على الديناميات الحرارية للأغلفة الكوكبية.
مناقشة
في هذا القسم، يناقش المؤلفون ديناميات الأغلفة الغازية المحيطة بالكواكب المدمجة، مع التركيز على تأثيرات عمليات التبريد والتسخين. يتم نمذجة تبريد الغاز كعملية استرخاء حراري، مع التعبير عن معدل التبريد كـ \( Q_{\text{cool}} = -e_0 e^{-t/t_{\text{cool}}} \). تفحص الدراسة التسخين الناتج عن تراكم الحصى، الممثل بإضاءة \( L(r) \approx GM_p \dot{M} r \)، وتستخرج معدل التسخين لكل وحدة حجم \( Q_{\text{acc}} \) من معادلة نقل الطاقة. تستخدم المحاكاة وحدات طبيعية وتستكشف معدلات تراكم مختلفة، مما يكشف أن زمن التبريد، الذي يتم تحديده كـ \( \beta \equiv t_{\text{cool}} \Omega \)، يؤثر بشكل كبير على الهيكل الحراري للغلاف.
يحدد المؤلفون ثلاثة أنظمة تبريد متميزة: التبريد السريع، المتوسط، والبطيء. في نظام التبريد السريع (\( \beta \lesssim 1 \))، يكون الغلاف تقريبًا متساوي الحرارة، مع تدرج حراري مستقر. يظهر النظام المتوسط (\( 1 \lesssim \beta \lesssim 300 \)) هيكلًا ثلاثي الطبقات يتكون من طبقات الحمل، الإشعاع، وإعادة التدوير، حيث تنتقل ملفات الحرارة والكثافة بين الحالات المتساوية الحرارة والحرارية. في نظام التبريد البطيء (\( \beta \gtrsim 1000 \))، يصبح الغلاف كامل الحمل، متطابقًا مع الحد الحراري. تبرز الدراسة أن ديناميات الغلاف تحكمها توازن التسخين والتبريد، مع زيادة إضاءة التراكم التي تؤدي إلى تدرجات حرارية أكثر حدة وزيادة الحمل. تشير النتائج إلى أن هيكل الغلاف قوي ضد التغيرات في الكتلة الكوكبية، ولكن زمن التبريد الحرج الذي يفصل بين الأنظمة حساس لإضاءة التراكم.
القيود
تنبع قيود الدراسة بشكل أساسي من فرضيتين رئيسيتين: نموذج تبريد β وإغفال تدرجات الضغط في غاز القرص. نموذج تبريد β، الذي يبسط التبريد الإشعاعي من خلال تعيين قيمة β ثابتة، لا يأخذ في الاعتبار التغير في الشفافية بسبب جزيئات الغبار الصغيرة في الغلاف الكوكبي. بينما يكون النموذج قابلًا للتطبيق في ظل ظروف الشفافية المنخفضة (κ ≲ $10^{-5}$ إلى $10^{-4}$ سم²/غ) أو على مسافات أكبر من حوالي 10 وحدة فلكية، توجد عدم يقين كبير بشأن الشفافية الفعلية، والتي يمكن أن تتأثر بنمو الغبار واستقراره. في المناطق الأقرب إلى الكوكب (داخل ≲ 0.1 $R_B$)، قد تنخفض الشفافية إلى κ ≲ $10^{-4}$ سم²/غ بسبب تكتل الغبار، مما قد يؤثر على دقة النموذج.
بالإضافة إلى ذلك، تفترض الدراسة دوران كيبلري خالص، متجاهلة تدرج الضغط لغاز القرص الخلفي. من المتوقع أن يكون لهذا التبسيط تأثير طفيف على النتائج، حيث يظهر غاز القرص الحقيقي دورانًا دون كيبلري، مما يقدم رياحًا معاكسة قد تقلل من حجم الطبقة الإشعاعية في الغلاف. يتم وصف العلاقة بين نصف قطر الغلاف وتأثير الرياح المعاكسة بالمعادلة \( R_{KK24}^{atm} = R_{atm,0} (1 – D_1 \frac{\eta v_K}{c_s} \frac{R_B}{H}) \)، حيث \( R_{atm,0} = C_1 R_B \) و\( C_1 \) و\( D_1 \) هما معاملات ملائمة. تحت الظروف النموذجية، يمكن أن تقلل هذه الرياح المعاكسة من حجم الغلاف الجوي بنسبة تصل إلى حوالي 7%.
DOI: https://doi.org/10.1051/0004-6361/202557579
Publication Date: 2026-01-22
Author(s): Ayumu Kuwahara et al.
Primary Topic: Astro and Planetary Science
Overview
In this research, the authors explore the dynamics of planetary envelopes formed through core accretion, emphasizing the interplay between heating from solid accretion and cooling processes. Traditional one-dimensional models have treated these envelopes as static, but recent three-dimensional hydrodynamical simulations reveal that gas exchange between the envelope and the surrounding disk is dynamic and influenced by the cooling rates, characterized by the dimensionless parameter $\beta$. The study identifies three distinct cooling regimes: fast cooling (with $\beta \lesssim 1$), where the envelope remains nearly isothermal; intermediate cooling ($1 \lesssim \beta \lesssim 300$), which features a three-layer structure comprising convective, radiative, and recycling layers; and slow cooling ($\beta \gtrsim 1000$), where the envelope becomes fully convective.
The findings indicate that envelopes in the inner regions of protoplanetary disks (within $\sim 1$ au) are likely to be fully convective, leading to efficient material exchange but potential depletion of volatiles. In contrast, envelopes in the outer disk (beyond $\sim 10$ au) are more prone to develop radiative layers, allowing for the accumulation of sublimated volatiles. The research underscores the importance of considering both the thermodynamic and dynamic aspects of envelope evolution in understanding planetary growth and composition, suggesting that future models should integrate dust evolution, radiative transfer, and multi-fluid dynamics to enhance our comprehension of planet formation processes.
Introduction
In the context of the core accretion model, this study investigates the formation and dynamics of planetary envelopes within protoplanetary disks. Initially, planetary embryos accrete planetesimals and pebbles, and once they reach a mass exceeding approximately 10 Earth masses ($M_\oplus$), they begin to capture gas, forming an envelope. Classical one-dimensional models treat these envelopes as static and isolated; however, recent three-dimensional hydrodynamical simulations indicate that they interact dynamically with the surrounding disk through a process known as recycling. This interaction facilitates continuous gas exchange, which can influence the cooling and contraction of the envelope, thereby affecting the timing of runaway gas accretion.
The growth of the envelope is governed by the balance between heating from accreted solids and cooling through radiation or convection, with the efficiency of recycling being significantly influenced by the radiative cooling efficiency of the gas. The study highlights that gas dynamics can also impact the accretion rates of solids, particularly pebbles, by reducing their capture due to gas drag and potentially leading to the sublimation of volatile species. To address these complexities, the authors present a comprehensive investigation of envelope dynamics across various heating and cooling rates, utilizing a simplified model to parameterize radiative cooling efficiency. The paper outlines its structure, detailing the numerical setup for simulations, categorization of cooling regimes, and modeling of recycling efficiency for particles and vapors within the envelope.
Methods
In this section, the authors describe the numerical methods employed to simulate the heating and cooling processes of a pebble-accreting planet within a non-self-gravitating disk using the Athena++ code (Stone et al. 2020). The simulations are conducted in spherical polar coordinates centered on the planet, with variables defined as $r$ (distance from the planet), $\theta$ (polar angle), and $\phi$ (azimuth angle). The code solves the equations of continuity, momentum, and energy for a compressible, inviscid fluid, incorporating terms for thermal relaxation and accretion heating.
The gravitational potential of the planet is modeled using a smoothed function to ensure hydrostatic equilibrium at the inner boundary, avoiding numerical heating that would arise from a typical gravitational potential formulation. The smoothing function is defined as $f_{sm} = \frac{(r – r_{in})^2}{(r – r_{in})^2 + r_{sm}^2}$, where $r_{sm} = 0.1 r_{in}$. Additionally, the gravitational influence of the planet is gradually introduced into the disk to mitigate shock formation, utilizing an injection function $f_{inj} = 1 – \exp\left(-\frac{1}{2} \frac{t}{t_{inj}}^2\right)$, with $t_{inj} = 0.5 \Omega^{-1}$, where $\Omega$ is the orbital frequency.
Results
The results of the study categorize the (thermo)dynamic states of planetary envelopes into three distinct regimes based on their cooling rates: the fast cooling regime ($\beta < 1$), the intermediate cooling regime ($1 \lesssim \beta < 1000$), and the slow cooling regime ($\beta \gtrsim 1000$). Each regime exhibits unique properties, which are elaborated upon in subsequent sections (3.1-3.3) for a fiducial case with a planetary mass of $m = 0.1$ and an accretion luminosity of $L_{\text{acc}} = 1.1 \times 10^{-5} (m/0.1)^{5/3} (t_{\text{acc}} / 10^5 \text{ yr})^{-1}$. Further analysis regarding the influence of planetary mass and accretion luminosity on these regimes will be discussed in Section 3.4, providing a comprehensive understanding of how these parameters affect the thermal dynamics of planetary envelopes.
Discussion
In this section, the authors discuss the dynamics of gas envelopes surrounding embedded planets, focusing on the effects of cooling and heating processes. The cooling of the gas is modeled as a thermal relaxation process, with the cooling rate expressed as \( Q_{\text{cool}} = -e_0 e^{-t/t_{\text{cool}}} \). The study examines the heating due to pebble accretion, represented by a luminosity \( L(r) \approx GM_p \dot{M} r \), and derives the heating rate per unit volume \( Q_{\text{acc}} \) from the energy transport equation. The simulations utilize normalized units and explore various accretion rates, revealing that the cooling timescale, parameterized as \( \beta \equiv t_{\text{cool}} \Omega \), significantly influences the envelope’s thermal structure.
The authors identify three distinct cooling regimes: fast, intermediate, and slow cooling. In the fast cooling regime (\( \beta \lesssim 1 \)), the envelope is nearly isothermal, with a stable temperature gradient. The intermediate regime (\( 1 \lesssim \beta \lesssim 300 \)) exhibits a three-layer structure comprising convective, radiative, and recycling layers, where the temperature and density profiles transition between isothermal and adiabatic states. In the slow cooling regime (\( \beta \gtrsim 1000 \)), the envelope becomes fully convective, matching the adiabatic limit. The study highlights that the dynamics of the envelope are governed by the balance of heating and cooling, with higher accretion luminosity leading to steeper temperature gradients and enhanced convection. The results suggest that the envelope structure is robust against variations in planetary mass, but the critical cooling timescale separating regimes is sensitive to the accretion luminosity.
Limitations
The limitations of the study primarily stem from two key assumptions: the β-cooling model and the neglect of pressure gradients in the disk gas. The β-cooling model, which simplifies radiative cooling by assigning a constant β value, does not account for the variability in opacity due to small dust grains in the planetary envelope. While the model is applicable under low opacity conditions (κ ≲ $10^{-5}$ to $10^{-4}$ cm²/g) or at distances greater than approximately 10 AU, significant uncertainties exist regarding the actual opacity, which can be influenced by dust growth and settling. In regions closer to the planet (within ≲ 0.1 $R_B$), opacity may decrease to κ ≲ $10^{-4}$ cm²/g due to dust coagulation, potentially affecting the model’s accuracy.
Additionally, the study assumes pure Keplerian rotation, disregarding the pressure gradient of the background disk gas. This simplification is expected to have a minor impact on the results, as real disk gas exhibits sub-Keplerian rotation, which introduces a headwind that could reduce the size of the radiative layer in the envelope. The relationship between atmospheric radius and the headwind effect is described by the equation \( R_{KK24}^{atm} = R_{atm,0} (1 – D_1 \frac{\eta v_K}{c_s} \frac{R_B}{H}) \), where \( R_{atm,0} = C_1 R_B \) and \( C_1 \) and \( D_1 \) are fitting coefficients. Under typical conditions, this headwind can reduce the atmospheric size by up to approximately 7%.