Molecular-Level Understanding of Polyethylene Glycol Adsorption on TiO2(001) at the Solid–Water Interfaces

기술 요약: 고체-물 계면에서 TiO₂(001) 상의 폴리에틸렌 글리콜 흡착에 대한 분자 수준의 이해

문제 정의
폴리에틸렌 글리콜(PEG) 개질은 TiO₂ 나노입자의 콜로이드 안정성, 수용성 및 계면 기능성을 향상시키기 위한 표준 전략이다. 그러나 산화물/액체 물 계면에서의 포화 PEG 피복도를 예측하는 것은 여전히 큰 과제로 남아 있다. 이러한 어려움은 PEG 흡착이 표면 결합, 사슬 간 패킹, 수화, 물 치환 및 형태적 엔트로피를 포함하는 복잡한 열역학적 균형에 의해 지배되기 때문에 발생한다. 기존 연구들이 아나타제 TiO₂(101) 면에 대한 PEG 흡착을 조사해 왔으나, 아나타제 TiO₂(001) 면의 뚜렷한 표면 원자 구조와 더 높은 반응성은 서로 다른 계면 거동을 초래할 수 있다. 수성 환경에서 포화 피복도를 합리적으로 제어하기 위해서는 TiO₂(001) 상의 PEG 흡착에 대한 분자 수준의 열역학적 이해가 필요하다.

방법론
본 연구는 수화된 아나타제 TiO₂(001) 상에 대한 PEG-5 (n=5) 흡착을 조사하기 위해 밀도 범함수 이론(DFT), 제일 원리 분자 동역학(AIMD), 가우시안 프로세스 회귀(GPR)를 결합한 계산 프레임워크를 채택하였다.

• 시뮬레이션 설정: 분산 상호작용을 고려하기 위해 vdW-DF2 범함수와 Ti 3d 오비탈에 대한 Hubbard U 보정(3.0 eV)을 사용한 VASP 패키지 기반의 스핀 편극 DFT 계산을 수행하였다. TiO₂(001) 표면은 세 개의 Ti–O 층을 가진 (2 × 2) 단위 셀로 모델링되었다.
• AIMD 프로토콜: 시뮬레이션은 298.15 K에서 액체 물 밀도(~1.00 g/cm³)를 갖는 정준 앙상블에서 수행되었다. PEG 분자는 표면에 수직인 브러시 형태의 구성으로 초기화되었다. 시스템을 평형화한 후, 10 ps의 AIMD 궤적을 실행하였다. 1 ps 간격으로 10개의 스냅샷을 추출하여 0 K에서 기하 구조를 최적화하고, 대표적인 에너지를 얻기 위해 평균을 산출하였다.
• 피복도 분석: PEG-5 흡착은 0.25, 0.50, 0.75, 1.00 단분자층(ML) 피복도에서 조사되었다.
• 자유 에너지 분해: 흡착 자유 에너지($\Delta G_{ads}$)는 다음 항목들로 분해되었다:
  1. 엔탈피 기여: 헤드 그룹(head-group) 흡착, 사슬 간 상호작나, PEG–물 용매화.
  2. 영점 에너지(ZPE).
  3. 엔트로피 기여: 엔트로피 손실이 면 특이적 결합보다는 형태적 제한에 의해 지배된다는 가정하에, 이전의 TiO₂(101) 연구로부터 채택된 1차 근사치.
  4. H₂O 탈착: AIMD에서 관찰된 분자적으로 흡착된 물 종(species)을 기반으로 계산됨.
• 통계 분석: 이산적인 AIMD 데이터로부터 연속적인 자유 에너지 프로파일을 보간하기 위해 radial basis function 커널을 가진 GPR이 사용되었다. 보간된 프로파일의 95% 신뢰 구간을 추정하기 위해 부트스트랩 리샘플링(1,000회 반복)이 적용되었다.

주요 결과

