Chiral p-wave pairing of ultracold fermionic atoms due to a quadratic band touching

We study the superfuild ground state of ultracold fermions in optical lattices with a quadratic band touching. Examples are a checkerboard lattice around half filling and a kagome lattice above one third filling. Instead of pairing between spin states, here we focus on pairing interactions between different orbital states. We find that our systems have only odd-parity(orbital) pairing instability while the singlet(orbital) pairing instability vanishes thanks to the quadratic band touching. In the mean field level, the ground state is found to be a chiral p-wave pairing superfluid(mixed with finite f-wave pairing order-parameters) which supports Majorana fermions.     

Quantum phase transition of cold atoms trapped in optical lattices

We review our recent theoretical advances in phase transition of cold atoms in optical lattices, such as triangular lattice, honeycomb lattice, and Kagomé lattice. By employing the new developed numerical methods called dynamical cluster approximation and cellular dynamical mean-field theory, the properties in different phases of cold atoms in optical lattices are studied, such as density of states, Fermi surface and double occupancy. On triangular lattice, a reentrant behavior of phase translation line between Fermi liquid state and pseudogap state is found due to the Kondo effect. We find the system undergoes a second order Mott transition from a metallic state into a Mott insulator state on honeycomb lattice and triangular Kagomé lattice. The stability of quantum spin Hall phase towards interaction on honeycomb lattice with spin-orbital coupling is systematically discussed. And we investigate the transition from quantum spin Hall insulator to normal insulator in Kagomé lattice which includes a nearest-neighbor intrinsic spin-orbit coupling and a trimerized Hamiltonian. In addition, we propose the experimental protocols to observe these phase transition of cold atoms in optical lattices.     

Ferromagnetism in the Hubbard model on the square lattice: Improved instability criterion for the Nagaoka state

The instability of the fully polarized ferromagnetic state (Nagaoka state) with respect to single spin flips is re-examined for the Hubbard model on the square lattice with a large family of variational wave functions which include correlation effects of the majority spins in the vicinity of the flipped spin. We find a critical hole density of δ_cr = 0.251 for U = ∞ and a critical coupling of U_cr = 77.7t. Both values improve previous variational results considerably.     

Comparison of diluted antiferromagnetic Ising models on frustrated lattices in a magnetic field

We study diluted antiferromagnetic Ising models on triangular and kagome lattices in a magnetic field, using the replica-exchange Monte Carlo method. We observe seven and five plateaus in the magnetization curve of the diluted antiferromagnetic Ising model on the triangular and kagome lattices, respectively, when a magnetic field is applied. These observations contrast with the two plateaus observed in the pure model. The origin of multiple plateaus is investigated by considering the spin configuration of triangles in the diluted models. We compare these results with those of a diluted antiferromagnetic Ising model on the three-dimensional pyrochlore lattice in a magnetic field pointing in the [111] direction, sometimes referred to as the “kagome-ice” problem. We discuss the similarity and dissimilarity of the magnetization curves of the “kagome-ice” state and the two-dimensional kagome lattice.     

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.     

Orbital order in Mott insulators of spinless p-band fermions

A gas of strongly interacting single-species (spinless) p-orbital fermionic atoms in 2D optical lattices is proposed and studied. Several interesting new features are found. In the Mott limit on a square lattice, the gas is found to be described effectively by an orbital exchange Hamiltonian equivalent to a pseudospin-1/2 XXZ model. For a triangular, honeycomb, or kagome lattice, the orbital exchange is geometrically frustrated and described by a new quantum 120 degrees model. We determine the orbital ordering on the kagome lattice, and show how orbital wave fluctuations select ground states via the order by disorder mechanism for the honeycomb lattice. We discuss experimental signatures of various orbital ordering.     

Ferromagnetism in the Hubbard model: instability of the Nagaoka state on the square lattice

