Superconductivity in the Two-Dimensional Hubbard Model

It is shown that the existence of d-wave superconductivity in the two-dimensional Hubbard model close to half-filling can be inferred from a renormalization group analysis at one-loop level.     

Unconventional Landau Levels of Ultracold Fermionic Atoms on Two-Dimensional Honeycomb Lattice

We investigate the unconventional Landau levels of ultracold fermionic atoms on the two-dimensional honeycomb optical lattice subjected to an effective magnetic field, which is created with optical means. In the presence of the effective magnetic field, the energy spectrum of the unconventional Landau levels is calculated. Furthermore, we propose to detect the unconventional Landau levels with Bragg scattering techniques.     

Electrical Control of Edge Magnetism in Two-Dimensional Buckled Honeycomb Lattice

We theoretically study indirect spin coupling strength between two magnetic impurities located on honeycomb Kane-Mele zigzag ribbon （KMZR） with periodic and open boundary （PB and OB）. We show that spin interaction J in PB ribbons displays an AFM-FM oscillating behavior with increasing the staggered potential and electron density, and approaches to maximum at the edges. While the spin coupling in OB KMZR shows a trivial smooth AFM coupling with varying staggered potential. Such a novel J（ A） behavior is the combining effect of finite size, topological edge states and inversion symmetry breaking induced by the staggered potential. We propose that one could control the edge magnetism electrically in two-dimensional buckled honeycomb materials, e.g., silicence, germanene and stanene.     

Finite-Temperature Optical Conductivity of the Two-Dimensional Hubbard Model

The finite-temperature optical conductivity of the two-dimensional Hubbard model is studied by means of the grand canonical quantum Monte Carlo method. Results are reported for the half-filled and nonhalf-filled band. The formation and variation of Hubbard gap are observed. For the half-filled band case, our results support that there exists a high-temperature insulator-metal transition.     

Effective-Field Theory for Kinetic Ising Model on Honeycomb Lattice

As an analytical method, the effective-field theory （EFT） is used to study the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field. The effective-field equations of motion of the average magnetization are given for the honeycomb lattice （Z = 3）. The Liapunov exponent A is calculated for discussing the stability of the magnetization and it is used to determine the phase boundary. In the field amplitude ho / ZJ-temperature T/ZJ plane, the phase boundary separating the dynamic ordered and the disordered phase has been drawn. In contrast to previous analytical results that predicted a tricritical point separating a dynamic phase boundary line of continuous and discontinuous transitions, we find that the transition is always continuous. There is inconsistency between our results and previous analytical results, because they do not introduce sufficiently strong fluctuations.     

A Monte Carlo Study of the Half-Filled Two-Dimensional Hubbard Model

Using the path integral representation,bulk thermodynamic properties of the two-dimensional Hubbard model on a square lattice are studied numerically using the hybrid Monte Carlo simulation algorithm.Results for averaged energy,spin-spin correlation functions and Cooper pair correlation functions are presented.All these results suggest that the system exhibits an anti-ferromagnetic correlation at half-filling for a wide range of temperature.``     

A Monte Carlo Study of the Half-Filled Two-Dimensional Hubbard Model

Using the path integral representation,bulk thermodynamic properties of the two-dimensional Hubbard model on a square lattice are studied numerically using the hybrid Monte Carlo simulation algorithm.Results for averaged energy,spin-spin correlation functions and Cooper pair correlation functions are presented.All these results suggest that the system exhibits an anti-ferromagnetic correlation at half-filling for a wide range of temperature.     

An Anyon Model in a Toric Honeycomb Lattice

We study an anyon model in a toric honeycomb lattice. The ground states and the low-lying excitations coincide with those of Kitaev toric code model and then the excitations obey mutual semionic statistics. This model is helpful to understand the toric code of anyons in a more symmetric way. On the other hand, there is a direct relation between this toric honeycomb model and a boundary coupled Ising chain array in a square lattice via Jordan-Wignertransformation. We discuss the equivalence between these two modelsin the low-lying sector and realize these anyon excitations in a conventional fermion system. The analysis for the ground state degeneracy in the last section can also be thought of as a complementarity of our previous work [Phys. A: Math. Theor. 43 (2010) 105306].     

