Manufaktur industri
Industri Internet of Things | bahan industri | Pemeliharaan dan Perbaikan Peralatan | Pemrograman industri |
home  MfgRobots >> Manufaktur industri >  >> Industrial materials >> bahan nano

Sifat Osilasi Elektromagnetik Longitudinal pada Logam dan Eksitasinya pada Permukaan Planar dan Bulat

Abstrak

Definisi umum dari permitivitas dispersi spasial direvisi. Respons gas elektron terdegenerasi pada medan listrik yang memenuhi persamaan vektor Helmholtz ditemukan dengan solusi persamaan Boltzmann. Fungsi dielektrik longitudinal yang dihitung bertepatan dengan yang diperoleh Klimontovich dan Silin pada tahun 1952 dan Lindhard pada tahun 1954. Namun, itu tergantung pada kuadrat bilangan gelombang, parameter vektor persamaan Helmholtz, tetapi bukan vektor gelombang dari gelombang elektromagnetik bidang. Konsep baru ini menyederhanakan simulasi efek nonlokal, misalnya, dengan teori Lorents-Mie yang digeneralisasi, karena tidak ada transformasi Fourier yang harus dibuat. Koefisien Fresnel digeneralisasi memungkinkan untuk eksitasi gelombang elektromagnetik longitudinal. Untuk memverifikasi teori, spektrum kepunahan untuk bola berukuran nanometer perak dan emas dihitung. Untuk partikel-partikel ini, teori Lorents-Mie yang digeneralisasikan memberikan pergeseran biru dan perluasan resonansi plasmon yang sangat sesuai dengan data eksperimen. Selain itu, teori nonlokal menjelaskan hilangnya resonansi plasmon yang diamati untuk bola emas dengan diameter kurang dari atau sama dengan 2 nm. Perhitungan menggunakan Klimontovich-Silin-Lindhard dan fungsi dielektrik hidrodinamik untuk perak ditemukan memberikan hasil yang mendekati pada energi foton dari 3 hingga 4 eV. Kami menunjukkan bahwa nilai absolut bilangan gelombang gelombang longitudinal dalam zat padat jauh lebih tinggi daripada nilai gelombang transversal.

Latar Belakang

Iradiasi permukaan logam bidang oleh pulsa laser femtosecond sering menghasilkan pembentukan struktur permukaan periodik yang diinduksi laser (LIPSSs) [1]. Selain LIPSS, riak hiperhalus yang disebut LIPSS frekuensi spasial tinggi (HSFL) diamati [1, 2]. Periode spasial HSFL secara signifikan lebih kecil dari panjang gelombang iradiasi λ 0 . Misalnya, untuk aluminium periode ini diperkirakan berkisar antara 20 hingga 200 nm pada λ 0 =0,8 μ m [2, 3]. Sementara orientasi riak dalam LIPSS biasa tegak lurus terhadap polarisasi sinar laser, orientasi HSFL seringkali tegak lurus dan terkadang sejajar dengan polarisasi. HSFL serupa terbentuk pada permukaan dielektrik transparan, semikonduktor, dan logam. Asal usul HSFL dijelaskan oleh mekanisme yang berbeda seperti generasi harmonik kedua, keterlibatan jenis mode plasmon tertentu, pengaturan sendiri, dan peningkatan medan lokal selama kerusakan homogen dalam bahan dielektrik [2, 3].

Tujuan dari penelitian ini adalah untuk mencari proses gelombang yang dapat menghasilkan suatu pola dengan periode yang pendek Λ λ 0 . Kami memeriksa sifat gelombang elektromagnetik longitudinal (L) dalam logam yang juga dikenal sebagai gelombang plasma. Studi kami terdiri dari langkah-langkah baru berikut. Pertama, kami memulai penelitian kami dengan definisi dispersi spasial permitivitas. Seperti yang ditunjukkan di bawah ini, definisi umum tidak berguna jika media yang diteliti tidak seragam dan tak terbatas. Oleh karena itu, kami mengusulkan konsep baru fungsi dielektrik dispersif spasial ε . Fungsi ini menetapkan proporsionalitas langsung antara dua bidang vektor, E (r ,ω ) dan D (r ,ω ), tetapi bukan amplitudo E (k ,ω ) dan D (k ,ω ) gelombang bidang. Akibatnya, kuantitas ε tergantung pada kuadrat dari bilangan gelombang, k 2 , parameter vektor persamaan Helmholtz untuk medan listrik E (r ,ω ), tetapi bukan vektor gelombang k dari gelombang pesawat. Kemudian, untuk menurunkan fungsi baru seperti itu, kami menentukan respons elektron konduksi pada mode elektromagnetik dengan menyelesaikan persamaan transpor Boltzmann yang ditulis dalam perkiraan waktu relaksasi. Yang disebut fungsi dielektrik Lindhard transversal dan longitudinal diperoleh. Selanjutnya, kami menemukan bahwa Lindhard longitudinal dan fungsi hidrodinamika yang lebih sederhana memiliki jarak yang dekat dalam berbagai parameter. Kepunahan cahaya oleh nanospheres perak dan emas dianggap untuk menggambarkan teori. Kami menunjukkan untuk pertama kalinya bahwa teori Mie nonlokal menjelaskan pergeseran biru, perluasan, dan akhirnya lenyapnya resonansi plasmon yang diamati dengan penurunan ukuran nanosfer logam mulia. Akhirnya, model teoretis yang baru dikembangkan diterapkan untuk memeriksa kemungkinan keterlibatan mode longitudinal dalam pembentukan struktur permukaan yang diinduksi laser. Untuk tujuan ini, kami memodifikasi teori Fresnel dengan memperhitungkan gelombang longitudinal yang ditransmisikan.