The instability of the fully polarized ferromagnetic ground state (Nagaoka state) of the Hubbard model on the square lattice is investigated. We use single spin flip variational wave functions including majority spin correlation effects and calculate spin flip energies in the thermodynamic limit. With very local wave functions and with moderate numbers of variational parameters we reproduce the best known estimate for the critical hole density δ_cr = 0.29 and we obtain an estimate of U_cr = 63 t for the critical coupling which is considerably better than the best estimate of U_cr = 42 t previously known. The simplicity of our wave functions makes the physical origin of the various aspects of the instability particularly transparent.     

U(1) gauge theory of the Hubbard model: spin liquid states and possible application to kappa-(BEDT-TTF)2Cu2(CN)3

We formulate a U(1) gauge theory of the Hubbard model in the slave-rotor representation. From this formalism it is argued that spin liquid phases may exist near the Mott transition in the Hubbard model on triangular and honeycomb lattices at half filling. The organic compound kappa-(BEDT-TTF)2Cu2(CN)3 is a good candidate for the spin liquid state on a triangular lattice. We predict a highly unusual temperature dependence for the thermal conductivity of this material.     

Ferromagnetism in the Hubbard model: Spin waves and instability of the Nagaoka state

We discuss two single spin flip variational wave functions describing spin wave excitations which were proposed earlier by Shastry, Krishnamurthy and Anderson (SKA) and by Basile and Elser (BE), respectively, in order to investigate the instability of the fully polarized ferromagnetic state (Nagaoka state) in the infinite U Hubbard model. We calculate the energy of these variational states for the square lattice and for multiple chains. At the zone boundary in the vicinity of the point (0, π) the spin wave energy is reduced substantially by the binding of the spin up hole to the flipped down spin. For the square lattice this leads to a critical hole density of δ_cr = 0.407 for the SKA spin wave and of δ_cr = 0.322 for the BE spin wave which implies remarkable improvements in comparison to the corresponding scattering states investigated previously.     

Magnetic and non-magnetic ground states of the Kondo lattice

We use the variational method to investigate the ground state phase diagram of the Kondo lattice Hamiltonian for arbitraryJ/W, and conduction electron concentrationn _c (J is the Kondo coupling andW the bandwidth). We are particularly interested in the question under which circumstances the globally singlet (collective Kondo) Fermi liquid type ground state becomes unstable against magnetic ordering. For the collective Kondo singlet we use the lattice generalization of Yosida's wavefunction which implies the existence of a large Fermi volume, in accordance with Luttinger's theorem. Using the Gutzwiller approximation, we derive closed-form results for the ground state energy at arbitraryJ/W andn _c, and for the Kondo gap atn _c=1. We introduce simple trial states to describe ferromagnetic, antiferromagnetic, and spiral ordering in the small-J (RKKY) regime, and Nagaoka type ferromagnetism at largeJ/W. We study three particular cases: a band with a constant density of states, and the (tight binding) linear chain, and square lattice periodic Kondo models. We find that the lattice enhancement of the Kondo effect, which is described in our theory of the Fermi liquid state, pushes the RKKY-to-nonmagnetic phase boundary to much smaller values ofJ/W than it was previously thought. In our study of the square lattice case, we also find a region of itinerant, Nagaoka-type ferromagnetism at largeJ/W forn _c 1/3.     

Degenerate ferrimagnetic ground state of frustrated kagome lattice

The frustrated Ising model on kagome lattice with nearest-neighboring antiferromagnetic interaction is investigated by using Monte Carlo simulation of the Wang-Landau algorithm and Glauber dynamics. The geometrical frustration leads to a particularly high degeneracy of ground states in this system. A small magnetic field applied can lift the degeneracy partially, and produce the magnetization plateau of 1/3 saturate value (Ms), which is analogous to the magnetic behavior in triangular antiferromagnetic system. However, different from the long-range ferrimagnetic state responsible for 1/3 Ms plateau in triangular lattice, the ferrimagnetic ground state corresponding to 1/3 Ms plateau in kagome lattice is short-ranged and still highly degenerate. Furthermore, the spin configuration of these degenerate ferrimagnetic ground states show an inherent characteristic that the spins along the magnetic field must be aligned on the closed loops, which can be well understood in terms of geometrical frustration.     

