صفحه اصلی | ان آر تی سی | صفحه اصلی<

09392522438  
   EN | FA
دوشنبه, 14 مهر 1399 ساعت 09:57

Hydrogenated Ψ-graphene as an ultraviolet optomechanical sensor

PSI (ψ)-graphene is a dynamically and thermally stable two-dimensional (2D) allotrope of carbon composed of 5-6-7 carbon rings. Herein, we study the opto/mechanical behavior of two graphene allotropes, Ψ-graphene and its hydrogenated form, Ψ-graphane under uniaxial and biaxial strain using density functional theory (DFT) calculations. We calculated the elastic constants and second Piola-Kirchhoff (PK2) stresses, in which both nanostructures indicate a similar elasticity behavior to graphene. Also, the plasmonic behavior of these structures in response to various strains has been studied. As a result, plasmonic peaks varied up to about 2 eV under strain. Our findings reveal that these two structures have a large peak in the ultraviolet (UV) region and can be tuned by different applied strain. In addition, Ψ-graphene has smaller peaks in the IR and UV regions. Therefore, both Ψ-graphene and Ψ-graphane can be used as UV optomechanical sensors, whereas Ψ-graphene could be used as an infrared (IR) and visible sensor.

 

CIF files:

- Ψ-graphene

- Ψ-graphane

 

Ref: https://pubs.rsc.org/en/content/articlepdf/2020/ra/d0ra03104f

 

 

 

 

منتشرشده در مقالات ما

 

محاسبات ساختاری انجام شده بر روی پروسکایت ها

پروسکایتها، با فرمول شیمیایی ABX3 (شکل-1)، بطور معمول سلول واحدی شامل 5 اتم در یک ساختار مکعبی دارند (فاز α)، که در آنها کاتیون B، 6 همسایه ی اول و آنیون X و کاتیون A هر کدام 12 همسایه اول دارند [1].

 

Perovskite unit cell

شکل-1: سلول واحد پروسکایت ABX3 [2].

در وضعیت ایده آل، این ساختار برای آنکه حداکثر تقارن مکعبی را حفظ کند، شعاع یون های A، B و X باید طوری باشند که ضریب تلورانس[1] t تقریبا برابر یک شود:

(‏2‑1)

t=(RA+RX)/[(√2)(RA+RX)]

که در این رابطه RA، RB و RX به ترتیب شعاع یون های متناظر ساختار ABX3 هستند. انحراف از مقدار 1 باعث می شود ساختار مکعبی واپیچیده[2] شده و بلور تقارن کمتری خواهد داشت. باید توجه داشت، برای آنکه t≈1 باشد یون A می بایست بسیار بزرگتر از یون B باشد. در پروسکایت های هالوژنی، جایگاه B معمولا با اتم های بزرگ Pb یا Sn اشغال شده و بنابر این اتم A می بایست بسیار بزرگتر باشد. Cs تقریبا بزرگترین اتم گروه I جدول تناوبی عناصر می باشد. اما با این حال Cs به اندازه ی کافی برای حفظ پایداری ساختار مکعبی پروسکایت بزرگ نیست، بنابراین می بایست توسط یک مولکول بزرگ جایگزین شود. شاید دلیل آنکه ترکیب CH3NH3PbI3 پایدارتر بوده و عملکرد بهتری نسبت به CsPbI3 دارد همین باشد.

در دماهای محدود اگر1<t<8/0 باشد ساختار مکعبی همچنان وجود خواهد داشت [3]کمترین مقدار t منجر به ساختار چهارگوش[3] با کمترین تقارن (فاز β) و یا لوزی رخ[4] (فاز γ) می گردد، در حالی که مقدار بزرگ t (t>1) پایداری شبکه ی سه بعدی B-X را برهم زده و منجر به یک ساختار دوبعدی و لایه ای می شود. لازم به یادآوری است که در محاسبات DFT در دمای صفر، فاز γ همیشه پایدارترین فاز بوده و فاز α ناپایدار می باشد، چرا که نمی تواند شرط t=1 را برآورده کند. با این وجود، در دمای محیط اغلب انتقال از یک فاز به فاز دیگر اتفاق می افتد [4, 5]. در حقیقت، مدهای فونونی آکوستیک عرضی به راحتی آنیون X را از وسط فاصله ی اتمهای B-B فاز مکعبی جابجا می کنند. به خاطر شعاع یونی متفاوت و استحکام متفاوت ساختارها، دمای گذار پروسکایت های مختلف با یکدیگر فرق می کنند. برای پروسکایت CH3NH3PbI3 گذار فاز α به β به γ، به ترتیب در دماهای K 330 و K 160 رخ میدهد [6]. برای برخی پروسکایت ها مانند HC(NH2)2PbI، FAPbI3، CsPbI3 و CsSnI3 یک فاز δ غیرپروسکایتی نیز یافت شده است [4, 5]. بر خلاف فاز β و γ، فاز δ نمی تواند از فاز α بر اثر واپیچش زاویه ی B-X-B پدید بیاید، بلکه در عوض با شکسته شدن پیوند B-X بوجود می آید.

Perovskite training CH3NH3PbI3

شکل-2: ساختار پروسکایت CH3NH3PbI3، الف) فاز α، ب) فاز β، ج) فاز γ، د) فاز δ ا[7].

 

