[고분자 물리] Polymers in dilute solution

jmt·2024년 6월 8일
0

고분자물리

목록 보기
11/13

Introduction

예전에 second virial coefficient로 A2A_2를 배웠었다. 이때 A2A_2가 증가하게 되면 good solvent를 의미하고 이는 solvent power를 증가시키기에 excluded volume도 마찬가지로 증가하게 된다. 그리고 Flory-Huggins Equation의 χ\chi값은 감소한다고 배웠다. 그런데 이런 경햐을 보이는 이유는 Flory-Huggins Theory가 통계역학을 기반으로 lattice model을 이용하여 계산을 유도했기 때문이다. 하지만 lattice model은 종에 따른 분자의 밀도가 동일하다고 보고, 분자가 동일한 크기의 lattice cell을 차지한다고 가정하게 된다. 하지만 dilute solution의 경우에서는 not homogeneous한 경우로 부분적으로 균일한 locally uniform 상태와 균일하지 않은 not uniform 상태로 lattice model을 적용시키기 어렵다. 이를 해결하기 위해 dilute solution에서는 Flory Krigbaum theory를 적용시킨다. 이는 mean-field approximation을 사용한다. lattice에 존재하게 되는 polymer molecules가 균일하게 분포되어 있다는 가정을 전제로 한다. 앞에서 이야기했듯 이 가정은 높은 농도(not dilute solution)에서만 성립하지만, 고분자 영역에서 원소 체적 내의 segment 밀도는 여전히 일정하다고 가정 할 수 있으며 이는 lattice model을 적용하게 한다.

chemical potential

우선 chemical potential을 고려해보자. chemical potential은 삼투압(osmotic pressure)로 측정할 수 있다고 재료물리화학 시간에 배웠다.

π=ρgh\pi = \rho gh

chemical potential의 정의는 ΔGM\Delta G^Mn1n_1에 대해 편미분한 값으로 다음과 같다.

ΔGMn1=ΔGMN1Navo=NavoΔGMN1\frac{\partial \Delta G^M}{\partial n_1} = \frac{\partial \Delta G^M}{\partial \frac{N_1}{N_{avo}}} = N_{avo} \frac{\partial \Delta G^M}{\partial N_1}

여기서 osmotic pressure로 chemical potential은 다음과 같이 구해짐을 이용하자.

μ1μ10=V1π\mu_1 - \mu_1^0 = - V_1 \pi

ΔGMn1\frac{\partial \Delta G^M}{\partial n_1}μ1μ10\mu_1 - \mu_1^0과 같기때문에, 다음이 성립한다.

NavoΔGMN1=V1πΔGMN1=V1πNavoN_{avo} \frac{\partial \Delta G^M}{\partial N_1} = -V_1 \pi \Rightarrow \frac{\partial \Delta G^M}{\partial N_1} = \frac{-V_1 \pi}{N_{avo}}

수평균 분자량 Mˉn\bar{M}_n을 구할 때 삼투압(π\pi)은 다음과 같이 정의됨을 배웠다.

πc2=RT(1Mn+A2c2+A3c22+)\frac{\pi}{c_2} = RT \left(\frac{1}{M_n} + A_2c_2 + A_3c_2^2 + \cdots \right)

이를 그대로 대입해주면, 다음과 같다. (RNavo=k\frac{R}{N_{avo}} = k임을 기억하자. kk는 볼츠만 상수이다.)

ΔGMN1=V1c2kT(1Mn+A2c2+A3c22+)\frac{\partial \Delta G^M}{\partial N_1} = -V_1c_2kT\left(\frac{1}{M_n} + A_2c_2 + A_3c_2^2 + \cdots \right)

이제 위 식의 좌변을 구해보자. 앞 시간에 배운 ΔGM\Delta G^M을 그대로 사용하여 이를 N1N_1에 대해 편미분해보자. 단, 여기서 주의해야 할 것은 다음과 같은 변수간의 상관관계이다.

N0=N1+rN2;ϕ1=N1N0;ϕ2=rN2N0.N_0 = N_1 + rN_2;\quad \phi_1 = \frac{N_1}{N_0};\quad \phi_2 = \frac{rN_2}{N_0}.

또한 지금 우리는 lattice model을 적용하고 있음을 기억하자!