Orbital ordering and frustration of p-band Mott insulators

We investigate the general structure of orbital exchange physics in Mott-insulating states of p-orbital systems in optical lattices. Orbital orders occur in both the triangular and kagome lattices. In contrast, orbital exchange in the honeycomb lattice is frustrated as described by a novel quantum 120 degrees model. Its classical ground states are mapped into configurations of the fully packed loop model with an extra U(1) rotation degree of freedom. Quantum orbital fluctuations select a six-site plaquette ground state ordering pattern in the semiclassical limit from the "order from disorder" mechanism. This effect arises from the appearance of a zero energy flat band of orbital excitations.     

Superconducting states of pure and doped graphene

We study the superconducting phases of the two-dimensional honeycomb lattice of graphene. We find two spin singlet pairing states; s wave and an exotic p+ip that is possible because of the special structure of the honeycomb lattice. At half filling, the p+ip phase is gapless and superconductivity is a hidden order. We discuss the possibility of a superconducting state in metal coated graphene.     

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.     

Partially disordered antiferromagnetic state in two-dimensional triangular-lattice antiferromagnet

The partially disordered antiferromagnetic (PDA) state, as an exotic phase peculiar to the antiferromagnet with Ising spin in triangular lattice, is investigated by using Monte Carlo simulations of Wang-Landau algorithm and Glauber algorithm. It is revealed that PDA state, as the ground state of the triangular antiferromagnetic system, presents the complicated spin configuration due to geometrical frustration. The formation of multi-domain structure within the framework of honeycomb antiferromagnet results in the high degeneracy of PDA state. And this degeneracy of ground state can be lifted by a small magnetic field. Consequently the system shows the ferrimagnetic state, and the magnetization plateau of 1/3 saturate value (Ms) is observed in many experiments. Moreover, due to the multi-domain structure of PDA state and those spins on domain walls, the metastable steps may manifest themselves superposed on the 1/3Ms plateau in some special cases.     

Vortex lattice transitions in cyclic spinor condensates

We study the energetics of vortices and vortex lattices produced by rotation in the cyclic phase of F=2 spinor condensates. In addition to the familiar triangular lattice predicted by Tkachenko for 4He, many more complex lattices appear in this system as a result of the spin degree of freedom. In particular, we predict a magnetic-field-driven transition from a triangular lattice to a honeycomb lattice. Other transitions and lattice geometries are driven at constant field by changes in the temperature-dependent ratio of charge and spin stiffnesses, including a transition through an aperiodic vortex structure. Finally, we compute the renormalization of the ratio of the spin and charge stiffnesses from thermal fluctuations using a nonlinear sigma model analysis.     

Supersolid hard-core bosons on the triangular lattice

We determine the phase diagram of hard-core bosons on a triangular lattice with nearest-neighbor repulsion, paying special attention to the stability of the supersolid phase. Similar to the same model on a square lattice we find that for densities rho<1/3 or rho>2/3 a supersolid phase is unstable and the transition between a commensurate solid and the superfluid is of first order. At intermediate fillings 1/3相似文献   

Diverse solid and supersolid phases of bosons in a triangular lattice

We investigate the ground state of bosons with long-range interactions in the large U limit on a triangular lattice. By mapping this system to the spin-1/2 XXZ model in a magnetic field, we can apply the spin wave theory to this study. We demonstrate how to construct the phase diagrams within the spin wave theory. The phase diagrams are given in an extensive parameter region, where, besides the superfluid phase, diverse solid and supersolid phases are shown to exist in this model. Especially, we find that the phase diagram obtained in this method is consistent with the one obtained previously using numerical techniques in the Ising limit. This confirms the effectiveness of our method. We analyze the stability of all the obtained supersolids and show that they will not be ruined by the quantum fluctuations. We observe that the quantum fluctuations in the stripe supersolid phase could be enhanced by the external field. We also discuss the relevance of our result with the experiment that may be realized with ultracold bosonic polar molecules in a triangular optical lattice.