# Impact of the transport of magnetospheric electrons on the composition of the Triton atmosphere

B Benne, B Benmahi, M Dobrijevic, T Cavalié, J-C Loison, K M Hickson, M

Barthélémy, J Lilensten

► To cite this version:

B Benne, B Benmahi, M Dobrijevic, T Cavalié, J-C Loison, et al.. Impact of the transport of magnetospheric electrons on the composition of the Triton atmosphere. Astronomy & Astrophysics - A&A, 2024, 686, 10.1051/0004-6361/202346699. hal-05357917

HAL Id: hal-05357917

https://hal.science/hal-05357917v1

Submitted on 10 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Impact of the transport of magnetospheric electrons on the composition of the Triton atmosphere

B. Benne<fnref n="1" />,<fnref n="2" />, B. Benmahi<fnref n="3" />, M. Dobrijevic<fnref n="1" />, T. Cavalié<fnref n="1" />,<fnref n="4" />, J.-C. Loison<fnref n="5" />, K. M. Hickson<fnref n="5" />, M. Barthélémy<fnref n="6" />, and J. Lilensten<fnref n="6" />

Received 19 April 2023 / Accepted 23 February 2024

## ABSTRACT

**Context.** Due to its inclined orbit and the complex geometry of the magnetic field of Neptune, Triton experiences a highly variable magnetic environment. As precipitation of magnetospheric electrons is thought to have a large impact on the Triton atmosphere, a better understanding of the interaction between its atmosphere and the magnetosphere of Neptune is important.

**Aims.** We aim to couple a model of the Triton atmosphere with an electron transport model to compute the impact of a varying electron precipitation on the atmospheric composition.

**Methods.** We coupled a recent photochemical model of the Triton atmosphere with the electron transport model TRANSPlanets. The inputs of this code were determined from Voyager 2 observations and previous studies. The main inputs were the electron precipitation flux, the orbital scaling factor, and the magnetic field strength. The electron-impact ionization and electron-impact dissociation rates computed by TRANSPlanets were then used in the photochemical model. We also analyzed the model uncertainties.

**Results.** The coupling of the two models enabled us to find an electron density profile, as well as N₂ and N number densities, that are consistent with the Voyager 2 observations. We found that photoionization and electron-impact ionization are of the same order, in contrast to the results of previous photochemical models. However, we emphasize that this result depends on the hypotheses we used to determine the input variables of TRANSPlanets. Our model would greatly benefit from new measurements of the magnetic environment of Triton, as well as of the electron fluxes in the Neptune magnetosphere.

**Key words.** astrochemistry – planets and satellites: atmospheres – planets and satellites: magnetic fields

## 1. Introduction

Since the flyby of Triton by Voyager 2 in 1989, we have known that the Triton ionosphere is surprisingly more dense than that of Titan. As Triton is three times farther from the Sun than Titan, it is difficult to explain this density. Thus, an additional ionization source seems to be required. This supplementary source was quickly thought to be electrons precipitating from the Neptune magnetosphere (Majeed et al. 1990; Strobel et al. 1990), as energetic particles were observed by Voyager 2 in the Neptune magnetosphere (Krimigis et al. 1989). In addition, these precipitating electrons could bring the additional energy needed to explain the observed thermospheric temperatures (Stevens et al. 1992; Krasnopolsky et al. 1993). This precipitation is considered in the majority of Triton atmospheric models (Summers & Strobel 1991; Krasnopolsky et al. 1993; Krasnopolsky & Cruikshank 1995; Strobel & Summers 1995; Strobel & Zhu 2017; Benne et al. 2022; Krasnopolsky 2023), with the exception of Lyons et al. (1992). The latter model did not consider electron precipitation, but was able to reproduce the electron profile by using an unrealistically high reaction rate for the charge exchange between N₂⁺ and C.

The magnetic environment of Triton is highly variable due to the combination of its inclined orbit and the complex configuration of the magnetic field of Neptune. As Voyager 2 only passed through the Neptune magnetosphere once, we lack crucial data to precisely describe this environment. Strobel et al. (1990) used measurements performed by the Plasma Science System (PLS) and by the Low-Energy Charged Particles (LECP) instruments in the vicinity of Triton to determine the electron precipitation flux and compute the electron-impact ionization profile and the power deposited in the Triton atmosphere. To take into account the variability of the environment along the Triton orbit, they also used an orbital scaling factor of 0.25, which reduced the total amount of power deposited in the atmosphere. Summers & Strobel (1991) recommended shifting the ionization profile of Strobel et al. (1990) upward by two scale heights (≈100 km) so that the maximum electron production rate would correspond to the electron density peak measured by Voyager 2. This modification was then used in the models of Stevens et al. (1992), Krasnopolsky et al. (1993), Krasnopolsky & Cruikshank (1995), Strobel & Summers (1995), and Benne et al. (2022). This shift of the ionization profile is equivalent to considering that high-energy electrons (E > 50 keV) do not precipitate in the Triton

<footnote n="1">¹ Laboratoire d'Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="2">² The University of Edinburgh, School of Geosciences, Edinburgh, UK
e-mail: benjamin.benne@ed.ac.uk</footnote>

<footnote n="3">³ Laboratoire de Physique Atmosphérique et Planétaire, B5C Quartier Agora, Allée du Six Août 19c, 4000 Liège, Belgium</footnote>

<footnote n="4">⁴ LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France</footnote>

<footnote n="5">⁵ Institut des Sciences Moléculaires, CNRS, Univ. Bordeaux, 351 Cours de la Libération, 33400 Talence, France</footnote>

<footnote n="6">⁶ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>
<table><caption>**Table 1.** Comparison between the different OTD models used to approximate the magnetic field of Neptune.</caption><thead><tr><th>Model</th><th>Dipole moment<br>(T.R<sub>N</sub><sup>3</sup>)</th><th>x<sub>c</sub><br>(R<sub>N</sub>)</th><th>y<sub>c</sub><br>(R<sub>N</sub>)</th><th>z<sub>c</sub><br>(R<sub>N</sub>)</th><th>Tilt angle<br>(°)</th><th>Direction<br>(°W)</th></tr></thead><tbody><tr><td>OTD1<sup>(1)</sup></td><td>1.33 × 10<sup>-5</sup></td><td>0.17</td><td>0.46</td><td>-0.24</td><td>46.8</td><td>79.5</td></tr><tr><td>OTD2<sup>(2)</sup></td><td>1.30 × 10<sup>-5</sup></td><td>0.19</td><td>0.48</td><td>-0.19</td><td>45.2</td><td>76.5</td></tr><tr><td>OTD-O8<sup>(3)</sup></td><td>1.42 × 10<sup>-5</sup></td><td>0.05</td><td>0.48</td><td>0.00</td><td>46.9</td><td>72.0</td></tr></tbody></table>

**Notes.** ($x_c$, $y_c$, $z_c$) are the coordinates of the center of the dipole, using the coordinate system defined by the Voyager Project Steering Group (Connerney et al. 1991). This coordinate system differs from that of the International Astronomical Union (see for example Archinal et al. 2018). The tilt angle corresponds to the inclination of the magnetic dipole axis with respect to the Neptune rotation axis, and the direction indicates toward which longitude the dipole axis is inclined. The west longitude is defined as in Connerney et al. (1991).

**References.** (1) Ness et al. (1989); (2) Ness et al. (1995); (3) Connerney et al. (1991).

atmosphere (Stevens et al. 1992). However, this hypothesis was questioned by Sittler & Hartle (1996), whose results suggest that high-energy electrons are the most likely to precipitate in the Triton atmosphere. In their own model, Krasnopolsky et al. (1993) considered the orbital scaling factor as a free parameter. They chose to use a factor of 0.162 in their nominal model, with which they multiplied the ionization profile from Strobel et al. (1990) to fit the electron density observations from Voyager 2. Therefore, ad hoc modifications were made in previous models of the Triton atmosphere to model the electron precipitation and find electron density profiles consistent with the Voyager 2 observations.

Benne et al. (2022) found that electron-impact ionization was the main source of ionization in the Triton atmosphere, with a ratio of the photoionization to electron-impact ionization of 3/8. However, their maximum electron density was 2.5–5 times higher than the measured peaks presented in Tyler et al. (1989), even when chemical uncertainties were considered. Thus, it appeared necessary to model the precipitation of magnetospheric electrons better in order to improve these results. This is the aim of the present work: we couple the photochemical model of Benne et al. (2022) with the electron transport model TRANS-Planets. This model is derived from the family of the TRANS* models that were used to compute electron transport in various planetary atmospheres (Gronoff et al. 2009a). In this paper, we adapt it to the Triton conditions. Its main input parameters (electron precipitation flux, magnetic field strength, and orbital scaling factor) are computed from a model of the Neptune-Triton system and from previous studies (Strobel et al. 1990; Sittler & Hartle 1996). TRANSPlanets then computes the electron-impact ionization and dissociation rates that are used in the photochemical model of the Triton atmosphere. We present how we treated the highly variable magnetic environment of Triton in order to determine the inputs of the electron transport model in Sect. 2. We detail TRANSPlanets, its inputs, and how we coupled it to our photochemical model in Sect. 3. We show the secondary ionization profiles in Sect. 4 and our results, both nominal and considering chemical uncertainties, in Sect. 5, and we also discuss them. We present a sensitivity study of the electron temperature in Sect. 6 and finally give our conclusions in Sect. 7.

## 2. Variable magnetic environment of Triton

The Neptune magnetic field is very complex, as it cannot be represented at close distances without considering quadrupole and octupole terms (Ness et al. 1989, 1995; Connerney et al. 1991). When it is approximated with a dipole (this approximation is valid from 4 to 15 Neptunian radii ($R_N$) (Ness et al. 1995),

with $R_N = 24765$ km), its center is offset from the center of the planet, and the dipole axis is not aligned with the rotation axis of the planet. Thus, the simple models of the magnetic field of Neptune are offset tilted dipole (OTD) models. The first model of this type was presented in Ness et al. (1989) and is referred to as OTD1 throughout this paper. An updated version of this model was provided in Ness et al. (1995) and is referred to as OTD2 in the following. Finally, Connerney et al. (1991) reported a more complex model of the magnetic field of Neptune with quadrupole and octupole moments. Because these terms are only predominant at close distances from Neptune, however, they also give a dipolar approximation of their model for computations at larger distances. We refer to this model as OTD-O8. The parameters of the three models are given in Table 1, which is taken from Ness et al. (1995).

The offset-tilted nature of this magnetic dipole, coupled to the inclination of the Triton orbit and to the fast rotation of the planet (16.1 h) in comparison to the orbital period of Triton (5.88 Earth days), causes the satellite distance to the center of the dipole, as well as its magnetic latitude, to vary strongly. This can be observed when computing the L-shell parameter of Triton. This parameter corresponds to the distance at which a given field line of a dipolar field crosses the magnetic equator (Ness et al. 1989). In our case, this corresponds to the L-shell of the line going through the position of Triton. For a dipole field, this parameter is obtained for the following formula:

