Recent ALMA (Atacama Large Millimeter/submillimeter Array) observations of young protostellar objects detected warm SO emission, which could be associated with a forming protostellar disk. In order to investigate if such warm gas can be produced by accretion shock onto the forming disk, we calculate the sputtering and thermal desorption of various grain surface species in one dimensional shock waves. We find that thermal desorption is much more efficient than the sputtering in the post-shock region. While H$_{2}$O can be thermally desorbed, if the accretion velocity is larger than 8 km s$^{-1}$ with the pre-shock gas number density of 10$^{9}$ cm$^{-3}$, SO is desorbed, if the accretion velocity $\gtrsim$ 2 km s$^{-1}$ and $\gtrsim$ 4km s$^{-1}$, with the pre-shock density of 10$^{9}$ cm$^{-3}$ and 10$^{8}$ cm$^{-3}$, respectively. We also find that the column density of hydrogen nuclei in warm post-shock gas is $N_{{\rm warm}} \sim 10^{21}$ cm$^{-2}$.

We performed one-dimensional hydrodynamic simulations with detailed cooling, heating and chemical processes to examine the thermal stability of shocked gas in cold neutral medium (CNM) and molecular clouds. We find that both CNM and molecular clouds can be thermally unstable in the cooling layer behind the shock wave. The characteristic wavelength of the thermal instability ranges from 10^-5 pc to 0.1 pc in the CNM, and from 10^-7 pc to 0.1 pc in the molecular clouds. This coincides with the size of observed tiny scale structures in the CNM and molecular clouds, indicating that the thermal instability in the shocked gas could be a formation mechanism of these tiny structures in the interstellar medium. We have also calculated the e-folding number of the thermal instability to estimate the amplification of the density fluctuation in the shocked gas. Density perturbations in the CNM grow by a factor of exp(5)~150, whereas the perturbations in the molecular clouds grow only by a factor of a few behind a high Mach number shock. The amplification factor is larger at lower densities and higher velocities. Formation of very small-scale structures by thermal instability in shocked gas is more effective in lower densities.

We study the evolution of phosphorous-bearing species in one-dimensional C-shock models. We find that the abundances of P-bearing species depend sensitively on the elemental abundance of P in the gas phase and on the abundance of N atoms in the pre-shock gas. The observed abundance of PN and the non-detection of PO towards L1157 B1 are reproduced in C-shock models with shock velocity v=20km s-1 and pre-shock density n(H2) =10^4 - 10^5, if the elemental abundance of P in the gas phase is 10^-9 and the N-atom abundance is n(N)/nH - 10^-5 in the pre-shock gas. We also find that P-chemistry is sensitive to O- and N-chemistry, because N atoms are destroyed mainly by OH and NO. We identify the reactions of O-bearing and N-bearing species that significantly affect P chemistry.