Metode

Untuk menentukan medan elektromagnetik dalam media homogen sepotong-sepotong, teori elektromagnetik klasik diterapkan. Medan listrik E di setiap domain seragam dari media heterogen diasumsikan sebagai solusi dari persamaan vektor Helmholtz (VHE):

$$ \Delta\,\mathbf{E} + k^{2}\, \mathbf{E}=0, $$ (1)

dimana Δ adalah operator Laplace.

Seperti biasa, komponen tangensial dari listrik E dan magnet H bidang terus menerus melintasi batas-batas media. Selain itu, kami memperhitungkan bahwa elektron terkurung dalam logam; oleh karena itu, kondisi batas tambahan (ABC) berikut untuk komponen normal rapat arus j di permukaan logam S digunakan:(j n )| r S =0.

Untuk menentukan arus konduksi dalam logam, kami memecahkan persamaan transportasi Boltzmann (BTE) yang ditulis dalam pendekatan waktu relaksasi:

$$ \frac{\partial f}{\partial t}+\mathbf{v}\,\frac{\partial f}{\partial \mathbf{r}}+ \frac{e}{m}\,\ kiri(\mathbf{E}+\mathbf{v}\times\mathbf{B} \kanan)\,\frac{\partial f}{\partial \mathbf{v}} =\frac{f_{0}- f}{\tau}, $$ (2)

dimana f adalah fungsi distribusi partikel tunggal dalam ruang fase (r , v ), v adalah kecepatan elektron mikroskopis, e dan m adalah muatan dan massa elektron berturut-turut, B adalah induksi magnetik, f 0 adalah fungsi distribusi kesetimbangan, dan τ adalah waktu relaksasi.

Di bawah ini, kami memperoleh rumus untuk fungsi dielektrik dispersif spasial. Kemudian, kami menggunakannya untuk mempelajari pantulan cahaya dari permukaan logam datar dan hamburan cahaya pada nanosfer logam mulia.

Hasil dan Diskusi

Dispersi Spasial ε dalam Media Heterogen

Dalam literatur, fungsi dielektrik dispersif spasial ε didefinisikan melalui relasi berikut [4–6]:

$$ \mathbf{D}(\omega,\, \mathbf{r}) =\epsilon_{0} \iiint\limits_{-\infty}^{\infty} \! \mathbf{d} \mathbf{r}^{\prime}\, \epsilon\left(\omega, \, \mathbf{r}-\mathbf{r}^{\prime}\right)\,\mathbf {E}\kiri(\omega,\, \mathbf{r}^{\prime}\kanan), $$ (3)

dimana ε 0 adalah konstanta listrik, D (ω , r ) adalah amplitudo osilasi vektor perpindahan dengan frekuensi sudut ω di titik r , dan E (ω , r ) amplitudo osilasi medan listrik di titik r . Transformasi Fourier dari Persamaan. (3) berikan persamaan

$$ \mathbf{D}(\omega,\, \mathbf{k}) =\epsilon_{0} \, \epsilon(\omega,\, \mathbf{k})\,\mathbf{E}(\ omega,\, \mathbf{k}) $$ (4)

di mana dispersi spasial ε (ω , k ) bergantung pada vektor gelombang k dari gelombang elektromagnetik bidang. Menurut pendapat kami, Persamaan. (3) tidak ambigu hanya dalam volume homogen tak terbatas tetapi kita berurusan dengan sistem heterogen sepotong-sepotong di mana batas-batas harus diperhitungkan dan k tidak sama di media yang berbeda.

Pendekatan kami tidak menggunakan perluasan gelombang elektromagnetik di atas gelombang bidang. Permitivitas dispersi spasial menentukan hubungan antara D (ω , r ) dan solusi khusus untuk vektor Helmholtz Persamaan. (1):

$$ \mathbf{D}(\omega,\, \mathbf{r}) =\epsilon_{0} \, \epsilon(\omega, \, k)\,\mathbf{E}(\omega,\, \mathbf{r}). $$ (5)