ساختار اتمی پروسکایت های آلی-معدنی پیچیده تر از پروسکایت های معدنی هستند. مولکول آلی، با تقارن غیر مرکزی، می تواند آرایش و جهتگیری های گوناگونی را به خود بگیرد (شکل-2). برای ترکیب CH3NH3PbI3 در فاز α مولکلول CH3NH3 در دماهای بالا آرایش تصادفی داشته و در نهایت بلور تقارن Oh به خود می گیرد. در دماهای معمولی (فاز β) آزادی مولکول CH3NH3کاهش یافته ولی هنوز جهتگیری فضایی متنوعی دارد. در دماهای پایین (فاز γ) مولکول CH3NH3 در جای خود ثابت بوده و جهتگیری های منظمی خواهند داشت [5].

بی نظمی جهتگیری مولکولها در فاز α و β را نمی توان در یک سلول واحد محدود مدل سازی کرد. برای فاز α با مقایسه ی انرژی های ساختار CH3NH3PbI3 در شرایطی که پیوند C-N در راستاهای [001]، [110] و [111] جهتگیری کرده اند، کمترین انرژی و پایدارترین حالت قرار گیری پیوند C-N در راستای [111] بدست آمده است. اما مطالعات محدودی نیز به پایدارتر بودن راستای [001] اشاره کرده اند، البته تفاوت چندانی وجود نداشته است. لازم به ذکر است که ساختار نواری و الکترونی CH3NH3PbI3 به جهتگیری CH3NH3 بستگی چشمگیری ندارد. نتایج محاسباتی و تجربی پارامتر شبکه ی فازهای α، β، γ و δ در جدول-1 گزارش شده است [1].

Calculated grid parameter and band gap

جدول-1: پارامتر شبکه و گاف نواری محاسبه شده برای سه فاز CH3NH3PbX3 (X=I, Br, Cl) و مقایسه آنها با نتایج تجربی. چون مولکول CH3NH3 تقارن مرکزی ندارد، برای فاز α پارامتر شبکه ی محاسبه شده برای شبه مکعب بصورت a=V1/3 و همچنین برای فاز β، a=1/2(a+b) تعریف شده است [8].

گزارش هایی که منتشر شده نشان می دهند که تقریب LDA پارامتر شبکه را کمتر از مقدار واقعی و تقریب GGA بیشتر از مقدار واقعی محاسبه می کنند. با این حال محققین، برای ترکیب CH3NH3PbI3، برآورد بیشتر پارامتر شبکه های محاسبه شده به روش GGA-(PBE)، را به برهمکنش های واندروالسی نادیده گرفته شده مابین CH3NH3 و شبکه ی [Pb-I] نسبت داده اند [9]. زمانی که برهمکنش دافعه ای وارد محاسبات می شود پارامترهای شبکه محاسبه شده با نتایج تجربی توافق خوبی خواهند داشت [9-11]. بر همین اساس برهمکنش های واندروالس در پروسکایت های هیبرید هالیدی به نظر بسیار مهم می آیند. البته باید اشاره کرد که اهمیت برهمکنش واندروالس به دو دلیل مهم زیر است:

اولاً، تقریب GGA/PBE نه فقط برای CH3NH3PbI3 بلکه برای سایر ترکیبات نیز منجر به بیشتر برآورد کردن پارامتر شبکه ها می شود و این اتفاق ناشی از مشکل ذاتی این تابع تبادلی-همبستگی است و نه نوع ماده. ثانیاً، اگر برای پروسکایت معدنی CsPbI3 برهمکنش دافعه ای را در نظر بگیریم، پارامتر شبکه ی محاسبه شده با نتایج تجربی یکی خواهد شد. یعنی استفاده از تقریب GGA است که منجر به کاهش طول پیوند ترکیبات می شود و با اضافه کردن برهمکنش دافعه ای به محاسبات، طول پیوندها افزایش یافته و خطای تقریب GGA اصلاح می شود. این نتایج نشان می دهند که CH3NH3 به مانند یک یون رفتار کرده و برهکنش CH3NH3 و [Pb-I] مشابه برهمکنش Cs و [Pb-I] بوده و ویژگی های یونی از خود نشان می دهد. تاکنون برای پروسکایتهای ترکیبی، نتایج بهینه سازی ساختاری با تقریب PBEsol

بهترین همخوانی را با نتایج تجربی حاصل از پراش نوترون و اشعه ی X داشته است (جدول-2).

 

 

Experimental and theoretical values of lattice parameter and bond length

جدول-2: مقادیر تجربی و نظری پارامتر شبکه و طول پیوند های CH3NH3PbI3DFT[11], X-ray[4], X-ray3[5], Neutron[12], X-ray5[13].

References

[1]        Yin W-J, Yang J-H, Kang J, Yan Y, & Wei S-H (2015) Halide perovskite materials for solar cells: a theoretical review. Journal of Materials Chemistry A 3(17):8926-8942.

[2]        Yin W-J, Shi T, & Yan Y (2014) Unusual defect physics in CH3NH3PbI3 perovskite solar cell absorber. Applied Physics Letters 104(6):063903.

[3]        Li C, Lu X, Ding W, Feng L, Gao Y, & Guo Z (2008) Formability of ABX3 (X= F, Cl, Br, I) Halide Perovskites. Acta Crystallographica Section B: Structural Science 64(6):702-707.

[4]        Baikie T, Fang Y, Kadro JM, Schreyer M, Wei F, Mhaisalkar SG, Graetzel M, & White TJ (2013) Synthesis and crystal chemistry of the hybrid perovskite (CH 3 NH 3) PbI 3 for solid-state sensitised solar cell applications. Journal of Materials Chemistry A 1(18):5628-5641.

[5]        Stoumpos CC, Malliakas CD, & Kanatzidis MG (2013) Semiconducting tin and lead iodide perovskites with organic cations: phase transitions, high mobilities, and near-infrared photoluminescent properties. Inorganic chemistry 52(15):9019-9038.

