Dimers on the kagome lattice I: Finite lattices

We report exact results on the enumeration of close-packed dimers on a finite kagome lattice with general asymmetric dimer weights under periodic and cylindrical boundary conditions. For symmetric dimer weights, the resulting dimer generating functions reduce to very simple expressions, and we show how the simple expressions can be obtained from the consideration of a spin-variable mapping.     

Floquet topological phase transitions and chiral edge states in a kagome lattice

The Floquet topological phases and chiral edge states in a kagome lattice under a circularly-polarized driving field are studied. In the off-resonant case, the system exhibits the similar character as the kagome lattice model with staggered magnetic fluxes, but the total band width is damped in oscillation. In the on-resonant case, the degeneracy splitting at the Γ point does not always result in a gap. The positions of the other two gaps are influenced by the flat band. With the field intensity increased, these two gaps undergo closing-then-reopening processes, accompanied with the changing of the winding numbers.     

Ferromagnetism in the Hubbard model: instability of the Nagaoka state on the triangular,honeycomb and kagome lattices

In order to analyse the lattice dependence of ferromagnetism in the two-dimensional Hubbard model we investigate the instability of the fully polarised ferromagnetic ground state (Nagaoka state) on the triangular, honeycomb and kagome lattices. We mainly focus on the local instability, applying single spin flip variational wave functions which include majority spin correlation effects. The question of global instability and phase separation is addressed in the framework of Hartree-Fock theory. We find a strong tendency towards Nagaoka ferromagnetism on the non-bipartite lattices (triangular, kagome) for more than half filling. For the triangular lattice we find the Nagaoka state to be unstable above a critical density of n = 1.887 at U = ∞, thereby significantly improving former variational results. For the kagome lattice the region where ferromagnetism prevails in the phase diagram widely exceeds the flat band regime. Our results even allow the stability of the Nagaoka state in a small region below half filling. In the case of the bipartite honeycomb lattice several disconnected regions are left for a possible Nagaoka ground state.     

Dimer statistics of honeycomb lattices on Klein bottle,Möbius strip and cylinder

Dimer statistics is a central problem in statistical physics. In this paper the enumerations of close-packed dimers of honeycomb lattices on Klein bottle, Möbius strip and cylinder are considered. By establishing a Pfaffian orientation or a crossing orientation, and then computing the determinants of the skew-symmetric matrices of the resulting orientation graphs, we obtain explicit expressions of the number of close-packed dimers of the Klein-bottle polyhex, the Möbius polyhex and the cylindrical polyhex.     

Influence of dilution on magnetization properties of geometrically frustrated magnetic systems: Effective-field theory cluster approximations on kagome lattice

The influence of the site dilution on the magnetization properties of the antiferromagnetic spin-1/2 Ising model on the kagome lattice is systematically investigated using various n-site cluster effective-field theory approximations up to the cluster size $n=12$. It is shown that, regardless of the cluster approximation, the site dilution of the system leads, in addition to the existence of the standard saturated ground state, to the formation of four nontrivial plateau ground states together with five single-point ground states that separate them in the zero temperature limit. It is also shown that, while magnetization properties of the saturated ground state and two single-point ground states are stable with respect to the used cluster approximation, the magnetization properties of all four nontrivial plateau ground states and three single-point ground states that separate them strongly depend on the used approximation and have tendency to approach each other with increasing of the cluster approximation.     

Interaction and spin–orbit effects on a kagome lattice at 1/3 filling

We investigate the competing effects of spin-orbit coupling and electron--electron interaction on a kagome lattice at 1/3 filling. We apply the cellular dynamical mean-field theory and its real-space extension combined with the continuous time quantum Monte Carlo method, and obtain a phase diagram including the effects of the interaction and the spin-orbit coupling at T = 0. 1t, where T is the temperature and t is the hopping energy. We find that without the spin-orbit coupling, the system is in a semi-metal phase stable against the electron--electron interaction. The presence of the spin-orbit coupling can induce a topological non-trivial gap and drive the system to a topological insulator, and as the interaction increases, a larger spin--orbit coupling is required to reach the topological insulating phase.     

Calculation of critical properties for the anisotropic two-layer Ising model on the Kagome lattice: Cellular automata approach

The critical point (K_c) of the two-layer Ising model on the Kagome lattice has been calculated with a high precision, using the probabilistic cellular automata with the Glauber algorithm. The critical point is calculated for different values of the inter- and intra-layer couplings (K₁≠K₂≠K₃≠K_z), where K₁, K₂ and K₃ are the nearest-neighbor interactions within each layer in the 1, 2 and 3 directions, respectively, and K_z is the intralayer coupling. A general ansatz equation for the critical point is given as a function of the inter- and intra-layer interactions, ξ=K₃/K₁,σ=K₂/K₁ and ω=K_z/K₁ for the one- and two-layer Ising models on the Kagome lattice.     

