Solar Wind Charge-Exchange X-ray Emissions from the O⁵⁺ Ions in the Earth’s Magnetosheath

1. Introduction

The Earth's magnetic field is compressed by solar wind to form a magnetosphere. The boundary layer between the magnetosphere and the magnetosheath flow is the magnetopause. The magnetosheath is confined between the most two important boundaries of the Earth’s magnetosphere, namely the bow shock and the magnetopause. Pioneer 1, launched in 1958, was the first to conduct in-situ detection of the Earth's magnetosheath [1]. With advances in aerospace and computer technology, in situ measurement data and computer numerical simulations are used to gain a more complete understanding of the Earth's magnetosheath. The magnetosheath plasma comes from the solar wind heated after the bow shock. The density decreases with the distance from the Earth's center, but it is always higher than the magnetospheric plasma density. Earth's particles reach the magnetosheath by magnetic reconnecting at the magnetopause. Since the ion-electron temperature ratio of the magnetosheath is 6 to 7 times that of the solar wind, the temperature of the solar wind plasma increases even more after entering the magnetosheath, while the electron temperature does not change significantly. The temperature anisotropy of the magnetosheath plasma becomes stronger near the top of the magnetosphere, and the ion temperature anisotropy is much stronger than the temperature anisotropy. Nabert et al. [2] used global MHD simulations to obtain the distribution of magnetosheath plasma density and velocity. They found that the magnetosheath plasma flows from the subsolar point to the two flanks, and the plasma density decreases as the distance from the center of the earth increases. Dimmock and Nykyri [3] used 5-year measurement data from THEMIS satellite to confirm the reliability of Nabert et al.'s MHD model and found that the magnetosheath plasma parameters were axially asymmetric. Using global hybrid simulations, Omidi et al. [4] found differences in the magnetosheath density perturbations for different Mach numbers. Interplanetary shocks and fractures keep the plasma and magnetic field in the magnetosheath in a changing state and are related to the type of bow shock perturbations. The magnetic field perturbation is severe downstream of the quasi-parallel bow shock but is smaller in the case of the quasi-perpendicular bow shock. There are various free energies (such as ion beam current, temperature anisotropy, etc.) in the magnetosheath, which can excite various plasma waves or instabilities [5,6]. The magnetic field in the magnetosheath profoundly affects the solar wind-magnetosphere interactions. For example, when the magnetosheath magnetic field has a southward component, magnetic field reconnection will occur at the top of the low-latitude magnetosphere. At this time, the energy and momentum transferred to the magnetosphere, plasma and magnetic flux increased significantly.

The SWCX [7] mechanism is that highly charged solar wind heavy ions (such as C⁵⁺, C⁶⁺, N⁷⁺, O⁶⁺, O⁷⁺ and O⁸⁺) collide with neutral particles from the Earth or comets, and the solar wind heavy ions obtain an electron. The electrons are at an excited state, and then the excited solar wind ions de-excite and transmit photons in extreme ultraviolet (EUV) [8] or soft X-rays (Soft X-ray, SXR) [9]. Regardless of the line of sight, X-ray emission produced by SWCX in interplanetary space contaminates every astrophysical observation. However, the main emission spectral lines of SWCX also happen to be an important diagnostic method for astrophysical plasmas [10].

Zhou et al. [11] used the Suzaku X-ray satellite to observe the SWCX event driven by the corona mass ejection (CME) and found that the arrival time of SWCX and CME magnetic clouds were consistent. They simulated the light curve changes of SWCX and concluded that the solar wind has anisotropy in the extreme tip region. Ishikawa et al. [12] analyzed the intense SWCX X-ray data obtained by the Suzaku satellite during geomagnetic storms and found the time-varying changes in diffuse SXR emission related to the solar wind proton flux. Ringuette et al. [13] successfully separated the SWCX emissions from the X-ray spectrum obtained by HaloSat at low ecliptic latitudes. Zhang et al. [14] compared the ROSAT all-sky survey (RASS) data and the Quiet method (a quiet observation data method, which uses the same satellite to observe the same target for a long time) using XMM-Newton when solar wind conditions change dramatically. XMM-Newton observed SWCX emissions from the magnetosheath. The results show that the SWCX spectral composition of the magnetosheath obtained is close for the two methods, and the temporal changes of the intensity are similar. However, significant differences exist at energies lower than 0.7 keV between the two methods. Connor et al. [15] used global magnetohydrodynamics model to simulate the strong SXR emissions in the magnetosheath and cusps regions and proposed that the magnetopause reconnection can be investigated with the X-ray images. Sun et al. [16] simulated different solar wind conditions, and the virtual camera observed X-ray images of the Earth's magnetosheath and cusps regions. These images can capture the main responses of the magnetopause and cusps regions. Therefore, soft X-ray large-scale imaging technology can be used to perform long-distance panoramic imaging of the magnetosphere, thereby understanding the basic mode of solar wind-magnetosphere interaction on a large scale. Based on these simulations, the Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) [17] will conduct X-ray imaging of the magnetosheath. Collier et al. [18] provided the development progress of the wide-field SXR imager (Soft X-ray Imager, SXI) for the ESA AXIOM mission. Somana et al. [19] gave the design, detection, optimization, performance prediction and further work of the CCD370 equipped with SXI on the SMILE satellite. The SMILE satellite is planned to be launched in 2025.

