Dayside ionosphere of Titan: Impact on calculated plasma densities due to variations in the model parameters

A one dimensional photochemical model for the dayside ionosphere of Titan has been developed for calculating the density profiles of ions and electrons under steady state photochemical equilibrium condition. We concentrated on the T40 flyby of Cassini orbiter and used the in-situ measurements from instruments onboard Cassini as input to the model. An energy deposition model is employed for calculating the attenuated photon flux and photoelectron flux at different altitudes in Titan's ionosphere. We used the Analytical Yield Spectrum approach for calculating the photoelectron fluxes. Volume production rates of major primary ions, like, N2+, N+, CH4+, CH3+, etc due to photon and photoelectron impact are calculated and used as input to the model. The modeled profiles are compared with the Cassini Ion Neutral Mass Spectrometer (INMS) and Langmuir Probe (LP) measurements. The calculated electron density is higher than the observation by a factor of 2 to 3 around the peak. We studied the impact of different model parameters, viz. photoelectron flux, ion production rates, electron temperature, dissociative recombination rate coefficients, neutral densities of minor species, and solar flux on the calculated electron density to understand the possible reasons for this discrepancy. Recent studies have shown that there is an overestimation in the modeled photoelectron flux and N2+ ion production rates which may contribute towards this disagreement. But decreasing the photoelectron flux (by a factor of 3) and N2+ ion production rate (by a factor of 2) decreases the electron density only by 10 to 20%. Reduction in the measured electron temperature by a factor of 5 provides a good agreement between the modeled and observed electron density. The change in HCN and NH3 densities affects the calculated densities of the major ions (HCNH+, C2H5+, and CH5+); however the overall impact on electron density is not appreciable ( < 20%). Even though increasing the dissociative recombination rate coefficients of the ions C2H5+ and CH5+ by a factor of 10 reduces the difference between modeled and observed densities of the major ions, the modeled electron density is still higher than the observation by  ∼ 60% at the peak. We suggest that there might be some unidentified chemical reactions that may account for the additional loss of plasma in Titan's ionosphere.