We consider a quadrature-based eigensolver to find eigenpairs of Hermitian matrices arising in lattice quantum chromodynamics. To reduce the computational cost for finding eigenpairs of such Hermitian matrices, we propose a new technique for solving
shifted linear systems with complex shifts by means of the shifted CG method. Furthermore using integration paths along horizontal lines corresponding to the real axis of the complex plane, the number of iterations for the shifted CG method is also reduced. Some numerical experiments illustrate the accuracy and efficiency of the proposed method by comparison with a conventional method.
We report on the status of our study towards the equation of state in 2+1 flavor QCD with improved Wilson quarks. To reduce the computational cost which is quite demanding for Wilson-type quarks, we adopt the fixed scale approach, i.e. the temperatur
e T is varied by N_t at fixed lattice spacing. Since the conventional integral method to obtain the pressure is inapplicable at a fixed scale, we adopt the T-integral method, to calculate the pressure non-perturbatively. Reduction of the computational cost of T=0 simulations thus achieved is indispensable to study EOS in QCD with dynamical quarks.
Dissociation temperatures of J/psi, psi, and chi_c states play key roles in the sequential J/psi suppression scenario for high energy heavy ion collisions. We report on a study of charmonium dissociation temperatures in quenched lattice QCD. On aniso
tropic lattices, we first subtract the effects of the constant mode in finite temperature meson correlators, which have lead to unphysical results in previous studies. We then extract ground and first exited state masses by diagonalizing correlation functions among different source and sink operators. To distinguish bound states from scattering states, we first compare the charmonium mass spectra under different spatial boundary conditions, and examine the shape and the volume-dependence of their Bethe-Salpeter wave functions. From these studies, we found so far no sign of scattering states up to about 2.3T_c.