[6]        Ball JM, Lee MM, Hey A, & Snaith HJ (2013) Low-temperature processed meso-superstructured to thin-film perovskite solar cells. Energy & Environmental Science 6(6):1739-1743.

[7]        Yin WJ, Shi T, & Yan Y (2014) Unique properties of halide perovskites as possible origins of the superior solar cell performance. Advanced Materials 26(27):4653-4658.

[8]        Yin W-J, Yan Y, & Wei S-H (2014) Anomalous alloy properties in mixed halide perovskites. The journal of physical chemistry letters 5(21):3625-3631.

[9]        Egger DA & Kronik L (2014) Role of dispersive interactions in determining structural properties of organic–inorganic halide perovskites: insights from first-principles calculations. The journal of physical chemistry letters 5(15):2728-2733.

[10]      Wang Y, Gould T, Dobson JF, Zhang H, Yang H, Yao X, & Zhao H (2013) Density functional theory analysis of structural and electronic properties of orthorhombic perovskite CH3NH3PbI3. Physical Chemistry Chemical Physics 16(4):1424-1429.

[11]      Menéndez-Proupin E, Palacios P, Wahnón P, & Conesa JC (2014) Self-consistent relativistic band structure of the CH 3 NH 3 PbI 3 perovskite. Physical Review B 90(4):045207.

[12]      Eames C, Frost JM, Barnes PRF, O’regan BC, Walsh A, & Islam MS (2015) Ionic transport in hybrid lead iodide perovskite solar cells. Nature communications 6:7497.

[13]      Weller MT, Weber OJ, Henry PF, Di Pumpo AM, & Hansen TC (2015) Complete structure and cation orientation in the perovskite photovoltaic methylammonium lead iodide between 100 and 352 K. Chemical Communications 51(20):4180-4183.


[1]Tolerance factor

[2]Distorted

[3]Tetragonal

[4]Orthorhombic

منتشرشده در مقاله
یکشنبه, 19 آبان 1398 ساعت 01:03

نانونوار یا ریبون های گرافن

نانونوار یا ریبون های گرافن را بشناسیم

 

دانشمندان معاصر، مهندسین و سرمایه گذاران پیش بینی می کنند که علم نانو مسیری است به سوی میدانی وسیع اما در ابعادی کوچک. تلاش محققین در هرچه کوچکتر کردن ابعاد دستگاه های میکروالکترونیک در سال های اخیر بسیار شدت یافته و در ابعاد بسیار پایین پدیده ای به نام محدودیت کوانتومی[1] الکترون ها در نانو ساختار ها، می تواند ابزار قدرتمندی برای کنترل خواص الکتریکی، نوری، مغناطیسی، ترموالکتریکی این مواد پیشرفته باشد.

 

انواع نانوساختارهایی که تا کنون سنتز شده اند عبارتند از:

نانولوله[2] ها

نانوریبون ها

نانو میله[3] ها

نانوسیم[4] ها

و نقاط کوانتومی.

در این میان، نانوریبون ها یکی از مهمترین اعضای خانواده ی نانومواد یک بعدی است که ساختارش بطور گسترده مورد مطالعه ی محققین قرار گرفته و انتظار می رود که کاندید خوبی برای حسگر، مبدل و تشدید کننده نانو مقیاس باشد. جذابیت های پژوهشی نانوریبون ها ناشی از مورفولوژی، خواص فیزیکی، الکترونیکی و نوری منحصر به فرد این ساختارها است.

 

در میان تمامی ساختارهای کربنی، یکی از انواع ساده ی ساختارهای مشابه گرافن، ریبون های گرافنی با عرضی به ابعاد چند نانومتر هستند که بسیار مورد مطالعه قرار گرفته اند [1]. اهمیت این ساختار در کاربردهای آن ها در قطعات الکترونیکی و ترانزیستور ها می باشد که به تازگی روش های متعددی برای استفاده از این ساختارها در این قطعات ابداع شده است. نانوریبون های گرافنی نیز با روش های متعددی مانند برش های مکانیکی گرافن های ورقه ورقه شده و یا با روش های برآرایی[5] با استفاده از ماسک های ویژه ای تولید می شوند.

از زمانی که ریبون های گرافنی مورد توجه قرار گرفته اند، مطالعات بسیاری بر روی ساختار الکترونی این مواد و تاثیر تغییر پهنا و عرض این ماده بر روی خواص مختلف آن با روش های مختلفی مانند تقریب بستگی قوی و تئوری تابعی چگالی صورت گرفته است. این محاسبات نشان دادند که [6]GNR ها با لبه های آرمچیر، هم می توانند فلز باشند و هم نیمه رسانا، و گذار میان این دو حالت با توجه به تغییر عرض GNR رخ می دهد. و GNR هایی با لبه های زیگزاگ همگی فلز هستند [2].

در سال 2006، یانگ و همکارانش نشان دادند که GNR هایی که لبه های آنها توسط هیدروژن اشباع شده اند، همگی گاف نواری مستقیم و غیر صفری دارند [3]. به علاوه آنها نشان دادند که اندازه ی گاف نواری تابعی از پهنای GNR ها می باشد.

این محاسبات به کمک رهیافت تئوری تابعی چگالی و تقریب چگالی (اسپینی) جایگزیده[7] انجام شده است.

 

 

Geometric shape of graphene nanoribbons

شکل-1: شکل هندسی نانوریبون های گرافنی: a) آرمچیر، b) زیگزاگ.

 

References

