A Revised Sensitivity Model for Cassini INMS: Results at Titan

Cassini Ion Neutral Mass Spectrometer (INMS) measurements from roughly a hundred Titan encounters over the Cassini mission yield neutral and ion densities systematically lower, by factors approximately 2 to 3, than estimates from several other spacecraft systems, including the Attitude and Articulation Control System, and Navigation system. In this paper we present a new INMS instrument sensitivity model, obtained by re-analyzing (1) the capture and transmission of neutral gas through the instrument, and (2) the detector gain reduction during pre-launch testing. By correcting for an under-estimation of gas leakage out of the instrument into space by the original calibration model, and adjusting for the gain change, the new model brings INMS densities into much closer agreement with the other Cassini systems. Accordingly, the INMS ion densities are revised upward by a constant detector sensitivity correction factor of 1.55±21 %, while the neutral sensitivities have a complex instrument pointing direction dependence, due (mostly) to the effect of the INMS vent and antechamber-to-closed source tube. In the special case of on-ram pointing the neutral densities are revised upward by a constant factor of 2.2±23 %. The corrected neutral and ion sensitivities given here are applicable to all previously published INMS results at Titan, Enceladus and elsewhere in the Saturn system. The new model gives reliable densities at high ram angles, in some cases above 90 degrees, thereby expanding the list of Titan flybys from which INMS densities may be extracted. We apply the model to obtain accurate densities from several off-ram Titan flybys which gave unusual neutral density vs. altitude profiles, or unreasonably high densities, with the original calibration.


Introduction
The primary scientific objectives of the Cassini Ion Neutral Mass Spectrometer (INMS) investigation are: (1) to characterize the composition, structure, and chemical behavior of Titan's upper atmosphere and its interaction with Saturn's magnetosphere; and (2) to investigate the neutral and plasma environments of the rings and icy satellites and their magnetospheric interactions (Waite et al. 2004). Careful determination of the gas densities of the environment is necessary for the proper execution of these objectives. The ion density extraction process has been discussed in detail by Mandt et al. (2012) for ram pointed flybys, and many other aspects of the INMS calibration and instrument response to ions and neutrals have now also been addressed in literature [see, e.g., Cui et al. 2009;Gell and Waite 2015;Vuitton et al. 2008;Waite et al. 2004], including (1) compositional extraction from mass spectra, (2) wall adsorption of 'sticky' gases such as H 2 O and NH 3 , (3) thruster firing contamination, (4) radiation background, (5) high-pressure background effects, (6) residual gas 'tails', (7) highsensitivity detector saturation, and (8) on-going work to refine the instrument calibration with the ground-based refurbished engineering unit. However the systematic difference between the INMS neutral and ion densities with other Cassini systems-which has been become apparent during the last decade over more than 100 Titan flybys-has not yet been addressed.
The INMS measurements of Titan's atmospheric neutral mass density were found to be systematically lower than that measured by the Cassini Attitude and Articulation Control System (AACS) and by Navigation (NAV), both of which use spacecraft drag data and aerodynamic models of the spacecraft to derive mass density [as reported by the Titan Atmospheric Working Group [TAMWG], Lee and Hanover 2005]. RPWS Langmuir probe electron densities can serve as a proxy for total ion densities, while ion densities can be retrieved from CAPS-IBS data (Crary et al. 2009). Comparisons of these sources of ion densities also show systematic differences with the ion densities measured by INMS. We have traced the differences between INMS neutral density measurements and other spacecraft systems in part to gas escape from the "closed" source (CS) through leakage pathways not accounted for in the original INMS calibration model. The escape competes with the ram dynamic pressure entering the INMS, yielding lower CS gas densities than indicated by the original calibration, and resulting in a lower INMS neutral gas sensitivity estimate. Therefore the INMS original calibration model that assumed a sealed CS-used until now to relate the CS density to the ambient density-must be replaced by a more accurate approach. This work describes in-flight density comparisons that lead to a new formulation of the ram enhancement factor in the INMS ion sources and, accordingly, a re-estimation of the instrument sensitivity versus spacecraft speed and pointing. We have also re-analyzed the INMS detector gain, and find that the gain reduction which occurred during pre-launch characterization testing (Waite et al. 2004), had not been accounted for in the original calibration model. This gain reduction contributes to the reduced INMS sensitivity estimate in the model presented here. The new calibration model brings Cassini INMS results into agreement with the other Cassini systems on the Titan ion and neutral atmospheric densities, and enables extraction of Titan neutral density profiles from several previous flybys that yielded unusual and unreasonable densities using the original calibration model.
The paper is organized as follows. In the main text we limit the discussion to the key issues, i.e., a basic explanation of the new calibration model, its origin and application, and key equations (Sect. 3), the discussion of the detector sensitivity correction (Sect. 4

Instrument Description
The INMS (Fig. 1) is a mass spectrometer described in detail by Waite et al. (2004). The instrument is a quadrupole mass spectrometer ( Fig. 2) with two source inlets: (1) a closed source (CS) and antechamber for measurement of neutral gas; and (2) an open source (OS) for the measurement of ions and reactive neutrals without wall collisions with the instrument  The ram factor model considers the sections above the quadrupole mass analyzer within the ion source enclosure (ISE not shown) interior (Fig. 3). The sources are housed in an Ion Source Enclosure (ISE). In flight, measurements alternate in a programmable fashion between the two sources. Gas enters the CS inlet through the entrance aperture as a well collimated molecular beam since the relative gas-spacecraft speed is typically much greater than the molecule thermal speeds. Molecule velocities are isotropized and slowed (thermally accommodated) by colliding (roughly 120 times) with the walls of the inlet's antechamber and antechamber-to-CS tube ) before passing through the tube to the CS ionization region (Fig. 3, labeled with density n s ). The resulting density enhancement in the antechamber and CS is a function of the ram angle between the spacecraft velocity vector and the entrance aperture normal. A typical Titan flyby at 6 km/s produces a ram factor of ∼71 relative to the density for a spacecraft at rest if one neglects the leakage paths out of the CS. Scaling the ram factor by √ T a /T s yields the CS to ambient density ratio at the CS and ambient gas temperatures T s and T a (see Eq. (A.17)). The side vent (Figs. 1 and 3) limits the accumulation of residual gas inside the instrument, and constitutes a major leakage pathway out of the instrument (Table 1). Neutral gas is ionized in the CS using an electron beam (commandable to ∼25 or 70 eV) generated by one of two redundant hot filaments. The OS inlet is located below the CS inlet, and consists of an entrance aperture and ion beam collimator. In-flight OS measurement of ambient ion densities is described by Mandt et al. (2012). Positive ions from the CS or OS are extracted into the electrostatic quadrupole lens (Fig. 3), which is used as a deflector to switch ion trajectories between the two sources (Mahhafy and Lai 1990). The quadrupole lens voltages are set to transfer the ions into the analyzer ion focus lenses and in turn into the quadrupole mass analyzer (QMA). The voltage setting on the ion focus lenses can increase or reduce the sensitivity by affecting how well the beam of ionized neutrals is focused into the mass analyzer.
The QMA selects the mass-to-charge ratio (m/q) in the range of 0.5 to 8.5 and 11.5 to 99.5 u per charge, by way of a quadrupole radio frequency mass filter. The detector fore  Fig. 1) with ion sources and gas inlet system shown, together with the azimuthal (φ), ram (θ ) and vent ram (θ ) angles, spacecraft coordinates axes, gas flows (red arrows: ram flux, blue arrows: thermalized flows), and the approximate locations of densities n s , n os , and n ISE . Conduction paths considered in the ram factor model are colored as indicated by the legend, with number-letter notation (Appendix A) in parentheses. Asterisks: effective area estimated by Monte Carlo simulations. Dagger: effective area estimated by fit to Cassini T85 and T88 data. For remaining pathways, we estimated the effective areas analytically. For clarity, leakage through both CS electron guns is shown here as a single conduction pathway: see Fig. 4 for multiple pathways considered in the CS e-gun conductance calculation. Grey colored components are neglected in the model, since these sections are treated as having either (i) negligible (e.g., CS enclosure around n s , ion focus lenses & QMA), or (ii) maximum (quadrupole lenses) conductance optics consists of three lenses: the Einzel, Mask I and Mask II, which together focus the ions into the detectors. The lens voltage settings can modify the INMS sensitivity depending on how well they focus the beam exiting the QMA into the detectors. Two secondary electron multipliers detect the resulting ion beam operating in a pulse counting mode. These detectors differ in effective sensitivity, with the lower sensitivity detector operating off of secondary electrons produced from the adjacent higher sensitivity detector (Waite et al. 2004). INMS neutral species-dependent sensitivity factors s f (for species f ) are calibrated on the ground by seven factors: (1) the electron impact efficiency in the sources as determined by the ionization cross sections of each neutral species; (2) the efficiency for extraction of ions from the CS or OS into the quadrupole lens; (3) the transmission bandpass of the quadrupole switching lens; (4) the degree of focusing in the ion focus lenses into the mass analyzer; (5) the QMA transmission; (6) the detector fore optics transmission; and (7) the detector efficiency.
The OS has three modes of operation: the ion (OSI), neutral beam (OSNB), and neutral thermal (OSNT) modes. In the OSI and OSNB modes the OS is used to analyze the ram flux (or molecular beam) of ambient ions and neutrals entering the OS at spacecraft speed, respectively. We do not address OSI and OSNB extraction of ion and neutral densities and velocities in this paper: for more details see Mandt et al. (2012) and Waite et al. (2004). Here we restrict our analysis of the OS to its significance as a neutral gas conduction pathway to/from the ambient environment, and the effect on the neutral gas densities in both sources. While the instrument sensitivity correction factor discussed in Sect. 4 applies equally to all four INMS modes (CS, OSI, OSNB, OSNT), the new ram factor modeling only affects the sensitivities of the CS and OSNT modes. The OSNT mode measures the thermalized neutral gas in the OS (Waite et al. 2004), i.e., residual neutral gas molecules which have already undergone collisions in the instrument (contrary to OSNB mode which analyzes the ram neutral beam that has not undergone wall collisions). The OSNT mode thereby corresponds to operating the OS in a manner analogous to the CS. While OSNT mode is very rarely used, we discuss the OSNT sensitivity in this paper because data acquired in this mode during the T85, T88 and T96 flybys is relevant for our analysis (Appendix B).

