## 3. Formulation of the problemWe are interested here in the distribution of oxygen of interstellar origin in the computational domain shown in Fig. 1. The flow of atoms must be described kinetically, because, as already said, the free path length of the neutral atoms is of the order of the size of the interaction region. In order to obtain the kinetic distribution function, Boltzman's equation must be solved: where is the oxygen atom distribution function, the local Maxwell distribution function of the protons with the hydrodynamical values , , from the model described above; is the hydrogen atom distribution function, the local Maxwell distribution function of oxygen ions with the hydrodynamic values described later, the individual velocities of oxygen atoms, protons, oxygen ions and hydrogen atoms, the charge exchange cross-section of an O-atom with a proton, the charge exchange cross-section of an hydrogen atom with an oxygen ion, the photoionization rate, the mass of the oxygen atom, the solar gravitational force. Eq. (1) takes into account the following processes in which oxygen is involved: a) the solar gravitation: the potential energy is , where , is the universal gravitational constant, and - the mass of the Sun. b) the charge exchange of O-atoms with protons:
, with a charge exchange rate
, given by formula:
,
where
is the relative atom-proton velocity, and
is the charge exchange cross section.
,
are constants,
,
(Stebbings & al, 1964, Banks &
Kockarts, 1973). The relative velocity is measured in
c) the reverse charge exchange of hydrogen atoms with oxygen ions: . The cross section for this reaction is connected with the cross section of the direct charge exchange by the following relation: d) the photoionization: the rate is given by: (Banks and Kockarts, 1973), where is the photoionization rate at the Earth's orbit and - is 1 A.U. Collisions between oxygen atoms and electrons are not taken into
account in the present model. These collisions may be significant in
the heliosheath where solar wind electrons are hot, and may ionize
substantially the oxygen. Following the arguments of Fahr, Osterbart
and Rucinski (1995) and the Lotz (1967) formula, we have estimated
this effect for an interstellar plasma density of 0.1 cm-3. An upper
limit of 9% extinction between the termination shock and the
heliopause can be derived, assuming a 2.10 Eq. (1) contains the oxygen ions distribution function . This distribution function obeys a relation of the type: which is expressed according to the plasma structure. Actually, we do not calculate explicitly the distribution function , but its first moments, i.e. the number density of oxygen atoms , the bulk velocity , the 'temperature' (average kinetic energy) , and the pick-up ions source term defined by: In the last definition, is the number of newly created ions per time unit, or equally the number of photoionized or charge- exchanged oxygen atoms per time unit. Assuming that plasma picks up the new ions instantly, ionized atoms acquire immediately the velocity and the temperature of the solar wind. In these conditions, the number density of ions obeys the continuity equation: The boundary condition for the distribution function is the Maxwell distribution: (where ). The number density of oxygen ions at the outer boundary is determined by the ionization balance which prevails in the unperturbed medium: © European Southern Observatory (ESO) 1997 |