The solar wind-magnetosphere interaction produces a variety of plasma physical processes, such as magnetic field reconnection, turbulence, particle acceleration and heating, etc., which subsequently triggers various space weather effects such as magnetic storms, substorms, and auroral brightening, ultimately affecting satellite navigation, and having a profound impact on terrestrial communications, power grid maintenance and human life. Therefore, strengthening the research on changes in the solar-terrestrial space environment is not only of scientific significance, but also of practical significance.

In this paper, a three-dimensional (3-D) global hybrid model will be used to obtain solar wind heavy ion parameters (including the O⁶⁺ ion number density $n_{O^{6+}}$, the solar wind speed V_SW, and temperature T_SW). The distribution of neutral hydrogen atoms near the Earth was obtained through the 3-D geocoronal hydrogen model. The SXR spectrum and corresponding efficiency factor α produced by the collision between O⁶⁺ ions and hydrogen atoms were obtained through the Spectral Analysis System for Astrophysical and Laboratory (SASAL). We then establish a spatial coordinate system and bring the above parameters into the X-ray emission integral equation, and finally simulated with a virtual X-ray camera.

We then extract the efficiency factor data of the spectral line with the strongest emission intensity and simulate the magnetosheath X-ray images at two special positions a: (0.0, 0.0, 60.0 R_E) and b: (0.0, 60.0 R_E, 0.0), at six stages (Stage 1-6) during the passage of a fast forward perpendicular interplanetary shock (IP shock) embedded with a tangential discontinuity (TD) with it normal along the X_GSM direction (Hereafter referred to as IP-TD). Finally, the positions of the bow shock and the magnetopause on the Sun-Earth line are extracted from the hybrid simulation results and compare with the SXR images to explore the feasibility of SXR imaging of the magnetosheath.

2. Methods

The SXR spectrum in the SWCX process is determined by the efficiency factor of each particle type, charge state, and energy level transition participating in the reaction. Therefore, this section will introduce the SWCX mechanism and the relevant models.

2.1. SWCX Mechanism

O⁶⁺ ions are the main charge state of oxygen ions in both the fast and slow solar winds [20], and the reaction rate of O⁶⁺ ions is significantly higher than other charge states. The collisions between O⁶⁺ ions and hydrogen atoms produce O^5+* ions (excited state) and hydrogen ions. The reaction equation is as follows:

$$O^{6+}+H\to O^{{5+}^{\mathrm{*}}}+H^{+}$$

(1)

During the decay of O^5+* to the low energy state O⁵⁺, the corresponding X-ray emission intensity (I, keV cm^-2 s^-2 sr^-1) can be calculated by the following integration along line of sight (LOS):

$$I=\frac{1}{4\pi }\underset{LOS}{\overset{\mathrm{\infty }}{\int }}{n}_{O^{6+}}\left(s\right)u_{O^{6+}}\left(s\right)n_H\left(s\right)\alpha \left(u_{O^{6+}}\right)ds,$$

(2)

where s=s(r,θ,φ) represents a point in LOS in the Geocentric Solar Magnetospheric (GSM) coordinate system, $n_{O^{6+}}$ is the density of O⁶⁺ ions, $u_{O^{6+}}$ is the relative collision speed between O⁶⁺ ions and hydrogen atoms, n_H is the exospheric hydrogen density, and α is the efficiency factor. Normally, the speed of exospheric hydrogen atoms is much slower than the solar wind ions, and the speed of O⁶⁺ ions can be represented by the solar wind speed V_SW. Besides, the thermal velocity V_th mush be considered in the magnetosheath region [21]. The formula for calculating the thermal velocity is

