Advanced
Discrete Dipole Approximation Calculation of the Extinction Spectra of Gold-Silver Alloy Nanoparticles
Discrete Dipole Approximation Calculation of the Extinction Spectra of Gold-Silver Alloy Nanoparticles
Journal of the Korean Chemical Society. 2007. Apr, 51(2): 136-140
Copyright © 2007, The Korean Chemical Society
  • Received : February 15, 2007
  • Published : April 20, 2007
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
희미 이
Saha Leton Chandra
준경 장

Abstract
이산 쌍극자 근사를 이용하여 수용액상의 금-은 합금 나노입자의 표면 플라즈몬 공명 스펙트럼을 계산하였다. 직경 10 나노미터의 금-은 합금 입자의 경우에 스펙트럼의 최대 소광 파장이 합금의 은 성분이 높아짐에 따라 선형적으로 짧은 파장대로 이동하며 최대 소광 세기는 지수적으로 증가함을 관측하였다. 이러한 계산결과는 실험 결과들과 잘 일치하는 것을 확인하였다.
Keywords
서 론
최근 금속 나노입자의 모양, 크기, 주변환경(입자의 표면 배치, 입자-유무기물 혼성체) 등을 조절하는 기술과 나노입자에 기능성을 부여하여 검출하고자 하는 분자(예: DNA, 단백질)와 선택적으로 반응하게 하는 기술이 크게 발전하였다. 1 , 2 용액상 합성 방법으로 막대 및 삼각형 등 비구형 나노입자를 합성할 수 있게 되었고 3 다양한 리소그래피(lithography) 기술을 이용하여 원하는 크기와 모양의 입자들을 기판에 배치할 수 있게 되었다. 4 이러한 기술들을 결합하여 만들어진 화학 및 생물 센서는 의료 진단, 환경 처방, 독극물 검출 등에 있어서 기존의 센서에 비해 높은 감도를 가지는 것으로 알려져 있다. 나노 센서는 작은 크기로 인한 높은 집적도 때문에 많은 분자들을 동시에 감지해 낼 수 있는 집적 센서로의 발전 가능성도 높다 하겠다.
다양한 금속들 중에서 특히 은(Ag) 나노입자는 센서로 사용되었을 경우 정확하고 높은 감도를 나타내는 것으로 알려져 있다. 하지만 은은 화학적으로 불안정하여 액체상에서의 센서로 쓰이기에는 적당치 못하다. 5 최근에는 은과 화학적으로 안정한 금을 합금(alloy)으로 만들어 생물 및 화학 센서로 사용될 수 있다는 연구결과들이 보고 되었다. 6 , 7 지금까지 다양한 조성을 가진 금-은 합금 입자들이 콜로이드 용액, 박막, 또는 매트릭스(matrix) 형태로 합성되었다. 8 - 19 콜로이드 상의 금-은 합금 나노 입자의 경우 스펙트럼의 최대 소광(extinction) 파장이 은 성분이 증가함에 따라 선형적으로 감소함이 관측되었다. 9 , 11 , 12 Papavassiliou 12 는 직경 10 nm의 콜로이드 합금 입자의 경우 은 함량이 높아짐에 따라 광학 스펙트럼의 최대 소광 파장이 약 540 nm(순수한 금)에서 400 nm(순수한 은)로 선형적으로 감소함을 보고하였다. Chen과 9 El-Sayed 11 도 최대 흡수 파장과 은 성분 사이에 비슷한 관계를 관측하였다. 한편 El-Sayed 11 는 스펙트럼의 최대 소광 세기(피크 높이)는 은 성분에 대해 지수적인 (exponential) 의존성을 가짐을 보고하였다.
금속입자의 합성 및 광학 실험에서는 흔히 불순물, 입자 크기 및 모양의 불균일 등이 나타난다. 예를 들면 실험에서 합금 성분을 바꾸면 입자의 크기도 함께 변화하여 순수한 성분 변화에 따른 소광 스펙트럼 변화를 관찰하기 힘들다. 17 이론 연구는 이러한 실험적 난관들을 없애고 입자 크기 및 모양을 완벽하게 제어할 수 있는 장점이 있다. 금속 입자의 광학적 성질은 표면 플라즈몬 공명(Surface Plasmon Resonance, SPR) 현상 때문인 것으로 잘 알려져 있다. SPR은 외부 전자기장에 의해 금속입자의 전도성 전자들이 원자핵으로부터 일정한 진동수를 가지고 집합적으로 움직이는 현상을 말한다. 20 최근 전산 전자기학의 발전으로 인해 나노 입자의 플라즈몬 공명에 대한 이론적 연구가 가능하게 되었다. 본 연구는 이산 쌍극자 근사(Discrete Dipole Approximation, DDA) 21 를 이용하여 금-은 합금 입자의 성분 변화에 따른 표면 플라즈몬 공명 스펙트럼 변화를 고찰하고자 한다. 다양한 전산전자기학 방법들 22 - 24 중 DDA를 선택한 것은 이 방법이 임의의 크기와 모양을 가진 나노 입자의 광학 성질 연구에 적합하여 향후 다양한 모양의 입자에 대해 적용이 용이하기 때문이다. 합금입자에 대한 DDA 계산을 실험과 비교하여 이 방법의 정확성을 가늠하고자 한다.
전산 방법
DDA에 대한 자세한 내용은 참고문헌 20 21 으로 미루고 여기서는 대략적으로만 설명하고자 한다. DDA는 금속입자를 N 개의 정육각형 격자로 나타내고 각 격자에는 외부 전자기장에 의해 유도된 쌍극자가 위치한다고 가정한다. 이 방법은 임의의 모양과 조성을 가진 입자를 자유롭게 구현할 수 있다는 장점을 지니고 있다.
PPT Slide
Lager Image
에 위치한 j 번째 쌍극자는 이 위치에서의 전기장이
PPT Slide
Lager Image
이라면 다음과 같은 분극(Polarization)
PPT Slide
Lager Image
를 가진다.
PPT Slide
Lager Image
. α j 는 분극성(Polarizability)이며 Draine과 Flatau는 주어진 진동수 ω에서의 나노 입자의 유전체 함수 ε(ω)(이 값은 복소수 값임)와 α j 사이의 관계를 유도하였다. 21 한편
PPT Slide
Lager Image
는 외부에서 걸어준 전기장
PPT Slide
Lager Image
j 번째 쌍극자를 제외한 나머지 N-1개의 쌍극자에 의한 유도 전기장의 합으로서 식 (1)과 같이 주어진다.
PPT Slide
Lager Image
식 (1)에서
PPT Slide
Lager Image
는 다음과 같이 주어지는 3×3 매트릭스이다.
PPT Slide
Lager Image
위 식에서 k = ω / c ,
PPT Slide
Lager Image
, 그리고
PPT Slide
Lager Image
이다. 만일 우리가
PPT Slide
Lager Image
= 1/α j 라고 정의하면
PPT Slide
Lager Image
에 관한 다음과 같은 선형방정식을 얻는다.
PPT Slide
Lager Image
위의 식을
PPT Slide
Lager Image
에 관해서 풀면 소광 단면적(extinction cross section) Cext 는 다음 식으로 주어진다.
PPT Slide
Lager Image
소광 스펙트럼 계산에 있어서는 금-은 합금 입자를 반지름 r 이 5 nm인 구라고 가정하였다. 외부 빛의 파장을 200 nm에서 700 nm까지 5 nm씩 증가시키면서 주어진 파장에서의 소광 단면적 Cext 을 계산하였다. 실제 스펙트럼은 Cext 를 구의 단면적 πr 2 으로 나눈 소광 효율(extinction efficiency)값을 나타내었다( . 1 참조). DDA 계산에는 총 739개의 쌍극자들이 포함되었고 DDA계산에 필요한 물과 합금 입자의 유전체 함수는 다음과 같이 구하였다. 우선 물의 유전체 함수는 Hale과 Querry가 보고한 값을 사용하였다. 25 그리고 합금 입자의 유전체 함수 ε(ω)(복소수 값을 가짐)는 다음과 같이 은의 몰성분(mole fraction) XAg 에 대해 선형적인 관계식을 이용하였다.
PPT Slide
Lager Image
은성분 변화에 따른 수용액 상의 금-은 합금 입자의 소광 스펙트럼 변화. 은의 몰성분을 0%, 10%, 30%, 50%, 70%, 90%, 100%로 증가시킴에 따라 스펙트럼이 짧은 파장쪽으로 이동하고, 스펙트럼 폭이 좁아지며, 피크의 높이가 올라감을 볼 수 있다.
PPT Slide
Lager Image
위 식에서 ε Ag (ω)와 ε Au (ω)는 각각 주어진 주파수 ω에서의 은과 금의 유전체 함수들이며 참고 문헌 26 에 보고된 값을 사용하였다. 식 (3)의 해를 구하기 위해 PBCGST를 27 이용한 반복적인(iterative) 방법을 사용하였고, 상대 오차 10 -5 이하를 수렴조건으로 사용하였다.
결과 및 토론
. 1 은 수용액상의 다양한 성분을 가진 금-은 합금(직경 10 nm)의 소광 스펙트럼을 위에 기술한 DDA 방법으로 계산한 결과이다. 우선 모든 스펙트럼은 두개가 아닌 하나의 플라즈몬 피크를 가짐을 알 수 있다. 그리고 은 성분이 증가할수록 소광 스펙트럼의 최대값이 나타나는 파장이 짧은 쪽으로 이동하는 것을 볼 수 있다. 또한 순수한 은에 가까울수록 스펙트럼의 폭이 좁아지면서 뚜렷한 플라즈몬 공명현상을 보임을 알 수 있다. 이로부터 합금입자의 은 성분이 증가하면 날카로운 소광 스펙트럼을 가져서 센서로서의 감도가 더욱 좋아질 것이라 예측할 수 있다. 그러나 앞서 지적한 대로 순수한 은은 화학적으로 불안정하기 때문에 금 성분을 첨가함으로써 화학적인 안정성을 개선할 수 있을 것이다. 그림에서 은의 몰성분이 90%인 경우 순수한 은에 비해 감도가 크게 떨어지지 않음을 볼 수 있다.
PPT Slide
Lager Image
금-은 합금 나노 입자의 최대 소광 파장(Peak position)과 은의 몰성분(%Ag)의 관계. 금-은 합금 입자에서 은성분을 변화시키면서 스펙트럼의 최대 소광 파장(원으로 표시됨)을 DDA 방법으로 계산하였다. 은의 성분이 높아짐에 따라 최대 소광 파장은 약 520 nm에서 약 410 nm로 감소함을 볼수 있다. 직선으로 표시된 것은 DDA에 대한 선형 fitting 함수이다.
. 2 는 금-은 합금입자의 플라즈몬 공명 파장이 은 성분의 변화에 따라 어떻게 변화하는지를 정량적으로 보여주고 있다. . 2 를 구하기 위해 우리는 . 1 의 각 스펙트럼에 대해 최대의 세기를 주는 파장을(그림에는 Peak position으로 나타남) 구하여 이를 은의 몰성분(%Ag)의 함수로 나타내었다. Papavassiliou 12 는 직경 10 nm의 금-은 합금 입자(콜로이드)에 대한 실험을 통해 은 성분이 증가함에 따라 피크 파장이 선형적으로 긴 파장에서 짧은 파장으로 이동한다는 것을 보고하였다. El-Sayed 11 도 직경이 약 20 nm의 금-은 합금 입자에 대해 비슷한 선형적 관계를 관측하였다. 이러한 실험 결과들과 유사하게 . 2 에 나타난 DDA 계산 결과도 스펙트럼의 피크 파장이 은 성분에 대해 선형에 가까운 의존성을 가짐을 보여준다. 그림에 직선으로 나타난 것은 원들로 표시된 DDA 계산 결과들을 선형 방정식으로 fitting한 것이며 Fitting 함수는 Peak position = 532 nm - 1.15×(%Ag)으로 주어진다.
. 1 에 볼 수 있는 또 한가지 특징은 은 성분이 증가함에 따라 스펙트럼의 피크 높이가 증가한다는 것이다. Elsayed는 11 직경 20 nm의 금-은 합금 나노입자의 경우 소광 스펙트럼의 최대 소광 세기는 은 성분이 증가함에 따라 지수적으로(exponentially) 증가한다는 것을 관측하였다. 이러한 실험 결과와 비교하기 위해 . 1 에 나타난 최대 소광 세기와 은 성분의 상관 관계를 . 3 에 정량적으로 나타내었다(원들로 표시됨). . 3 에 곡선으로 나타난 것은 Levenverg-Marquardt 28 방법을 이용한 비선형 fitting의 결과이다. Fitting 함수식은 Peak intensity = 0.56 + 0.0237*exp[0.0428 × (%Ag)]으로 주어진다. 그림에서 볼 수 있듯이 El-Sayed의 보고와 비슷하게 은 성분과 최대 소광 세기는 지수함수에 가까운 비선형적인 관계를 보이고 있다. 그러므로 우리의 DDA 계산이 El-Sayed의 실험의 주요 결론과도 일치함을 알 수 있다.
. 3
PPT Slide
Lager Image
은 성분 변화에 따른 최대 소광 세기 변화(원들로 표시됨). 직경 10 nm의 금-은 합금 입자에서 은의 몰 성분(%Ag)을 0%, 10%, 30%, 50%, 70%, 90%, 100%로 바꾸면서 최대 소광 세기(Peak intensity)를 계산하였다. 곡선으로 나타난 것은 지수함수 형태의 fitting 결과이다.
결 론
본 연구에서 우리는 금-은 합금 나노 입자에서 순수한 조성 변화에 따른 플라즈몬 공명 스펙트럼의 추이를 이론적으로 살펴보았다. 금과 은의 유전체 함수를 조성비율로 평균한 값(식 (5))을 합금 입자의 유전체 함수로 쓰고 DDA 방법을 이용하여 직경 10 나노미터의 금-은 합금 입자의 수용액상의 소광 스펙트럼을 계산하였다. 이러한 간단한 DDA 계산 결과들은 El-Sayed 11 와 Papavassiliou 12 가 실험으로 얻은 주요 결론들과 잘 일치하는 것을 확인하였다. 즉, 은의 성분이 증가함에 따라 플라즈몬 피크(최대 소광) 파장이 약 540 nm에서 400 nm로 선형적으로 감소하고 최대 소광 세기(피크의 높이)가 지수적으로 증가하는 현상을 이론적으로 재현할 수 있었다. DDA 방법은 임의의 모양 및 크기를 가진 나노입자의 스펙트럼 계산을 하는데 용이하므로 향후 비구형 및 새로운 형태의 합금 나노입자에 대한 표면 플라즈몬 공명 연구에 유용하게 적용될 수 있을 것으로 기대된다.
이 논문은 부산대학교 자유과제 학술연구비(2년)에 의하여 연구되었음.
References
Mirkin C. A. , Letsinger R. L. , Mucic R. C. , Storhoff J. J. 1996 Nature (London) 382 607 -    DOI : 10.1038/382607a0
Boal A. K. , Ilhan F. , DeRoucher J. E. , Thurn-Albrech T. , Rotello V. M. 2000 Nature (London) 404 746 -    DOI : 10.1038/35008037
Hao E. , Bailey R. C. , Schatz G. C. , Hupp J. T. , Li S. 2004 Nano Lett. 4 327 -    DOI : 10.1021/nl0351542
Zhang X. , Yonzon C. , Van Duyne P. 2006 J. Mat. Res. 21 1083 -    DOI : 10.1557/jmr.2006.0136
Sharma A. , Gupta B. D. 2006 Nanotechnology 17 124 -    DOI : 10.1088/0957-4484/17/1/020
Ehler T. T. , Neo L. J. 1995 Langmuir 11 4177 -    DOI : 10.1021/la00010a088
Sharma A. K. , Gupta B. D. 2005 Opt. Commun. 245 159 -    DOI : 10.1016/j.optcom.2004.10.013
Moskovits M. , Srnova-Sloufova I. , Vlckova B. 2002 J. Chem. Phys. 116 10435 -    DOI : 10.1063/1.1449943
Chen D.-H. , Chen C.-J. 2002 J. Mater. Chem. 12 1557 -    DOI : 10.1039/b110749f
Teo B. K. , Keating K. , Kao Y.-H. 1987 J. Am. Chem. Soc. 109 3494 -    DOI : 10.1021/ja00245a070
Link S. , Wang Z. L. , El-Sayed M. A. 1999 J. Phys. Chem. B 103 3529 -    DOI : 10.1021/jp990387w
Papavassiliou G. C. 1976 J. Phys. F: Metal Phys. 6 L103 -    DOI : 10.1088/0305-4608/6/4/004
Mulvaney P. 1996 Langmuir 12 788 -    DOI : 10.1021/la9502711
Kreibig U. 1974 J. Phys. F: Metal Phys. 4 999 -    DOI : 10.1088/0305-4608/4/7/007
Hornyak G. L. , Patrissi C. J. , Oberhauser E. B. , Martin C. R. , Valmalette J-C , Lemaire L. , Dutta J. , Hofmann H. 1997 NanoStructured Materials 9 571 -    DOI : 10.1016/S0965-9773(97)00127-X
Hubenthal F. , Ziegler T. , Hendrich C. , Alschinger M. , Trager F. 2005 Eur. Phys. J. D 34 165 -    DOI : 10.1140/epjd/e2005-00138-1
Roy R. K. , Mandal S. K. , Pal A. K. 2003 Eur. Phys. J. B 33 109 -    DOI : 10.1140/epjb/e2003-00147-x
Gaudry M. , Lerme J. , Cottacin E. , Pellarin M. , Vialle J.-L. , Broyer M. , Prevel B. , Treilleux M. , Melinon P. 2001 Phys. Rev. B 64 085407 -    DOI : 10.1103/PhysRevB.64.085407
Shi H. , Zhang L. , Cai W. 2000 J. Appl. Phys. 87 1572 -    DOI : 10.1063/1.372053
Kelly K. L. , Coronado E. , Zhao L. L. , Schatz G. C. 2003 J. Phys. Chem. B 107 668 -    DOI : 10.1021/jp026731y
Draine B. T. , Flatau P. 1994 J. Opt. Soc. Am. A 11 1491 -    DOI : 10.1364/JOSAA.11.001491
Moreno E , Erni D. , Hafner C. , Vahldieck R. 2002 J. Opt. Soc. A 19 101 -    DOI : 10.1364/JOSAA.19.000101
Novotny L. , Bian R. X. , Xie X. S. 1997 Phys. Rev. Lett. 79 645 -    DOI : 10.1103/PhysRevLett.79.645
Bian R. X. , Dunn R. C. , Xie X. , Leung P. T. 1995 Phys. Rev. Lett. 75 4772 -    DOI : 10.1103/PhysRevLett.75.4772
Hale G. M. , Querry M. R. 1973 Appl. Opt. 12 555 -    DOI : 10.1364/AO.12.000555
Palik E. D. 1985 Handbook of Optical Constants of Solids Academic New York
Flatau P. 1997 J. Opt. Lett. 22 1205 -    DOI : 10.1364/OL.22.001205
Press W. H. , Teukolsky S. A. , Vetterling W. T. , Flannery B. P. 1986 Numerical Recipes in Fortran Cambridge University Press Cambridge