Next Article in Journal
Highly Conductive Al/Al Interfaces in Ultrafine Grained Al Compact Prepared by Low Oxygen Powder Metallurgy Technique
Next Article in Special Issue
Energy Stored and Capacitance of a Circular Parallel Plate Nanocapacitor
Previous Article in Journal
The Effect of Quasi-Spherical Gold Nanoparticles on Two-Photon Induced Reactive Oxygen Species for Cell Damage
Previous Article in Special Issue
Scanning Thermal Microscopy of Ultrathin Films: Numerical Studies Regarding Cantilever Displacement, Thermal Contact Areas, Heat Fluxes, and Heat Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Konrad Jerzy Kapcia
Faculty of Physics, Adam Mickiewicz University in Poznań, ulica Uniwersytetu Poznańskiego 2, PL-61614 Poznań, Poland
Nanomaterials 2021, 11(5), 1181; https://doi.org/10.3390/nano11051181
Submission received: 25 March 2021 / Revised: 27 April 2021 / Accepted: 29 April 2021 / Published: 30 April 2021
(This article belongs to the Special Issue Novel Research in Low-Dimensional Systems)

Abstract

:
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 W 1 and W 2 (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 W 1 > 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 W 2 0 as well as finite temperature ( T > 0 ) phase diagrams for W 2 = 0 are presented. Two different types of charge order within 3 × 3 unit cell can occur. At T = 0 , for W 2 = 0 phase separated states are degenerated with homogeneous phases (but T > 0 removes this degeneration), whereas attractive W 2 < 0 stabilizes phase separation at incommensurate fillings. For U / W 1 < 0 and U / W 1 > 1 / 2 only the phase with two different concentrations occurs (together with two different phase separated states occurring), whereas for small repulsive 0 < U / W 1 < 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.

Graphical Abstract

1. Introduction

It is a well known fact that the classical lattice gas model is useful phenomenological model for various phenomena. It has been studied in the context of experimental studies of adsobed gas layers on crystaline substrates (cf., for example pioneering works [1,2,3,4]). For instance, the adsorbed atoms exhibit tendency to occupy a triangular lattice formed by periodic potential of the underlying crystal surface. This lattice is shown in Figure 1a. Such a lattice is formed by, e.g., a single layer of graphane or the graphite surface [i.e., the honeycomb lattice; (0001) hexagonal closed-packed (hcp) surface], and (111) face-centered cubic (fcc) surface. Atoms from (111) fcc surface are organized in the triangular lattice, whereas the triangular lattice is a dual lattice for the honeycomb lattice [5]. Note also that arrangements of atoms on (110) base-centered cubic (bcc) surface as well as on (111) bcc surface (if one neglects the interactions associated with other layers under surface) are quite close to the triangular lattice. One should mention that the triangular lattice and the honeycomb lattice are two examples of two-dimensional hexagonal Bravais lattices. Formally, the triangular lattice is a hexagonal lattice with a one-site basis, whereas the honeycomb lattice is a hexagonal lattice with a two-site basis. The classical lattice gas model is equivalent with the S = 1 / 2 Ising model in the external field [1,6,7,8,9] (the results for this model on the triangular lattice will be discussed in more details in Section 2).
In the present work, an extension of the lattice gas model to S = 1 / 2 fermionic particles is analyzed. Such a model has a form of the atomic limit of the extended Hubbard model [10], cf. (1). In this model, each lattice site can be occupied not by only one particle as in the model discussed in previous paragraph, but also by two particles. In addition to long-range (i.e., intersite) interactions between fermions, the particles located at the same site can also interact with each other via onsite Hubbard U interaction. For a description of the interacting fermionic particles on the lattice, the single-orbital extended Hubbard model with intersite density-density interactions has been used widely [11,12,13,14,15,16,17,18,19]. It is one of the simplest model capturing the interplay between the Mott localization (onsite interactions) and the charge-order phenomenon [15,16,17,18,19,20,21,22,23,24,25]. However, in some systems the inclusion of other interactions and orbitals is necessary [11,12,13,14,26,27].
This work can be palced among recent theoretical and experimental studies of adsorption of various atoms on (i) the (0001) hcp surface of the graphite [28,29,30,31,32] or of other materials [33,34] and (ii) the (111) fcc surface of metals and semimetals [34,35,36,37,38,39]. Although in the mentioned works the adsorbed particles on surface are rather classical and the analysis of classical lattice gas can give some predictions, taking into account of the quantum properties of adsorbed particles is necessary for, e.g., a description of experiments with He 4 and He 3 [28,40,41]. Moreover, there is plethora of recent experimental and theoretical studies of quasi-two-dimensional systems, e.g., Na x CoO 2 [42], NbSe 2 [43,44,45,46,47], TiSe 2 [48], TaSe 2 [49], VSe 2 [50], TaS 2 [51], and other transition metals dichalcogenides [52] as well as organic conductors [53,54], where various charge-ordered patterns have been observed on the triangular lattice. However, for such phenomena the atomic limit of the model studies is less reliable and more realistic description includes also electron hoping term as in the extended Hubbard model [15,16,17,18,19] or coupling with phonons as in the Holstein-Hubbard model [55]. In such cases, results obtained for atomic limit can be treated as a benchmark for models including the itinerant properties of fermionic particles.
The present work is organized as follows. In Section 2 the model and the methods (together with the most important equations) are presented. Section 3 is devoted to the discussion of ground state phase diagrams of the model with non-zero next-nearest neighbor interactions. Next, the finite temperature properties of the model with only the nearest-neighbor interactions are presented in Section 4. Finally, the most important conclusions and supplementary discussion are included in Section 5.

2. The Model and the Method

The extended Hubbard model in the zero-bandwidth limit (i.e., in the atomic limit) with interactions restricted to the second neighbors (or, equivalently, the next-nearest neighbors) can be expressed as:
H ^ = U i n ^ i n ^ i + 1 2 W 1 z 1 i , j 1 n ^ i n ^ j + 1 2 W 2 z 2 i , j 2 n ^ i n ^ j μ i n ^ i ,
where n ^ i = σ n ^ i σ , n ^ i σ = c ^ i σ c ^ i σ , and c ^ i σ ( c ^ i σ ) denotes the creation (annihilation) operator of an electron with spin σ at the site i. U is the onsite density interaction, W 1 and W 2 are the intersite density-density interactions between the nearest neighbors (NN) and the next-nearest neighbors (NNNs), respectively. z 1 and z 2 are numbers of NN and NNNs, respectively. μ is the chemical potential determining the total concentration n of electrons in the system by the relation n = ( 1 / L ) i n ^ i , where 0 n 2 and L is the total number of lattice sites. In this work phase diagrams emerging from this model are inspected. The analyses are performed in the grand canonical ensemble.
In this work the mean-field decoupling of the intersite term is used in the following form
n ^ i n ^ j = n ^ i n ^ j + n ^ i n ^ j n ^ i n ^ j ,
which is an exact treatment only in the limit of large coordination number ( z n ; or limit of infinite dimensions) [10,11,12,13,56,57,58]. Thus, for the two-dimensional triangular lattice (with z 1 = z 2 = 6 ) it is an approximation in the general case. It should be underlined that the treatment of the onsite term is rigorous in the present work. Please note that that the interactions U, W 1 and W 2 should be treated as effective parameters for fermionic particles including all possible contributions and renormalizations originating from other (sub-)systems.
Model (1) for W 2 0 has been intensively studied on the hypercubic lattices (see, e.g., [59,60,61,62,63,64,65,66,67,68] and references therein). Also the case of two-dimensional square lattice was investigated in detail for W 2 = 0 [61,62,63,64] as well as for W 2 0 [65,66,67,68,69]. There are also rigorous results for one-dimensional chain for W 2 = 0 [70,71] and W 2 0 [72].
In [73] the model with W 2 = 0 was investigated on the triangular lattice at half-filling by using a classical Monte Carlo method, and a critical phase, characterized by algebraic decay of the charge correlation function, belonging to the universality class of the two-dimensional X Y model with a Z 6 anisotropy was found in the intermediate-temperature regime. Some preliminary results for model (1) on the triangular lattice and for large attractive U < 0 and W 2 = 0 within the mean-field approximation were presented in [74].
The model in the limit U is equivalent with the S = 1 / 2 Ising model with antiferromagnetic (ferromagnetic) J n interactions if W n interaction in model (1) are repulsive, i.e., W n > 0 (attractive, i.e., W n < 0 , respectively). The relation between interaction parameters in both models is very simple, namely J n = W n . There is plethora of the results obtained for the Ising model on the triangular lattice. One should mention the following works (not assuming a comprehensive review): (a) exact solution in the absence of the external field H, i.e., for H = 0 (only with NN interactions, at arbitrary temperature) [5,75,76,77,78,79]; (b) for the model with NNN interactions included: ground state exact results [2], Bethe-Peierls approximation [1], Monte Carlo simulation both for H = 0 [80] and H 0 [3] (and other methods, e.g., [81,82]); (c) exact ground state results for the model with up to 3rd nearest-neighbor interactions for both H = 0 case [83] and H 0 case [4,84]. The most important information arising from these analyses is that only for W 2 0 (and arbitrary W 1 ) one can expect that consideration of 3 × 3 unit cells (i.e., tri-subblatice orderings) is enough to find all ordered states (particle arrangements) in the model. The reason is that the range of W 2 interaction is larger than the size of the unit cell. Thus, this is the point for that the present analysis of the model including only 3 × 3 unit cell orderings with restriction to W 2 0 is justified. One should not expect occurrence of any other phases beyond the tri-sublattice assumption in the studied range of the model parameters.
Please note that for W 2 > 0 it is necessary to consider a larger unit cell to find the true phase diagram of the model even in the U limit (cf., e.g., [2,4,81]). This is a similar situation as for model (1) on the square lattice, where for W 2 > 0 and any U not only checker-board order occurs (the two-sublattice assumption), but also other different arrangements of particles are present (the four-sublattice assumption, e.g., various stripes orders) [67,68].

2.1. General Definitions of Phases Existing in the Investigated System

In the systems analyzed only three nonequivalent homogeneous phases can exist (within the tree-sublattice assumption used). They are determined by the relations between concentrations n α ’s in each sublattice α ( n α = ( 3 / L ) i α n ^ i ), but a few equivalent solutions exist due to change of sublattice indexes. For intuitive understanding of rather complicated phase diagrams each pattern is marked with adequate abbreviation. The nonordered (NO) phase is defined by n A = n B = n C (all three n α ’s are equal), the charge-ordered phase with two different concentrations in sublattices (DCO phase) is defined by n A = n B n C , n B = n C n A , or n C = n A n B (two and only two out of three n α ’s are equal, 3 equivalent solution), whereas in the charge-ordered phase with three different concentrations in sublattices (TCO phase) n A n B , n B n C , and n A n C (all three n α ’s are different, 6 equivalent solutions). All these phases are schematically illustrated in Figure 1b. These phases exist in several equivalent solutions due to the equivalence of three sublattices forming the triangular lattice. Each of these patterns can be realized in a few distinct forms depending on specific electron concentrations on each sublattice (cf. Table 1 and Table 2 for T = 0 ). In addition, the degeneracy of the ground state solutions is contained in Table 1 (including charge and spin degrees of freedom).

2.2. Expressions for the Ground State

In the ground state (i.e., for T = 0 ), the grand canonical potential ω 0 per site of model (1) can be found as
ω 0 = H ^ / L = E D + E W + E μ ,
where contributions associated with the onsite interaction, the intersite interactions, and the chemical potential, respectively, has the following forms
E D = U 6 n A ( n A 1 ) + n B ( n B 1 ) + n C ( n C 1 ) ,
E W = 1 6 W 1 ( n A n B + n B n C + n C n A ) + 1 6 W 2 ( n A 2 + n B 2 + n C 2 ) ,
E μ = 1 3 μ ( n A + n B + n C ) .
In the above expressions, concentrations n α at T = 0 take the values from { 0 , 1 , 2 } set (cf. also Table 1). Please note that the above equations are the exact expressions for ω 0 of model (1) on the triangular lattice.
The free energy per site of homogeneous phases at T = 0 within the mean-field approximation is obtained as
f 0 = H ^ + μ i n ^ i / L = U D o c c + E W ,
where E W is expressed by (5). D o c c = ( 1 / L ) i n ^ i n ^ i denotes the double occupancy and this quantity is found to be exact, cf. Table 2. One should underline that above expression for f 0 is an approximate result for model (1) on the triangular lattice. Here, it is assumed that concentration n α are as defined in Table 2 and they are the same in each 3 × 3 unit cell in the system. Formally, it could be treated as exact one only if the numbers z n ( n = 1 , 2 ) goes to infinity.
The expressions presented in this subsection (for W 2 = 0 ) can be obtained as the T 0 limit of the equations for T > 0 included in Section 2.3.

2.3. Expressions for Finite Temperatures

For finite temperatures ( T > 0 ), the expressions given in [59] for the three-sublattice assumption takes the following forms (cf. also these in [68] given for the four-sublattice assumption). In approach used, the onsite U term is treated exactly and for the intersite W 1 term the mean-field approximation (2) is used. For a grand canonical potential ω (per lattice site) in the case of the lattice presented in Figure 1 one obtains
ω = 1 6 α Φ α n α 1 3 β α ln Z α .
where β = 1 / ( k B T ) is inverted temperature, coefficients Φ α are defined as Φ α = μ μ α ,
Z α = 1 + 2 exp β μ α + exp β 2 μ α U ,
and μ α is a local chemical potential in α sublattice ( α { A , B , C } )
μ A = μ 1 2 W 1 ( n B + n C ) , μ B = μ 1 2 W 1 ( n A + n C ) , μ C = μ 1 2 W 1 ( n A + n B ) .
For electron concentration n α in each sublattice in arbitrary temperature T > 0 one gets
n α = 2 Z α exp β μ α + exp β 2 μ α U ( for α { A , B , C } ) .
The set of three Equations (11) for n A , n B , and n C determines the (homogeneous) phase occurring in the system for fixed model parameters U, W 1 , and μ . If n = ( 1 / 3 ) ( n A + n B + n C ) is fixed, one has also set of three equations, but it is solved with respect to μ , n A , and n B (the third n α is obviously found as n C = 3 n n A n B ).
The free energy f per site is derived as
f = ω + 1 3 μ n A + n B + n C ,
where ω and n α ’s are expressed by (8)–(11).

2.4. Macroscopic Phase Separation

The free energy f P S of the (macroscopic) phase separated state (as a function of total electron concentration n; and at any temperature T 0 ) is calculated from
f P S ( n ) = n n n + n f + ( n + ) + n + n n + n f ( n ) ,
where f ± ( n ± ) are free energies of separating homogeneous phases with concentrations n ± . The factor before f ± ( n ± ) is associated with a fraction of the system, which is occupied by the phase with concentration n ± . Such defined phase separated states can exist only for n fulfilling the condition n < n < n + . For n ± only the homogeneous phase exists in the system (one homogeneous phase occupies the whole system). Concentrations n ± are simply determined at the ground state, whereas for T > 0 they can be found as concentrations at the first-order (discontinuous) boundary for fixed μ or by minimizing the free energy f P S [i.e., (13)] with respect to n + and n (for n fixed). For more details of the so-called Maxwell’s construction and macroscopic phase separations see, e.g., [59,68,85,86]. The interface energy between two separating phases is neglected here.

3. Results for the Ground State ( W 1 > 0 and W 2 0 )

3.1. Analysis for Fixed Chemical Potential μ

The ground state diagram for model (1) as a function of (shifted) chemical potential μ ¯ = μ W 1 W 2 is shown in Figure 2. The diagram is determined by comparison of the grand canonical potentials ω 0 ’s of all phases collected in Table 1 [cf. (3)]. It consists of several regions, where the NO phase occurs (3 regions: NO 0 , NO 1 and NO 2 ), the DCO phase occurs (6 regions: DCO 1 , DCO 2 , DCO 3 , DCO 1 * , DCO 2 * , and DCO 3 * ) and the TCO phase occurs (1 region).
All boundaries between the phases in Figure 2 are associated with a discontinuous change of at least one of the n α . The only boundaries associated with a discontinuous jump of two n α ’s are: DCO 2 –DCO 3 (DCO 2 * –DCO 3 * ) and TCO–NO 1 . At the boundaries ω 0 ’s of the phases are the same. It means that both phases can coexist in the system provided that a formation of the interface between two phases does not require additional energy. For W 2 = 0 , only the boundaries DCO 2 –DCO 3 (DCO 2 * –DCO 3 * ) and TCO–NO 1 have finite degeneracy (6 and 7, respectively, modulo spin degrees of freedom) and the interface between different types of 3 × 3 unit cells increases the energy of the system. Thus, the mentioned phases from neighboring regions cannot coexist at the boundaries. The other boundaries exhibit infinite degeneracy (it is larger than 3 · 2 L / 3 modulo spin) and entropy per site in the thermodynamic limit is non-zero. It means that at these boundaries both types of unit cells from neighboring regions can mix with any ratio and the formation of the interface between two phases does not change energy of the system. However, some conditions for arrangement of the cells can exist. For example, the DCO 2 phase with ( 0 , 0 , 2 ) can mix with the DCO 2 * phase with ( 0 , 2 , 2 ) or ( 2 , 0 , 2 ) , but not with the DCO 2 * phase with ( 2 , 2 , 0 ) . Please note that it is also possible to mix all three unit cells: ( 0 , 0 , 2 ) , ( 0 , 2 , 2 ) , and ( 2 , 0 , 2 ) . In such a case, ( 0 , 2 , 2 ) and ( 2 , 0 , 2 ) cells of the DCO 2 * phase cannot be located next to each other, i.e., they need to be separated by ( 0 , 0 , 2 ) unit cells of the DCO 2 phase. Thus, the degeneracy of the DCO 2 –DCO 2 * boundary is indeed larger than 3 · 2 L / 3 modulo spin. This is so-called macroscopic degeneracy, cf. [68]). In such a case, we say that the microscopic phase separation occurs. For W 2 < 0 these degeneracies are removed and all boundaries exhibit finite degeneracy (neglecting spin degrees of freedom). In this case the phases cannot be mixed on a microscopic level.
Please also note that for W 2 = 0 as well as for W 2 < 0 inside the regions shown in Figure 2, the 3 × 3 unit cells of the same type with different orientation cannot mix. It denotes that orientation of one type of the unit cell determines the orientation of other unit cells (of the same type). Thus, the degeneracy of the state of the system is finite (modulo spin) and the system exhibits the long-range order at the ground state inside each region of Figure 2. This is different from the case of two dimensional square lattice, where inside some regions different unit cells (elementary blocks) of the same phase can mix with each other [67,68].
One should underline that the discussed above ground state results for fixed chemical potential are the exact results for model (1) on the triangular lattice. This is due to the fact that the model is equivalent with a classical spin model, namely the S = 1 Blume-Cappel model with two-fold degenerated value of S = 0 (or the S = 1 classical Blume-Cappel with temperature-dependent anizotropy without degeneration), cf. [10,60,63]. For such a model, the mean-field approximation is an exact theory at the ground state and fixed external magnetic field (which corresponds to the fixed chemical potential in the model investigated).

