Quantum mechanical calculation is performed for the $O(^1D)$ + HCl ${ ightarrow}$OH + Cl reaction using Reactive Infinite Order Sudden Approximation. Shifting approximation is also employed for the l ${ eq}$ 0 partial wave contributions. Various dynamical quantities are calculated and compared with available experimental results and quasiclassical trajectory results. Vibrational distributions agree well with experimental results i.e. product states mostly populated at $v_f$ = 3, 4. Our results also show small peak at $v_f$ = 0, which indicates bimodal vibrational distribution. The results show two significant broad peaks in ${gamma}_i$ dependence of the cross section, one is at ${gamma}_i$ = $15^{circ}-35^{circ}$ and the another is at ${gamma}_i$= $55^{circ}-75^{circ}$ which can be explained as steric effects. At smaller gi, the distribution is peaked only at higher state ($v_f$ = 3, 4) while at the larger gi, both lower state ($v_f$ = 0) and higher state ($v_f$ = 3, 4) are significantly populated. Such two competing contributions (smaller and larger ${gamma}_i$) result in the bimodal distribution. From these points we suggest two mechanisms underlying in current reaction system: one is that reaction occurs in a direct way, while the another is that reaction occurs in a indirect way.