[논문 리뷰] Hartree-Fock Model of a Self-Avoiding Flexible Polymer

C. Bouchiat|arXiv (Cornell University)|2002. 01. 24.

Electrostatics and Colloid Interactions참고 문헌 1인용 수 3

한 줄 요약

이 논문은 짧은 범위의 반발성 단위-단위 상호작용을 가진 연속적인 자유결합 사슬 모델을 사용하여, 자기방치성 유연 고분자의 힘-연장 거동을 모델링하기 위해 하트리-폭크(Hartree-Fock, HF) 장 이론 접근법을 개발한다. 주요 결과는 자가일관성 적분방정식에 대한 해석적 해를 구한 것으로, 이는 확장 곡선을 정확하게 예측하며, 반발 강도를 조정할 경우 N=100 단위에서 몬테카를로 시뮬레이션과 양호한 일치를 보인다. 이는 d ≤ b일 경우 힘-연장 형태가 잠재력 범위에 민감하지 않음을 시사한다.

ABSTRACT

Recent measurements of the force versus extension curves in stretched single stranded DNA, under conditions where the hydrogen bonding between complementary bases is inhibited, provide a new handle for the study of self-avoiding effects in a flexible polymer. We report in this paper upon analytic computations of the force versus extension curves within a continuous version of the freely joining chain model, with monomer-monomer repulsive interactions. The problem is formulated as a Statistical Field Theory model, endowed with a fixed cutoff associated with the curvature of the extension versus force curve, in the free polymer limit. Using a Field Theory version of the Hartree-Fock approximation, the self-avoiding single polymer problem reduces to the solving of a one-dimension self-consistent integral equation. Taking the short range potential limit, we obtain an explicit analytical solution, involving the roots of an algebraic equation. In the low force regime, we find that the slope of the relative extension at zero force increases steadily with the total monomer number N, while it stays constant for a free polymer. We compare, for N=100, our analytic results with those obtained by Monte Carlo simulations of a discrete freely joining chain, in presence of excluded volume effects. Despite a difference between the repulsive potential ranges, an appropriate choice of our single free parameter leads to a good agreement between the two computations.

연구 동기 및 목표

• 자유결합 사슬 모델의 연속적 버전을 제시하고, 단위-단위 반발을 포함한 연속적 장 이론 프레임워크를 개발하여, 수소결합이 억제되는 조건에서 단일사슬 DNA와 같은 유연 고분자에서 자기방치 효과를 모델링하는 데 목적이 있다.
• 자유결합 사슬 모델의 연속적 형태를 수립하고, 단위-단위 반발 상호작용을 포함하여, 탄성 고분자 사슬에서 배제 부피 상호작용을 다루는 데 목적이 있다.
• 다체 문제를 해결 가능한 일차원 자가일관성 적분방정식으로 줄이기 위해 하트리-폭크 근사를 적용한다.
• 해석적 HF 결과를 몬테카를로 시뮬레이션과 비교하여, 현실적인 고분자 거동을 예측할 수 있는 모델의 유효성을 검증한다.
• 반발 잠재력의 범위에 따른 변화에 대해 힘-연장 곡선 형태의 강건성(로버스트니)을 탐구한다.

제안 방법

• 자유결합 사슬 모델의 연속적 형태를 수립하고, 고분자 분할함수를 이미너리 시간에서의 양자역학 문제로 매핑한다.
• 영역 한계를 통한 짧은 범위의 반발 잠재력을 도입하여, 장 이론 수식에서 델타 함수 상호작용을 유도한다.
• 다체 상호작용을 분리하기 위해 하트리-폭크 근사를 적용하여, 일차원 자가일관성 적분방정식으로 문제를 축소한다.
• 짧은 범위 잠재력 근사에서 유도된 적분방정식을 해석적으로 풀어, 대수방정식의 근을 포함하는 명시적 해를 도출한다.
• 해결된 해를 사용하여 힘에 따른 상대적 연장도를 계산하고, 저력 범위를 초월하여 시뮬레이션과의 비교를 위해 외삽한다.
• N=100일 때 d=b 조건에서 몬테카를로 시뮬레이션 데이터와 해석 곡선을 일치시켜 단일 자유 매개변수(결합 상수 λ)를 校정한다.

실험 결과

연구 질문

• RQ1 자기방치 효과를 포함할 경우, 이상적인 자유결합 사슬에 비해 유연 고분자의 힘-연장 반응은 어떻게 변화하는가?
• RQ2 d ≤ b 조건을 만족할 경우, 반발 단위-단위 잠재력의 범위 d에 따라 힘-연장 곡선 형태가 얼마나 민감하지 않은가?
• RQ3 하트리-폭크 근사는 유연 고분자에서 배제 부피 효과를 정확하고 해석적으로 기술할 수 있는가?
• RQ4 해석적 HF 해가 동일한 시스템에서 배제 부피 제약 조건을 가진 몬테카를로 시뮬레이션 결과를 얼마나 잘 재현하는가?
• RQ5 이 모델을 더 현실적인 단위-단위 잠재력, 예를 들어 요쿠다 또는 지수형 형태로 확장하여 생물학적 고분자, 예를 들어 ssDNA를 더 잘 기술할 수 있는가?

주요 결과

• 자기방치 경우, 상대적 연장도는 총 단위 수 N이 증가함에 따라 증가하지만, 이상적 고분자 경우는 일정하게 유지된다.
• N=100일 때 λ=0.7인 해석적 하트리-폭크 해는 딱딱한 코어 반발(d=b)을 사용한 몬테카를로 시뮬레이션과 양호한 정량적 일치를 보인다.
• 힘-연장 곡선 형태는 잠재력의 범위 d에 대해 강건함을 발견하였으며, d ≤ b 조건을 만족할 경우 특정 짧은 범위 반발의 형태에 민감하지 않음을 시사한다.
• 모델은 영력 상태에서도 비영이 연장도를 예측하여, 낮은 힘에서도 시스템이 확장된 상태를 유지함을 나타내며, 실험 관측 결과와 일치한다.
• 해석적 해와 시뮬레이션 결과 간의 일치는, 서로 다른 근사(영역 한계 vs. 유한 범위 잠재력)를 사용하고 있음에도 불구하고, 물리적 거동이 잠재력의 범위보다는 배제 부피 상호작용의 강도에 의해 지배됨을 시사한다.
• 결과는 생물학적 고분자에 대한 하트리-폭크 방법의 추가 탐색을 장려하며, 비평면 보정 및 더 현실적인 상호작용 잠재력의 포함이 가능하다.