Di sini E (ω , r ) menunjukkan distribusi medan listrik tetapi bukan hanya vektor E di titik r .

Fungsi Dielektrik Membujur dan Transversal

Permitivitas logam biasanya dinyatakan melalui konduktivitas σ [4]:

$$ \epsilon=\epsilon_{\mathrm{g}}+\frac{i\,\sigma}{\omega\,\epsilon_{0}}, $$ (6)

dimana ε g adalah bagian dari fungsi dielektrik yang memungkinkan terjadinya polarisasi padatan; ε g =1 untuk logam sederhana. Untuk menentukan σ , kami menghitung kerapatan arus

$$ \mathbf{j}=e \iiint\limits_{-\infty}^{\infty} \! \mathbf{v}\, f\, \mathrm{d}\/ \mathbf{v}=\sigma\,\mathbf{E}, $$ (7)

di mana \(\mathrm {d} \mathbf {v}=\frac {v}{m}\,\mathrm {d} \epsilon \,\mathrm {d}\,\Omega,\) d ​​Ω =dosaθ dθ dϕ , v , θ , ϕ adalah koordinat bola dari kecepatan. Tidak seperti penelitian sebelumnya, kami tidak memperkenalkan vektor gelombang k tetapi menemukan solusi BTE dalam bentuk deret tak hingga yang mengandung operator v bertindak atas v E :

$$ f=f_{0} + \frac{e}{-i \omega+\Gamma}\,\frac{\partial f_{0}}{\partial \epsilon} \left[ 1+\frac{\mathbf {v}\,{\mathbf{\nabla}}}{-i \omega+\Gamma} \right]^{-\,1}\! \mathbf{v}\,\mathbf{E}, $$ (8)

dimana Γ =1/τ . Lalu, f 0 didekati dengan distribusi Fermi-Dirac suhu nol dan, setelah integrasi lebih dari ε dalam Persamaan. (7), kita mendapatkan

$$ \mathbf{j}=\frac{\omega_{\mathrm{p}}^{2}\,\varepsilon_{0}}{-i \omega+\Gamma}\,\frac{3}{4 \ pi} \iint \mathbf{u}\,\left(1+ l\,\mathbf{u} \mathbf{\nabla} \right)^{-\,1} (\mathbf{u}\,\mathbf {E})\, \mathrm{d}\,\Omega, $$ (9)

di mana \(\omega _{\mathrm {p}}^{2}=\frac {e^{2}\,n_{e}}{m\,\varepsilon _{0}},\) p adalah frekuensi plasma, \(\mathbf {u}=\frac {\mathbf {v}}{v}\) adalah vektor satuan dalam arah v ,\(l=\frac {v_{\mathrm {F}}}{-i \omega +\Gamma },\) v B adalah kecepatan Fermi. Selanjutnya, kami menghitung integral

$$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\, (\mathbf{u}\,\mathbf{E})\,\mathbf{d}\,\ Omega =\frac{4\/\pi}{3}\,\mathbf{E} \end{array} $$ (10) $$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\,(\mathbf{u}\,\mathbf{\nabla})^{2 n-1} (\mathbf{u}\,\mathbf{E})\,\mathbf{d} \,\Omega =0 \end{array} $$ (11) $$\begin{array}{*{20}l} &\textstyle \iint \mathbf{u}\,(\mathbf{u}\, \mathbf{\nabla})^{2 n} (\mathbf{u}\,\mathbf{E})\,\mathbf{d}\,\Omega =\frac{4 \pi}{2 n+3 } \\ &\times \, \Delta^{n-1} \left[ \mathbf{\nabla}\, (\mathbf{\nabla}\cdot\mathbf{E})-\frac{1}{2 n+1}\, \mathbf{\nabla}\times \mathbf{\nabla}\times\mathbf{E}\right] \end{array} $$ (12)

dimana n adalah bilangan asli. Ketergantungan berikut dari j dengan sewenang-wenang medan listrik E akhirnya diperoleh

$$\begin{array}{*{20}l} \mathbf{j}&=\frac{\omega_{\mathrm{p}}^{2}\,\varepsilon_{0}}{-i \omega+ \Gamma} \left\{\mathbf{E} + 3\, \sum\limits_{n=1}^{\infty} l^{\,2 n} \kanan. \\ &\quad\left.\times \frac{\Delta^{n-1}}{2 n+3} \left[ \mathbf{\nabla}\, (\mathbf{\nabla}\cdot\mathbf {E})-\frac{\mathbf{\nabla}\times \mathbf{\nabla}\times\mathbf{E}}{2 n+1}\right]\right\}. \end{array} $$ (13)

Ada dua jenis solusi untuk Persamaan. (1), bebas divergensi yang memenuhi persamaan ·E =0 dan tanpa rotasi yang memenuhi persamaan

