Zero-Length Column에서 Intracrystalline Diffusivity를 추출하는 역문제

확산계수를 직접 볼 수 있다면 실험은 간단할 것입니다. 하지만 제올라이트 결정 내부에서 분자가 움직이는 모습을 ZLC 실험이 직접 관찰하는 것은 아닙니다.

실험에서 실제로 얻는 것은 퍼지 가스로 전환한 뒤 시간에 따라 감소하는 출구 농도입니다. 우리가 알고 싶은 것은 결정 내부의 확산계수인데, 측정되는 것은 장치 밖으로 빠져나온 분자가 만든 농도 응답인 셈입니다.

그래서 Zero-Length Column(ZLC)은 전형적인 역문제(inverse problem)입니다. 관측된 탈착곡선에서 출발해 그 곡선을 만들어낸 결정 내부의 확산시간과 확산계수를 거꾸로 추정해야 합니다.

‘Zero-Length’인데 정말 길이가 0인 것은 아니다

ZLC라는 이름만 보면 컬럼의 길이가 실제로 0인 것처럼 느껴집니다. 의미는 조금 다릅니다.

일반적인 크로마토그래피 컬럼처럼 긴 충전층을 사용하는 대신 매우 적은 양의 흡착제를 얇게 배치해 컬럼 내부의 축방향 농도분포나 긴 충전층 자체가 만드는 영향을 최소화합니다.

작은 흡착제 시료를 먼저 흡착질과 평형시킨 다음, 시점 t = 0에서 흡착질이 없는 퍼지 가스로 전환하고 출구 농도의 감소를 기록합니다. 매우 적은 양의 시료를 사용함으로써 긴 충전층에서 발생하는 복잡한 농도 전파보다 흡착제 내부의 수송 특성을 더 직접적으로 관찰하려는 방법입니다.

측정곡선은 보통 초기 농도 c₀로 정규화한 c/c₀의 시간 변화로 나타냅니다.

처음에는 흡착되어 있던 분자가 빠르게 방출되고, 시간이 지나면 신호가 점점 낮아집니다. 하지만 이 곡선의 감소 속도를 그대로 확산계수라고 읽을 수는 없습니다.

유량과 흡착평형, 결정 크기, 장치 내부의 기체 체류 등이 같은 응답에 함께 영향을 줄 수 있기 때문입니다.

탈착곡선에서 확산시간을 거꾸로 찾는다

결정 내부 확산이 지배적인 경우 ZLC 응답은 하나의 중요한 시간척도와 연결됩니다.

구형 결정의 반지름을 Rₚ, 결정내 확산계수를 D라고 하면 확산시간은 대략 다음과 같은 비율로 결정됩니다.

τ_D = Rₚ² / D

같은 크기의 결정이라면 D가 작을수록 분자가 결정 밖으로 빠져나오는 데 오래 걸리고 탈착곡선의 꼬리도 길어집니다. 반대로 D가 크면 결정 내부의 농도분포가 더 빠르게 완화됩니다.

여기서 중요한 것은 ZLC가 D를 직접 측정하는 장치가 아니라는 점입니다. 먼저 실험에서 c/c₀ 대 시간의 응답을 얻고, 이를 확산모델이 예측하는 곡선과 비교해 Rₚ²/D 또는 D/Rₚ²와 연결된 시간척도를 추정합니다.

따라서 결정 크기를 정확히 알지 못하면 D 역시 정확하게 결정하기 어렵습니다. 같은 탈착시간이라도 결정이 크기 때문에 오래 걸린 것인지, 결정 내부의 확산이 실제로 느리기 때문인지 구분해야 하기 때문입니다.

이 과정이 역문제인 이유도 여기에 있습니다. 결과인 탈착곡선은 하나지만 비슷한 곡선을 만들어낼 수 있는 원인은 여러 개일 수 있습니다.

같은 꼬리곡선이 같은 확산계수를 뜻하지 않을 수 있다

ZLC 해석에서 조심해야 할 상황은 느린 탈착을 발견하고 곧바로 ‘결정내 확산이 느리다’고 결론내리는 것입니다.

