با استفاده از نرم افزار گوسین ساختارهای ممکن برای برهمکنش های بین مولکولی، مولکول های hclcs و c6h5x(x=f,cl,br,ch3,h) در نظر گرفته شدودرادامه اثر تغییر استخلاف x در مولکول c6h5x بر قدرت پیوند در کمپلکس های به وجود آمده مورد مطالعه قرار گرفت سپس هر یک از ساختارهای پیش بینی شده فوق در چهارسطح محاسباتی mp2/6-311++g(d,p)،mp2/6-311++g(2d,2p) ، b3lyp/6-311++g(d,p) و b3lyp/6-311++g(2d,2p)بهینه گردید و بر روی ساختارهای بهینه شده محاسبه فرکانس مربوطه انجام گرفت. در کمپلکس های c6h5x…hclcs انرژی پایداری برای هر x به ترتیب زیر افزایش می یابد: x=f s11< s14< s13<s12 x=cl s21, s23< s24<s22 x=br s32< s34 <s33< s31 x=ch3 s41 < s43< s5< s44<s42 در نهایت با استفاده از تئوریaim نوع برهمکنش ها مشخص گردید.