$$V_{th}=\sqrt{\frac{3k_B{T}_{sw}}{m_O}},$$

(3)

where k_B is Boltzmann's constant, and m_O is the mass of O⁶⁺ ions. The relative collision speed of O⁶⁺ ions and hydrogen atoms is

$$u_{O^{6+}}=\sqrt{V_{sw}^2+V_{th}^2}.$$

(4)

2.2.3-D Global Hybrid Model

In the 3-D global hybrid model treats the ions as discrete particles and the electrons as massless fluids [22,23,24]. Both the dayside and nightside magnetosphere are included in the model. By assuming quasi-electric neutrality, the motion equation of particles in the simulation unit is:

$$\frac{{d\mathit{v}}_p}{dt}=\mathit{E}+{\mathit{v}}_p\times \mathit{B}-\nu \left({\mathit{V}}_p-{\mathit{V}}_e\right),$$

(5)

where v_p is the ion velocity, B is the magnetic field, E is the electric field, ν is the collision frequency related to the current, and V_p and V_e are the bulk flow velocities of ions and electrons, respectively. In the inner magnetosphere region (r<6.5 R_E) the number density of cold ionic fluid is

$$N_f=\left(1000/r^3\right)[1-\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{h}(r-6.5\left)\right],$$

(6)

where r is the distance to the center of the Earth in R_E.

The input parameters for the hybrid model are V_SW, N_SW, T_SW, interplanetary magnetic field (IMF) components Bx, By and Bz, and the ratio of O⁶⁺ ions in the solar wind $f_{O^{6+}}$. According to previous studies, the He²⁺/H⁺ ratio is ~5% [25], and the He²⁺/O⁶⁺ ratio is ~50 [26], and the O/H flux ratio observed by Ulysses is ~0.001 [27]. Therefore, throughout the study, we adopt the value of $f_{O^{6+}}$=0.001.

2.3. 3-D Geocoronal Hydrogen Model

The Monte Carlo simulation method of planetary exosphere was developed from early studies of the lunar atmosphere. The regolith displayed by the moon provides a clear exosphere [28]. Assuming that particles in the planet's exosphere do not collide with each other, the atoms leaving the exosphere obey the Maxwell velocity distribution. Particles distributed in the high-energy tail will break away from the planet and escape from the planet. The remaining particles pass through the exosphere on ballistic trajectories and then return to the exobase. Brinkmann [29] studied the collision between the neutral particles at the bottom of the exosphere and the escaping neutral particles and found that since the downward velocity of H or He atoms when they returned to the distribution at the bottom of the exosphere did not exceed the escape velocity, This shows that there is energy dissipation after the neutral particles escape. Cole [30] believed that the charge exchange between hot hydrogen atoms and hot protons in the ionosphere would increase the escape of hydrogen atoms. Subsequent calculations by Hodges et al. [31] provided quite strong arguments for charge exchange being the cause of hydrogen escaping from the planets.

This article adopts the Hodges model [28]. Compared with other models, Hodges' 3-D Monte Carlo model considered the interaction process between neutral particles and the interaction process between ions and neutral particles. At a certain position (r, θ, φ), the density of hydrogen atoms is

$$n_H\left(r,\theta ,\phi \right)=N\left(r\right)\sqrt{4\pi }\sum _{l=0}^3\sum _{m=0}^l\left\{A_{lm}\cos \left(m\varphi \right)+B_{lm}\sin \left(m\varphi \right)\right\}{Y}_{lm}\left(\theta \right),$$

(7)

where N(r), A_lm and B_lm are given by Hodges [28], and Y_lm(θ) is the spherical harmonic Legendre function.

2.4. SASAL

The SASAL [32] adopts the optically thin assumption (when light passes through a thin transparent medium, its propagation path can be approximately regarded as linear propagation) to describe the atomic model of solar wind heavy ion spectral line emission, including collision excitation of solar wind heavy ions, electron/photon ionization, dielectronic/radiative recombination (DR/RR), charge transfer caused by collision with neutral ions, etc.

The SASAL uses the nuclear number from H to Zn ion data in Chianti version 7 as the baseline data for this model. In terms of calculating photoionization, the data provided by Verner et al. [33] were added to SASAL, and the analysis method used was consistent with the non-relativistic calculation of the ground state of atoms and ions. For the cross section of charge exchange, the multichannel Landau–Zener theory with rotational coupling [34] was adopted in the SASAL.

