Improved lattice QCD with quarks: the 2 dimensional case

QCD in two dimensions is investigated using the improved fermionic lattice Hamiltonian proposed by Luo, Chen, Xu, and Jiang. We show that the improved theory leads to a significant reduction of the finite lattice spacing errors. The quark condensate and the mass of lightest quark and anti-quark bound state in the strong coupling phase (different from t'Hooft phase) are computed. We find agreement between our results and the analytical ones in the continuum.Comment: LaTeX file (including text + 10 figures

Thermodynamical quantities of lattice full QCD from an efficient method

I extend to QCD an efficient method for lattice gauge theory with dynamical fermions. Once the eigenvalues of the Dirac operator and the density of states of pure gluonic configurations at a set of plaquette energies (proportional to the gauge action) are computed, thermodynamical quantities deriving from the partition function can be obtained for arbitrary flavor number, quark masses and wide range of coupling constants, without additional computational cost. Results for the chiral condensate and gauge action are presented on the $10^4$ lattice at flavor number $N_f=0$, 1, 2, 3, 4 and many quark masses and coupling constants. New results in the chiral limit for the gauge action and its correlation with the chiral condensate, which are useful for analyzing the QCD chiral phase structure, are also provided.Comment: Latex, 11 figures, version accepted for publicatio

Bound States and Critical Behavior of the Yukawa Potential

We investigate the bound states of the Yukawa potential $V(r)=-\lambda \exp(-\alpha r)/ r$, using different algorithms: solving the Schr\"odinger equation numerically and our Monte Carlo Hamiltonian approach. There is a critical $\alpha=\alpha_C$, above which no bound state exists. We study the relation between $\alpha_C$ and $\lambda$ for various angular momentum quantum number $l$, and find in atomic units, $\alpha_{C}(l)= \lambda [A_{1} \exp(-l/ B_{1})+ A_{2} \exp(-l/ B_{2})]$, with $A_1=1.020(18)$, $B_1=0.443(14)$, $A_2=0.170(17)$, and $B_2=2.490(180)$.Comment: 15 pages, 12 figures, 5 tables. Version to appear in Sciences in China

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Li et al., Phys. Rev. E. 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/sqrt(a). Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.Comment: 9 figures, 2 tables, accepted by Physical Review E (in press

Matrix convex functions with applications to weighted centers for semidefinite programming

In this paper, we develop various calculus rules for general smooth matrix-valued functions and for the class of matrix convex (or concave) functions first introduced by Loewner and Kraus in 1930s. Then we use these calculus rules and the matrix convex function -log X to study a new notion of weighted convex centers for semidefinite programming (SDP) and show that, with this definition, some known properties of weighted centers for linear programming can be extended to SDP. We also show how the calculus rules for matrix convex functions can be used in the implementation of barrier methods for optimization problems involving nonlinear matrix functions.matrix convexity;matrix monotonicity;semidefinite programming