덴드리머 간 유효 상호작용의 본질: 전자화 수준에 따른 전기·반데르발스 조절

초록

본 연구는 전산분자동역학(MD)과 우산 샘플링(US) 기법을 이용해 폴리아미도아민(PAMAM) 덴드리머 두 개 사이의 평균 힘 퍼텐셜(PMF)을 계산하였다. 프로톤화 정도와 세대(크기)에 따라 상호작용이 순수한 반발에서 부분적인 인력으로 전환됨을 보였으며, PMF는 지수함수와 가우시안 함수의 합으로 잘 맞는다. 이는 기존 이론·코스게인 모델이 예측한 순수 가우시안 형태와는 차이가 있다.

상세 분석

이 논문은 전통적인 코스게인 모델이 간과해 온 원자 수준의 전하 분포와 수화층 효과를 전면에 반영한 전산 실험을 수행하였다. 우산 샘플링을 통해 중심‑대‑중심 거리(R)를 반응 좌표로 설정하고, 50~70개의 윈도우를 2 ns씩 평형화한 뒤 WHAM으로 PMF를 재구성하였다. 비프로톤화(고 pH) 조건에서는 R이 약 2 nm 이상일 때 약한 인력이 나타나며, 이는 말단 아민의 수소결합 및 반데르발스 상호작용이 주도함을 에너지 분해 분석으로 확인하였다. 반면 프로톤화(저 pH)에서는 전하가 전부 양전하로 변해 전기적 반발이 지배적이며, 전체 PMF가 전 구간에서 양(반발) 값을 유지한다.

PMF 형태를 정량화하기 위해 가우시안(두 개)과 지수·가우시안 혼합 모델을 비교하였다. 가우시안만으로는 강한 겹침 영역(R < R_g)에서 실험 데이터와 큰 차이를 보였으며, 특히 고세대(G4)에서 오차가 두드러졌다. 반면 지수함수(α₁·R/R_g)와 가우시안(α₂·R²/R_g²)의 합은 전체 거리 구간에서 높은 결정계수를 기록했다. 파라미터 ε₁이 ε₂보다 크게 나타난 점은 지수항이 주된 기여를 함을 의미한다. 이는 전하가 높은 경우 전기적 포텐셜이 거리와 거의 선형적으로 감소하는 반데르발스‑전기 복합 효과와 일맥상통한다.

구조적 분석에서는 두 덴드리머가 겹칠 때 가지가 얽히는 현상을 관찰했으며, 접촉 원자 수와 겹침 함수 O(r)를 통해 구형 입자 모델이 부적합함을 입증했다. 비구형성 지표 δ가 R이 감소함에 따라 증가하는데, 이는 전하와 스테릭 효과가 덴드리머 형태를 크게 변형시킨다는 증거이다. 또한 염 농도(0.1 M NaCl) 실험에서는 전기적 차폐가 강화돼 ε₁이 감소하고, 가우시안 항의 비중이 약간 늘어나는 경향을 보였다.

결과적으로, 이 연구는 (1) 프로톤화 수준에 따라 덴드리머 간 상호작용이 반발‑인력으로 전이될 수 있음을, (2) PMF가 순수 가우시안이 아니라 지수·가우시안 혼합 형태가 더 정확함을, (3) 고농도·고전하 조건에서 덴드리머의 형태 변형과 수화층 재배열이 상호작용에 중대한 영향을 미친다는 점을 명확히 했다. 이러한 발견은 코스게인 시뮬레이션에 새로운 파라미터화 지침을 제공하고, 나노복합재료 설계 시 덴드리머 간 거리와 전하 조절을 통한 구조 제어 전략에 직접적인 시사점을 준다.