We studied the critical ionization velocity (CIV) process by first investigating a coupled system of equations describing the production of several ion species and electrons by impact ionization, their collisions with neutrals, and the heating of electrons. Analytic relations derived from this were tested with the help of a particle simulation including collisonal processes between neutrals and plasma particles. It was found that resistive heating of electrons plays an important role when the density of the neutrals is high and that electron heating due to lower hybrid waves is significant when the neutral density is low. In both cases we verified the control of the plasma production rate by the ratio of the beam velocity to the critical velocity. Himreel et al., 1976; Axniis, 1978; Mb'bius, 1979; Venkatararnani and Mattoo, 1980;Brenning, 1981], and related phenomena were reported in a rocket experiment [Haerendel, 1982] in relation with the space shuttle [Papadopoulos, 1984; Sasaki et al., 1986] and in some other phenomena in space [Lindeman et al., 1974; Cloutier et al., 1978; Petelski, 1981; Forrnisano et al., 1982; Galeev and Chabibrachmanov, 1983-]. However, the critical ionization velocity (CIV) process itself is not easily explained, since the neutral particles cannot be ionized by collisions with either background ions or electrons when the relative velocity is only slightly larger than V• [Danielsson, 1973].Early attempts to understand this phenomena focused on the interpretation of laboratory experiments and did not specify the basic mechanism underlying the phenomena, although the importance of some anomalous electron heating was recognized. A review of these theoretical considerations is given by Sherman [1973]. Relevance of the lower hybrid wave to the CIV process was postulated and studied by several authors [Raadu, 1978;Galeev, 1981;Abe, 1984;Papadopoulos, 1985]. According to them, it is possible to heat the electrons to energies roughly equal to the kinetic energy of the flowing neutrals. These energized electrons ionize the neutrals, producing new plasma to act as an additional generator of waves and thus forming a positive feedback loop.Investigations of the electron heating in the nonlinear stage of the lower hybrid wave instability were performed by computer simulations [McBride et al., 1972; Tanaka and Papadopoulos, 1983;Abe and Machida, 1985]. Moreover, in order to study the feedback mechanism directly, Machida et al. [1984] In this study we are concerned with the CIV process in both high and low neutral density cases. It is necessary to include several collisional processes between plasma and neutral particles in a high neutral density case in addition to the impact ionization collisions treated previously by several researchers [Raadu, 1978;Galeev, 1981;Papadopoulos, 1985]. In general, elastic collisions between electrons and neutrals occur most frequently when the electron energy is in a range below several tens of electron volts, which is the energy consideration in this work. It shou...