Advanced
Study on the Micellization of CPC/Brij 35 Mixed Surfactant Systems in Water
Study on the Micellization of CPC/Brij 35 Mixed Surfactant Systems in Water
Journal of the Korean Chemical Society. 2009. Apr, 53(2): 118-124
Copyright © 2009, The Korean Chemical Society
  • Received : February 25, 2009
  • Published : April 20, 2009
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
한내 길
병환 이

Abstract
CPC(1-hexadecylpyridinium chloride)와 Brij 35(polyoxyethylene(23) lauryl ether)의 혼합계면활성제의 미셀화에 대한 임계미셀농도(CMC)와 반대이온 결합상수값( B )을 순수 물에서 CPC의 겉보기몰분율(α 1 )의 함수로서 전도도법으로 측정하였다. 이와 같이 측정한 CMC값에 비이상적 혼합미셀 모델을 적용함으로써 여러 가지 열역학적 함수값( Xi , γi , C i , α i M , β 및 ΔH mix )을 계산하고 분석하였다. 또한 온도에 따른 CMC값의 변화로부터 CPC/Brij 35 혼합계면활성제의 미셀화에 대한 열역학 함수값(ΔG o m , ΔH o m 및 ΔS o m )을 계산하였다. 그 결과, ΔG o m 값은 모두 음의 값을 나타내었지만, ΔH o m 와 ΔS o m 의 값은 측정한 온도와 α 1 의 값에 따라 음 혹은 양의 값을 나타내었다.
Keywords
서 론
계면활성제의 생산과정에는 생성물인 계면활성제뿐만 아니라 동시에 미반응의 유기물질들이 다량 함유되어서 실생활에서는 순수 계면활성제보다는 혼합계면활성제가 널리 이용되고 있다. 혼합계면활성제는 대체적으로 단일 성분의 계면활성제가 나타낼 수 없는 특수한 물리적 특성을 나타냄으로써 경제적으로 활발한 연구의 대상이 되고 있다. 1 5 혼합계면활성제들의 미셀화 현상을 고찰하기 위하여 다양한 이론들이 제안되고 있으며, 6 8 또한 열역학적인 고찰을 위하여 온도에 따른 임계미셀농도(CMC)와 반대이온결합상수( B )값의 변화를 측정하고 있다. 10 11
섬유, 염색 및 정밀화학 산업분야에서는 기름 혹은 고분자 성분을 물에 유화시키기 위하여 다양한 종류의 혼합계면활성제가 이용되고 있다. 일반적으로 음이온/비이온 혼합계면활성제가 많이 사용되고 있으며, 특수한 경우에는 양이온/비이온 혼합계면활성제가 이용되고 있다. 양이온 계면활성제인 CPC(1-hexadecylpyridinium chloride)와 비이온계면활성제인 Brij 35(polyoxyethylene(23) lauryl ether)의 혼합계면활성제에 대한 연구는 지금까지 진행된 바 없다. 이들 두 성분의 혼합계면활성제에 대한 연구는 화장품, 페인트, 식품, 의약품과 같은 산업현장에서 안정한 에멀젼의 제조를 위하여 응용될 수 있다. 따라서 본 연구에서는 CPC/Brij 35 혼합계면활성제의 비이상적 혼합 미셀화 현상을 분석하였으며 또한 이 혼합계면활성제의 미셀화 현상에 대한 열역학적인 연구를 수행하고자 한다.
실 험
양이온성 계면활성제인 CPC와 비이온성 계면활성제인 Brij 35는 순도가 98%이상 인 Aldrich 제품을 더 이상 정제하지 않고 바로 사용하였으며, 실험방법은 전편의 논문에서 사용한 방법을 이용하였다. 11 14 즉, 용액 제조를 위하여 우선 농도가 각각 0.01 M인 CPC와 Brij 35의 저장용액들을 만들었으며, 이들 저장용액들을 정해진 CPC의 겉보기 몰분율조성( α 1 )에 따라 일정한 부피비로 혼합함으로써 조성이 고정된 혼합계면활성제의 저장용액들을 만들었다. 제조한 혼합계면활성제의 저장용액을 순수 물로 묽혀서 12개의 용액들을 만들었으며, 각 용액의 전도도 값을 측정하여 CMC와 B 값을 결정하였다. 15 혼합계면활성제의 겉보기몰분율조성( α 1 )은 1(순수 CPC), 0.8, 0.6, 0.4, 0.2 및 0(순수 Brij 35)으로 하였으며, 온도는 288K에서 308K까지 5K 간격으로 변화시키면서 측정하였다.
결과 및 고찰
순수 물에서 CPC/Brij 35 혼합계면활성제의 CMC와 B 값을 온도와 α 1 의 변화에 따라 측정하였으며, 그 결과를 1 에 모두 나타내었다. 각 온도에서 측정한 혼합계면활성제의 CMC값은 순수 CPC와 순수 Brij 35의 CMC값보다 훨씬 큰 값을 나타내었다. 또한 각 몰분율조성에서 CMC값은 온도의 증가에 따라 감소하다가 증가하는 경향을 보였으며, 293K에서 최소값을 나타내었다. 한편 혼합미셀 상에서 반대이온(Cl - )들의 결합상수( B )값은 α 1 의 증가에 따라 증가하는 경향을 보였으며, 온도의 증가에 따라 B 의 값은 거의 변화가 없거나 감소하는 경향을 보였다. 이러한 현상들은 α 1 이 증가하면 미셀 상에서 이온성 계면활성제인 CPC가 많아짐으로써 반대이온이 미셀과 더욱 많이 결합하기 때문이며, 또한 온도의 증가에 따라 반대이온의 활동도가 증가함으로써 미셀과의 결합이 감소하기 때문이다.
Values of CMC(±0.1, mM) andB(±0.02) for the micellization of CPC/Brij 35 mixed surfactant systems in water at various temperatures
PPT Slide
Lager Image
Values of CMC(±0.1, mM) and B(±0.02) for the micellization of CPC/Brij 35 mixed surfactant systems in water at various temperatures
임의의 α 1 에서 측정한 혼합계면활성제의 CMC 값과 순수 CPC와 순수 Brij 35에 대한 임계미셀 농도값인 CMC 1 과 CMC 2 를 비이상적 혼합미셀모델인 식 (1)에 대입하면 혼합미셀 상에서 CPC 분자의 몰분율(X 1 )의 값을 계산할 수 있다. 298K 에서 CPC/Brij 35 혼합계면활성제의 미셀화에 대하여 이와 같이 임의의 α 1 에서 X 1 을 계산하였으며, 그 결과를 2 에 나타내었다. 여기서 보듯이 CPC/Brij 35 혼합계면활성제의 미셀화에서는 X 1 값이 α 1 값보다 훨씬 큰 값을 나타내었으며, 이것은 CMC값에서 유추할 수 있듯이 X 1 값이 이상적 미셀화로부터 큰 양의 벗어남을 나타낸다.
Values of the measured CMC(±0.01) and calculated thermodynamic parameters from the nonideal mixed micellar model for the micellization of CPC/Brij 35 mixed surfactant systems in pure water at 298K
PPT Slide
Lager Image
Values of the measured CMC(±0.01) and calculated thermodynamic parameters from the nonideal mixed micellar model for the micellization of CPC/Brij 35 mixed surfactant systems in pure water at 298K
PPT Slide
Lager Image
계면활성제분자들이 혼합미셀을 이룰 때 두종류의 계면활성제분자들 사이에 소수성 인력 혹은 정전기적 반발력이 생기게 되며, 이러한 상호작용의 세기를 나타내는 상호인력상수(β)값은 X 1 과 CMC값에 의하여 계산할 수 있다. 16 , 17 2 의 X 1 과 CMC값들을 이용하여 임의의 α 1 에서 β값들을 계산하였으며, 그 결과를 2 에 함께 나타내었다. 이들 β값들은 모두 음의 값을 나타내었으며, α 1 의 값이 0.6일 때 β값은 최소값을 나타내었다. 이와 같이 CPC/Brij 35 혼합계면활성제의 미셀화에 대한 β값이 음의 값을 나타내는 것은 두 성분들 사이에 강한 인력이 작용함을 나타내며, 특히 두 성분의 head-그룹들 사이에 이온-쌍극자 상호인력으로 인하여 혼합계면활성제의 미셀화가 이상적인 미셀화로부터 크게 벗어나게 된다.
비이상적 혼합미셀 모델에 의하여 혼합미셀상에서 CPC분자와 Brij 35분자들의 활동도계수를 나타내는 γ 1 γ 2 의 값을 2 의 X 1 과 β값을 이용하여 계산하였으며, 14 16 그 결과를 2 에 함께 나타내었다. 그리고 α 1 의 변화에 따른 이들 활동도계수 값들의 변화를 . 1 에 함께 나타내었다. γ 1 의 값은 대각선(이상적 혼합미셀화)으로부터 큰 양의 벗어남을 보였으며, γ 2 의 값은 그래프가 대각선으로부터 큰 음의 벗어남을 나타내었다. 또한 γ 1 γ 2 의 값을 이용하면 단량체상으로 존재하는 각 성분들의 몰농도를 나타내는 C 1 과 C 2 의 값을 계산할 수 있다. 14 , 15 각 용액에 대하여 이와 같이 계산한 C 1 과 C 2 의 값을 2 에 함께 나타내었다. 여기서 알 수 있듯이 C 1 과 C 2 의 값은 γ 1 γ 2 의 값처럼 각각 이상적 혼합미셀화로부터 양의 벗어남과 음의 벗어남을 보였다.
한편 미셀 상을 열역학적으로 분리된 하나의 상으로 취급하고 순수 성분 상태를 미셀 상의 표준 상태로 정하면 혼합미셀 상에서 각 성분들의 활동도를 나타내는 α 1 M α 2 M 는 각 성분의 활동도 계수와 몰분율로부터 계산할 수 있으며, 16 계산한 결과를 2 에 함께 나타내었다. 미셀 상에서 CPC분자의 활동도를 나타내는 α 1 M α 1 에 비하여 훨씬 큰 양의 벗어남을 나타내었으며, Brij35의 활동도를 나타내는 α 2 M α 2 보다 훨씬 작은 음의 벗어남을 나타내었다. 계산한 α 2 M 의 값을 α 1 M 에 대하여 도시하였으며, 그 결과를 . 2 에 나타내었다. 여기서 그래프의 위쪽은 혼합 미셀을 형성하는 영역이며 그래프의 아래쪽 영역은 두 성분의 농도가 매우 묽기 때문에 계면활성제분자들이 미셀을 이루지 못한 상태의 영역을 나타낸다. 이 그래프에서 알 수 있듯이 CPC/Brij 35 혼합계면활성제의 미셀화에서는 미셀을 이루는 영역이 미셀을 이루지 않는 영역보다 훨씬 넓음을 알 수 있다.
PPT Slide
Lager Image
Plots of γ1(●) and γ2(▼) versus α1 for the micellization of CPC/Brij 35 mixed surfactant systems in water at 298K.
PPT Slide
Lager Image
Plots of α1M versus α1M for the micellization of CPC/Brij 35 mixed surfactant systems in water at 298K.
혼합미셀 상에서 각 성분들이 서로 혼합됨으로써 생기는 엔탈피 변화값(ΔH mix )은 비이상적혼합미셀 모델에 의하면 X 1 과 β값의 함수로서 주어지며, 17 계산한 ΔH mix 값들을 2 에 함께 나타내었다. 여기서 ΔH mix 값은 β값처럼 모두 음의 값이었으며, α 1 값이 0.2일 때 최소값을 나타내었다. 이것은 혼합미셀 상에서 두 종류의 계면활성제분자들 사이에 상호인력이 발생했음을 의미하며 또한 미셀 상에서 두 성분들이 서로 혼합됨으로써 열역학적으로 더욱 안정화됨을 의미한다.
Values of ΔGom(±0.05, kcal/mol), ΔHom(±0.05, kcal/mol), and ΔSom(±0.5, cal/mol K) for the micellization of CPC/Brij 35 mixed surfactant systems in water at various temperatures
PPT Slide
Lager Image
Values of ΔGom (±0.05, kcal/mol), ΔHom (±0.05, kcal/mol), and ΔSom (±0.5, cal/mol K) for the micellization of CPC/Brij 35 mixed surfactant systems in water at various temperatures
CPC/Brij 35 혼합계면활성제의 미셀화에 대하여 1 의 CMC와 B에 대한 데이터를 식 (2)에 대입하여 ΔG o m 값을 계산하였으며, 18 그 결과를 3 에 나타내었다. 각 온도에서 측정한 ΔG o m 의 값은 모두 큰 음의 값을 나타내었으며, 측정한 온도와 α 1 의 값에 따라 서로 다른 값을 나타내었다. 계면활성제분자들이 회합하여 미셀을 이룸으로써 계의 전체 자유에너지가 감소하게 되며, 그 결과 ΔG o m 의 값은 음의 값을 나타내게 된다. 이것은 계면활성제분자가 어떤 농도 이상에서는 미셀화가 자발적으로 형성하게 하는 이유라고 할 수 있다. ΔG o m 의 값은 ΔH o m 와 ΔS o m 값의 크기와 부호에 따라 결정되며, 이들 값의 크기와 기여도에 따라 미셀화 현상이 엔탈피지배반응 혹은 엔트로피 지배반응으로 나누어진다. 19
PPT Slide
Lager Image
CPC/Brij 35 혼합계면활성제의 미셀화에 대한 ΔH o m 와 ΔS o m 의 값은 각각 식 (3)과 (4)에 의하여 계산할 수 있으며, 계산한 결과를 3 에 함께 나타내었다. 순수 Brij 35를 제외한 순수 CPC와 모든 혼합계면활성제의 미셀화에 대하여 계산한 ΔS o m 값은 측정한 모든 영역에서 양의 값을 나타내었다. 또한 ΔH o m 값은 측정한 온도가 293K이하일 때에는 양의 값을 그리고 온도가 298K이상일 때에는 음의 값을 나타내었다. 그 결과 293K이하의 낮은 온도에서는 ΔG o m 값이 양의 값을 나타내는 ΔH o m 보다 양의 값인 ΔS o m 에 의하여 음의 값을 나타내게 되며, 298K이상의 높은 온도에서는 음의 값인 ΔH o m 와 양의 값인 ΔS o m 가 공동으로 기여함으로써 ΔG o m 의 값이 더욱 음의 값을 나타내게 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
3 에 있는 열역학 함수값(ΔG o m , ΔH o m 및 ΔS o m )과 온도와의 관계를 조사하기 위하여 α 1 이 0.4일 때 계산한 열역학 함수값들을 온도에 대하여 . 3 에 도시하였다. 여기서 보듯이 온도의 증가에 따라 모든 열역학 함수값들은 거의 직선적으로 정비례하면서 감소하는 경향을 보였다. 3 에 의하면 ΔH o m 와 ΔS o m 의 값은 측정한 모든 몰분율조성에서 온도의 증가에 따라 감소하였지만, ΔG o m 의 경우는 측정한 α 1 의 값에 따라 서로 다른 경향을 보였다. 즉, α 1 이 0.8 미만일 때 ΔG o m 의 값은 온도의 증가에 따라 감소하는 경향을 보였지만, α 1 이 0.8 이상일 때에는 오히려 약간 증가하는 경향을 보였다. 이것은 순수 CPC의 미셀화에서는 온도의 증가에 따라 엔탈피와 엔트로피의 감소가 작은데 반하여 순수 Brij 35의 미셀화에서는 온도변화에 따른 이들 열역학 함수값들이 크게 감소하기 때문에 나타나는 현상이다. Brij 35의 경우는 CPC와는 다르게 부피가 매우 큰 머리-그룹(polyoxyethylene)이 물 분자와 약하게 수소결합을 이룸으로써 친수성을 띠게 된다. 온도가 올라감에 따라 이 친수성 그룹과 결합하고 있는 물 분자들이 떨어져 나오기 때문에 Brij 35의 경우 온도의 변화에 따른 엔트로피와 엔탈피가 크게 변하게 된다.
PPT Slide
Lager Image
Plots of thermodynamic parameters versus temperature for the micellization of CPC/Brij 35 mixed surfactant system (α1=0.4) in water: △Hom(■); △Som(▲); △Gom(●).
한편 α 1 의 증가에 따른 열역학 함수값들의 변화를 조사하기 위하여 293K에서 측정한 열역학 함수값들을 α 1 에 대하여 도시하여 보았으며, 그 결과를 . 4 에 나타내었다. 여기서 알 수 있듯이 α 1 의 증가에 따라 ΔG o m 의 값은 증가하다가 감소하였으며, α 1 이 0.2일 때 최대값을 나타낸다. 그러나 ΔH o m 와 ΔS o m 의 값은 α 1 의 증가에 따라 감소하다가 증가하는 경향을 보였으며, α 1 의 값이 0.2일 때 최소값을 나타내었다. 3 에 의하면 α 1 에 따른 이러한 열역학 함수값들이 변하는 경향은 측정한 온도에 따라 약간 다르게 나타나고 있다. 즉, 298K 이하의 낮은 온도에서는 이들 열역학 함수값들이 모두 . 4 와 같은 경향을 나타내었지만, 303K 이상의 높은 온도에서는 . 4 와는 다른 경향을 보였다.
PPT Slide
Lager Image
Plots of thermodynamic parameters versus α1 for the micellization of CPC/Brij 35 mixed surfactant systems in water at 293K; △Hom(■); △Som(▲); △Gom(●).
CPC/Brij 35 혼합계면활성제의 미셀화에 대한 ΔH o m 와 ΔS o m 의 상관성을 조사하기 위하여 . 5 에 ΔH o m 을 ΔS o m 에 대하여 도시하였다. 그 결과, 모든 몰분율조성에서 직선을 이루었다. 따라서 ΔH o m 와 ΔS o m 사이에는 식 (5)와 같이 일차함수의 식으로 나타낼 수 있으며, ΔG o m 의 값은 식 (6)과 같이 나타낼 수 있다. 여기서 직선의 기울기( p )는 등구조온도(iso-structural temperature)로서 식 (6)에 의하면 ΔG o m 에 대한 ΔS o m 의 기여도가 0이 되는 온도를 나타낸다. 20 22 즉, 측정한 온도가 등구조온도( p )보다 높을 때에는 ΔS o m 에 의하여 미셀화가 더욱 촉진되지만, 반대로 등구조온도 보다 낮을 때에는 ΔS o m 에 의하여 미셀화는 더욱 억제된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Plots of △Hom versus △Som for the micellization of CPC/Brij 35 mixed surfactant systems in pure water: α1=0(●); α1=0.2(■); α1=0.4(▲); α1=0.6(○); α1=0.8(□); α1=1(△).
Derived least square parameters (from equation 5) and root mean square deviation (RMSD) for the micellization of CPC/Brij 35 mixed surfactant systems in water
PPT Slide
Lager Image
Derived least square parameters (from equation 5) and root mean square deviation (RMSD) for the micellization of CPC/Brij 35 mixed surfactant systems in water
. 5 에 나타낸 각 몰분율조성에 대한 직선의 기울기( p )와 절편( q )의 값을 최소자승법으로 계산하였으며, 그 결과를 4 에 나타내었다. . 5 4 에서 알 수 있듯이 각 몰분율조성에 대하여 계산한 p 의 값은 α 1 의 증가에 따라 감소하였으며, 특히 α 1 이 1(순수 CPC)일 때 최소값을 나타낸다. 이것은 순수 Brij 35의 경우보다 순수 CPC의 미셀화에서 엔트로피 효과가 더욱 크게 작용하고 있음을 나타낸다. 4 에 의하면 q 의 값은 △Gom값처럼 α 1 의 증가에 따라 변화하는 경향을 보였으며, α 1 이 1 일 때 최소값인 -8.65kcal/mol을 나타내었다.
결 론
CPC/Brij 35 혼합계면활성제의 미셀화에 대한 CMC와 B 값을 측정하여 여러 가지 열역학 함수 값을 계산하여 분석하였으며, 그 결과 혼합미셀상에서 Brij 35분자보다 CPC분자가 더욱 많이 존재하며 또한 여러 가지 열역학 함수값들이 비이상적 혼합미셀화 현상을 나타냄을 알 수 있었다. 즉, 계산한 X 1 의 값은 α 1 값과 비교하여 큰 양의 벗어남을 보였으며, α 1 M , γ 1 , C 1 의 값도 이상적 혼합미셀 모델로부터 크게 양의 벗어남을 보였다. 이에 반하여 α 2 M , γ 2 및 C 2 의 값은 큰 음의 벗어남을 보였다. 미셀 상에서 두 성분들 사이의 상호작용의 세기를 나타내는 β값은 모두 음의 값을 나타내었으며, 그 결과 ΔH mix 값도 음의 값을 나타내었다. 이것은 미셀 상에서 두 성분의 head-그룹들 사이에 이온-쌍극자 상호작용이 미치고 있음을 나타낸다.
한편 CPC/Brij 35 혼합계면활성제의 미셀화에 대한 CMC값은 온도의 증가에 따라 감소하다가 증가하는 경향을 보였다. 이러한 변화로부터 계산한 ΔG o m 값은 -7.5 kcal/mol에서 -11 kcal/mol까지 모두 음의 값을 나타내었다. 또한 ΔS o m 값은 순수 Brij 35의 경우를 제외한 모든 혼합계면활성제와 순수 CPC의 경우에 모두 양의 값을 나타내었다. ΔH o m 값은 -12 kcal/mol에서 5.7 kcal/mol까지 변하였으며, 측정한 온도에 따라 부호가 변화하였다. 즉, 293K 이하의 온도에서는 양의 값을 그리고 298K 이상의 온도에서는 모두 음의 값을 나타내었다. ΔH o m 을 ΔS o m 에 대하여 도시한 결과 직선이 얻어졌으며, 그 직선의 기울기( p )로 부터 구한 등구조온도는 244.9K에서 295.7K까지 변하였다. α 1 의 증가에 따라 p의 값은 감소하는 경향을 보였으며, α 1 이 1 일 때 최소값을 나타내었다.
References
Sierra M. B. , Morini M. A. , Schulz P. C. , Junguera E. , Aicart E. 2007 J. Phys. Chem. B 111 11692 -    DOI : 10.1021/jp072766s
Shang B. Z. , Wang Z. , Larson R. G. 2008 J. Phys. Chem. B 112 2888 -    DOI : 10.1021/jp0773841
Kim H. U. , Lee J. K. , Lim K. H. 2005 J. Korean Ind. Eng. Chem. 16 231 -
Peyre V. , Patil S. , Durand G. , Pucci B. 2007 Langmuir 23 11465 -    DOI : 10.1021/la701579e
Zhang G. , Chen X. , Zhao Y. , Xie Y. , Qiu H. 2007 J. Phys. Chem. B 111 11708 -    DOI : 10.1021/jp074945f
Deb N. , Shannigrahi M. , Bagchi S. 2008 J. Phys. Chem. B 112 2868 -    DOI : 10.1021/jp710299t
Garcia-Rio L. , Mendez M. , Paleo M. R. , Sardina F. J. 2007 J. Phys. Chem. B 111 12756 -    DOI : 10.1021/jp073510p
Imanishi K. , Einaga Y. 2007 J. Phys. Chem. B 111 62 -    DOI : 10.1021/jp065317l
Hildebrand A. , Garidel P. , Neubert R. , Blume A. 2004 Langmuir 20 320 -    DOI : 10.1021/la035526m
Lu S. , Somasundaran P. 2008 Langmuir 24 3874 -    DOI : 10.1021/la703233d
Kim H. U. , Lee J. K. , Lim K. H. 2005 J. Korean Ind. Eng. Chem. 16 231 -
Szymczyk K. , Janczuk B. 2007 Langmuir 23 8740 -    DOI : 10.1021/la7008495
Gil H. N. , Lee B. H. 2007 J. Korean Ind. Eng. Chem. 18 267 -
Park I. J. , Lee B. H. 2006 J. Kor. Chem. Soc. 50 190 -    DOI : 10.5012/jkcs.2006.50.3.190
Gil H. N. , Lee B. H. 2008 J. Kor. Chem. Soc. 52 111 -    DOI : 10.5012/jkcs.2008.52.2.111
Gil H. N. , Lee B. H. 2007 J. Kor. Chem. Soc 51 129 -    DOI : 10.5012/jkcs.2007.51.2.129
Holland P. M. , Rubingh D. N. 1983 J. Phys. Chem. 87 1984 -    DOI : 10.1021/j100234a030
Singh J. , Unlu Z. , Ranganathan R. , Griffiths P. 2008 J. Phys. Chem. B 112 3997 -    DOI : 10.1021/jp077380w
Lim K. H. , Kang K. H. , Lee M. J. 2006 J. Korean Ind. Eng. Chem. 17 625 -
Gil H. N. , Lee B. H. 2008 J. Kor. Chem. Soc. 52 461 -    DOI : 10.5012/jkcs.2008.52.5.461
Khatua D. , Ghosh S. , Dey J. , Ghosh G. , Aswal V. K. 2008 J. Phys. Chem. B 112 5374 -    DOI : 10.1021/jp7096636
Miyake M. , Einaga Y. 2007 J. Phys. Chem. B 111 535 -    DOI : 10.1021/jp0664465