[1] M. Ezawa, Physical Review B, 73 (2006) 045432.

[2] M. Fujita, K. Wakabayashi, K. Nakada, K. Kusakabe, Journal of the Physical Society of Japan, 65 (1996) 1920-1923.

[3] Y.-W. Son, M.L. Cohen, S.G. Louie, Physical review letters, 97 (2006) 216803.

 


[1]Quantum Confinement

[2]Nanotube

[3]Nanorod

[4]Nanowire

[5]Epitaxy

[6]Graphene Nanoribbons

[7] Local (Spin) Density Approximation )LSDA(

منتشرشده در مقاله
یکشنبه, 19 آبان 1398 ساعت 00:21

آلایش در گرافن (doping)

آلایش در گرافن و تاثیر آن بر خواص این ماده

گرافن داپت شده با نیتروژن

با توجه به تنوع بسیار زیاد در فرهای مختلف ساختار گرافن و خواص منحصر به فرد این ماده، به تازگی مورد توجه محققین جهت استفاده در ابرخازن های قرار گرفته است. نیتروژن، در جدول تناوبی عناصر در کنار کربن قرار داشته و الکترون نگاتیوی آن از کربن بیشتر است، به همین دلیل جایگزین کردن یکی از اتم های کربن در گرافن منجر به تغییر در ساختار الکترونی آن می گردد. بر همین اساس، گرافن داپت شده با نیتروژن بر اساس ساختار ویژه اش کاندید امیدبخشی جهت ارتقاء عملکرد ابر خازن ها، باتری های لیتیومی، و ترانزیستورها می باشد [1].

ترکیب نیتروژن با گرافن می تواند سه نوع ساختار متفاوت را شکل دهد:

گرافیتی-[1]N

پیریدینیک-N[2]

و پیرولیک-[3]N(شکل-1).

در نوع گرافیتی-N اتم نیتروژن داپت شده جایگزین یکی از اتم های حلقه ی هگزاگونال گرافن می شود.

اما در ساختارهای پیریدینیک-N و پیرولیک-N، الکترون های اوربیتال π پیوندهایی با هیبریداسیون، به ترتیب، sp2 و sp3 تشکیل می دهند [2].

Different structures induced by nitrogen atom dopants in graphene

شکل-1: ساختارهای مختلف ناشی از داپت اتم نیتروژن در گرافن.

 

بر اثر داپینگ نیتروژن، تراز فرمی به بالا ی نقطه ی دیراک جابجا شده و سپس چگالی حالت ها در اطراف تراز فرمی به یک حالت ثابتی منتهی می شود که در نتیجه میان نوار رسانش و نوار ظرفیت گاف نواری ای تقریبا به اندازه ی 0.2 eV پدید می آید (شکل-2) [3].

به همین دلیل گرافن داپت شده با نیتروژن به عنوان کاندید مناسبی جهت استفاده در ادوات نیمه هادی بسیار مناسب است.

 

Band gap created by doping of nitrogen to graphene at point K

شکل-2: گاف نواری ایجاد شده بر اثر داپ نیتروژن به گرافن، در نقطه ی K.

 

گرافن داپت شده با برم

عنصر برم (B)، خواص منحصر بفرد و غیر قابل مقایسه ای دارد. به این دلیل که برم در جدول تناوبی عناصر در کنار کربن قرار دارد، محققین بسیاری امکان داپ برم را در گرافن بررسی کرده اند [4]. ساختار گرافن داپت شده با برم، به روش جالبی سنتز می شود.

در این روش اکسیدگرافیت در حضور اکسید برم (به عنوان منبع برم) بازپخت شده و گرافن داپت شده با برم تولید می گردد. به علاوه الکترون نگاتیوی برم از کربن کمتر است و از همین روی B-گرافن از خود پتانسیل خوبی جهت استفاده در واکنش کاهش اکسیژن (ORR[4]) نشان می دهد.

همچنین به دلیل پایداری و طول عمر بالا، به عنوان ماده ی اصلی کاتالیزور های مبتنی بر پلاتین (Pt) مورد استفاده قرار می گیرد. علاوه بر این، مطالعات مبتنی بر رهیافت DFT نشان داده اند که داپت برم در گرافن منجر به القاء چگالی اسپینی بالاییدر ساختار گرافن می گردد که نقش بسیار تاثیرگذاری بر جذب اکسیژن و مولکول OOH در مقایسه با گرافن خالص دارد [5].

همانطوری که در شکل-3 نشان داده شده است، از آنجایی که اتم برم تنها سه الکترون در لایه ی ظرفیت دارد، اوربیتال pz آن خالی خواهد بود، در مقابل تک الکترون درون اوربیتال نیمه پر pz اتم کربن همسایه در آن یک چگالی اسپینی جایگزیده القاء می کند. از آنجایی که الکترون نگاتیوی اتم برم از اتم کربن کمتر است، پیوند تشکیل شده بیشتر به سمت کربن متمایل می شود و در نتیجه این جایگاه قطبیده ی پدید آمده در ساختار گرافن داپت شده با برم، به عنوان سایت فعال کاتالیز می تواند عمل کند. 

 

How to fill carbon and bromine pz orbitals in B-graphene

شکل-3: نحوه ی پر شدن اوربیتال های pz کربن و برم در B-گرافن.

 

گرافن داپت شده با گوگرد

نخستین بار ساختار الکترونی گرافن داپت شده با گوگرد بصورت تئوری محاسبه شد، که نشان داد این ترکیب با توجه به درصد سولفور داپت شده در گرافن، می تواند نیمه رسانایی با گاف نواری کوچک و یا یک فلز باشد [6]. دنیس[5] و همکارانش در این تحقیق نشان دادند که حتی یک تک اتم گوگرد در یک صفحه ی گرافنی 5×5 می تواند تاثیر جدی در ساختار نواری این ماده پدید بیارورد. محاسبات آنان که با نرم افزار محاسباتی SIESTA انجام شد نشان داد که گوگرد در این صفحه ی گرافنی گاف نواری ای به اندازه ی 0.3 eV و یک پیک بسیار تیز در نزدیکی تراز فرمی به وجود می آورد (شکل-4).

 

The density of graphene states a before sulfur doping b after sulfur dubbing

شکل-4: چگالی حالت های گرافن، a) قبل از داپ گوگرد، b) بعد از داب گوگرد.

 