3.2. Analysis for Fixed Particle Concentration n

The ground state diagram as a function of particle concentration n is shown in Figure 3. The rectangular regions are labeled by the abbreviations of homogeneous phases (cf. Table 2). At commensurate filling, i.e., i / 3 ( i = 0 , 1 , 2 , 3 , 4 , 5 , 6 ; but only on the vertical boundaries indicated in Figure 3) the homogeneous phase occurs, which can be found in Table 1 and Figure 2. On the horizontal boundaries the phases from both neighboring regions have the same energies.
For W 2 = 0 phase separated states (mentioned in the last column of Table 2) are degenerated with the corresponding homogeneous phases inside all regions of the phase diagram. This degeneracy can be removed in finite temperatures and in some regions the phase separated states can be stable at T > 0 (such regions are indicated by slantwise patter in Figure 3, cf. also Section 4). E.g., for W 2 = 0 , the TCO phase can exist only in the range of 0 < U / W 1 < 1 / 2 at T 0 . For W 2 < 0 the phase separated states have lower energies and they occur on the phase diagram (inside the rectangular regions of Figure 3). Obviously, at commensurate filling and for any W 2 0 , the homogeneous states can only occur (i.e., solid vertical lines in Figure 3). Please note that the following boundaries between homogeneous states (obtained by comparing only energies of homogeneous phases): (i) the DCO A and DCO B phases, (ii) the DCO A and DCO C phases, and (iii) the TCO A and TCO B phases are located at U / W 1 = 0 (and these corresponding for n > 1 ; the dashed line in Figure 3). For W 2 < 0 these lines do not overlap with the boundaries between corresponding phase separated states at U / W 1 | k | = 0 (or U / | W 2 | = 1 ), but in such a case the homogeneous states have higher energies than the phase separated states. In fact, the homogeneous states for W 2 < 0 are unstable (i.e., μ / n < 0 ) inside the regions of Figure 1. For W 2 < 0 they are stable only for commensurate fillings (solid lines in Figure 3).
For the system on the square lattice the similar observation can be made (Figure 1 from [59])—compare HCO A –LCO A and HCO A –HCO B boundaries at U / W 1 = 0 with PS1 A –PS1 B and PS1 A –PS1 B boundaries at U / | W 2 | = 1 , respectively. In [68] the boundaries between homogeneous phases for W 2 < 0 are not shown in Figure 3. Only boundaries between corresponding phase separated states are correctly presented in that figure for W 2 < 0 . For U / W 1 > 0 , the CBO A phase (corresponding to the HCO A phase from [59]) is not the phase with the lowest energy (among homogeneous phases) in any range of n (but for U / W 1 < 0 it has the lowest energy among all homogeneous states). However, the corresponding phase separated state NO 0 /CBO 2 (i.e., PS1 A from [59]) can occur for U / W 1 > 0 (and for U / | W 2 | < 1 ) as shown in Figure 3 of [68].
The vertical boundaries for homogeneous phases (i.e., the transitions with changing n) are associated with continuous changes of all n α ’s and D o c c , but the chemical potential μ (calculated as μ = f / n ) changes discontinuously. Boundaries DCO A –DCO B , DCO A –DCO C , and TCO A –TCO B (and other transitions for fixed n at U / W 1 = 0 ) between homogeneous phases are associated with discontinuous change of only D o c c . One should note that it is similar to transition between two checker-board ordered phases on the square lattice, namely CBO A –CBO B and CBO A –CBO C boundaries, cf. [68] (or the HCO A –LCO A and HCO A –HCO B boundaries, respectively, from [59]). At the other horizontal boundaries (i.e., transitions for fixed n at U / W 1 | k | = 1 / 2 in Figure 3) two of n α ’s and D o c c change discontinuously. At commensurate fillings transitions with changing U / W 1 occur only at points indicated by squares in Figure 3.
All horizontal boundaries between phase separated states (which are stable for W 2 < 0 ) are connected with discontinuous changes of D o c c . These boundaries located at U / W 1 | k | = 0 are also associated to a discontinuous change of particle concentration in one of the domains.
The diagram presented in Figure 3 is constructed by the comparison of (free) energies of various homogeneous phases and phase separated states collected in Table 1. The energies of homogeneous phases are calculated from (7), whereas energies of phase separated states are calculated from (13). Please note that it is easy to calculate energies of f ± ( n ± ) of separating homogeneous phase (with commensurate fillings) at the ground state by just taking μ = 0 in ω 0 ’s collected in Table 1. Obviously, one can also calculate energies of the phases collected in Table 2 at these fillings (from both neighboring regions). For example, the DCO B phase and the DCO C phase at n = 1 / 3 reduce to DCO 1 phase.