ΔGM=kT(N1lnϕ1+N2lnϕ2+N1ϕ2χ)ΔGMN1=kT(lnϕ1+N1ϕ1N1ϕ1+N2ϕ2N1ϕ2+ϕ2χ+N1ϕ2N1χ)=kT(lnϕ1+N1ϕ2N0ϕ1+N2ϕ2N0ϕ2+ϕ2χN1ϕ2N0χ)=kT(lnϕ1+ϕ2N1N0ϕ1(=N1)N2N0rr+ϕ2χ(1N1N0))=kT(lnϕ1+ϕ21rϕ2+χϕ22)=kT(lnϕ1+ϕ2(11r+χϕ22))=kT(ln(1ϕ2)+ϕ2(11r)+χϕ22)\begin{aligned} \Delta G^M &= kT(N_1 \ln \phi_1 + N_2 \ln \phi_2 + N_1 \phi_2 \chi)\\ \frac{\partial \Delta G^M}{\partial N_1} &= kT \left(\ln \phi_1 + N_1 \frac{\frac{\partial \phi_1}{\partial N_1}}{\phi_1} + N_2 \frac{\frac{\partial \phi_2}{\partial N_1}}{\phi_2} + \phi_2 \chi + N_1 \frac{\partial \phi_2}{\partial N_1}\chi \right)\\ &= kT \left(\ln \phi_1 + N_1 \frac{\frac{\phi_2}{N_0}}{\phi_1} + N_2 \frac{-\frac{\phi_2}{N_0}}{\phi_2} + \phi_2 \chi - N_1 \frac{\phi_2}{N_0} \chi \right)\\ &= kT\left(\ln\phi_1 + \frac{\phi_2 N_1}{N_0 \phi_1 (=N_1)} - \frac{N_2}{N_0}\frac{r}{r} + \phi_2\chi\left(1- \frac{N_1}{N_0} \right) \right)\\ &= kT\left( \ln \phi_1 + \phi_2 - \frac{1}{r}\phi_2 + \chi \phi_2^2\right)\\ &= kT\left(\ln \phi_1 + \phi_2 \left(1- \frac{1}{r} + \chi\phi_2^2 \right) \right)\\ &= kT \left(\ln(1-\phi_2) + \phi_2\left(1-\frac{1}{r}\right) + \chi\phi_2^2 \right) \end{aligned}

A2A_2 and μ\mu

힘들게 구한 편미분의 결과는 chemical potential(μ\mu)를 의미한다. 반드시 기억해두자. 이제 μ\muA2A_2간의 관계를 알아보자. 결국 혼합용액의 깁스 자유 에너지를 N1N_1에 대해 편미분한 결과와 삼투압을 통해 유도한 비리얼 계수를 포함한 식은 같다. 그러므로 다음이 성립한다.

V1πNavo=ΔGMN1=kT(ln(1ϕ2)+ϕ2(11r)+χϕ22)=kT(ϕ2ϕ222ϕ233+ϕ2(11r)+χϕ22)=kT(ϕ2(1r)+(12χ)ϕ22+ϕ233+)\begin{aligned} \frac{-V_1 \pi}{N_{avo}} &= \frac{\partial \Delta G^M}{\partial N_1}\\ &=kT \left(\ln(1-\phi_2) + \phi_2\left(1-\frac{1}{r}\right) + \chi\phi_2^2 \right)\\ &=kT\left(-\phi_2 - \frac{\phi_2^2}{2} - \frac{\phi_2^3}{3} - \cdots + \phi_2 \left(1-\frac{1}{r} \right) + \chi \phi_2^2 \right)\\ &=-kT\left(\phi_2 \left(\frac{1}{r} \right) + \left(\frac{1}{2} - \chi \right)\phi_2^2 + \frac{\phi_2^3}{3} + \cdots \right) \end{aligned}

위 유도과정에서 아래와 같은 talyor expansion을 적용했다.

Talyor expansion ln(1x)=ϕ2ϕ222ϕ233\text{Talyor expansion } \Rightarrow \ln(1-x) = -\phi_2 - \frac{\phi_2^2}{2} - \frac{\phi_2^3}{3} - \cdots

이제 A2A_2와 연관지어야 하는데 V1ϕNavo-\frac{V_1 \phi}{N_{avo}}는 농도 c2c_2의 다항식으로 표현이 된다. 여기서 농도는 mass fraction이다. 하지만 우리가 혼합용액의 깁스 자유 에너지를 편미분한 다항식의 항에는 mass fraction이 아닌 volume fraction인 ϕ\phi값으로 표현된다. 그렇기에 volume fraction을 mass fraction으로 변환해줘야 한다. 이는 다음과 같은 과정을 따른다.