برای سنتز S-گرافن، از هر دو ماده ی سولفید هیدروژن و دی سولفید بنزیل به عنوان مواد تامین کننده ی گوگرد، به همراه اکسید گرافن استفاده می شود.

 

گرافن داپت شده با سیلیسیم

تحقیقات تئوری و تجربی بسیاری جهت اصلاح عملکرد سنسوری گرافن جهت شناسایی مولکولهای مختلف انجام شده است. مکانیزم این سنسورهای گازی بر مبنای تغییر رسانندگی الکتریکی گرافن به واسطه ی انتقال بار میان گرافن و ماده ی جذب شونده می باشد. اکثر مطالعات نشان می دهند که گازهای NO، NO2 و ... که به عنوان گازهای آلاینده می شناسیم، بصورت فیزیکی بر روی گرافن خالص جذب می شوند. داپ سیلیسیم (Si) در گرافن بسیار کمیاب است، اما به تازگی محققین جهت استفاده در کاتالیزورهای مبتنی بر گرافن به سراغ این ماده نیز رفته اند. تا کنون گزارشی مبنی بر استفاده از Si-گرافن در در کاتالیست ها ارائه نشده است.

زاهو و همکارانش [7] به کمک رهیافت تئوری تابعی چگالی (DFT) به بررسی جذب گازهای NO، N2O و NO2 بر روی Si-گرافن پرداختند و نشان دادند که این ماده می تواند سنسور بسیار خوبی برای شناسایی گازهای NO و NO2 باشد. داپ اتم Si در گرافن باعث جابجایی نوار رسانش گرافن شده و آن را کمی به بالاتر منتقل می کند در نتیجه گاف نواری بسیار کوچکی به اندازه ی 0.054 eV در این ماده به وجود می آید (شکل-5). با وجود آنکه Si چهار ظرفیتی است اما در این ساختار به دلیل طول پیوند های بلند تری که با اتم کربن تشکیل می دهد (Si-C) ساختار هرم شکلی پیدا می کند و به همین دلیل هیبریداسیون sp3 به خود می گیرد، که همین هیبریداسیون منجر به ایجاد ترازهای جایگزیده در اطراف تراز فرمی شده و تاثیر کوچکی بر ویژگی های شبه فلزی گرافن می گذارد.

 

Geometric shape and band structure of graphene doped with silicon

شکل-5: شکل هندسی و ساختار نواری گرافن داپت شده با سیلیسیم.

 

 

References

[1] Y. Lu, Y. Huang, M. Zhang, Y. Chen, Journal of nanoscience and nanotechnology, 14 (2014) 1134-1144.

[2] X.-K. Kong, C.-L. Chen, Q.-W. Chen, Chemical Society Reviews, 43 (2014) 2841-2857.

[3] D. Usachov, O. Vilkov, A. Gruneis, D. Haberer, A. Fedorov, V.K. Adamchuk, A.B. Preobrajenski, P. Dudin, A. Barinov, M. Oehzelt, Nano letters, 11 (2011) 5401-5407.

[4] A. Lherbier, X. Blase, Y.-M. Niquet, F. Triozon, S. Roche, Physical Review Letters, 101 (2008) 036808.

[5] X. Kong, Q. Chen, Z. Sun, ChemPhysChem, 14 (2013) 514-519.

[6] P.A. Denis, R. Faccio, A.W. Mombru, ChemPhysChem, 10 (2009) 715-722.

[7] Y. Chen, B. Gao, J.-X. Zhao, Q.-H. Cai, H.-G. Fu, Journal of molecular modeling, 18 (2012) 2043-2054.

 


[1]Graphitic N

[2]Pyridinic N

[3]Pyrrolic N

[4] Oxygen Reduction Reaction

[5]Denis

منتشرشده در مقاله
دوشنبه, 06 آبان 1398 ساعت 15:17

نانو محاسبات

نانومحاسبات چیست؟
تئوری تابعی چگالی (Density Functional Theory - DFT) یکی از کارآمدترین روش های محاسباتی سیستم های بس ذره ای است و این به دلیل کارایی محاسباتی بالای آن، کاهش حجم محاسبات و دقت خوب آن در مقایسه با سایر روش ها است. اهمیتDFT هنگامی روشن تر می شود که بدانیم به خاطر آن جایزه نوبل شیمی در 1998 به کوهن یکی از بنیان گذاران آن داده شده است. کوهن و شم یک دستگاه فرضی شامل الکترون هایی که هیچ برهم کنشی برهم ندارند را در نظر گرفتند و در نهایت به معادلاتی رسیدند که با آن ها به صورت خود-سازگار خواص حالت پایه دستگاه را محاسبه کردند. پیدا کردن این انرژی به صورت کلی امکان ندارد بنابراین باید از روشهای تقریبی استفاده کرد که بعضی از این تقریب ها، تقریب چگالی موضعی (LDA) و تقریب شیب تعمیم یافته (GGA) است.
همانگونه که در شکل (1) دیده می شود DFT انقلابی در عرصه ی مطالعه ی تئوری جامدات و درشت مولکول ها پدید آورده و در سالهای اخیر رشد چشمگیر و فزاینده ای در بررسی خواص مواد و طراحی مواد جدید داشته است. تلاش های بسیاری جهت افزایش دقت محاسباتی مبتنی بر DFT انجام شده و به مدد توسعه ی سخت افزارهای کامپیوتری و افزایش توان پردازشی کامپیوترهای محاسباتی در دو دهه ی گذشته، امکان مطالعه ی سیستم ها و ساختار های بزرگتر و پیچیده تر را فراهم کرده است.

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

 