The New INMS Sensitivity Model: In Summary
We have re-analysed the ram factor and detector gain of the Cassini's spacecraft's Ion Neutral Mass Spectrometer, and have identified an effect of (1) gas leakage pathways out of the instrument, and (2) a reduction of the overall instrument calibration sensitivity s f factors, which largely explains previous systematic differences between the Titan atmospheric neutral and ion density measurements of INMS and other Cassini systems. We carried out the ram factor modelling by solving the gas flow equations through the INMS, considering the gas entry and exit pathways, and internal instrument gas conduction pathways as detailed in Table 1, Figs. 3 and 4, and Appendix A. Since the ram flux arriving at the instrument has a non-isotropic molecular velocity distribution, we performed Monte Carlo simulations (Appendix F) versus spacecraft speed and pointing to estimate gas transmission into the instrument through individual conduction pathways in terms of the 'effective areas' of these pathways (Appendix A). With two exceptions (Table 1), we estimated the effective areas of internal pathways (i.e., not exposed to the ram flux) with standard shape-dependent analytical expressions (Dushman and Lafferty 1962) that assume an isotropic 'thermalized' velocity distribution. Since the gas flow time constants are milliseconds for volatile species , we solved the steady-state problem, with a set of flow balance equations between the CS, OS, ISE and the ambient environment (Appendix A). We obtain a new spacecraft speed and pointing dependent expression for the ambient neutral gas densities n af encountered by the spacecraft, in terms of the count rate X f for each species measured in CSN mode (see also Eqs. (A.18) and (A.22a), (A.22b)): where F and F are given in terms of spacecraft speed and pointing by Eqs. (A.13) and (A.16), T 0 is room temperature (293 K), R 1 is a constant 0.6254, and the estimated (speed and pointing-dependent) R 0 , R 2 and R 3 factors can be (i) obtained from lookup tables (to be made available on the PDS, pending publication of this work), or (ii) approximated by the empirical fitting functions given in Appendix C with the parameter values in Table 2. Speed and pointing are expressed in terms of the speed ratio S 0 (Eq. (A.14)) and the ram (θ ) and azimuthal (ϕ) angles shown in Fig. 3. We calibrated the model against the T85 and T88 Titan flybys during which the spacecraft pointing swept over a wide range of angles (ideal for constraining the model) aiming toward, and then away from the vent (Appendix B). During the inbound portion of these two flybys the INMS operates 'in reverse' as discussed in Appendix B, with the ram flux entering the instrument through the vent, and escaping via the CS and OS. In Fig. 5 we show the pointing dependence of the CS ram factor D s with respect to the spacecraft surfaces in view, and in Fig. 6 we give the range of pointing directions not shadowed by the spacecraft (Appendix D) where the model is applicable. The model is valid in the 'collisionless' regime below ∼4 × 10 −9 kg/m 3 (∼8 × 10 16 molecules/m 3 ) for which the molecular mean free path in all directions (including the ram direction) exceeds the spacecraft size, such that the ram enhancement around the spacecraft is insufficient to scatter significant flux toward or away from the INMS. To date this density limit has not been exceeded during any Titan encounters (Figs. 7,8 and 9). In addition to the revised ram factor modeling, the species-dependent sensitivity factors s f used to date, must be lowered by a (species-independent) correction factor β = 1.55 ± 21 % (again, with the revised values to be made available on the PDS) as explained in Sect. 4. The revised s f values apply equally to the neutral (i.e., the CSN and OSNT modes) and the ion densities (i.e., the OSI and OSNB modes).

Fig. 5 Color:
The estimated CSN ram factor D s for S 0 = 20, vs. θ and ϕ. Sensitivity enhancements can be seen for ram pointing into (i) the antechamber-to-CS tube and (ii) the vent. Brown shaded region: ram directions shadowed by the Cassini spacecraft with 0 degrees CAPS angle. The cases of 90 degrees CAPS angle, and the Huygens probe, are shown in light blue and purple shades, respectively. Here D s does not include the effect of the spacecraft, which is only important for ram directions approaching within a few degrees of the shadow zone (Fig. 6). Lines show the ram direction traces from −400 to +400 sec from closest approach, for the flybys shown in Figs. 12, 13 and 19 (the flybys all had S 0 close to 20). Each trace terminates at the dot (+400 sec)

Fig. 6
Range of pointing directions outside the spacecraft shadow zone (blue) where our revised INMS ram factor modeling is applicable. Contours show the cutoff for S 0 = 5, 10, 20 and 50, within which more than 5% of the ram flux is shadowed by the spacecraft. With increasing S 0 the contours get closer to the shadow zone edge as the ram flux thermal angular spread narrows (in the spacecraft reference frame) with increasing spacecraft speed The original INMS calibration model has assumed as the basis for estimating the ram factor, an 'ideal' closed source with no leakage pathways other than the CS inlet (Fig. 3), resulting in an under-estimation of the ambient gas density required to pressurize the closed source. Contrary to the original model, we find here (Eq. (A.7)) that the total effective area of the leakage paths out of the CS (3.32 mm 2 ) actually exceeds that of the CS inlet (∼2 mm 2 ) as shown in Table 1. After escaping the CS, the vent serves as the primary (though not the only) means of escape of the gas into space.
to the original INMS calibration model (Fig. 10). In the special case of zero degrees INMS ram angle pointing (assuming suprathermal spacecraft speeds, as is typical), the CSN correction factor is approximately a constant: ∼2.2 ± 23 % (Appendix E) and approximately independent of the gas species. These corrections supersede the factor three used by several previous INMS based studies (in anticipation of the more detailed determination reported here).