A Lattice Boltzmann Model for Two-Dimensional Magnetohydrodynamics

A lattice Boltzmann model （LBM） has been developed for simulating magnetohydrodynamics （MHD） along the line of Dellar [J. Comput. Phys. 179 （2002）95]. In this model the magnetic field is presented by a vector valued magnetic distribution function which obeys a vector Boltzmann equation. The truncated error of the equilibrium distribution in the present model is up to order O（u^4） in velocity u rather than the usual 0（u^3）. For verification, the model is applied to solve the shock tube problem and the main features of the flow predicted by the model are found to compare well with the corresponding results obtained with high-order semi-discrete schemes [J. Comput. Phys. 201 （2004） 261]. The numerical experiments have also shown that the present LBM model with the equilibrium distribution truncated at O（u^4） performs much better in terms of numerical stability than those truncated at O（u^3）.     

The Colored Hofstadter Butterfly for the Honeycomb Lattice

We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number \(\sigma _\mathrm{H}\) is given as the winding number of an eigenvector of a \(2 \times 2\) transfer matrix, as a function of the quasi-momentum \(k\in (0,2\pi )\) . This method is computationally efficient (of order \(\mathcal {O}(n^4)\) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for \(\sigma _\mathrm{H}\) for flux \(p/q\) in the \(r\) -th gap conforms with the Diophantine equation \(r=\sigma _\mathrm{H}\cdot p+ s\cdot q\) , which determines \(\sigma _\mathrm{H}\mod q\) . A window such as \(\sigma _\mathrm{H}\in (-q/2,q/2)\) , or possibly shifted, provides a natural further condition for \(\sigma _\mathrm{H}\) , which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition \(\sigma _\mathrm{H}\in (-q,q)\) .     

Exact Potts Model Partition Functions for Strips of the Honeycomb Lattice

We present exact calculations of the Potts model partition function Z(G,q,v) for arbitrary q and temperature-like variable v on strip graphs G of the honeycomb lattice for a variety of transverse widths equal to L _y vertices and for arbitrarily great length, with free longitudinal boundary conditions and free and periodic transverse boundary conditions. These partition functions have the form , where m denotes the number of repeated subgraphs in the longitudinal direction. We give general formulas for N _Z,G,j for arbitrary L _y . We also present plots of zeros of the partition function in the q plane for various values of v and in the v plane for various values of q. Plots of specific heat for infinite-length strips are also presented, and, in particular, the behavior of the Potts antiferromagnet at is investigated.     

Theory of Photo-emission from the Two-Dimensional t-J Model on a Triangular Lattice

The electron spectral function, the electron dispersion, and the electron density of states of the two-dimensional t-J model on a triangular lattice are studied within the selfconsistent fermion-spin mean-field theory. It is shown that there is a gap in the electron dispersions, and the global features of the electron spectral function is qualitatively similar to the result of the t-J model on a square lattice.     

Weak Universality in the Disordered Two-Dimensional Antiferromagnetic Potts Model on a Triangular Lattice

The critical behavior of the disordered two-dimensional antiferromagnetic Potts model with the number of spin states q= 3 on a triangular lattice with disorder in the form of nonmagnetic impurities is studied by the Monte Carlo method. The critical exponents for the susceptibility γ, magnetization β, specific heat α, and correlation radius ν are calculated in the framework of the finite-size scaling theory at spin concentrations p = 0.90, 0.80, 0.70, and 0.65. It is found that the critical exponents increase with the degree of disorder, whereas the ratios and do not change, thus holding the scaling equality \(\frac{{2\beta }}{\nu } + \frac{\gamma }{\nu } = d\). Such behavior of the critical exponents is related to the weak universality of the critical behavior characteristic of disordered systems. All results are obtained using independent Monte Carlo algorithms, such as the Metropolis and Wolff algorithms.     

Drude Weight,Optical Conductivity of Two-Dimensional Hubbard Model at Half Filling

We study the Drude weight D and optical conductivity of the two-dimensional (2D) Hubbard model at half filling with staggered magnetic flux (SMF). When SMF being introduced, the hopping integrals are modulated by the magnetic flux. The optical sum rule, which is related to the mean kinetic energy of band electrons, is evaluated for this 2D Hubbard Hamiltonian. Our present result gives the dependence of the kinetic energy, D and the optical conductivity on SMF and U. At half filling D vanishes exponentially with system size. We also find in the frequency dependence of the optical conductivity, there is δ-function peak at ω≈2|m|U and the incoherent excitations begin to present themselves extended to a higher energy region.     

Drude Weight, Optical Conductivity of Two-Dimensional Hubbard Model at Half Filling

We study the Drude weight D and optical conductivity of the two-dimensional （2D） Hubbard model at half filling with staggered magnetic flux （SMF）. When SMF being introduced, the hopping integrals are modulated by the magnetic flux. The optical sum rule, which is related to the mean kinetic energy of band electrons, is evaluated for this 2D Hubbard Hamiltonian. Our present result gives the dependence of the kinetic energy, D and the optical conductivity on SMF and U. At half filling D vanishes exponentially with system size. We also find in the frequency dependence of the optical conductivity, there is δ-function peak at ω ≈ 2｜m｜U and the incoherent excitations begin to present themselves extended to a higher energy region.     

Berezinskii-Kosterlitz-Thouless Transition in a Two-Dimensional Random-Bond XY Model on a Square Lattice

We perform Monte Carlo simulations to study the two dimensional random-bond XY model on a square lattice. Two kinds of bond randomness with the coupling coefficient obeying the Gaussian or uniform distribution are discussed. It is shown that the two kinds of disorders lead to similar thermodynamic behaviors if their variances take the same value. This result implies that the variance can be chosen as a characteristic parameter to evaluate the strength of the randomness. In addition, the Berezinskii-Kosterlitz-Thouless transition temperature decreases as the variance increases and the transition can even be destroyed as long as the disorder is strong enough.     

二维正方晶格Compass-XY模型的方向序和拓扑序

采用Monte Carlo方法研究二维正方晶格Compass-XY模型,通过调节过渡参量α,计算了降低交换相互作用的阻挫对方向序和拓扑序的影响.结果表明,方向序和拓扑序之间的转变是一个逐渐过渡的过程,并不存在明显的界限.在某一范围内,这两种序是相互重叠的,难以区分.阻挫的降低,易于形成拓扑序,同时会抑制方向序.当阻挫变得足够弱时,方向序被破坏.     

三角晶格二维光子晶体的横磁模特性

基于平面波展开法,以介电常数分别为1F/m和12F/m的物质构成三角晶格二维光子晶体,改变空气孔半径r与晶格常数a之间的大小,数值模拟得到了三角晶格二维光子晶体横磁模带隙,当r=.4α时形成三处三角晶格二维光子晶体横磁模的带隙,其中最大三角晶格二维光子晶体横磁模带隙出现在0.6766—0.8000Hz,差值为0.1234Hz。研究结论为光子晶体器件的制作提供参考。     

Highly Efficient Lattice Boltzmann Model for Compressible Fluids:Two-Dimensional Case

We present a highly efficient lattice Boltzmann model for simulatingcompressible flows. This model is based on the combination of an appropriatefinite difference scheme, a 16-discrete-velocity model [Kataoka andTsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion anddissipation terms. The dispersion term effectively reduces the oscillationat the discontinuity and enhances numerical precision. The dissipation termmakes the new model more easily meet with the von Neumann stabilitycondition. This model works for both high-speed and low-speed flows witharbitrary specific-heat-ratio. With the new model simulation results for thewell-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation ofstate or free-energy functional, the new model has the potential tostudythe complex procedure of shock wave reaction on porous materials.