ϕ2=N2V2V=N2M2V×V2V=c2×v2\phi_2 = \frac{N_2 V_2}{V} = \frac{N_2M_2}{V} \times \frac{V_2}{V} = c_2 \times v_2

위 관계식을 앞에서 구한 식에 그대로 대입해주면,

V1πNavo=ΔGMN1=kT(ϕ2(1r)+(12χ)ϕ22+ϕ233+)=kT(c2v2(1r)+(12χ)(c2v2)2+(c2v2)33+)=V1c2kT(c2v2V1c21r+(12χ)(v2)2V1c2c22+(v2)33V1c2c23+)=V1c2kT(v2V11r+(12χ)(v2)2V1c2+(v2)33V1c22+)\begin{aligned} \frac{-V_1 \pi}{N_{avo}} &= \frac{\partial \Delta G^M}{\partial N_1}\\ &= -kT\left(\phi_2 \left(\frac{1}{r} \right) + \left(\frac{1}{2} - \chi \right)\phi_2^2 + \frac{\phi_2^3}{3} + \cdots \right)\\ &=-kT\left(c_2v_2 \left(\frac{1}{r} \right) + \left(\frac{1}{2} - \chi \right)(c_2 v_2)^2 + \frac{(c_2v_2)^3}{3} + \cdots \right)\\ &=-V_1c_2kT\left(\frac{c_2v_2}{V_1c_2}\frac{1}{r} + \left(\frac{1}{2} -\chi \right)\frac{(v_2)^2}{V_1c_2}c_2^2 + \frac{(v_2)^3}{3V_1 c_2}c_2^3 + \cdots \right)\\ &= -V_1c_2kT \left(\frac{v_2}{V_1}\frac{1}{r} + \left(\frac{1}{2} - \chi \right)\frac{(v_2)^2}{V_1}c_2 + \frac{(v_2)^3}{3V_1}c_2^2 + \cdots \right) \end{aligned}

여기서 다음을 이용하자.

r=V2V1=V2/M2V1/M2=v2V1/M2v2r=V1M2r = \frac{V_2}{V_1} = \frac{V_2 /M_2}{V_1/M_2} = \frac{v_2}{V_1/M_2} \rightarrow \frac{v_2}{r} = \frac{V_1}{M_2}

위 식을 대입해주면, 다음과 같다.

V1πNavo=ΔGMN1=V1c2kT(v2V11r+(12χ)(v2)2V1c2+(v2)33V1c22+)=V1c2kT(1M2+(12χ)(v2)2V1c2+(v2)33V1c22+)\begin{aligned} \frac{-V_1 \pi}{N_{avo}} &= \frac{\partial \Delta G^M}{\partial N_1}\\ &= -V_1c_2kT \left(\frac{v_2}{V_1}\frac{1}{r} + \left(\frac{1}{2} - \chi \right)\frac{(v_2)^2}{V_1}c_2 + \frac{(v_2)^3}{3V_1}c_2^2 + \cdots \right)\\ &= -V_1c_2kT \left(\frac{1}{M_2} + \left(\frac{1}{2} - \chi \right)\frac{(v_2)^2}{V_1}c_2 + \frac{(v_2)^3}{3V_1}c_2^2 + \cdots \right)\\ \end{aligned}

그러므로 다음이 성립한다.

V1πNavo=ΔGMN1V1c2kT(1Mn+A2c2+A3c22+)=V1c2kT(1M2+(12χ)(v2)2V1c2+(v2)33V1c22+)\begin{aligned} \frac{-V_1 \pi}{N_{avo}} &= \frac{\partial \Delta G^M}{\partial N_1}\\ -V_1c_2kT\left(\frac{1}{M_n} + A_2c_2 + A_3c_2^2 + \cdots \right)&=-V_1c_2kT \left(\frac{1}{M_2} + \left(\frac{1}{2} - \chi \right)\frac{(v_2)^2}{V_1}c_2 + \frac{(v_2)^3}{3V_1}c_2^2 + \cdots \right)\\ \end{aligned}

위 식의 결과로 알 수 있는 것은 χ\chiA2A_2의 관계이다. 괄호 안의 두번째항을 보면 다음과 같은 관계식을 알 수 있다.

A2=(12χ)(v2)2V1A_2 = \left(\frac{1}{2} - \chi \right)\frac{(v_2)^2}{V_1}

Excluded Volume

A2A_2는 2차 비리얼 계수로 excluded volume으로 근사할 수 있다. 그럼 다음과 같은 결론을 도출할 수 있다.

χ>1/2 poor solventA2<0 excluded volume χ=1/2 theta solventA2=0 excluded volume =0χ<1/2 good solventA2>0 excluded volume \begin{aligned} \chi > 1/2 \rightarrow \text{ poor solvent} &\Rightarrow A_2 < 0 \rightarrow \text{ excluded volume }\downarrow\\ \chi = 1/2 \rightarrow \text{ theta solvent} &\Rightarrow A_2 = 0 \rightarrow \text{ excluded volume }=0\\ \chi < 1/2 \rightarrow \text{ good solvent} &\Rightarrow A_2 > 0 \rightarrow \text{ excluded volume }\uparrow\\ \end{aligned}

여기서 또 내릴 수 있는 결론은 혼합 용액의 깁스 자유 에너지를 N1N_1에 대해 편미분한 결과식에서 entropic effect와 enthalphic effect를 의미하는 term으로 나눌 수 있다는 것이다. 이때 A2A_2χ\chi의 관계를 나타내는 관계식에서 (12χ)\left(\frac{1}{2} -\chi\right)이라는 계수를 주목하자. 12\frac{1}{2}은 entropic term에서 도출된 값이고 χ\chi는 enthalphic term에서 도출된 값임을 알 수 있다.

Modification of FH theroy : Temperature dependency

Excess Gibbs Free Energy

우선 과잉 깁스 자유 에너지(Excess Gibbs Free Energy)의 정의를 알아보자. 과잉 깁스 자유 에너지는 실제 혼합물의 열역학적 특성을 설명하는 중요한 개념이다. 이는 이상 혼합물과 실제 혼합물 간의 차이를 나타내며, 혼합물 내에서 발생하는 비이상적인 상호작용을 반영한다.

GE=GrealGidealG^E = G^{real} - G^{ideal}

그리고 앞에서 언급한 A2A_2항에 해당하는 RT(12χ)ϕ22-RT\left( \frac{1}{2}-\chi\right) \phi_2^2 해당하는 항을 ΔGE\Delta G^E라 정의하자. 여기서 Flory와 Krigbaum은 enthalpy parameter인 κ\kappa와 entropy parameter Ψ\Psi를 새롭게 정의하였다. 새롭게 정의한 parameter로 excess gibbs free energy인 ΔGE\Delta G^E를 다음과 같이 나타낼 수 있다.

ΔGE=ΔHETΔSEwhere ΔHE=RTκϕ22  ΔSE=RΨϕ22.=RT(κΨ)ϕ22\begin{aligned} \Delta G^E &= \Delta H^E - T\Delta S^E\\ \text{where }\Delta H^E &= RT\kappa \phi_2^2\; \quad \Delta S^E = R\Psi \phi_2^2.\\ &= -RT( \kappa - \Psi ) \phi_2^2 \end{aligned}

그러므로 위 경우에서 κ=Ψ\kappa = \Psi인 경우에 ΔGE=0\Delta G^E = 0이 되고, 결국 Greal=GidealG^{real} = G^{ideal}한 경우로, 실제 혼합 용액의 거동이 이상적인 경우와 동일하다고 볼 수 있다. 이는 theta stae, Θ\Theta-state라고 볼 수 있다.

Excluded volume effect at theta temperature

그리고 우리는 A2A_2항과 ΔGE\Delta G^E항이 같다고 보기로 했으므로,

RT(12χ)ϕ22=RT(κΨ)ϕ22-RT\left(\frac{1}{2}-\chi\right)\phi_2^2 = -RT(\kappa - \Psi)\phi_2^2

가 성립한다. 즉, (12χ)=κΨ(\frac{1}{2}-\chi) = \kappa - \Psi이다. 여기서 우리는 theta state, 즉 이상적인 거동을 하기 위해선 κ=Ψ\kappa = \Psi인 상태여야 한다. 이 개념을 이용해 다음과 같이 온도 κ\kappaΨ\Psi의 관계를 theta state일때의 온도를 θ\theta, 현재 polymer solution의 온도를 TT라고 했을 때 다음과 같이 정의할 수 있다.

θTκΨκ=θTΨ.\frac{\theta}{T} \equiv \frac{\kappa}{\Psi} \rightarrow \kappa = \frac{\theta}{T} \Psi.

그러면 우리는 excluded volume A2Ψ(1θT)\approx A_2 \approx \Psi\left(1-\frac{\theta}{T} \right) 로 근사할 수 있게 된다. TT에 대한 항으로 나타낸 이유는 temperature dependency를 알아보기 위함이다.

profile
고분자/컴공

0개의 댓글

관련 채용 정보