منتشرشده در مقاله
پنج شنبه, 02 آبان 1398 ساعت 13:57

DFT الکترونی

 

بررسی خواص الکترونی ساختارها به کمک DFT

میزان رسانایی الکترونی در مواد جامد بسیار متنوع است. بر اساس میزان مقاومت مواد در عبور جریان الکتریکی، مواد مختلف را می‌توان به دسته‌های مختلف رسانا، نیمه رسانا و عایق دسته‌بندی کرد.

ساختار نواری یعنی ترازهایی از انرژی که الکترون می‌تواند در آنها حضور داشته باشد و با حضور نداشته باشد.

به مناطقی که الکترون نمی‌تواند حضور داشته باشد منطقه ی ممنوعه، گاف نواری یا Band Gap می گویند. با محاسبه ی ساختار نواری و گاف نواری مواد می توان به فلز (Metal)، شبه فلز (Semimetal)، نیمه رسانا (Semiconductor) و یا نارسانا (Insulator) بودن مواد پی برد. حتی با محاسبات اسپینی می توان به ساختار هایی که نیمه فلز (Half-metal) هستند نیز پی برد.

امروزه تئوری تابعی چگالی (DFT) به یکی از رایجترین روش های بررسی خواص الکتریکی نانوساختار ها و ساختار های بلوری تبدیل شده است و توسعه ی روش ها و تقریب های محاسباتی منجر به افزایش دقت این محاسبه ها شده است. برخی از اطلاعاتی را که می توان از این رهیافت استخراج کرد در زیر مشاهده می کنید:

    • رسم ساختار نواری (BandStructure)
    • محاسبه ی گاف نواری (Band gap)
    • رسم چگالی حالت ها (DOS)
    • سم چگالی حالت های موضعی (PDOS)
    • محاسبه ی جرم موئثر حامل ها
    • محاسبه ی تحرک پذیری حامل ها
    • تاثیر ناخالصی بر خواص الکتریکی
    • تاثیر نقص بلوری بر خواص الکتریکی.

IMAGE-xqy

 

 

 

در حال حاضر پروژه های بسیاری در قالب پایانامه های دکتری و ارشد برای داشنجویان تحصیلات تکمیلی رشته های فیزیک، شیمی، برق و مهندسی مواد جهت بررسی خواص الکتریکی مواد بلوری و نانوساختار های تعریف شده و به یکی از زمینه های بسیار محبوب در میان اساتید و دانشجویان این رشته در آمده که منجر به هجوم محققین به این علم نوین گردیده است.

به دلیل وسعت بسیار زیاد این حوزه، انتشار مقاله در ژورنال های بسیار معتبر به سادگی امکانپذیر است. نرم افزارهای مختلفی می توانند این پارامترها را محاسبه کنند که از پرکاربردترین آنها Quantum ESPRESSO و Siesta می باشند.

 

 

نمونه مقالاتی که در این حوزه منتشر شده اند:

 

    • Silva, Juliana CM, Heitor A. De Abreu, and Hélio A. Duarte. "Electronic and structural properties of bulk arsenopyrite and its cleavage surfaces–a DFT study." RSC Advances 5.3 (2015): 2013-2023.
    • Park, Changwon, et al. "Electronic Properties of Bilayer Graphene Strongly Coupled to Interlayer Stacking and an External Electric Field." Physical review letters 115.1 (2015): 015502.
    • Rossi, Antonio, et al. "Nano-Scale Corrugations in Graphene: A Density Functional Theory Study of Structure, Electronic Properties and Hydrogenation." The Journal of Physical Chemistry C 119.14 (2015): 7900-7910
منتشرشده در مقاله
دوشنبه, 29 مهر 1398 ساعت 19:50

LAMMPS دانشگاه خواجه نصیر


دوره مقدماتی نرم افزار محاسباتی لمپس - LAMMPS
دانشگاه صنعتی خواجه نصیرالدین طوسی
تاریخ: 30 اردیبهشت 1395

 

LAMMPS

 

 

LAMMPS 2

 

 

LAMMPS 3

 

 

LAMMPS 4

 

 

LAMMPS 5

 

 

LAMMPS 7

 

 

LAMMPS 8

 

 

LAMMPS 9

 

 

منتشرشده در پیش رو

 Nano plat elongation

 

 

بررسی خواص مکانیکی ساختارها (نانولوله- نانوریبون- صفحات- ساختارهای بالک)

 