$$ \mathbf{\nabla} \times \mathbf{E}=0. $$ (14)

Untuk gelombang bidang, dengan E exp[i (k r ω t )], Persamaan. (14) ditransformasikan ke dalam relasi k ×E =0 yang menunjukkan bahwa gelombangnya membujur (L). Untuk mensimulasikan proses dalam tubuh bola, akan lebih mudah untuk menggunakan harmonik bola vektor L , M , dan N sebagai satu set lengkap fungsi ortogonal. Dalam hal ini, Persamaan. (14) menentukan harmonik L . Bilangan gelombang dari gelombang L dan L mode ditentukan oleh hukum dispersi berikut

$$ \epsilon^{\mathrm{L}}\left(\omega, \, k^{\mathrm{L}}\right)=0. $$ (15)

Dari Persamaan. (6) dan (13) kami menemukan bahwa solusi untuk Persamaan. (1) memenuhi kendala Persamaan. (14) berikan permitivitas longitudinal berikut

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{2} \, \Phi \left(a^{2},\, 1,\,\frac{3}{2} \right) $$ (16)

dimana Φ adalah fungsi Phi Lerch,

$$ \frac{3}{2} \, \Phi \left(a^{2},\, 1,\,\frac{3}{2} \right) =\sum\limits_{n=0} ^{\infty} \frac{3}{2 n+3}\, a^{2\,n}, ​​$$ (17)

\(a=\frac {k v_{\mathrm {F}}}{\omega +i \Gamma }\).

Permitivitas yang diperoleh berbeda dari yang didefinisikan oleh Kliewer dan Fuchs [7] hanya dalam notasi:

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}+ \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{a^{2}}\left[1-\frac{1}{ia} \tan^{-1}(ia) \right] $$ (18)

Identitas

$$ \frac{1}{i a}\tan^{-1}(i a)=\frac{1}{2}\ln\frac{1+a}{1-a} $$ (19)

memungkinkan seseorang untuk menulis ulang Persamaan. (18) sebagai berikut

$$ \epsilon^{\mathrm{L}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega\,(\omega+i \Gamma)}\,\frac{3}{a^{2}}\left[1-\frac{1}{2 a} \, \frac{\ln(1+a)}{\ln(1 -a)} \kanan]. $$ (20)

Dalam kasus Γ =0, rumus ini mengambil bentuk persamaan yang diturunkan oleh Klimontovich dan Silin [8] yang mempelajari dumping Landau dalam plasma degenerasi (lihat [9], [10, Persamaan (40.17)], dan [11]). Permitivitas dari Persamaan setara. (16), (18), dan (20) biasa disebut fungsi dielektrik Lindhard (dengan mengacu pada [12]) meskipun fungsi ini pertama kali diperoleh oleh Klimontovich dan Silin [8].

Permitivitas Lindhard melintang [7] dapat ditemukan dengan Persamaan. (13) ketika ·E =0. Dalam kasus sebenarnya v B k ω , ini mengurangi fungsi dielektrik Drude

$$ \epsilon^{\mathrm{T}}=\epsilon_{\mathrm{g}}- \frac{\omega_{\mathrm{p}}^{2}}{\omega^{2}+i\ ,\Gamma\/\omega}. $$ (21)

Fungsi ini sesuai dengan data eksperimen pada banyak logam [13]. Jika |a |<1, permitivitas longitudinal (16) disederhanakan menjadi fungsi dielektrik hidrodinamik:

$$ \epsilon^{\mathrm{L}}\left(\omega,\,k^{\mathrm{L}}\right)=\epsilon_{\mathrm{g}}- \frac{\omega_{\ mathrm{p}}^{2}}{\omega^{2}+i\,\Gamma \omega-\frac{3}{5}\,\left(v_{\mathrm{F}}\, k ^{\mathrm{L}}\kanan)^{2}}. $$ (22)

Pemantulan Gelombang Elektromagnetik Bidang dari Permukaan Logam Datar

Kondisi Batas

Pada bagian ini, kita menentukan arah vektor gelombang k L dan amplitudo gelombang L yang tereksitasi dalam logam selama pantulan gelombang elektromagnetik bidang dari permukaan logam datar.

Pertimbangkan insiden gelombang bidang pada antarmuka dielektrik-logam z =0 dengan vektor gelombang terletak di xz pesawat. Medan listrik dalam media dielektrik 1 terdiri dari kejadian E i dan mencerminkan E r gelombang, medan pada logam 2 memiliki transversal E t dan, dalam beberapa kasus, memanjang E L komponen. Menurut kondisi batas Maxwell, komponen transversal dari vektor medan listrik dan magnet kontinu dalam bidang z =0. Selain itu, elektron tidak dikeluarkan dari logam; oleh karena itu, komponen normal rapat arus listrik adalah nol pada z =0,