Correction to the Calibration Sensitivity Factors s f
The calibration sensitivity factors for CH 4 and N 2 were measured pre-launch by exposing the INMS to gas leaked into a laboratory calibration chamber, and comparing (1) the instrument count rate for each species, to (2) the ambient gas density measurement from a separate chamber-mounted pressure sensor. Following a subsequent pre-launch characterization testing phase [detailed by Waite et al. 2004] the decision was made to operate the INMS primary (high-sensitivity) secondary electron multiplier (SEM) detector at 2701 V, rather than the 2500 V initially used for calibration. The voltage increase resulted in an increase of SEM gain, which was estimated at the time to increase the sensitivity by a factor of 1.27. Accordingly, the s f factors for multiple different gas species, measured post-launch in the laboratory, with the INMS refurbished engineering unit, or REU, have been scaled to the original CH 4 and N 2 flight unit calibration, times 1.27. However the characterization tests also subjected the SEM to sufficient counts to degrade its performance, by lowering the gain as shown in Waite et al. (2004). The detector gain has been estimated during pre-launch testing, and in flight, by fitting (by linear regression) the proportionality constant of the SEM anode current to the count rate registered by the pulse counter circuitry. As part of this work, we have re-investigated the measurements performed during the characterization period, and discovered that the gain reduction was, erroneously,

Fig. 10
Correction factor to densities estimated from the original calibration model, vs. ram and azimuthal angle. The correction has an uncertainty of ±23 %. In this example S 0 = 20 typical of a Titan flyby, i.e. m = 28 u per charge for N 2 , u = 6 km/s, T a = 150 K (Eq. (A.14)). Near θ = 15 deg, ϕ = 180 deg the ram flux is aimed into the antechamberto-CS tube, resulting in a high sensitivity spot. At higher θ the factor is highly ϕ dependent due to the effect of the vent. Blue shaded region at right is the spacecraft shadowing zone with CAPS at 90 deg. Note at 0°ram that the factor is ∼2.2. (a) Using Eq. (E.5). The 'noise' in the contour lines originates from statistical error in effective areas estimated by Monte Carlo simulations. (b) Correction factor resulting from the Eqs. (C.1a)-(C.4j) empirical fitting functions not included in the 1.27 sensitivity adjustment factor. Additionally, the factor 1.27 itself appears to have over-estimated the effect of the SEM voltage change because only two test runs were compared (at the end of the characterization period), leading to a large margin of error. As shown in Fig. 11a [based on the data shown in Waite et al. 2004], the pre-launch characterization tests produced a gain decrease of ∼82 % from a starting value of ∼7 × 10 7 over the 77 test runs performed, even after accounting for the SEM voltage change from 2500 to 2701 V. Fortunately most of the voltage pulses (produced by ion strikes on the SEM) are still sufficiently in excess of the discriminator threshold to be registered by the pulse counter, despite the gain loss, and therefore the effect of the gain loss on the count rate is not very drastic. In Fig. 11b we plot the count rate versus gain from the characterization tests: a linear fit (on logarithmic scale) yields: over the relevant gain range, with the non-zero exponent a consequence of the finite pulse height distribution width. By Eq. (4), an 82 % gain loss translates to a factor 1.21 count rate reduction during the characterization tests. Since the species calibration sensitivity factors s f used to date have assumed a factor 1.27 increase, rather than a reduction, we find here that these sensitivity factors should be reduced by 1.21 ± 21 % times 1.27, i.e., β = 1.55 ± 21 %. Gain measurements on the flight unit during the T23, 25, 29, 40, 41, 49, 51, 55, 77, 84, and 98 Titan flybys reveal no statistically discernible trend, implying that β has remained

Explaining 'anomalous' Titan Flybys
The ram factor modeling also provides an explanation for the T11, 20, 46, 47, 62, 85, 88 and 96 Titan encounters, where the original calibration [Appendix E] (1) over-estimates Titan's N 2 density (sometimes drastically), and/or (2) yields an unreasonably large difference between the inbound and outbound measured densities. As shown in Fig. 5, all of the relevant encounters had spacecraft pointing far off ram, and in all cases the vent ram flux (not included in the original calibration model) played a major role in the signal. As shown in Fig. 12 for T11, 20, 46, and 62 (and Fig. 13 for the T47, 85, 88, and 96 flybys discussed in Appendix B), the new ram factor modelling largely corrects the problems, with (aside for spacecraft shadowing effects) the only uncorrected discrepancy being observed at T46 where the inbound/outbound difference appears to be real. Hence the new ram factor modelling not only increases the accuracy of the INMS neutral densities estimates, but also expands the list of Titan encounters from which reliable densities may be extracted.

