Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0807.3615 · Journal of the Physical Society of Japan

Variational Monte Carlo Study of a Spinless Fermion t – V Model on a Triangular Lattice: Formation of a Pinball Liquid

Mitake Miyazaki, Chisa Hotta, Shin Miyahara, Keisuke Matsuda, Nobuo Furukawa

Published Jan 12, 2009• 17 citations
arXivHTMLPDF

Abstract

We analyze a model of spinless fermions on a triangular lattice at half-filling interacting via strong nearest-neighbor repulsive interactions, V, using the variational Monte Carlo simulation technique. The existence of three-sublattice long-range order is confirmed by the finite-size scaling analysis of the charge structural factor at V_c/t > 12. This ordered phase shows characteristics expected for a so called "pinball liquid" state, which has the spontaneous separation of fermionic degrees of freedom into coexisting Wigner crystal-like charge order (pin) and a metal (ball). The pins are fixed in order to maximize the kinetic energy gain of balls which move almost freely. The Fermi surface is reconstructed at V=V_c and remains towards the strong coupling limit. These features reminiscent of the strong correlation together with the large value of V_c/t distinguishes the pinball liquid from the conventional charge-density-wave.

DOI
10.1143/jpsj.78.014707

Related papers

Momentum distributions of a spinless fermiont-Vmodel on a triangular lattice in the strong coupling region

K Matsuda, Shin Miyahara, Nobuo Furukawa

Spinless fermion t-V model on a triangular lattice has been studied. In the model, the presence of a novel metallic charge-ordered ground state named a pinball liquid state has been proposed. To clarify the metallic features in the pinball state, we calculate the momentum distribution function of the model. For comparison, we also analyze a corresponding model where charge ordering is strongly pinned. The results are consistent with the picture expected for the pinball liquid, where an insulating ordered state (pin) and a metallic state on the honeycomb lattice (ball) coexist.

Journal of Physics Conference Series•Mar 1, 2009•0 citations

Ground state and spectral properties across a charge density wave transition in a triangular lattice spinless fermion model

Trithep Devakul, Rajiv R. P. Singh

We study ground state properties and particle excitation spectra across a commensurate charge density wave transition in a system of strongly interacting fermions, using series expansion methods and mean-field theory. We consider a 1/3-filled system of spinless fermions on a triangular lattice, with hopping parameter $t$, nearest-neighbor repulsion $V$, and a sublattice dependent chemical potential ${\ensuremath{\mu}}_{s}$. The phase transition is found to be first order for ${\ensuremath{\mu}}_{s}=0$, but becomes continuous with increasing ${\ensuremath{\mu}}_{s}$. The particle and hole excitation spectra exhibit dramatic changes in the vicinity of the phase transitions and in the charge density wave ordered phase. We discuss the relevance of this study to the pinball Fermi liquid phase postulated theoretically in earlier studies as well as to various strongly correlated triangular lattice materials.

arXiv:1409.5063•Physical Review B•Sep 22, 2014•2 citations

Particle-hole asymmetry and pinball liquid in a triangular-lattice extended Hubbard model within the mean-field approximation

Aleksey Alekseev, Agnieszka Cichy, Konrad Jerzy Kapcia

Recently, triangular lattice models have received a lot of attention since they can describe a number of strongly-correlated materials that exhibit superconductivity and various magnetic and charge orders. In this research we present an extensive analysis of the charge-ordering phenomenon of the triangular-lattice extended Hubbard model with repulsive onsite and nearest-neighbor interaction, arbitrary charge concentration, and $\sqrt{3}\times\sqrt{3}$ supercell (3-sublattice assumption). The model is solved in the ground state with the mean-field approximation which allowed to identify $8$ charge-ordered phases and a large variety of phase transitions. An exotic pinball-liquid phase was found and described. Moreover, strong particle-hole asymmetry of the phase diagram is found to play an important role for triangular lattices. The detailed analysis of band structures, unavailable for more advanced methods, such as dynamical mean-field theory, allowed us to interpret the found triangular-lattice phases and provided a great insight into the mechanisms behind the phase transitions that can also be met when correlation effects are taken into account. The complexity of the mean-field phase diagram showed the importance and usefulness of the results for the further research with correlation effects included. Together with atomic-limit approximation it can serve them as both a starting point, and a tool to interpret results.

arXiv:2504.20719•Physical review. B./Physical review. B•Sep 12, 2025•3 citations

Charge-Order on the Triangular Lattice: A Mean-Field Study for the Lattice S = 1/2 Fermionic Gas

Konrad Jerzy Kapcia

