Evolusi Morfologi yang Digali Difusi Permukaan Substrat Si(001) Berpola Lubang: Studi Lapangan Fase

Latar Belakang

Heteroepitaksi logam tidak cocok dari semikonduktor seperti Ge/Si atau InGaAs/GaAs sering menghasilkan pembentukan pulau tiga dimensi melalui mode pertumbuhan Stranski–Krastanow (SK). Meskipun pengorganisasian sendiri dapat menghasilkan titik kuantum, nukleasi acak yang inherent pada proses ini sering menghasilkan distribusi ukuran dan bentuk yang luas yang membatasi kinerja perangkat.

Selama beberapa dekade, berbagai teknik telah dikembangkan untuk mengarahkan pertumbuhan heteroepitaksial menuju struktur terurut. Di antaranya, substrat berpola lubang telah muncul sebagai salah satu pendekatan paling serbaguna untuk mencapai pengurutan samping yang tinggi dan kontrol ukuran yang tepat pulau.

Pola lubang dibuat menggunakan metode top-down seperti nanoimprint litografi, litografi sinar elektron digabungkan dengan etsi reaktif-ion atau etsi kimia basah, dan nanoindentation. Teknik ini memungkinkan desain array lubang terurut dengan presisi nanometer, yang, dalam kondisi pertumbuhan yang dioptimalkan, menghasilkan pengurutan samping pulau yang hampir sempurna.

Karena bentuk lubang mempengaruhi total energi permukaan dan, akibatnya, nukleasi pulau, kontrol akurat morfologi lubang sangat penting. Pada suhu tinggi, kapilaritas mengarahkan evolusi lubang menuju konfigurasi yang energis menguntungkan, kadang-kadang menghasilkan penyembuhan lengkap. Oleh karena itu, aniling atau depoisi material tambahan secara rutin digunakan untuk memstabilkan bentuk lubang metastabil. Bahkan setelah stabilisasi, evolusi lebih lanjut dapat terjadi selama fase pertumbuhan heteroepitaksial.

Dalam studi ini, kami menggambarkan evolusi substrat berpola lubang yang digali oleh pengurangan energi permukaan melalui difusi permukaan. Kami mengadopsi formulasi lapangan fase yang menangkap skala panjang dan waktu yang relevan dengan eksperimen dan telah sebelumnya diterapkan pada kinetika terbatas difusi dalam sistem heteroepitaksial. Model ini mampu mereproduksi bentuk seimbang dengan energi permukaan anizotropik yang realistis dan telah divalidasi terhadap data eksperimental.

Kami berfokus pada permukaan Si(001) berpola lubang, sistem yang telah banyak diselidiki dalam literatur. Kertas ini diatur sebagai berikut: model lapangan fase diperkenalkan pada bagian berikutnya; perataan lubang Si(001) dan jalur kinetik menuju keseimbangan dibahas pada Bagian 3; aplikasi pada pertumbuhan Ge, yang menginduksi rotasi lubang, disajikan pada Bagian 4; dan kertas ini diakhiri dengan diskusi tentang implikasi temuan kami.

Metode

Model Lapangan Fase

Model lapangan fase memperkenalkan parameter urutan kontinu φ yang bervariasi halus dari φ = 1 (padat) ke φ = 0 (vakuum). Evolusi φ diatur oleh persamaan Cahn–Hilliard degeneratif yang menangkap difusi permukaan sambil menghormati anizotropi kuat energi permukaan. Fungsi energi bebas total adalah

F = ∫Ω γ(n̂) \,[ ε/2 |∇φ|² + 1/ε B(φ) ] d𝑟 + ∫Ω β/(2ε) [ -ε∇²φ + 1/ε B′(φ) ]² d𝑟,

di mana Ω adalah domain simulasi, ε ketebalan antarmuka, B(φ)=18φ²(1−φ)² potensi double-well, dan β parameter regulasi Willmore yang membulatkan sudut tajam dalam regime anizotropi kuat. Potensial kimia μ=δF/δφ mengarahkan evolusi melalui

∂φ/∂t = D ∇·[ M(φ) ∇μ ],

dengan D koefisien difusi dan M(φ)=36/ε φ²(1−φ)² mobilitas permukaan. Ekspresi untuk μ mencakup densitas energi permukaan anizotropik γ(n̂) dan istilah regulasi Willmore, memastikan deskripsi yang realistis secara fisik dari evolusi facies.

Energi Permukaan Anizotropik

Kami memodelkan densitas energi permukaan sebagai

γ(n̂) = γ₀ [ 1 − Σᵢ αᵢ (n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ) ],

di mana vektor mᵢ sesuai dengan arah kristalografi energi minimum (〈001〉, 〈113〉, 〈110〉, 〈111〉 untuk silikon). Kedalaman setiap minimum αᵢ berasal dari nilai energi permukaan relatif γᵢ yang dilaporkan dalam literatur, sedangkan parameter lebar wᵢ mengontrol jangkauan sudut facies. Untuk silikon kami menetapkan γ₀ = 1, α_{〈001〉}=0.15, dan nilai w 50 (kecuali w_{〈113〉}=100). Parameter regulasi Willmore β dipilih sebagai 0.005 untuk menyelesaikan sudut yang dibulatkan sambil mempertahankan kelayakan komputasi.