$$L(\lambda, r) = \frac{r}{\cos^2(\lambda)}, \qquad (1)$$

where $r$ is the distance between Triton and the center of the dipole in $R_N$, and $\lambda$ is the magnetic latitude.

To compute the different parameters needed to describe the magnetic environment of Triton, we developed a Python code. This code allows us to derive the position of Triton with respect to Neptune and to the magnetic dipole, using the ephemeris from the JPL Horizons System<fnref n="1" />. We found that $\lambda$ ranges from roughly -70° to +70°, which causes $L$ to vary from 14 to more than 100 for every OTD model (128 for the OTD1 and OTD-O8 models, and 109 for the OTD2 model), as shown in panel a of Fig. 1. As the inner magnetosphere of Neptune seems to be delimited by the minimum L-shell of Triton (Mauk et al. 1995) and as the electron number density is only assumed to be important near the magnetic equator (Strobel et al. 1990), the electron precipitation into the Triton atmosphere may vary strongly along its orbit.

In panel b of Fig. 1, Triton remains within the distance interval in which the dipolar approximation is valid. Thus, knowing the magnetic latitude of Triton, the magnetic moment of the

<footnote n="1">¹ https://ssd.jpl.nasa.gov/horizons/app.html#/</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. Examples of the variation in some parameters used to describe the magnetic environment of Triton during one orbit around Neptune. These results are obtained with the OTD-O8 model from Connerney et al. (1991), starting from August 21, 1989. (a) Variation in the L-shell parameter of Triton. (b) Variation in the distance between Triton and the center of the dipole.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 2. Variation in the magnetic field at Triton. These results were obtained with the OTD-O8 model from Connerney et al. (1991). (a) Variation in the magnetic field norm over one orbit around Neptune, starting from August 21, 1989. (b) Histogram of the value of the norm of the magnetic field at Triton when the Triton L-shell is lower than 15.5. These results were computed over one period of variation in the magnetic field at Triton, i.e., 35 days, from August 21, 1989.</figcaption></figure>

dipole, and the distance between its center and the satellite, we can compute the norm of the magnetic field as

$$ \|\vec{B}\| = \frac{M}{r^3} \times \sqrt{1 + 3 \sin^2(\lambda)}, \qquad (2) $$

with M as the magnetic moment in $T \cdot R_N^3$ (r is then expressed in $R_N$). When we use the same time period as in Fig. 1, we obtain the magnetic field at Triton presented in panel a of Fig. 2, which shows a factor of about two between extrema.

### 3. The TRANSPlanets model

TRANS is an electron transport model that was first developed to model the effects of electron precipitation into the Earth atmosphere (Lilensten 1989). It was then adapted to a variety of Solar System bodies such as Venus (Gronoff et al. 2007, 2008), Mars (Witasse et al. 2002, 2003; Simon et al. 2009; Nicholson et al. 2009), Jupiter (Menager 2011), or Titan (Lilensten et al. 2005a,b; Gronoff et al. 2009a,b, 2011).

The TRANS model computes the production of ions and electrons following the propagation of a flux of suprathelial electrons in the atmosphere. This flux is composed of photoelectrons produced by the photoionization of atmospheric species by solar extreme-UV (EUV) photons and of precipitating magnetospheric electrons. This primary ionization (i.e., photoionization) is also computed in TRANS. The code can be used with photoelectrons only, magnetospheric electrons only, or with both, as in this work. Then, the code computes the secondary ionization (ionization by the suprathelial electrons from primary ionization and electron precipitation) and dissociation through electron-impact reactions between this propagating flux and the atmospheric species. We describe the flow of operations of the code in Fig. 3. The secondary production rates of the ions and electrons computed by TRANS are of interest in our case, as they can be used in our photochemical model. We consequently modified the original code used for Jupiter from Menager (2011) into a generic Fortran 90 version, which we call TRANSPlanets. This can be used for any planetary body and was adapted here to Triton.
<figure><img src="image_4.png" /><figcaption>Fig. 3. TRANSPlanets model workflow. The inputs are shown in the yellow boxes, the computational parts in the red boxes, the intermediate variables in the blue box, and the final outputs in the green box. ME stands for magnetospheric electrons.</figcaption></figure>

In this section, we describe the fundamental equations used in the code, the inputs we used for Triton, and how we coupled it to our photochemical model. We point out that complete descriptions of the fundamental equations used in TRANS and TRANSPlanets have been presented by Lilensten (1989), Simon (2006), Gronoff (2009), Menager (2011), and Benmahi (2022).

## 3.1. Fundamental equations

To compute the secondary ionization of the neutral species following the propagation of a suprathermal electron flux, TRANSPlanets solves the nonconservative Boltzmann equation, which describes the variation in the electron distribution function $f(\vec{r}, \vec{v}, t)$ in the phase space (Stamnes & Rees 1983b; Karamcheti & Sentman 1968),

$$ \frac{\partial f}{\partial t} + \vec{v} \cdot \frac{\partial f}{\partial \vec{r}} + \frac{\partial}{\partial \vec{v}} \left( \frac{f \vec{X}}{m_e} \right) = S, \quad (3) $$

where $\vec{r}$ is the position of the electrons, and $\vec{v}$ is their speed at time $t$. $\vec{X}$ represents the friction force between the suprathermal electrons and the thermalized electrons that are already part of the atmosphere, $m_e$ is the electron mass, and $S(\vec{r}, \vec{v}, t)$ is the source term describing the electron production. We express the friction force as in Menager (2011),

$$ \vec{X} = -n_e L(E) \frac{\vec{v}}{v}, \quad (4) $$

with

$$ L(E) = \frac{3.37 \times 10^{-12}}{E^{0.94} n_e^{0.03}} \left( \frac{E - E_{th}}{E - 0.53 E_{th}} \right)^{2.36}, \quad (5) $$

where $E$ is the energy of the suprathermal electron, $n_e$ is the number density of the thermal electrons, and $E_{th}$ is their energy. All the energies are expressed in eV. $L(E)$ is the friction function expressed in eV cm$^2$ as given in Swartz et al. (1971) and recommended by Stamnes & Rees (1983a).

Equation (3) can be expressed depending on the suprathermal electron flux $\Phi$ instead of on their distribution function by making the following variable change:

$$ \Phi(\vec{r}, E, \vec{u}, t) = \frac{v^2}{m_e} f(\vec{r}, \vec{v}, t), \quad (6) $$

with $E$ the kinetic energy of the electrons and $\vec{u}$ their direction of propagation, such that $\vec{u} = \frac{\vec{v}}{v}$, giving

$$ \frac{1}{v} \frac{\partial \Phi}{\partial t} + \frac{\vec{v}}{v} \cdot \frac{\partial \Phi}{\partial \vec{r}} - n_e \frac{\partial [L(E)\Phi]}{\partial E} = \frac{1}{v} \hat{S} \quad (7) $$

with $\hat{S}(\vec{r}, E, \vec{u}, t) = \frac{v^2}{m_e} S(\vec{r}, \vec{v}, t)$.

Considering steady state and a plane-parallel atmosphere, we obtain the equation describing the propagation of a flux $\Phi$ of suprathermal electrons along a magnetic field line projected on the vertical axis,

$$ \mu \frac{\partial \Phi(\tau, \mu, E)}{\partial \tau(z, E)} = -\Phi(\tau, \mu, E) + \frac{n_e(z)}{\sum_k n_k(z) \sigma_k^{\text{tot}}(E)} \frac{\partial [L(E)\Phi(\tau, \mu, E)]}{\partial E} + D(z, \mu, E) + P(z, \mu, E), \quad (8) $$

where $\mu$ is the cosine of the pitch angle, $\tau$ is the collision depth, which is similar to an optical depth in radiative transfer, $z$ is the altitude of the considered level, $n_k(z)$ is the number density of species $k$ at this level, and $\sigma_k^{\text{tot}}(E)$ is the total electron-impact cross section of species $k$ at energy $E$. The collision depth is computed as

$$ \tau(z, E) = \int_{z}^{z_{\text{max}}} d\tau \quad \text{with} \quad d\tau = \sum_{k} n_{k}(z) \sigma_{k}^{\text{tot}}(E) dz. \quad (9) $$

The suprathermal electrons center of mass propagates along the magnetic field line. The projection on the vertical axis is made by simply dividing by the cosine of the magnetic dip angle. This is incorrect for large angles (typically larger than 45°) because the cosine tends toward zero. Rigorously, the most accurate results are obtained when the local magnetic field line coincides with the vertical dimension.

The source term is decomposed into two terms $P$ and $D$. $P$ is the electron production term, corresponding to the sum of the electron production from photoionization and the precipitation flux. $D$ is the diffusion term, resulting from the diffusion of suprathermal electrons to lower energies through collisions with atmospheric species. The code is multistream, meaning that $\mu$ is discretized on a certain number of angles chosen by the user. This allows computing the different output fluxes on various angles and not just the downward and upward fluxes. In our case, we chose to run our calculations over 16 angles over 360° because this discretization provided the best compromise in terms of accuracy and computational time.

## 3.2. TRANSPlanets inputs

As shown in Eq. (8), TRANSPlanets needs some additional data to solve the Boltzmann equation. It needs the electron-impact ionization and dissociation cross sections for the considered atmospheric species, an input flux, an energy range, an energy grid, an altitude grid, the number density profiles of the considered species and of the thermal electrons, the temperature
profiles for these species, and a solar flux with the associated photolysis cross sections. In the following sections, we first describe the inputs that are not related to the photochemical model. These inputs are detailed in Sect. 3.4.

### 3.2.1. Electron-impact cross sections and energy grid

In our model, the electron-impact cross sections come from different sources, depending on the species. The elastic and inelastic (electron-impact ionization and dissociation) cross sections of N₂ were taken from the Atomic and Molecular Cross section for Ionization and Aurora Database (AtMoCiad; Gronoff et al. 2021). Those of CO were taken from Itikawa (2015). For C, we took the elastic cross sections from the NIST database, while the electron-impact ionization cross section was taken from Jonauskas (2018) and the electron-impact excitation cross sections from Suno & Kato (2006). For the N(⁴S) and N(²D) cross sections, we used inelastic cross sections from Jonauskas (2022) and Kato (1994). The elastic cross sections below 50 eV were taken from Wang et al. (2014), and they were completed above this energy with data from the NIST database. As the maximum energy used in the elastic cross section of N₂ in the AtMoCiad database is 10⁵ eV, we used this energy as the upper limit of our energy range. The lower boundary was set to 1 eV to take into account rotational and vibrational excitation of N₂. The energy grid consisted of 600 points spaced logarithmically between these two energies, such that $E(i) = \exp[E(i-1) + \frac{\ln(E_{max})-\ln(E_{min})}{600}]$.

### 3.2.2. Variable parameters