Fig. 12
Estimated N 2 density vs. altitude from the mass 28 channel for the Cassini T11, 20, 46, and 62 Titan encounters between −400 and 400 sec from CA, using the original 'ideal' calibration (red circles), and the new ram factor model (blue squares) with s f = 4.05 × 10 −10 m 3 /s for N 2 (the factor β is applied to s f ). The arrows indicate the inbound branch of the data. The original calibration yields an inbound/outbound density difference as the instrument pointing changed over the flybys (Fig. 5). However the new model anticipates and largely cancels the effect. High-altitude outbound data (not shown) over-estimate the densities due a significant contribution from residual gas in the instrument. On the outbound portion of T20 the ram pointing is briefly shadowed by the tip of CAPS ELS (at 0 deg rotation angle) as can be seen in Fig. 5, resulting in a brief dip in signal which is uncorrected by the new model. Density measurements in the dip are therefore not valid and are shown only to illustrate the effect. The T11 and 62 flyby pointing (nearly steady through the flyby) was greater than 90 deg ram but close to vent ram. Hence the original calibration drastically under-estimates the instrument sensitivity, and thereby over-estimates the T11 and 62 densities. The new model gives much lower, more reasonable densities where M denotes the spacecraft mass, C the drag coefficient, V = VV the spacecraft velocity in the gas reference frame, and A the spacecraft cross section area projected along V. C can be estimated using formulae given in Stalder and Zurick (1951). We use C = 2.1 ± 0.1, which is the value calculated for a cylinder in free molecular flow, and also that found in Monte Carlo simulations of the gas flow interaction with the Cassini spacecraft (Liechty 2006). This is a reasonable drag coefficient value when compared with results determined using orbital data of Earth-orbiting satellites . AACS measures the acceleration in two ways: (1) directly through an accelerometer (powered on only for selected Titan flybys), or (2) through the atmospheric drag torque T Atm on the spacecraft, determined by way of the torques and duty cycles of the thrusters (most flybys) which fire automatically to stabilize the spacecraft orientation. T Atm satisfies the spacecraft equation of rotational motion (Lee 1999;Lee and Hanover 2005): or, where ↔ I SC denotes the spacecraft's inertia tensor, ω the spacecraft's angular velocity,˙ ω the spacecraft's angular acceleration, T RWA and H RWA the torque and angular momentum of the reaction wheels (zero if the wheels are powered off), and T Thrstr the torque exerted by the thrusters. The spacecraft's inertia tensor is estimated by ground software, and has previously been confirmed by an in-flight calibration technique (Lee and Wertz 2002). An onboard attitude estimator (the Kalman-Bucy filter) provides ω at 125 ms time intervals. T Thrstr is extrapolated by flight software from the estimated force impulse due to all thruster firings, including effects due to the thruster rise and tail-off dynamics. The impulses are converted into the per-axis torque impulses using the known moment arms of all the thrusters. ε contains the torques from solar radiation, Titan's gravity gradient, etc., but these torques are negligible (<1.1 mNm) to first order (Lee and Hanover 2005). The drag acceleration is related to the atmospheric torque by (Feldman et al. 2007;Lee and Hanover 2005;Sarani 2007): or, where r = r CP − r CM and r ⊥ = r − ( r ·V)V, with r CP and r CM the positions of the spacecraft centers of pressure and mass as estimated by ground software (Feldman et al. 2007; Lee and Hanover 2005;Lee and Lim 2013;Sarani 2007Sarani , 2009). We therefore obtain (by substituting Eq. (6b) into (7b), then (7b) into (5b)) the density in terms of thruster torque: where we have dropped the negligible ε term. Similar methodologies have provided atmospheric density estimates of Venus (Espiritus and Tolson 1995) and Mars (Spencer et al. 1999), and of Titan using the Huygens HASI accelerometer data (Fulchignoni 2007;Fulchignoni et al. 2005;Zarnecki et al. 2003). NAV uses radiometric Doppler tracking data to estimate the spacecraft velocity, which requires HGA pointing and transmission to earth at the time of the measurement. All the forces acting on the spacecraft are modeled (such as thruster impulses, gravitational harmonics and atmospheric drag, see Fig. 14) and integrated with the equations of motion along the computed spacecraft trajectory. The calculation is iterated to fit the model parameters to the tracking data by a least squares method. Using this approach a model of the atmospheric density (an exponential profile vs altitude) is fit to yield the drag acceleration (via Eq. (5a)) implied by the tracking data. The total one-sigma (random) error of the result (typically a few percent) includes uncertainty in C (±1.6 %), V (±0.005 %), A (±0.7 %), M (±0.73 %), and the Doppler-based measurement and estimation error (varies between flybys, of order ±1 %). Opportunities for earthward HGA pointing are scarce due to the pointing requirements of other spacecraft systems. Therefore, apart from the T70, T87, and T107 encounters, every low altitude Titan flyby of the Cassini mission is "in the blind", as far as NAV is concerned, meaning that no tracking data is captured during the time Cassini dips into the atmosphere and experiences drag. All NAV can do in these cases is to measure (after the flyby) the net acceleration of drag.
Only the T70, T87 and T107 flybys had continuous HGA pointing to earth, yielding altitude resolution in the Doppler data. To retrieve the neutral mass density during these flybys, the model atmosphere is divided into isothermal altitude layers between 950 and 1300 km. The density and (using Eq. (5a)) the drag acceleration at altitude h are modeled as where ρ i is the density at a reference altitude h i and H i is the scale height for the ith layer of the atmosphere: The layer thickness is defined such that the accumulated V from drag is ∼10 times the noise level of the Doppler during Cassini's traversal of the layer. The layer placement chosen for T87 is shown in Fig. 14. Note that different layer placements were allowed between the inbound and outbound legs to capture variations in the atmosphere. In Figs. 7-9 we compare the INMS Titan closest approach densities for multiple flybys to those of AACS and NAV, with the original and new calibration models. As shown the original calibration model yields systematic deviations of INMS from AACS and NAV, with AACS/INMS and NAV/INMS ratios of 3.30 and 2.56. However the new model yields much better agreement, with ratios of 1.47 and 1.16, respectively. Therefore INMS shows better agreement with NAV than AACS, reflecting the small factor of 1.27 disagreement between NAV and AACS (Fig. 9). The one-sigma errors on the ratios are ± σ 2 INMS + σ 2 AACS , ± σ 2 INMS + σ 2 NAV and ± σ 2 AACS + σ 2 NAV , where σ is the systematic (i.e. model-dependent) standard error of the densities. The one-sigma AACS random error is estimated to be ±5.8 % (Feldman et al. 2007), and those for NAV are usually a few percent shown as in Figs. 8-9. However, we ignore the random measurement error since it contributes negligibly due to the large number of Titan flybys used to compute the ratios. In contrast to the random error, the AACS and NAV systematic errors are presently unknown, and we therefore include only the INMS contributions in Figs. 7-9 (cross-hatched regions). As shown in Fig. 8 σ INMS = 23 % yields agreement (i.e. within error) between INMS and NAV irrespective of any (systematic) NAV error. Interestingly the INMS error by itself is not sufficient to yield agreement with AACS (Fig. 7). This remaining INMS/AACS disagreement, at least 24 % outside INMS error, is currently unexplained. At present we can say only to our knowledge that the disagreement cannot be entirely attributed to uncertainty in INMS modeling and measurements.
The total electron density is measured by RPWS and is a reliable proxy for the total ion density in Titan's atmosphere between 1100-1250 km; a region where high mass negative ions do not constitute a significant fraction of the negative charge density (Well-brock et al. 2013). Above 1250 km the RPWS density at times is higher than the INMS ion density because ion drifts at higher altitudes shift the ions velocities outside of the INMS field of view. Below ∼1100 km a significant contribution of heavy negative ions to the negative charge, and the presence of heavy positive ions above the 100 amu INMS mass range, can also effect the ratio of RPWS to INMS (Crary et al. 2009;Mandt et al. 2012;. Therefore, the ratio of RPWS electron densities to INMS total ion density was evaluated versus altitude between 1100 and 1250 km for 10 flybys at 0°ram angle, and found to be approximately constant in this altitude range. The profile comparison provided an average ratio of 1.37 ± 0.6 using the original INMS calibration model, but the ratio reduces to 0.88 ± 0.39 with the new model. Therefore both models yield agreement within error between INMS and RPWS, but the agreement is better with the new model. Charge balance considerations require close agreement of the ion and electron densities in the 1100-1250 km altitude range, and hence INMS and RPWS are, in effect, making independent measurements of the same quantity (i.e., ion/electron density). Accordingly, when both INMS ion and RPWS electron measurements are available and agree within error, we recommend that the ion/electron density should be evaluated as the mean value of the overlap of the INMS and RPWS error bars (ideally including random and systematic errors of both systems). CAPS-IBS also measures the ion densities. In Titan's atmosphere, the majority of the ion energy of the thermal ionosphere is kinetic energy derived as a result of the velocity of the spacecraft through the near stationary atmosphere, which allows the high energy resolution of the IBS energy spectrometer to be used as a low resolution mass spectrometer (Crary et al. 2009). The energy range of IBS overlaps with the mass range of the major components of Titan's ionosphere allowing total ion densities to be derived from IBS and compared to INMS. We compared the IBS ion densities below 100 amu (the INMS mass range) to INMS total ion densities for 7 flybys. This comparison gives a ratio of 2.48 ± 0.68 with the original INMS calibration, which reduces to 1.60 ± 0.44 with the new calibration. We note that all errors given here are under-estimates since we have not included RPWS and IBS systematic errors.
Later in the Cassini mission the project allocated the T87 (13 Nov 2012) and T107 (10 Dec 2014) flybys during the Solstice Tour to compare the AACS and NAV densities to INMS. These two flybys were selected because the geometry allowed the first density measurements with altitude resolution to be simultaneously acquired from all three systems. The opportunity was also taken to power on the AACS accelerometer (ACC) during T87, yielding four independent density measurements versus altitude (Thrusters, ACC, NAV, and INMS). The NAV measurement required continuous HGA Earth pointing during the flyby, ∼52 degrees from the ram direction. In Fig. 15 we compare the T87 densities vs altitude measured by the AACS thrusters and accelerometer, and NAV, to the INMS mass densities with the new calibration (as estimated from the 28 u per charge N 2 signal). With the exception of a partial inbound INMS-AACS disagreement, all the systems, including INMS, agree within errors on the density at all altitudes.