تولید فلزات و آلیاژهایی با اندازه کم‌تر از 100 نانومتر باعث دستیابی به موادی با استحکام بسیار بالا شده است. در واقع کوچک کردن دانه‌ها در مواد، ابزار قدرتمندی است تا ساختارهایی با خواص مکانیکی عالی تولید گردد. در بررسی خواص مکانیکی مواد نانوساختار مشکلات زیادی از جمله عدم امکان تهیه نمونه مطلوب، وجود تخلخل و میکروترک، تنش‌های داخلی شدید، وجود ناخالصی‌ها و گازهای حبس شده و نیز عدم امکان ارزیابی برخی کمیت‌ها، نظیر اندازه‌گیری کرنش به دلیل کوچک بودن نمونه‌ها وجود دارد. وجود چنین مشکلاتی باعث شده تا داده‌های آزمایشگاهی مربوط به خواص مکانیکی برای این گروه از مواد محدود باشد. از همین روی محاسبه ی دقیق تئوری این خواص با روش DFT در بررسی خواص اینگونه مواد می تواند بسیار حائز اهمیت باشد.

پایه ی اصلی محاسبه ی خواص مکانیکی و الاستیکی مواد تغییر انرژی و به تبع آن تغییر نیروهای بین اتمی در ساختار می باشد. بدین ترتیب با اعمال کرنش بر ماده و محاسبه ی تغییر انرژی و یا تغییر تنش در ساختار و رسم نمودارهای انرژی-کرنش و یا تنش-کرنش می توان بسیاری از خواص الاستیکی و مکانیکی ساختارها (نانولوله- نانوریبون- صفحات- ساختارهای بالک) را محاسبه کرد.

 

 