The L-shell of Triton varies strongly with time, as shown in Fig. 1. Krimigis et al. (1989) have shown that the electron flux decreases strongly outside of the minimum L-shell of Triton, this value presumably defining the limit of the inner magnetosphere of Neptune (Mauk et al. 1995). This led Strobel et al. (1990) to consider that Triton experiences constant electron precipitation if L(Triton) < 15.5, with no precipitation for higher L values. With this assumption, they computed an average precipitation flux of 25% of the flux at the minimum L-shell of Triton. Thus, their orbital scaling factor was $m_{orb} = 0.25$.

In their model, Krasnopolsky et al. (1993) considered $m_{orb}$ as a free parameter. Their best value to match the electron density profiles from Voyager 2 was $m_{orb} = 0.162$. This value was used in their nominal model and then in Krasnopolsky & Cruikshank (1995).

With our modeling of the Neptune-Triton system, we are able to compute the value of the Triton L-shell at any moment. We could then have used a relation between the value of L and the magnetospheric electron flux, but no such relation exists due to the scarcity of the Voyager 2 electron flux measurements. Thus, we decided to keep the criterion used in Strobel et al. (1990) and considered that electron precipitation only occurs when the L-shell of Triton is lower than 15.5. Therefore, the orbital scaling factor corresponds to the ratio of the time when L is lower than 15.5 to the time interval considered. We consistently find for all three magnetospheric models presented in Table 1 that $L < 15.5$ for 27–28.5% of the time. Therefore, we used $m_{orb} = 0.27$.

The value of the magnetic field strength B at Triton varies significantly with time, as shown in panel a of Fig. 2. When we plot the histogram of the value of $\|\vec{B}\|$ when $L < 15.5$, we obtain panel b of Fig. 2 (with the OTD-O8 model). The mean value of the magnetic field norm is $(5.07 \pm 0.35) \text{ nT}$. This value

is consistent with Fig. 3 of Ness et al. (1989) and with the value found with Eq. (2) using the values of r and $\lambda$ from Fig. 7 of Mauk et al. (1995). This value is lower than the value used in Strobel et al. (1990) and Sittler & Hartle (1996), which is 8 nT. This difference impacts the power deposited in the atmosphere by curvature drift, as shown in Sect. 3.2.3. When we assume that the electrons deposit their energy by curvature drift (as in Strobel et al. 1990), using 5 nT instead of 8 nT leads to a higher deposited energy (see Eqs. (13) and (14) in Sect. 3.2.3). In our model, we used $B = 5 \text{ nT}$.

### 3.2.3. The electron precipitation flux

In this section, we determine the electron precipitation flux, that is, the flux that reaches the atmosphere of Triton and therefore interacts with atmospheric species. This flux is different from the electron flux that was measured in the Neptune magnetosphere with the LECP instrument of Voyager 2 (Krimigis et al. 1989) (called magnetospheric flux in the following). The main problem here is that the transition from the magnetospheric flux to the precipitation flux is unknown. The total energy carried by the magnetospheric electrons with an energy higher than 20 keV in the Neptune magnetosphere is 0.42 erg cm⁻² s⁻¹ (Krimigis et al. 1989; Sittler & Hartle 1996). This flux is more than 250 times higher than the flux needed to explain the Triton thermospheric temperature, which is $1.6 \times 10^{-3} \text{ erg cm}^{-2} \text{ s}^{-1}$ (Broadfoot et al. 1989). Moreover, the magnetospheric flux was measured at $12 R_N$ from Neptune (Krimigis et al. 1989), which is not in the direct vicinity of Triton.

In our work, we determined the electron precipitation flux reaching the Triton atmosphere from the work of Strobel et al. (1990) and Sittler & Hartle (1996). We recall, however, that the steps we took to compute the precipitation flux from the magnetospheric flux may not result in a flux that is representative of the real precipitation flux at Triton. Only an orbiter measuring the magnetic environment of Triton over many orbits can reveal the precipitation flux at Triton. Here, we tried to determine the flux based on the few observations we have, physics-based considerations, and following what was assumed in previous studies.

The magnetospheric flux was taken from Strobel et al. (1990), based on LECP and PLS measurements from Voyager 2 near the Triton flyby. It corresponds to a flux $j(E)$ (in cm⁻² s⁻¹ sr⁻¹ keV⁻¹) of

$$j(E) = 2.10^4 \times \left( \frac{E}{28 \text{ keV}} \right)^{-2.7} \quad \text{for } E \ge 28 \text{ keV} \qquad (10)$$

$$j(E) = 2.10^4 \times \left( \frac{E}{28 \text{ keV}} \right)^{-0.7} \quad \text{for } E < 28 \text{ keV}, \qquad (11)$$

where E is the energy of the electrons in keV.

According to Sittler & Hartle (1996), not all the electrons contained in the magnetospheric flux reach the Triton atmosphere. Their results show that the higher the electron energy, the more likely the precipitation. As stated in their article, this contradicts Summers & Strobel (1991) (and therefore Krasnopolsky et al. 1993; Krasnopolsky & Cruikshank 1995; Strobel & Summers 1995; Benne et al. 2022), who inferred from the shift in the ionization profile of Strobel et al. (1990) toward higher altitudes that the precipitation of high-energy electrons is inhibited. More precisely, Sittler & Hartle (1996) computed that fewer than 3% of the electrons with energies lower than 5 keV precipitate, fewer than 15% of electrons with energies between 5 and 50 keV,
<figure><img src="image_5.png" /><figcaption>Fig. 4. Comparison of the magnetospheric electron flux $j(E)$ from Strobel et al. (1990) (dotted blue line), derived from Voyager 2 observations, and the modified flux $j^*(E)$ (dash-dotted red line), computed following the recommendations of Sittler & Hartle (1996). The flux $j^*(E)$ is computed from the flux $j(E)$ multiplied by a precipitation probability $p(E)$.</figcaption></figure>

and half of the higher-energy electrons. Therefore, we modified the magnetospheric flux $j(E)$ from Strobel et al. (1990) to account for these probabilities. We thus computed a flux $j^*(E)$, such that

$$j^*(E) = j(E) \times p(E), \qquad (12)$$

where $p(E)$ is the precipitation probability. We considered that 3 and 15% of the electrons with energies of 5 and 50 keV precipitate, respectively, and 50% of the electrons with the highest energy considered in TRANSPlanets, that is, 100 keV. To compute the precipitation probability $p(E)$ for the electrons with intermediate energies, we fit a second-order polynomial using these three values. The $j(E)$ and $j^*(E)$ fluxes are plotted in Fig. 4.

Knowing $j^*$, we can estimate the power that is deposited by the electrons in the Triton atmosphere. Strobel et al. (1990) claimed that the power is deposited by curvature drift. In this case, the deposited power $P_{CD}(E)$ is computed with formulas (13) and (14), taken from Strobel et al. (1990),

$$P_{CD}(E) = j^*(E) \times \pi \times R_i \times \delta(E), \qquad (13)$$

where $R_i$ is the radius of the ionopause, and $\delta$ is the inward drift per pass through the ionosphere, such as

$$\delta = \frac{cm_e v_{||} \pi}{eB}, \qquad (14)$$

with $c$ the speed of light, $m_e$ the electron mass, $v_{||}$ the electron speed parallel to the magnetic field, $B$ the strength of this field, and $e$ the elementary charge.

Sittler & Hartle (1996) defined the radius of the ionopause as $1/\sigma n \approx R_i$, with $\sigma$ the ion-neutral cross section, and $n$ the neutral number density. In our case, this gives $R_i \approx 2300$ km, which nearly corresponds to their value of $\approx 2350$ km, while Strobel et al. (1990) took $\approx 2000$ km.

Using the flux $j$ in Eq. (13), considering $B = 8$ nT in Eq. (14), and using the flux between $E_{min} = 200$ eV and $E_{max} = 10^5$ eV, we find $P(E) = 8 \times 10^8$ W. This power can be compared to the total power carried by the electron flux, $E_{flux}$, which is computed

<figure><img src="image_6.png" /><figcaption>Fig. 5. Comparison of the electron flux $j^*(E)$ (dash-dotted red line) derived from the flux $j(E)$ from Strobel et al. (1990) and Eq. (12) following the results of Sittler & Hartle (1996), and the renormalized flux $j_{CD}^*(E)$ (solid orange line). This flux is computed from $j^*(E)$ with Eq. (16), so that the electrons deposit the correct amount of power in the Triton atmosphere, with the hypothesis that this power is deposited by curvature drift.</figcaption></figure>

as

$$P_{\text{flux}} = 4\pi R_i^2 \times 1.602176634 \times 10^{-19} \times \int_{E_{\text{min}}}^{E_{\text{max}}} j(E)EdE \ [\text{W}]. \qquad (15)$$

With $E_{min} = 200$ eV and $E_{max} = 10^5$ eV, we find $P_{flux} = 3, 4 \times 10^{10}$ W. Therefore, the power deposited by curvature drift is lower than the total power carried by the electron flux $j(E)$.

The same can be found for $j^*$: the total power carried by the electrons is $3.8 \times 10^9$ W, whereas the power that should be deposited by curvature drift is only $1.8 \times 10^8$ W (using $B = 5$ nT, as computed in Sect. 2). Thus, to deposit the correct amount of power in the atmosphere, that is, to mimic a power deposition by curvature drift, we chose to renormalize the precipitation flux using Eq. (16),

$$j_{CD}^*(E) = j(E) \times p(E) \times \frac{P_{CD}(E)}{P_{\text{flux}}(E)}. \qquad (16)$$

The flux $j_{CD}^*(E)$ is compared to the flux $j^*(E)$ in Fig. 5. Ultimately, the power carried by the flux $j_{CD}^*$ is $2.7 \times 10^{-3}$ erg cm$^{-2}$ s$^{-1}$, which is more consistent with the value of $1.6 \times 10^{-3}$ erg cm$^{-2}$ s$^{-1}$ computed by Broadfoot et al. (1989) to explain the thermospheric temperature.

We clarify that we performed a test using the nonrenormalized flux $j^*$ in TRANSPlanets. In this case, all the power carried by the flux was classically deposited in the atmosphere, and therefore $3.8 \times 10^9$ W were deposited instead of $1.8 \times 10^8$ W. This resulted in high electron-impact ionization rates (see Fig. A.1), electron densities, and N number densities. These densities are inconsistent with Voyager 2 observations, as shown in Fig. A.2. This demonstrates the need to renormalize the precipitation flux so that the deposited power matches the deposit by curvature drift.

Finally, we took the orbital scaling factor into account. Strobel et al. (1990) multiplied the power deposited by curvature drift by $m_{orb}$ to allow for the intermittence of the precipitation. Strobel et al. (1990) did not indicate whether the orbital scaling factor needed to be applied to the ionization profile, but this
was done by Krasnopolsky & Cruikshank (1995). In their model, they multiplied the ionization profile computed by Strobel et al. (1990) by their adjusted orbital scaling factor of 0.162 to take into account the variability of the electron precipitation over one orbit. In a similar way, we chose to multiply $j_{\text{CD}}^*$ by $m_{\text{orb}}$ in our model. We used $m_{\text{orb}} = 0.27$, as computed in Sect. 3.2.2, with our model of the Neptune-Triton system.