The adsorbed atoms exhibit tendency to occupy a triangular lattice formed by periodic potential of the underlying crystal surface. Such a lattice is formed by, e.g., a single layer of graphane or the graphite surfaces as well as (111) surface of face-cubic center crystals. In the present work, an extension of the lattice gas model to S=1/2 fermionic particles on the two-dimensional triangular (hexagonal) lattice is analyzed. In such a model, each lattice site can be occupied not by only one particle, but by two particles, which interact with each other by onsite U and intersite W1 and W2 (nearest and next-nearest-neighbor, respectively) density-density interaction. The investigated hamiltonian has a form of the extended Hubbard model in the atomic limit (i.e., the zero-bandwidth limit). In the analysis of the phase diagrams and thermodynamic properties of this model with repulsive W1>0, the variational approach is used, which treats the onsite interaction term exactly and the intersite interactions within the mean-field approximation. The ground state (T=0) diagram for W2≤0 as well as finite temperature (T>0) phase diagrams for W2=0 are presented. Two different types of charge order within 3×3 unit cell can occur. At T=0, for W2=0 phase separated states are degenerated with homogeneous phases (but T>0 removes this degeneration), whereas attractive W2<0 stabilizes phase separation at incommensurate fillings. For U/W1<0 and U/W1>1/2 only the phase with two different concentrations occurs (together with two different phase separated states occurring), whereas for small repulsive 0<U/W1<1/2 the other ordered phase also appears (with tree different concentrations in sublattices). The qualitative differences with the model considered on hypercubic lattices are also discussed.

arXiv:2105.01472•Nanomaterials•Apr 30, 2021•10 citations

Phase stability, charge ordering, and charge-liquid formation in the Falicov–Kimball model on triangular and kagome lattices

Ammar Nejati, Marvin Lenk

Abstract We present a systematic investigation of the stability of the phase diagram of the spinless Falicov–Kimball model on two-dimensional non-bipartite lattices, the triangular and kagome lattices, using Markov-chain Monte Carlo simulations. Three filling or chemical potential conditions are examined: fixed chemical potentials ( $$\mu _f = \mu _c = U/2$$ μ f = μ c = U / 2 ), “generalized” half-filling ( $$\bar{n}_f = 1/3$$ n ¯ f = 1 / 3 , $$\bar{n}_c = 2/3$$ n ¯ c = 2 / 3 ), and conventional half-filling ( $$\bar{n}_f = \bar{n}_c = 1/2$$ n ¯ f = n ¯ c = 1 / 2 ). For each condition, the effects of a perturbatively small next-nearest-neighbor hopping are studied. On the triangular lattice, Anderson insulator, Mott insulator, and charge-density wave (CDW) phases are found in all cases, with the CDW driven by Coulomb interaction-mediated nesting. Evidence for a charge-liquid regime is found in the region between the weakly correlated regime and the CDW, characterized by competing charge-ordering wavevectors and the absence of universal scaling behavior, whose extent shrinks progressively from the grand-canonical to the conventional half-filling condition. On the kagome lattice, the CDW phase is completely absent; instead, an insulating ground state, plausibly of the Mott type, is found, with signatures of a possible quantum phase transition at zero temperature. Our results establish the crucial role of lattice geometry and particle-hole symmetry breaking in shaping the phase diagram of correlated electron systems. Graphical abstract The 2d triangular lattice (left) and the 2d kagome lattice (right) with their unit cells highlighted. Bottom: The corresponding phase diagram of the Falicov-Kimball model (FKM) on each lattice, illustrating the non-Fermi liquid (nFL), charge-density wave (CDW), Anderson insulator (AI) and Mott insulator (MI) regimes.

arXiv:2608.28898•The European Physical Journal B•Sep 1, 2026•0 citations

Quantum Monte Carlo study of the one-dimensional Holstein model of spinless fermions

Ross H. McKenzie, Chris J. Hamer, D. W. Murray

