Penelitian Dinamika Molekuler tentang Pemotongan Berlian Ultra Presisi Cerium

Abstrak

Pengkaitan antara transformasi fase struktural dan dinamika dislokasi menyulitkan pemahaman terhadap deformasi skala nano cerium. Di sini kami menggunakan dinamika molekuler (MD) untuk mengurai mekanisme cerium selama pemotongan berlian ultra presisi. Model MD yang mengandung potensi Ce–Ce dan Ce–C yang diperoleh secara empiris dibangun, dan respons elastis dua fase fcc cerium dievaluasi. Simulasi menunjukkan bahwa geser dislokasi mendominasi deformasi plastis, sementara analisis fungsi distribusi radial (RDF) mengungkap perubahan fase γ→δ yang minor pada permukaan yang dimesin dan chip. Studi parameter mengidentifikasi kondisi pemesinan optimal untuk mencapai permukaan berkualitas tinggi.

Latar Belakang

Cerium (Ce, Z=58) adalah logam lanthanida yang dihargai karena properti mekanis, fisik, dan kimiawi unik. Morfologi permukaan secara kritikal mempengaruhi performa, daya tahan korosi, dan perilaku penyimpanan hidrogen bagian cerium [1–3]. Pemotongan berlian ultra presisi dapat menghasilkan permukaan ultra halus dengan kerusakan subsurface minimal [6,7], namun sedikit pekerjaan eksperimental atau teoretis yang membahas pemotongan cerium. Karena radius ujung alat dalam pemotongan ultra presisi dapat dibandingkan dengan kedalaman pemotongan, properti material mendominasi proses, membuat mekanisme pemotongan cerium sangat menantang.

Ce menunjukkan diagram fase tekanan-suhu yang kaya yang didorong oleh delokalisasi elektron 4f. Pada tekanan ambien dan <110 K, fase α-Ce (fcc) stabil, bertransformasi menjadi β-Ce (dhcp) antara 45–275 K, γ-Ce (fcc) antara 270–999 K, dan δ-Ce (bcc) antara 999–1071 K [8–11]. Secara menarik, transformasi γ→α pada 295 K dan 8 kbar menurunkan volume sekitar 20% [12–14], sementara transisi γ→δ dapat dipicu oleh suhu dan tekanan tinggi. Transformasi ini mengubah ikatan dan properti elastis, mempengaruhi deformasi. Ce logam juga duktil, diatur oleh aktivitas dislokasi [15]. Namun, interaksi antara perubahan fase dan dislokasi selama pemotongan berlian tetap sebagian besar belum terjelajahi.

Simulasi MD menawarkan jalur yang kuat untuk meneliti mekanisme pemesinan fundamental pada skala atom. Studi sebelumnya telah memeriksa nanogores dan pemotongan ortogonal Cu, Al, Fe, dan alloy entropi tinggi, mengungkap wawasan tentang nukleasi dislokasi, tumpukan permukaan, dan perilaku gesekan [16–27]. Pekerjaan tersebut menyoroti pentingnya orientasi kristal, sudut rake alat, dan potensi interaksi. Namun, tidak ada studi MD yang telah menyelidiki pemesinan mekanik cerium, yang mendorong pekerjaan saat ini.

Metode

Model MD Pemotongan Berlian

Gambar 1 menggambarkan pengaturan simulasi: workpiece Ce kristal tunggal (41 × 25 × 31 nm) yang berisi sekitar 1 × 10^6 atom γ-Ce, dan alat pemotongan berlian (0.1 × 10^6 atom C) dengan sudut relief 9°. Workpiece dibagi menjadi lapisan bawah tetap 2 nm dan lapisan mobile yang mengikuti dinamika Newton dengan langkah waktu 1 fs. Batas periodik diterapkan hanya sepanjang arah longitudinal. Tiga orientasi kristalografis—(010), (110), dan (111)—dipertimbangkan. Tujuh sudut rake (−30°, −20°, −10°, 0°, 10°, 20°, 30°) diuji. Alat berlian diperlaku sebagai bodi kaku; sehingga interaksi C–C dihilangkan.

Interaksi atom dimodelkan sebagai berikut: Gaya Ce–Ce menggunakan parameterisasi embedded atom method (EAM) Sheng et al. [29], yang secara akurat mereproduksi konstanta elastis fcc Ce. Interaksi Ce–C menggunakan potensi Morse:

$$E_{\text{tot}}=\sum_{ij}D_0\left[e^{-2\alpha(r-r_0)}-2e^{-\alpha(r-r_0)}\right]$$ dengan D_0=0.087 eV, α=5.14, dan r_0=2.93 Å; batas potong 1.0 nm diterapkan [30].

Sistem pertama kali diseimbangkan pada 30 K, 0 bar dalam ensemble NPT, lalu ditetapkan pada pemotongan berlian dengan kecepatan konstan 100 m/s dan kedalaman pemotongan 4 nm dalam ensemble NVT. Gaya pemotongan dipecahkan sepanjang arah pemotongan, dengan arah pemotongan ditunjukkan oleh panah merah dalam gambar. Meskipun dimensi yang disimulasikan dan kecepatan pemotongan berorde magnitudo lebih kecil dari kondisi pemotongan ultra presisi eksperimental, model MD menangkap fisika esensial proses. Analisis tetangga umum (CNA) membedakan atom fcc, hcp, bcc, dan cacat [31]. Simulasi dilakukan dengan LAMMPS (langkah 1 fs) dan divisualisasikan menggunakan OVITO [32,33].

Karakterisasi Fase Cerium

Lima fase cerium—γ, α, β, ε, dan δ—dimodelkan menggunakan parameter jaringan dari literatur [8–11]. Konfigurasi bulk mengalami tarikakse uniaksial, geser, dan kompresi untuk mengekstrak konstanta elastis dan properti mekanis (Tabel 2). Modulus Young yang dihitung γ-Ce (24.17 GPa) sejalan dengan data nanoindentation eksperimental (36.7 GPa) [10], dan konstanta elastis sejalan dengan hasil neutron inelastik [34]. α-Ce yang lebih padat menunjukkan kekakuan yang jauh lebih tinggi dari γ-Ce.

Analisis fungsi distribusi radial (RDF) mengkarakterisasi perbedaan struktural antar fase. Gambar 2a menunjukkan puncak RDF yang berbeda untuk setiap fase, yang sesuai dengan jarak tetangga terdekat (misalnya, γ-Ce: 3.64, 5.13, 6.30 Å; α-Ce: 3.41, 4.85, 5.92 Å). Kompresi seragam γ-Ce untuk mencapai kerusutan volume 20% menghasilkan RDF yang cocok dengan α-Ce (Gambar 2b), mengonfirmasi transisi γ→α [12–14].

Hasil dan Pembahasan

Mekanisme Pemesinan Cerium

Simulasi MD pemotongan berlian pada Ce(010) dengan sudut rake 0° mengungkap tiga regime gaya (Gambar 3). Awalnya, gaya adheksi menghasilkan gaya pemotongan dan normal negatif; deformasi elastis mengikutinya, yang mengakibatkan kenaikan cepat gaya. Pada panjang pemotongan sekitar 2.3 nm, penurunan tajam menunjukkan awal plastisitas, yang didominasi oleh dislokasi parsial Shockley 1/6<112> yang terbentuk di dekat permukaan bebas dan meluncur pada bidang {111}. Gaya pemotongan dan normal stabil antara 10–35 nm, sementara dislokasi menumpuk di depan dan di bawah alat. Di luar 35 nm, alat terpisah dari workpiece, menyebabkan penurunan gaya. Sepanjang waktu, gaya normal lebih rendah dari gaya pemotongan.

Gambar 4 menggambarkan evolusi struktur cacat dan morfologi permukaan pada panjang pemotongan kunci. Parsial Shockley menghasilkan tumpukan bersimpang yang dibatasi oleh inti dislokasi. Seiring berjalannya pemotongan, awan padat parsial muncul dari permukaan depan, membentuk chip di sisi rake (Gambar 4f). Dislokasi di belakang alat naik ke permukaan atas, mempromosikan pemulihan permukaan. Saat alat mendekati batas kiri, propagasi dislokasi terhalang, meningkatkan volume chip (Gambar 4g). Setelah pemisahan chip, densitas dislokasi menurun karena anihilasi di permukaan bebas (Gambar 4h).

