[논문 리뷰] Maximum Entropy and the Variational Method in Statistical Mechanics: an Application to Simple Fluids

Chih‐Yuan Tseng, Ariel Caticha|ArXiv.org|2004. 11. 25.

Statistical Mechanics and Entropy참고 문헌 28인용 수 3

한 줄 요약

이 논문은 통계역학에서 복잡한 확률분포를 최적의 근사로 생성하기 위해 최대 엔트로피(ME) 프레임워크를 제안한다. 특히 단순 유체에 대해 적용한다. 상대 엔트로피를 사용하여 경량화된 하드스피어 직경에 대한 가중 평균을 포함한 시험 분포의 집합에서 최적의 분포를 선택함으로써, 전통적인 변분 방법을 개선하고 레나드존스 유체에 대해 BH 및 WCA와 같은 최첨단 섭동 이론과 경쟁 가능한 결과를 도출한다.

ABSTRACT

We develop the method of Maximum Entropy (ME) as a technique to generate approximations to probability distributions. The central results consist in (a) justifying the use of relative entropy as the uniquely natural criterion to select a "best" approximation from within a family of trial distributions, and (b) to quantify the extent to which non-optimal trial distributions are ruled out. The Bogoliuvob variational method is shown to be included as a special case. As an illustration we apply our method to simple fluids. In a first use of the ME method the "exact" canonical distribution is approximated by that of a fluid of hard spheres and ME is used to select the optimal value of the hard-sphere diameter. A second, more refined application of the ME method approximates the "exact" distribution by a suitably weighed average over different hard-sphere diameters and leads to a considerable improvement in accounting for the soft-core nature of the interatomic potential. As a specific example, the radial distribution function and the equation of state for a Lennard-Jones fluid (Argon) are compared with results from molecular dynamics simulations.

연구 동기 및 목표

• 통계역학에서 비가역적인 확률분포 근사치를 생성하기 위한 일반적인 방법을 개발하는 것.
• 시험 분포의 집합에서 최적의 근사치를 선택하는 데 상대 엔트로피가 유일하게 자연스러운 기준임을 정당화하는 것.
• 단일 분포 근사치를 넘어서 시험 분포의 가중 평균을 포함시킴으로써 변분 방법을 확장하는 것.
• 클래식한 단순 유체, 특히 레나드존스 유체를 대상으로 ME 프레임워크를 시험하여 정확성과 강건성을 기준으로 평가하는 것.
• ME 기반 근사치가 BH 및 WCA와 같은 기존 섭동 이론에 비해 열역학적 및 구조적 성질 예측에서 경쟁력 있거나 이를 초월할 수 있음을 보여주는 것.

제안 방법

• 상대 엔트로피로 정보 손실을 정량화함으로써 새로운 제약 조건을 만족하면서도 사전 분포를 후행 분포로 업데이트하는 최대 엔트로피(ME) 원리를 사용한다.
• 다양한 직경을 가진 하드스피어 시험 분포의 집합을 사용하여 유체의 정확한 캐논리컬 분포를 근사한다.
• 다양한 하드스피어 직경에 대한 가중 분포를 ME 원리에 따라 최적화함으로써 소프트 코어 상호작용을 더 잘 기술하는 가중 평균 분포를 도출한다.
• ME 최적화 분포에서 반지름 분포 함수(RDF)와 상태 방정식을 유도하고 분자 동역학 시뮬레이션 및 섭동 이론과 비교한다.
• 단일 하드스피어 직경만을 사용할 경우 보골리우보 변분 원리가 ME 방법의 특수한 경우임을 밝힌다.
• 상대 엔트로피 분석과 비최적 분포가 얼마나 많이 배제되는지 분석함으로써 ME 방법이 표준 변분 방법보다 얼마나 향상되었는지 정량화한다.

실험 결과

연구 질문

• RQ1 최대 엔트로피 방법을 통계역학에서 다루기 어려운 확률분포 근사치를 체계적으로 생성하기 위해 어떻게 적용할 수 있는가?
• RQ2 왜 시험 분포 집합에서 최적의 근사치를 선택하는 데 상대 엔트로피가 유일하게 자연스러운 기준인가?
• RQ3 시험 분포의 가중 평균을 허용함으로써 단일 분포가 아닌 ME 방법이 표준 변분 방법보다 얼마나 향상되는가?
• RQ4 ME 방법이 레나드존스 유체의 반지름 분포 함수와 상태 방정식에 대해 정확한 예측을 낼 수 있는가? 기존 섭동 이론과 비교해보면 어떻게 되는가?
• RQ5 다양한 밀도와 온도에서 ME 접근법이 분자 동역학 시뮬레이션 및 BH 및 WCA와 같은 고급 이론과 비교해 유체 거동을 얼마나 잘 캐릭터라이즈하는가?

주요 결과

• ME 방법은 서로 다른 직경을 가진 하드스피어 분포에 대한 가중 평균인 확률분포를 성공적으로 생성하여 소프트 코어 원자 간 상호작용을 더 잘 기술한다.
• ME 방법에서 유도된 반지름 분포 함수(RDF)는 하드스피어 직경을 조정하지 않더라도 베르젯 시뮬레이션 결과의 기능 형태를 매우 잘 재현한다.
• ME 방법을 사용해 계산한 상태 방정식은 단순한 변분 방법보다 명확한 개선을 보이며, 특히 중간 밀도에서 분자 동역학 시뮬레이션과 양호한 일치를 보인다.
• 저밀도에서는 모든 방법(ME, 변분, BH, WCA)이 시뮬레이션과 일치하지만, 장거리 두 입자 상관관계 처리의 한계로 인해 차이가 발생한다.
• 고밀도에서는 WCA 섭동 이론이 가장 우수한 성능을 보이지만, ME 방법도 경쟁력 있고 광범위한 조건에서 뛰어난 성능을 보인다.
• 저온에서 ME 방법은 압력-밀도 이소스터머에 반데르발스 루프를 생성하여 액체-기화 전이의 정확한 정성적 행동을 보이며, 이는 시뮬레이션과 실험 결과와 일치한다.