Advanced
Theoretical Study on the Effects of Particle Size Distribution on the Optical Properties of Colloidal Gold
Theoretical Study on the Effects of Particle Size Distribution on the Optical Properties of Colloidal Gold
Journal of the Korean Chemical Society. 2007. Oct, 51(5): 407-411
Copyright © 2007, The Korean Chemical Society
  • Received : August 10, 2007
  • Published : October 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에서 28 나노미터인 경우 투과전자현미경 분석 결과와 잘 일치하였다.
Keywords
서 론
외부 전자기장으로 인해 나노 입자의 전자들이 집합적으로 진동하는 표면 플라즈몬 공명(Surface Plasmon Resonance, SPR) 1 현상은 센서나 분자전자소자 등 나노 과학기술 전반에서 폭넓게 응용되고 있다. 나노 입자의 SPR 측정은 대부분의 경우 많은 개수의 입자들의 집합체에 대해서 수행된다. 예를 들면, 나노 입자 콜로이드 용액의 소광 스펙트럼(extinction spectrum)을 얻는 경우, 스펙트럼은 여러가지 크기를 가지는 입자들의 집합적인 광학 현상 발현이라고 보는 것이 타당할 것이다. 그럼에도 불구하고 많은 경우의 스펙트럼 해석은 나노 입자가 단 하나의 크기만을 가진다는 가정에 바탕을 두고 있다. 2 이러한 가정의 타당성을 조사하기 위해 우리는 입자크기 분포를 고려한 경우와 고려하지 않은 경우에 대해 소광 스펙트럼을 계산하고 이를 Joo et al . 3 의 실험과 비교하고자 한다. 이를 통해 나노 입자의 크기분포가 콜로이드 용액의 광학 성질에 미치는 영향을 연구한다.
나노 입자의 광학, 촉매, 전기 등 다양한 성질들은 그 크기에 민감하게 의존하므로 입자 크기를 정확히 결정하는 것은 매우 중요한 일이다. 현재 입자 크기를 측정하는 대표적인 비 파괴 방법은 동적 빛 산란법(Dynamic Light Scattering, DLS) 4 이다. DLS법은 나노 입자에 의해 산란된 빛으로부터 입자의 확산 계수를 구하고 스톡스-아인시타인(Stokes-Einstein) 관계식으로부터 입자 크기를 구하는 방법이다. 흥미롭게도 최근에 Mori et al. 5 은 소광 스펙트럼과 X-선 산란의 두 가지 실험을 병행하면 콜로이드 나노 입자의 크기 분포를 얻을 수 있다고 보고하였다. 우리는 본 연구에서 자외-가시광(UV-VIS) 소광 스펙트럼 하나만으로 입자크기의 분포를 도출하는 것이 가능한지 이론적으로 검토하고자 한다. 이러한 연구가 성공적으로 수행된다면 입자 크기를 측정하는 아주 편리한 신기술 개발로 이어질 수 있을 것이다.
이론 및 전산 방법
나노 입자의 소광 스펙트럼 계산은 미(Mie) 이론 6 을 사용하였다. 미 이론은 구형 입자에 대한 맥스웰(Maxwell) 방정식의 정확한 해를 제공하는데 이에 대한 자세한 내용은 참고문헌 6 으로 미루겠다. 미 이론에 따르면 외부 빛 파장이 λ 로 주어지면 반지름 r 을 가진 구형 입자의 소광 단면적 Cext ( r , λ )은 다음과 같이 주어진다.
PPT Slide
Lager Image
식 (1)에서 α = 2 πr / λ 이고, Re( x )는 복소수 x 의 실수부분을 의미한다. 그리고 aj bj 는 각각 산란되는 전기파와 자기파의 진폭을 나타내며 다음과 같이 표현된다.
PPT Slide
Lager Image
여기에서 β = 이고 m 은 금 나노입자의 굴절률(refractive index) nm 과 용매의 굴절률 n 0 의 비율( m = nm / n 0 , nm 은 복소수)을 나타낸다. ψj ζj ψj ζj 의 일차 미분 값을 의미하고 함수 ψj ζj 는 다음과 같이 주어진다.
PPT Slide
Lager Image
식 (3)에서 Jj 는 베셀(Bessel) 함수이며 Hj (1) 는 1차 행켈(Hankel) 함수이다. 7
PPT Slide
Lager Image
금에 대한 다양한 굴절률 데이터를 사용하여 계산된 소광 스펙트럼들. 금 나노 입자의 반지름이 49 nm인 경우에 대해 Palik(실선), Johnson-Christy(끊어진 선), Morris-Collins(실선 및 점선)이 보고한 굴절률을 사용하여 소광단면적 Cext을 계산하고 이를 파장의 함수로 나타내었다.
식 (1)-(3)을 이용하여 소광 스펙트럼을 계산하기 위해서는 각 파장에 따른 입자(금)와 용매(물)의 굴절률이 필요하다. 물의 굴절률은 Hale과 Querry가 보고한 데이터를 사용 하였다. 8 금의 굴절률은 Morris와 Collins, 9 Johnson과 Christy, 10 Palik 11 이 보고한 데이터들을 모두 사용하여 각각 스펙트럼을 계산한 다음 서로 비교해 보았다( . 1 ). Morris-Collins 데이터의 경우 일부 파장에서 굴절률이 존재하지 않았고 Palik의 굴절률을 사용한 경우 . 1 에 나타난 것처럼 스펙트럼에 인위적인 이중 피크가 나타났다. 그러므로 우리는 이후의 모든 계산에서 Johnson-Christy의 굴절률을 사용하였다.
우리는 금 나노 입자의 반지름 r 을 1 nm부터 1000 nm까지 그리고 빛의 파장 λ 를 300 nm에서 800 nm까지 변화시키면서 주어진 반지름과 파장에서의 소광 단면적 Cext ( λ , r )을 계산하였다. 금 나노 입자의 크기분포는 Joo et al .이 3 투과전자현미경(Tunneling Electron Microscope, TEM) 이미지 분석으로부터 얻은 데이터를 이용하였다. 1 은 7개의 서로 다른 금 콜로이드 용액들에 대해 입자 반지름의 평균
PPT Slide
Lager Image
과 표준편차 σ 를 TEM 이미지들로부터 측정한 값들이다. 3
TEM 이미지로부터 얻은 7개의 금 콜로이드 용액들의 입자 크기 분포
PPT Slide
Lager Image
TEM 이미지로부터 얻은 7개의 금 콜로이드 용액들의 입자 크기 분포
1 에서 입자 반지름이 커질수록 대체로 입자크기 분포의 폭(표준 편차, σ )이 넓어짐을 볼 수 있다. 콜로이드 용액 1-5의 경우 입자크기의 표준 편차가 평균값에 비해 15% 이하이지만 용액 6과7의 경우 입자 반지름 표준 편차가 평균값의 25% 이상임을 알 수 있다.
콜로이드 용액의 입자 반지름 분포 P ( r )은 다음과 같은 가우시안(Gaussian) 분포로 주어진다고 가정하였다.
PPT Slide
Lager Image
식 (4)의
PPT Slide
Lager Image
과 표준편차 σ 를 TEM 이미지들로부터 측정한 값들이다. 3 σ 1 의 값들로 주어진다. 입자크기 분포를 고려한 소광 스펙트럼 Cext ( λ )은 다음과 같이 주어진다.
PPT Slide
Lager Image
식 (5)에서 Cext ( λ , r )은 입자 반지름에 따른 소광 스펙트럼이며, 식 (5)의 수치적분을 위해 사다리 법칙(trapezoidal rule) 이나 심슨(Simpson) 법칙 12 을 사용하였다. 두 가지 수치 적분법들 중 어떤 방법을 쓰든지 비슷한 계산 결과를 얻었다.
실험에서 얻은 소광 스펙트럼 Cexp ( λ )으로부터 입자 크기 분포 P ( r )을 아래와 같이 도출하였다. 우선식 (5)의 좌변 Cext ( λ )을 실험 스펙트럼 Cexp ( λ )으로 대체하고 반지름 r 에 대한 적분을 다음과 같이 격자화 하였다.
PPT Slide
Lager Image
식 (6)에서 Δ r 은 반지름의 격자 간격이고 N ×Δ r 은 계산에서 고려한 최대 입자 반지름이다. 우리는 Δ r =1 nm, N = 1000을 계산에 사용하였다. 주어진 실험 스펙트럼 Cexp ( λ )와 이론으로 계산된 입자 크기에 따른 스펙트럼 Cext ( λ , r )을 이용하여 식 (6)을 P ( r )에 관해서 푸는 것은 매우 어려운 수학적 문제로 알려져 있다. 12 그 이유는 주어진 Cext ( λ , r )에 대해서 식 (6)을 좌변을 만족시키는 P ( r )의 해가 무수히 많이 존재하고 좌변의 Cexp ( λ )가 조금만 변해도(노이즈 등으로 인해) P ( r )의 해가 큰 폭으로 변할 수 있기 때문이다.
우리는 일반적인 P ( r )에 대해 식 (6)을 푸는 대신 P ( r )이 식 (4)와 같은 가우시안 분포로 주어진다고 가정하고 최소자승법(least squares method)을 12 써서 P ( r )을 구하였다. 즉, 다음과 같이 주어지는 자승차를 최소화하는 σ
PPT Slide
Lager Image
을 수치적으로 구하였다.
PPT Slide
Lager Image
식 (7)에서 λ max λ min 은 각각 실험 스펙트럼에서 측정한 최대와 최소 파장을 말한다. 실제 계산에서는 식(7)에 나타난 파장에 관한 적분을 식 (6)과 같이 격자화하고 파장의 격자 간격은 0.5 nm로 잡았다.
결과 및 토론
입자 크기 분포가 소광 스펙트럼에 어떤 영향을 미치는지 알아보기 위하여 1 에 나타난 7개의 콜로이드 용액들에 대해 소광 스펙트럼을 계산하였다( . 2 ). 입자 평균 반지름이 22 nm 이상인 경우, 입자 크기의 분포를 고려한 이론 스펙트럼(실선 및 점선)이 크기 분포가 없는 경우의 이론 스펙트럼(끊어진선)과 그 모양이나 피크 위치에 있어서 확연히 다름을 알 수 있다. 그리고 크기 분포를 고려한 이론 스펙트럼이 실험(실선)에 더 근접함을 알 수 있다.
입자가 작은 경우(반지름 13 nm 미만)에는 입자 크기 분포가 있는 경우와 없는 경우의 스펙트럼이 큰 차이가 없으나 스펙트럼 피크의 위치가 조금 씩 다른 것은 . 2 에서 쉽게 확인할 수 있다. 대체로 입자 크기 분포가 있는 경우의 피크 위치가 크기 분포가 없는 경우 보다 상대적으로 긴 파장에 있는 것도 알 수 있다.
PPT Slide
Lager Image
입자 크기 분포가 소광 스펙트럼에 미치는 영향. Table 1에 나타난 금 나노입자 크기 분포로부터 자외-가시 소광 스펙트럼을 이론적으로 계산하였다(실선 및 점선). 끊어진 선으로 나타난 것은 입자 크기 분포가 없는 경우의 이론 스펙트럼이며 실선으로 나타난 것은 실험 스펙트럼이다.
PPT Slide
Lager Image
금 나노 입자의 평균 반지름 변화에 따른 스펙트럼피크 위치 변화. 표 1에 나타난 7개의 콜로이드 용액에 대해 입자의 평균 반지름과 소광 스펙트럼의 피크 파장 사이의 관계를 나타냈다. 입자 크기 분포를 고려한 계산(X표)과 고려하지 않은 계산(동그라미) 결과를 실험 결과(네모)와 비교하였다.
스펙트럼의 피크 위치와 입자 크기 분포와의 관계를 정량적으로 보여주기 위해 . 3 에 실험(네모), 입자 크기 분포를 고려한 계산(X표), 크기 분포를 고려치 않은 계산(동그라미)으로부터 얻은 피크 위치를 평균 입자 반지름의 함수로 나타내었다. 대체로 입자 크기 분포를 고려한 계산에서 나온 피크 파장이 분포를 고려치 않은 계산에서 나온 피크 파장보다 상대적으로 긴 파장을 가지며 실험에 더 가까운 것을 보여주고 있다.
. 3 에서 입자의 평균 반지름이 17 nm 이상인 경우, 입자 크기분포를 고려한 계산과 고려치 않은 계산이 스펙트럼의 피크 위치를 예측하는데 상당한 차이를 나타낸다는 것을 알 수 있다. 그리고 입자크기의 분포를 고려하더라도 이론적인 계산이 실험과 완벽하게 일치하지는 않는다는 것도 알 수 있다. 실험과의 더 나은 일치를 위해서는 아마도 우리는 입자 반지름의 분포를 가우시안(식 (4))으로 가정하지 않고 다른 종류의 분포를 이용해야 될 것으로 보인다. 그러나 입자 평균 반지름이 49 nm인 경우를 제외하고는 가우시안 분포 만으로도 입자 크기를 고려하지 않은 경우에 비해 실험에 상당히 가까워진다는 것을 알 수 있다.
. 4 는 위에서 설명한 최소자승법을 써서 얻은 입자 반지름의 평균[위 그림]과 표준 편차[아래 그림]를 보여준다. 그림에는 TEM 이미지 분석을 통해 얻은 입자 반지름 평균과 표준 편차도 함께 나타내었다( 1 의 값들). 그림에서 보듯이 입자 반지름이 5.4 nm에서 13.5 nm 사이의 경우 최소자승법으로 얻은 입자 크기 분포가 실험과 잘 일치하는 경향을 보인다. 그러나 입자 반지름이 아주 작은 경우(3.2 nm)와 반지름이 큰 경우(17.4 nm 이상) 최소자승법으로 얻은 값이 실험값과 상당한 차이를 보임을 알 수 있다. 이것은 입자가 매우 작은 경우, 입자 반지름 분포의 폭이 매우 좁아 가우시안 분포 가정이 성립하지 않기 때문인 것으로 생각된다. 그리고 입자 반지름이 17 nm 이상인 경우 입자 크기 분포의 폭이 너무 넓어서 또다시 가우시안 분포 가정이 성립하지 않아 최소자승법과 실험값이 큰 차이를 보인다고 생각된다. 하지만 평균 반지름이 5 nm에서 14 nm 사이의 용액은 입자 크기가 가우시안 분포로 잘 근사됨을 알 수 있다.
PPT Slide
Lager Image
소광 스펙트럼으로부터 이론적으로 도출된 입자크기 분포. 입자 크기에 따라 계산된 이론 스펙트럼을 이용하여 입자 반지름의 평균(위 그림)과 표준 편차(아래 그림)를 최소 자승법으로 얻었다. 최소자승법으로 얻은 이론값(X표)과 TEM 이미지 분석으로 얻은 값(동그라미)을 함께 나타냈다.
입자 크기가 가우시안 분포로 주어지지 않는 경우(입자 반지름이 3.2 nm 이하 또는 17.4 nm 이상인 경우) 우리는 일반적인 크기 분포에 대해서 식 (6)을 풀어야 할 것이다. 그러나 앞에서 지적한 대로 이것은 수학적으로 매우 어려운 문제일 수 있다. 이러한 수학적 어려움을 극복하기 위해 우리는 식 (7)에 P ( r )의 해를 안정화시키는 규칙화(regularization) 항을 더하여 이들 두 항의 합을 최소화시키는 방법을 생각해 볼 수 있다. 13 향후 우리는 이러한 규칙화 방법을 이용하여 가우시안이 아닌 일반적인 크기 분포를 실험 스펙트럼으로부터 얻는 방법을 연구할 것이다.
결 론
금 콜로이드 용액의 입자 크기 분포를 고려하여 자외-가시광 소광 스펙트럼을 미 이론으로 계산하고 실험과 비교하였다. 입자 크기 분포를 고려하기 위해 금 나노 입자의 반지름을 1 nm부터 1000 nm까지 변화시키며 각 입자 크기에 따른 소광 스펙트럼을 계산하였다. 그런 다음 TEM 이미지 분석으로 얻은 입자 크기 분포를 이용하여 금 콜로이드의 소광 스펙트럼을 계산하였다. 이러한 입자크기 분포를 고려한 계산은 크기 분포를 무시한 계산보다 실험에 더 가까운 것을 확인하였다. 입자 반지름 분포를 가우시안으로 가정하고 최소자승법을 사용하여 실험 스펙트럼으로부터 입자 반지름 분포를 도출하였다. 이러한 최소자승법은 나노입자 평균 반지름이 5 nm에서 14 nm 사이인 경우 TEM 이미지와 잘 일치함을 확인하였다. 그러나 입자 평균 반지름이 5 nm보다 작거나 14 nm보다 큰 경우에는 실험과 잘 일치하지 않았다. 이는 입자 크기 분포가 너무 좁거나 넓어서 최소자승법의 가우시안 분포 가정이 성립하지 않기 때문인 것으로 생각된다.
이 논문은 부산대학교 자유과제 학술연구비(2년)에 의하여 연구되었음.
References
Kelly K. L. , Coronado E. , Zhao L. L. , Schatz G. C. 2003 J. Phys. Chem. B 107 668 -    DOI : 10.1021/jp026731y
Link S. , El-Sayed M. A. 1999 J. Phys. Chem. B 103 4212 -    DOI : 10.1021/jp984796o
Jang S. , Park J. , Shin S. , Yoon C. , Choi B. K. , Gong M. , Joo S. W. 2004 Langmiur 20 1922 -    DOI : 10.1021/la0348665
Berne B. J. , Pecora R. 2000 Dynamic Light Scattering: With Applications to Chemistry, Biology, and Physics Dover NY
Nakamura K. , Kawabata T. , Mori Y. 2003 Powder Technology 131 120 -    DOI : 10.1016/S0032-5910(02)00340-6
Bohren C. F. , Huffman D. F. 1983 Absorption and Scattering of Light by Small Particles Wiley New York
Abramowitz M. , Stegun I. A. 1970 Handbook of Mathematical Functions Dover New York
Hale G. M. , Querry M. R. 1973 Appl. Opt. 12 555 -    DOI : 10.1364/AO.12.000555
Morris R. H. , Collins L. F. 1964 J. Chem. Phys. 41 3357 -    DOI : 10.1063/1.1725733
Johnson P. B. , Christy R. W. 1972 Phys. Rev. B 12 4370 -    DOI : 10.1103/PhysRevB.6.4370
Palik E. D. 1985 Handbook of Optical Constants of Solids Academic New York
Press W. H. , Teukolsky S. A. , Vetterling W. T. , Flannery B. P. 1986 Numerical Recipes in Fortran Cambridge University Press Cambridge
Tikhonov A. N. , Arsenin V. Y. 1977 Solutions of Ill-Posed Problems Winston & Sons Washington, D. C.