.برخی از مهمترین خواص مکانیکی (نانولوله- نانوریبون- صفحات- ساختارهای بالک) در زیر لیست شده است.

    • نمودار انرژی-کرنش
    • نمودار تنش کرنش
    • تمامی ضرایب الاستیک
    • مدول یانگ (Young's modulus)
    • مدول بالک (Bulk modulus)
    • مدول برشی (Shear modulus)
    • ضریب پوآسون (Poisson's ratio)
    • تراکم پذیری
    • ناحیه الاستیک-پلاستیک (Elastic-Plastic region)

 

lagrangian strain

 

 

 

لازم به ذکر است که این اطلاعات در قالب پروژه ها ارشد و دکتری دانشجویان رشته های مکانیک، فیزیک، متالوژی و شیمی قابل تعریف و پیاده سازی بوده و به لحاظ جدید بودن و به روز بودن این دانش، به راحتی امکان چاپ مقالات متعدد در این حوزه وجود دارد. برای محاسبه ی این پارامتر های می توان از بسته های محاسباتی بسیار پرقدرت و دقیقی به نام های Quantum ESPRESSO و Siesta بهره برد.

 

 

نمونه مقالاتی که در این حوزه منتشر شده اند:

 

 Miller, David C., Mauricio Terrones, and Humberto Terrones. "Mechanical properties of hypothetical graphene foams: Giant Schwarzites." Carbon 96 (2016): 1191-1199. 

 Wu, Guo Xun, et al. "Mechanical and Electrical Properties of TiN with Stacking Fault: A DFT Study." Applied Mechanics and Materials. Vol. 727. Trans Tech Publications, 2015.

 Peng, Qing, Zhongfang Chen, and Suvranu De. "A density functional theory study of the mechanical properties of graphane with van der Waals corrections." Mechanics of Advanced Materials and Structures 22.9 (2015): 717-721.

 

 

 

 

 

منتشرشده در مقاله

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

ایمیل مرکز محاسباتی NRTC:

این آدرس ایمیل توسط spambots حفاظت می شود. برای دیدن شما نیاز به جاوا اسکریپت دارید

 

 

منتشرشده در اخبار مرکز
سه شنبه, 18 دی 1397 ساعت 15:41

Quantum ESPRESSO

 

Quantum ESPRESSO :

بسته ی نرم افزاری کوانتوم اسپرسو، نرم افزاری چند منظوره برای محاسبات ابتدا به ساکن در سیستم های ماده چگال (تناوبی و بی نظم) می باشد که بر پایه ی نظریه ی تابعی چگالی (DFT) کار می کند. از این نرم افزار بصورت گسترده در بررسی و شبیه سازی نانو ساختار ها استفاده می شود.

کلمه اسپرسو به معنای قهوه فشرده است، که مخفف کلمات زیر می باشد:

opEn Source Package for Research in Electronic Structure, Simulation and Optimaization

 

ین نرم افزار برای کار در محیط لینوکس نوشته شده و دارای سه کد مجزای کامپیوتری PWscf، FPMD و CP می‌باشد که هر سه‌ی آن‌ها بر پایه‌ی نظریه‌ی تابعی چگالی کار می‌کنند و برای انجام محاسبات از امواج تخت به همراه شبه پتانسیل استفاده می‌کنند. همچنین این بسته‌ی نرم افزاری قابلیت محاسبه‌ی اکثر تقریب‌های انرژی تبادلی- همبستگی مانند LDA، GGA، LSDA و LDA+U را دارد.


قابلیت هایQE

در بسته محاسباتی Quantum ESPRESSO سه نوع محاسبات Pwscf، Phonon و PostProc انجام می پذیرد.
1- Pwscf

Pwscf می تواند محاسبات زیر را انجام دهد:

1. محاسبه نیروهای بین اتمی، تنش، و بهینه سازی ساختار یونی.
2. محاسبه دینامیک مولکولی در حالت پایه سطح بورن-اوپنهایمر(Born-Oppenheimer) ، با تغییرات سلول.
3. قطبش ماکروسکوپی و میدان الکتریکی محدود از طریق تئوری مدرن قطبش(Berry Phases) .
4. محاسبه انرژی آزاد سطح در سلول ثابت از طریق متا دینامیک. ( به کمک بسته PLUMED)
5. محاسبه‌ی انرژی حالت پایه‌ی سیستم و استخراج ویژه تابع‌ها و مقدارهای معادلات کوهن- شم با استفاده از حل خودسازگار سیستم.
6. محاسبه‌ی تنش وارد بر دیوارهای سلول واحد و نیروی وارد بر اتم‌ها.
7. محاسبه‌ی واهلش سیستم یا دینامیک مولکولی، که به دو روش قابل بررسی است. در حالت اول با ثابت نگه داشتن دیواره‌های سلول واحد و جابجا کردن اتم‌ها در فضای آن تا سیستم به حالت پایه برسد و حالت دوم حرکت دادن دیواره‌ها که ابعاد سلول واحد اجازه‌ی تغییر دارند.
8. محاسبه‌ی فرکانس فونون ها با استفاده از نظریه اختلالی تابعی چگالی.
9. محاسبه‌ی قطبیدگی ماکروسکوپیک از طریق فاز بِری.
10. بررسی برهمکنش اسپین-مدار.

 

همه محاسبات فوق برای هر دو نوع عایق و فلزات، در هر ساختار بلوری، برای بسیاری از توابع همبستگی-تبادلی (exchange correlation) از جمله قطبش اسپین، DFT + U، توابع غیرموضعیVdW، توابع هیبریدی، شبهه پتانسیلهای بار پایسته (Hamann-Schluter-Chiang PPs (NCPPs)) یا فوق نرم (Vanderbilt, PPs, USPPs) مورد استفاده قرار می گیرد. درمحاسبات، مغناطیس غیرهمراستا و برهمکنش اسپین-مدار نیز لحاظ شده است. همچنین میدان الکتریکی محدود با یک پتانسیل الکتریکی دندان اره ای در یک سلول در بسته پیاده سازی شده است (NEB).


2- Phonon

پکیج Phonon در بسته ی Quantum ESPRESSO می تواند محاسبات زیر را انجام دهد:

1. محاسبه فرکانسهای فونونی و ویژه بردارهای آنها با استفاده از نظریه تابعی چگالی اختلالی Density-Functional perturbation Theory .
2. بار موثر و تانسورهای دی الکتریک.
3. ضریب برهمکنش الکترون-فونون در فلزات.
4. برهمکنشهای ثابت نیرو در فضای مستقیم.
5. طول عمر فونونهای غیرهماهنگ مرتبه سوم.
6. سطح مقطع مادون قرمز و رامان.

 

Phonon هم مانند Pwscf میتواند از DFT+U وVdW استفاده کند. اما USPP و PAW در محاسبات فونونی پیاده سازی نشده اند. محاسبات تقریب شبه هارمونیک Quasi-Harmonic approximations با استفاده از بسته QHA انجام میشود که نیاز به محاسبات چگالی حالتهای فونونی دارد.
3- PostProc

پکیج PostProc در Quantum ESPRESSO میتواند محاسبات زیر را انجام دهد:

1. اسکن تصاویر میکروسکوپ تونل زنی (STM).
2. رسم توابع محلی سازی الکترونی (ELF).
3. چگالی حالتها (DOS) و چگالی حالتهای تصویر شده (PDOS).
4. میانگین گیری سطحی و کروی.


دوره ی آموزشی QE

دوره ی آموزشی Quantum ESPRESSO در دو سطح مقدماتی و پیشرفته برگزار می شود. مباحث آموزش داده شده در هر سطح به شرح زیر می باشد:


سطح مقدماتی:

1. اهمیت و جایگاه فیزیک/شیمی محاسباتی.
2. آشنایی با مبانی نظری Quantum ESPRESSO از جمله DFT و شبه پتانسیل.
3. نصب نرم افزار Quantum ESPRESSO.
4. آشنایی با محیط نرم افزار Quantum ESPRESSO و بخش های مختلف آن.
5. نحوه ی ورودی نویسی و معرفی ساختار مورد نظر به QE و انجام تنظیمات اولیه.
6. بهینه سازی ساختار مورد مطالعه و محاسبه طول و زوایای پیوندی و شکل کلی ساختار.
7. آشنایی با نرم افزار کمکی XcrysDen جهت بررسی شماتیک جزئیات ساختار مورد مطالعه.
8. رسم بانداستراکچر، DOS و PDOS.
9. محاسبه و رسم توپولوژی بار (Charge Density- ELF).
10. مطالعه ی جذب مولکول بر روی سطوح متنوع تحت فشارهای مختلف.


سطح پیشرفته:

1. بررسی خواص مکانیکی ساختارها (تئوری/محاسباتی) به کمک Quantum ESPRESSO.
2. بررسی خواص اپتیکی ساختارها (تئوری/محاسباتی) با استفاده از Quantum ESPRESSO.
3. بررسی خواص فونونی ساختارها و محاسبه ی توابع ترمودینامیکی و شیمیایی آنها (تئوری/محاسباتی) در QE.
4. بررسی خواص رسانندگی ساختارها (تئوری/محاسباتی) در QE.
5. اعمال میدان الکتریکی خارجی بر ساختارها (تئوری/محاسباتی) در QE.
6. بررسی خواص گوناگون ساختارها با اضافه کردن برهمکنش اسپین- مداری (تئوری/محاسباتی) در QE.
7. آشنایی با محاسبات مبتنی بر نظریه ی تابعی چگالی وابسته به زمان، TDDFT (تئوری/محاسباتی) در QE.

 

 

منتشرشده در مقاله
صفحه1 از2

7 روز هفته، 24 ساعته پاسخگوی شما هستیم

social 16social 13social 09 social 05