### 3.3. Caveats

Some assumptions and simplifications in our work have a significant impact on the results. They are listed below.

– The magnetic field line is considered to be vertical. As mentioned in Sect. 3.1, this configuration is the most accurate when considering the projection of the axis parallel to the magnetic field line on the vertical axis. However, other configurations with inclined lines would be easily and accurately solved as well. Moreover, this allows modeling mean conditions, as was done for Titan in Dobrijevic et al. (2016). However, as considered in Strobel et al. (1990) and Sittler & Hartle (1996), the magnetic field lines are likely to drape themselves around the Triton atmosphere. Inclined lines change the thickness of the atmosphere that is crossed by the electrons, thus modifying the altitude at which they deposit their energy. Here, we lack knowledge of the magnetic environment of Triton. A systematic study of any field configuration is therefore meaningless, and we have to reduce our approach, keeping in mind that this is a severe limitation.

– As shown in the previous section, we used a magnetospheric electron flux that was measured at some distance from Triton to compute the precipitation flux used in TRANSPlanets. We also renormalized this flux to deposit the correct amount of power, considering that this power is deposited by curvature drift. Therefore, our input flux is dependent on the initial measurements, and we have no guarantees that this flux is representative of the real precipitation flux at the top of the Triton atmosphere. Until new in situ measurements of both the magnetospheric and ionospheric electron fluxes are made, however, we are limited by this assumption in our model. This can only be achieved with an orbiter performing several flybys of Triton for different magnetospheric configurations, and/or a probe/lander sent into the Triton atmosphere.

– We limited our calculations to 100 keV, despite the measurements of higher-energy electrons in the Neptune magnetosphere and the results of Sittler & Hartle (1996) showing that higher-energy electrons are more likely to precipitate in the Triton atmosphere. We made this choice because TRANSPlanets solves the nonrelativistic Boltzmann equation and because our electron-impact cross sections for $N_2$ are only known until 100 keV. There is currently no globally accepted way of taking into account relativistic electrons in electron transport models. Some solutions are to use the relativistic Boltzmann equation or relativistic cross sections. In this work, we chose to stay at relatively low energies to avoid modeling problems that are related to relativistic effects.

– Liuzzo et al. (2021) modeled the interaction of Triton with the Neptune magnetosphere and predicted it to be complex. According to this model, the incident electron flux impacting the Triton ionosphere may strongly vary spatially and temporally. This variability is not taken into account in our model, as we chose to develop our approach based on Voyager 2 data, as was done in other photochemical models of the Triton atmosphere (Strobel et al. 1990; Strobel & Summers 1995; Krasnopolsky & Cruikshank 1995; Benne et al. 2022; Krasnopolsky 2023).

### 3.4. Coupling TRANSPlanets with the photochemical model

The remainder of the TRANSPlanets inputs were provided by our photochemical model, such as the atmospheric number density profiles of the neutral species and the thermal electrons. The altitude grid of the photochemical model was also used in TRANSPlanets. This grid was sampled with $H/5$ steps, where $H$ is the atmospheric scale height, giving 96 altitude levels between the surface and the top of the atmosphere (i.e., 1026 km). We chose to consider five species in TRANSPlanets, namely $N_2$, $N(^4S)$, $N(^2D)$, C, and CO, as this allowed us to consider the three most abundant species in most of the atmosphere. The temperature profile for the neutral species was the same as in the photochemical model and was taken from the Triton-3 model of Strobel & Zhu (2017). As in Benne et al. (2022), the electron temperature was taken to be equal to the neutral temperature. We adopted this hypothesis as no measurements of this physical parameter were made during the Voyager flyby. The electron temperature is likely to be higher than the neutral temperature, however, at least in the upper atmosphere. Using pressure balance arguments, Sittler & Hartle (1996) computed that the electron temperature at the ionopause should be between 308 and 1230 K, depending on the hypotheses used. This is significantly higher than the neutral temperature in the thermosphere in our model, which is 92 K. We ran some simulations with higher electron temperatures up to 1230 K (see Sect. 6). They produced a higher electron density because a higher electron temperature lowers the radiative and dissociative recombination rates. Therefore, to maintain a nominal electron profile consistent with Voyager 2 measurements, we would need to add another loss process for ions such as, for example, ion escape, as used in Krasnopolsky & Cruikshank (1995) and Krasnopolsky (2023). However, we show in Sect. 6 that the use of an electron temperature between 300 and 1230 K leads to electron density profiles that are not significantly different from those presented in Tyler et al. (1989) when we consider the chemical uncertainties and that the Tyler et al. profiles were given without any uncertainty. Thus, we chose to keep the electron temperature equal to the neutral temperature in our main model.

TRANSPlanets was used to compute the secondary ionization and dissociation of atmospheric species by a suprathreshold electron flux. This flux consists of both precipitating electrons and photoelectrons. Therefore, the photoionization rates of the atmospheric species by EUV photons need to be computed first because they give the photoelectron flux as a function of energy. These photoionization rates were already computed by the photochemical model, with a high-resolution solar flux and high-resolution cross sections for $N_2$ and CO photolysis. However, these calculations only give the electron and ion production rates, thus losing all the information about the energy of the photoelectrons. In order to keep the energy dependence of the photoelectron flux, we calculated the photoionization rates within TRANSPlanets, adopting the solar flux and cross sections of the photochemical model.

We used the same solar flux in both codes, corresponding to maximum solar activity, as was the case at the time of the Voyager flyby. The flux had a resolution of 1 nm between 1 and 730 nm and was taken from Thuillier et al. (2004). Of the five species considered in TRANSPlanets, the photochemical model currently only considers the photoionization of $N_2$, $N(^4S)$, and CO. Therefore, TRANSPlanets only computes the primary ionization of these three species and the photodissociations of $N_2$ and CO. The maximum wavelength absorbed by at least one of these species being 163 nm, we consequently restrained the
<table><caption>**Table 2.** Electron-impact ionization and electron-impact dissociation reactions considered in the photochemical model and in TRANSPlanets.</caption><tr><td>N<sub>2</sub> + e<sup>-</sup></td><td>→</td><td>N<sub>2</sub><sup>+</sup> + 2e<sup>-</sup></td></tr><tr><td>N<sub>2</sub> + e<sup>-</sup></td><td>→</td><td>N<sub>2</sub><sup>+</sup> + N(<sup>2</sup>D) + 2e<sup>-</sup></td></tr><tr><td>N<sub>2</sub> + e<sup>-</sup></td><td>→</td><td>N(<sup>4</sup>S) + N(<sup>2</sup>D) + e<sup>-</sup></td></tr><tr><td>N(<sup>4</sup>S) + e<sup>-</sup></td><td>→</td><td>N<sup>+</sup> + 2e<sup>-</sup></td></tr><tr><td>N(<sup>2</sup>D) + e<sup>-</sup></td><td>→</td><td>N<sup>+</sup> + 2e<sup>-</sup></td></tr><tr><td>C + e<sup>-</sup></td><td>→</td><td>C<sup>+</sup> + 2e<sup>-</sup></td></tr><tr><td>CO + e<sup>-</sup></td><td>→</td><td>CO<sup>+</sup> + 2e<sup>-</sup></td></tr><tr><td>CO + e<sup>-</sup></td><td>→</td><td>C + O(<sup>3</sup>P) + e<sup>-</sup></td></tr></table>

calculations of TRANSPlanets to the 1–163 nm interval. We also used the solar optical depths computed by the photochemical model in TRANSPlanets and a solar zenith angle of 50° in both codes, corresponding to mean global conditions at equinox<fnref n="2" />.

Finally, to be consistent between the photochemical model and TRANSPlanets, we considered the same electron-impact ionization and electron-impact dissociation reactions in both codes. These reactions are given in Table 2. The chemical scheme we used is that of Benne et al. (2022), to which we added the last five reactions of Table 2, as well as the reactions (N(<sup>4</sup>S); N(<sup>2</sup>D)) + H<sub>2</sub>CN → (NH + HCN; <sup>3</sup>CH<sub>2</sub> + N<sub>2</sub>), whose rates were taken from Hébrard et al. (2012) with updated calculations.

Our workflow was the following: We first ran the photochemical model without the electron-impact dissociation and electron-impact ionization reactions. Then, the N<sub>2</sub>, N(<sup>4</sup>S), N(<sup>2</sup>D), C, CO and electron number density profiles, the neutral and electron temperature profiles, and the solar optical depths were used as inputs for TRANSPlanets. Then, TRANSPlanets computed the electron-impact ionization and electron-impact dissociation rates, which we subsequently used back in the photochemical model to calculate the atmospheric composition. We iterated this process until steady state was reached, that is, until the variations in the mole fractions between two successive iterations were small compared to the chemical uncertainties. This corresponded to three iterations. The method of the coupling of the two codes is presented in Fig. 6.

### 3.5. Global mean conditions

We used our photochemical model to compute the steady-state composition of the Triton atmosphere. It thus ran over a long time interval, forcing us to consider global mean conditions. To verify whether global mean conditions were appropriate, we checked that the main ions of the Triton atmosphere had sufficiently long lifetimes to maintain an ionosphere between precipitation events. Based on previous modeling results (Majeed et al. 1990; Strobel et al. 1990; Strobel & Summers 1995; Krasnopolsky & Cruikshank 1995; Benne et al. 2022; Krasnopolsky 2023), the ionosphere may be mainly composed of atomic ions with slow recombination rates. We find that the lifetime of the main ion, C<sup>+</sup>, varies between 4.8×10<sup>6</sup> s and 7.5×10<sup>7</sup> s in the ionosphere. The one of N<sup>+</sup> ranges between 3.9 s at 204 km and 5.6×10<sup>5</sup> s at 1026 km, while the lifetime of H<sup>+</sup> lies between 3.1×10<sup>6</sup> s and 5.3×10<sup>7</sup> s. With our model of the Neptune-Triton system, we computed the time interval between precipitation events (see Fig. 7). The mean value

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Method used to couple our photochemical model of the Triton atmosphere and the electron transport code TRANSPlanets. R<sub>i</sub> is the ionosphere radius, ||B|| is the strength of the magnetic field at Triton, and m<sub>orb</sub> is the orbital scaling factor.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Histogram of the time interval between precipitation events computed with our model of the Neptune-Triton system. These values are computed over one period of variation in the magnetic field at Triton, i.e., 35 days, from August 21, 1989.</figcaption></figure>

is 5.2 h (1.9×10<sup>4</sup> s), and the time interval that occurs most frequently is 7.2 h (2.6×10<sup>4</sup> s). This value is consistent with the time interval between each Neptune magnetic equator crossing by Triton given in Strobel & Summers (1995) and Krasnopolsky & Cruikshank (1995), which is 7 h.

Therefore, the lifetimes of C<sup>+</sup> and H<sup>+</sup> are significantly longer than the time interval between precipitation events. The lifetime of N<sup>+</sup> is longer than 5.2 h above 660 km and longer than 7.2 h above 676 km. As a consequence, we find it unlikely that the abundances of these ions are significantly affected by the short-term variability in the magnetosphere. Global mean conditions for the photochemical model are thus valid.