Magnetization and Lyapunov exponents on a kagome chain with multi-site exchange interaction

The Ising approximation of the Heisenberg model in a strong magnetic field, with two, three and six spin exchange interactions is studied on a kagome chain. The kagome chain can be considered as an approximation of the third layer of ³He absorbed on the surface of graphite (kagome lattice). By using dynamical approach we have found one- and multi-dimensional mappings (recursion relations) for the partition function. The magnetization diagrams are plotted and they show that the kagome chain is separating into four sublattices with different magnetizations. Magnetization curves of two sublattices exhibit plateaus at zero and 2/3 of the saturation field. The maximal Lyapunov exponent for multi-dimensional mapping is considered and it is shown that near the magnetization plateaus the maximal Lyapunov exponent also exhibits plateaus.     

A hyperbolic function approach to constructing exact solitary wave solutions of the Hybrid lattice and discrete mKdV lattice

Some new exact solitary wave solutions of the Hybrid lattice and discrete mKdV lattice are obtained by using a hyperbolic function approach. This approach can also be applied to other nonlinear differential-difference equations.     

Green's function approach to the impurity state in the Lieb lattice

Wavefunction of an impurity attached to the Lieb lattice is considered by directly computing the Green's function (GF) for the nearest-neighbor tight-binding model. By replacing the lattice GF with a GF array including all the nine GFs defining on the three-atom unit cell of the Lieb lattice, an accurate and efficient numerical technique is developed. Agreement of both the real and imaginary components of the GF between numerical simulation of the lattice GF and continuum-space Fourier transformation is achieved in the whole resonant-energy range. Both results demonstrate that the wavefunction amplitude decays in a power-law pattern while the resonant energy is small and it decays in a pattern stronger than the power law and weaker than exponentially while the resonant energy is large. The decaying exponents depend on the adatom type and location in the unit cell, which directly modifies the temperature-dependent conductance under the variable-range hopping theory.     

Hard-core bosons on the kagome lattice: valence-bond solids and their quantum melting

Using large scale quantum Monte Carlo simulations and dual vortex theory, we analyze the ground state phase diagram of hard-core bosons on the kagome lattice with nearest-neighbor repulsion. In contrast with the case of a triangular lattice, no supersolid emerges for strong interactions. While a uniform superfluid prevails at half filling, two novel solid phases emerge at densities rho=1/3 and rho=2/3. These solids exhibit an only partial ordering of the bosonic density, allowing for local resonances on a subset of hexagons of the kagome lattice. We provide evidence for a weakly first-order phase transition at the quantum melting point between these solid phases and the superfluid.     

Quantum dimer model on the kagome lattice: solvable dimer-liquid and ising gauge theory

We introduce quantum dimer models on lattices made of corner-sharing triangles. These lattices include the kagome lattice and can be defined in arbitrary geometry. They realize fully disordered and gapped dimer-liquid phase with topological degeneracy and deconfined fractional excitations, as well as solid phases. Using geometrical properties of the lattice, several results are obtained exactly, including the full spectrum of a dimer liquid. These models offer a very natural-and maybe the simplest possible-framework to illustrate general concepts such as fractionalization, topological order, and relation to Z2 gauge theories.     

On the generalized Hermite-based lattice Boltzmann construction,lattice sets,weights, moments,distribution functions and high-order models

The influence of the use of the generalized Hermite polynomial on the Hermite-based lattice Bofiz- mann （LB） construction approach, lattice sets, the thermal weights, moments and the equilibrium distribution function （EDF） are addressed. A new moment system is proposed. The theoretical possibility to obtain a unique high-order Hermite-based singel relaxation time LB model capable to exactly inatch some first hydrodynamic inoments thermally i） on-Cartesian lattice, ii） with thermal weights in the EDF, iii） whilst the highest possible hydrodynamic moments that are exactly reatched are obtained with the shortest on-Cartesian lattiee sets with some fixed real-valued temperatures, is also analyzed.     

Correlation functions of odd numbers of spins with finite separations on the Onsager-Ising lattice

A method is developed for making exact evaluations of correlation functions of odd numbers of spins on the Onsager-Ising lattice, applicable to cases in which the separations between the spins are finite. The method is based on an identity which permits the reduction of determinants of infinite-dimensional matrices to those of finite dimension. Particularly simple results are obtained when all spins are on a straight line. Numerical calculations are carried out for a few cases.     