Analisis RDF zona cacat mengonfirmasi transformasi fase γ→δ yang terbatas. RDF wilayah cacat di bawah permukaan cocok dengan puncak δ-Ce, menunjukkan transisi γ→δ yang kecil yang didorong oleh pemanasan lokal (Gambar 5a). Tanda tangan δ-Ce yang sama muncul dalam chip (Gambar 5b), tetapi jumlah secara keseluruhan tetap remeh, mengindikasikan bahwa perubahan fase tidak dominan selama pemotongan.

Pengaruh Sudut Rake

Gambar 6 menunjukkan bahwa baik gaya pemotongan maupun normal menurun seiring meningkatnya sudut rake, sejalan dengan teori bidang geser Merchant [37]. Sudut rake yang lebih tinggi juga mengurangi selisih gaya dan meningkatkan stabilitas pemesinan.

Analisis cacat (Gambar 7a,b) mengungkapkan bahwa sudut rake negatif (−30°) menghasilkan densitas dislokasi yang lebih tinggi daripada sudut positif (30°). Jumlah segmen dislokasi (Gambar 7c) menurun monoton dengan sudut rake, mencerminkan deformasi plastis yang kurang kompleks dan work hardening yang lebih rendah.

Morfologi permukaan (Gambar 8) membenarkan temuan ini: sudut rake negatif menghasilkan tumpukan permukaan yang asimetris dan terucuk, sedangkan sudut rake 30° menghasilkan volume tumpukan terendah dan profil simmetris. Akibatnya, sudut rake 30° optimal untuk mencapai gaya pemotongan rendah, densitas dislokasi minimal, dan finishing permukaan superior.

Pengaruh Orientasi Kristal

Di bawah sudut rake optimal 30°, kami membandingkan orientasi Ce(010), Ce(110), dan Ce(111) (Gambar 9). Ce(010) menunjukkan gaya pemotongan dan jumlah dislokasi terendah. Meskipun Ce(111) menawarkan sedikit dislokasi, bidang {111} yang sejajar dengan permukaan mempromosikan tumpukan permukaan, meningkatkan resistensi pemotongan. Morfologi permukaan yang dimesin (Gambar 10) menunjukkan bahwa Ce(010) menghasilkan volume tumpukan terkecil dan profil paling simmetris, sedangkan Ce(110) menderita tumpukan terbesar.

Kesimpulan

Kami telah mengembangkan kerangka MD yang menggabungkan potensi EAM dan Morse untuk memodelkan pemotongan berlian ultra presisi cerium. Konstanta elastis, properti mekanis, dan perilaku transisi fase fase Ce divalidasi, mengonfirmasi kemampuan prediktif model. Simulasi mengungkapkan bahwa nukleasi dan geser dislokasi mendominasi plastisitas, dengan hanya perubahan fase γ→δ yang minor yang diamati. Kondisi pemesinan optimal—sudut rake 30° dan orientasi kristal (010)—menghasilkan gaya pemotongan terendah, densitas dislokasi minimal, dan kualitas permukaan tertinggi.

Referensi

[1]–[38] (lihat artikel asli untuk daftar lengkap kutipan)

Gambar

Gambar 1. Model MD pemotongan berlian Ce. (a) Tampilan depan, (b) tampilan atas. Merah = atom Ce bawah, biru = atom Ce mobile, abu-abu = atom C.

Gambar 2. (a) RDF fase cerium. (b) RDF sebelum dan setelah kompresi γ-Ce hingga kerusutan volume 20%.

Gambar 3. Variasi gaya pemotongan dan normal dengan panjang pemotongan. Subgambar menampilkan konfigurasi cacat per zona.

Gambar 4. Struktur cacat dan morfologi permukaan yang dimesin pada berbagai panjang pemotongan.

Gambar 5. Analisis RDF zona cacat (a) dan chip (b) setelah pemotongan.

Gambar 6. Ketergantungan sudut rake pada gaya pemotongan dan normal rata-rata.

Gambar 7. Struktur cacat untuk sudut rake −30° (a) dan 30° (b); jumlah dislokasi (c).

Gambar 8. Morfologi permukaan yang dimesin untuk sudut rake −10°, −20°, −30°, 10°, 20°, 30°.

Gambar 9. Gaya pemotongan dan jumlah dislokasi versus orientasi kristal.

Gambar 10. Morfologi permukaan yang dimesin untuk Ce(010), Ce(110), Ce(111).