Stabilitas Bentuk Nanopiring Logam: Wawasan dari Simulasi Dinamika Molekuler

Abstrak

Nanopiring logam menunjukkan kemampuan fungsional yang luar biasa yang erat kaitannya dengan morfologinya. Menggunakan simulasi dinamika molekuler atomistik, kami memeriksa stabilitas bentuk nanopiring kubus tubuh-tengah (bcc) dengan permukaan indeks rendah. Simulasi kami mengungkapkan bahwa (1) permukaan (110) secara intrinsik paling stabil, sementara nanopiring (111) dan (001) mengalami deformasi yang dipicu suhu. Papan (001) berubah menjadi struktur bertopang sediuluh yang terdiri dari faceted (110) sebelum runtuh menjadi partikel tidak teratur. Evolusi ini dikendalikan oleh distribusi tegangan geser, yang bergantung pada orientasi faceted. Dengan memodulasi orientasi faceted, ukuran papan (diameter dan ketebalan), dan komposisi material, kami menunjukkan stabilitas bentuk yang dapat disesuaikan di seluruh beberapa logam bcc (Fe, W, Nb, Mo, Cr). Temuan ini menyediakan kerangka kerja atomistik untuk mengontrol morfologi nanopiring selama sintesis dan aplikasi.

Pendahuluan

Nanopiring logam adalah nanomaterial fungsional kunci dengan aplikasi yang bervariasi dari katalisis hingga optoelektronika. Performa mereka ditentukan oleh luas permukaan, kimia faceted, dan bentuk secara keseluruhan. Meskipun termodinamika memprediksi kecenderungan menuju peminimuman energi permukaan berbentuk bola, hambatan kinetik sering mengstabilkan platelet anisotropi. Papan tipis (beberapa lapisan atom) memiliki rasio permukaan-volume yang tinggi, sehingga sangat sensitif terhadap suhu dan rangsangan lingkungan. Studi in-situ TEM sebelumnya telah mengamati fragmentasi tepi dan reorientasi faceted saat dipanaskan, namun wawasan atomistik kuantitatif masih terbatas. Di sini kami menggunakan simulasi dinamika molekuler (MD) skala besar untuk menangkap evolusi morfologi yang bergantung pada suhu secara lengkap dari nanopiring logam bcc, dengan fokus pada peran orientasi kristal, ukuran, dan komposisi.

Metode

Nanopiring Fe dibangun dengan konstanta jaringan a = 2,8665 Å dan diameter d = 32a, masing-masing terdiri dari tiga lapisan atom. Permukaan awal diatur ke orientasi (111), (001), atau (110). Struktur analog dihasilkan untuk W, Nb, Mo, dan Cr. Model dibangun dalam LAMMPS dan di-relax melalui minimisasi conjugate-gradient. Pemanasan selanjutnya dilakukan dalam ensembel NVT, menaikkan suhu dari 1 K hingga 300 K (atau lebih) dalam kenaikan 1 K, dengan timestep 2 fs dan seimbang 200 ps per langkah. Energi potensial dan tensor tegangan direkam; 8 ps terakhir dari setiap langkah dirata-rata. Gaya antaratom menggunakan potensi EAM yang divalidasi dalam studi MD sebelumnya. Tegangan lokal untuk atom i dihitung menggunakan ekspresi virial standar:

σ_{αβ}=\frac{1}{2Ω_i}\sum_{j\ne i}F_{ij}^{α}R_{ij}^{β}

di mana Ω_i adalah volume Voronoi atom i.

Hasil dan Pembahasan