4. Results for Finite Temperatures ( W 1 > 0 and W 2 = 0 )

One can distinguish four ranges of U interaction, where the system exhibits qualitatively different behavior, namely: (i) U / W 1 < 0 , (ii) 0 < U / W 1 < ( 1 / 3 ) ln ( 2 ) , (iii) ( 1 / 3 ) ln ( 2 ) < U / W 1 < 1 / 2 , and (iv) U / W 1 > 1 / 2 . In Figure 4, Figure 5, Figure 6 and Figure 7, the exemplary finite temperature phase diagrams occurring in each of these ranges of onsite interaction are presented. All diagrams are found by investigation of the behavior of n α ’s determined by (11) in the solution corresponding to the lowest grand canonical potential [Equation (8), when μ is fixed] or to the lowest free energy [Equations (12) and (13) if n is fixed]. The set of three nonlinear Equations (11) has usually several nonequivalent solutions and thus it is extremely important to find a solution, which has the minimal adequate thermodynamic potential. In Figure 8 the behavior of n α ’s as a function of temperature or chemical potential is shown for some representative model parameters. Figure 9 presents the phase diagram of the system for half-filling.
For U / W 1 < 0 and U / W 1 > 1 / 2 the phase diagrams of the model are similar and the DCO phase is only ordered homogeneous one occurring on the diagrams. In the first range, there are two regions of ordered phase occurrence (cf. Figure 4 and [74]), whereas in the second case one can distinguish four regions of the DCO phase stability (cf. Figure 5). The NO–DCO transitions for fixed μ are discontinuous for any values of onsite interaction and chemical potential in discussed range of model parameters and thus phase separated state PS 1 :NO/DCO occurs in define ranges of n. In this state domains of the NO and the DCO phases coexist.
For U / W 1 < 0 the temperature of NO–DCO transition is maximal for μ ¯ = 0 (i.e., at half-filling)–Figure 4a. Its maximal value T M monotonously decreases with increasing of U from k B T M / W 1 = 1 / 2 for U and at U = 0 it is equal to 1 / 4 . This transition exhibits re-entrant behavior (for fixed | μ ¯ | > 1 ). At T = T M and μ ¯ = 0 and at only this point, this transition exhibits properties of a second order transition (cf. Figure 8a). In particular, with increasing T for fixed μ ¯ = 0 n α ’s changes continuously at T M , but two equivalent solutions still exist for any T < T M (similarly as in the ferromagnetic Ising model at zero field [9]). At μ ¯ = 0 and T < T M the discontinuous transition between two DCO phases occurs. In the DCO phase for μ ¯ < 0 ( n < 1 ) [connecting with the DCO 1 (DCO A ) region at T = 0 ] the relation n A = n B < n C is fulfilled, whereas in the DCO phase for μ ¯ > 0 ( n > 1 ) [connecting with the DCO 1 * (DCO A * ) region at T = 0 ] the relation n A < n B = n C occurs ( n C can be larger than 1 for some temperatures), cf. also Figure 8g,h as well as [74]. Both discontinuous transitions for fixed chemical potential are associated with occurrence of phase separated states. On the diagrams obtained for fixed n three region of phase separated states occurs (Figure 4b). For W 2 = 0 the PS 1 :NO/DCO phase separated state occurs only for T > 0 . For T 0 the concentrations in both domains of the PS 1 state approach 0 (or 2), whereas for T T M they approach to 1. Near n = 1 the PS 2 :DCO/DCO state is stable for 0 T < T M . In this state domains of two DCO phases (with different particle concentrations) coexist in the system.
For U / W 1 > 1 / 2 the diagrams are similar, but the double occupancy of sites is strongly reduced due to repulsive U (Figure 5). Thus, their structure exhibits two lobs of the DCO phase occurrence in cotrary to the case of U / W 1 < 0 , where a single lob of the DCO phase is present (as expected from previous studies of the model, cf. [10,59,60]). The maximal value k B T M / W 1 of NO–DCO transition occurs for μ ¯ / W 1 corresponding approximately quarter fillings (i.e., near n = 1 / 2 and n = 3 / 2 ). With increasing U it decreases and finally in the limit U + it reaches 1 / 8 . At this point DCO–NO boundary exhibits features of continuous transition as discussed previously. In this range, the phase diagrams are (almost) symmetric with respect to these fillings (when one considers only one part of the diagram for 0 < n < 1 or for 1 < n < 2 ).
The most complex diagrams are obtained for 0 < U / W 1 < 1 / 2 , where the TCO phase appears at T = 0 and for finite temperatures near half-filling. For 0 < U / W 1 < ( 1 / 3 ) ln ( 2 ) the region of the TCO phase is separated from the NO phase by the region of DCO phase, Figure 6a. The TCO–DCO transition is continuous (cf. Figure 8g,h for U / W 1 = 0.35 ) and its maximal temperature is located for half-filling (at μ ¯ = 0 or n = 1 ). At this point two first-order NO–DCO and two second-order TCO–DCO boundaries merge (for fixed chemical potential). It is the only point for fixed U / W 1 in this range of model parameters, where a direct continuous transition from the TCO phase to the NO phase is possible (Figure 8b). The continuous TCO–DCO transition temperature can be also found as a solution of (11) and (A4) as discussed in Appendix A. Similarly as for U / W 1 < 0 , the temperature of NO–DCO transition is maximal at half-filling. For fixed n, the narrow regions of PS 1 :NO/DCO states are present between the NO region and DCO regions. Please note that for T > 0 there is no signatures of the discontinuous DCO 1 –DCO 2 (DCO 1 * –DCO 2 * ) boundary occurring at T = 0 . It is due to the fact that the discontinuous jumps of n α ’s occurring for T = 0 at these boundaries are changed into continuous evolutions of sublattice concentrations at T > 0 and there is no criteria for distinction of these two DCO phases at finite temperatures (cf. also [59,60,61]). From the same reason, there is no boundary at T > 0 for fixed n associated to the DCO B –DCO C (DCO B * –DCO C * ) line occurring at T = 0 (Figure 6b). However, strong reduction of one n α from the case where n α 2 to the case of n α 1 is visible (some kind of a smooth crossover inside the DCO region), cf. Figure 8f–h for U / W 1 = 0.35 .
For ( 1 / 3 ) ln ( 2 ) < U / W 1 < 1 / 2 , the maximum of the NO–DCO transition temperature is shifted towards larger | μ ¯ | / W 1 (or smaller | 1 n | ). This is associated with forming of the two-lob structure of the diagram found for U / W 1 > 1 / 2 . Inside the regions of the DCO phase occurrence discontinuous transitions between two DCO phases appear—See Figure 7a as well as Figure 8e,i. These new regions of the DCO phase at T > 0 [with n A < n B = n C (for μ ¯ < 0 or n < 1 ); cf. Figure 8e,i] are connected with the DCO 3 and DCO 3 * regions occurring at the ground state. The boundaries DCO–DCO weakly dependent on μ ¯ are associated with occurrence of phase separated PS 2 :DCO/DCO states (at high temperatures) in some ranges of n, cf. Figure 7b. The other DCO–DCO transitions (which are almost temperature-independent) are not connected with phase separated states. Also the first-order TCO–NO line is present near half-filling, cf. Figure 8d. One should underline that all four lines (three first-order boundaries: DCO–NO, DCO–DCO, TCO–NO and the second-order TCO–DCO boundary) merge at single point with numeric accuracy. However, it cannot be excluded that the DCO–NO and TCO–DCO boundaries connect with the temperature-independent line in slightly different points, what could result in, e.g., the TCO–DCO–NO sequence of transition with increasing temperature for small range of chemical potential μ ¯ . All of these almost temperature-independent boundaries (i.e., the DCO–DCO and the TCO–NO lines) are located at temperature, which decreases with increasing U / W 1 and approaches 0 at U / W 1 = 1 / 2 [i.e., they connect with the DCO 2 –DCO 3 (DCO 2 * –DCO 3 * ) and TCO–NO 1 boundaries at T = 0 for fixed μ or with the TCO–DCO D (TCO * –DCO D * ) lines at T = 0 for fixed n]. From the analysis of (11) similarly as it was done in the case of the square lattice [10] (see also Appendix A) one obtains that the point, where the TCO–NO transition changes its order at half-filling, is k B T / W 1 = 1 / 6 and U / W 1 = ( 1 / 3 ) ln ( 2 ) .
For better overview of the system behavior, the phase diagram of the model for half-filling ( μ ¯ = 0 or n = 1 ) is presented in Figure 9. The temperature of the order-disorder transition decreases with increasing U / W 1 . In low temperatures and for U / W 1 < 0 , the DCO phases exist in the system (precisely, if μ is fixed—at μ ¯ = 0 the DCO–DCO discontinuous boundary occurs; whereas if n is fixed—the PS 2 :DCO/DCO state is stable at n = 1 ), cf. also Figure 4. For 0 < U / W 1 < 1 / 2 the TCO phase is stable below the order-disorder line, but for ( 1 / 3 ) ln ( 2 ) < U / W 1 < 1 / 2 and k B T / W 1 < 1 / 6 the TCO–NO phase transition is discontinuous (cf. also Figure 7). For U / W 1 < ( 1 / 3 ) ln ( 2 ) the order-disorder boundary presented in Figure 9 is a merging point of several boundaries as presented in Figure 4 and Figure 6, and discussed previously. Thus, formally this order-disorder boundary for U / W 1 < ( 1 / 3 ) ln ( 2 ) occurring at half-filling is a line of some critical points of a higher order.
Please note that the order-disorder transition is discontinuous for any value of onsite interaction and chemical potential [excluding only the TCO–NO boundary for half-filling and 0 < U / W 1 < ( 1 / 3 ) ln ( 2 ) ] in contrast to the case of two- [10,59,60] or four-sublattice [67,68] assumptions, where it can be continuous one for some range of model parameters). In [74] also metastable phases have been discussed in detail for the large onsite attraction limit and the triangular lattice.