In this article, SASAL can be simplified accordingly, that is, only the collision of heavy ions and neutral atoms is considered. Electrons are usually captured to an energy level with a given quantum number nq, which can be expressed as [35]:

$$n_q=q{\left(2\left|I_H\right|\left(1+\frac{q-1}{2\sqrt{q}+1}\right)\right)}^{-\frac{1}{2}},$$

(8)

where I_H is the ionization potential of hydrogen atoms. Liang et al. [35] gave a detailed calculation process. Then, by solving the following rate equation [36], the O⁵⁺ density of the i-th energy state ($N_i^{5+}$) is obtained through the equation:

$$\begin{gathered} \frac{d}{dt}{N}_i^{5+}\left(r\right)=\sum _{j>i}{N}_i^{5+}\left(r\right)A_{ji}-\sum _{j<i}{N}_i^{5+}{A}_{ij}+n_H\left(r\right)n_{O^{6+}}{C}_{0i}\left(v\right) \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}+n_H\sum _{j\ne i}\left[N_i^{5+}\left(r\right)Q_{ji}\left(T_e\right)-N_j^{5+}\left(r\right)Q_{ij}\left(T_e\right)\right]=0, \end{gathered}$$

(9)

where $C_{ij}\equiv 〈v{\sigma }_{cx}\left(v\right)〉$ is the single electron capture rate coefficient, v is the relative collision speed of O⁶⁺ ions and H, ${\sigma }_{cx}\left(v\right)={\sigma }_{O^{6+}}\left(u_{O^{6+}}\right)$ is the cross section of the single electron transfer process, A_ij is the radiative attenuation rate of a specific transition i→j, and Q_ij(T_e) is the proton impact excitation rate at a given temperature T_e. Once $N_i^{5+}$ of a given i energy state is determined, the efficiency factor corresponding to a given transition i→j can be determined as

$${\alpha }_{ij}={\sigma }_{O^{6+}}\left(u_{O^{6+}}\right)\frac{N_i^{5+}}{n_{O^{6+}}}{∆E}_{ij}$$

(10)

where ${∆E}_{ij}$ is the transition energy.

3. Results

3.1. Distributions of O⁶⁺ Ions and Hydrogen Atoms

An example of the distributions of $n_{O^{6+}}$, V_SW, and temperature T_SW for solar wind conditions of N_SW =6 cm^-3, V_SW =700 km s^-1, T_SW =10 eV, IMF magnitude B_T=8 nT, IMF clock angle 𝜃_c=0° (northward), and $f_{O^{6+}}$=0.001 are shown in Figure 1 a-i. The distributions of hydrogen atoms obtained from the 3-D geocoronal hydrogen model are shown in Figure 1 j-l.

In Figure 1 a-c, the O⁶⁺ number density reaches 0.3 cm^-3 in the dayside magnetosheath. Clear bow shock and magnetopause structures are shown. There are obvious O⁶⁺ ions transported along the magnetic field lines into the magnetosphere, reaching a density peak of 0.5 cm^-3 at about 3 R_E. The solar wind speed is significantly slowed down downstream of the bow shock, especially in the subsolar region (Figure 1 d-f). From Figure 1d, a region with fast O⁶⁺ ion speed can be seen at X_GSM ≤ -15 R_E on the night side of the Earth. This is due to the dawn-dusk asymmetry of the plasma sheet. More high-speed reconnection flows occur on the dusk side than on the dawn side, and O⁶⁺ ions are possibly accelerated through reconnection [22].

3.2. Spectrum of O⁶⁺ Ions during SWCX Process

Through the SASAL model mentioned above, the emission spectrum of O^5+* ion can be calculated as shown in Table 1. The excited state O^5+* ions are generated by collisions between O⁶⁺ ions and hydrogen atoms through the SWCX process, and then the O^5+* ions de-excite from one higher energy state to lower (indicated by → in the table) to generate the SXR photons with corresponding energy. It can be seen from the table that the spectral lines of O^5+* are in a narrow energy range between 105.607~118.291 eV. The 106.561 eV SXR produced by the O^5+* ion deexcitation process 5 in the table is consistent with the SXR produced on comets through the SWCX process of O⁶⁺ [38].

Figure 2 shows the curves of the efficiency factor 𝛼 for various spectral lines as a function of collision speed. The 𝛼 value for 11.582 nm is the highest compared with other spectral lines. The 𝛼 values for 11.582 nm and 11.583 nm lines are very sensitive to the collision speed. As the collision speed increases, their 𝛼 values can be reduced to up to 1/3 of the original values. For other spectral lines, no obvious correlation is found between the 𝛼 values and the collision speed, and their 𝛼 values are always at a low level.