$$ \hat{\mathbf{z}}\,\mathbf{j}|_{z=0}=0. $$ (23)

adalah \(\hat {\mathbf {z}}\) adalah vektor satuan dalam arah z sumbu.

Semua suku dalam kondisi batas Maxwell harus memiliki ketergantungan yang sama pada x dan y . Persyaratan ini memiliki beberapa konsekuensi. Pertama, dapat ditentukan bahwa gelombang L dapat tereksitasi hanya dalam kasus polarisasi p ketika vektor listrik dari gelombang datang E (i ) adalah sejajar dengan bidang datang. Dengan kata lain, polariton plasmon dapat dibangkitkan oleh gelombang magnetik transversal (TM). Efeknya hampir sama seperti di bola logam [14]. Kedua, rumus yang mirip dengan hukum Snell dapat diturunkan dari kondisi

$$ k_{1x}=k_{2x}=k_{2x}^{\mathrm{L}}=k_{1}\,\sin\theta_{1} $$ (24)

di mana indeks 1x dan 2x menunjukkan x -proyeksi vektor di media 1 dan 2, masing-masing, θ 1 adalah sudut datang.

Koefisien Refleksi dan Transmisi

Mari kita tentukan medan yang dibentuk oleh insiden gelombang elektromagnetik terpolarisasi bidang p pada permukaan logam datar. Lebih mudah untuk mengekspresikan komponen medan listrik dan magnet melalui x komponen E (i ) , yaitu E x (r)=− r E x (i) untuk gelombang pantul, E x (t)=t E x (i) untuk gelombang transversal yang ditransmisikan, dan

$$ E^{\mathrm{(a)}}_{x}=\delta\,E^{\mathrm{(t)}}_{x}=t_{\mathrm{L}}\,E^ {\mathrm{(i)}}_{x} $$ (25)

untuk gelombang longitudinal yang ditransmisikan, di sini r adalah koefisien refleksi, t dan t L adalah koefisien transmisi.

Dari kondisi batas Maxwell dan ABC dari Persamaan. (23) ditulis dalam bentuk berikut

$$ \hat{\mathbf{z}}\,(\mathbf{D}-\epsilon_{0}\epsilon_{\mathrm{g}}\,\mathbf{E})|_{z=0}=0, $$ (26)

kami mendapat

$$\begin{array}{*{20}l} r&=- \frac{(1+\delta)\,\epsilon_{1}\,k_{2z}-\epsilon_{2}\,k_{1z }}{(1+\delta)\,\epsilon_{1}\,k_{2z}+\epsilon_{2}\,k_{1z}}=1-(1+\delta)\,t \end{ array} $$ (27) $$\begin{array}{*{20}l} t&=\frac{2\,\epsilon_{1}\,k_{2z}}{\epsilon_{2}\,k_ {1z}+(1+\delta)\,\epsilon_{1}\,k_{2z}}, \end{array} $$ (28) $$\begin{array}{*{20}l} \ delta&=\frac{\epsilon_{\mathrm{g}}-\epsilon}{\epsilon_{\mathrm{g}}}\,\frac{k_{2x}^{2}}{k_{2z}\, k_{2z}^{\mathrm{L}}} \end{array} $$ (29)

Di δ =0, koefisien r menjadi koefisien Fresnel refleksi dari gelombang terpolarisasi p (lihat, misalnya, Persamaan (2.49) dari [4]). Dalam kondisi yang sama, t bukan koefisien transmisi Fresnel karena definisi kami tentang t dan r berbeda dari Fresnel.

Kepunahan Cahaya oleh Nanosfer Logam

Dalam makalah sebelumnya, salah satu penulis menggeneralisasi teori Lorentz-Mie yang memungkinkan persamaan ABC. (23). Sebuah analog dari koefisien Fresnel r , koefisien Mie b l untuk mode TM yang dipantulkan dari l urutan th ditemukan

$$ b_{l}=- \frac{(1+\delta_{l})\,\epsilon_{1}\, \frac{k_{2}\,\psi_{l}^{\prime}(k_ {2} R)}{\psi_{l}(k_{2} R)} - \epsilon_{2}\, \frac{k_{1}\,\psi_{l}^{\prime}(k_{ 1} R)}{\psi_{l}(k_{1} R)}}{(1+\delta_{l})\,\epsilon_{1}\, \frac{k_{2}\,\psi_ {l}^{\prime}(k_{2} R)}{\psi_{l}(k_{2} R)} - \epsilon_{2}\, \frac{k_{1}\,\zeta_{ l}^{\prime}(k_{1} R)}{\zeta_{l}(k_{1} R)}}, $$ (30)

dimana

