Global potential energy surfaces (PESs) of the [H,C,N,O] system in singlet and tripler states have been investigated using the hybrid density functional B3LYP/6-311G(d,p) method. Isocyanic acid, HNCO 1, has been found to be the most stable isomer for both multiplicities. The adiabatic singlet-triplet splitting for 1 is 82.6 kcal/mol. In the singlet state, HNCO is energetically followed by cyanic acid, HOCN 2, 28.7 kcal/mol higher than 1, fulminic acid, HCNO 3 (67.9 kcal/mol), and isofulminic acid, HONC 4 (87.1 kcal/mol). In the triplet state, the branched NC(H)O isomer (3)7 is 0.3 kcal/mol higher than (3)1, followed by HOCN (3)2 (27.9 kcal/mol relative to tripler HNCO) and HCNO (3)3 (40.6 kcal/mol). The barriers for intramolecular rearrangements within singlet and triplet [H,C,N,O] system have been calculated to be high, and the isomerization processes in most cases are not expected to compete with fragmentations. Several minima on the singlet-triplet seam of crossing, relevant to the singlet [H,C,N,O] decomposition reactions, have also been found. The global features of the singlet and tripler PES have been applied to several important reactions, such as NH((3) Sigma(-))+CO, thermal decomposition of HNCO, O(P-3)+HCN, O(P-3)+HNC, and CH((2) Pi)+NO((2) Pi). For these reactions, major product channels have been speculated and their activation energies have been reported. Adiabatic ionization potentials for singlet and tripler [H,C,N,O] have been found to be high, in the range of 180-270 kcal/mol.