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

Localized fermions on the triangular lattice with Ising-like interactions

Basic Information

  • Paper ID: 2412.10808
  • Title: Localized fermions on the triangular lattice with Ising-like interactions
  • Authors: Lubomíra Regeciová (Slovak Academy of Sciences), Konrad Jerzy Kapcia (Adam Mickiewicz University)
  • Classification: cond-mat.stat-mech cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.comp-ph
  • Publication Date: October 17, 2025
  • Paper Link: https://arxiv.org/abs/2412.10808

Abstract

This study analyzes the localized fermion model on a triangular lattice in the canonical ensemble through Monte Carlo simulations. The system's Hamiltonian takes the form of the extended Hubbard model (atomic limit), containing nearest-neighbor Ising-like magnetic interactions J and on-site Coulomb interactions U. The investigation examines both signs of J, arbitrary U interactions, and arbitrary chemical potential μ (or equivalently, arbitrary particle density n). Based on analysis of specific heat and sublattice magnetization, the phase diagram of the model is determined. For the ferromagnetic case (J<0), the transition from the ordered phase can be either a second-order phase transition (at sufficiently high temperatures kBT/|J|≳0.2) or a first-order phase transition (at half-filling when -1<U/|J|<-0.65). For the J>0 case, ordered phases appear in the range -1/2<U/|J|<0 (at n=1), and short-range ordered states are also found for larger U values.

Research Background and Motivation

  1. Core Problem: Investigating the magnetic phase transition behavior of localized fermions in geometrically frustrated systems, particularly the phase diagram structure of the extended Hubbard model on the triangular lattice in the atomic limit.
  2. Significance:
    • The triangular lattice possesses inherent geometric frustration, preventing antiferromagnetic interactions from simultaneously satisfying all interaction constraints
    • This model represents a non-trivial extension of the classical Ising model to fermionic systems
    • It is crucial for understanding strongly correlated electron systems and magnetic materials
  3. Limitations of Existing Methods:
    • Research on hypercubic lattices is well-developed, but exact analytical results for the triangular lattice are lacking
    • Mean-field approximations overestimate phase transition temperatures
    • Detailed numerical studies of geometric frustration effects are scarce
  4. Research Motivation:
    • Filling the gap in numerical studies of this model on the triangular lattice
    • Exploring the influence of geometric frustration on phase transition behavior
    • Providing theoretical foundations for experimental materials (such as Li2MnTeO6, KCeO2, etc.)

Core Contributions

  1. First Systematic Monte Carlo Study: Comprehensive numerical simulations of the extended Hubbard model with Ising-like magnetic interactions on the triangular lattice
  2. Complete Phase Diagram Construction: Determination of the three-dimensional temperature-chemical potential-interaction strength phase diagram for both ferromagnetic (J<0) and antiferromagnetic (J>0) cases
  3. Phase Transition Analysis: Precise determination of first-order and second-order phase transition boundaries, revealing that all antiferromagnetic phase transitions are second-order
  4. Frustration Effects Revelation: Discovery that geometric frustration significantly affects the phase diagram structure, with ordered phase stability temperatures in the antiferromagnetic case being orders of magnitude lower than in the ferromagnetic case
  5. Short-Range Order Identification: Identification and characterization of short-range ordered state regions in the antiferromagnetic case

Methodology Details

Task Definition

Investigating the phase transition behavior of localized fermion systems on the triangular lattice, with model parameters (J, U, μ, T) as input and the system's phase diagram and phase transition properties as output.

Model Architecture

Hamiltonian: H^=U∑in^i↑n^i↓+2Jz∑⟨i,j⟩s^is^j−(μ+U2)∑in^i\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

where:

  • s^i=(n^i↑−n^i↓)/2\hat{s}_i = (\hat{n}_{i\uparrow} - \hat{n}_{i\downarrow})/2: local spin operator
  • n^i=n^i↑+n^i↓\hat{n}_i = \hat{n}_{i\uparrow} + \hat{n}_{i\downarrow}: particle number operator
  • z=6z = 6: coordination number of the triangular lattice

Lattice Site States: Each lattice site can occupy four states:

  • ∣0⟩|0\rangle: vacancy
  • ∣↑⟩|\uparrow\rangle: singly occupied (spin up)
  • ∣↓⟩|\downarrow\rangle: singly occupied (spin down)
  • ∣↑↓⟩|\uparrow\downarrow\rangle: doubly occupied

Technical Innovations

  1. Canonical Ensemble Simulation: Employing the canonical ensemble with particle number fluctuations for more accurate description of chemically controlled systems
  2. Finite-Size Analysis: Systematic investigation of finite-size effects for systems ranging from L = 12×12 to 36×36
  3. Phase Transition Identification Methods:
    • Using the energy distribution function P(E) method of Challa et al. to distinguish first-order and second-order phase transitions
    • Combining scaling behavior analysis of specific heat and sublattice magnetization
  4. Sublattice Decomposition: Decomposing the triangular lattice into three equivalent sublattices for precise analysis of magnetic ordering structures

Experimental Setup

Dataset

  • System Sizes: Triangular clusters with L = 12×12, 24×24, 36×36 lattice sites
  • Parameter Ranges:
    • Temperature: kBT/|J| = 0.005 - 1.0
    • Interactions: U/|J| = -1.0 - 5.0
    • Chemical potential: μ/|J| = -0.6 - 0 (utilizing particle-hole symmetry)