5. Final Remarks

In this work, the mean-field approximation was used to investigate the atomic limit of extended Hubbard model [hamiltonian (1)] on the triangular lattice. The phase diagram was determined for the model with intersite repulsion between the nearest neighbors ( W 1 > 0 ). The effects of attractive next-nearest-neighbor interaction ( W 2 < 0 ) were discussed in the ground state. The most important findings of this work are that (i) two different arrangements of particles (i.e., two different charge-ordered phases: the DCO and TCO states) can occur in the system and (ii) attractive W 2 < 0 or finite T > 0 removes the degeneration between homogeneous phases and phase separated states occurring at T = 0 for W 2 = 0 . It was shown that TCO phase is stable in intermediate range of onsite repulsion 0 < U / W 1 < 1 / 2 (for W 2 = 0 ). All transition from the ordered phases to the NO are discontinuous for fixed chemical potential (apart from TCO–NO boundary at half-filling for 0 < U / W 1 < ( 1 / 3 ) ln ( 2 ) ) and the DCO–NO boundaries at single points corresponding to n = 1 / 2 , 1 , 3 / 2 as discussed in Section 4), thus the phase separated states occur on the phase diagram for fixed particle concentration.
One should stress that hamiltonian (1) is interesting not only from statistical point of view as a relatively simple toy model for phase transition investigations. Although it is oversimplified for quantitative description of bulk condensed matter systems, it can be useful in qualitative analysis of, e.g., experimental studies of adsorbed gas layers on crystalline substrates.
Additionally, one notes that the mean-field results for model (1) with attractive W 1 < 0 and W 2 0 are the same for both two-sublattice and tri-sublattice assumptions. In such a case, three different nonordered phases exist with the discontinuous first-order transition between them (at μ ¯ = 0 for U < 0 or for | μ ¯ | 0 for U / ( | W 1 | + | W 2 | ) > 1 ), and thus for fixed n, several so-called electron-droplet states (phase separations NO/NO) exist (cf. [60,68,87,88], particularly Figure 2 of [60]).
Notice that the mean-field decoupling of the intersite term is an approximation for purely two-dimensional model investigated, which overestimates the stability of ordered phases. For example, the order-disorder transition for the ferromagnetic Ising model is overestimated by the factor two (for the honeycomb, square and triangular lattices rigorous solution gives k B T c / | J | as 0.506 , 0.568 , 0.607 , respectively, whereas the mean-field approximation gives k B T c / | J | = 1 ) [76]. Moreover, the results for the antiferromagnetic Ising model on the triangular lattice [the limit U ± of model (1)] do not predict long-range order at zero field [1,3,76] and T > 0 [corresponding to n = 1 or n = 1 / 2 , 3 / 2 , respectively, in the case of model (1)]. However, longer-range interactions [3] or weak interactions between adsorbed particles and the adsorbent material occurring in realistic systems could stabilize such an order (such systems are rather quasi-two-dimensional). It should be also mentioned that the charge Berezinskii-Kosterlitz-Thouless-like phase was found in the intermediate-temperature regime between the charge-ordered phase (with long-range order, coresponding to the TCO phase here) and disordered phases in the investigated model [73].
The recent progress in the field of optical lattices and a creation of the triangular lattice by laser trapping [89,90] could enable testing predictions of the present work. The fermionic gases in harmonic traps are fully controllable systems. Note also that the superconductivity in the twisted-bilayer graphene [91,92,93,94,95,96] is driven by the angle between the graphene layers. It is associated with an occurrence of the Moiré pattern (the triangular lattice with very large supercell). Hetero-bilayer transition metals dichalcogenides system is the other field where this pattern appears [97,98]. This makes further studies of properties of different models on the triangular lattice desirable.

Author Contributions

Conceptualization, K.J.K.; methodology, K.J.K.; software, K.J.K.; validation, K.J.K.; formal analysis, K.J.K.; investigation, K.J.K.; resources, K.J.K.; data curation, K.J.K.; writing—original draft preparation, K.J.K.; writing—review and editing, K.J.K.; visualization, K.J.K.; supervision, K.J.K.; project administration, K.J.K.; funding acquisition, K.J.K. The author has read and agreed to the published version of the manuscript.

Funding

The support from the National Science Centre (NCN, Poland) under Grant SONATINA 1 no. UMO-2017/24/C/ST3/00276 is acknowledged. Founding in the frame of a scholarship of the Minister of Science and Higher Education (Poland) for outstanding young scientists (2019 edition, no. 821/STYP/14/2019) is also appreciated.

Institutional Review Board Statement

Not applicable. The study does not involve humans or animals.

Informed Consent Statement

Not applicable. The study does not involve humans.

Data Availability Statement

The data presented in this study are available on request from the author. All data presented in Section 4 have been obtained by numerical solving of self-consistent equations given in Section 2.3.

Acknowledgments

The author expresses his sincere thanks to J. Barański, R. Lemański, R. Micnas, P. Piekarz, and A. Ptok for very useful discussions on some issues raised in this work. The author also thanks R. Micnas and I. Ostrowska for careful reading of the manuscript.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
NONon-ordered phase with the same concentrations in all three sublattices
DCOCharge-ordered phase with two different concentrations in sublattices
TCOCharge-ordered phase with three different concentrations in sublattices

Appendix A. Analytic Expressions for Continuous Transition Temperatures

Equations (11) can be written in a different form, namely n α = f α , where f α 2 ( t α + t α 2 a ) / ( 1 + 2 t α + t α 2 a ) , t α exp ( β μ α ) and a exp ( β U ) . One can define Δ ( n A n B ) / 2 and χ ( n B n C ) / 2 . From (10) one gets:
μ A = μ W 1 n 1 3 2 Δ + χ ,
μ B = μ W 1 n + 1 3 Δ χ ,
μ C = μ W 1 n + 1 3 Δ + 2 χ .
Taking the limit χ 0 of both sides of the equation ( f B f C ) / ( 2 χ ) = 1 (using de l’Hospital theorem) one gets ( g B g C ) / 2 = 1 , where g α f α χ = f α t α t α μ α μ α χ . One easily finds that f α / t α = 2 ( 1 + 2 t α a + t α 2 a ) / ( 1 + 2 t α + t α 2 a ) 2 , t α / μ α = β t α as well as μ A / χ = W 1 / 3 , μ B / χ = W 1 / 3 , μ C / χ = 2 W 1 / 3 . Finally, the equation determining temperature T c of a continuous transition (at which n B n C ) has the form
1 β c W 1 = 1 + 2 t B C a ¯ + t B C 2 a ¯ t B C 1 + 2 t B C + t B C 2 a ¯ 2 ,
where t B C exp ( β c μ B C ) , μ B C = μ W 1 ( n + Δ / 3 ) (in the considered limit μ B = μ C and n B = n C ), a ¯ exp ( β c U ) , β c 1 / ( k B T c ) . Concentrations n A and n B C n B = n C are calculated from (11) for β c self-consistently. Thus, for fixed μ (or n) one has a set of three equation which is solved with respect to β c , n (or μ ) and Δ .
The solutions of (A4) and (11) with Δ 0 (i.e., n A n B ) correspond to the TCO–DCO boundaries. Such determined temperatures coincide with those found from the analysis of (11) and (8) or (12) and presented in Figure 6 and Figure 7, what supports the findings that the TCO–DCO boundaries are indeed continuous.
The solutions of (A4) and (11) with Δ = 0 (i.e., n A = n B ) correspond to the continuous DCO–NO boundaries. On the diagrams presented in Section 4 such solutions for T c are located inside the regions of the DCO phase occurrence (and they correspond to the transitions between metastable phases [74] or to vanishing of the NO metastable solution, cf. [88,99]). In the present case of model (1) studied, they coincide with the DCO–NO transitions presented in Figure 4, Figure 5, Figure 6 and Figure 7 only at T = 0 (i.e., for n = 0 , 2 as well as for n = 1 and U / W 1 > 1 / 2 ; or corresponding μ ¯ ) and at T = T M (i.e., maximal temperature of the DCO–NO transition, occurring for U / W 1 < ( 1 / 3 ) ln ( 2 ) and n = 1 or μ ¯ = 0 , as well as for U / W 1 > 1 / 2 and n 1 / 2 , 3 / 2 or corresponding μ ¯ ; for ( 1 / 3 ) ln ( 2 ) < U / W 1 < 1 / 2 it is located for some intermediate concentrations 1 / 2 < n < 1 and 1 < n < 3 / 2 ). For Δ = 0 , (A4) and (11) give the following results: (i) for U : k B T c / W 1 = n ( 2 n ) / 2 ; (ii) for U = 0 : k B T c / W 1 = n ( 2 n ) / 4 ; and (iii) for U + : k B T c / W 1 = n ( 1 n ) / 2 (if n < 1 ) and k B T c / W 1 = ( 2 n ) ( n 1 ) / 2 (if n > 1 ). Please note that such determined T c for Δ = 0 is two times smaller than corresponding continuous transitions for the model considered on the hypercubic lattice within the mean-field aprroximation for the intersite term (for the same U / W 1 and n) [10,59,60].