Finally, the plasma parameters in the magnetosheath and the efficiency factors are input into the integration equation (Equation 2) to calculate the line-of-sight emission intensity. A virtual SXR camera is set at position b: (0.0, 60.0 R_E, 0.0). The field of view of the camera is 40º and the spatial resolution is 0.2º in the projection plane of GSM xz. The simulated images for the nine spectral lines are shown in Figure 3.

Taking Figure 3 c as an example, it can be clearly seen that there is a red broadband area (magnetosheath) that is symmetrical to the X_GSM axis, and its left boundary is the bow shock. The strongest emission intensity on the dayside is generated near the cusp region, and the strongest emission intensity on the nightside is located in the area near the magnetopause in the magnetosheath. The simulation results of SXR emission intensity in the magnetosheath of different O^5+* ion energy states have almost the same distribution pattern. The small differences may be due to the relationship of the efficiency factor to the collision speed.

The strongest emission intensities for the nine spectral lines on the dayside (red) and nightside (blue) are shown in Figure 4. When the O^5+* ion energy state is deexcited from 1s² 4d 2P1/2 to 1s² 2p 2D3/2, the emission intensity is the strongest, and the spectral line energy is 107.047 eV. The second largest emission intensity is the deexcitation of O^5+* ions from 1s² 4p 2P1/2 to 1s² 2p 2P3/2, with a spectral line energy of 107.040 eV. The strongest emission intensity is 6.07 keV cm^-2 s^-1 sr^-1, and the difference between the strongest emission intensity between day and night is basically maintained at half an order of magnitude.

3.3. Effects of Interplanetary Shock Waves on Soft X-ray Emission

In this section, we simulate the effect of a forward interplanetary shock passing through the Earth with a virtual camera at two special positions a and b. During this process, six stages are chosen (shown by the dashed vertical lines in Figure 5). All images are simulated at the strongest spectral line (energy 107.047 eV).

First, the initial solar wind parameters (N_SW=6 cm^-3, V_SW=700 km s^-1, T_SW=10 eV, B_T=8 nT, θ_c=0°, and $f_{O^{6+}}=0.001$, the same as those adopted in Figure 1) were input into the hybrid model and run until an equilibrium state was achieved. Then, the solar wind IP-TD structure was inserted into the model run. The IMF was turned southward (θ_c=180°) with B_T increased to 16 nT, N_SW was doubled to 12 cm^-3, and T_SW was enhanced to 15 eV. The value of $f_{{\mathrm{O}}^{6+}}$ in solar wind was kept at 0.001. During the following evolution of the solar wind-magnetosphere interaction, five stages after the IP-TD structure arrived the dayside bow shock were selected to investigate the temporal variations of the X-ray emissions.

In the first stage, when the IP-TD structure does not reach the bow shock, X-ray emission intensity is low (Figure 6a,b, similar to Figure 3c). In Stage 2, when the IP-TD structure reaches about 10 R_E upstream of the Earth and collides with the bow shock, the dayside maximum emission intensity occurs in the area where the IP-TD structure contacts the magnetosheath, i.e., subsolar region. As can be seen in Stages 3-6 in Figure 6, after IP-TD structure passing through the cusp, the cusp is always the area with the highest emission intensity. An obvious signature in Figure 6 is that the X-ray emission intensity in the subsolar and low latitude magnetosheath region is generally lower than the intensities in the flank region or the high latitude region, possibly because the ion velocities in the subsolar and low latitude region are significantly lower than those in other regions (Figure 1d-e). One can also see that a clear shock structure passed from the upstream towards the Earth and becomes not clear after State 5.

3.4. Emission Intensity Enhancement

The dayside and nightside SXR emission intensity peak values in Figure 6 are further plotted in Figure 7. During the passing of the IP-TD structure, the emission intensity of X-ray increases with time, and the maximum value achieves 27.11 keV cm^-2 s^-2 sr^-1 on the dayside, which are 4.5 times that before the IP-TD structure reaches the bow shock. This means that when the solar wind density increased by 2 times with southward IMF turning, the X-ray emission intensity can enhance by at least 4 times. Such an emission intensity enhancement is significantly smaller than the EUV emissions which could increase 40 times (38). The maximum emission intensity imaged at location b is always higher than that at location a, especially under conditions of enhanced solar wind conditions. The increase in the nightside is small, only about twice the initial. This may be because heavier solar wind ions (O⁶⁺) are difficult to be transported to the nightside magnetosphere, while He²⁺ ions are lighter and can be transported to the nightside magnetosphere [38].

