Vacancy defects in square-triangle tilings and their implications for quasicrystals formed by square-shoulder particles

기술 요약: 정사각형-삼각형 타일링에서의 공공 결함과 정사각형-숄더 입자로 형성된 준결정의 함의

문제 정의
정사각형-삼각형 준결정(QC)은 연성 물질 시스템에서 널리 관찰되지만, 일반적으로 추가적인 육각형 타일로 나타나는 점 형태의 높은 밀도 결함을 포함하는 특성이 있다. 이러한 결함의 열역학적 역할은 여전히 명확히 이해되지 않았다. 즉, 높은 결함 농도가 자기 조립 과정 중의 동역학적 트래핑(kinetic trapping)의 결과인지, 아니면 준결정의 안정성에 기여하는 평형상의 고유한 특징인지가 불분분하다. 본 연구는 12회 대칭성을 가진 정사각형-삼각형 준결정의 안정성에 대한 이러한 결함의 열역학적 기여를 조사한다.

방법론
본 연구는 격자 기반 타일 모델과 미시적 입자 모델을 결합한 두 부분의 접근 방식을 채택한다:

1. 격자 기반 몬테카를로 시뮬레이션:

   • 모델: 저자들은 정사각형, 삼각형, 그리고 두 종류의 불규칙한 육각형 결함 타일인 "실드(shield)"와 "에그(egg)"(후자는 $Egg_L$과 $Egg_R$로 존재하는 카이랄 구조임)로 구성된 이상화된 타일 기반 모델을 활용한다. 이러한 결함은 완벽한 정사각형-삼각형 타일링에서 정점(vertex)이 제거될 때 자연스럽게 발생한다.
   • 알고리-즘: 저자들은 "이동 가능한 에지(mobile edges)"(경계에 있거나 결함 타일의 일부인 에지)를 식별하고 이를 $\pi/6$만큼 회전시킴으로써 구성 공간을 샘플링하는 몬테카를로(MC) 알고리즘을 사용한다. 이를 통해 구멍이나 겹침을 생성하지 않고도 국부적인 재배열 및 타일 유형 변화가 가능하다.
   • 앙상블: 시뮬레이션은 $N\gamma T$ 앙상블(고정된 정점 수 $N$, 선 장력 $\gamma$, 온도 $T$)에서 수행된다. 결함 농도를 제어하기 위해 각 결함 타일에 자유 에너지 페널티 $\epsilon_D$가 할당된다.
   • 엔트로피 계산: 열역학적 적분을 사용하여 결함이 있는 타일링과 결함이 없는 타일링 사이의 구성 엔트로피 차이($\Delta S$)를 계산한다. 연구는 엔트로피가 시스템 크기 및 결함 농도에 따라 어떻게 스케일링되는지 분석하며, 개별 결함의 기여와 조합적 혼합 엔트로피를 구분한다.

2. 미시적 입자 모델:

   • 시스템: 정사각형-숄더 포텐셜(하드 코어 $\sigma$, 높이 $\epsilon$ 및 폭 $\delta=1.4\sigma$의 연성 척력 숄더)로 상호작용하는 코어-코로나(core-corona) 입자 모델을 사용한다.
   • 자유 에너지 계산: 저자들은 에인슈타인 적분과 열역학적 적분을 사용하여 준결정 및 경쟁 주기적 상(정사각형 및 육각형 결정)의 진동 자유 에너지를 계산한다.
   • 결함 비용: 준결정에서의 공공(vacancy) 생성에 따른 진동 자유 에너지 비용을 결정하기 위해, 주기적 결정에 대한 방법론을 응용한다. 저자들은 단일 공공이 준결정 내에서 한 쌍의 독립적인 결함(실드 또는 에그)으로 완화된다는 점을 고려한다. 자유 에너지 계산 중 결함 확산을 방지하기 위해, 입자의 운동은 국부적인 타일링 기하학에 의해 정의된 "가상" 보로노이 유사 셀(Voronoi-like cells) 내로 제한된다.
   • 평형 농도: 저자들은 입자 모델로부터 얻은 진동 자유 에너지와 타일 모델로부터 얻은 구성 엔트로피를 결합하여 전체 깁스 자유 에너지를 구축한다. 평형 결함 농도는 이 자유 에너지를 최소화함으로써 도출된다.

주요 기여 및 결과

• 결함의 본질: 정사각형-삼각형 타일링에서 단일 입자(정점)를 제거하는 것은 단순한 공공을 만드는 것이 아니라, 하나의 실드와 하나의 에그(왼손잡이 또는 오른손잡이)로 이루어진 한 쌍의 육각형 타일을 생성하는 국부적 재배열을 유발한다.
• 구성 엔트로피 이득:
  • 결함의 도입은 상당한 구성 엔트로피 이득을 가져온다. 열역학적 극한에서, 결함 비용이 0($\beta\epsilon_D = 0$)인 시스템의 정점당 구성 엔트로피는 $0.554(1) k_B$로 발견되었으며, 이는 결함이 없는 정사각형-삼각형 타일링($0.120 k_B$)보다 약 5배 높다.
  • 결정적으로, 엔트로피 이득은 단순히 개별 결함 엔트로피의 합이 아니다. 개별 실드 및 에그 결함은 단순한 점 결함에 비해 음의 엔트로피 상수를 가지지만(이는 고립된 상태에서 이들이 덜 유리함을 시사함), 세 가지 결함 유형(Shield, $Egg_L$, $Egg_R$)의 조합적 혼합이 상당한 추가 엔트로피 항을 생성한다. 이 혼합 엔트로피는 총 엔트로피 이득을 순 양(+)의 값으로 높여 결함의 존재를 안정화시킨다.
• 평형 결함 농도:
  • 이러한 발견을 정사각형-숄더 입자 모델에 적용했을 때, 저자들은 주기적 결정에 비해 비정상적으로 높은 평형 결함 농도를 예측한다.
  • 저온($k_B T/\epsilon = 0.1$)에서는 결함 농도가 무시할 만한 수준($\sim 10^{-8}$)이다. 그러나 온도가 높아짐에 따라 진동 자유 에너지 비용은 감소하고 엔트로피 이득이 지배적이 된다.
  • $k_B T/\epsilon = 0.3$에서 총 결함 농도는 약 $1\%$ ($\sim 0.01$ per lattice site)에 도달하며, 이는 하드 스피어 결정에서 전형적인 공공 농도인 $\sim 10^{-4}$보다 현저히 높다.
• 상 안정성:
  • 구성 엔트로피를 포함하지 않을 경우, 정사각형-삼각형 준결정은 정사각형 및 육각형 결정의 공존 상태에 비해 열역학적으로 안정하지 않다.
  • 결함의 구성 엔트로피를 포함하면, 12회 대칭 준결정을 위한 안정 영역이 나타난다. 따라서 높은 결함 농도는 이 상의 안정성에 필요한 내재적인 열역학적 특징이다.
  • 결함의 존재는 공존 화학 퍼텐셜과 압력을 약간 변화시키지만, 준결정의 평형 밀도 범위를 약 7% 넓힌다.

의의
본 논문은 연성 물질 준결정에서 관찰되는 실드 및 에그 결함의 높은 빈도가 동역학적 트래핑이 아닌 평형 열역학의 결과라고 결론짓는다. 본 연구는 준결정에서의 결함이 건설적인 역할을 한다는 것을 보여준다: 개별 타일의 자유도와 다중 결함 유형의 혼합을 통해 접근 가능한 구성을 증가시킴으로써, 결함은 주기적 결정이 달성할 수 없는 상당한 엔트로피적 안정화를 제공한다. 결과적으로, 연성 물질 준결정을 정확하게 기술하기 위해서는 이러한 점 결함의 거동과 높은 농도를 반드시 고려해야 하며, 이 결함들은 이들 상의 역학에도 깊은 역할을 한다.