References

  1. Campbell, C.E.; Schick, M. Triangular Lattice Gas. Phys. Rev. A 1972, 5, 1919–1925. [Google Scholar] [CrossRef]
  2. Kaburagi, M.; Kanamori, J. Ordered Structure of Adatoms in the Extended Range Lattice Gas Model. Japan. J. Appl. Phys. 1974, 13 (Suppl. S2), 145–148. [Google Scholar] [CrossRef] [Green Version]
  3. Mihura, B.; Landau, D.P. New Type of Multicritical Behavior in a Triangular Lattice Gas Model. Phys. Rev. Lett. 1977, 38, 977–980. [Google Scholar] [CrossRef]
  4. Kaburagi, M.; Kanamori, J. Ground State Structure of Triangular Lattice Gas Model with up to 3rd Neighbor Interactions. J. Phys. Soc. Jpn. 1978, 44, 718–727. [Google Scholar] [CrossRef]
  5. Wannier, G.H. The Statistical Problem in Cooperative Phenomena. Rev. Mod. Phys. 1945, 17, 50–60. [Google Scholar] [CrossRef]
  6. Ising, E. Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 1925, 31, 253–258. [Google Scholar] [CrossRef]
  7. Onsager, L. Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition. Phys. Rev. 1944, 65, 117–149. [Google Scholar] [CrossRef]
  8. Binney, J.J.; Dowrick, N.J.; Fisher, A.J.; Newman, M.E.J. The Theory of Critical Phenomena: An Introduction to the Renormalization Group; Oxford University Press: Oxford, UK, 1992. [Google Scholar]
  9. Vives, E.; Castán, T.; Planes, A. Unified Mean-Field Study of Ferro-and Antiferromagnetic Behavior of the Ising Model with External Field. Amer. J. Phys. 1997, 65, 907–913. [Google Scholar] [CrossRef]
  10. Micnas, R.; Robaszkiewicz, S.; Chao, K.A. Multicritical Behavior of the Extended Hubbard Model in the Zero-Bandwidth Limit. Phys. Rev. B 1984, 29, 2784–2789. [Google Scholar] [CrossRef]
  11. Micnas, R.; Ranninger, J.; Robaszkiewicz, S. Superconductivity in Narrow-Band Systems with Local Nonretarded Attractive Interactions. Rev. Mod. Phys. 1990, 62, 113–171. [Google Scholar] [CrossRef]
  12. Georges, A.; Kotliar, G.; Krauth, W.; Rozenberg, M.J. Dynamical Mean-Field Theory of Strongly Correlated Fermion Systems and the Limit of Infinite Dimensions. Rev. Mod. Phys. 1996, 68, 13–125. [Google Scholar] [CrossRef] [Green Version]
  13. Imada, M.; Fujimori, A.; Tokura, Y. Metal-Insulator Transitions. Rev. Mod. Phys. 1998, 70, 1039–1263. [Google Scholar] [CrossRef] [Green Version]
  14. Kotliar, G.; Savrasov, S.Y.; Haule, K.; Oudovenko, V.S.; Parcollet, O.; Marianetti, C.A. Electronic Structure Calculations with Dynamical Mean-Field Theory. Rev. Mod. Phys. 2006, 78, 865–951. [Google Scholar] [CrossRef] [Green Version]
  15. Davoudi, B.; Hassan, S.R.; Tremblay, A.M.S. Competition Between Charge and Spin Order in the t-U-V Extended Hubbard Model on the Triangular Lattice. Phys. Rev. B 2008, 77, 214408. [Google Scholar] [CrossRef] [Green Version]
  16. Cano-Cortés, L.; Ralko, A.; Février, C.; Merino, J.; Fratini, S. Geometrical Frustration Effects on Charge-Driven Quantum Phase Transitions. Phys. Rev. B 2011, 84, 155115. [Google Scholar] [CrossRef] [Green Version]
  17. Merino, J.; Ralko, A.; Fratini, S. Emergent Heavy Fermion Behavior at the Wigner-Mott Transition. Phys. Rev. Lett. 2013, 111, 126403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Tocchio, L.F.; Gros, C.; Zhang, X.F.; Eggert, S. Phase Diagram of the Triangular Extended Hubbard Model. Phys. Rev. Lett. 2014, 113, 246405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Litak, G.; Wysokiński, K.I. Evolution of the Charge Density Wave Order on the Two-Dimensional Hexagonal Lattice. J. Magn. Magn. Mater. 2017, 440, 104–107. [Google Scholar] [CrossRef]
  20. Aichhorn, M.; Evertz, H.G.; von der Linden, W.; Potthoff, M. Charge Ordering in Extended Hubbard Models: Variational Cluster Approach. Phys. Rev. B 2004, 70, 235107. [Google Scholar] [CrossRef] [Green Version]
  21. Tong, N.H.; Shen, S.Q.; Bulla, R. Charge Ordering and Phase Separation in the Infinite Dimensional Extended Hubbard Model. Phys. Rev. B 2004, 70, 085118. [Google Scholar] [CrossRef] [Green Version]
  22. Amaricci, A.; Camjayi, A.; Haule, K.; Kotliar, G.; Tanasković, D.; Dobrosavljević, V. Extended Hubbard Model: Charge Ordering and Wigner-Mott Transition. Phys. Rev. B 2010, 82, 155102. [Google Scholar] [CrossRef] [Green Version]
  23. Ayral, T.; Biermann, S.; Werner, P.; Boehnke, L. Influence of Fock Exchange in Combined Many-Body Perturbation and Dynamical Mean Field Theory. Phys. Rev. B 2017, 95, 245130. [Google Scholar] [CrossRef] [Green Version]
  24. Kapcia, K.J.; Robaszkiewicz, S.; Capone, M.; Amaricci, A. Doping-Driven Metal-Insulator Transitions and Charge Orderings in the Extended Hubbard Model. Phys. Rev. B 2017, 95, 125112. [Google Scholar] [CrossRef] [Green Version]
  25. Terletska, H.; Chen, T.; Paki, J.; Gull, E. Charge Ordering and Nonlocal Correlations in the Doped Extended Hubbard Model. Phys. Rev. B 2018, 97, 115117. [Google Scholar] [CrossRef] [Green Version]
  26. Freericks, J.K.; Zlatić, V. Exact Dynamical Mean-Field Theory of the Falicov-Kimball Model. Rev. Mod. Phys. 2003, 75, 1333–1382. [Google Scholar] [CrossRef] [Green Version]
  27. Kapcia, K.J.; Lemański, R.; Zygmunt, M.J. Extended Falicov–Kimball Model: Hartree–Fock vs DMFT approach. J. Phys. Condens. Matter 2020, 33, 065602. [Google Scholar] [CrossRef]
  28. Aziz, R.A.; Buck, U.; Jónsson, H.; Ruiz-Suárez, J.; Schmidt, B.; Scoles, G.; Slaman, M.J.; Xu, J. Two- and Three-Body Forces in the Interaction of He Atoms with Xe Overlayers Adsorbed on (0001) Graphite. J. Chem. Phys. 1989, 91, 6477–6493. [Google Scholar] [CrossRef]
  29. Caragiu, M.; Finberg, S. Alkali Metal Adsorption on Graphite: A Review. J. Phys. Condens. Matter 2005, 17, R995–R1024. [Google Scholar] [CrossRef]
  30. Petrović, M.; Lazić, P.; Runte, S.; Michely, T.; Busse, C.; Kralj, M. Moiré-Regulated Self-Assembly of Cesium Adatoms on Epitaxial Graphene. Phys. Rev. B 2017, 96, 085428. [Google Scholar] [CrossRef]
  31. Dimakis, N.; Valdez, D.; Flor, F.A.; Salgado, A.; Adjibi, K.; Vargas, S.; Saenz, J. Density Functional Theory Calculations on Alkali and the Alkaline Ca Atoms Adsorbed on Graphene Monolayers. Appl. Sur. Sci. 2017, 413, 197–208. [Google Scholar] [CrossRef]
  32. Zhour, K.; El Haj Hassan, F.; Fahs, H.; Vaezzadeh, M. Ab Initio Study of the Adsorption of Potassium on B, N, and BN-Doped Graphene Heterostructure. Mater. Today Commun. 2019, 21, 100676. [Google Scholar] [CrossRef]
  33. Huang, Y.C.; Zhao, K.Y.; Liu, Y.; Zhang, X.Y.; Du, H.Y.; Ren, X.W. Investigation on Adsorption of Ar and N2 on α-Al2O3(0001) Surface from First-Principles Calculations. Vacuum 2020, 176, 109344. [Google Scholar] [CrossRef]
  34. Xing, H.; Hu, P.; Li, S.; Zuo, Y.; Han, J.; Hua, X.; Wang, K.; Yang, F.; Feng, P.; Chang, T. Adsorption and Diffusion of Oxygen on Metal Surfaces Studied by First-Principle Study: A Review. J. Mater. Sci. Technol. 2021, 62, 180–194. [Google Scholar] [CrossRef]
  35. Profeta, G.; Ottaviano, L.; Continenza, A.   3 × 3 R30→3×3 Distortion on the C/Si(111) Surface. Phys. Rev. B 2004, 69, 241307. [Google Scholar] [CrossRef]
  36. Tresca, C.; Brun, C.; Bilgeri, T.; Menard, G.; Cherkez, V.; Federicci, R.; Longo, D.; Debontridder, F.; D’angelo, M.; Roditchev, D.; et al. Chiral Spin Texture in the Charge-Density-Wave Phase of the Correlated Metallic Pb/Si(111) Monolayer. Phys. Rev. Lett. 2018, 120, 196402. [Google Scholar] [CrossRef] [Green Version]
  37. Rodríguez, B.C.R.; Santana, J.A. Adsorption and Diffusion of Sulfur on the (111), (100), (110), and (211) Surfaces of FCC Metals: Density Functional Theory Calculations. J. Chem. Phys. 2018, 149, 204701. [Google Scholar] [CrossRef]
  38. Patra, A.; Peng, H.; Sun, J.; Perdew, J.P. Rethinking CO Adsorption on Transition-Metal Surfaces: Effect of Density-Driven Self-Interaction Errors. Phys. Rev. B 2019, 100, 035442. [Google Scholar] [CrossRef] [Green Version]
  39. Menkah, E.S.; Dzade, N.Y.; Tia, R.; Adei, E.; de Leeuw, N.H. Hydrazine Adsorption on Perfect and Defective FCC Nickel (100), (110) and (111) Surfaces: A Dispersion corrected DFT-D2 study. Appl. Sur. Sci. 2019, 480, 1014–1024. [Google Scholar] [CrossRef] [Green Version]
  40. Bretz, M.; Dash, J.G. Ordering Transitions in Helium Monolayers. Phys. Rev. Lett. 1971, 27, 647–650. [Google Scholar] [CrossRef]
  41. Bretz, M.; Dash, J.G.; Hickernell, D.C.; McLean, E.O.; Vilches, O.E. Phases of He3 and He4 Monolayer Films Adsorbed on Basal-Plane Oriented Graphite. Phys. Rev. A 1973, 8, 1589–1615. [Google Scholar] [CrossRef]
  42. Zhou, S.; Wang, Z. Charge and Spin Order on the Triangular Lattice: NaxCoO2 at x = 0.5. Phys. Rev. Lett. 2007, 98, 226402. [Google Scholar] [CrossRef] [Green Version]
  43. Soumyanarayanan, A.; Yee, M.M.; He, Y.; van Wezel, J.; Rahn, D.J.; Rossnagel, K.; Hudson, E.W.; Norman, M.R.; Hoffman, J.E. Quantum Phase Transition from Triangular to Stripe Charge Order in NbSe2. Proc. Natl. Acad. Sci. USA 2013, 110, 1623–1627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ugeda, M.M.; Bradley, A.J.; Zhang, Y.; Onishi, S.; Chen, Y.; Ruan, W.; Ojeda-Aristizabal, C.; Ryu, H.; Edmonds, M.T.; Tsai, H.Z.; et al. Characterization of Collective Ground States in Single-Layer NbSe2. Nat. Phys. 2016, 12, 92–97. [Google Scholar] [CrossRef]
  45. Xi, X.; Wang, Z.; Zhao, W.; Park, J.H.; Law, K.T.; Berger, H.; Forro, L.; Shan, J.; Mak, K.F. Ising Pairing in Superconducting NbSe2 Atomic Layers. Nat. Phys. 2016, 12, 139–143. [Google Scholar] [CrossRef] [Green Version]
  46. Ptok, A.; Głodzik, S.; Domański, T. Yu-Shiba-Rusinov States of Impurities in a Triangular Lattice of NbSe2 with Spin-Orbit Coupling. Phys. Rev. B 2017, 96, 184425. [Google Scholar] [CrossRef] [Green Version]
  47. Lian, C.S.; Si, C.; Duan, W. Unveiling Charge-Density Wave, Superconductivity, and Their Competitive Nature in Two-Dimensional NbSe2. Nano Lett. 2018, 18, 2924–2929. [Google Scholar] [CrossRef]
  48. Kolekar, S.; Bonilla, M.; Ma, Y.; Diaz, H.C.; Batzill, M. Layer- and Substrate-Dependent Charge Density Wave Criticality in 1T-TiSe2. 2D Mater. 2018, 5, 015006. [Google Scholar] [CrossRef]
  49. Ryu, H.; Chen, Y.; Kim, H.; Tsai, H.Z.; Tang, S.; Jiang, J.; Liou, F.; Kahn, S.; Jia, C.; Omrani, A.A.; et al. Persistent Charge-Density-Wave Order in Single-Layer TaSe2. Nano Lett. 2018, 18, 689–694. [Google Scholar] [CrossRef] [Green Version]
  50. Pásztor, A.; Scarfato, A.; Barreteau, C.; Giannini, E.; Renner, C. Dimensional Crossover of the Charge Density Wave Transition in Thin Exfoliated VSe2. 2D Mater. 2017, 4, 041005. [Google Scholar] [CrossRef] [Green Version]
  51. Zhao, J.; Wijayaratne, K.; Butler, A.; Yang, J.; Malliakas, C.D.; Chung, D.Y.; Louca, D.; Kanatzidis, M.G.; van Wezel, J.; Chatterjee, U. Orbital Selectivity Causing Anisotropy and Particle-Hole Asymmetry in the Charge Density Wave Gap of 2H-TaS2. Phys. Rev. B 2017, 96, 125103. [Google Scholar] [CrossRef] [Green Version]
  52. Chhowalla, M.; Shin, H.S.; Eda, G.; Li, L.J.; Loh, K.P.; Zhang, H. The Chemistry of Two-Dimensional Layered Transition Metal Dichalcogenide Nanosheets. Nat. Chem. 2013, 5, 263–275. [Google Scholar] [CrossRef]
  53. Kaneko, R.; Tocchio, L.F.; Valentí, R.; Gros, C. Emergent Lattices with Geometrical Frustration in Doped Extended Hubbard Models. Phys. Rev. B 2016, 94, 195111. [Google Scholar] [CrossRef] [Green Version]
  54. Kaneko, R.; Tocchio, L.F.; Valenti, R.; Becca, F. Charge Orders in Organic Charge-Transfer Salts. New J. Phys. 2017, 19, 103033. [Google Scholar] [CrossRef] [Green Version]
  55. Han, Z.; Kivelson, S.A.; Yao, H. Strong Coupling Limit of the Holstein-Hubbard Model. Phys. Rev. Lett. 2020, 125, 167001. [Google Scholar] [CrossRef] [PubMed]
  56. Müller-Hartmann, E. Correlated Fermions on a Lattice in High Dimensions. Z. Phys. B Condens. Matter 1989, 74, 507–512. [Google Scholar] [CrossRef]
  57. Pearce, P.A.; Thompson, C.J. The Anisotropic Heisenberg Model in the Long-Range Interaction Limit. Commun. Math. Phys. 1975, 41, 191–201. [Google Scholar] [CrossRef]
  58. Pearce, P.A.; Thompson, C.J. The High Density Limit for Lattice Spin Models. Commun. Math. Phys. 1978, 58, 131–138. [Google Scholar] [CrossRef]
  59. Kapcia, K.; Robaszkiewicz, S. The Effects of the Next-Nearest-Neighbour Density-Density Interaction in the Atomic Limit of the Extended Hubbard Model. J. Phys. Condens. Matter 2011, 23, 105601. [Google Scholar] [CrossRef] [Green Version]
  60. Kapcia, K.J.; Robaszkiewicz, S. On the Phase Diagram of the Extended Hubbard Model with Intersite Density-Density Interactions in the Atomic Limit. Phys. A 2016, 461, 487–497. [Google Scholar] [CrossRef]
  61. Borgs, C.; Jedrzejewski, J.; Koteckỳ, R. The Staggered Charge-Order Phase of the Extended Hubbard Model in the Atomic Limit. J. Phys. A Math. Gen. 1996, 29, 733–747. [Google Scholar] [CrossRef] [Green Version]
  62. Fröhlich, J.; Rey-Bellet, L.; Ueltschi, D. Quantum Lattice Models at Intermediate Temperature. Commun. Math. Phys. 2001, 224, 33–63. [Google Scholar] [CrossRef] [Green Version]
  63. Pawłowski, G. Charge Orderings in the Atomic Limit of the Extended Hubbard Model. Eur. Phys. J. B 2006, 53, 471–479. [Google Scholar] [CrossRef]
  64. Ganzenmüller, G.; Pawłowski, G. Flat Histogram Monte Carlo Sampling for Mechanical Variables and Conjugate Thermodynamic Fields with Example Applications to Strongly Correlated Electronic Systems. Phys. Rev. E 2008, 78, 036703. [Google Scholar] [CrossRef] [PubMed]
  65. Jędrzejewski, J. Phase Diagrams of Extended Hubbard Models in the Atomic Limit. Phys. A 1994, 205, 702–717. [Google Scholar] [CrossRef]
  66. Rademaker, L.; Pramudya, Y.; Zaanen, J.; Dobrosavljević, V. Influence of Long-Range Interactions on Charge Ordering Phenomena on a Square Lattice. Phys. Rev. E 2013, 88, 032121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Kapcia, K.J.; Barański, J.; Robaszkiewicz, S.; Ptok, A. Various Charge-Ordered States in the Extended Hubbard Model with On-Site Attraction in the Zero-Bandwidth Limit. J. Supercond. Nov. Magn. 2017, 30, 109–115. [Google Scholar] [CrossRef] [Green Version]
  68. Kapcia, K.J.; Barański, J.; Ptok, A. Diversity of Charge Orderings in Correlated Systems. Phys. Rev. E 2017, 96, 042104. [Google Scholar] [CrossRef] [Green Version]
  69. Lee, S.J.; Lee, J.R.; Kim, B. Patterns of Striped Order in the Classical Lattice Coulomb Gas. Phys. Rev. Lett. 2001, 88, 025701. [Google Scholar] [CrossRef] [Green Version]
  70. Mancini, F.; Mancini, F.P. One-Dimensional Extended Hubbard Model in the Atomic Limit. Phys. Rev. E 2008, 77, 061120. [Google Scholar] [CrossRef] [Green Version]
  71. Mancini, F.; Mancini, F.P. Extended Hubbard Model in the Presence of a Magnetic Field. Eur. Phys. J. B 2009, 68, 341–351. [Google Scholar] [CrossRef] [Green Version]
  72. Mancini, F.; Plekhanov, E.; Sica, G. Exact Solution of the 1D Hubbard Model with NN and NNN Interactions in the Narrow-Band Limit. Eur. Phys. J. B 2013, 86, 408. [Google Scholar] [CrossRef] [Green Version]
  73. Kaneko, R.; Nonomura, Y.; Kohno, M. Thermal Algebraic-Decay Charge Liquid Driven by Competing Short-Range Coulomb Repulsion. Phys. Rev. B 2018, 97, 205125. [Google Scholar] [CrossRef] [Green Version]
  74. Kapcia, K.J. Charge Order of Strongly Bounded Electron Pairs on the Triangular Lattice: The Zero-Bandwidth Limit of the Extended Hubbard Model with Strong Onsite Attraction. J. Supercond. Nov. Magn. 2019, 32, 2751–2757. [Google Scholar] [CrossRef]
  75. Houtappel, R.M.F. Statistics of Two-Dimensional Hexagonal Ferromagnetics with “Ising”-Interaction Between Nearest Neighbours Only. Physica 1950, 16, 391–392. [Google Scholar] [CrossRef]
  76. Houtappel, R.M.F. Order-Disorder in Hexagonal Lattices. Physica 1950, 16, 425–455. [Google Scholar] [CrossRef]
  77. Wannier, G.H. Antiferromagnetism. The Triangular Ising Net. Phys. Rev. 1950, 79, 357–364. [Google Scholar] [CrossRef] [Green Version]
  78. Wannier, G.H. Antiferromagnetism. The Triangular Ising Net (erratum). Phys. Rev. B 1973, 7, 5017. [Google Scholar] [CrossRef]
  79. Schick, M.; Walker, J.S.; Wortis, M. Phase Diagram of the Triangular Ising Model: Renormalization-Group Calculation with Application to Adsorbed Monolayers. Phys. Rev. B 1977, 16, 2205–2219. [Google Scholar] [CrossRef]
  80. Metcalf, B.D. Ground State Spin Orderings of the Triangular Ising Model with the Nearest and Next Nearest Neighbor Interaction. Phys. Lett. A 1974, 46, 325–326. [Google Scholar] [CrossRef]
  81. Oitmaa, J. The Triangular Lattice Ising Model with First and Second Neighbour Interactions. J. Phys. A Math. Gen. 1982, 15, 573–585. [Google Scholar] [CrossRef]
  82. Saito, Y.; Igeta, K. Antiferromagnetic Ising Model on a Triangular Lattice. J. Phys. Soc. Jpn. 1984, 53, 3060–3069. [Google Scholar] [CrossRef]
  83. Tanaka, Y.; Uryû, N. Ground State Spin Configurations of the Triangular Ising Net with the First, Second and Third Nearest Neighbor Interactions. Prog. Theor. Phys. 1976, 55, 1356–1372. [Google Scholar] [CrossRef] [Green Version]
  84. Kudo, T.; Katsura, S. A Method of Determining the Orderings of the Ising Model with Several Neighbor Interactions under the Magnetic Field and Applications to Hexagonal Lattices. Prog. Theor. Phys. 1976, 56, 435–449. [Google Scholar] [CrossRef] [Green Version]
  85. Arrigoni, E.; Strinati, G.C. Doping-Induced Incommensurate Antiferromagnetism in a Mott-Hubbard Insulator. Phys. Rev. B 1991, 44, 7455–7465. [Google Scholar] [CrossRef]
  86. Bąk, M. Mixed Phase and Bound States in the Phase Diagram of the Extended Hubbard Model. Acta Phys. Pol. A 2004, 106, 637–646. [Google Scholar] [CrossRef]
  87. Bursill, R.J.; Thompson, C.J. Variational Bounds for Lattice Fermion Models II. Extended Hubbard Model in the Atomic Limit. J. Phys. A Math. Gen. 1993, 26, 4497–4511. [Google Scholar] [CrossRef]
  88. Kapcia, K.; Robaszkiewicz, S. Stable and Metastable Phases in the Atomic Limit of the Extended Hubbard Model with Intersite Density-Density Interactions. Acta. Phys. Pol. A. 2012, 121, 1029–1031. [Google Scholar] [CrossRef]
  89. Becker, C.; Soltan-Panahi, P.; Kronjäger, J.; Dörscher, S.; Bongs, K.; Sengstock, K. Ultracold Quantum Gases in Triangular Optical Lattices. New J. Phys. 2010, 12, 065025. [Google Scholar] [CrossRef] [Green Version]
  90. Struck, J.; Ölschläger, C.; Le Targat, R.; Soltan-Panahi, P.; Eckardt, A.; Lewenstein, M.; Windpassinger, P.; Sengstock, K. Quantum Simulation of Frustrated Classical Magnetism in Triangular Optical Lattices. Science 2011, 333, 996–999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Cao, Y.; Fatemi, V.; Fang, S.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; Jarillo-Herrero, P. Unconventional Superconductivity in Magic-Angle Graphene Superlattices. Nature 2018, 556, 43–50. [Google Scholar] [CrossRef] [PubMed]
  92. Cao, Y.; Fatemi, V.; Demir, A.; Fang, S.; Tomarken, S.L.; Luo, J.Y.; Sanchez-Yamagishi, J.D.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; et al. Correlated Insulator Behaviour at Half-Filling in Magic-Angle Graphene Superlattices. Nature 2018, 556, 80–84. [Google Scholar] [CrossRef]
  93. Yankowitz, M.; Chen, S.; Polshyn, H.; Zhang, Y.; Watanabe, K.; Taniguchi, T.; Graf, D.; Young, A.F.; Dean, C.R. Tuning Superconductivity in Twisted Bilayer Graphene. Science 2019, 363, 1059–1064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Wu, F.; MacDonald, A.H.; Martin, I. Theory of Phonon-Mediated Superconductivity in Twisted Bilayer Graphene. Phys. Rev. Lett. 2018, 121, 257001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Xu, C.; Balents, L. Topological Superconductivity in Twisted Multilayer Graphene. Phys. Rev. Lett. 2018, 121, 087001. [Google Scholar] [CrossRef] [Green Version]
  96. Lian, B.; Wang, Z.; Bernevig, B.A. Twisted Bilayer Graphene: A Phonon-Driven Superconductor. Phys. Rev. Lett. 2019, 122, 257002. [Google Scholar] [CrossRef] [Green Version]
  97. Wang, G.; Chernikov, A.; Glazov, M.M.; Heinz, T.F.; Marie, X.; Amand, T.; Urbaszek, B. Colloquium: Excitons in Atomically Thin Transition Metal Dichalcogenides. Rev. Mod. Phys. 2018, 90, 021001. [Google Scholar] [CrossRef] [Green Version]
  98. Xu, Y.; Liu, S.; Rhodes, D.A.; Watanabe, K.; Taniguchi, T.; Hone, J.; Elser, V.; Mak, K.F.; Shan, J. Correlated Insulating States at Fractional Fillings of Moiré Superlattices. Nature 2020, 587, 214–218. [Google Scholar] [CrossRef]
  99. Kapcia, K. Metastability and phase separation in a simple model of a superconductor with extremely short coherence length. J. Supercond. Nov. Magn. 2014, 27, 913–917. [Google Scholar] [CrossRef] [Green Version]

