[논문 리뷰] Solidification of small para-H2 clusters at zero temperature

E. Sola, J. Boronat|arXiv (Cornell University)|2011. 03. 04.

Quantum, superfluid, helium dynamics참고 문헌 18인용 수 17

한 줄 요약

이 연구는 영온도에서 소형 페라-헤륨 클러스터의 기본상태 위상(액체 또는 고체)을 결정하기 위해 단계 제어된 시험 파동함수를 사용한 확산 몽테카를로(DMC) 방법을 사용한다. N = 55 이하에서는 균일한 상이 없이 번갈아가며 액체상과 고체상이 존재하는 비균일한 결정화가 나타나며, 이후에는 모든 클러스터가 고체상으로 전이된다. 이는 두 상 사이의 에너지 차이가 몇 킬빈 내외로 극히 작음에도 불구하고, 밀도 프로파일이 액체와 고체의 구조를 명확히 구분함을 보여준다.

ABSTRACT

We have determined the ground-state energies of para-H$_2$ clusters at zero temperature using the diffusion Monte Carlo method. The liquid or solid character of each cluster is investigated by restricting the phase through the use of proper importance sampling. Our results show inhomogeneous crystallization of clusters, with alternating behavior between liquid and solid phases up to N=55. From there on, all clusters are solid. The ground-state energies in the range N=13--75 are established and the stable phase of each cluster is determined. In spite of the small differences observed between the energy of liquid and solid clusters, the corresponding density profiles are significantly different, feature that can help to solve ambiguities in the determination of the specific phase of H$_2$ clusters.

연구 동기 및 목표

• 영온도에서 페라-헤륨 클러스터의 기본상태 위상(액체 또는 고체)을 결정하는 것.
• 에너지 및 구조적 차이를 비교하여 헤륨 클러스터의 위상 식별에 대한 모호함을 해결하는 것.
• 에너지 및 화학적 포텐셜 분석을 통해 마법의 클러스터 크기와 그 안정성 존재 여부를 조사하는 것.
• 중요도 샘플링을 통한 위상 제어가 DMC 시뮬레이션에서 정확한 위상 식별에 기여하는 바를 평가하는 것.
• N = 13에서 75까지의 클러스터에 대한 기본상태 에너지 및 구조적 성질을 확립하는 것.

제안 방법

• 보스성 페라-헤륨 클러스터의 다체 슈뢰딩거 방정식을 정확히 해결하기 위해 확산 몽테카를로(DMC) 방법을 사용하였다.
• 액체 또는 고체 구형으로 제한하기 위해 위상 특화 시험 파동함수와 중요도 샘플링을 사용하였다.
• 상호작용으로 실버라-골드먼 잠재력 사용 및 잘 깊이 약 -37 K.
• 일차 및 이차 상관 인자와 함께 가우시안 항을 포함한 제스트로 유형의 시험 파동함수 사용.
• 각 N에 대해 액체 및 고체 상에 대해 병렬 DMC 시뮬레이션 수행하여 에너지 및 구조적 성질 비교.
• 이阶 시간단계 전개 및 그린 함수 근사를 적용하여 통계적 수렴성 확보 및 시간단계 편향 감소.

실험 결과

연구 질문

• RQ1 N = 13에서 75까지의 페라-헤륨 클러스터가 영온도에서 기본상태 위상(액체 또는 고체)은 무엇인가?
• RQ2 클러스터 크기 N에 따라 액체상과 고체상 간의 에너지 차이가 어떻게 변화하는가?
• RQ3 소규모 에너지 차이에도 불구하고 밀도 프로파일이 액체와 고체 헤륨 클러스터를 명확히 구분할 수 있는가?
• RQ4 특정 클러스터 크기(마법의 수)가 액체 또는 고체 상에서 강화된 안정성을 보이는가?
• RQ5 화학적 포텐셜 µ(N)은 N 범위 전체에서 어떻게 변화하며, 이는 위상 번갈아가름을 어떻게 드러내는가?

주요 결과

• N ≤ 55일 경우, 페라-헤륨 클러스터는 균일한 상이 없이 액체상과 고체상이 번갈아가며 존재하며, 상에 대한 선호도가 없다.
• N = 55 이후로는 모든 클러스터가 고체상에서 에너적으로 안정되며, 완전한 결정화로의 전이를 나타낸다.
• 액체상과 고체상 간의 에너지 차이는 입자당 1 K 미만으로, 기본상태 에너지가 거의 degenerate함을 시사한다.
• 소규모 에너지 차이에도 불구하고 액체 클러스터와 고체 클러스터는 밀도 프로파일에서 현저한 차이를 보인다. 액체 클러스터는 중심 밀도 피크를 보이며, 고체 클러스터는 격자 구조로 인해 중심 빈공간을 나타낸다.
• 두 번째 에너지 차이 ∆2(N)는 N = 13, 19, 23, 25, 28, 30, 34, 37, 40, 43, 47, 49, 51, 53에서 명확한 마법의 클러스터 크기를 보이며, 이는 안정성 향상을 시사한다.
• 동일한 N에 대해 밀도 프로파일이 현저히 다르게 나타난다: 액체 클러스터는 강한 H2-H2 상호작용으로 인해 이중 껍질 구조를 보이며, 고체 클러스터는 중심 빈공간과 추가 껍질 구조를 유지하지만, N이 증가함에 따라 이 구조가 약화된다.