[논문 리뷰] Theory of complex fluids in the warm-dense-matter regime, and application to an unusual phase-transitions in liquid carbon

M. W. C. Dharma‐wardana|arXiv (Cornell University)|2016. 07. 25.

High-pressure geophysics and materials인용 수 4

한 줄 요약

이 논문은 온난 밀도 물질(WDM) 영역에서 복잡한 유체를 모델링하기 위한 일급 원리적 접근으로 중성-가짜원자(NPA) + 초망원체체계(HNC) 방법을 제안하며, 액체 탄소에서 일시적인 공유 결합을 정확히 포착한다. 약 7 eV 및 1.0 g/cm³에서 액체-기화 및 금속-반도체 상전이가 동시에 발생하며, 평균 이온화도 $ Z $ 에서 급격한 불연속성이 나타나며, DFT+MD 및 PIMC 시뮬레이션과의 비교를 통해 검증된다.

ABSTRACT

Data from recent laser-shock experiments, density-functional theory (DFT) with molecular-dynamics (MD), and path-integral Monte Carlo (PIMC) simulations on carbon are compared with predictions from the neutral-pseudo-atom (NPA)+ hyper-netted-chain (HNC) approach for carbon, a complex liquid in the warm-dense matter regime. The NPA results are in good agreement, not only with high-density regimes that have been studies via PIMC, but even at low densities and low temperatures where transient covalent bonding dominates ionic correlations. Thus the `pre-peak' due to the C-C bond at $\\sim$1.4-1.6 \\AA$\\,$ and other features found in the pair-distribution function from DFT+MD simulations at 0.86 eV and 3.7 g/cm$^3$ etc., are recovered accurately in the NPA+HNC calculations. Such C-C bonding peaks have not been captured via average-atom ion-sphere (IS) models. Evidence for an unusual liquid $\ o$ vapor and metal$\ o$ semi-metal transition occurring simultaneously is presented. Here a strongly correlated metallic-liquid with transient C-C bonds, i.e., carbon at density $\\sim$ 1.0 g/cm$^3$ and mean ionization $Z=4$ transits abruptly to a disordered mono-atomic vapour at 7 eV, with $Z\\simeq$ 3. Other cases where $Z$ drops abruptly are also noted. The nature of $Z$, its discontinuities, and the role of exchange correlation, are reviewed. The limitations of IS models in capturing the physics of transient covalent bonding in warm dense matter are discussed.

연구 동기 및 목표

• 기존 모델이 실패하는 온난 밀도 물질(WDM) 영역에서 탄소와 같은 복잡한 액체를 정확하게 모델링할 수 있는 일급 원리적 방법을 개발하기 위해.
• WDM에서 일시적인 공유 결합과 전자 상관관계를 포착하지 못하는 평균 원자 이온-구역(IS) 모델의 한계를 해결하기 위해.
• 강한 상관관계를 가진 금속 액체에서 $ Z $, 구조 인자, 쌍임계함수, 전자 스펙트럼을 계산하기 위한 통합 이론적 프레임워크를 제공하기 위해.
• 탄소와 같은 저원자번호 원소에서 $ Z $ 의 불연속적 전이의 성격과 상 행동과의 연관성을 규명하기 위해.
• 넓은 밀도 및 온도 범위에서 고정밀 DFT+MD 및 PIMC 시뮬레이션과의 비교를 통해 NPA+HNC 접근법을 검증하기 위해.

제안 방법

• NPA 모델은 전자를 한 입자 밀도 $ n(r) $ 와 이온을 $ \rho(r) $ 로 처리하여, 효과적인 전자-이온 잠재력과 함께 코른-샤움 DFT 계산을 가능하게 한다.
• 초망원체계(HNC) 적분방정식을 사용하여 NPA에서 유도된 잠재력으로부터 쌍분포함수(PDFs)를 계산하며, 이온 상관관계를 포착한다.
• 자기일관성 있는 방법으로 프리델 합과 전하중성 조건을 통해 평균 이온화도 $ Z $ 를 결정하며, 열역학적 제약 조건과 일치함을 보장한다.
• 이 방법은 저온 영역(공유 결합이 지배하는 영역)에서부터 DFT+MD가 접근할 수 없는 고온 WDM 조건까지 원활하게 확장 가능하다.
• 유한 온도에서 코른-샤움-머민 형식을 사용하여 전자 구조를 계산하며, 점유수는 페르미-디랙 통계에서 유도된다.
• 모델은 결합된 전자와 자유 전자의 기여를 구분하며, $ Z = Z_n - n_b $ 로 표현되며, 여기서 $ Z_n = 6 $ 이고 $ n_b $ 는 결합 전자의 수이다.

실험 결과

연구 질문

• RQ1 NPA+HNC 방법은 넓은 밀도 및 온도 범위에서 액체 탄소의 쌍분포함수(PDFs)를 정확히 재현할 수 있는가, 특히 일시적인 C-C 공유 결합에 기인한 프리피크를 포함하여?
• RQ2 약 7 eV 및 1.0 g/cm³에서 관측된 액체-기화 및 금속-반도체 상전이의 성격은 무엇인가?
• RQ3 평균 원자 이온-구역(IS) 모델은 왜 WDM에서 공유 결합 특성을 포착하지 못하며, NPA 모델은 이러한 한계를 어떻게 극복하는가?
• RQ4 특히 약 7 eV 근처에서 관측되는 평균 이온화도 $ Z $ 의 불연속성은 전자 구조 변화와 탄소의 상 전이와 어떻게 관련되어 있는가?
• RQ5 NPA+HNC 방법은 탄소와 같은 경량 원소의 WDM 연구에서 계산 비용이 큰 DFT+MD 및 PIMC 시뮬레이션을 어느 정도 대체할 수 있는가?

주요 결과

• NPA+HNC 방법은 이온-이온 PDF에서 약 1.4–1.6 Å 부근에 위치한 C-C 결합 프리피크를 정확히 복원하며, 이는 IS 모델에는 존재하지 않지만 DFT+MD 및 PIMC 시뮬레이션에는 존재하는 특징이다.
• 0.86 eV 및 3.7 g/cm³ 조건에서 NPA+HNC PDF는 DFT+MD 결과와 높은 정밀도로 일치하며, 일시적인 공유 결합 특징까지 잘 반영한다.
• 약 7 eV 및 1.0 g/cm³에서 액체-기화 및 금속-반도체 상전이가 동시에 발생하며, $ Z $ 가 약 4에서 약 3으로 급격히 감소함으로써 전자 구조의 불연속적 변화가 나타난다.
• 평균 이온화도 $ Z $ 는 약 7 eV 근처에서 불연속성을 보이며, T = 7 eV 조건에서 두 개의 해를 가진다: 하나는 $ Z = 4 $ (1s는 완전 충만, 2s는 빈 상태)이고, 다른 하나는 $ Z = 2.78 $ 로, 강한 상관관계 상태를 나타낸다.
• NPA 모델은 $ Z_{cn} $, $ Z_b $ 와 프리델 합 간의 일致성을 보이며, 조건이 불연속적 전이 지점에서도 약 93% 이내로 일관되고 명확한 $ Z $ 값을 제공한다.
• 이 방법은 일시적인 C-C 결합이 존재하는 강한 상관관계 금속 액체에서의 전이와 고체가 아닌 단일 원자 기체 상태로의 전이를 성공적으로 포착하며, 압축성 및 전도도의 변화가 뚜렷하게 나타난다.