After the steady-state composition of the atmosphere was determined, our model could compute the atmospheric composition under the specific conditions of the Voyager 2 flyby. However, we lack several critical inputs to do this, such as the electron precipitation flux and the magnetic field at Triton.

<footnote n="2">² Using this solar zenith angle does not produce significant differences from the case using a classic solar zenith angle of 60°, as the uncertainties on our results are large.</footnote>
<figure><img src="image_9.png" /><figcaption>Fig. 8. Secondary electron production profiles from Strobel et al. (1990) (dotted black line), Benne et al. (2022) (dash-dotted red line), and from this work (solid orange line).</figcaption></figure>

## 4. Secondary ionization profiles

In this section, we present the secondary electron production profile computed with TRANSPlanets. To do this, we used the flux $j_{CD}^*(E)$ presented in Fig. 5, an orbital scaling factor $m_{orb} = 0.27$, and a magnetic field $B = 5$ nT. As done for Titan in Dobrijevic et al. (2016), we considered a vertical precipitation, even though magnetic field lines are expected to drape themselves around Triton (Strobel et al. 1990; Sittler & Hartle 1996) (see Sect. 3.3). We emphasize that the following results were obtained considering secondary ionization and dissociation from both photoelectrons and magnetospheric electrons, which is not the case for the results presented in Strobel et al. (1990), where only the secondary ionization by magnetospheric electrons was considered.

The secondary electron production profile obtained with our model is compared to those from Benne et al. (2022) and Strobel et al. (1990) in Fig. 8.

Figure 8 shows that we obtain a profile with two peaks. The first peak, around 400 km, results from the ionization by photoelectrons. The second peak, at 71 km, is due to electron-impact ionization by magnetospheric electrons. The maximum production rate of secondary electrons is $1.82 \text{ cm}^{-3} \text{ s}^{-1}$.

In addition to this profile, TRANSPlanets allowed us to obtain the reaction rates for every reaction of electron-impact ionization and electron-impact dissociation separately. This is the main difference between this work and the model of Benne et al. (2022). In the latter work, the electron-impact ionization and electron-impact dissociation rates were computed from the ionization profile shown in Fig. 8 (dash-dotted red curve) and constant branching ratios from Fox & Victor (1988):

$$0.8 \text{ for the reaction } N_2 + e^- \rightarrow N_2^+ + 2e^-,$$

$$0.2 \text{ for } N_2 + e^- \rightarrow N^+ + N(^2D) + 2e^-$$

$$0.6 \text{ for } N_2 + e^- \rightarrow N(^4S) + N(^2D) + e^-.$$

Using TRANSPlanets, we did not need to use branching ratios because the reaction rates were computed independently for each reaction. Therefore, these profiles, one per reaction presented in Table 2, were used in the photochemical model. They are presented in Fig. 9. We describe the results from the photochemical model using these rates in the next section.

## 5. Results and discussion

In the sections that follow, we compare some of the results from this work to results from a model that did not use TRANSPlanets, similar to Benne et al. (2022). This model is a slightly modified version of the model from Benne et al. (2022), in which we improved the way in which the interstellar medium flux is taken into account. This slightly modifies the nominal results presented in Benne et al. (2022), but the changes remain small compared to the chemical uncertainties. Moreover, the key chemical reactions and the main chemical pathways remain the same. Thus, we refer to this model as Benne et al. (2022). We recall that in this model, the electron-impact dissociation and electron-impact ionization rates were not computed with TRANSPlanets, but derived from the profile of Strobel et al. (1990), as described in Sect. 4.

As in Benne et al. (2022), we performed an evaluation of the uncertainties on the model results that are caused by the uncertainties in the chemical reaction rates. The method used here is the same as was presented in Benne et al. (2022). The Monte Carlo procedure was applied using 250 sets of the chemical reaction rates.

### 5.1. Ionization sources

The main result from this work is a significant decrease in the electron density in most of the atmosphere compared to Benne et al. (2022), as shown in panel a of Fig. 10, resulting from a decrease of 66% in the electron-impact ionization rate, as listed in Table B.1. For the nominal version of the model, ionization through photoionization is slightly more important than electron-impact ionization, with a ratio of 6/5. This is in contrast to the results of Benne et al. (2022), where this ratio was 3/8, similar to 1/2 found in Krasnopolsky & Cruikshank (1995), even though we should underline that it is now close to unity. To confirm this result, we computed the photoionization to electron-impact ionization ratio for the 250 runs of the Monte Carlo procedure. The histogram of the ratios is given in panel b of Fig. 10. We observe that electron-impact ionization is the main ionization source in the Triton atmosphere in only 16% of the 250 runs of the Monte Carlo procedure, and the mean value of the photoionization to electron-impact ionization ratio is close to the nominal value of 6/5. This indicates that photoionization is slightly more efficient than electron-impact ionization in the Triton atmosphere. This result depends on the magnitude of the precipitation flux, however, and thus on the hypotheses we made to compute $j_{CD}^*(E)$. This point is discussed in Sect. 5.5.

### 5.2. Electron density

With our model, the significant decrease in the electron-impact ionization rate results in a nominal electron peak density of $3.9 \times 10^4 \text{ cm}^{-3}$ at 367 km, which agrees well with the Voyager 2 observations of a peak of $(3.5 \pm 1.0) \times 10^4 \text{ cm}^{-3}$ at (340–350) km (Tyler et al. 1989; Krasnopolsky & Cruikshank 1995). Therefore, the electron peak from this work is significantly smaller than that from Benne et al. (2022), where the peak electron density was $9.8 \times 10^4 \text{ cm}^{-3}$ at 345 km. Our results also agree better with the Voyager 2 observations than the results from Strobel & Summers (1995) and Krasnopolsky & Cruikshank (1995), who find their peaks at 275 and 320 km, with peak number densities of $3.4 \times 10^4$ and $3.5 \times 10^4 \text{ cm}^{-3}$, respectively.

The electron density profile is strongly influenced by the reaction rate of the charge-exchange reaction between $N_2^+$ and C. This reaction is crucial for $C^+$ because it contributes 81% of
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Electron-impact ionisation and electron-impact dissociation rates computed with TRANSPlanets using the precipitation flux $j_{CD}^*$ multiplied by an orbital scaling factor of 0.27. (a) Reactions involving N$_2$, N($^4$S), and N($^2$D). (b) Reactions involving C and CO.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Results from this work about the electron density and ionization processes. (a) Nominal electron density profiles from Benne et al. (2022) (dash-dotted green line, using the electron production profile presented in Fig. 8 and not TRANSPlanets) and from this work (red, using TRANSPlanets to compute the electron-impact ionization and electron-impact dissociation rates). The red triangles and blue circles give the electron density profiles measured by Voyager 2 at the winter and summer occultation points, respectively (data from Tyler et al. 1989). The gray profiles are the 250 electron profiles obtained from this work with the Monte Carlo procedure. (b) Histogram of the ratio of photoionization to electron-impact ionization obtained from the 250 runs of the Monte Carlo procedure.</figcaption></figure>

its integrated production. Since the main ion in our model is C$^+$, this reaction also impacts the electron density. The nominal rate constant of this reaction is $10^{-10} \text{ cm}^3 \text{ s}^{-1}$, taken from the KIDA database (Wakelam et al. 2012). The uncertainty factor for this reaction is $F(300 \text{ K}) = 3$. The rate of this reaction is thus significantly lower than the rate used in Lyons et al. (1992), who took $10^{-9} \text{ cm}^3 \text{ s}^{-1}$, and explained the electron profiles with photoionization alone. This reaction is also a key uncertainty reaction, which means that it significantly affects the uncertainties on the model results. It is thus necessary to measure this reaction rate at temperatures representative of the Triton atmospheric conditions to improve Triton photochemical models.

The results obtained with the Monte Carlo simulation are also consistent with the observations, as shown in panel a of Fig. 10. The electron density at 367 km, the altitude of the nominal peak, is $4.2 \times 10^4 \text{ cm}^{-3}$, with an uncertainty factor of 1.53. The altitude of the peak is $(372 \pm 8) \text{ km at } 1\sigma$, which is slightly

above the (340–350) km measured by Voyager 2 (Tyler et al. 1989). However, we note that the altitude of the electron peak depends on the inclination of the magnetic field line because it affects the atmospheric layers in which the electron-impact ionization and dissociation reactions take place. Thus, the altitude of the peak might change for a more inclined magnetic field line. For now, as shown in Fig. 8, the secondary electron production from electron-impact ionization by the precipitating electrons is maximum at 71 km. A more inclined line should bring this peak to a higher altitude, which would shift the electron peak number density, possibly toward lower altitudes. We also note that in the end, the difference in the altitude of the electron peak may not be that significant because the data presented in Tyler et al. (1989) were given without any uncertainty. Voyager 2 RSS data was reanalyzed with modern tools by Togni et al. (2023), who found that the electron peak may be located at higher altitudes than previously presented by Tyler et al. (1989).
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Nominal simulated mole fraction profiles of the main species of the Triton atmosphere. (a, b) Main neutral species. (c) Main hydrocarbons and HCN. (d) Main ions. The bars on the N₂ and N(⁴S) profiles are the measurements from the UVS instrument of Voyager 2 (Broadfoot et al. 1989; Krasnopolsky et al. 1993).</figcaption></figure>

## 5.3. Results on the neutral atmosphere

The nominal mole fraction profiles of the main species of the Triton atmosphere are given in Fig. 11. We obtain nominal N₂ and N number densities consistent with Voyager Ultraviolet Spectrometer (UVS) observations (Broadfoot et al. 1989; Krasnopolsky et al. 1993), as shown in panels a and b of Fig. 11. We find a N₂ number density of 4.11×10⁸ cm⁻³ at 575 km (with an uncertainty factor of 1.0), in agreement with (4±0.4)×10⁸ cm⁻³ from the observations. For atomic nitrogen, we can compare our results to measurements at 400 and 200 km, as given in Krasnopolsky et al. (1993); Krasnopolsky & Cruikshank (1995). We note that for this species, we computed the uncertainty factor over the mole fraction of N(⁴S) at the corresponding levels, as its abundance is more than two orders of magnitude higher than that of N(²D) at 400 km and six orders of magnitude higher at 200 km. At 400 km, we find a nominal number density of 8.38×10⁷ cm⁻³. Considering chemical uncertainties, we have a mean value of 6.42×10⁷ cm⁻³ with an uncertainty factor of 1.31. Thus, these values are consistent with the measurements of (1±0.25)×10⁸ cm⁻³ at this level. In the 196–204 km range, we have nominal values of 7.65–6.93×10⁸ cm⁻³ and mean values of 6.33–5.70×10⁸ cm⁻³ considering the chemical

uncertainties, with uncertainty factors of 1.37 and 1.36, respectively, in agreement with (5±2.5)×10⁸ cm⁻³ at 200 km measured from UVS.

