تابع پتانسیل جفتی هیدروژن را با استفاده از انتگرال های برخورد کاهش یافته ویسکوزیته در فشار صفر محاسبه کرده و با پتانسیل شبه hfd منطبق کرده ایم. با استفاده از این پتانسیل ضریب دوم ویریال و خواص انتقالی شامل ویسکوزیته و هدایت گرمایی هیدروژن را در گستره ی وسیعی از دما و چگالی در توافق خوبی با تجربه به دست آورده ایم. ما هم چنین شبیه سازی دینامیک مولکولی (md) را برای به دست آوردن فشار، انرژی درونی و ضریب خود نفوذی سیال هیدروژن در دماها و چگالی های مختلف با استفاده از پتانسیل جفتی تکمیل شده با تصحیحات کوانتومی روش فیمن-هیبس(fh) انجام داده ایم. نتایج نشان داد که در محاسبه ی فشار در ? t?^*?4 در نظرگرفتن اثرات کوانتومی ضرورتی ندارد. نکته قابل توجه در این کار این است که عبارت سه ذره ای هاوسچیلد و و پراوسنیتز به عنوان تابعی از دما و چگالی گسترش داده شد و پتانسیل سه ذره ای ونگ و سادوس با پتانسیل شبه hfd دوذره ای برای پیش بینی بهتر خواص انتقالی و ترمودینامیکی هیدروژن بدون نیاز به محاسبات سه ذره ای هزینه-بر استفاده شده است. هم چنین شبیه سازی دینامیک مولکولی هیدروژن برای تعیین معادله حالت (eos) جدید برای هیدروژن در توافق خوب با تجربه(nist eos استفاده شده است.

نتایج شبیه سازی دینامیک مولکولی برای به دست آوردن فشار کاهش یافته سیال کربن دی اکسید با استفاده از پتانسیل شبه hfd دو ذره ای (معادله 3-19) و پتانسیل های کل که از طریق معادلات 3-34 ، 3-37 و 3-46 به ترتیب با استفاده از پتانسیل های سه ذره ای هاوسچیلد و پراسنیتیز [42]، وانگ وسادوس [43]، و گیوزمن و همکارانش [46]، به دست آمده انجام شده است. قردادهای معمول با استفاده از پارامترهای شبه hfd پتانسیل co2 – co2 ، برای دانسیته کاهش یافته (??*= ????3 )، دمای کاهش یافته (t*= kt/?? )، و فشار کاهش یافته (p*= p??3/?? ) استفاده شده اند (جدول 3-1). نتایج فشار کاهش یافته برای سیال دی اکسید کربن در انسامبل nvt ، در 1t*= (در حالت مایع) و در 3 ، 9/1 ، 6/1 t*= (در منطقه بخار و فوق بحرانی) و دانسیته های مختلف با داده های تجربی[53]، در شکل 4-1 مقایسه شده اند. همانطور که این شکل نشان می دهد، هر دو پتانسیل سه ذره ای وانگ و سادوس [43] و گیوزمن و همکارانش [46] اثرات کوچکی بر روی نتایج دو ذره ای دارند و انحراف آنها را از مقدار تجربی افزایش می دهد؛ اما پتانسیل چند ذره ای هاوسچیلد و پراسنیتیز [42] نتایج را کاهش داده و آنها را به تجربه نزدیکتر کرده است. همچنین انتظار می رفت و دیده شد که، پتانسیل های مختلف در دانسیته های پایین و دماهای بالا انحراف کمتری را نشان می دهند. ما همچنین مقادیر فشار محاسبه شده خود را با مقدار به دست آمده آن از پتانسیل مکان – مکان هریس و یانگ [48] مقایسه کردیم. این برهمکنش مکان- مکان با روش بار الکتروستاتیک ، برای دی اکسید کربن که بیش از ده سال قبل توسعه یافته است، روش قابل فهمتری است . این پتانسیل بین مولکولی به علت توافق بهتر آن با داده های تجربی [55] هنوز به طور وسیعی استفاده می شود. طبق شکل4-1، اگرچه مقادیر فشار محاسبه شده با استفاده از روش مکان- مکان هریس و یانگ [48] در دمای پایین (1t*=) در توافق بهتری با تجربه است اما مقادیر ما با استفاده از پتانسیل چند ذره ای هایوسچیلد و پرایوسنیتیز در دمای بالاتر توافق بهتری را نشان می دهد. طبق شکل 4-1، نتایج فشار با استفاده از پتانسیل کل با استفاده از پتانسیل چند ذره ای هاوسچیلد و پراسنیتیز (معادله 3-34 ) در توافق خیلی خوب با داده های تجربی است. بنابر این شبیه سازی دینامیک مولکولی برای تعیین معادله حالت دی اکسید کربن با استفاده از پتانسیل کل انجام شده. مقادیر فشار کاهش یافته دی اکسید کربن برای دماهای کاهش یافته از 1 تا 10 ( k2 /241 - 2410) و دانسیته کاهش یافته از 1/0 تا 0/1 (kg/m3 140 - 1400) شبیه سازی شده. بررسی دقیق تر داده های دینامیک مولکولی نشان داد که معادله حالت جدید از طریق معادله 4-1 (با در صد انحراف از میانگین 5%) بیان می شود که ممکن است به عنوان یک مرجع برای سیال دی اکسید کربن استفاده شود. که در آن gi وابسته به دانسیته است و از معادله زیر پیروی می کند: که در آن z نمایانگر gi است . مقادیر پارامترهای این معادله در جدول 4-1 داده شده است . پارامتر های h , k وl نیز وابسته به چگالی هستند و از معادله زیر پیروی می کنند: که در آن z نمایانگر h , kو l است. مقادیر پارامترهای این معادله نیز در جدول 4-1 داده شده است. معادله حالت جدید (معادله 4-1) با معادله حالت nist [53]، و معادلات حالت محسن نیا و همکارانش [49]، و دوان و هیو [50]، در چهار دمای کاهش یافته مختلف، در شکل 4-2 ، مقایسه شده اند. مشخص شده است که معادله حالت جدید توافق بهتری را با معادله حالت nist نسبت به بقیه نشان می دهد . 4- 2-1-2 ضریب خودنفوذی با استفاده از فرمول انیشتن [29]، ضریب خود نفوذی ،d ، از جابجایی مجذور میانگین وابسته به زمان به دست آمده است، که انتظار می رود در زمان های طولانی به رابطه خطی تبدیل شود. که در آن ri(0) و ri(t) موقعیت ذرات به ترتیب در زمان های t=0 و t هستند. ضریب خود نفوذی کاهش یافته محاسبه شده ، دی اکسید کربن از جابجایی مجذور میانگین وابسته به زمان با استفاده از پتانسیل های مختلف در 44/1t*= در دانسیته های مختلف با مقادیر تجربی اتس و همکارانش [10] که از طریق تکنیک انعکاس اسپین شیب پالس به دست آمده، در شکل 4-3 مقایسه شده است. همانطور که این شکل نشان می دهد، مقادیر خود نفوذی کاهش یافته شبیه سازی شده ما در توافق خوبی با تجربه است. همچنین مشخص می شود که پتانسیل های سه ذره ای سهم خیلی کوچکی در خودنفوذی دارند اما مقادیر شبیه سازی شده را برای توافق بهتر با تجربه افزایش می دهد. 4 – 2- 2 مخلوط های co2 -ar و co2 -ch4 شبیه سازی دینامیک مولکولی برای به دست آوردن فشارهای کاهش یافته مخلوط های ar- co2وch4 - co2با استفاده از پتانسیل کل و شبه hfd دو ذره ای انجام شده است. قراردادهای معمول همچنین برای دانسیته ، دما و فشار کاهش یافته با استفاده از پارامترهای پتانسیل شبه hfd ،co2 - co2 (جدول 3-1) استفاده شده است. برای مقایسه نتایج شبیه سازی شده ما با تجربه، برای فشارهای پایین و متوسط از بسط ویریال قطع شده بعد از جمله سوم استفاده کرده ایم [51]: با که در آن xi کسر مولی ترکیب در مخلوط و bii* و ciii* به ترتیب ضرایب دوم و سوم ویریال جزء خالص i هستند. در این معادلات ، 1 بجای آرگون یا متان قرار گرفته و 2 بجای دی اکسید کربن قرار گرفته است. پارامترهای c112* و c122* با استفاده از معادلات زیر محاسبه می شوند: مقادیر تجربی پارامترهای, b22*, b12* , b11* c111* و c222* برای مخلوط co2 -ar در 34/1t*= و مخلوط co2 -ch4 در 20/1t*= در جدول 4-2 [52]، ارائه شده اند. نتایج فشار شبیه سازی شده مخلوط های co2 -ar و co2 -ch4 در دماهای ثابت و در انسامبل nvt و سه کسر مولی و دانسیته های مختلف با استفاده از پتانسیل های کل و دو ذره ای با داده های تجربی (معادله 4-5) به ترتیب در شکل 4-4 و 4-5 مقایسه شده است. همانطور که این شکل نشان می دهد ، اگر چه هر دو پتانسیل های سه ذره ای هاوسچیلد و پراسنیتیز [42] و وانگ و سادوس [45] نتایج دو ذره ای را افزایش داده اند و آنها را به تجربه نزدیکتر کرده اند اما پتانسیل چند ذره ای هاوسچیلد و پراسنیتیز اثرات بزرگتری روی فشار های دو ذره ای دارند (باید یادآوری شود که پتانسیل چند ذره ای هاوسچیلد و پراسنیتیز فشارهای دو ذره ای co2 خالص را کاهش داد). برای مولکول های غیر قطبی انتظار می رود (نه الزاما همیشه) که پتانسیل های سه ذره ای سه دو قطبی (از قبیل پتانسیل وانگ و سادوس) نتایج دو ذره ای را افزایش دهد. به عبارت دیگر، پتانسیل های سه ذره ای معمولا پتانسیل های دو ذره ای را بهبود می دهند که نتایج را در کمتر از مقدار تجربی به دست می آورد که بستگی به موثر بودن پتانسیل دو ذره ای دارد (از قبیل مخلوط های co2 -ar و co2 -ch4 در شکل 4-4 و 4-5 ). اما مشخص شده است که پتانسیل هاوسچیلد و پراسنیتیز نتایج دو ذره ای سیستم های خالص را کاهش و سیستم های مخلوط را افزایش می دهد و آنها را به تجربه نزدیکتر می کند (شکل 4-1 ، 4-4 و 4-5). بنابر این پتانسیل وابسته به دانسیته ایزوتروپیک هاوسچیلد و پراسنیتز نه تنها برهمکنش سه دو قطبی را بلکه اثرات دیگر همچون جملات دافعه ی سه ذره ای یا اثرات چند ذره ای را نیز شامل می شود. این اساسا بخاطر پارامتر تنظیم پذیر ?? است که از طریق مقایسه با داده های تجربی به دست آمده است. بنابراین، شاید بهتر باشد که عبارت سه ذره ای هاوسچیلد و پراسنیتیز به عنوان پتانسیل چند ذره ای [42] بیان شود. همچنین انتظار می رفت که با استفاده از پتانسیل سه ذره ای وانگ و سادوس در غلظت های بالای co2 در مخلوط co2 - ar و در دانسیته های بالای آن میزان بهبود در توافق با تجربه بیشتر شود. اگر چه اثرات چند ذره ای از پتانسیل هاوسچیلد و پراسنیتیز در دانسیته های بالا افزایش یافته اما چنین اثر غلظتی در شکل 4-4 و 4-5 را نمی توان مشاهده کرد؛ که می تواند بخاطر این واقعیت باشد که اثرات چند ذره ای از پتانسیل هاوسچیلد و پراسنیتز (مخصوصا پارامتر ?? ) تنها به تفاوت بین نتایج دو ذره ای و تجربی مربوط می شود. ما همچنین نتایج فشار سیستم co2 -ch4 را با استفاده از عبارت سه ذره ای اوکلی و همکارانش (معادله 3-44) با مقادیر تجربی در شکل 4-5 مقایسه کردیم. برای استفاده از عبارت سه ذره ای اوکلی و همکارانش در شبیه سازی، در ابتدا مقادیر گزارش شده k به وسیله اوکلی و همکارانش، که با منطبق کردن kدر کسرهای مولی مختلف [20] به دست آمده را استفاده کردیم . همانطور که شکل 4-5 نشان می دهد، اگرچه پتانسیل سه ذره ای اوکلی و همکارانش نتایج را بهبود داده اند اما تاثیر خیلی کوچکی روی فشارهای دو ذره ای دارد. این ممکن است بستگی به نوع پتانسیل های دو ذره ای، یا رنج دما یا فشار استفاده شده برای تعیین مقادیر k ، که با کار ما فرق دارد، داشته باشد. بنابراین ما مقادیر k را برای سیستم co2 -ch4 با استفاده از پتانسیل های خود در دانسیته های مختلف در 34/1t*= ( به وسیله مقایسه بین مقادیر پیش بینی شده و مقدار فشار تجربی) محاسبه کردیم و در جدول 4-3 ارائه دادیم که در معادله زیر فیت کرده ایم: که در آن ci از معادله زیر پیروی می کند: که در آن y بجای ci و x1 کسر مولی متان است. مقادیر پارامترهای این معادله در جدول 4-4 داده شده است. همانطور که در شکل 4-5 دیده می شود، مقادیر فشار شبیه سازی شده ما با استفاده از پتانسیل سه ذره ای اوکلی و همکارانش [20] با مقادیر k از معادله (4-10) در توافق خوبی با تجربه است.

در بخش اول پایان نامه، ابتدا سنتز استر متیل پرفلوئورواکتانوات انجام شده است.در ابتدا،اسید پرفلوئورواکتانوات در حضور مقدار کاتالیستی از پیریدین، وارد واکنش با تیونیل کلرید شده و پس از جداسازی و خالص سازی لایه پایینی، پرفلوئورواکتانوئیل کلرید بدست می آید.سپس این محصول, در حضور مقدار کاتالیستی از پیریدین، وارد واکنش با متانول شده و پس از جداسازی و خالص سازی لایه پایینی، متیل پرفلوئورواکتانوات بدست می آید. در بخش دوم پایان نامه, ابتدا ساختار بهینه متیل اکتانوات و متیل پرفلوئورواکتانوات توسط نرم افزار gaussian 09w بدست آمده و توسط نرم افزار dl_poly 4.03 در دو دمای 300 و 350 کلوین، شبیه سازی دینامیک مولکولی متیل اکتانوات و متیل پرفلوئورواکتانوات انجام و مقادیر محاسبه شده دانسیته و آنتالپی این دو ترکیب با یکدیگر و با مقادیر تجربی مقایسه شده است.

سولفورهگزافلورید،دارای ساختاراکتاهدرال،غیرقابل انفجار و غیرقابل اشتعال می باشد.در حالت خالص غیرسمی است.بسیار پایدارو به شدت الکترونگاتیو وجاذب الکترون است. یکی از سنگین ترین گازهای شناخته شده وگازی بی رنگ و بی بو است. ویژگی دی الکتریکی فوق-العاده آن نتیجه ی خصلت الکترونگاتیوی مولکول های آن می باشد که تمایل مشخصی به گرفتن الکترون های آزاد دارد و به تشکیل یون های سنگین با تحرک کم می انجامد. قدرت دی الکتریکی آن 5/2 برابر بیشتر از هوا در شرایط یکسان می باشد. به علت دمای تفکیک پایین و انرژی تفکیک بالا به عنوان یک گاز سردکننده قوس بسیار عالی است. سرعت صوت در آن یک سوم هواست وآن را عایق صوتی خوبی می سازد. به دلیل امتیازاتی مثل قدرت دی الکتریک بالا با هزینه کمتر نسبت به عایق های مایع دیگر، توانایی خوداحیایی و افزایش فشار کمتر طی تفکیک، به طور وسیعی در صنایع الکتریکی (به عنوان قطع کننده مداربا ولتاژبالا ، ایستگاه های مخابرات، ترانسفورماتورها ، شتاب دهنده ها و خطوط انتقال و کابل هایولتاژ بالاوغیره) کاربرد دارد. همچنین از آن در مصارف غیرالکتریکی به طور وسیعی استفاده می شود. برای مثال در ریخته-گری منیزیم و آلومینیوم به ترتیب از آن به عنوان گاز محافظ و گاززدا، در صنعت نیمه هادی ها برای تولید تراشه ها و مدارهای الکتریکی، به عنوان ضربه گیر در کفش های آموزشی به دلیل اندازه بزرگ مولکول های آن، در پنجره ها به عنوان عایق صدا و گرما، در پزشکی، در میکروسکوپ های الکترونی، ژنراتورهای واندرگراف، رادارهای هواپیما و گاز ردیاب نیز کاربرد دارد. آلوده کننده های اتمسفری تولیدشده توسط فعالیت های انسانی با توجه به اثراتشان به دو دسته طبقه بندی می شوند: کاهش دهنده اوزون استراتوسفری و گرم کننده سطح زمین (اثرگلخانه-ای). سولفورهگزافلوئورید سهم عمده ای در کاهش اوزون استراتوسفری ندارد چون شامل کلر که عامل اصلی فعل و انفعالات در اوزون می باشد، نیست. همچنین با وجود اینکه عمر آن در جوطولانیست، به دلیل غلظت کم آن در اتمسفر، سهم آن در گرم کردن کره زمین خیلی ناچیز می باشد.[6-3] برخلاف اهمیت سولفورهگزافلوئورید و مخلوط های آن در بسیاری از زمینه های علوم و تکنولوژی،توجه کمتری به برهمکنش بین مولکولی مولکول های سولفورهگزافلوئورید و ساختاروخواص آن شده است و بررسی های بیشتری براساس مدل های پتانسیل مناسب برای آن لازم است. بنابراین هدف تحقیق حاضر، ارائه یک تابع پتانسیل جفتی دقیق برایسولفورهگزافلوئورید با روش وارونگی انتگرال های برخورد کاهش یافته در فشار صفر می-باشد. همچنین ضریب دوم ویریالسولفورهگزافلوئورید با استفاده از پتانسیل جفتی بدست آمده و در دماهای مختلف با تجربه مقایسه شده است. برخی خواص انتقالی سولفورهگزافلوئوریداز قبیل ویسکوزیته، هدایت گرمایی و ضریبخودنفوذی با استفاده از پتانسیل جفتی بدست آمده، در گستره وسیعی از دماها محاسبه شد. همچنین ازشبیه سازی md برای بدستآوردن بعضی خواص ترمودینامیکی(فشار)، ساختاری (تابع توزیع شعاعی)و انتقالی (ضریب خود نفوذی) سولفورهگزافلوئورید در دماها و چگالی های مختلف با استفاده از پتانسیل جفتی محاسبه شده، استفاده شد. همچنین فشار مخلوط های سیالsf6-ar وsf6-kr در دمای ثابت ودرچگالی های مختلف با استفاده ازمدل هایlike- hfd دوذره ای بدست آمد. به منظور در نظر گرفتن نیروهای سه ذره ای یک عبارت نیروی میانگین توسعه یافته هاسچایلد و پروسنیتذ [7]و پتانسیل سه ذره-ای وانگ و سادوس [8و9] همراه با پتانسیل های دوذره ای برای سیستم های سیال sf6وsf6-arوsf6-krبکار گرفته شد. همچنین شبیه سازی md سولفورهگزافلوئورید برای تعیین یک معادله حالت جدید استفاده شده است.

از دیدگاه تئوری، برای مولکول edpb چهل و شش پیکربندی پایدار انولی و نوزده پیکربندی پایدار کتو می توان در نظر گرفت. در میان پیکربندی های انولی، یازده پیکربندی درگیر در حلقه کی لیتی شش عضوی و پیوند هیدروژنی درون مولکولی هستند. حالت ترانس یک پیکربندی نیز درگیر در حلقه کی لیتی پنج عضوی و پیوند هیدروژنی ضعیفی در آن مشاهده می شود. بر طبق محاسبات نظری انجام شده، در تمام سطوح محاسباتی b3lyp و mp2 ، پیکربندیe1 همواره پایدارترین پیکربندی می باشد. همچنین مشخص شده است که بین پیکربندی های edpb اختلاف ناچیزی از نظر انرژی وجود دارد. با مقایسه مقادیر انرژی انول ها با فرم های کتو آن ها می توان نتیجه گرفت پیوند هیدروژنی باعث پایداری فرم انولی می شود. مقدار تفاوت انرژی بین پایدارترین پیکربندی با سایر پیکربندی ها (انرژی نسبی) در سطح محاسباتی b3lyp با افزایش تابع پایه بیشتر می شود. مقدار انرژی پیوند هیدروژنی محاسبه شده در سطح mp2 کمتر از مقدار محاسبه شده در سطح b3lyp است. همچنین این مقدار با تغییر تابع پایه از 6-311g بهg** ++6-311 کمتر می شود. این نتایج در توافق با فاصله o…o محاسبه شده است. با توجه به اطلاعات بیان شده می توان نتیجه گرفت دقت در انجام محاسبات با افزایش سطح و پایه بیشتر می شود. (سطح محاسباتی mp2 به خاطر درنظر گرفتن تصحیح انرژی تا مرتبه دوم نتایج دقیق تری از سطح محاسباتیb3lyp به دست می دهد.) قدرت پیوند هیدروژنی در edpb به اندازه kj/mol296/11 بیشتر از dmoa و kj/mol480/25 کمتر از ba است. مقایسه طول o…o در ترکیب edpb (535/2) و dmoa (621/2) و ba (520/2) همچنین بررسی زاویه oho در ترکیب edpb (722/147) با dmoa (981/144) و ba (997/148) نشان دهنده افزایش قدرت پیوند هیدروژنی edpb نسبت به dmoa و کاهش آن نسبت به ba است. بنابراین قدرت پیوند هیدروژنی ترکیبات را می توان به صورت زیر ترسیم نمود: edpb > dmoa<ba پارامترهای ساختاری edpbنیز پیوند هیدروژنی قویتری را در مقایسه با dmoa و پیوند هیدروژنی ضعیف تری را در مقایسه با ba نشان می دهند.

: ساختار مولکولی باتجزیه و تحلیل کامل کنفورمرهای انول و کتو و پیوند هیدروژنی درون مولکولی و نیز فرکانس های ارتعاشی مولکول (z)-4،6-دی اُکسو-6-فنیل هگز-2-اِنوئیک اسید توسط روش تئوری تابعی چگال (dft) مورد بررسی قرار گرفت. بعلاوه پارامترهای کنفورمرهای سیس- انول در سطوح محاسباتی mp2/6-31g**و b3lyp/6-311++g** و b3lyp /6-311g** بهینه شد. نتایج با مولکول بنزوئیل استن به عنوان مولکول مادر مقایسه گردید. علاوه بر ده کنفورمر انولی که بوسیله پیوند هیدروژنی درون مولکولی پایدار می گردند تعداد سه کنفورمر کتو نیز مشاهده گردید. محاسبات تئوری نشان می دهند که پیوند هیدروژنی در zdpea نسبتاَ ضعیفتر از بنزوئیل استن می باشد.

فشار و ضریب خودنفوذی سیال تترا فلورید کربن با استفاده از پتانسیل hfd-like دو ذره ای که از وارونگی انتگرال برخورد چگالی کاهش یافته در فشار صفر بدست آمده و همچنین فشار مخلوط های سیال cf_4-arcf_4-ch_4 و در دمای ثابت و در چگالی های مختلف با استفاده از توابع پتانسیل دقیق و جدید hfd-like دو ذره ای با استفاده از شبیه سازی دینامیک مولکولی محاسبه شده است. برای داشتن نیروهای چند ذره ای دقیق، پتانسیل های سه ذره ای هاوسچایلد و پراسنیتز(1993)، وانگوسادوس (2006)، و گیوزمن (2011) همراه پتانسیل های دو ذره ای مورد استفاده قرار گرفتند. اهمیت این کار در این است که پتانسیل چند ذره ای هاوسچایلد و پراسنیتز به عنوان یک تابع چگالی، دما، و کسرمولی به کار می رود و با پتانسیل های hfd-like دو ذره ای سیستم های cf_4-ch_4 ,cf_4-ar ,cf_4 برای بهبود بخشیدن پیشگویی مقادیر فشار بدون نیاز به محاسبات پر هزینه سه ذره ای استفاده می شود . شبیه سازی دینامیک مولکولی تترا فلورید کربن همچنین برای تعیین یک معادله حالت جدید مورد استفاده قرار گرفته است. نتایج توافق خوبی با مقادیر تجربی دارند.

نتایج ما نشان داد که پتانسیل جند ذره ای توسعه یافته ما از پتانسیل های سه ذره ای گیوزمن و همکاران(2011)و ونگ و سادوس (2006) بهتر است

تصحیحات کوانتومی با استفاده از روش فیمن- هیبس برای نئون جامد در نظر گرفته شده است. یک عبارت ساده و جدید با پتانسیل hfd-like بدون نیاز به محاسبات وقت گیر سه ذره ای به کار برده ایم. نتایج به دست آمده نشان می دهد که پتانسیل چند ذره ای و تصحیح یافته کوانتومی hfd-like پیش بینی نتایج دوذره ای کلاسیکی را بهبود می بخشد.

: در این تحقیق ما با استفاده از شبیه سازی دینامیک مولکولی،جذب مخلوط های مختلفی (کسرهای مولی مختلف)از آرگون وهلیم را برنانو خوشه های نقره،با شعاع های مختلف و در دماها وفشارهای متفاوت مورد بررسی قرار دادیم.در این کار ما از پتانسیل ساتن-چن برای توصیف برهمکنش های فلز-فلز وبرای سایربرهمکنش ها از پتانسیل لنارد جونز استفاده کرده ایم.به منظور ثابت نگه داشتن دما وفشاردر طول شبیه سازی از امکانات نرم افزار دی ال پلی بهره گرفته شده است.ترموستات برندسون در تمامی شبیه سازی ها با زمان آسایش1/0 برای ثابت نگه داشتن دما در سل شبیه سازی استفاده شده است. ما میانگین انرژی برهمکنش بین اتم های گاز ونانوخوشه ها، ثابت جذب وآنتالپی جذب را برای همه نانو خوشه های نقره بدست آورده ایم.نتایج نشان می دهد که ثابت جذب، میانگین انرژی برهمکنش وآنتالپی جذب با افزایش اندازه نانو ذرات کاهش می یابد.همچنین با کاهش اندازه نانو ذرات کسر بیشتری از اتم های گاز بر سطح نانو خوشه ها جذب می شوند.زیرا اتم های سطح در نانو خوشه ها در مقایسه با اتم های بالک انرژی پیوندی کمتری دارند بنابر این با کاهش اندازه نانو ذره انرژی چسبندگی آن کاهش می یابد ودر نتیجه ساختار ناپایدارتر می شود.پس اتم های سطح تمایل بیشتری برای جذب اتم های گاز و افزایش انرژی چسبندگی دارند.به همین دلیل با افزایش اندازه نانو خوشه ها، کاهش ثابت جذب متوقع خواهد بود

تابع انرژی پتانسیل جفتی جدید برای فلوئور از طریق وارونگی انتگرال های برخورد کاهش یافته ی ویسکوزیته تعیین گردید. سپس ضریب دوم ویریال، ویسکوزیته وهدایت گرمایی با استفاده از پتانسیل بدست آمده، محاسبه گردید که توافق خوبی با داده های تجربی ظاهر شد. ما همچنین پتانسیل سیستم f2-f2 رابااستفاده از متدهای ab initio درجهت گیری های متفاوت تعیین کردیم. همچنین شبیه سازی دینامیک مولکولی برای تعیین بعضی خواص سیال فلوئورازطریق وارونگی پتانسیل جفتی تکمیل شده بوسیله تصحیح کوانتومی فیمن- هیبس وپتانسیل ab initioومتد ccsd(t)/cc-pvdz انجام شد. برای بررسی برهمکنش های چند ذره ای، بدون نیاز به محاسبات سنگین، ازپتانسیل های هاوسچیلد وپراوسنیتز(1993) و ونگ وسادوس (2006) همراه وارونگی پتانسیل استفاده گردید. ما همچنین برهمکنش های سه ذره ای تریمر f2 را ازطریق محاسبات ab initio در جهت گیری های مختلف محاسبه کردیم. نتایج بدست آمده حاکی از آن است که پتانسیل های سه ذره ای خواص انتقالی وترمودینامیکی دو ذره ای را بهبود بخشید.

سقط خود بخودی طبق تعریف سازمان بهداشت جهانی ختم حاملگی یا تولد نوزاد کمتر از 500 گرم قبل از هفته بیستم حاملگی تعریف شده است .از جمله فاکتور های شناسایی شده در سقط آنومالی¬های کروموزومی والدین، فاکتورهای فیزیولوژیکی رحم مادر،عفونتها،فاکتورهای محیطی،مثل استعمال الکل و سیگار و علل ژنتیکی را میتوان نام برد. از جمله علل ژنتیکی میتوان به تاثیر برخی از پلی¬مورفیسم های ژنتیکی در این بیماری اشاره کرد. ژن vegf از جمله ژن¬های کلیدی و مهم در رگ¬زایی است. از طرف دیگر در طی تشکیل جفت رگ¬زایی انجام می¬گیرد. بنابراین تغییرات در فاکتور رشد رگ¬زایی اندوتلیالی می¬تواند عاملی در سقط باشد. اخیراً در بررسی¬ها نشان داده¬اند که میزان vegf در طی بارداری در زنان افزایش می-یابد و بلافاصله پس از اتمام بارداری کاهش می¬یابد. هدف از این تحقیق بررسی شش پلی¬مورفیسم ژن vegf و ارتباط آن دو با سقط جنین در زنان استان آذربایجان¬شرقی می¬باشد. به ¬این منظور از 50 زن مبتلا به بیماری سقط مکرر و 50 زن سالم بدون سابقه سقط نمونه خون جمع آوری شد. dna ژنومی از لوکوسیت¬های خون محیطی استخراج گردید. از روش sequencing محصولاتpcr جهت تعیین ژنوتیپ افراد بیمار و کنترل مورد استفاده قرار گرفت. نتایج بدست آمده نشان¬دهنده عدم تفاوت معنی¬دار در فراوانی ژنوتیپ rs3025033 ([ci], 0.296-1.661 p=0.43) و ([ci],0.769-5.419 p=0.14) rs3025035 وrs114074122 ([ci], 0.984-1.102 p=0.15) و ([ci], 0.942-1.020 p=0.31) rs3025034 وrs3025032 ([ci], 0.954-3.633 p=0.06 وrs182288251 در بین دو گروه بود . بنابراین یافته¬های ما از این فرضیه که ژنوتیپ rs3025033 rs3025035 ،rs3025032 ،rs3025034 ،rs114074122و rs182288251با افزایش سقط همراه است حمایت نمی¬کند. اگر چه مطالعات بیشتر در مورد ژن vegf و عملکرد بیولوژی آن در جهت درک بهتر نقش این ژن در سقط مورد نیاز است .

در این پروژه به بررسی رطوبت نسبی که در آن رطوبت، نانو ذرّات نمکی موجود در جو از جمله nacl، kcl،cscl نم پذیر می شوند وهمین طور وابستگی کشش سطحی در سیستم های نانوذرّات - بخار آب به اندازه با استفاده از شبیه سازی پرداخته شده است. شبیه سازی های دینامیک مولکولی در انسامبل کانونیکال و با استفاده از ترموستات نوز- هوور و با زمان آسایشps 1.0 انجام شده است. سیستم ها بعد ازns 5.0 به تعادل رسیدند و زمان کل هر شبیه سازی ns3 به طول انجامید. برهم کنش های یون-یون در شبیه سازی ها با استفاده از تابع پتانسیل لنارد-جونز و مولکول-های آب با استفاده از پتانسیل spc/e شرح داده شده است. شبیه سازی های دینامیک مولکولی در دمای 300 کلوین و فشار کل 1 اتمسفر و رطوبت نسبی های 70%-98% بررسی شدند. شبکه ی نانوذرّات با اندازه های مختلف، سدیم کلراید: 3.2، 1.3، 9.3، 7.4و5.5 نانومتر، سزیم کلراید: 9.2، 4.3، 6.4، 7.5و9.6 نانومتر و پتاسیم کلراید: 6.2، 5.3، 4.4، 2.5و1.6 نانومتر ساخته شدند. از آنجا که هیچ بر هم کنشی بین آب و سطح یون های سدیم کلراید کاملاً مسطح مشاهده نشده است بنابراین در ساختار اولیه ی نانوذرّات 10% حفره در نظر گرفته شده است. برای انتقال مبدأ نانوذرات به مرکز جعبه از یک تبدیل خطی استفاده می شود. این حوزه ای به قطر 3 برابر طول لبه ی نانوذرات است و سپس در مبدأ قرار داده و حجم موجود در داخل این حوزه در مجاورت نانوذرات تعریف می شود. سپس، ما تعداد مولکول های آب موجود در داخل این حوزه به عنوان تابعی از زمان در طول چند شبیه سازی در رطوبت های مختلف رسم شده است. برای محاسبه ی کشش سطحی از روش ناحیه ی مورد تست استفاده شده است. محدوده ی بالا و پایین drh بر اساس مقادیر محدود رطوبت نسبی شبیه سازی شده برای آنچه که جذب آب مشاهده شده است تعیین شد. مقادیر رطوبت نسبی نم پذیری نانوذرات nacl، kclوcscl افزایش می یابد، زمانیکه اندازه ی ذرات به ترتیب 5.5 تا 3.2، 1.6تا7.2و 8.6تا8.2 کاهش می یابند

شبیه سازی دینامیک مولکولی روشی مناسب برای مدلسازی میکروسکوپی در مقیاس اتمی و مولکولی فراهم می کند. در این پایان نامه با استفاده از شبیه سازی دینامیک مولکولی، کشش سطحی نانوذرات نمکی nh4cl محاسبه شده است. برای محاسبه کشش سطحی از روش جدیدی استفاده شده است. در این روش، با استفاده از انرژی جذب سیستم nh4cl در تماس با آب در قطرهای مختلف نانوذره، می توان drh را در rhهای مختلف به دست آورد. در واقع نم پذیری زمانی اتفاق می افتد که به یکباره با افزایش رطوبت نسبی نمودار انرژی جذب بر حسب قطر نانوذره افزایش یابد. هنگامی که اندازه نانوذره کاهش می یابد، بخش بزرگتری از اتمها روی سطح خوشه قرار دارند. از آنجاکه اتمهای سطح در مقایسه با اتمهای توده دارای انرژی پیوند کمتری هستند، همچنین با کاهش تعداد ذرات انرژی پیوستگی کاهش می یابد و سبب ایجاد ساختار ناپایدارتری می شود. همچنین با کاهش قطر نانوذره، drh افزایش می یابد. بنابراین با استفاده از وابستگی drh و همچنین کشش سطحی به قطر ذره می توان نمودارهای کشش سطحی را رسم کرد.

چکیده: گاز طبیعی یکی از اعضای خانواده سوخت های فسیلی است که استفاده از آن ها برای تامین انرژی اخیرا افزایش یافته است. ارزش گاز طبیعی با کاهش مقداردی اکسید کربن افزایش می یابد. بناابراین جدا سازی کربن دی اکسید از متان یک عملیات رایج درصنعت گاز طبیعی است. از زمان کشف نانولوله های کربنی تعداد زیادی ااز مطالعات آزمایشگاهی به جذب گاز روی نانو لوله کربنی اختصاص یافت. در این مطالعه ما از یک نانولوله کربنی با 416 اتم کربن در داخل یک جعبه شبیه سازی استفاده کردیم. اخیرا شبیه سازی های مولکولی به ابزار قدرتمندی برای مطالعه ی رفتار های جذب سطحی تبدیل شده است. شبیه سازی دینامیک مولکولی برای پیش بینی جذب هم دمای خالص ، و همچنین مخلوط در کسر های مولی مختلف از متان بر روی نانو لوله ی کربنی در دمای 300 کلوین انجام گرفت . ما شبیه سازی دینامیک مولکولی را با استفاده از نرم افزار dl-poly4.03 انجام دادیم. ساختار اولیه توسط توزیع تصادفی مخلوط دی اکسید کربن و متان روی سطح نانولوله ی کربنی ساخته شد. سیستم در طول 1 نانوثانیه توسط انسامبل nvt با ترموستات برندسن ثابت شده بود . معادلات حرکت با استفاده از الگوریتم ورلت جهش قورباغه ای با time step=1fsساخته شدند. شبیه سازی ها بعد از 1 نانوثانیه به تعادل می رسند و زمان مورد نیاز برای محاسبه خواص 2 نانو ثانیه است. همه ی بر هم کنش های بین اتمی در جعبه شبیه سازی با 10 =cutoffمحاسبه شد در طول تمام شبیه سازی ها نانو لوله فیکس بود . از نتایج گاز مخلوط می توان به عنوان یک تکنولوژی مناسب برای جداسازی کربن دی اکسید از متان استفاده کرد . به منظور بررسی تاثیر فشار گاز بر میزان جذب تمام شبیه سازی ها برای تعداد مختلفی از مولکول های گاز انجام شد . نتایج نشان دادند میزان جذب دی اکسید کربن خالص ،متان خالص و مخلوط دی اکسید کربن و متان با افزایش فشار گاز افزایش می یابد و فاکتور گزینش پذیری با افزایش فشار گاز کاهش می یابد.

در این مطالعه: ما تاثیر برهمکنش های سه ذره ای جذب سطحی گازها روی گرافیت را بررسی کردیم. اصلاح پتانسیل برهمکنش های بین مولکولی با استفاده از وارد کردن سهم مثبت برهمکنش های سه ذره ای در پتانسیل دو ذره ای انجام شد. یک ضریب تنظیم پذیر(?) برای پتانسیل دو ذره ای در نظر گرفته شد، این پارامتر نشان می دهد که بر همکنشهای سه ذره ای چه درصدی از پتانسیل دو ذره ای را شامل می شوند. مقدار سهم برهمکنش های سه ذره ای با پتانسیل دو ذره ای جمع شده تا پتانسیل اصلاح شده حاصل گردد.

به منظور بدست آوردن تابع توزیع شعاعی یک سیال لنارد-جونز محصور درون نانولوله کربنی، شبیه سازی روی 160 تا 320 اتم آرگون محصور در نانولوله کربنی با کایرالیته (33،0) از دمای اتاق (298 درجه کلوین) تا دمای 600 درجه کلوین انجام شد. پس از بدست آوردن مقادیر تابع توزیع شعاعی، آن ها را با معادلات ریاضی فیت کردیم که مقادیر بدست آمده از این معادلات توافق خوبی با مقادیر شبیه سازی داشتند. همچنین با استفاده از این تابع، انرژی پیکربندی، فشار محوری و فشار شعاعی سیال درون نانولوله را محاسبه نمودیم که مقادیر حاصل از این محاسبات نیز توافق خوبی با مقادیر شبیه سازی شده نشان داده اند. در ادامه نیز با استفاده از شبیه سازی دینامیک مولکولی و درنظرگرفتن بردارهای فشار به صورت ناهمسانگرد، دو معادله حالت متفاوت در دو جهت شعاعی و محوری برای سیال لنارد-جونز محصور شده درون نانولوله های (11،11)،(12،10)،(19،0)،(19،19)،(18،0)،(33،0)و دماهای 120 تا 600 درجه کلوین بدست آمده است. همچنین پژوهش هایی در رابطه با اثر شعاع نانولوله، کایرالیته نانولوله، چگالی ذرات درون نانولوله و دما بر روی برخی از ویژگی های ساختاری و دینامیکی سیال لنارد-جونز محصور درون نانولوله، مانند rdf و ضریب خودنفوذی انجام گرفته است.