$$ \delta_{l}=\frac{\epsilon^{\mathrm{T}}- \epsilon_{\mathrm{g}}}{\epsilon_{\mathrm{g}}}\, \frac{l\ ,(l+1)\,j_{l}(k_{2} R)\,j_{l}(k_{2}^{\mathrm{L}} R)}{\psi_{l}^{\ prime}(k_{2} R)\,k_{2}^{\mathrm{L}} R \,j_{l}^{\prime}\left(k_{2}^{\mathrm{L}} R\kanan)}, $$ (31)

ψ l dan ζ l adalah fungsi Riccati-Bessel dan Riccati-Hankel dari ordo l , masing-masing; j l adalah fungsi Bessel berbentuk bola, prima menunjukkan turunan dari suatu fungsi sehubungan dengan argumennya.

Mari kita bandingkan prediksi teori Lorentz-Mie klasik dan umum dengan data eksperimen. Dalam [15], Hilger, Tenfelde, dan Kreibig mempelajari spektrum kepunahan nanopartikel perak yang disimpan pada permukaan dielektrik. Pada penelitian tahap pertama, peneliti membuat berkas partikel perak dengan diameter rata-rata 2, 3,5, dan 4 nm, menentukan distribusi ukuran partikel untuk salah satu berkas, merekam spektrum kepunahan, dan menaksir parameter A =0,25 rumus fenomenologis Γ =Γ b +A v B /R , di mana Γ b adalah tingkat relaksasi logam curah, untuk bola perak dalam ruang hampa. Pertama, kami menghitung spektrum kepunahan untuk seberkas bola perak dengan diameter rata-rata D =2 nm dan distribusi ukuran eksperimental yang membentang wilayah dari D =1 sampai D =4nm. Teori kami tidak mengandung parameter yang dapat disesuaikan. Untuk menentukan fungsi dielektrik, kami menggunakan tabulasi indeks bias perak curah yang diusulkan oleh Lynch dan Hunter [16] (lihat Gambar 1). Kami juga menerapkan Persamaan. (16), (21), dan (22) dengan ω p =9,17 eV, Γ b =0,021 eV, v B =1,39×10 6 m/s, dan A =0,25. Hasil perhitungan dan spektrum eksperimen disajikan pada Gambar 2.

Nyata (a ) dan imajiner (b ) bagian dari fungsi dielektrik perak menurut Johnson dan Christy (◇) [20], Lynch dan Hunter (+) [16], Weber (△) [21], Hao dan Nordlander (garis putus-putus ) [22], dan Drachev et al. (garis padat ) [23]

Spektrum kepunahan cahaya oleh partikel berukuran nanometer perak diamati pada [15] dan dihitung dengan model lokal dan nonlokal. Semua spektrum teoretis disajikan dalam unit relatif umum

Spektrum teoritis pada Gambar. 2 dihitung menggunakan Klimontovich-Silin-Lindhard dan fungsi dielektrik hidrodinamika yang lebih sederhana. Mengejutkan bahwa kedua perhitungan memberikan hasil yang mendekati meskipun |a |>1 di wilayah resonansi plasmon.

Untuk bola perak berukuran nanometer, maksimum dalam spektrum kepunahan, yang disebut resonansi Fröhlich [17], plasmon, dan plasmon polariton permukaan (SPP) [15], diketahui bergeser dari 3,5 menjadi 3,65 eV [18]. Model nonlokal sangat sesuai dengan data eksperimen, sedangkan teori lokal (Mie) memberikan nilai maksimum pada ω 3.5 eV (lihat Gambar 2 dan Tabel 1).

Perhitungan pergeseran biru dari resonansi plasmon dapat didukung oleh pertimbangan berikut. Dalam pendekatan elektrostatik, hanya b 1 berkontribusi pada penampang kepunahan Q eks dan Persamaan. (30) dapat disederhanakan dengan menggunakan pendekatan berikut

$$ \frac{k_{2} R\,\psi_{l}^{\prime}(k_{2} R)}{\psi_{l}(k_{2} R)}\simeq l+1; \,\,\, \frac{k_{1} R\,\zeta_{l}^{\prime}(k_{1} R)}{\zeta_{l}(k_{1} R)} \simeq -\,l. $$ (32)

Jadi, Q eks memiliki maksimum pada

$$ \Re [2\,(1+\delta_{1})\,\epsilon_{1}+ \epsilon_{2}]=0. $$ (33)

Kondisi yang diperoleh (33) memperhitungkan eksitasi L mode (dengan istilah δ 1 ) dan, oleh karena itu, berbeda dari kondisi resonansi Fröhlich [17]:.

$$ \Re (2\,\epsilon_{1}+ \epsilon_{2})=0. $$ (34)