An important outcome of this work is that we find nominal results consistent with the observations of Voyager 2 for the electron density profiles and the N₂ and N number densities, without having to apply ad hoc modifications to the electron-impact ionization profile (e.g., a shift toward higher altitudes), as was done in previous work (Summers & Strobel 1991; Krasnopolsky & Cruikshank 1995; Strobel & Summers 1995; Benne et al. 2022). In the latter article, the electron-impact ionization was computed using the ionization profile computed by Strobel et al. (1990). This profile was shifted by 100 km to higher altitudes and multiplied by an orbital scaling factor, to take the intermittence of the precipitation into account. This orbital scaling factor was taken from Krasnopolsky et al. (1993), who considered this factor as a free parameter of their model to be able to fit the observed electron density. In this work, we derived the orbital scaling factor from the modeling of the Neptune-Triton system by considering that this factor is the ratio of the time when Triton is near the magnetic equator (and therefore, where the precipitation occurs) over the total time considered. Therefore, this approach is only based on physical considerations.
## 5.4. Model result uncertainties and main chemical pathways

The study of the uncertainties on the model results shows us that despite the changes made in the modeling of the electron-impact ionization and dissociation reactions, the uncertainty factors remain large for the main species at the altitude where their uncertainty is maximum, even though we note some variations in comparison to Benne et al. (2022). This is shown in Table B.2, which gives the uncertainty factors *F* for all the main species of the Triton atmosphere at the level where the uncertainty on their mole fraction is maximum. We recall that the uncertainty factor *F* was used to find the $1\sigma$ interval around the mean mole fraction of a species, $\bar{y}_i$, such that this interval is $|\bar{y}_i/F, \bar{y}_i \times F|$ ($\bar{y}_i$ is not necessarily equal to the nominal value of the mole fraction at the same altitude level). Finding large uncertainty factors is logical because the electron-impact ionization and dissociation of $N_2$ were only a small part of the key uncertainty reactions identified by Benne et al. (2022), and we used the same reaction rates for the other reactions. This confirms the need for new measurements at low temperatures for the key uncertainty reactions, as presented in Benne et al. (2022).

The main chemical pathways mostly remain the same between Benne et al. (2022) and this work. We only observe some modifications of the contribution of the key chemical reactions to the chemical production and loss of the main atmospheric species (key chemical reactions are the crucial reactions producing or destroying a given species, as defined in Benne et al. 2022). Because the main chemical pathways were already detailed in Benne et al. (2022), we only give the key chemical reactions from this work in Table B.3.

## 5.5. Caveats

The current model has some limitations, however. The first limitation is the geometry used for the electron precipitation: We considered it to be vertical, which led to a secondary electron production by these electrons that is maximum below 100 km. The magnetic field lines are expected to drape themselves around Triton (Strobel et al. 1990; Sittler & Hartle 1996), as is the case for Titan (Gronoff et al. 2009b; Modolo & Chanteur 2008). Strobel et al. (1990) therefore considered a magnetic field line inclined by 9° and Sittler & Hartle (1996) worked with a horizontal line. The number of different configurations with a draping field is so large that a systematic study would be meaningless. We chose the configuration that was most accurate and prevented any additional uncertainty due to the projection toward the vertical. This configuration is found on Triton like any other and is similar to that used to study Titan in Dobrijevic et al. (2016). In the future, the study of other magnetic configurations as in Gronoff et al. (2009b) will allow us to address how the energy is deposited in the atmosphere. For now, we did not compute the deposited energy flux, and we therefore cannot give any constraints on the heating of the atmosphere by magnetospheric electrons from our model. However, this was a critical argument supporting the electron precipitation hypothesis because the solar flux seemed too weak to explain the observed thermospheric temperature (Stevens et al. 1992; Krasnopolsky et al. 1993; Strobel & Zhu 2017). Again, as the electron-impact ionization and electron-impact dissociation reactions resulting from electron precipitation mostly take place below 100 km, we expect the energy to be deposited in this area as well, which would certainly not allow us to explain the thermospheric temperature (if a secondary energy source is needed). However, considering a more inclined magnetic field line as described in Sect. 5.2 should

shift the region in which the electrons deposit their energy toward a higher altitude. Therefore, it should help to reproduce the thermospheric temperature.

Moreover, better constraints on the electron densities inside the Neptune magnetosphere and better constraints of the magnetic field strength at Triton depending on its position in the magnetosphere would allow us to improve our model significantly. This could be achieved with an orbiter equipped with a magnetometer and a plasma/energetic particles experiment that would make multiple flybys of Triton, as Cassini did for Titan.

## 5.6. Electron precipitation is still needed to explain the Voyager 2 observations

We tested our photochemical model without considering electron precipitation and found an electron density profile consistent with the Voyager 2 observations presented in Tyler et al. (1989), as shown in Fig. C.1. However, even though a consistent electron density profile can be found without electron precipitation, it is more likely that this precipitation should be implemented in models of the Triton atmosphere (Lyons et al. 1992). This additional source of energy is needed to explain the thermospheric temperature according to Stevens et al. (1992) and Krasnopolsky et al. (1993). In addition, it is needed to obtain a number density of atomic nitrogen at 400 km that matches the data presented in Krasnopolsky et al. (1993).

# 6. Sensitivity study of the electron temperature

As explained in previous sections, the electron temperature $T_e$ was not explicitly provided in the majority of the articles about the photochemistry of the Triton atmosphere because no measurements of this physical quantity were performed during the Voyager 2 flyby. We therefore used an electron temperature equal to the neutral temperature in our work. Because the atmosphere of Triton is tenuous, however, we can expect $T_e$ to be higher than the neutral temperature $T_n$. Using various hypotheses to compute the pressure balance, Sittler & Hartle (1996) found that $T_e$ at the Triton ionopause should be between 308 and 1230 K. These values are higher than the maximum neutral temperature reached in the thermosphere, which is $102\pm3$ K according to Krasnopolsky et al. (1993). In their photochemical model, Krasnopolsky (2023) assumed $T_e = 300$ K, but did not justify this choice.

As the electron temperature remains unknown, we performed multiple simulations with electron temperatures higher than the neutral temperature. Based on the results of Sittler & Hartle (1996), our three main simulations considered $T_e = 300, 10 \times T_n$ (resulting in a maximum electron temperature of 920 K) and 1230 K. The results from these simulations are shown in Fig. 12 and Table 3.

A higher electron temperature logically leads to a higher electron peak density because a higher electron temperature lowers the rates of dissociative and radiative recombination reactions. We observe in Table 3 that the simulation with $T_e = 300$ K is consistent with Voyager 2 data from Tyler et al. (1989) at $1\sigma$, which gives a maximum electron density of $(3.5\pm1) \times 10^4$ cm$^{-3}$. The simulations with $T_e = 10 T_n$ and 1230 K give marginally higher electron peak densities at $1\sigma$ because the values from Tyler et al. (1989) are given without any uncertainties. Moreover, the nominal electron density profiles remain in the uncertainty range of our reference case where $T_e = T_n$.
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Nominal electron density profiles for $T_e = T_n$ (solid red line), $T_e = 300$ K (dash-dotted green line), $T_e = 10 \times T_n$ (dashed blue line) and $T_e = 1230$ K (dotted black line). The gray profiles are the electron density profiles obtained considering chemical uncertainties for the $T_e = T_n$ case. The red triangles and blue circles give the electron density profiles measured during the Voyager 2 flyby from Tyler et al. (1989).</figcaption></figure>

<table><caption>**Table 3.** Comparison of the maximum electron density and its altitude for various cases using a different electron temperature profile.</caption><thead><tr><th>$T_e$</th><th>Max e<sup>-</sup> density (cm<sup>-3</sup>)</th><th>Max e<sup>-</sup> density at 1-σ (cm<sup>-3</sup>)</th><th>e<sup>-</sup> density peak (km)</th></tr></thead><tbody><tr><td>$T_n$</td><td>4.2 × 10<sup>4</sup></td><td>[2.7–6.4] × 10<sup>4</sup></td><td>367</td></tr><tr><td>300 K</td><td>6.0 × 10<sup>4</sup></td><td>[3.9–9.2] × 10<sup>4</sup></td><td>367</td></tr><tr><td>10 $T_n$</td><td>8.0 × 10<sup>4</sup></td><td>[5.2–12.2] × 10<sup>4</sup></td><td>367</td></tr><tr><td>1230 K</td><td>1.0 × 10<sup>5</sup></td><td>[6.5–15.3] × 10<sup>4</sup></td><td>367</td></tr></tbody></table>

**Notes.** The $1\sigma$ interval is computed with the uncertainty factor of 1.53 that was calculated at the electron peak in Sect. 5.2 for the $T_e = T_n$ case.

Therefore, increasing the electron temperature in our model leads to a higher electron density, but this increase is not significant considering the chemical uncertainties and that the profiles from Tyler et al. (1989) are given without uncertainties. In addition, recent reprocessing of Voyager 2 data by Togni et al. (2023) indicates that the electron density measured during the flyby is likely to be higher than what is presented in Tyler et al. (1989). To make our nominal electron density profiles consistent with the profiles from Tyler et al. (1989), we could also have added another unconstrained loss process in our model such as, for example, ion escape at the top of the atmosphere, as was done in Krasnopolsky & Cruikshank (1995); Krasnopolsky (2023).

## 7. Conclusion

We presented a model that was improved over the model of Benne et al. (2022) to better account for the interaction between the Triton atmosphere and the Neptune magnetosphere. To do this, we first modeled the Neptune-Triton system to study the variation in the magnetic environment of Triton depending on its position in the Neptune magnetosphere. The L-shell parameter of Triton varies strongly with time from 14 to more than 100. Thus, Triton is not always near the magnetic equator of Neptune, where

the magnetospheric electrons are located. Therefore, we computed the time ratio when Triton is in this area, called the orbital scaling factor, and applied it to compute the mean precipitation, as done in Strobel et al. (1990). The initial input precipitation flux was taken from the latter article, but was then modified following the work of Sittler & Hartle (1996), in which the electron precipitation rate depends on their energy. We also computed the mean magnetic field using OTD models of the Neptune magnetic field that were developed after the flyby of Voyager 2. This gave a mean field of 5 nT, slightly lower than the value of 8 nT used by Strobel et al. (1990) and Sittler & Hartle (1996). This difference impacts the energy deposited by the electrons. We considered that the electron flux propagated vertically, that is, perpendicular to the atmospheric layers.

