در این پایان نامه، جزییات سازوکار واکنش های o2h + so1، 3o + s3، 3o + 2s3، 3o + 2so و 2so + s3 روی سطح های انرژی پتانسیل حالت پایه و (یا) حالت برانگیخته از دیدگاه ترمودینامیکی و سینتیکی بررسی شده است. ساختار همه ی گونه های درگیر در این واکنش ها با روش dft در سطح محاسبه ای b3lyp/6-311+g(3df) یا روش mp با سطح محاسبه ای mp2/6-31+g(d)با نرم افزار گوسین 2003 بهینه شده است. برای افزایش دقت، انرژی تمام گونه ها در این واکنش ها در سطح های محاسبه ای همانند، ccsd(t)،g3b3 یا g3mp2b3 با استفاده از ساختارهای بهینه شده ی مدل های b3lyp یا mp2 به دست آمده است. برای هر واکنش، در ابتدا کمپلکس های حاصل از برخورد مناسب بین واکنش دهنده ها روی سطح های انرژی پتانسیل حالت پایه و (یا) حالت برانگیخته تعیین شد. سپس، با شروع از این کمپلکس ها ی برخورد و عبور از حدواسط ها و حالت های گذار، مسیر هایی برای انجام هر واکنش و رسیدن به محصول های ممکن مورد بررسی قرار گرفت. با استفاده از این حدواسط ها و حالت های گذار، نمودار انرژی-مسیر انجام هر واکنش روی سطح های انرژی پتانسیل حالت پایه و (یا) حالت برانگیخته رسم شد. پارامتر های ترمودینامیکی برای کانال های هر واکنش در دمای معمولی محاسبه شده و در نهایت، با استفاده از نرم افزار هایgpop وssumes ، ثابت سرعت برای مطلوب ترین مسیر هر یک از واکنش های یاد شده در یک گستره ی دمایی محاسبه شده است.