Dalam percobaan [15], frekuensi puncak ω m dan lebar resonansi Δ ω spektrum kepunahan hampir tidak tergantung pada D . Fitur Δ . ini ω tampaknya tidak setuju dengan teori Mie klasik. Sungguh, teori lokal memprediksi perluasan resonansi plasmon dengan penurunan D (di A =0,25) seperti yang ditunjukkan pada Tabel 1. Pada saat yang sama, teori nonlokal memberikan lebar resonansi yang kira-kira sama tetapi posisi puncaknya berbeda. Superposisi kontribusi dari semua partikel memberikan nilai Δ ω yang sangat sesuai dengan data eksperimen. Sangat menarik bahwa teori nonlokal memprediksi perluasan resonansi plasmon sinar bahkan pada A =0.

Di ω>4 eV, kurva teoritis halus pada Gambar. 2 terletak lebih tinggi dari tumbukan puncak eksperimental sempit yang terletak dekat. Penyerapan antar pita mendominasi dalam rentang spektral ini seperti yang dapat dikonfirmasikan oleh Gambar 1. Kekhususan spektrum yang diamati kemungkinan besar merupakan konsekuensi dari transisi dari pita kontinum ke struktur tingkat diskrit. Efek ukuran kuantum seperti itu ditemukan sebelumnya dalam studi tentang sifat optik nanosfer emas [19]. Ketika ukuran bola perak ditingkatkan menjadi D =3,5 nm, penyerapan pertama meningkat relatif terhadap maksimum dan membentuk dataran tinggi dengan serangkaian dips kecil yang berjarak sama. Kemudian, penyerapan sedikit menurun pada D =4 nm.

Untuk mempelajari pembentukan sayap biru resonansi plasmon, kami menghitung spektrum kepunahan partikel perak ultra-timah dan menyajikannya pada Gambar. 3. Fitur luar biasa pada Gambar 3 adalah hilangnya resonansi plasmon sepenuhnya pada D =1nm. Sebelumnya, efek ini diamati dalam studi eksperimental nanospheres emas [19]. Secara khusus, pada Gambar. 9 dari [19], spektrum eksperimental partikel dengan diameter 1.7, 1.9, 2.0, 2.1, 2.3, dan 2,5 nanometer dibandingkan dengan spektrum dihitung dengan teori Mie lokal. Kesepakatan itu buruk, gagal untuk menggambarkan perluasan resonansi plasmon dan posisinya [19]. Upaya untuk meningkatkan kecocokan dengan memvariasikan ukuran partikel dan modifikasi fungsi dielektrik tidak berhasil. Menurut penulis [19], pita osilasi kolektif yang lebar atau tertekan yang diamati secara abnormal menolak untuk dipasangi dengan koreksi yang diusulkan dari teori Mie lokal. Seperti dapat dilihat dari Gambar 4, situasi berubah secara dramatis jika teori Mie nonlokal diterapkan. Perhatikan bahwa kami tidak menggunakan parameter yang dapat disesuaikan. Tabulasi indeks bias kompleks oleh Johnson dan Christy [20] digunakan untuk menentukan fungsi dielektrik emas. Parameter lainnya, termasuk A =1 dan indeks bias toluena (1,37) diambil dari [19].

Penampang kepunahan ternormalisasi partikel perak dengan diameter 2,2, 1,8, 1,4, dan 1,0 nm dihitung dengan lokal (garis putus-putus ) dan nonlokal (kurva padat ) Teori mie. Semakin kecil partikelnya, semakin rendah kurvanya. Semua penampang teoretis disajikan dalam satuan relatif umum

Spektrum serapan dihitung dengan lokal (garis putus-putus ) dan nonlokal (garis padat ) Teori mie dan data eksperimen (titik ) diekstraksi dari Gambar 9 dari [19] untuk bola emas dengan D =2.5, 2.1, dan 1.7 nm dalam toluena (kurva 1 dan lingkaran , kurva 2 dan persegi , dan kurva 3 dan segitiga , masing-masing). Semua spektrum teoretis dinormalisasi menjadi satu pada 4,12 eV dan dipindahkan secara vertikal

Jumlah Gelombang dari Gelombang Membujur

Modus longitudinal berbeda dari yang transversal dengan nilai bilangan gelombang yang jauh lebih tinggi. Misalnya, untuk perhitungan yang disajikan pada Gambar 2, bagian real \(k_{2}^{\mathrm {L}}\) sesuai dengan periode spasial \(\Lambda =2 \pi /\Re k_{ 2}^{\mathrm {L}}\) menurun dari 9 menjadi 2 nm di ω meningkat dari 3 menjadi 4 eV. Dalam ω . ini interval, nilai absolut rasio \(k_{2}^{\mathrm {L}}/k_{2}\) menurun dari 130 menjadi 100 dan parameter δ dari Persamaan. (27) menurun dari 0,01 menjadi 0,005 di θ 1 =π /4. Oleh karena itu, kami menyimpulkan bahwa eksitasi gelombang L pada permukaan perak datar dapat diabaikan. Namun, L mode telah ditemukan penting dalam cluster perak berukuran nanometer.