Impact of New Calibration Model on Previous Titan Studies
Numerous publications using INMS thermal neutral measurements include absolute densities and mixing ratios of neutral species in Titan's atmosphere and the plume of Enceladus. INMS neutral densities have primarily been used for dynamical models of Titan's atmosphere (Bell et al. 2010a(Bell et al. , 2010b; Yelle et al.  Vuitton et al. 2007Vuitton et al. , 2006. Temperature derivations based on the scale height use the ratio of densities at two points in the atmosphere . For θ = 0 fixed flybys the temperature derivations are not affected by the (on ram) 2.2 instrument sensitivity correction constant, which cancels from the ratio. However, other publications have used the INMS density vs. altitude profiles to constrain the dynamics of the upper atmosphere (Bell et al. 2010a(Bell et al. , 2010b. Ideally, these dynamical models [e.g., Bell et al. 2010aBell et al. , 2010b] use densities and mixing ratios at their lower boundaries that are constrained by lower atmospheric observations from the Cassini Composite Infrared Spectrometer (CIRS) and the Huygens Atmospheric Science Investigation (HASI). For a factor 2.2 increase in INMS upper atmospheric densities these models require a warmer thermosphere to match the scaled INMS densities-assuming that these models maintain constant lower boundary conditions. These higher thermospheric densities (and simultaneously warmer temperatures) will alter key aspects of these dynamical simulations, such as the absolute value of the eddy diffusion and binary molecular diffusion coefficients at a given altitude. Some studies have found that matching the scaled up INMS densities alters the altitude of peak methane photodestruction, the implied thermal structure of the background atmosphere, and the methane vertical dynamics. These aspects in turn have been found to alter the inferred escape rates from non-hydrostatic models (Bell et al. 2010a(Bell et al. , 2010b. These models are consistent with the kinetic treatments by Tucker and Johnson (2009) and Tucker et al. (2013) using the same scaled up INMS densities. Alternatively, some studies have matched the scaled up INMS densities but have not found differences in their inferred escape rates of methane [cf., Strobel 2010Strobel , 2012, meaning that the topic of methane escape is still a subject of active investigation and debate.
Increasing the INMS densities by a factor of 2.2 will alter all of the major photochemical dissociation and heating rates. The optical depth calculations increase by the same factor, moving the peak of solar heating and dissociation up by almost a scale height (∼70 km). Chemical reaction rates will be modified as well, but the impact is more difficult to constrain because both production and loss rates are affected.  found that the peak ion densities calculated in their photochemical model shifted to a higher altitude (∼60 km) when the INMS neutral densities were increased by a factor of 2.5, bringing their model results into better agreement with the INMS ion data. These findings have been confirmed by recent modeling (Galand et al. 2010;Richard et al. 2014;Vigren et al. 2013).
Acknowledgements The INMS team acknowledges support from NASA and the Jet Propulsion Laboratory under SwRI subcontract 1405853. Support for this work was provided by a Centre National d'Etudes Spatiales Cassini Participating Scientist grant (to VV). Titan atmospheric densities estimated by the Attitude Control team represent work done by many individuals. T. Barber, T. Burk, A. Feldman, R.S. Lim, Siamak Sarani (formerly with JPL), E.K. Wang, J.L. Webster, and Genevie Yang, our colleagues at JPL, have provided us with invaluable support. We are especially grateful to co-author and Cassini INMS Co-Investigator Hasso Niemman, who passed away July 11, 2013. Hasso leaves a decades long legacy in the field of spaceflight mass spectrometry, and his insight, analysis and detailed knowledge of the INMS were essential to this work. He will be missed by the team and the community.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Appendix A: Neutral Density Ram Factor Modeling
The gas conduction pathways through the INMS are indicated in Table 1 and Figs. 3 and 4 by number-letter pairs. Ambient gas flows to the closed source (CS) via the CS inlet (subscript 0) − consisting of the antechamber aperture + antechamber + tube in series. As shown in Figs. 3 and 4, gas flows between the CS and the ISE through openings in the electron gun structure (3a, 3b, 3c, and 3d, see Fig. 4b), and separately (i.e. in parallel), through the focus lenses (2a) in series with the ion nozzle (2b). The ISE also communicates with the ambient environment through the OS pathways and the vent (5). The OS exchanges gas with the ambient environment through the OS collimator (6), and with the ISE by leakage (i) through the OS assembly (7), i.e., where the OS ionization region resides, or (ii) through the OS collimator (8) via slits. Ambient gas can also bypass the OS assembly by flowing directly through the collimator slits to/from the ISE (4). Gas can flow between the ISE and QMA through the ion focus lenses (Fig. 3), but this pathway is not relevant for a steady state model since there are no other gas conduction pathways out of the QMA enclosure. The conduction pathways considered in the model are summarized in Table 1 and Fig. 3, and the labeled flow diagrams in Fig. 4 show schematically how the closed source, open source and ISE densities couple to the ambient density through the various pathways. The gas density and temperature inside the sources and ISE are approximately spatially homogenous, since the dimensions of these structures are much larger than those of the gas escape pathways. Additionally, for volatile species we estimate the time-constants in these regions to be of the order of milliseconds , and much less than the INMS integration period (IP) of 31 ms. Therefore the gas flows and resulting densities are always in approximate steady state during the measurements.
Considering first the CS, we write the steady-state expression equating ingoing and outgoing gas flows as follows: where n s is the molecule number density in the source and n i is the density outside the source at pathway i. Here we model every conduction pathway as consisting of (i) two apertures with areas a i and a i opposite and adjacent to the closed source and (ii) internal structure with transmission probabilities t i and t i for particles to pass through the pathway. v i and v s are the average magnitude of the molecule velocity components projected into the apertures of i from outside and inside the source. We consider v s to be independent of i due to the approximate spatial homogeneity of source density and temperature. We estimate transmittances of internal conduction pathways, i.e. the lens system and electron gun structures, with the standard analytical expressions for slits and cylinders (Dushman and Lafferty 1962), and that of the CS inlet exposed directly to the ram flux, with Monte Carlo calculations (Appendix E). Starting from the flow balance expression Eq. (A.1) with A = at we solve for the source density: where A is the 'effective' area; related to gas conductance C by a factor of the average aperture projected speed, i.e., C = vA. We note that v is (i) 1/4 of the average molecule speed for internal components exposed to an isotropic velocity distribution, and (ii) a function of spacecraft speed and pointing and ambient gas temperature T a for components exposed to the ram flux. The denominator i A i is the total effective area to exit the close source, which contains the contribution of the CS inlet (A 0 ), the lens system plus ion nozzle (A 2a & A 2b ), and both electron gun structures (2 where & denotes effective areas in series, i.e., We include subscripts '0' and '1', e.g. A 0 → A 00 or A 01 where necessary to indicate (0) ram flux entering the conduction pathway, or (1) thermal gas escaping through the pathway to space, respectively. The Eq. (A.2) numerator is the sum of all ingoing fluxes, consisting of one contribution from the ambient environment n a v 00 A 00 (i.e. via the CS inlet) and another from the ISE where n ise denotes the ISE density. We may therefore write Eq. (A.2) in the form and the total effective area A CS = 3.32 mm 2 between the closed source and ISE is with the effective areas on the RHS given in Table 1. R 1 is a constant number (R 1 = 0.6254) due to its dependence solely on internal conductances, and is unaffected by spacecraft speed and pointing. Notably, A CS exceeds A 01 = 1.989 mm 2 (Table 1), and therefore neutral gas is more likely to "leak" from the CS into the rest of the instrument, than the escape back out through the CS inlet. The ISE steady-state density also satisfies a flow balance expression of the same form as Eq. (A.2): and like the expression for n s given in Eq. (A.5), similar considerations yield an analogous expression for n ise : Substituting Eq. (A.9) into (A.5) we obtain Considering the ambient gas to have a Maxwell-Boltzmann velocity distribution with temperature T a , the speed ratios v 00 /v s and v 00 Here T s and T ise are the CS and ISE gas temperatures, and where s x , s y and s z are the molecule velocity components in units of u 0 = 2kT a /m f (the most probable speed of the ambient gas molecules in the ambient frame) and the integral F is evaluated over the velocity half space pointed toward the conduction pathway. Here θ is the ram angle ( Fig. 3) with respect to the rest frame of the ambient gas, and S 0 is the spacecraft speed u in units of u 0 , i.e., where m f is the molecular mass of species f . Note that we do not assume an equivalence of F with the CS ram enhancement factor, as was done in the original INMS calibration model (Appendix E). We have distinguished the average projected velocity v 00 into the antechamber and open source from the projected velocity into the vent v 50 , which points at 90 degrees to the closed and open sources (Figs. 1 & 3). An expression of the same form as Eq. (A.13) can be written for the speed ratios v 50 /v s and v 50 /v ise , by replacing θ with the angle θ = cos −1 (cos φ sin θ) of the vent normal to the spacecraft direction of motion in the ambient frame (the vent ram angle): where we define the azimuthal angle φ as the left-handed angle of the spacecraft velocity about +X from −Y in spacecraft coordinates (Fig. 3). where is the closed source ram enhancement factor (contrary to the original calibration where D s = F , see Appendix E). Likewise, substitution of Eqs. (A.12b) and (A.15b) into (A.9) gives The detector count rate X [the number of counts divided by the 31 ms IP minus background, after correcting for instrumental effects, ] in CSN (closed source neutral) mode for species f is where s f is the CSN sensitivity for f . The ratio √ T s /T 0 appears in Eq. (A.21) because laboratory s f calibrations are performed with respect to the equilibrium density at room temperature T 0 (∼293 K) surrounding the instrument, which is different from n s by a factor √ T s /T 0 . The combination of Eqs. (A.17) and (A.21) yields the final expressions for count rate and ambient density: Similarly in open source neutral thermal (OSNT) mode, which measures the OS neutral density (Waite et al. 2004), we have with X os , s os , T os , and n os the OSNT count rate, sensitivity, temperature and density, respectively. Writing an expression of the same form as Eq. (A.2) for the OS density, we have since the OS communicates with the ambient density through the OS collimator assembly (path 6 in Table 1, Fig. 4a), and with the ISE through (i) the OS assembly (path 7), and (ii) the OS collimator assembly (path 8). Substituting Eq. (A.9) for the ISE density into Eq. (A.24) above, we obtain From Eqs. (A.23) and (A.27), we obtain the OSNT count rate in terms of the ambient density: Though the INMS rarely operates in OSNT mode, data acquired in this mode during the T85 and T88 Titan encounters are useful to the present analysis as we discuss below.

Appendix B: In Flight Model Calibration: Path 7
The OS assembly conductance (path 7) is cumbersome to model owing to the complexity of this section and limited detail of available drawings. Rather than estimate this pathway analytically or by Monte Carlo simulations, we used in flight data from the T85 and T88 flybys on 24 July and 29 Nov, 2012 during which the vent ram angle θ was swept from below to above 90 deg (Fig. 16). This spacecraft pointing provided a unique opportunity to calibrate the model against its two limiting cases: i.e., entry of gas flux predominantly into (i) the CS inlet and OS collimator, or (ii) the vent. During both encounters INMS measured the CS and OS densities, using the CSN and OSNT modes, in only three channels to allow high (0.1 sec) time resolution: 2, 16 and 28 u per charge; corresponding mostly to Titan Fig. 17 (a) OSNT counts at 2, 16 and 28 u per charge during the T85 flyby, corresponding mostly to Titanian H 2 , CH 4 and N 2 , respectively. The ambient density for each species can be obtained from the count rate X osf (counts divided by the 0.031 sec integration period) by Eq. (A.30) (and, for CSN data, by Eq. (A.22a), (A.22b)). Exhaust from thruster firings is the source of the positive H 2 spikes. Negative N 2 and CH 4 spikes are due to scattering of the ambient atmosphere by the exhaust plume, as shown in Fig. 16 (by comparison N 2 from the thrusters contributes negligibly to signal spikes). (b) Same as (a) for CSN signals atmospheric H 2 , CH 4 and N 2 (Fig. 17). The Titan densities are also uncertain by a factor of (roughly) three due to atmospheric scale height differences between flybys, depending on time and positional variability of solar and/or plasma heating (Cui et al. 2009;Westlake et al. 2011). We therefore cancel the effect of Titan's densities by plotting the ratio of CSN and OSNT, obtained by dividing Eq. (A.30) by (A.22a): where s f = 4.05 × 10 −10 m 3 /s for 14 N 2 . Additionally X f and X osf are affected equally by spacecraft shadowing effects which are significant at T85 and T88 (Appendix D), and therefore this source of uncertainty is also canceled by the ratio. We have fit Eq. (B.1) to the ratio of OSNT to CSN counts at 28 u per charge for T85 and T88, with A 7 and the OSNT sensitivity factor s osf as the fitting parameters. We fit s osf since-by comparison to the s f factors-the OSNT sensitivity is much less well characterized, with only one measurement available for N 2 (Waite et al. 2004). Our fit yields s osf = 1.33(±0.12) × 10 −10 m 3 /s for Fig. 18 Ratio of the OSNT to CSN N 2 signal in the mass 28 channel, during the T85 (circles) and T88 (squares) flybys, together with the optimal model (Eq. (B.1)) fit with OSNT sensitivity s osf = 1.33 × 10 −10 m 3 /s and OS-ISE effective area A 7 = 1.55 × 10 −7 m 2 (solid lines). Dashed line: model prediction with A 7 three times optimal. Dotted line: same with A 7 = 0. The step like decrease in the ratio occurs, as anticipated by the model, when the vent ram angle exceeds 90 degrees. One can see that the success in fitting the data is not very sensitive to uncertainty in A 7 14 N 2 , which agrees within error with the (Waite et al. 2004) value of 1.23 × 10 −10 m 3 /s (their Table IV). We note that the s f and s osf values given here have been reduced by the correction factor by β = 1.55, as discussed in Sect. 4. In Fig. 18 we compare the ratio of OSNT to CSN counts at 28 u per charge at T85 to Eq. (B.1) for different values of A 7 , including the best fit of 1.5(±0.8) × 10 −7 m 2 . As expected both flybys exhibit a step like decrease in the ratio when the vent ram angle approaches 90 deg. Above 90 deg the dominant gas source is through the CS inlet and OS collimator, rather than the vent, resulting in a different distribution of densities between the CS, OS and ISE. The fit is shifted up/down on Fig. 18 by s osf since this factor is just a multiplicative constant in Eq. (B.1), while the magnitude of the step is controlled by A 7 . As shown in Fig. 18 the success of the fit is in fact rather insensitive to uncertainty in A 7 , since most gas flow between the ISE and ambient environment is through pathways 4 and 5 (bypassing path 7) due to their large effective areas (Table 1). Two other flybys: T47 and T96, had similar spacecraft pointing (Fig. 5) to T85 and T88. T96 was at higher altitude (1400 km closest approach) than T47, 85 and 88 (1023, 1012 and 1014 km CA, respectively), and therefore the density was lower, and the data noisier (Fig. 19). Nevertheless, we found that the T96 OSNT/CSN ratio behaves as predicted by the model. A fit of the ratio to T47 data was not possible since OSNT data was not acquired during this flyby.

Appendix C: Ram Factor Empirical Function Fits
Ambient neutral density estimates by way of Eq. (A.22b) require knowledge of D versus θ , φ and S 0 . Since the ram factor modeling depends in part on the results of Monte Carlo simulations, D s must be estimated from lookup tables. The INMS team plans in the near future to provide the lookup tables on the Planetary Data System. However we have also fit the R 0 , R 2 , and R 3 values from the tables with the empirical fitting functions given below, from which D s can be directly calculated (Eq. (A.18)). The empirical functions can serve as a convenient alternative to the tables. Except for R 1 which is constant, the functions are parameterized with respect to the components S 0 , θ , φ of the spacecraft velocity vector in the Here the constant value of R 1 is estimated from the ram factor modeling to be 0.6254. The empirical functions for R 0 , R 2 , and R 3 , which feed into the above expressions, are as follows: R 0 = Q 1 for 0 < φ < π/2 or 3π/2 < φ < 2π (C.1b) where the values of the fitting constants, written in uppercase, are listed in Table 2, and the lowercase quantities are given by Note that we did not fit R 3 for π/2 ≤ φ ≤ 3π/2 since the behavior is complex (requiring lengthy fitting functions) while the actual contribution of R 3 to D s is negligible in this angular range. Though not strictly correct, it is convenient to set R 3 = 0.55 at these angles to ensure continuity of R 3 at φ = π/2 and 3π/2. Finally, we have also fit A 60 : which (together with R 0 , R 1 , R 2 and R 3 given above) is required to calculate R 4 and R 5 by Eqs. (A.29a), (A.29b) and, in turn, D os by Eq. (A.28).

Appendix D: Spacecraft Shadowing
One can forward model the expected N 2 count rate by entering a reasonable density vs. altitude function into Eqs. (A.22a) and (A.30). We use here n(x) = C × exp(−x/H ) with fitting parameters C and H , where the scale height H is approximated as constant [as adequate for our purposes, Westlake et al. 2011]. However using the T85, 88 and 96 inbound densities to constrain the fit of X s and X os , we found that the model over-estimates the signal after θ exceeds ∼90 deg by a factor of ∼6 in the worst case (Fig. 19). For T47 the deviation from the model is comparatively minor (Fig. 19), and only after θ exceeds ∼110 deg. The explanation is the actuator angle of the rotatable Cassini Plasma Spectrometer (CAPS), which lies beneath INMS on the spacecraft deck. Contrary to T47 when CAPS was actuated to a rotation angle of 0 deg for most of the encounter (+/−10 min), the rotation angle was near 90 deg (Fig. 20) during the T85, T88 and T96 flybys, which was the final position when CAPS went offline on June 2, 2012 prior to T85. At 90 deg rotation angle the CAPS instrument protrudes farther into the INMS ram line of sight shortly after θ passes 90 deg (Fig. 20), i.e., approximately the same time as the N 2 signal drops significantly below the expected level.
In Fig. 5 we show the θ , ϕ ram directions shadowed by the spacecraft, overlaid onto an estimate of the CSN ram factor D s for S 0 = 20 (typical of a Titan flyby). One can see for the most part the spacecraft shadow zone lies (by design) in the areas of low INMS sensitivity, i.e., far off ram. We do not currently include the effect of spacecraft shadowing in our ram factor modelling, i.e., in the estimates of the D factors. This would require multiplication of F and F by spacecraft correction factors, obtained from modeling of (1) direct shadowing of the ram flux, and (2) secondary scattering from other spacecraft surfaces into the INMS. Direct shadowing becomes important for pointing within a few degrees of the spacecraft shadow zone, due to the small thermal spread of the incoming ambient molecule velocity vectors about the ram direction (assuming S 0 1). Inside the shadow zone (i.e., ram from behind the spacecraft), secondary scattering from the spacecraft becomes the dominant contribution to the signal. Inclusion of the spacecraft correction factors is of course a possibility for future work, which would extend the range of validity of the ram factor modelling to all 4π steradians of θ and ϕ. However, with the caveat that S 0 should be significantly greater than unity (essentially always true), we recommend for now that the ram factor modelling be used for angles for which no more than 5 % of the ram flux is shadowed by the spacecraft. These angles are shown in Fig. 6. For a typical Titan flyby S 0 of 20 this requirement translates to a minimum angle of ∼3 degrees of the pointing from the edge of spacecraft.

Appendix E: Comparison to Original 'Ideal' Calibration Model
The 'ideal' closed source that was the basis of the original calibration model was approximated as completely sealed, i.e., with the exception of the CS inlet as the only gas conduction pathway. With this assumption the analysis of Appendix A yields which is independent of the flux into the vent or OS collimator. The approximation A 00 ≈ A 01 is made to the CS inlet since the combined shape is conductance limited mainly by the tube, which receives an isotropic molecular velocity distribution at both ends. This approximation breaks down for ram directions near the special angle [at θ = 13°, ϕ = 180°, see  where the ram molecular beam is aimed from the entrance aperture directly into the tube (Fig. 5), i.e., without antechamber wall collisions to pre-isotropize the molecule velocities. However the original calibration model did not consider this case. Therefore Eq. (E.1) becomes with the CSN detector count rate X for species f (Eq. (A.21)) given by and the correction factor to the ambient neutral densities from the original calibration is F/D s times the reduction factor β = 1.55 ± 21 % (Sect. 4): with D s from Eq. (A.18). One can see by comparing Eqs. (E.2) and (A.17) that the ram factor D s is equivalent to F in the original calibration, unlike the revised ram factor model presented in this paper. The angular θ and φ dependence of the correction at S 0 = 20 is shown in Fig. 10. The factor exceeds unity because the new INMS calibration model is less sensitive than the original: the result of (i) the s f sensitivity reduction, and (ii) the gas leakage out of the CS now accounted for. Analyzing 53 Titan flybys, we find that the closest approach densities are on average a factor of 2.2 greater with the new INMS calibration than with the original calibration. This value is equal (within error) to the limiting S 0 1 (supra-thermal spacecraft speed) value of the θ = 0 (on ram) correction factor: also 2.2 ± 23 %. This constant is useful (in CSN mode) for re-interpreting densities evaluated with the original calibration for flybys fixed near 0 degrees ram. Unlike the neutral densities obtained in CSN and OSNT modes (subject to the ram factor modeling discussed above), we emphasize that the INMS ion and neutral beam measurements (OSI and OSNB modes) do not involve wall collisions and thermalization of the molecules in the instrument, and are therefore only affected by the (angle independent) s f sensitivity correction constant discussed in Sect. 4. Many recent INMS publications have adopted a constant correction factor of 3 for the CSN neutral densities, in anticipation of the more detailed assessment given in this paper. This factor is now superseded by the more accurate value of 2.2 ± 23 % (which includes ram factor and s f sensitivity corrections) for on-ram flybys, or, as appropriate, the angle dependent correction. Likewise, a constant factor 1.55 ± 21 % (i.e., s f sensitivity correction only, no angle dependence) should now be adopted for OSI and OSNB. In the S 0 1 limit considered here the correction factors are independent of the gas species since all gases enter the INMS with equivalent (near perfect) collimation. The factors are applicable in the 'collisionless' regime below ∼4 × 10 −9 kg/m 3 (∼8 × 10 16 molecules/m 3 ) in which ambient mean free paths are larger than the spacecraft in all directions (including the direction of motion); a condition satisfied at all Titan altitudes traversed by Cassini.

Appendix F: Monte Carlo Modeling
Here we give a more detailed description of the Monte Carlo modeling of the transmission factors used to calculate effective areas in the INMS ram factor modeling. We estimate the CS and vent transmittances t 0 and t 5 , as well as transmittances t 4 between the ambient environment and ISE through OS collimator slits, t 6 between ambient and the OS, and t 8 between the OS and ISE through the OS-aperture and OS collimator slits, with separate Monte Carlo simulations (Table 1). All pathways except 8 are exposed to the ram flux, and therefore we estimate the angle/speed dependent ingoing and outgoing transmittances (i.e., t i and t i , respectively) with separate simulations. Only one simulation was required for pathway 8 since the effective areas are equal in both directions, as required for all internal pathways in thermodynamic equilibrium: hence A 8 = A 8 , i.e., t 8 = (a 8 /a 8 )t 8 (the same equality holds for internal pathways 2, 3 and 7).
To estimate the CS inlet and OS collimator inward (unprimed) transmission we initialized particles from randomly distributed points on their entrance aperture planes, with an isotropic Maxwell-Boltzmann velocity distribution. We shifted the distribution in velocity space by the spacecraft velocity vector, and scaled it by velocity, i.e., f (υ) → υf (υ) such that particles are initialized according to the arrival rate at the aperture (Greenwood 2002). By contrast we subdivided the vent aperture into 20 flat sub-planes to approximate its curved surface, evaluating the number of molecules and the velocity vector distributions for each individual plane according to the portion υf (υ) within the velocity space hemisphere seen by that plane. The ingoing directional distribution is completely characterized by S 0 , θ and φ, and therefore t i = t i (S 0 , θ, φ), where i = 0, 4, 5 or 6 ( Table 1). We allowed molecules to bounce randomly inside the structure in question, and molecules striking the walls were re-emitted according to a cosine angular distribution with respect to the surface normal (Steckelmacher 1986;Wenaas 1971). We did not consider intermolecular collisions, since the size scales are much smaller than the mean free paths at typical densities encountered by INMS. We estimated the transmission probability as the ratio of the number of particles transmitted through the exit areas (i.e. slits or apertures or both), to the number created at the initialization area. The statistical error is the inverse square root of the total number of transmitted particles.
Likewise we estimate the outgoing transmissions by initializing molecules at random positions on the exit area(s) with a cosine directional distribution, as appropriate for an isotropic flux transiting a plane (Greenwood 2002). Unlike t i , the t i are constant numbers, with no dependence on S 0 , θ or φ, since we are now considering only the escape of the already thermally equilibrated gas. We follow the same procedure as for t i , estimating t i as the ratio of molecules transmitted through the entrance aperture back into space, to the number entering the structure from the exit area.
In the CS inlet simulation we consider the antechamber and the tube together, as shown in . We simulate the OS collimator as a short tube representing the entrance aperture, in series with a tube representing the collimator assembly. The exit areas consist of the four axial collimator slits, the circular inner and outer ring slits, and the exit aperture (see Fig. 3 detail on OS collimator). We consider the slit depth by assigning a transmission factor to each slit, given by the Clausing factor for the slit of the given width and depth (Dushman and Lafferty 1962). Particles that fail to transmit through the slit are re-directed isotropically back into the collimator assembly from the location where they struck the slit (we ignore 'lateral drift' of the particles in the slit, since for a narrow slit the effect is small). Particles re-entering the collimator assembly from the ISE must also pass through the slits, with the possibility of being rejected back into the ISE.
An equally valid alternative method of estimating transmission factors is a surface mesh approach, i.e. to define a shape's two-dimensional transfer matrix M kl describing the vapor transfer rate between the mesh elements which depends on element area, inter-element distance and relative orientation. The instantaneous molecule surface coverage σ k is modified iteratively with a time step t : σ k → σ k + t l M kl σ l . The transmission factor is estimated from the ratio of the outgoing to ingoing flux after the number of iterations is sufficient to reach steady state, i.e., σ steady = lim N→∞ ( ↔ I + t ↔ M) N σ . A detailed explanation is given by , who applied the surface mesh approach to simulate the propagation of the 'sticky' molecule H 2 O through the CS inlet. However the model is equally applicable to 'non-sticky' volatiles, and for example our calculations of CS inlet transmission t 0 with the surface mesh approach are in excellent agreement with the Monte Carlo simulations.

Appendix G: Thruster Spikes and Inter-molecular Collisions
Negative N 2 spikes during the T85 and T88 encounters (Figs. 16 & 17 for T85) were measured simultaneously with the firing times of Cassini's reaction control thrusters, which were used (instead of Cassini's three reaction wheels) to stabilize the spacecraft orientation against atmospheric drag. Each thruster had a different effect on the signal, with the most intense negative spikes produced by the Y3 and Z3 thrusters. As shown in Figs. 16 & 17 for T85, the negative spikes gained intensity until θ exceeded 90 deg (at ∼−140 and +60 sec for T85 & 88, respectively), whereupon the spikes abruptly became positive. We performed collisional Monte Carlo simulations of the exhaust plume to investigate the cause of the negative spikes, taking into account the plume composition (assuming as an approximation complete stoichiometric conversion of the N 2 H 4 propellant to ∼2:1 H 2 :N 2 ), mass (∼0.5 g/s expelled from the thrusters), and thrust [∼0.69 N Rizvi 2014], and intermolecular collisions with the ambient atmosphere streaming past the spacecraft. We simulated the ambient atmosphere as a drifted 150 K Maxwellian distribution of N 2 molecules, the spacecraft body as cylindrical with the conical high-gain antenna (HGA) at one end, and set the thruster positions and exhaust plume directions accordingly with respect to the HGA and INMS position on the spacecraft body. Thermal accommodation coefficients near unity are routinely observed by Earth-orbiting satellites , and we therefore assumed as an approximation full thermal accommodation of the impacting molecules, and scattered them from the spacecraft according to a 'diffuse' cosine ejection angle distribution. Additionally, we monitored the arrival angles of molecules at the INMS instrument position, and convolved onto the angle-dependent modelling described in the preceding sections to predict the change of CSN and OSNT count rates due to thruster firings.
The simulations demonstrated that the exhaust plume N 2 density (e.g. ∼7 × 10 17 N 2 /m 3 , or ∼3 × 10 −8 kg/m 3 , two meters from the thruster) is sufficient for scattering of a major Here the simulation assumes a 3000 K thruster exhaust temperature with a 0.5 g/s burn rate and 1500 m/s ejection speed (yielding 0.69 N of thrust) fraction of the incoming molecules away from the INMS by thrusters Y3 and Z3. The explanation is the near alignment of INMS's location along the Z axis with thruster block 3, which fires directly into the ram flux line of sight of the INMS vent (Figs. 16 & 21). However all three inlet pathways (CS, OS, and vent) are pointed obliquely to the ram direction for oblique θ , and hence the instrument sees only a minor fraction of the ram flux after θ exceeds 90 deg. Therefore the thruster spikes become positive after θ passes 90 deg since, as shown in Fig. 16, obstruction of the (small) direct flux by the exhaust plume is less significant than scattering of ambient N 2 toward the inlet pathways.
The effectiveness of scattering (as much as 40 % by the Y3 & Z3 thrusters, Fig. 22) raises the question of whether the drag and/or torque on the spacecraft is significantly reduced by intermolecular collisions during thruster firing events. However our simulations indicate that the 'shadow' cast by the thruster plume (Fig. 23) covers only a small fraction of the spacecraft in the worst case, and often the ram direction does not aim the scattering shadow Fig. 23 Estimate of the ram enhancement factor (cross-section shown) vs. position (grid lines: 1 meter) from the thruster exhaust simulation, with the thruster plume parameters of Fig. 22, for the T85 ram direction (from below) and estimated ambient density (∼6.7 × 10 14 N 2 /m 3 ) at −150 sec from closest approach to Titan, i.e., around the time when INMS observed negative N 2 spikes. Cassini is simulated as a cylinder plus a cone for the high gain antenna. The ram factor shown here includes incident and scattered ambient N 2 , but excludes thruster exhaust. Backscattered N 2 molecules produce drastic density ram enhancements near the spacecraft, except for the density wake downstream of the high gain antenna. Diamonds: thruster and INMS positions. The reduction in density near the thruster occurs as ambient molecules are scattered away by the exhaust plume. A slight "shadow" is also cast downstream of the thruster plume toward the INMS at the spacecraft at all. Moreover, the shadowing tends to be largely mitigated since every scattering event produces a recoil N 2 which has a good chance of striking the spacecraft. Clear evidence for recoils in the exhaust plume is seen in the positive H 2 spikes (Fig. 17), which increase from <10 counts per IP (below the noise) to several hundred near 0 sec, as the higher atmospheric densities at closest approach scatter more thruster H 2 back toward the INMS. During thruster firings we find that the thruster collisional correction to drag and torque is at most a few percent in the worst case, and therefore not a serious problem for AACS-based density estimates, nor a significant contributor to the difference with INMS.
A related consideration is the effect of inter-molecular collisions generally, i.e., without thruster plumes, since the mean free path as viewed in the spacecraft reference frame shortens to only ∼10 meters at Titan closest approach densities (∼10 16 molecules/m 3 , or ∼5 × 10 −10 kg/m 3 , at 1000 km altitude) for molecules forward-directed into the incoming ram flow. For example ambient N 2 forward scattered off the spacecraft can re-collide in the spacecraft vicinity at closest approach densities, thereby slowing or scattering away a fraction of the incoming N 2 . However our simulations showed that the effect is small even at closest approach, and largely canceled by recoiling N 2 re-striking the spacecraft. The simulations show a drop in drag coefficient at high density through conversion of the ram dynamic pressure into heat as the incoming gas slows in front of the spacecraft, but the effect is small, e.g., with only a ∼3 % drop predicted at 2 × 10 −9 kg/m 3 (∼4 × 10 16 molecules/m 3 ): approximately the maximum Titan density encountered in the mission (Figs. 7-9), imply-ing that the gas flow around the spacecraft is essentially collisionless. Moreover, significant collisional flow stagnation in front of the spacecraft would change the derived INMS and AACS densities in different proportions if intermolecular collisions were not accounted for, since the gas flux and drag force seen, respectively, by INMS and AACS, differ by a factor of (gas) velocity. Accordingly, the proportionality of the INMS and AACS densities over the entire range of Titan flyby densities (Fig. 7) indicates a constant drag coefficient within error at the relevant densities. We therefore conclude as expected that inter-molecular collisions cannot be a significant source of error between AACS and INMS.