3.5. Bow Shock and Magnetopause Compressions

The data of the plasma temperature on the Earth's dayside along the X_GSM axis are extracted from the hybrid simulation result (Figure 8). The corresponding color of each stage is displayed in the upper right corner of the figure. In the first and second stages (green and purple lines in the figure), the temperature rises rapidly from ~1.0 eV at 11.5 R_E, and quickly reaches ~1.0 keV at 10.9 R_E. This sharp boundary corresponds to the bow shock. Downstream of the bow shock, there is a boundary with sharp temperature decrease at ~8.0 R_E, corresponding to the magnetopause. Therefore, the positions where the temperature dramatically increases and decreases are identified as the positions of bow shock and magnetopause, respectively.

From State 3 on, when the IP-TD structure enters the magnetosheath, both the bow shock and the magnetopause are significantly compressed and move towards the Earth. We extract the position data of the bow shock and magnetopause in the six stages in Figure 8 and draw a line chart of the position changes in Figure 9. It can be clearly seen that the bow shock moves from 11.4 R_E to 9.3 R_E, and the magnetopause moves from 7.8 R_E to 6.6 R_E. This is a very strong compression of the magnetosphere. The compression is the most significant during the period when the IP-TD contact the dayside bow shock and pass through the dayside magnetopause (State 2 to Stage 3).

Figure 10 shows the simulated SXR emission profiles along the X_GSM axis from position a. The vertical solid and dashed lines show the positions of bow shock and magnetopause in Figure 8, respectively. It is shown that the position of the bow shock corresponds well with the right edge of the radial profile of SXR emission, that is, the position of the bow shock is consistent with the coordinates of the first sudden increase in SXR emission. Therefore, the SXR emission can be reliably utilized to identify the location and shape of the bow shock in future X-ray imaging. However, the magnetopause structure in SXR images may be not so clear due to the LOS overlapping effect. Other method should be investigated in future to extract the magnetopause structures in the X-ray images.

5. Conclusions

This paper uses the global hybrid simulation model, the geocoronal hydrogen model and the SASAL to simulate the spectral lines of solar wind O⁶⁺ ions to study the impact of interplanetary shock on SXR emission in the Earth's magnetosheath. The main results are summarized as follows:

• A global hybrid model is used to simulate the number density, velocity, and temperature distribution of O⁶⁺ in the magnetosheath. Then the SASAL is used to calculate the state distributions of excited O^5+* ions, which is generated in the collisions between O⁶⁺ ions and hydrogen atoms, to obtain the efficiency factors for the X-ray emissions. Then the X-ray emission intensities are simulated at 60.0 R_E at dusk side and at north pole.
• It is found that there are two closed strong emission lines for O^5+* at 10.582 nm and 10.583 nm, corresponding to energy of 107.047 eV and 107.040 eV. The efficiency factors of both lines vary significantly with collision speed. All energy states of O^5+* ions have almost the same distribution pattern in the simulate images. The strongest emission intensity on the day side is generated in the cusp region, and the strongest emission intensity on the night side is in the area near the magnetopause in the magnetosheath.
• During the passage of an IP-TD structure, significant emission intensity enhancement and compressions of both bow shock and magnetopause are observed in the simulations. The passage of the shock in the magnetosheath is clearly shown in the simulated images. The emission intensities simulated at dusk are always stronger than those simulated at north pole, possibly because that the solar wind heavy ions in the magnetosheath are mainly confined near the equatorial region.
• The emission intensity profiles can reliably be used to extract the location of bow shock, while the location of the magnetopause is difficult to be identified possibly due to the LOS overlapping effect.

During the passing of solar wind structures (such as CME, discontinuity, IP shock, etc.) in the magnetosphere, the influence in the magnetosheath can be clearly observed in the X-ray images, such as the earthward motions of bow shock and magnetopause. Thus, imagers with large field of view (e.g., the SMILE SXI and other EUV imagers in future) will be powerful to monitor the large-scale dynamic evolution of the magnetosheath and reveal the mass and energy transportation from solar wind into the magnetosphere. If the imagers are capable of distinguishing X-ray spectral lines, dynamics of specific heavy ions in the magnetosheath can be traced. The emission properties for other solar wind heavy ions will be investigated in future.