Evaluation Metrics

  1. Specific Heat: C=⟨ε2⟩−⟨ε⟩2/T2C = \langle\varepsilon^2\rangle - \langle\varepsilon\rangle^2/T^2
  2. Sublattice Magnetization: mα=(6/L)∑i∈α⟨s^i⟩m_\alpha = (6/L)\sum_{i\in\alpha}\langle\hat{s}_i\rangle
  3. Particle Density: n=⟨n^i⟩n = \langle\hat{n}_i\rangle
  4. Energy Distribution Function: P(E)P(E) for determining phase transition properties

Implementation Details

  • Algorithm: Metropolis Monte Carlo algorithm
  • Equilibration Steps: 10610^6 steps
  • Statistical Steps: 10610^6 steps
  • Temperature Sweep: Gradual cooling from high to low temperatures
  • Boundary Conditions: Periodic boundary conditions

Experimental Results

Main Results

Ferromagnetic Case (J < 0):

  • Highest phase transition temperature: kBT/|J| ≈ 0.605 (approaching the classical Ising model value of 0.607 as U → ∞)
  • Tricritical point location: kBT/|J| ≈ 0.18-0.20, nearly independent of U and μ
  • Phase transition nature: First-order at low temperatures, second-order at high temperatures
  • Ordered phase: Standard ferromagnetic phase with equal sublattice magnetizations

Antiferromagnetic Case (J > 0):

  • Highest phase transition temperature: kBT/|J| ≈ 0.06 (orders of magnitude lower than ferromagnetic case)
  • Ordered phase existence range: -1/2 < U/|J| < 0 (at half-filling)
  • Phase transition nature: All transitions are second-order
  • Ordered phase structure: Two sublattices form antiferromagnetic order, third sublattice has zero magnetization

Ablation Studies

Verification of thermodynamic limit behavior through system size variation:

  • First-order transitions: Bimodal structure in energy distribution
  • Second-order transitions: λ-type singularity in specific heat with peak growing with system size
  • Short-range order: Broad peaks in specific heat that do not scale with system size

Case Studies

Typical Phase Transition Behaviors:

  1. U = -0.7, μ = -0.10: First-order F-NO transition with discontinuous jumps in magnetization and particle number
  2. U = -0.6, μ = -0.12: Second-order F-NO transition with continuous order parameter evolution
  3. U = -0.2, μ = -0.15: Second-order AF-NO transition with antiferromagnetic ordering of two sublattices

Experimental Findings

  1. Strong Impact of Geometric Frustration: Phase transition temperatures in the antiferromagnetic case are orders of magnitude lower than in the ferromagnetic case
  2. Asymmetry in Phase Transition Nature: Ferromagnetic case exhibits both first-order and second-order transitions, while antiferromagnetic case shows only second-order transitions
  3. Ubiquity of Short-Range Order: Short-range ordered states are widely observed in the antiferromagnetic case

Related Work

  1. Classical Ising Model: Exact solutions on the triangular lattice (Houtappel, 1950) provide benchmarks for this study
  2. Extended Hubbard Model: Variational and numerical studies on hypercubic lattices provide methodological foundations
  3. Frustrated Magnetic Systems: Theoretical and experimental studies of triangular lattice antiferromagnets provide background for result interpretation
  4. Cold Atom Experiments: Realization of ultracold atomic gases in optical lattices provides possibilities for experimental verification of the model

Conclusions and Discussion

Main Conclusions

  1. Geometric frustration of the triangular lattice significantly affects phase transition behavior in localized fermion systems
  2. Ferromagnetic and antiferromagnetic interactions lead to completely different phase diagram structures and phase transition properties
  3. Ordered phases in the antiferromagnetic case possess unique two-sublattice antiferromagnetic structures
  4. All antiferromagnetic transitions are continuous, while ferromagnetic cases exhibit first-order to second-order transition crossovers

Limitations

  1. System Size Constraints: Computational resources limit maximum system size to 36×36
  2. Atomic Limit Approximation: Neglecting particle hopping terms may affect quantitative predictions for real materials
  3. Classical Spin Approximation: May be insufficient in regions with strong quantum fluctuations

Future Directions

  1. Investigation of longer-range interactions
  2. Application of quantum Monte Carlo methods
  3. Quantitative comparison with real materials
  4. Experimental verification in cold atom systems

In-Depth Evaluation

Strengths

  1. Rigorous Methodology: Employs mature Monte Carlo methods and established phase transition identification techniques
  2. Comprehensive Results: Constructs complete three-dimensional phase diagrams covering multiple parameter regions
  3. Deep Physical Insights: Reveals the fundamental impact of geometric frustration on phase transition behavior
  4. Strong Experimental Relevance: Connects to multiple real materials and cold atom systems

Weaknesses

  1. Insufficient Theoretical Analysis: Lacks deep theoretical explanation of phase transition mechanisms
  2. Missing Dynamical Properties: Does not address system dynamics and transport properties
  3. Neglected Quantum Effects: Atomic limit approximation ignores important quantum fluctuation effects

Impact

  1. Academic Value: Provides important numerical benchmarks for frustrated magnetic system theory
  2. Practical Value: Offers guidance for understanding and designing related materials
  3. Reproducibility: Detailed method descriptions ensure good reproducibility of results

Applicable Scenarios

  1. Theoretical modeling of triangular lattice magnetic materials
  2. Theoretical guidance for cold atom optical lattice experiments
  3. Phase transition studies in strongly correlated electron systems
  4. Benchmark testing for frustrated quantum many-body systems

References

This paper cites 71 important references covering Ising model theory, Hubbard model research, Monte Carlo methods, frustrated magnetic systems, and related experimental materials, providing a solid theoretical foundation and experimental background for the research.