Abstract:Observations of transmission spectra reveal that hot Jupiters and Neptunes are likely to possess escaping atmospheres driven by stellar radiation. Numerous models predict that magnetic fields may exert significant influences on the atmospheres of hot planets. Generally, the escaping atmospheres are not entirely ionized, and magnetic fields only directly affect the escape of ionized components within them. Considering the chemical reactions between ionized components and neutral atoms, as well as collision processes, magnetic fields indirectly impact the escape of neutral atoms, thereby influencing the detection signals of planetary atmospheres in transmission spectra. In order to simulate this process, we developed a magneto-hydrodynamic multi-fluid model based on MHD code PLUTO. As an initial exploration, we investigated the impact of magnetic fields on the decoupling of H$^+$ and H in the escaping atmosphere of the hot Neptune GJ436 b. Due to the strong resonant interactions between H and H$^+$, the coupling between them is tight even if the magnetic field is strong. Of course, alternatively, our work also suggests that merging H and H$^+$ into a single flow can be a reasonable assumption in MHD simulations of escaping atmospheres. However, our simulation results indicate that under the influence of magnetic fields, there are noticeable regional differences in the decoupling of H$^+$ and H. With the increase of magnetic field strength, the degree of decoupling also increases. For heavier particles such as O, the decoupling between O and H$^+$ is more pronounced. Our findings provide important insights for future studies on the decoupling processes of heavy atoms in the escaping atmospheres of hot Jupiters and hot Neptunes under the influence of magnetic fields.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the decoupling phenomenon of the escaping atmosphere of the hot Neptune GJ436 b under the influence of the magnetic field. Specifically, the author simulated the interaction and decoupling process of hydrogen atoms (H) and hydrogen ions (H⁺) under the action of the magnetic field by establishing a multi - fluid magnetohydrodynamic (MHD) model. The following are the core problems and background of the paper:
### Research Background
1. **Observed Phenomena**:
- Through the observation of transmission spectra, it is found that hot Jupiters and hot Neptunes may have escaping atmospheres driven by stellar radiation.
- Observations also indicate that there may be energetic neutral atoms (ENAs) in the atmospheres of these planets, which are produced by charge exchange between the stellar wind and the planetary wind.
2. **Theoretical Models**:
- Many models predict that the magnetic field has a significant impact on the atmospheric escape of hot planets.
- Since the escaping atmosphere is not fully ionized, the magnetic field only directly affects the charged components (such as H⁺) in it, and indirectly affects neutral atoms (such as H) through chemical reactions and collisions.
### Research Problems
1. **The Influence of the Magnetic Field on Atmospheric Escape**:
- How does the magnetic field affect the decoupling process of hydrogen atoms (H) and hydrogen ions (H⁺)?
- What are the effects of different magnetic field intensities on the decoupling of H and H⁺?
2. **The Mechanism of the Decoupling Phenomenon**:
- How does the decoupling degree of H and H⁺ change under different conditions (such as temperature, density, magnetic field intensity, ionization rate, etc.)?
- Does the decoupling phenomenon depend on the physical conditions in specific regions (such as the equatorial dead zone, polar region, wind zone)?
### Model and Method
In order to study these problems, the author developed a multi - fluid magnetohydrodynamic model based on the PLUTO code, focusing on simulating the decoupling process of H and H⁺ in the escaping atmosphere of GJ436 b under the influence of the magnetic field. The model considers two fluids: one represents H, and the other represents the charged fluid composed of H⁺ and electrons. The main equations of the model are as follows:
- **Continuity Equation**:
\[
\frac{\partial n_t}{\partial t}+\nabla\cdot(n_t u_t) = 0,\quad t = n, s
\]
- **Momentum Equation**:
\[
\frac{\partial(m_n n_n u_n)}{\partial t}+\nabla\cdot(m_n n_n u_n u_n + I p_n)=C_{ns}m_a m_u n_n n_s(u_s - u_n)+m_n n_n a_{ext}
\]
\[
\frac{\partial(m_s n_s u_s)}{\partial t}+\nabla\cdot(m_s n_s u_s u_s + I(p_s+\frac{B^2}{8\pi})-\frac{1}{4\pi}B B)=C_{ns}m_a m_u n_n n_s(u_n - u_s)+m_s n_s a_{ext}
\]
- **Pressure Evolution Equation**:
\[
\frac{\partial p_t}{\partial t}+(u_t\cdot\nabla)p_t+\gamma p_t(\nabla\cdot u_t) = 0,\quad t = n, s
\]
- **Magnetic Field Evolution Equation**:
\[
\frac{\partial B}{\partial t}=\nabla\times(u_s\times B)
\]
### Results and Discussion
Through numerical simulation, the author found:
1. **The Influence of Magnetic Field Intensity**:
- As the magnetic field intensity increases, the decoupling degree of H and H⁺ also increases.
- In the polar region, as the magnetic field intensity increases, the fluid velocity gradually decreases and finally forms a dead zone.
2. **The Influence of Temperature**:
- When the temperature rises, the gas pressure gradually cancels out the magnetic field pressure, leading...