With these inputs, we used the electron transport model TRANSPlanets to compute the electron-impact ionization and electron-impact dissociation rates. TRANSPlanets is a modified version of the TRANS model that has been used to study electron transport in various planetary atmospheres Gronoff et al. (2009a). These rates were then used in our photochemical model after coupling it with TRANSPlanets. This significantly changed the results of the photochemical model compared to the case where the electron-impact ionization and electron-impact dissociation rates were computed as in Benne et al. (2022). The electron impact-ionization and dissociation integrated column rates decreased by 66% between the two models, thus affecting the atmospheric composition. Consequently, electron-impact ionization is no longer the dominant ionization source in the Triton atmosphere because its rate is slightly lower than the rate of photoionization. This is opposite to what was found in the previous photochemical models of the Triton atmosphere (Strobel & Summers 1995; Krasnopolsky & Cruikshank 1995; Benne et al. 2022). However, we note that this result depends on the input flux used in TRANSPlanets, which depends on some of our assumptions because we lack knowledge of the electron fluxes at Triton, the time variation of its magnetic environment, and of the magnetosphere-atmosphere interaction. Maybe more importantly, these results also depend on the shape of the magnetic field around Triton. We made the strong assumption of a vertical magnetic field line. Considering the multiplicity of the different configurations of a satellite orbiting a complex magnetized planet, there is no doubt that this configuration may be found, and no doubt either that it is not a general rule. Therefore, this work must be considered as a first approach, showing that we are now able to address a detailed and deep study of the upper atmosphere of Triton. Observations are still missing, however. We strongly need a dedicated mission to study the Neptune magnetosphere in depth, with several flybys of Triton when it is at different places in the magnetosphere.

Even though the atmospheric composition was altered, the main chemical pathways remained the same as in Benne et al. (2022). The changes come from the varying contribution of the different key chemical reactions. With this model, we find a nominal electron peak number density consistent with the Voyager observations presented in Tyler et al. (1989), at an altitude that is roughly consistent with observations. The nominal number densities of $N_2$ and N are also consistent with the observations reported in Broadfoot et al. (1989); Krasnopolsky et al. (1993); Krasnopolsky & Cruikshank (1995). We performed a Monte Carlo simulation with 250 runs of the model by varying the chemical reaction rates. We found that some uncertainty factors were lower in this work compared to Benne et al. (2022), but the overall uncertainties are still large. We also showed that using increased electron temperature leads to higher electron
densities, but these densities remain consistent with the Voyager 2 observations presented in Tyler et al. (1989) when we consider the chemical uncertainties and that the profiles from Tyler et al. (1989) did not display any uncertainty.

The current model can be improved in the following ways. The possibility of using a horizontal magnetic field line could have a significant impact on the results by shifting the electron-impact ionization profile resulting from the electron precipitation upward. This case would also model the Triton-magnetosphere interaction better because it is thought that the magnetic field lines drape themselves around Triton (Strobel et al. 1990; Sittler & Hartle 1996). The computation of the heat deposited by the electron precipitation would also help to constrain the importance of the geometry of the magnetic field lines in the upper atmosphere of Triton because magnetospheric electrons may explain the thermospheric temperature (Stevens et al. 1992; Krasnopolsky et al. 1993; Strobel & Zhu 2017). This model would also benefit from a better knowledge of the Neptune magnetosphere, such as the electron number density depending on the L-shell to better model the input precipitation flux. The acquisition of this type of data requires an orbiter in the Neptunian system.

*Acknowledgements.* B. Benne, M. Dobrijevic and T. Cavalié acknowledge funding from CNES. B. Benne, M. Dobrijevic, T. Cavalié and K.M. Hickson acknowledge funding from the "Programme National Planétologie" (PNP) of CNRS/INSU. K.M. Hickson acknowledges support from the French program "Physique et Chimie du Milieu Interstellaire" (PCMI) of the CNRS/INSU with the INC/INP co-funded by the CEA and CNES.

## References

Archinal, B. A., Acton, C. H., A'Hearn, M. F., et al. 2018, Celest. Mech. Dyn. Astron., 130, 22

Benmahi, B. 2022, PhD thesis, Université de Bordeaux, France

Benne, B., Dobrijevic, M., Cavalié, T., Loison, J.-C., & Hickson, K. M. 2022, A&A, 667, A169

Broadfoot, A. L., Atreya, S. K., Bertaux, J. L., et al. 1989, Science, 246, 1459

Connerney, J. E. P., Acuña, M. H., & Ness, N. F. 1991, J. Geophys. Res.: Space Phys., 96, 19023

Dobrijevic, M., Loison, J., Hickson, K., & Gronoff, G. 2016, Icarus, 268, 313

Fox, J., & Victor, G. 1988, Planet. Space Sci., 36, 329

Gronoff, G. 2009, PhD thesis, Université Joseph Fourier, France

Gronoff, G., Lilensten, J., Simon, C., et al. 2007, A&A, 465, 641

Gronoff, G., Lilensten, J., Simon, C., et al. 2008, A&A, 482, 1015

Gronoff, G., Lilensten, J., Desorgher, L., & Flückiger, E. 2009a, A&A, 506, 955

Gronoff, G., Lilensten, J., & Modolo, R. 2009b, A&A, 506, 965

Gronoff, G., Mertens, C., Lilensten, J., et al. 2011, A&A, 529, A143

Gronoff, G., Hegyi, B., Wedlund, C. S., & Lilensten, J. 2021, https://doi.org/10.5281/zenodo.4632426

Hébrard, E., Dobrijevic, M., Loison, J. C., Bergeat, A., & Hickson, K. M. 2012, A&A, 541, A21

Itikawa, Y. 2015, J. Phys. Chem. Ref. Data, 44, 013105

Jonauskas, V. 2018, A&A, 620, A188

Jonauskas, V. 2022, A&A, 659, A11

Karamcheti, K., & Sentman, L. H. 1968, J. Chem. Phys., 48, 2181

Kato, T. 1994, Atom. Data Nucl. Data Tables, 57, 181

Krasnopolsky, V. A. 2023, Icarus, 406, 115741

Krasnopolsky, V. A., & Cruikshank, D. P. 1995, J. Geophys. Res., 100, 21271

Krasnopolsky, V. A., Sandel, B. R., Herbert, F., & Vervack, R. J. 1993, J. Geophys. Res., 98, 3065

Krimigis, S. M., Armstrong, T. P., Axford, W. I., et al. 1989, Science, 246, 1483

Lilensten, J. 1989, PhD thesis, Institut National Polytechnique de Grenoble, France

Lilensten, J., Simon, C., Witasse, O., et al. 2005a, Icarus, 174, 285

Lilensten, J., Witasse, O., Simon, C., et al. 2005b, Geophys. Res. Lett., 32, L03203

Liuzzo, L., Paty, C., Cochrane, C., et al. 2021, J. Geophys. Res.: Space Phys., 126, e2021JA029740

Lyons, Yung, Y., & Allen, M. 1992, Science, 256, 204

Majeed, T., McConnell, J. C., Strobel, D. P., & Summers, M. E. 1990, Geophys. Res. Lett., 17, 1721

Mauk, B. H., Krimigis, S. M., Cheng, A. F., & Selesnick, R. S. 1995, Neptune and Triton, 169

Menager, H. 2011, PhD thesis, Université de Grenoble, France

Modolo, R., & Chanteur, G. M. 2008, J. Geophys. Res.: Space Phys., 113, A01317

Ness, N. F., Acuña, M. H., Burlaga, L. F., et al. 1989, Science, 246, 1473

Ness, N. F., Acuña, M. H., & Connerney, J. E. P. 1995, Neptune and Triton, 141

Nicholson, W. P., Gronoff, G., Lilensten, J., Aylward, A. D., & Simon, C. 2009, MNRAS, 400, 369

Simon, C. 2006, PhD thesis, Université Joseph Fourier, France

Simon, C., Witasse, O., Leblanc, F., Gronoff, G., & Bertaux, J. L. 2009, Planet. Space Sci., 57, 1008

Sittler, E. C., & Hartle, R. E. 1996, J. Geophys. Res.: Space Phys., 101, 10863

Stamnes, K., & Rees, M. H. 1983a, Geophys. Res. Lett., 10, 309

Stamnes, K., & Rees, M. H. 1983b, J. Geophys. Res.: Space Phys., 88, 6301

Stevens, M. H., Strobel, D. F., Summers, M. E., & Yelle, R. V. 1992, Geophys. Res. Lett., 19, 669

Strobel, D. F., & Summers, M. E. 1995, Neptune and Triton, 1107

Strobel, D. F., & Zhu, X. 2017, Icarus, 291, 55

Strobel, D. F., Cheng, A. F., Summers, M. E., & Strickland, D. J. 1990, Geophys. Res. Lett., 17, 1661

Summers, M. E., & Strobel, D. F. 1991, Geophys. Res. Lett., 18, 2309

Suno, H., & Kato, T. 2006, Atom. Data Nucl. Data Tables, 92, 407

Swartz, W. E., Nisbet, J. S., & Green, A. E. S. 1971, J. Geophys. Res., 76, 8425

Thuillier, G., Floyd, L., Woods, T., et al. 2004, Adv. Space Res., 34, 256

Togni, A., Caruso, A., Buccino, D., et al. 2023, 25th EGU General Assembly Conference Abstracts, online at https://egu23.eu/, EGU-9714

Tyler, G. L., Sweetnam, D. N., Anderson, J. D., et al. 1989, Science, 246, 1466

Wakelam, V., Herbst, E., Loison, J.-C., et al. 2012, ApJS, 199, 21

Wang, Y., Zatsarinnyo, O., & Bartschat, K. 2014, Phys. Rev. A, 89, 062714

Witasse, O., Dutuit, O., Lilensten, J., et al. 2002, Geophys. Res. Lett., 29, 104

Witasse, O., Dutuit, O., Lilensten, J., et al. 2003, Geophys. Res. Lett., 30, 1360
**Appendix A: Results obtained with a not renormalized precipitation flux**

We present results obtained with TRANSPlanets and the photochemical model when the input precipitation flux used in the former model is *j*'(E), that is, the flux that is not renormalized. Thus, this flux carries significantly more energy that can be deposited by curvature drift, as computed with Eqs. (13) and (14), which were taken from Strobel et al. (1990). Fig. A.1 shows the secondary electron production rate, and Fig. A.2 plots the mole fractions of the main neutral species and the main ions of the Triton atmosphere.

Fig. A.2 shows that both the mole fraction of atomic nitrogen and the electron density profiles are inconsistent with Voyager 2 data. They were taken from Krasnopolsky et al. (1993) and Tyler et al. (1989).

<figure><img src="image_14.png" /><figcaption>Fig. A.1. Secondary electron production profile computed with TRANSPlanets using a not renormalized precipitation flux *j*'(E).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. A.2.** Results from the photochemical model using electron-impact ionization and electron-impact dissociation rates computed by TRANS-Planets with a not renormalized precipitation flux *j*'(E). (a) Mole fractions of the main neutral species. The horizontal bars represent Voyager 2 measurements as presented in Krasnopolsky et al. (1993). (b) Number density profiles of the main ions. The electron profile can be compared with the profiles presented in Tyler et al. (1989).</figcaption></figure>
## Appendix B: Outputs from the nominal model using TRANSPlanets

We present in Table B.1 the integrated reaction rates of the electron-impact ionization and electron-impact dissociation reactions in Benne et al. (2022) and in this work. Table B.2 presents the altitude at which the uncertainty on the mole fraction of the main atmospheric species is maximum for this work, along with the uncertainty factor on this mole fraction. Finally, Table B.3 compares the key chemical reactions of Benne et al. (2022) and this work.