Hierarki energi dan evolusi bentuk. Energi potensial awal per atom adalah −2,833, −3,457, dan −3,668 eV/atom untuk papan Fe (111), (001), dan (110), masing-masing, mencerminkan energi permukaan 2,58, 2,47, dan 2,37 J/m². Seiring meningkatnya suhu, papan (111) dan (001) dengan cepat mengadopsi bentuk lengkung yang metastabil, sedangkan papan (110) mempertahankan geometri datarnya hingga suhu simulasi tertinggi. Papan (001) mengembangkan bentuk bertopang sediuluh yang didominasi oleh faceted (110), seperti yang dikonfirmasi dengan memasang lapisan menengahnya ke permukaan kuadrik z = ax² + by² + cxy + d. Transisi kritis terjadi pada 8, 129, dan 205 K, dengan perubahan yang sesuai dalam koefisien pemasangan dan nilai R² > 0,8 yang menunjukkan geometri bertopang sediuluh yang kuat hingga 205 K.

Deformasi yang dipicu tegangan. Peta tegangan geser mengungkap distribusi anisotropi yang terkait dengan orientasi faceted: tegangan positif pada faceted merah dan negatif pada faceted biru menciptakan gaya yang berlawanan yang menggerakkan lenturan. Pada 9 K, lenturan yang meningkat menghilangkan kesenjangan tegangan, menurunkan energi potensial. Transisi sekunder antara 129–134 K sesuai dengan lenturan pada daerah pusat, yang menyebabkan penurunan tiba-tiba R² dan runtuhan akhirnya menjadi partikel tidak teratur.

Modulasi faceted dan ketergantungan material. Papan Fe "(110) yang dimodulasi", yang terdiri dari empat faceted dengan orientasi berbeda, menunjukkan pembentukan bertopang sediuluh yang serupa pada 1 K tetapi bertransisi kembali ke cakram pada 179 K dan akhirnya runtuh pada 277 K. Sebaliknya, papan (110) reguler tetap stabil hingga 552 K, menekankan pengaruh desain faceted. Memperluas studi ke papan W, Nb, Mo, dan Cr dengan d = 32a menunjukkan bahwa stabilitas bertopang sediuluh berkisar dari 62 K (Cr) hingga 582 K (W), berkaitan dengan energi potensial awal dan stabilitas struktural setiap elemen.

Efek ukuran dan ketebalan. Papan Fe (001) dengan diameter 12a, 40a, dan 50a menunjukkan jendela stabilitas yang berbeda: papan 12a runtuh pada 84 K, papan 40a mempertahankan bertopang sediuluh hingga 190 K, sementara papan 50a hanya hingga 134 K. Variasi ketebalan menunjukkan bahwa papan 4 lapis mempertahankan bentuknya hingga 98 K sebelum berubah menjadi faceted (110); papan yang lebih tipis runtuh lebih awal. Secara keseluruhan, papan yang lebih besar meningkatkan stabilitas dasar, namun kontribusi kinetik dan entropi dapat mengurangi stabilitas bertopang sediuluh pada diameter yang sangat besar.

Kesimpulan

Studi MD kami menguraikan mekanisme atomistik yang mengatur evolusi bentuk dalam nanopiring logam bcc. Permukaan (110) menawarkan stabilitas struktural tertinggi, sedangkan papan (111) dan (001) berubah menjadi bentuk bertopang sediuluh yang terdiri dari faceted (110) sebelum runtuh. Tegangan geser, yang terkait dengan orientasi faceted, menggerakkan deformasi ini. Dengan merancang orientasi faceted, dimensi papan, dan memilih komposisi material, stabilitas morfologi yang diinginkan dapat disesuaikan secara halus. Wawasan ini menuntaskan dasar teoretis untuk kontrol morfologi dalam sintesis nanomaterial logam.

Ketersediaan Data dan Bahan

Semua data simulasi dan skrip analisis termasuk dalam publikasi ini.

Singkatan

a

Konstanta jaringan

bcc

Kubus tubuh-tengah

d

Diameter

EAM

Metode atom terbenam

E_p

Energi potensial

MD

Dinamika molekuler

R²

Koefisien determinasi

TEM

Mikroskopi transmisi elektron