We present a theoretical analysis of formation and kinetics of hot OH molecules in the upper atmosphere of Mars produced in reactions of thermal molecular hydrogen and energetic oxygen atoms. Two major sources of energetic O considered are the photochemical production, via dissociative recombination of O$_{2}^{+}$ ions, and energizing collisions with fast atoms produced by the precipitating Solar Wind (SW) ions, mostly H$^+$ and He$^{2+}$, and energetic neutral atoms (ENAs) originating in the charge-exchange collisions between the SW ions and atmospheric gases. Energizing collisions of O with atmospheric secondary hot atoms, induced by precipitating SW ions and ENAs, are also included in our consideration. The non-thermal reaction O + H$_2(v,j) \rightarrow$ H + OH$(v',j')$ is described using recent quantum-mechanical state-to-state cross sections, which allow us to predict non-equilibrium distributions of excited rotational and vibrational states $(v',j')$ of OH and expected emission spectra. A fraction of produced translationally hot OH is sufficiently energetic to overcome Mars' gravitational potential and escape into space, contributing to the hot corona. We estimate the total escape flux from dayside of Mars for low solar activity conditions at about $5\times10^{22}$ s$^{-1}$, or about 0.1\% of the total escape rate of atomic O and H. The described non-thermal OH production mechanism is general and expected to contribute to the evolution of atmospheres of the planets, satellites, and exoplanets with similar atmospheric compositions.

The precipitation of energetic neutral atoms, produced through charge exchange collisions between solar wind ions and thermal atmospheric gases, is investigated for the Martian atmosphere. Connections between parameters of precipitating fast ions and resulting escape fluxes, altitude-dependent energy distributions of fast atoms and their coefficients of reflection from the Mars atmosphere, are established using accurate cross sections in Monte Carlo simulations. Distributions of secondary hot atoms and molecules, induced by precipitating particles, have been obtained and applied for computations of the non-thermal escape fluxes. A new collisional database on accurate energy-angular dependent cross sections, required for description of the energy-momentum transfer in collisions of precipitating particles and production of non-thermal atmospheric atoms and molecules, is reported with analytic fitting equations. 3D Monte Carlo simulations with accurate energy-angular dependent cross sections have been carried out to track large ensembles of energetic atoms in a time-dependent manner as they propagate into the Martian atmosphere and transfer their energy to the ambient atoms and molecules. Results of the Monte Carlo simulations on the energy-deposition altitude profiles, reflection coefficients, and time-dependent atmospheric heating, obtained for the isotropic hard sphere and anisotropic quantum cross sections, are compared. Atmospheric heating rates, thermalization depths, altitude profiles of production rates, energy distributions of secondary hot atoms and molecules, and induced escape fluxes have been determined.