Polymer Models
Overview
세포 핵 내 크로마틴/DNA의 구조와 동역학을 기술하기 위해 사용되는 coarse-grained 폴리머 모델들을 정리한 노트다. 가장 단순한 Gaussian (Ideal) Chain에서 시작해, excluded volume을 고려한 Self-Avoiding Walk, 동역학을 부여한 Rouse Model, 그리고 그 일반화 버전들(Freely-Joined-Chain, β-polymer, Worm-Like Chain, Rod-like)까지 다룬다. 모든 모델은 결국 같은 질문으로 수렴한다 — monomer 간 상호작용(또는 그 부재)이 폴리머의 정적 통계량(
내용은 Amitai & Holcman, Physics Reports 678 (2017) 리뷰의 Section 3을 기반으로 하며, 📑 Journal reading - Polymer physics of nuclear organization and function에서 분리되어 독립 개념 노트로 정리되었다.
Key Points
공통 기호
| 기호 | 의미 |
|---|---|
| 모노머 수 | |
| Kuhn length | |
| 스프링 상수 | |
| 확산계수 | |
| 이상확산 지수 | |
| end-to-end 벡터 | |
| 회전반경 | |
| Flory exponent | |
| dimension of space |
기본 폴리머 (Gaussian Chain, Ideal Chain)
기본 아이디어: 폴리머를 일정한 통계적 규칙을 따르는 random walk trajectory로 본다. 즉 한 monomer에서 다음 monomer로의 변위
이때 self-avoidance를 무시한다 (monomer끼리 통과 가능). 그래서 Ideal Chain 혹은 Phantom Chain이라고도 부른다.
왜 monomer 간 변위를 Gaussian으로 두는가?
여기서 ‘monomer’는 실제 화학적 단위가 아니라, 여러 개의 미시적 결합을 묶어 만든 coarse-grained segment (Kuhn segment) 다. 실제 폴리머(예: dsDNA)에서 인접한 화학 결합들은 방향이 서로 강하게 상관되어 있다 (backbone이 뻣뻣함). 그러나 persistence length보다 충분히 긴 스케일로 묶으면, 그 너머의 방향 상관은 사라진다.
하나의 Kuhn segment 변위(step)
이 미시 벡터들은 (거의) 독립이고 동일 분포를 따르므로, 중심극한정리(Central Limit Theorem)에 의해 그 합
- Kuhn length
: 방향 상관이 소멸하는 길이 스케일. 이 단위로 coarse-grain하면 각 step이 독립이 된다. - 결과적으로 Gaussian chain은 각 step이 독립인 random walk이며, 전체 통계도 Gaussian이 된다.
정적 통계량
End-to-end 벡터:
Step들이 독립이고 등방적이라(
여기서
모든 monomer 질량을 1로 두면, center of mass:
Radius of gyration
에너지: Entropic Spring
평형 상태에서 어떤 configuration이 존재할 확률은 Boltzmann 분포
이 “스프링”은 실제 화학 결합의 강성이 아니라, 각 step 변위가 Gaussian이라는 사실의 다른 표현이다. Gaussian 분포
스프링 상수는
두 monomer
요약: Gaussian chain은 폴리머를 독립 step의 random walk로 보는 정적(equilibrium) 구조 모델이다. CLT가 step과 전체 통계 모두를 Gaussian으로 만들고, 그 결과 entropic spring 에너지와
스케일링이 따라온다. (여기에 시간축을 얹어 동역학을 부여한 것이 Rouse Model이다.)
Self-Avoiding Random Walk Model (SAW)
Gaussian chain은 random walk의 궤적이기 때문에, walker가 한번 지나왔던 장소를 다시 지나가는 것 처럼, monomer의 위치가 겹칠 수도 있는 (혹은 아주 가까워질 수 있는) 모델이다.
실제 polymer는 그렇지 않다. momomer가 부피를 가지기 때문이다.
이러한 점을 고려한 게 self-avoiding random walk이다.
Excluded volume interaction 에너지를 추가한다. 디락델타를 써서 monomer 위치가 겹치는 순간 에너지가 무한해진다. 같은 위치에 있는 걸 절대 금지시키는 것이다.
이 효과 때문에 동일 조건의 gaussian chain보다 polymer가 부풀어 오른다.
AI 보충 (논문 외 내용) — 왜 excluded volume 포텐셜이 디랙델타로 주어지는가?
아래 내용은 원 논문에 없는, Glia(AI)가 보충한 설명이다.
의문: 진짜 반발은
인 모든 순간 무한대여야 하는데, 왜 식 (9)는 한 점에서만 터지는 디랙델타일까? 그리고 왜 가 앞에 곱해져 있을까? 핵심: 식 (9)는 진짜 hard-core 포텐셜이 아니라, 그것을 큰 스케일에서 coarse-grain한 유효(effective) 표현이다.
1. 진짜 포텐셜과 Boltzmann factor
진짜 hard-core 반발은
에서 , 그 밖에서 . 그러나 통계역학에서 실제로 중요한 건 자체가 아니라 Boltzmann factor 다. 2. Mayer function
두 monomer가 겹치는 작은 구 영역에서만
인 함수. 3.
는 어디서 오는가 — 적분하면 부피 상호작용의 통계적 효과는
을 공간에 적분한 양으로 측정되며, 이게 바로 excluded volume의 정의다:
즉
는 “유한 반지름 안에서 무한대”라는 hard-core 정보가 적분을 통해 부피 하나로 압축된 결과다. 반지름과 “무한대”는 사라지고, 그 영향력의 총량만 부피 로 남는다. 4. 그래서 디랙델타가 등장한다 (coarse-graining)
폴리머 통계는 monomer 크기보다 훨씬 큰 스케일(
)에서 일어난다. 이 큰 스케일에서 폭이 인 좁은 반발 영역은 사실상 점처럼 보인다. 따라서 같은 적분값을 갖는 디랙델타로 대체한다:
모든 쌍에 대해 합하고 double-counting 보정
을 붙이면 정확히 식 (9)가 된다. 5. 그럼 “겹치면 무한대”는 어디 갔나?
디랙델타 포텐셜은 한 점에서 무한대지만 적분 가능해서 겹침을 완전히 금지하진 않는다. 대신 겹침 configuration의 통계적 가중치를 정확히
만큼 깎는다. 큰 스케일 통계량( 등)을 계산할 때는 진짜 hard-core와 디랙델타 버전이 같은 답을 주기 때문에 이 교체가 정당화된다.
End to end distance가
Flory exponent는 차원으로 결정된다.
: :
AI 보충 (논문 외 내용) — Flory exponent의 유도와 0.588의 정체
아래 내용은 원 논문에 없는, Glia(AI)가 보충한 설명이다.
1. Flory의 mean-field 유도
폴리머 크기
에 대해 경쟁하는 두 자유에너지 항을 균형시킨다. ① Excluded volume 반발 (swelling): 부피
안에 개 monomer가 있으면 농도 . 충돌 쌍 수 에 excluded volume 를 곱하면
② Entropic spring 복원 (shrinking): ideal 크기보다 늘이면 엔트로피 손실이 생긴다.
③ 균형:
을 에 대해 최소화 ( ):
따라서 Flory 지수는
차원 비고 1 1 정확 (완전히 펴짐) 2 3/4 = 0.75 정확 3 3/5 = 0.6 거의 정확 4 1/2 정확 (ideal chain으로 환원) 2. 그럼 왜 0.6이 아니라 0.588인가?
Flory 공식은 mean-field 근사라 정확하지 않다. 반발 에너지 과대평가와 Gaussian 단순화의 오차가
에서 우연히 거의 상쇄되어 이라는 좋은 값을 줄 뿐이다. 정확한 값은 renormalization group (RG) 이론에서 나온다. self-avoiding walk를
극한의 field theory로 매핑하고 전개를 하면
고차항 + 수치적 재합(resummation), conformal bootstrap, Monte Carlo를 종합하면
이것이 여러 논문에서 보이는 0.588의 정체다.
3. 왜 이 숫자를 여기저기서 마주치는가 — Universality
self-avoiding walk는
인 universality class에 속한다. 그래서 같은 임계지수를 공유하는 모든 시스템(자성 물질의 임계현상, 이론 등)이 동일한 을 보인다. polymer physics에서 RG가 빛을 발하는 가장 깔끔한 입구. 요약:
= Flory mean-field, . / = RG · -전개로 얻은 정확한 임계지수 ( universality class).
Rouse Model
독립된 노트로 작성: Rouse Model
특성: 인접 모노머만 스프링으로 연결된 가장 기본적인 동적 폴리머 모델. Gaussian chain의 에너지(식 5)에 Langevin 동역학을 얹은 것. 스프링 상수는
동역학의 출발점: Langevin 방정식
질량
폴리머 동역학에서는 관성항이 마찰항에 비해 무시되는 overdamped(Smoluchowski) 극한
여기서 확산계수
모노머 동역학
Rouse 포텐셜(식 5)에서
각 항의
Rouse 모드 (정규 좌표)
식 (14)는 모노머들이 서로 결합되어 있어 직접 풀기 어렵다. Rouse 변환으로 모드를 분리(decouple)한다:
변환 계수 (식 18):
(식 17)은 descrete Fourier transform이며
고유값 (식 20):
이제 식 (14)가 분리된
여기서
단일 모노머의 MSD와 이상확산
모노머
각 OU 모드의 분산 (식 62):
완화 시간 (식 63):
이를 합치면 모노머 MSD의 완전한 표현 (식 64):
여기서 세 가지 체제가 나온다:
- 짧은 시간 (
): 지수항을 Taylor 전개하면 (p 무관). 합하면 (식 65):
- 중간 시간 (
): 합을 적분으로 근사 (식 66):
모든 내부 완화 시간이 기여하여 이상확산
- 긴 시간 (
): 지수항이 모두 포화되고 무게중심 확산만 남아 정상 확산, 단 확산계수는 .
핵심: Rouse 모델의 이상확산(
)은 외부 매질의 점탄성이 아니라, 내부 모드들의 완화 시간이 넓게 퍼져 있다는 사실에서 순수하게 기하학적으로 나온다. 이것이 이상확산에서 다룰 다른 메커니즘(FBM 등)과 구별되는 점이다. 이 섹션의 전체 유도 과정(friction의 미시적 기원, normal mode 분리, stress tensor)은 Rouse Model 노트에 자세히 정리되어 있다.
Freely-Joined-Chain Model (Generalized Rouse)
기본 아이디어: Rouse 모델은 인접 monomer 사이 스프링 포텐셜의 최소점이
에너지 (식 16):
- 정상상태에서 인접 monomer 간 평균 거리는
. 으로 두면 정확히 고전적 Rouse 모델로 환원된다.
왜
AI 보충 (논문 외 내용) — "FJC = ideal chain"과의 충돌 정리
아래 내용은 원 논문에 없는, Glia(AI)가 보충한 설명이다.
위키피디아 등에서는 FJC를 ideal chain과 같다고 소개하는데, 이 논문은 FJC를 generalized Rouse라고 한다. 둘 다 맞다 — 이름이 서로 다른 대상을 가리키기 때문이다.
- 정통 FJC (위키): 결합 길이를
로 엄격히 고정하고 각도만 무작위. 큰 에서 CLT에 의해 end-to-end 분포가 Gaussian으로 수렴 → 큰 스케일에서 ideal chain과 같은 통계. “FJC = ideal chain”은 이 의미. - 논문의 FJC (식 16): 결합 길이를
주위로 스프링으로 묶어 요동을 허용한 후 Langevin 동역학을 얹은 버전. Rouse가 인 특수 케이스이므로 저자 입장에서는 “최소점을 옮긴 generalized Rouse”.
β-폴리머 (Generalized Rouse)
독립된 노트: β-polymer (beta-polymer)
목적: 이상확산 지수
Rouse 고유값을 수정:
중간 시간대 이상확산 지수:
따라서:
(Rouse): : ← 실험값과 부합 : (완전 응축)
Worm-Like Chain (WLC)
특성: Semi-flexible polymer. Rouse/Gaussian chain은 구부림 강성이 없어 완전히 유연하지만, 실제 dsDNA·크로마틴 파이버는 일정 길이까지 뮣뮣하다. WLC는 이 stiffness(구부림 강성) 를 도입한 모델이다.
접선 벡터 상관함수
WLC의 정의적 특성은 접선 벡터
이산 모델에서 단위 접선 벡터로 평균을 내면 같은 지수 감쇄가 나온다 (식 37):
.
구부림 에너지 (Bending Energy)
stiffness는 인접한 결합 벡터 사이의 각도에 에너지 비용을 부과하여 구현한다. 단위 결합 벡터
두 형태가 같은 이유:
지속 길이와 구부림 강성의 관계
구부림 강성
동역학 시뮬레이션
WLC의 시간 발전은 총 포텐셜에 대한 overdamped Langevin 방정식으로 시뮬레이션한다 (식 38, 39):
여기서
에너지 스케일 비교
논문의 coarse-grained 크로마틴 시뮬레이션(
- 탄성 에너지:
N·m (DNA 최대 신장 ~10%, nm) - 구부림 에너지:
N·m ( 일 때)
비교하면, 같은 조건에서 순수 Rouse 폴리머는 구부림 항이 없어 구면 경계 근처에서도 구부러지지 않고 더 응축된 (
핵심: WLC는 Rouse/FJC에 구부림 강성
를 더한 것으로, 지속 길이 스케일 이하에서는 뻣뻣한 막대처럼, 그 이상에서는 유연한 사슬처럼 거동한다. 이 stiffness는 뒤의 NETP(좁은 탈출 시간)에서 탐색 효율에 직접 영향을 준다 (뻣뻣한 사슬이 경계 근처 타깃을 더 빨리 찾음).
Rod-like 모델
각도
장거리 모노머 상관을 가지며, Rouse와 다른 루핑 공식 유도 가능.
Questions & Insights
- 현재까지 기록된 질문 없음 — 학습 중 궁금한 점이 생기면 이곳에 추가.
Related Concepts
- Central Limit Theorem
- Entropic Spring Constant
- MSD exponent from scale-dependent Flory exponent
- Contact probability exponent and polymer scaling
- 이상확산 (Anomalous Diffusion)
References
- Amitai, A. & Holcman, D., “Polymer physics of nuclear organization and function,” Physics Reports 678, 1–83 (2017) — Section 3
- 📑 Journal reading - Polymer physics of nuclear organization and function
Notes from Claude
이 노트는 원래 📑 Journal reading - Polymer physics of nuclear organization and function의 Section 3 발췌본으로 시작됐던 것을, 독립적인 개념 노트 형식(frontmatter + Overview/Key Points/Related Concepts/References 구조)으로 재구성한 것이다.
- 본문에 있던 Glia(AI)의 보충 설명 콜아웃(
[!note] AI 보충)은 해당 수식 직후 원래 위치에 그대로 유지했다 — excluded volume 포텐셜이 디랙델타로 근사되는 이유, Flory exponent의 mean-field 유도와 0.588의 정체, FJC와 ideal chain 명명 충돌 정리 세 가지다. 내용 손실 없이 그대로 보존했다. - 원문에서
[[#3-이상확산|3. 이상확산...]]처럼 같은 파일 내부 앵커로 걸려 있던 링크는 실제로는 모(母) 저널 노트의 섹션을 가리키고 있었기 때문에, 정확한 cross-note 링크로 수정했다. - Related Concepts에 사용자의 현재 핵심 연구 주제와 직결되는 MSD exponent from scale-dependent Flory exponent를 추가했다 — 이 노트의 Rouse/β-polymer 이상확산 지수가 바로 그 연구의 이론적 출발점이기 때문이다.