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

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

Exotic mesons from quantum chromodynamics with improved gluon and quark actions on the anisotropic lattice

Hybrid (exotic) mesons, which are important predictions of quantum chromodynamics (QCD), are states of quarks and anti-quarks bound by excited gluons. First principle lattice study of such states would help us understand the role of ``dynamical'' color in low energy QCD and provide valuable information for experimental search for these new particles. In this paper, we apply both improved gluon and quark actions to the hybrid mesons, which might be much more efficient than the previous works in reducing lattice spacing error and finite volume effect. Quenched simulations were done at $\beta=2.6$ and on a $\xi=3$ anisotropic $12^3\times36$ lattice using our PC cluster. We obtain $2013 \pm 26 \pm 71$ MeV for the mass of the $1^{-+}$ hybrid meson ${\bar q}qg$ in the light quark sector, and $4369 \pm 37 \pm 99$Mev in the charm quark sector; the mass splitting between the $1^{-+}$ hybrid meson ${\bar c}c g$ in the charm quark sector and the spin averaged S-wave charmonium mass is estimated to be $1302 \pm 37 \pm 99$ MeV. As a byproduct, we obtain $1438 \pm 32 \pm 57$ MeV for the mass of a P-wave $1^{++}$ ${\bar u}u$ or ${\bar d}d$ meson and $1499 \pm 28 \pm 65$ MeV for the mass of a P-wave $1^{++}$ ${\bar s}s$ meson, which are comparable to their experimental value 1426 MeV for the $f_1(1420)$ meson. The first error is statistical, and the second one is systematical. The mixing of the hybrid meson with a four quark state is also discussed.Comment: 12 pages, 3 figures. Published versio

Algorithm for Computing Excited States in Quantum Theory

Monte Carlo techniques have been widely employed in statistical physics as well as in quantum theory in the Lagrangian formulation. However, in the conventional approach, it is extremely difficult to compute the excited states. Here we present a different algorithm: the Monte Carlo Hamiltonian method, designed to overcome the difficulties of the conventional approach. As a new example, application to the Klein-Gordon field theory is shown.Comment: 3 pages, uses Latex and aipproc.cl

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

Comment on "Time-Dependent Density-Matrix Renormalization Group: A Systematic Method for the Study of Quantum Many-Body Out-of- Equilibrium Systems"

In a recent Letter [Phys. Rev. Lett. 88, 256403(2002), cond-mat/0109158] Cazalilla and Marston proposed a time-dependent density- matrix renormalization group (TdDMRG) algorithm for the accurate evaluation of out-of-equilibrium properties of quantum many-body systems. For a point contact junction between two Luttinger liquids, a current oscillation develops after initial transient in the insulating regime. Here we would like to point out that (a) the observed oscillation is an artifact of the method; (b) the TdDMRG can be significantly improved by incorporating the non-equilibrium evolution of the goundstate into the density matrix.Comment: 1 page, 2 figure