The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (${\mathit{g}}^{2}$\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9512046•Physical review. B, Condensed matter•Apr 15, 1996•62 citations
We analyze a model of spinless fermions on a triangular lattice at half-filling interacting via strong nearest-neighbor repulsive interactions, V, using the variational Monte Carlo simulation technique. The existence of three-sublattice long-range order is confirmed by the finite-size scaling analysis of the charge structural factor at V_c/t > 12. This ordered phase shows characteristics expected for a so called "pinball liquid" state, which has the spontaneous separation of fermionic degrees of freedom into coexisting Wigner crystal-like charge order (pin) and a metal (ball). The pins are fixed in order to maximize the kinetic energy gain of balls which move almost freely. The Fermi surface is reconstructed at V=V_c and remains towards the strong coupling limit. These features reminiscent of the strong correlation together with the large value of V_c/t distinguishes the pinball liquid from the conventional charge-density-wave.
Spinless fermion t-V model on a triangular lattice has been studied. In the model, the presence of a novel metallic charge-ordered ground state named a pinball liquid state has been proposed. To clarify the metallic features in the pinball state, we calculate the momentum distribution function of the model. For comparison, we also analyze a corresponding model where charge ordering is strongly pinned. The results are consistent with the picture expected for the pinball liquid, where an insulating ordered state (pin) and a metallic state on the honeycomb lattice (ball) coexist.
We study ground state properties and particle excitation spectra across a commensurate charge density wave transition in a system of strongly interacting fermions, using series expansion methods and mean-field theory. We consider a 1/3-filled system of spinless fermions on a triangular lattice, with hopping parameter ttt, nearest-neighbor repulsion VVV, and a sublattice dependent chemical potential ${\ensuremath{\mu}}_{s}$. The phase transition is found to be first order for ${\ensuremath{\mu}}_{s}=0$, but becomes continuous with increasing ${\ensuremath{\mu}}_{s}$. The particle and hole excitation spectra exhibit dramatic changes in the vicinity of the phase transitions and in the charge density wave ordered phase. We discuss the relevance of this study to the pinball Fermi liquid phase postulated theoretically in earlier studies as well as to various strongly correlated triangular lattice materials.
Recently, triangular lattice models have received a lot of attention since they can describe a number of strongly-correlated materials that exhibit superconductivity and various magnetic and charge orders. In this research we present an extensive analysis of the charge-ordering phenomenon of the triangular-lattice extended Hubbard model with repulsive onsite and nearest-neighbor interaction, arbitrary charge concentration, and 3×3\sqrt{3}\times\sqrt{3}3​×3​ supercell (3-sublattice assumption). The model is solved in the ground state with the mean-field approximation which allowed to identify 888 charge-ordered phases and a large variety of phase transitions. An exotic pinball-liquid phase was found and described. Moreover, strong particle-hole asymmetry of the phase diagram is found to play an important role for triangular lattices. The detailed analysis of band structures, unavailable for more advanced methods, such as dynamical mean-field theory, allowed us to interpret the found triangular-lattice phases and provided a great insight into the mechanisms behind the phase transitions that can also be met when correlation effects are taken into account. The complexity of the mean-field phase diagram showed the importance and usefulness of the results for the further research with correlation effects included. Together with atomic-limit approximation it can serve them as both a starting point, and a tool to interpret results.
The adsorbed atoms exhibit tendency to occupy a triangular lattice formed by periodic potential of the underlying crystal surface. Such a lattice is formed by, e.g., a single layer of graphane or the graphite surfaces as well as (111) surface of face-cubic center crystals. In the present work, an extension of the lattice gas model to S=1/2 fermionic particles on the two-dimensional triangular (hexagonal) lattice is analyzed. In such a model, each lattice site can be occupied not by only one particle, but by two particles, which interact with each other by onsite U and intersite W1 and W2 (nearest and next-nearest-neighbor, respectively) density-density interaction. The investigated hamiltonian has a form of the extended Hubbard model in the atomic limit (i.e., the zero-bandwidth limit). In the analysis of the phase diagrams and thermodynamic properties of this model with repulsive W1>0, the variational approach is used, which treats the onsite interaction term exactly and the intersite interactions within the mean-field approximation. The ground state (T=0) diagram for W2≤0 as well as finite temperature (T>0) phase diagrams for W2=0 are presented. Two different types of charge order within 3×3 unit cell can occur. At T=0, for W2=0 phase separated states are degenerated with homogeneous phases (but T>0 removes this degeneration), whereas attractive W2<0 stabilizes phase separation at incommensurate fillings. For U/W1<0 and U/W1>1/2 only the phase with two different concentrations occurs (together with two different phase separated states occurring), whereas for small repulsive 0<U/W1<1/2 the other ordered phase also appears (with tree different concentrations in sublattices). The qualitative differences with the model considered on hypercubic lattices are also discussed.
Abstract We present a systematic investigation of the stability of the phase diagram of the spinless Falicov–Kimball model on two-dimensional non-bipartite lattices, the triangular and kagome lattices, using Markov-chain Monte Carlo simulations. Three filling or chemical potential conditions are examined: fixed chemical potentials ( μf=μc=U/2\mu _f = \mu _c = U/2μf​=μc​=U/2 μ f = μ c = U / 2 ), “generalized” half-filling ( nˉf=1/3\bar{n}_f = 1/3nˉf​=1/3 n ¯ f = 1 / 3 , nˉc=2/3\bar{n}_c = 2/3nˉc​=2/3 n ¯ c = 2 / 3 ), and conventional half-filling ( nˉf=nˉc=1/2\bar{n}_f = \bar{n}_c = 1/2nˉf​=nˉc​=1/2 n ¯ f = n ¯ c = 1 / 2 ). For each condition, the effects of a perturbatively small next-nearest-neighbor hopping are studied. On the triangular lattice, Anderson insulator, Mott insulator, and charge-density wave (CDW) phases are found in all cases, with the CDW driven by Coulomb interaction-mediated nesting. Evidence for a charge-liquid regime is found in the region between the weakly correlated regime and the CDW, characterized by competing charge-ordering wavevectors and the absence of universal scaling behavior, whose extent shrinks progressively from the grand-canonical to the conventional half-filling condition. On the kagome lattice, the CDW phase is completely absent; instead, an insulating ground state, plausibly of the Mott type, is found, with signatures of a possible quantum phase transition at zero temperature. Our results establish the crucial role of lattice geometry and particle-hole symmetry breaking in shaping the phase diagram of correlated electron systems. Graphical abstract The 2d triangular lattice (left) and the 2d kagome lattice (right) with their unit cells highlighted. Bottom: The corresponding phase diagram of the Falicov-Kimball model (FKM) on each lattice, illustrating the non-Fermi liquid (nFL), charge-density wave (CDW), Anderson insulator (AI) and Mott insulator (MI) regimes.
The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (g2{\mathit{g}}^{2}g2\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.