Penggantian istilah \(-\,\omega ^{2}_{\mathrm {p}}/(\omega ^{2}+i \Gamma \omega)\) dalam Persamaan. (16) oleh ε B ε g menurut Persamaan. (21) memungkinkan kita untuk menulis ulang Persamaan dispersi. (15) dalam bentuk berikut

$$ 1+\frac{3}{5}\,a^{2}+\frac{3}{7}\,a^{4}+\frac{3}{9}\,a^{6 }+\dots=\frac{1}{1-\epsilon^{\mathrm{T}}/\epsilon_{\mathrm{g}}}. $$ (35)

Dalam kasus paling sederhana ε g =1 dan Γ =0, Persamaan. (35) memprediksi bahwa logam transparan untuk gelombang transversal dan L di ω>ω p tapi keduanya k L dan k B kompleks di ω <ω p .

Jika padatan transparan, gelombang longitudinal dapat tereksitasi di bawah insiden miring gelombang terpolarisasi p pada permukaan bidang. Ada beberapa fitur yang berbeda dari efek ini. Pertama, gelombang longitudinal dapat dihasilkan pada permukaan datar, sedangkan upaya khusus harus dilakukan untuk membangkitkan polariton plasmon permukaan [4, 5]. Kedua, dalam pola interferensi, intensitas medan elektromagnetik dimodulasi tidak sepanjang tetapi tegak lurus terhadap antarmuka. Oleh karena itu, rongga dapat muncul pada bidang yang sejajar dengan permukaan akibat spalasi padatan. Menurut definisi ω p , kondisi ω>ω p dapat dipenuhi dalam padatan (misalnya, semikonduktor) dengan kepadatan pembawa arus yang rendah. Kami tidak memeriksa kasus ini di sini karena rumus ε L diturunkan untuk gas elektron yang terdegenerasi.

Kesimpulan

Untuk mendefinisikan fungsi dielektrik yang bergantung secara spasial, semua peneliti sebelumnya mempertimbangkan interaksi materi dengan gelombang elektromagnetik bidang. Pendekatan ini tidak konstruktif dan ketat dalam nano-optik ketika medan dilokalisasi dalam rongga dan kondisi batas harus diperhitungkan. Kami telah memecahkan masalah ini dengan menghitung respons media pada medan listrik yang memenuhi persamaan vektor Helmholtz. Fungsi dielektrik dispersif spasial yang diturunkan bergantung pada kuadrat bilangan gelombang, parameter persamaan Helmholtz, tetapi bukan vektor gelombang dari gelombang bidang.

Kami melaporkan koefisien refleksi Fresnel yang dimodifikasi karena eksitasi gelombang longitudinal dalam logam. Generalisasi serupa dibuat sebelumnya untuk koefisien Mie. Di sini, teori tersebut telah diverifikasi dengan simulasi kepunahan cahaya oleh kluster perak dan emas berukuran nanometer. Pergeseran yang dihitung dari 3,5 ke 3,65 eV dan lebar resonansi plasmon permukaan dari berkas partikel perak sangat sesuai dengan data eksperimen. Selain itu, model nonlokal menjelaskan hilangnya resonansi plasmon bola emas dengan diameter sekitar 2 nm. Adalah penting bahwa gelombang L dapat tereksitasi pada permukaan datar oleh gelombang datang bidang. This is the main difference of the plasmon polaritons from the surface plasmon polaritons.

The properties of the electromagnetic oscillations in metals have been examined. It has been found that the absolute values of the wavenumbers of the longitudinal waves are much larger than those of the transverse waves. For example, in silver at a photon energy of 3.5 eV, the ratio of the absolute values of the wavenumbers is equal to 130. There, the real part of the wavenumber of the longitudinal wave corresponds to a wavelength of 7 nm. The large difference in the wavenumbers prevents excitation of the L waves at a planar surface. However, the L modes have been shown to be excited in silver and gold nanometer-sized particles.


bahan nano

  1. Jenis Logam Merah dan Perbedaannya
  2. Preparasi dan Sifat Magnetik dari Nanopartikel Spinel FeMn2O4 Kobalt-Doped
  3. Sifat-sifat Logam, Nonlogam, dan Metaloid
  4. 20 Berbagai Jenis Logam Dan Sifatnya
  5. Panduan Definitif untuk Rem Elektromagnetik dan Penggunaannya dalam Peralatan Manufaktur
  6. Logam Tahan Api:Sifat, Jenis dan Aplikasi
  7. Bagaimana Logam Paduan Meningkatkan Sifatnya
  8. Berbagai Jenis Logam Nonferrous dan Kegunaannya
  9. Berbagai jenis logam dan klasifikasinya
  10. Jenis-jenis logam dan sifat-sifatnya