Boltzmann equation on a lattice: Existence and uniqueness of solutions

Cercignani, Greenberg, and Zweifel proved the existence and uniqueness of solutions of the Boltzmann equation on a toroidal lattice under the assumption that the collision kernel is bounded. We give an alternative, considerably simpler, proof which is based on a fixed point argument.     

Diffusion in Lorentz lattice gas cellular automata: The honeycomb and quasi-lattices compared with the square and triangular lattices

We study numerically the nature of the diffusion process on a honeycomb and a quasi-lattice, where a point particle, moving along the bonds of the lattice, scatters from randomly placed scatterers on the lattice sites according to strictly deterministic rules. For the honeycomb lattice fully occupied by fixed rotators two (symmetric) isolated critical points appear to be present, with the same hyperscaling relation as for the square and the triangular lattices. No such points appear to exist for the quasi-lattice. A comprehensive comparison is made with the behavior on the previously studied square and triangular lattices. A great variety of diffusive behavior is found, ranging from propagation, superdiffusion, normal, quasi-normal, and anomalous, to absence of diffusion. The influence of the scattering rules as well as of the lattice structure on the diffusive behavior of a point particle moving on the all lattices studied so far is summarized.     

The effect of a constant number of particles on the pair correlation function and the density profile in a one-dimensional lattice gas

A one-dimensional lattice gas (Ising model) of lengthL and with nearest-neighbor couplingJ is considered in a canonical ensemble with fixed number of particlesN=L/2. Exact expressions and asymptotic forms for largeL are derived for the density-density correlation function, using periodic boundary conditions, and for the density (magnetization) profile, using antisymmetric boundary conditions. The density-density correlation function,g, assumes for temperaturesT> T, withT = 2J(_BlnL)^–1 and forL large, the formg(x) =g ^gc(x) +BL ^–1 +a(x)L ^–1 +O(L^–2) wherex is a distance between considered lattice sites,B is known from earlier work of Lebowitz and Percus,(1b) anda(x) decays exponentially forx . For TT, the correlation function and the density profile behave differently, the latter exhibiting a step in the middle of the interface.     

Spin fluctuations and unconventional pairing on the Lieb lattice

The spin fluctuations and superconducting pairing symmetries in the dispersive band of Lieb lattice are studied by fluctuation exchange approximation. The antiferromagnetic spin density wave is found to exist on the A sublattice (the lattice sites with four nearest neighbors) at half filling. When slightly doped away from half filling, a balance between the combined effects of the (π,π$\pi ,\pi$) and (0.4π,0$0.4\pi ,0$) spin fluctuations and the gaining of the condensation energy leads to the nearly degenerate _dx2−y2$d_{x^2-y^2}$- and _gxy(x2−y2)$g_{xy(x^2-y^2)}$-wave pairing states. After further doped, the _dxy$d_{xy}$-wave state is favored via the intra-sublattice spin fluctuations with a wave vector (π,0$\pi ,0$). We emphasize that the sublattices' contribution and the renormalization of the spectral function play a crucial role on the spin fluctuations and the pairing symmetry. The effect of the imbalance of the on-site energy at different sublattices is also discussed.     

Bethe lattice spin glass: The effects of a ferromagnetic bias and external fields. II. Magnetized spin-glass phase and the de Almeida-Thouless line

In this and the companion paper, we analyze the ±J Ising spin-glass model on the Bethe lattice with fixed uncorrelated boundary conditions. Phase diagrams are derived as a function of temperature vs. concentration of ferromagnetic bonds and, for a symmetric distribution of bonds, external field vs. temperature. In this part we characterize magnetized spin-glass (MSG) phases by divergence of an appropriate susceptibility: at zero field this signals the existence of an intermediate MSG phase; at nonzero field, this is used to identify the de Almeida-Thouless line.     

Explicit relations for the radial distribution functions for one-dimensional lattice (arbitrary spacing) fluids and solutions

It is pointed out that the size of the matrix required to formulate the grand partition function for a one-dimensional lattice fluid for a fixed and finite range of the interatomic potential varies linearly with the density of lattice points used and hence is much smaller and more manageable than the expected size (which varies exponentially with the same quantity) and thus allows very fine grids to be examined. Using the matrix treatment of the grand partition function, it is shown that the radial distribution function for a one-dimensional fluid or solution can be formulated as an explicit matrix product which is simply performed by computer. The resulting distribution functions (which can be extrapolated to the continuum by varying the lattice spacing) are useful as starting solutions for the iterative solution of integral equations for three-dimensional fluids.