Short Biography of Authors

Nanomaterials 11 01181 i001Konrad Jerzy Kapcia is an assistant professor (adjunct) at Condensed Matter Theory Division of Faculty of Physics of Adam Mickiewicz University (AMU) in Poznań (Poland). He graduated in physics (M.Sc.) in 2009 and defended his dissertation (PhD) in theoretical physics in 2014 (both at AMU). Next, he worked in the Institute of Physics of the Polish Academy of Sciences (PAS) in Warsaw and in the Henryk Niewodniczański Institute of Nuclear Physics of PAS in Kraków. He was a visiting scientist at SISSA in Trieste (Italy) and CFEL-DESY in Hamburg (Germany). He is interested in theoretical physics, mainly in the field of strongly correlated systems (including cold gases on optical lattices) and superconductivity as well as in ab initio modeling of real materials (electron and phonon properties). He was the head and principal investigator of several projects financed by National Science Centre (Poland) and a recipient of a few scholarships of Ministry of Science and Higher Education (Poland). He is a member of the Polish Physical Society, the Institute of Physics (IoP, UK), and the American Physical Society.
Figure 1. (a) The schema of the triangular lattice on which the extended Hubbard model in the atomic limit is studied in the present work. The lattice is divided into three equivalent sublattices ( α = A , B , C ) denoted by different symbols. The dash-dotted line denotes the boundaries of 3 × 3 unit cell. By solid and dashed lines all nearest neighbors and all next-nearest neighbors of a chosen site from sublattice A are indicated, respectively. (b) There different types of particle arrangements in 3 × 3 unit cells (i.e., the tri-sublattice assumption) corresponding to NO, DCO, TCO phases (as labeled). Symbol shapes on each panel correspond to respective concentrations at the lattice sites.
Figure 1. (a) The schema of the triangular lattice on which the extended Hubbard model in the atomic limit is studied in the present work. The lattice is divided into three equivalent sublattices ( α = A , B , C ) denoted by different symbols. The dash-dotted line denotes the boundaries of 3 × 3 unit cell. By solid and dashed lines all nearest neighbors and all next-nearest neighbors of a chosen site from sublattice A are indicated, respectively. (b) There different types of particle arrangements in 3 × 3 unit cells (i.e., the tri-sublattice assumption) corresponding to NO, DCO, TCO phases (as labeled). Symbol shapes on each panel correspond to respective concentrations at the lattice sites.
Nanomaterials 11 01181 g001
Figure 2. Ground state phase diagram of the model on the triangular lattice as a function of shifted chemical potential μ ¯ = μ W 1 W 2 for W 1 > 0 and W 2 0 ( | k | = | W 2 | / W 1 ). The regions are labeled by the names of the phases defined in Table 1 and numbers corresponding to concentrations in each sublattice n A , n B and n C .
Figure 2. Ground state phase diagram of the model on the triangular lattice as a function of shifted chemical potential μ ¯ = μ W 1 W 2 for W 1 > 0 and W 2 0 ( | k | = | W 2 | / W 1 ). The regions are labeled by the names of the phases defined in Table 1 and numbers corresponding to concentrations in each sublattice n A , n B and n C .
Nanomaterials 11 01181 g002
Figure 3. Ground state phase diagram of the model as a function of particle concentration n for W 1 > 0 and W 2 0 ( | k | = | W 2 | / W 1 ). The regions are labeled by the names of the homogeneous phases (cf. Table 2). For W 2 = 0 all homogeneous phases are degenerated with macroscopic phase separated states indicated in the last column of Table 2. In regions filled by slantwise pattern the phase separated states occurs at infinitesimally T > 0 for W 2 = 0 . For W 2 < 0 the phase separated states occur inside the regions, whereas at the vertical boundaries for commensurate filling the homogeneous states (defined in Table 1) still exist. The boundary at U / W 1 = 0 (schematically indicated by dashed green line) denotes the boundaries between homogeneous phases, which do not overlap with the boundaries between phase separated states for W 2 < 0 . Squares denote transitions for fixed n between homogeneous phase at commensurate fillings.
Figure 3. Ground state phase diagram of the model as a function of particle concentration n for W 1 > 0 and W 2 0 ( | k | = | W 2 | / W 1 ). The regions are labeled by the names of the homogeneous phases (cf. Table 2). For W 2 = 0 all homogeneous phases are degenerated with macroscopic phase separated states indicated in the last column of Table 2. In regions filled by slantwise pattern the phase separated states occurs at infinitesimally T > 0 for W 2 = 0 . For W 2 < 0 the phase separated states occur inside the regions, whereas at the vertical boundaries for commensurate filling the homogeneous states (defined in Table 1) still exist. The boundary at U / W 1 = 0 (schematically indicated by dashed green line) denotes the boundaries between homogeneous phases, which do not overlap with the boundaries between phase separated states for W 2 < 0 . Squares denote transitions for fixed n between homogeneous phase at commensurate fillings.
Nanomaterials 11 01181 g003
Figure 4. Phase diagrams of the model for U / W 1 = 1.00 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). All transitions are first order and regions of phase separated state (PS 1 :NO/DCO and PS 2 :DCO/DCO) occurrence are present on panel (b). NO and DCO denote homogeneous phases defined in Figure 1b.
Figure 4. Phase diagrams of the model for U / W 1 = 1.00 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). All transitions are first order and regions of phase separated state (PS 1 :NO/DCO and PS 2 :DCO/DCO) occurrence are present on panel (b). NO and DCO denote homogeneous phases defined in Figure 1b.
Nanomaterials 11 01181 g004
Figure 5. Phase diagrams of the model for U / W 1 = 0.75 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). All transitions are first order. Other denotations as in Figure 4.
Figure 5. Phase diagrams of the model for U / W 1 = 0.75 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). All transitions are first order. Other denotations as in Figure 4.
Nanomaterials 11 01181 g005
Figure 6. Phase diagrams of the model for U / W 1 = 0.20 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). The boundary TCO–DCO is second order, the remaining are first order. Other denotations as in Figure 4. The diagrams are shown only for μ ¯ 0 and n 1 , but they are symmetric with respect to μ ¯ = 0 and n = 1 , respectively.
Figure 6. Phase diagrams of the model for U / W 1 = 0.20 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). The boundary TCO–DCO is second order, the remaining are first order. Other denotations as in Figure 4. The diagrams are shown only for μ ¯ 0 and n 1 , but they are symmetric with respect to μ ¯ = 0 and n = 1 , respectively.
Nanomaterials 11 01181 g006
Figure 7. Phase diagrams of the model for U / W 1 = 0.35 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). The boundary TCO–DCO is second order, the remaining are first order. Other denotations as in Figure 4. The diagrams are shown only for μ ¯ 0 and n 1 , but they are symmetric with respect to μ ¯ = 0 and n = 1 , respectively.
Figure 7. Phase diagrams of the model for U / W 1 = 0.35 as a function of (a) chemical potential μ ¯ / W 1 and (b) particle concentration n ( W 1 > 0 , W 2 = 0 ). The boundary TCO–DCO is second order, the remaining are first order. Other denotations as in Figure 4. The diagrams are shown only for μ ¯ 0 and n 1 , but they are symmetric with respect to μ ¯ = 0 and n = 1 , respectively.
Nanomaterials 11 01181 g007
Figure 8. Dependencies of particle concentrations n α ’s in the sublattices (red dotted, blue dashed, and green dot-dashed lines) as a function of k B T / W 1 [(af)] and μ ¯ / W 1 [(gh)] for W 2 = 0 . Black solid lines denote total particle concentration n = ( n A + n B + n C ) / 3 . They are obtained for: (a) U / W 1 = 1.00 , μ ¯ / W 1 = 0 ; (b) U / W 1 = 0.20 , μ ¯ / W 1 = 0 ; (c) U / W 1 = 0.35 , μ ¯ / W 1 = 0 ; (d) U / W 1 = 0.35 , μ ¯ / W 1 = 0.15 ; (e) U / W 1 = 0.35 , μ ¯ / W 1 = 0.5 ; (f) U / W 1 = 0.35 , μ ¯ / W 1 = 0.8 ; (g) U / W 1 = 0.35 , k B T / W 1 = 0.40 ; (h) U / W 1 = 0.35 , k B T / W 1 = 0.80 ; (i) U / W 1 = 0.35 , k B T / W 1 = 0.11 . Vertical solid and dashed lines indicate points of continuous and discontinuous transitions, respectively.
Figure 8. Dependencies of particle concentrations n α ’s in the sublattices (red dotted, blue dashed, and green dot-dashed lines) as a function of k B T / W 1 [(af)] and μ ¯ / W 1 [(gh)] for W 2 = 0 . Black solid lines denote total particle concentration n = ( n A + n B + n C ) / 3 . They are obtained for: (a) U / W 1 = 1.00 , μ ¯ / W 1 = 0 ; (b) U / W 1 = 0.20 , μ ¯ / W 1 = 0 ; (c) U / W 1 = 0.35 , μ ¯ / W 1 = 0 ; (d) U / W 1 = 0.35 , μ ¯ / W 1 = 0.15 ; (e) U / W 1 = 0.35 , μ ¯ / W 1 = 0.5 ; (f) U / W 1 = 0.35 , μ ¯ / W 1 = 0.8 ; (g) U / W 1 = 0.35 , k B T / W 1 = 0.40 ; (h) U / W 1 = 0.35 , k B T / W 1 = 0.80 ; (i) U / W 1 = 0.35 , k B T / W 1 = 0.11 . Vertical solid and dashed lines indicate points of continuous and discontinuous transitions, respectively.
Nanomaterials 11 01181 g008
Figure 9. Phase diagram of the model for half-filling ( μ ¯ = 0 or n = 1 ) as a function of onsite interaction U / W 1 ( W 1 > 0 and W 2 = 0 ). The order-disorder boundary for U / W 1 < ( 1 / 3 ) ln ( 2 ) is a line consisting of some higher-order critical points as discussed in the text.
Figure 9. Phase diagram of the model for half-filling ( μ ¯ = 0 or n = 1 ) as a function of onsite interaction U / W 1 ( W 1 > 0 and W 2 = 0 ). The order-disorder boundary for U / W 1 < ( 1 / 3 ) ln ( 2 ) is a line consisting of some higher-order critical points as discussed in the text.
Nanomaterials 11 01181 g009
Table 1. Homogeneous phases ( z n , n = 1 , 2 ) or 3 × 3 unit cells (the triangular lattice) at T = 0 (for fixed μ ). Star “*” in superscript indicates that the phase is obtained by the particle-hole transformation (i.e., n α 2 n α ; the NO 1 and TCO phases are invariant under this transformation). In the brackets also an alternative name is given. The degeneration d c × d s of the unit cells (equal to the degeneration of the ground state for z n limit) and degeneration D c × D s of the ground state phases constructed from the corresponding unit cells for the triangular lattice is given (with respect to charge and spin degrees of freedom).
Table 1. Homogeneous phases ( z n , n = 1 , 2 ) or 3 × 3 unit cells (the triangular lattice) at T = 0 (for fixed μ ). Star “*” in superscript indicates that the phase is obtained by the particle-hole transformation (i.e., n α 2 n α ; the NO 1 and TCO phases are invariant under this transformation). In the brackets also an alternative name is given. The degeneration d c × d s of the unit cells (equal to the degeneration of the ground state for z n limit) and degeneration D c × D s of the ground state phases constructed from the corresponding unit cells for the triangular lattice is given (with respect to charge and spin degrees of freedom).
Phase n A n B n C d c × d s D c × D s ω 0
NO 0 (NO 2 * )000 1 × 1 1 × 1 0
NO 1 (NO 1 * )111 1 × 8 1 × 2 L ( 2 μ + W 1 + W 2 ) / 2
NO 2 (NO 0 * )222 1 × 1 1 × 1 2 μ + U + 2 W 1 + 2 W 2
DCO 1 001 3 × 2 3 × 2 L / 3 ( 2 μ + W 2 ) / 6
DCO 1 * 122 3 × 2 3 × 2 L / 3 ( 10 μ + 4 U + 8 W 1 + 9 W 2 ) / 6
DCO 2 002 3 × 1 3 × 1 ( 2 μ + U + 2 W 2 ) / 3
DCO 2 * 022 3 × 1 3 × 1 ( 4 μ + 2 U + 2 W 1 + 4 W 2 ) / 3
DCO 3 011 3 × 4 3 × 4 L / 3 ( 4 μ + W 1 + 2 W 2 ) / 6
DCO 3 * 112 3 × 4 3 × 4 L / 3 ( 8 μ + 2 U + 5 W 1 + 6 W 2 ) / 6
TCO (TCO * )012 6 × 2 6 × 2 L / 3 ( 6 μ + 2 U + 2 W 1 + 5 W 2 ) / 6
Table 2. Homogeneous phases at T = 0 (for fixed n) defined by n α ’s and D o c c . n s and n f define the range [ n s , n f ] of n, where the phase is correctly defined. In the last column, the phase separated state degenerated with the homogeneous phase in range ( n s , n f ) for W 2 = 0 is mentioned. Star “*” in superscript indicates that the phase is obtained by the particle-hole transformation (i.e., n α 2 n α ; TCO A , TCO A * , TCO B , and TCO B * phases are invariant under this transformation).
Table 2. Homogeneous phases at T = 0 (for fixed n) defined by n α ’s and D o c c . n s and n f define the range [ n s , n f ] of n, where the phase is correctly defined. In the last column, the phase separated state degenerated with the homogeneous phase in range ( n s , n f ) for W 2 = 0 is mentioned. Star “*” in superscript indicates that the phase is obtained by the particle-hole transformation (i.e., n α 2 n α ; TCO A , TCO A * , TCO B , and TCO B * phases are invariant under this transformation).
Phase n A n B n C D occ n s n f PS
DCO A 00 3 n n / 2 0 2 / 3 NO 0 /DCO 2
DCO B 00 3 n 00 1 / 3 NO 0 /DCO 1
DCO C 00 3 n n 1 / 3 1 / 3 2 / 3 DCO 1 /DCO 2
DCO D 3 n 2 110 2 / 3 1DCO 3 /NO 1
TCO A 0 3 n 2 2 n / 2 2 / 3 4 / 3 DCO 2 /DCO 2 *
TCO B 0 3 n 2 2 1 / 3 2 / 3 1DCO 2 /TCO
TCO C 0 3 n 1 10 1 / 3 2 / 3 DCO 1 /DCO 3
DCO A * 3 n 4 22 n / 2 4 / 3 2DCO 2 * /NO 2
DCO B * 3 n 4 22 n 1 5 / 3 2DCO 1 * /NO 2
DCO C * 3 n 4 22 2 / 3 4 / 3 5 / 3 DCO 2 * /DCO 1 *
DCO D * 11 3 n 2 n 1 1 4 / 3 NO 1 /DCO 3 *
TCO A * 0 3 n 2 2 n / 2 2 / 3 4 / 3 DCO 2 /DCO 2 *
TCO B * 0 3 n 2 2 n 2 / 3 1 4 / 3 TCO * /DCO 2 *
TCO C * 1 3 n 3 2 n 1 4 / 3 5 / 3 DCO 3 * /DCO 1 *
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kapcia, K.J. Charge-Order on the Triangular Lattice: A Mean-Field Study for the Lattice S = 1/2 Fermionic Gas. Nanomaterials 2021, 11, 1181. https://doi.org/10.3390/nano11051181

AMA Style

Kapcia KJ. Charge-Order on the Triangular Lattice: A Mean-Field Study for the Lattice S = 1/2 Fermionic Gas. Nanomaterials. 2021; 11(5):1181. https://doi.org/10.3390/nano11051181

Chicago/Turabian Style

Kapcia, Konrad Jerzy. 2021. "Charge-Order on the Triangular Lattice: A Mean-Field Study for the Lattice S = 1/2 Fermionic Gas" Nanomaterials 11, no. 5: 1181. https://doi.org/10.3390/nano11051181

APA Style

Kapcia, K. J. (2021). Charge-Order on the Triangular Lattice: A Mean-Field Study for the Lattice S = 1/2 Fermionic Gas. Nanomaterials, 11(5), 1181. https://doi.org/10.3390/nano11051181

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop