2025-11-13T11:31:10.165583

Localized fermions on the triangular lattice with Ising-like interactions

RegeciovÃ¡, Kapcia

The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic $J$ interactions and onsite Coulomb $U$ interactions. The model is investigated for both signs of $J$, arbitrary $U$ interaction and arbitrary chemical potential $Î¼$ (or, equivalently, arbitrary particle concentration $n$). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case ($J<0$), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to $k_{B}T/|J| \approx 0.61$) is found to be second-order (for sufficiently large temperatures $k_{B}T/|J| \gtrsim 0.2$) or first-order (for $-1<U/|J|<-0.65$ at the half-filling, i.e., $n=1$). In the case of $J>0$, the ordered phase occurs in a range of $-1/2<U/|J|<0$ (for $n=1$), while for larger $U$ the state with short-range order is also found (also for $n \neq 1$). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for $Î¼\neq 0$ ($n \neq 1$) and $U/|J|>-1/2$ is always second-order for any model parameters. The ordered phase for $J>0$ can be stable up to $k_{B}T/|J| \approx 0.06$.

academic

삼각 격자 위의 국소화된 페르미온과 이징형 상호작용

기본 정보

논문 ID: 2412.10808
제목: Localized fermions on the triangular lattice with Ising-like interactions
저자: Lubomíra Regeciová (슬로바키아 과학 아카데미), Konrad Jerzy Kapcia (아담 미츠키에비치 대학교)
분류: cond-mat.stat-mech cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.comp-ph
발표 시간: 2025년 10월 17일
논문 링크: https://arxiv.org/abs/2412.10808

초록

본 연구는 몬테카를로 시뮬레이션을 통해 정준 앙상블에서 삼각 격자 위의 국소화된 페르미온 모델을 분석했습니다. 시스템의 해밀토니안은 확장된 허바드 모델(원자 극한)의 형태를 가지며, 최근접 이징형 자기 상호작용 J와 온사이트 쿨롱 상호작용 U를 포함합니다. 연구는 J의 두 가지 부호, 임의의 U 상호작용, 그리고 임의의 화학 포텐셜 μ(또는 동등하게 임의의 입자 농도 n)를 조사했습니다. 비열과 부분 격자 자화 강도 분석을 바탕으로 모델의 상도를 결정했습니다. 강자성 경우(J<0)에서, 정렬 상에서의 전이는 충분히 높은 온도(k_BT/|J|≳0.2)에서는 2차 상전이이거나 반채움 상태에서(-1<U/|J|<-0.65) 1차 상전이일 수 있습니다. J>0인 경우, 정렬 상은 -1/2<U/|J|<0 범위 내에서 나타나며(n=1일 때), 더 큰 U 값에 대해서는 단거리 정렬 상태도 발견되었습니다.

연구 배경 및 동기

핵심 문제: 기하학적으로 좌절된 시스템에서 국소화된 페르미온의 자기 상전이 거동을 연구하며, 특히 원자 극한에서 삼각 격자 위의 확장된 허바드 모델의 상도 구조를 규명합니다.
중요성:
- 삼각 격자는 고유한 기하학적 좌절 특성을 가지며, 이로 인해 반강자성 상호작용이 모든 상호작용을 동시에 만족할 수 없습니다
- 본 모델은 고전 이징 모델에서 페르미온 시스템으로의 비자명한 확장입니다
- 강상관 전자 시스템과 자기 재료의 이해에 중요한 의미를 가집니다
기존 방법의 한계:
- 초입방 격자에 대한 연구는 충분하지만, 삼각 격자의 정확한 해석 결과는 부족합니다
- 평균장 근사는 상전이 온도를 과대평가합니다
- 기하학적 좌절 효과에 대한 상세한 수치 연구가 부족합니다
연구 동기:
- 삼각 격자 위의 본 모델에 대한 수치 연구 공백을 채웁니다
- 기하학적 좌절이 상전이 거동에 미치는 영향을 탐색합니다
- 실험 재료(Li₂MnTeO₆, KCeO₂ 등)에 대한 이론적 기초를 제공합니다

핵심 기여

최초의 체계적 몬테카를로 연구: 삼각 격자 위의 이징형 자기 상호작용을 가진 확장된 허바드 모델에 대한 포괄적인 수치 시뮬레이션을 수행했습니다
완전한 상도 구성: 강자성(J<0)과 반강자성(J>0) 경우에 대한 온도-화학 포텐셜-상호작용 강도의 3차원 상도를 결정했습니다
상전이 성질 분석: 1차 및 2차 상전이의 경계를 정확히 결정하고, 반강자성 경우의 모든 상전이가 2차임을 발견했습니다
좌절 효과 규명: 기하학적 좌절이 상도 구조에 상당한 영향을 미치며, 반강자성 경우의 정렬 상 안정 온도가 강자성 경우보다 훨씬 낮음을 발견했습니다
단거리 정렬 상태 식별: 반강자성 경우에서 단거리 정렬 상태의 존재 영역을 발견하고 특성화했습니다

방법론 상세 설명

작업 정의

삼각 격자 위의 국소화된 페르미온 시스템의 상전이 거동을 연구하며, 입력은 모델 매개변수(J, U, μ, T)이고 출력은 시스템의 상도 및 상전이 성질입니다.

모델 구조

해밀토니안: $\hat{H} = U\sum_i \hat{n}_{i\uparrow}\hat{n}_{i\downarrow} + \frac{2J}{z}\sum_{\langle i,j\rangle} \hat{s}_i\hat{s}_j - \left(\mu + \frac{U}{2}\right)\sum_i \hat{n}_i$

여기서:

$\hat{s}_i = (\hat{n}_{i\uparrow} - \hat{n}_{i\downarrow})/2$ : 국소 스핀 연산자
$\hat{n}_i = \hat{n}_{i\uparrow} + \hat{n}_{i\downarrow}$ : 입자 수 연산자
$z = 6$ : 삼각 격자의 배위수

격점 상태: 각 격점은 네 가지 상태 중 하나를 가질 수 있습니다:

$|0\rangle$ : 공위
$|\uparrow\rangle$ : 단일 점유(스핀 위)
$|\downarrow\rangle$ : 단일 점유(스핀 아래)
$|\uparrow\downarrow\rangle$ : 이중 점유

기술적 혁신점

정준 앙상블 시뮬레이션: 입자 수 변동의 정준 앙상블을 채택하여 화학 포텐셜로 제어되는 시스템을 더 정확히 기술합니다
유한 크기 분석: L = 12×12에서 36×36까지의 유한 크기 효과를 체계적으로 연구했습니다
상전이 판정 방법:
- Challa 등의 에너지 분포 함수 P(E) 방법을 사용하여 1차 및 2차 상전이를 구분합니다
- 비열과 부분 격자 자화 강도의 스케일링 거동 분석을 결합합니다
부분 격자 분해: 삼각 격자를 세 개의 동등한 부분 격자로 분해하여 자기 정렬 구조를 정확히 분석합니다

실험 설정

데이터 집합

시스템 크기: 삼각형 클러스터, L = 12×12, 24×24, 36×36 격점
매개변수 범위:
- 온도: k_BT/|J| = 0.005 - 1.0
- 상호작용: U/|J| = -1.0 - 5.0
- 화학 포텐셜: μ/|J| = -0.6 - 0 (입자-공위 대칭성 활용)