Morfologi Awal dan Pengaturan Simulasi

Kami membangun geometri lubang awal pada permukaan (001) yang didefinisikan oleh depresi lingkaran dengan radius L dan kedalaman H, yang terhubung halus ke bidang di sekitarnya menggunakan ekor Gaussian dengan lebar σ=L/2. Jarak bertanda d(r) dari antarmuka dimasukkan ke dalam profil tangen hiperbolik φ(r)=½[1−tanh(3 d(r)/ε)] untuk menghasilkan kondisi awal. Ukuran domain tipikal L = 1 (dalam satuan arbitrer) dan ketebalan antarmuka ε = 0.2 digunakan. Kotak elemen hingga AMDiS digunakan dengan skema integrasi waktu semi-implicit dan penyempurnaan jaring di antarmuka. Kondisi batas periodik diterapkan di arah samping, dan kondisi tanpa fluks digunakan di sepanjang sumbu vertikal. Koefisien difusi D ditetapkan ke satuan, jadi waktu simulasi dilaporkan dalam satuan arbitrer.

Hasil dan Pembahasan

Perataan Lubang Si(001)

Gambar 2 menunjukkan evolusi awal waktu lubang dengan rasio aspek R=0.25 dan R=0.5. Profil yang awalnya tanpa fitur secara spontan mengembangkan facies yang sesuai dengan minimum densitas energi permukaan silikon. Untuk lubang dangkal (R=0.25), dasar persegi 〈001〉 dikelilingi oleh tepi 〈113〉 dan segitiga 〈110〉 yang sempit. Untuk lubang yang lebih dalam (R=0.5), facies 〈111〉 tambahan muncul di dekat dasar, dan facies 〈110〉 menjadi lebih lebar. Prediksi ini setuju secara kuantitatif dengan pengamatan eksperimental permukaan Si(001) berbentuk lubang.

Gambar 3 melacak evolusi jangka waktu panjang lubang yang lebih dalam (R=0.5). Urutan tampilan perspektif dan atas menunjukkan hilangnya bertahap facies 〈111〉 yang curam, kasarnya facies 〈113〉 yang bersebelahan, dan hilangnya facies 〈110〉, yang bersama-sama menghasilkan permukaan yang rata. Total energi permukaan menurun monotonik, seperti yang ditunjukkan pada Gambar 4, dan kinetika mengikuti jalum umum untuk semua R≤0.5. Skala waktu karakteristik untuk lubang puluhan nanometer pada 1100–1200 °C adalah pada urutan jam, sesuai dengan nilai literatur untuk difusi permukaan silikon.

Meniru Perubahan Bentuk karena Pertumbuhan Ge

Untuk mengevaluasi dampak depoisi materi dengan himpunan facies seimbang yang berbeda, kami mensimulasikan pertumbuhan Ge pada substrat Si(001) berpola lubang dengan R≈0.1. Memulai dari morfologi pada t≈5.0 (Gambar 3), kami mengganti densitas energi permukaan silikon dengan satu yang mencakup minimum dalam arah 〈105〉, karakteristik sistem Ge/Si(001). Parameter lebar w_{〈105〉} ditingkatkan ke 500 untuk menangkap facies 〈105〉 yang dangkal.

Gambar 5 menunjukkan bahwa permukaan secepatnya terorganisir ulang: facies 〈105〉 memulai nukleasi di antara facies 〈113〉 yang tersisa, memperluas pada bekas yang terakhir. Lubang akhir dibatasi secara eksklusif oleh facies 〈105〉, dan tampilan atas menunjukkan rotasi 45° garis lubang relatif terhadap orientasi awal. Rotasi ini telah dilaporkan secara eksperimental selama depoisi Ge dan kini ditunjukkan berasal sepenuhnya dari pengurangan energi permukaan.

Kesimpulan

Kami telah mengembangkan kerangka lapangan fase kontinu yang menangkap evolusi yang digali difusi permukaan substrat Si(001) berpola lubang. Dengan secara eksplisit memasukkan anizotropi kuat energi permukaan silikon, model ini mereproduksi pembentukan facies metastabil, perataan lubang yang akhirnya terjadi, dan jalur kinetik yang mengarah ke keadaan ini. Memperluas deskripsi energi permukaan untuk mencakup facies 〈105〉 menunjukkan bahwa rotasi garis lubang yang diamati selama pertumbuhan Ge dapat dijelaskan sepenuhnya oleh pengurangan energi permukaan. Model ini umum dan dapat dengan mudah disesuaikan dengan substrat lain dengan re-parameterisasi anizotropi, menawarkan alat prediktif untuk merancang geometri lubang yang memandu penempatan pulau heteroepitaksial dengan presisi nanometer.