퍼지 유량이 충분하지 않으면 외부의 기체상과 흡착상 사이의 평형 과정이 전체 응답에 크게 반영될 수 있습니다. 장치 내부의 빈 공간이나 배관에 남은 기체도 낮은 농도 영역의 신호에 영향을 줄 수 있습니다. 초기 평형이 완전하지 않거나 결정 크기 분포가 넓어도 응답곡선의 모양이 변합니다.

표면 장벽(surface barrier)이 존재하는 경우는 더 까다롭습니다. 분자가 결정 내부에서는 충분히 빠르게 움직이더라도 결정 표면을 통과하는 단계가 느리면 전체 탈착응답은 느려집니다.

따라서 느린 ZLC 응답 하나만으로 작은 D를 확정하는 것은 안전하지 않습니다. 어떤 저항이 실제 시간척도를 만들고 있는지를 실험조건과 모델을 함께 이용해 구분해야 합니다.

유량을 바꾸는 실험이 중요한 이유

이 문제를 확인하는 강력한 방법 가운데 하나가 퍼지 유량을 바꿔 ZLC 실험을 반복하는 것입니다.

평형의 영향이 큰 조건과 결정 내부 확산이 지배하는 조건은 유량 변화에 같은 방식으로 반응하지 않습니다. 따라서 여러 유량에서 얻은 곡선을 함께 비교하면 한 조건의 탈착곡선 하나만 맞추는 것보다 지배 메커니즘을 구별하는 데 도움이 됩니다.

ZLC에서는 확산시간과 대류·흡착 과정의 상대적인 크기를 나타내는 무차원 변수 L이 중요한 역할을 합니다. 확산의 영향이 충분히 강한 영역에서는 장시간 응답에서 얻는 확산시간의 해석이 유량 변화에 상대적으로 덜 민감해질 수 있습니다. 반대로 평형 과정의 영향이 큰 조건에서는 응답곡선이 유량에 더 민감하게 변할 수 있습니다.

또 서로 다른 유량의 응답을 유량과 시간의 곱인 Ft를 이용해 비교하면 단순히 곡선 하나의 모양만 보는 것보다 질량수지와 지배 메커니즘을 점검하는 데 도움이 됩니다.

그래서 ZLC에서 유량 변화는 단순한 반복 측정이 아닙니다. 현재 얻은 확산계수가 정말 결정 내부 확산에서 나온 값인지 확인하기 위한 실험적인 검증 과정입니다.

잘 맞는 피팅이 항상 유일한 답은 아니다

역문제에서는 실험곡선과 모델이 잘 맞는다는 사실만으로 충분하지 않습니다.

예를 들어 D와 결정 크기의 불확실성이 동시에 존재하거나, 표면 장벽과 내부 확산이 비슷한 시간척도를 가지면 서로 다른 매개변수 조합이 비슷한 탈착곡선을 만들 수 있습니다. 흡착등온선의 비선형성 역시 응답곡선에 영향을 줄 수 있습니다.

초기와 중간시간 영역에는 많은 정보가 들어 있지만 초기 평형, 결정 크기 분포, 등온선 비선형성 같은 조건에도 민감합니다. 반대로 장시간 점근영역(long-time asymptote)은 확산시간을 추정하기에는 비교적 단순하지만 농도가 매우 낮아지면서 사용할 수 있는 신호 자체가 작아지는 문제가 있습니다.

따라서 가장 예쁘게 맞는 곡선 하나를 찾았다고 해서 확산계수가 유일하게 결정됐다고 볼 수는 없습니다. 유량을 바꿔도 비슷한 확산시간이 얻어지는지, 장치 자체의 응답은 충분히 빠른지, 표면 장벽을 내부 확산으로 잘못 읽고 있지는 않은지까지 확인해야 합니다.

ZLC가 까다로운 이유는 실험장치 자체가 복잡해서가 아닙니다. 간단해 보이는 탈착곡선 하나에서 직접 보이지 않는 결정 내부의 수송 특성을 거꾸로 찾아야 하기 때문입니다.

참고자료

Eic & Ruthven, Zeolites (1988)

Brandani & Ruthven, Adsorption (1996)

Brandani, Adsorption (2021)

Brandani, Chemical Engineering & Technology (2016)