<table><caption>**Table B.1.** Comparison of the integrated column rates of the electron-impact ionization and electron-impact dissociation reactions</caption><thead><tr><th>Electron-impact ionization or<br>Electron-impact dissociation reaction</th><th>Integrated column<br>rates (cm<sup>-2</sup>.s<sup>-1</sup>)<br>Benne et al. (2022)</th><th>Integrated column<br>rates (cm<sup>-2</sup>.s<sup>-1</sup>)<br>This work</th><th>Δ<br>relative</th></tr></thead><tbody><tr><td>N<sub>2</sub> + e<sup>-</sup> → N<sub>2</sub><sup>+</sup> + 2e<sup>-</sup></td><td>6.6×10<sup>7</sup></td><td>2.4×10<sup>7</sup></td><td>-63.6%</td></tr><tr><td>N<sub>2</sub> + e<sup>-</sup> → N<sup>+</sup> + N(<sup>2</sup>D) + 2e<sup>-</sup></td><td>1.6×10<sup>7</sup></td><td>3.9×10<sup>6</sup></td><td>-75.6%</td></tr><tr><td>N<sub>2</sub> + e<sup>-</sup> → N(<sup>4</sup>S) + N(<sup>2</sup>D) + e<sup>-</sup></td><td>4.9×10<sup>7</sup></td><td>2.1×10<sup>7</sup></td><td>-57.1%</td></tr><tr><td>N(<sup>4</sup>S) + e<sup>-</sup> → N<sup>+</sup> + 2e<sup>-</sup></td><td>-</td><td>1.0×10<sup>5</sup></td><td>-</td></tr><tr><td>N(<sup>2</sup>D) + e<sup>-</sup> → N<sup>+</sup> + 2e<sup>-</sup></td><td>-</td><td>2.3×10<sup>3</sup></td><td>-</td></tr><tr><td>C + e<sup>-</sup> → C<sup>+</sup> + 2e<sup>-</sup></td><td>-</td><td>9.2×10<sup>3</sup></td><td>-</td></tr><tr><td>CO + e<sup>-</sup> → CO<sup>+</sup> + 2e<sup>-</sup></td><td>-</td><td>2.0×10<sup>4</sup></td><td>-</td></tr><tr><td>CO + e<sup>-</sup> → C + O(<sup>3</sup>P) + e<sup>-</sup></td><td>-</td><td>6.5×10<sup>3</sup></td><td>-</td></tr><tr><td>Sum</td><td>1.31×10<sup>8</sup></td><td>4.9×10<sup>7</sup></td><td>-62.6%</td></tr><tr><td>Sum for electron-impact ionization</td><td>8.2×10<sup>7</sup></td><td>2.8×10<sup>7</sup></td><td>-65.8%</td></tr></tbody></table>

<table><caption>**Table B.2.** Uncertainty factor on the mean abundance $F(\bar{y}_i)$ of the main atmospheric species.</caption><thead><tr><th>Species</th><th>Altitude (km)</th><th>F(\bar{y}_i)</th></tr></thead><tbody><tr><td>H<sub>2</sub></td><td>943</td><td>1.14</td></tr><tr><td>C</td><td>46</td><td>2.74</td></tr><tr><td>H</td><td>31</td><td>9.17</td></tr><tr><td>CH<sub>4</sub></td><td>153</td><td>9.28</td></tr><tr><td>C<sub>2</sub>H<sub>2</sub></td><td>127</td><td>3.14</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub></td><td>41</td><td>10.05</td></tr><tr><td>C<sub>2</sub>H<sub>6</sub></td><td>86</td><td>6.98</td></tr><tr><td>HCN</td><td>438</td><td>5.39</td></tr><tr><td>O(<sup>3</sup>P)</td><td>31</td><td>52.72</td></tr><tr><td>N<sub>2</sub></td><td>1026</td><td>1.01</td></tr><tr><td>CO</td><td>1026</td><td>1.07</td></tr><tr><td>N(<sup>4</sup>S)</td><td>31</td><td>26.44</td></tr><tr><td>N(<sup>2</sup>D)</td><td>115</td><td>42.01</td></tr><tr><td>H<sup>+</sup></td><td>127</td><td>5.36</td></tr><tr><td>C<sup>+</sup></td><td>121</td><td>5.60</td></tr><tr><td>N<sup>+</sup></td><td>923</td><td>1.59</td></tr><tr><td>N<sub>2</sub><sup>+</sup></td><td>1026</td><td>1.90</td></tr><tr><td>e<sup>-</sup></td><td>146</td><td>2.00</td></tr></tbody></table>

**Notes.** These values were computed at the level at which the uncertainty on the abundance of the studied species is maximum. The uncertainty factor $F$ gives the interval $[\bar{y}_i, \bar{y}_i \times F]$ at $1\sigma$.
<table><caption>**Table B.3.** Comparison of the key chemical reactions obtained in this work with a slightly modified version of the model of Benne et al. (2022), based on the nominal results.</caption><thead><tr><th rowspan="2">Reaction</th><th colspan="2">This work</th><th colspan="2">Benne et al. (2022)</th></tr><tr><th>Species (Production)</th><th>Species (Loss)</th><th>Species (Production)</th><th>Species (Loss)</th></tr></thead><tbody><tr><td>CH<sub>4</sub>+hν→<sup>1</sup>CH<sub>2</sub>+H<sub>2</sub></td><td>H<sub>2</sub>(33%)</td><td>CH<sub>4</sub>(34%)</td><td>H<sub>2</sub>(34%)</td><td>CH<sub>4</sub>(35%)</td></tr><tr><td>CH<sub>4</sub>+hν→CH<sub>3</sub>+H</td><td>H(27%)</td><td>CH<sub>4</sub>(29%)</td><td>H(30%)</td><td>CH<sub>4</sub>(30%)</td></tr><tr><td>H+<sup>3</sup>CH<sub>2</sub>→CH+H<sub>2</sub></td><td>H<sub>2</sub>(47%)</td><td>H(53%)</td><td>H<sub>2</sub>(48%)</td><td>H(56%)</td></tr><tr><td>CH+CH<sub>4</sub>→C<sub>2</sub>H<sub>4</sub>+H</td><td>C<sub>2</sub>H<sub>4</sub>(73%); H(21%)</td><td>CH<sub>4</sub>(23%)</td><td>C<sub>2</sub>H<sub>4</sub>(73%); H(23%)</td><td>CH<sub>4</sub>(22%)</td></tr><tr><td>H+HCNN→<sup>1</sup>CH<sub>2</sub>+N<sub>2</sub></td><td>N<sub>2</sub>(23%)</td><td>H(32%)</td><td>N<sub>2</sub>(17%)</td><td>H(34%)</td></tr><tr><td>N<sub>2</sub>+hν→N<sub>2</sub><sup>+</sup>+e<sup>-</sup></td><td>e<sup>-</sup>(48%); N<sub>2</sub><sup>+</sup>(55%)</td><td></td><td>e<sup>-</sup>(26%); N<sub>2</sub><sup>+</sup>(31%)</td><td></td></tr><tr><td>N<sub>2</sub>+hν→N<sup>(4</sup>S)+N<sup>(2</sup>D)</td><td>N<sup>(4</sup>S)(11%); N<sup>(2</sup>D)(20%)</td><td></td><td></td><td></td></tr><tr><td>N<sub>2</sub>+e<sup>-</sup>→N<sub>2</sub><sup>+</sup>+2e<sup>-</sup></td><td>e<sup>-</sup>(39%); N<sub>2</sub><sup>+</sup>(45%)</td><td></td><td>e<sup>-</sup>(57%); N<sub>2</sub><sup>+</sup>(69%)</td><td></td></tr><tr><td>N<sub>2</sub>+e<sup>-</sup>→N<sup>(4</sup>S)+N<sup>(2</sup>D)+e<sup>-</sup></td><td>N<sup>(4</sup>S)(12%); N<sup>(2</sup>D)(21%)</td><td></td><td>N<sup>(4</sup>S)(13%); N<sup>(2</sup>D)(24%)</td><td></td></tr><tr><td>N<sub>2</sub><sup>+</sup>+H<sub>2</sub>→N<sub>2</sub>H<sup>+</sup>+H</td><td></td><td>H<sub>2</sub>(51%); N<sub>2</sub><sup>+</sup>(37%)</td><td></td><td>H<sub>2</sub>(52%); N<sub>2</sub><sup>+</sup>(15%)</td></tr><tr><td>N<sub>2</sub><sup>+</sup>+e<sup>-</sup>→N<sup>(4</sup>S)+N<sup>(2</sup>D)</td><td>N<sup>(2</sup>D)(16%)</td><td>e<sup>-</sup>(26%); N<sub>2</sub><sup>+</sup>(29%)</td><td>N<sup>(4</sup>S)(11%); N<sup>(2</sup>D)(20%)</td><td>e<sup>-</sup>(35%); N<sub>2</sub><sup>+</sup>(42%)</td></tr><tr><td>N<sub>2</sub><sup>+</sup>+e<sup>-</sup>→N<sup>(2</sup>D)+N<sup>(2</sup>D)</td><td>N<sup>(2</sup>D)(27%)</td><td>e<sup>-</sup>(22%); N<sub>2</sub><sup>+</sup>(25%)</td><td>N<sup>(2</sup>D)(33%)</td><td>e<sup>-</sup>(30%); N<sub>2</sub><sup>+</sup>(36%)</td></tr><tr><td>N<sup>(4</sup>S)+CN→N<sub>2</sub>+C</td><td>C(72%)</td><td>N<sup>(4</sup>S)(20%)</td><td>C(69%); N<sub>2</sub>(17%)</td><td>N<sup>(4</sup>S)(37%)</td></tr><tr><td>N<sup>(2</sup>D)+CO→N<sup>(4</sup>S)+CO</td><td>N<sup>(4</sup>S)(56%)</td><td>N<sup>(2</sup>D)(93%)</td><td>N<sup>(4</sup>S)(48%)</td><td>N<sup>(2</sup>D)(77%)</td></tr><tr><td>H<sup>+</sup>+HCN→HNC<sup>+</sup>+H</td><td></td><td>H<sup>+</sup>(22%); HCN(87%)</td><td></td><td>H<sup>+</sup>(28%); HCN(70%)</td></tr></tbody></table>

**Notes.** The key chemical reactions are the reactions that contribute more than 10% of the total integrated production or loss of at least two of the main atmospheric species from this work.

## Appendix C: Results from our photochemical model without electron precipitation

In Fig. C.1 we present the number densities of the main ions of the Triton atmosphere obtained with our photochemical model without electron precipitation from the Neptune magnetosphere. Therefore, the rates of the electron-impact ionization and electron-impact dissociation reactions are zero.

<figure><img src="image_16.png" /><figcaption>**Fig. C.1.** Number density profiles of the main ions of the Triton atmosphere obtained with the photochemical model without the electron precipitation. The electron profile is compared to Voyager 2 data from Tyler et al. (1989).</figcaption></figure>