• 구조적 진화: 모든 피복도에서 PEG 사슬은 브러시 형태의 구성($R_F/D_m > 1.0$)을 취한다. 피복도가 0.25에서 1.00 ML로 증가함에 따라 사슬 간 중첩이 증가하고, 표면으로부터 물 분자가 점진적으로 치환된다. 0.75 및 1.00 ML에서는 조밀한 브러시 구조가 PEG 사슬에 대한 물의 접근성을 크게 감소시킨다.
• 엔탈피 기여:
  • 헤드 그룹 흡착: 분자당 정규화된 에너지는 0.25–0.50 ML에서 가장 유리하지만, 입체적 혼잡과 Ti 사이트에서의 국부적 전자 환경 변화로 인해 높은 피복도에서는 약화된다.
  • 사슬 간 상호작용: 인력적 상호작용은 피복도가 증가함에 따라 증가하며, 0.75–1.00 ML 근처에서 포화된다.
  • PEG–물 용매화: 이 항은 높은 피복도에서 점점 더 불리해진다. 총 PEG–물 수소 결합 수는 0.50 ML에서 정점에 도달한 후 0.75 및 1.00 ML에서 급격히 감소하며, 이는 조밀한 브러시로부터의 물 배제를 나타낸다.
• 열역학적 안정성:
  • 총 흡착 엔탈피($\Delta H_{PEG}$)는 모든 피복도에서 유리하게 유지되지만, 사슬 간 상호작용과 수화가 최적으로 균형을 이루는 중간 피복도(0.50 ML)에서 최대가 된다.
  • 총 흡착 자유 에너지($\Delta G_{ads}$)는 비단조적이다. 가장 유리한 값($\Delta G_{ads} = -7.85$ eV)은 0.50 ML에서 나타난다.
  • 0.75 및 1.00 ML에서 $\Delta G_{ads}$는 주로 형태적 제한 및 용매화 감소와 관련된 엔트로피 페널티로 인해 현저히 덜 유리해진다(1.00 ML에서 -1.30 eV까지 상승).
• 차등 흡착: 차등 흡착 자유 에너지($\Delta G_{diff}$)는 0.50 ML까지 음수이며, 이는 더 많은 PEG를 추가하는 것이 열역학적으로 유리함을 의미한다. 그러나 0.50 ML을 넘어서면 $\Delta G_{diff}$가 양수(0.75 ML에서 2.48 eV)가 되어, 추가적인 흡착이 불리함을 나타낸다.
• GPR 보간: GPR로 보간된 자유 에너지 프로파일은 0.51 ML 근처에서 최솟값(보간값 -8.11 eV)을 확인하였으며, 0.50 ML 데이터 포인트 주변의 좁은 95% 신뢰 구간을 통해 이산적인 AIMD 결과를 뒷받침한다.

의의 및 주장
본 논문은 수성 환경에서 TiO₂(001) 상의 PEG 포화에 대한 분자 수준의 열역학적 기초를 제공한다고 주장한다. 주요 발견은 포화된 PEG-5 피복도가 단순히 PEG–표면 결합 강도에 의해서만 결정되는 것이 아니라, 유리한 상호작용(표면 결합, 사슬 간 패킹, 수화)과 불리한 페널티(물 치환, 형태적 엔트로피 손실) 사이의 경쟁에 의해 결정된다는 것이다.

구체적으로, 본 연구는 수화된 TiO₂(001) 상에서 0.50 ML(표면 밀도 3.49 molecules/nm²에 해당)의 포화 PEG-5 피복도를 예측한다. 이 한계는 피복도를 넘어설 때 PEG–물 수소 결합의 상실, 용매화 감소, 그리고 조밀한 브러시 패킹과 관련된 엔트로피 페널티 증가로 인해 열역학적 구동력이 상쇄되기 때문에 도달된다. 이 연구는 PEG 피복도를 합리적으로 제어하려면 단순히 표면 결합 부위를 최대화하는 것이 아니라, 이러한 상충하는 엔탈피 및 엔트로피 요인들 사이의 균형을 맞추어야 함을 보여준다.